| 23 |
|
#include "MitAna/DataTree/interface/DecayParticleCol.h" |
| 24 |
|
#include "MitAna/DataTree/interface/PileupInfoCol.h" |
| 25 |
|
#include "MitAna/DataTree/interface/MCParticleCol.h" |
| 26 |
+ |
#include "MitAna/DataTree/interface/PFCandidateCol.h" |
| 27 |
|
|
| 28 |
|
#include "MitPhysics/Utils/interface/MVATools.h" |
| 29 |
+ |
#include "MitPhysics/Utils/interface/PhotonTools.h" |
| 30 |
+ |
|
| 31 |
+ |
class TRandom3; |
| 32 |
|
|
| 33 |
|
namespace mithep |
| 34 |
|
{ |
| 86 |
|
} |
| 87 |
|
void SetTriggerObjectsName(const char *n) { fTrigObjectsName = n; } |
| 88 |
|
|
| 89 |
+ |
|
| 90 |
+ |
void SetPhotonsFromBranch(bool b) { fPhotonsFromBranch = b; } |
| 91 |
+ |
|
| 92 |
|
void SetPVName(const char *n) { fPVName = n; } |
| 93 |
|
void SetPVFromBranch(bool b) { fPVFromBranch = b; } |
| 94 |
|
void SetIsData (Bool_t b) { fIsData = b;}; |
| 95 |
+ |
|
| 96 |
+ |
|
| 97 |
+ |
void SetShowerShapeType(const char *type) { fShowerShapeType = type; } |
| 98 |
+ |
|
| 99 |
+ |
// methods to set the MC smearing/energy correction values |
| 100 |
+ |
void AddEnCorrPerRun( UInt_t sRun, UInt_t eRun, |
| 101 |
+ |
Double_t corr_EBlowEta_hR9central, |
| 102 |
+ |
Double_t corr_EBlowEta_hR9gap, |
| 103 |
+ |
Double_t corr_EBlowEta_lR9, |
| 104 |
+ |
Double_t corr_EBhighEta_hR9, |
| 105 |
+ |
Double_t corr_EBhighEta_lR9, |
| 106 |
+ |
Double_t corr_EElowEta_hR9, |
| 107 |
+ |
Double_t corr_EElowEta_lR9, |
| 108 |
+ |
Double_t corr_EEhighEta_hR9, |
| 109 |
+ |
Double_t corr_EEhighEta_lR9) { |
| 110 |
+ |
|
| 111 |
+ |
fDataEnCorr_EBlowEta_hR9central.push_back(corr_EBlowEta_hR9central); |
| 112 |
+ |
fDataEnCorr_EBlowEta_hR9gap.push_back(corr_EBlowEta_hR9gap); |
| 113 |
+ |
fDataEnCorr_EBlowEta_lR9.push_back(corr_EBlowEta_lR9); |
| 114 |
+ |
fDataEnCorr_EBhighEta_hR9.push_back(corr_EBhighEta_hR9); |
| 115 |
+ |
fDataEnCorr_EBhighEta_lR9.push_back(corr_EBhighEta_lR9); |
| 116 |
+ |
fDataEnCorr_EElowEta_hR9.push_back(corr_EElowEta_hR9); |
| 117 |
+ |
fDataEnCorr_EElowEta_lR9.push_back(corr_EElowEta_lR9); |
| 118 |
+ |
fDataEnCorr_EEhighEta_hR9.push_back(corr_EEhighEta_hR9); |
| 119 |
+ |
fDataEnCorr_EEhighEta_lR9.push_back(corr_EEhighEta_lR9); |
| 120 |
+ |
fRunStart.push_back (sRun); |
| 121 |
+ |
fRunEnd.push_back (eRun); |
| 122 |
+ |
}; |
| 123 |
+ |
|
| 124 |
+ |
void SetMCSmearFactors(Double_t _EBlowEta_hR9central, |
| 125 |
+ |
Double_t _EBlowEta_hR9gap, |
| 126 |
+ |
Double_t _EBlowEta_lR9, |
| 127 |
+ |
Double_t _EBhighEta_hR9, |
| 128 |
+ |
Double_t _EBhighEta_lR9, |
| 129 |
+ |
Double_t _EElowEta_hR9, |
| 130 |
+ |
Double_t _EElowEta_lR9, |
| 131 |
+ |
Double_t _EEhighEta_hR9, |
| 132 |
+ |
Double_t _EEhighEta_lR9) { |
| 133 |
+ |
fMCSmear_EBlowEta_hR9central = _EBlowEta_hR9central; |
| 134 |
+ |
fMCSmear_EBlowEta_hR9gap = _EBlowEta_hR9gap; |
| 135 |
+ |
fMCSmear_EBlowEta_lR9 = _EBlowEta_lR9; |
| 136 |
+ |
fMCSmear_EBhighEta_hR9 = _EBhighEta_hR9; |
| 137 |
+ |
fMCSmear_EBhighEta_lR9 = _EBhighEta_lR9; |
| 138 |
+ |
fMCSmear_EElowEta_hR9 = _EElowEta_hR9; |
| 139 |
+ |
fMCSmear_EElowEta_lR9 = _EElowEta_lR9; |
| 140 |
+ |
fMCSmear_EEhighEta_hR9 = _EEhighEta_hR9; |
| 141 |
+ |
fMCSmear_EEhighEta_lR9 = _EEhighEta_lR9; |
| 142 |
+ |
}; |
| 143 |
+ |
|
| 144 |
|
void SetGoodElectronsFromBranch(Bool_t b) { fGoodElectronsFromBranch = b; } |
| 145 |
|
void SetGoodElectronName(TString name) { fGoodElectronName = name; } |
| 146 |
|
|
| 148 |
|
void SetBdtCutEndcap(double a) {fbdtCutEndcap = a; } |
| 149 |
|
|
| 150 |
|
|
| 151 |
+ |
void DoDataEneCorr(bool a) { fDoDataEneCorr = a; } |
| 152 |
+ |
void DoMCSmear(bool a) { fDoMCSmear = a; } |
| 153 |
+ |
void SetDoShowerShapeScaling(Bool_t b) { fDoShowerShapeScaling = b; } |
| 154 |
+ |
|
| 155 |
+ |
// replaced by ShowerShapeScaling (fab) |
| 156 |
+ |
/* void SetDoMCR9Scaling(Bool_t b) { fDoMCR9Scaling = b; } */ |
| 157 |
+ |
/* void SetMCR9Scale(Double_t ebscale, Double_t eescale) { fMCR9ScaleEB = ebscale; fMCR9ScaleEE = eescale; } */ |
| 158 |
+ |
/* void SetDoMCSigIEtaIEtaScaling(Bool_t b) { fDoMCSigIEtaIEtaScaling = b; } */ |
| 159 |
+ |
/* void SetDoMCWidthScaling(Bool_t b) { fDoMCWidthScaling = b; } */ |
| 160 |
+ |
|
| 161 |
+ |
void SetDoMCErrScaling(Bool_t b) { fDoMCErrScaling = b; } |
| 162 |
+ |
void SetMCErrScale(Double_t ebscale, Double_t eescale) { fMCErrScaleEB = ebscale; fMCErrScaleEE = eescale; } |
| 163 |
+ |
|
| 164 |
|
enum EPhIdType { |
| 165 |
|
kIdUndef = 0, //not defined |
| 166 |
|
kTight, //"Tight" |
| 167 |
|
kLoose, //"Loose" |
| 168 |
|
kLooseEM, //"LooseEM" |
| 169 |
< |
kBaseLineCiC, //"2011" Hgg BaseLine CiC |
| 170 |
< |
kMITMVAId, // MingMing MVA ID |
| 171 |
< |
kCustomId //"Custom" |
| 169 |
> |
kBaseLineCiC, //"2011" Hgg BaseLine CiC |
| 170 |
> |
kBaseLineCiCPF, //"2012" Hgg BaseLine CiC |
| 171 |
> |
kMITMVAId, // MingMing MVA ID |
| 172 |
> |
kMITPhSelection, //MIT loose preselection (for mva) |
| 173 |
> |
kMITPFPhSelection, //MIT loose preselection (for mva) |
| 174 |
> |
kVgamma2011Selection, // Vgamma 2011 Photon ID |
| 175 |
> |
kCustomId //"Custom" |
| 176 |
|
}; |
| 177 |
|
|
| 178 |
|
enum EPhIsoType { |
| 187 |
|
void Process(); |
| 188 |
|
void SlaveBegin(); |
| 189 |
|
|
| 190 |
+ |
Int_t FindRunRangeIdx(UInt_t run); |
| 191 |
+ |
Double_t GetDataEnCorr(Int_t runRange, PhotonTools::eScaleCats cat); |
| 192 |
+ |
Double_t GetMCSmearFac(PhotonTools::eScaleCats cat); |
| 193 |
+ |
|
| 194 |
|
TString fPhotonBranchName; //name of photon collection (input) |
| 195 |
|
TString fGoodPhotonsName; //name of exported "good photon" collection |
| 196 |
|
TString fTrackBranchName; // name of the track collection (only needed for PU corrected isolation) |
| 203 |
|
TString fPVName; |
| 204 |
|
TString fMCParticleName; |
| 205 |
|
TString fPileUpName; |
| 206 |
+ |
TString fPFCandsName; |
| 207 |
|
TString fPhotonIDType; //type of photon identification we impose |
| 208 |
|
TString fPhotonIsoType; //type of photon isolation we impose |
| 209 |
|
Double_t fPhotonPtMin; //min pt cut |
| 219 |
|
EPhIdType fPhIdType; //!identification scheme |
| 220 |
|
EPhIsoType fPhIsoType; //!isolation scheme |
| 221 |
|
Bool_t fFiduciality; //=true then apply fiducual requirement |
| 222 |
+ |
|
| 223 |
|
Double_t fEtaWidthEB; //max Eta Width in ECAL Barrel |
| 224 |
|
Double_t fEtaWidthEE; //max Eta Width in ECAL End Cap |
| 225 |
|
Double_t fAbsEtaMax; //max Abs Eta |
| 240 |
|
const VertexCol* fPV; //! |
| 241 |
|
const MCParticleCol *fMCParticles;//! |
| 242 |
|
const PileupInfoCol *fPileUp; //! |
| 243 |
+ |
const PFCandidateCol *fPFCands; //! |
| 244 |
|
|
| 245 |
|
Double_t fbdtCutBarrel; |
| 246 |
|
Double_t fbdtCutEndcap; |
| 249 |
|
TString fBarrelWeights; |
| 250 |
|
MVATools fTool; |
| 251 |
|
|
| 252 |
< |
bool fPVFromBranch; |
| 252 |
> |
Bool_t fDoMCR9Scaling; |
| 253 |
> |
Double_t fMCR9ScaleEB; |
| 254 |
> |
Double_t fMCR9ScaleEE; |
| 255 |
> |
|
| 256 |
> |
Bool_t fDoMCSigIEtaIEtaScaling; |
| 257 |
> |
Bool_t fDoMCWidthScaling; |
| 258 |
> |
|
| 259 |
> |
Bool_t fDoMCErrScaling; |
| 260 |
> |
Double_t fMCErrScaleEB; |
| 261 |
> |
Double_t fMCErrScaleEE; |
| 262 |
> |
|
| 263 |
> |
Bool_t fPhotonsFromBranch; |
| 264 |
> |
Bool_t fPVFromBranch; |
| 265 |
|
Bool_t fGoodElectronsFromBranch; |
| 266 |
|
Bool_t fIsData; |
| 267 |
|
|
| 268 |
|
|
| 269 |
+ |
// showershape |
| 270 |
+ |
TString fShowerShapeType; |
| 271 |
+ |
PhotonTools::ShowerShapeScales fSSType; |
| 272 |
+ |
|
| 273 |
+ |
bool fDoDataEneCorr; |
| 274 |
+ |
bool fDoMCSmear; |
| 275 |
+ |
Bool_t fDoShowerShapeScaling; |
| 276 |
+ |
|
| 277 |
+ |
|
| 278 |
+ |
// Vectroes to hols smeraring/correction factors |
| 279 |
+ |
std::vector<Double_t> fDataEnCorr_EBlowEta_hR9central; |
| 280 |
+ |
std::vector<Double_t> fDataEnCorr_EBlowEta_hR9gap; |
| 281 |
+ |
std::vector<Double_t> fDataEnCorr_EBlowEta_lR9; |
| 282 |
+ |
std::vector<Double_t> fDataEnCorr_EBhighEta_hR9; |
| 283 |
+ |
std::vector<Double_t> fDataEnCorr_EBhighEta_lR9; |
| 284 |
+ |
std::vector<Double_t> fDataEnCorr_EElowEta_hR9; |
| 285 |
+ |
std::vector<Double_t> fDataEnCorr_EElowEta_lR9; |
| 286 |
+ |
std::vector<Double_t> fDataEnCorr_EEhighEta_hR9; |
| 287 |
+ |
std::vector<Double_t> fDataEnCorr_EEhighEta_lR9; |
| 288 |
+ |
|
| 289 |
+ |
std::vector<UInt_t> fRunStart; |
| 290 |
+ |
std::vector<UInt_t> fRunEnd; |
| 291 |
+ |
|
| 292 |
+ |
Double_t fMCSmear_EBlowEta_hR9central; |
| 293 |
+ |
Double_t fMCSmear_EBlowEta_hR9gap; |
| 294 |
+ |
Double_t fMCSmear_EBlowEta_lR9; |
| 295 |
+ |
Double_t fMCSmear_EBhighEta_hR9; |
| 296 |
+ |
Double_t fMCSmear_EBhighEta_lR9; |
| 297 |
+ |
Double_t fMCSmear_EElowEta_hR9; |
| 298 |
+ |
Double_t fMCSmear_EElowEta_lR9; |
| 299 |
+ |
Double_t fMCSmear_EEhighEta_hR9; |
| 300 |
+ |
Double_t fMCSmear_EEhighEta_lR9; |
| 301 |
+ |
|
| 302 |
+ |
TRandom3* fRng; |
| 303 |
+ |
|
| 304 |
|
|
| 305 |
|
ClassDef(PhotonIDMod, 1) // Photon identification module |
| 306 |
|
}; |
