ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Mods/interface/PhotonIDMod.h
Revision: 1.28
Committed: Wed Jul 25 15:00:42 2012 UTC (12 years, 9 months ago) by fabstoec
Content type: text/plain
Branch: MAIN
Changes since 1.27: +7 -6 lines
Log Message:
added Vgamma Photon ID

File Contents

# User Rev Content
1 sixie 1.1 //--------------------------------------------------------------------------------------------------
2 fabstoec 1.28 // $Id: PhotonIDMod.h,v 1.27 2012/07/24 11:41:22 fabstoec Exp $
3 sixie 1.1 //
4     // PhotonIDMod
5     //
6     // This module applies photon identification criteria and exports a pointer to a collection
7     // of "good photons" according to the specified identification scheme.
8     //
9     // Authors: S.Xie, C.Loizides
10     //--------------------------------------------------------------------------------------------------
11    
12     #ifndef MITPHYSICS_MODS_PHOTONIDMOD_H
13     #define MITPHYSICS_MODS_PHOTONIDMOD_H
14    
15     #include "MitAna/TreeMod/interface/BaseMod.h"
16 loizides 1.10 #include "MitAna/DataTree/interface/PhotonFwd.h"
17 fabstoec 1.15 #include "MitAna/DataTree/interface/TrackCol.h"
18     #include "MitAna/DataTree/interface/BeamSpotCol.h"
19     #include "MitAna/DataTree/interface/PileupEnergyDensityCol.h"
20 bendavid 1.16 #include "MitAna/DataTree/interface/DecayParticleCol.h"
21     #include "MitAna/DataTree/interface/ElectronCol.h"
22 fabstoec 1.20 #include "MitAna/DataTree/interface/VertexCol.h"
23 bendavid 1.22 #include "MitAna/DataTree/interface/DecayParticleCol.h"
24     #include "MitAna/DataTree/interface/PileupInfoCol.h"
25     #include "MitAna/DataTree/interface/MCParticleCol.h"
26 bendavid 1.26 #include "MitAna/DataTree/interface/PFCandidateCol.h"
27 sixie 1.1
28 fabstoec 1.23 #include "MitPhysics/Utils/interface/MVATools.h"
29 fabstoec 1.27 #include "MitPhysics/Utils/interface/PhotonTools.h"
30    
31     class TRandom3;
32 fabstoec 1.23
33 sixie 1.1 namespace mithep
34     {
35     class PhotonIDMod : public BaseMod
36     {
37     public:
38     PhotonIDMod(const char *name="PhotonIDMod",
39 loizides 1.3 const char *title="Photon identification module");
40 sixie 1.1
41 sixie 1.13 Bool_t GetApplySpikeRemoval() const { return fApplySpikeRemoval; }
42 loizides 1.6 Bool_t GetApplyPixelSeed() const { return fApplyPixelSeed; }
43     const char *GetGoodName() const { return GetGoodPhotonsName(); }
44     const char *GetGoodPhotonsName() const { return fGoodPhotonsName; }
45     Double_t GetHadOverEmMax() const { return fHadOverEmMax; }
46     const char *GetIDType() const { return fPhotonIDType; }
47     const char *GetInputName() const { return fPhotonBranchName; }
48     const char *GetIsoType() const { return fPhotonIsoType; }
49     const char *GetOutputName() const { return GetGoodPhotonsName(); }
50     Double_t GetPtMin() const { return fPhotonPtMin; }
51 ceballos 1.11 Bool_t GetApplyFiduciality() const { return fFiduciality; }
52     Double_t GetEtaWidthEB() const { return fEtaWidthEB; }
53     Double_t GetEtaWidthEE() const { return fEtaWidthEE; }
54 ceballos 1.12 Double_t GetAbsEtaMax() const { return fAbsEtaMax; }
55 sixie 1.13 void SetApplySpikeRemoval(Bool_t b) { fApplySpikeRemoval = b; }
56 loizides 1.6 void SetApplyPixelSeed(Bool_t b) { fApplyPixelSeed = b; }
57 bendavid 1.16 void SetApplyElectronVeto(Bool_t b) { fApplyElectronVeto = b; }
58 bendavid 1.22 void SetInvertElectronVeto(Bool_t b) { fInvertElectronVeto = b; }
59 bendavid 1.16 void SetApplyElectronVetoConvRecovery(Bool_t b) { fApplyElectronVetoConvRecovery = b; }
60     void SetApplyConversionId(Bool_t b) { fApplyConversionId = b; }
61     void SetApplyTriggerMatching(Bool_t b) { fApplyTriggerMatching = b; }
62 loizides 1.6 void SetGoodName(const char *n) { SetGoodPhotonsName(n); }
63     void SetGoodPhotonsName(const char *n) { fGoodPhotonsName = n; }
64     void SetHadOverEmMax(Double_t hoe) { fHadOverEmMax = hoe; }
65     void SetIDType(const char *type) { fPhotonIDType = type; }
66     void SetInputName(const char *n) { fPhotonBranchName= n; }
67 fabstoec 1.15 void SetTrackName(const char *n) { fTrackBranchName = n; }
68     void SetBeamspotName(const char *n) { fBeamspotBranchName = n; }
69 loizides 1.6 void SetIsoType(const char *type) { fPhotonIsoType = type; }
70     void SetOutputName(const char *n) { SetGoodPhotonsName(n); }
71     void SetPtMin(Double_t pt) { fPhotonPtMin = pt; }
72 ceballos 1.11 void SetR9Min(Double_t x) { fPhotonR9Min = x; }
73     void SetEtaWidthEB(Double_t x) { fEtaWidthEB = x; }
74     void SetEtaWidthEE(Double_t x) { fEtaWidthEE = x; }
75 ceballos 1.12 void SetAbsEtaMax(Double_t x) { fAbsEtaMax = x; }
76 bendavid 1.14 void SetApplyR9Min(Bool_t b) { fApplyR9Min = b; }
77 bendavid 1.18 void SetApplyFiduciality(Bool_t b) { fFiduciality = b; }
78 fabstoec 1.15 void SetEffAreas(Double_t ecal, Double_t hcal, Double_t track) {
79 fabstoec 1.17 fEffAreaEcalEE = ecal; fEffAreaHcalEE = hcal; fEffAreaTrackEE = track;
80     fEffAreaEcalEB = ecal; fEffAreaHcalEB = hcal; fEffAreaTrackEB = track;
81     }
82     void SetEffAreasEEEB(Double_t ecalEE, Double_t hcalEE, Double_t trackEE,
83     Double_t ecalEB, Double_t hcalEB, Double_t trackEB) {
84     fEffAreaEcalEE = ecalEE; fEffAreaHcalEE = hcalEE; fEffAreaTrackEE = trackEE;
85     fEffAreaEcalEB = ecalEB; fEffAreaHcalEB = hcalEB; fEffAreaTrackEB = trackEB;
86     }
87 bendavid 1.16 void SetTriggerObjectsName(const char *n) { fTrigObjectsName = n; }
88 fabstoec 1.15
89 fabstoec 1.25
90     void SetPhotonsFromBranch(bool b) { fPhotonsFromBranch = b; }
91    
92 fabstoec 1.20 void SetPVName(const char *n) { fPVName = n; }
93     void SetPVFromBranch(bool b) { fPVFromBranch = b; }
94 bendavid 1.21 void SetIsData (Bool_t b) { fIsData = b;};
95 fabstoec 1.27
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 bendavid 1.22 void SetGoodElectronsFromBranch(Bool_t b) { fGoodElectronsFromBranch = b; }
145     void SetGoodElectronName(TString name) { fGoodElectronName = name; }
146 fabstoec 1.20
147 fabstoec 1.23 void SetBdtCutBarrel(double a) {fbdtCutBarrel = a; }
148     void SetBdtCutEndcap(double a) {fbdtCutEndcap = a; }
149    
150 fabstoec 1.27
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 fabstoec 1.25 void SetDoMCErrScaling(Bool_t b) { fDoMCErrScaling = b; }
162     void SetMCErrScale(Double_t ebscale, Double_t eescale) { fMCErrScaleEB = ebscale; fMCErrScaleEE = eescale; }
163 fabstoec 1.23
164 sixie 1.1 enum EPhIdType {
165     kIdUndef = 0, //not defined
166     kTight, //"Tight"
167     kLoose, //"Loose"
168 ceballos 1.2 kLooseEM, //"LooseEM"
169 fabstoec 1.28 kBaseLineCiC, //"2011" Hgg BaseLine CiC
170     kBaseLineCiCPF, //"2012" Hgg BaseLine CiC
171     kMITMVAId, // MingMing MVA ID
172     kMITPhSelection, //MIT loose preselection (for mva)
173 bendavid 1.26 kMITPFPhSelection, //MIT loose preselection (for mva)
174 fabstoec 1.28 kVgamma2011Selection, // Vgamma 2011 Photon ID
175     kCustomId //"Custom"
176 sixie 1.1 };
177 fabstoec 1.19
178 sixie 1.1 enum EPhIsoType {
179     kIsoUndef = 0, //not defined
180     kNoIso, //"NoIso"
181 ceballos 1.2 kCombinedIso, //"CombinedIso"
182 fabstoec 1.15 kCustomIso, //"Custom"
183     kMITPUCorrected //PileUp Corrected Hgg Isolation
184 sixie 1.1 };
185    
186     protected:
187 loizides 1.5 void Process();
188     void SlaveBegin();
189    
190 fabstoec 1.27 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 loizides 1.6 TString fPhotonBranchName; //name of photon collection (input)
195     TString fGoodPhotonsName; //name of exported "good photon" collection
196 fabstoec 1.15 TString fTrackBranchName; // name of the track collection (only needed for PU corrected isolation)
197     TString fBeamspotBranchName; //name of the Beamspot collection (only needed for PU corrected isolation)
198 bendavid 1.16 TString fPileUpDenName; //name of the PU density collection
199     TString fConversionName; //name of conversion branch
200     TString fElectronName;
201 bendavid 1.22 TString fGoodElectronName;
202 fabstoec 1.20 TString fTrigObjectsName; //name of trigger object collection
203     TString fPVName;
204 bendavid 1.21 TString fMCParticleName;
205     TString fPileUpName;
206 bendavid 1.26 TString fPFCandsName;
207 loizides 1.6 TString fPhotonIDType; //type of photon identification we impose
208     TString fPhotonIsoType; //type of photon isolation we impose
209 loizides 1.4 Double_t fPhotonPtMin; //min pt cut
210 loizides 1.6 Double_t fHadOverEmMax; //maximum of hadronic/em energy
211 sixie 1.13 Bool_t fApplySpikeRemoval; //whether apply spike removal
212 loizides 1.6 Bool_t fApplyPixelSeed; //=true then apply pixel seed constraint
213 bendavid 1.16 Bool_t fApplyElectronVeto; //=true then apply electron veto (with no conversion recovery)
214 bendavid 1.22 Bool_t fInvertElectronVeto; //=true then invert electron veto (for cic selection only atm)
215 bendavid 1.16 Bool_t fApplyElectronVetoConvRecovery; //=true then apply electron veto with conversion recovery
216     Bool_t fApplyConversionId; //=true then apply conversion id cuts
217     Bool_t fApplyTriggerMatching; //match to hlt photon (default=0)
218 ceballos 1.8 Double_t fPhotonR9Min; //min R9 value
219 loizides 1.4 EPhIdType fPhIdType; //!identification scheme
220     EPhIsoType fPhIsoType; //!isolation scheme
221 ceballos 1.11 Bool_t fFiduciality; //=true then apply fiducual requirement
222 fabstoec 1.27
223 ceballos 1.11 Double_t fEtaWidthEB; //max Eta Width in ECAL Barrel
224     Double_t fEtaWidthEE; //max Eta Width in ECAL End Cap
225 ceballos 1.12 Double_t fAbsEtaMax; //max Abs Eta
226 bendavid 1.14 Bool_t fApplyR9Min; //apply R9 min
227 fabstoec 1.17 Double_t fEffAreaEcalEE;
228     Double_t fEffAreaHcalEE;
229     Double_t fEffAreaTrackEE;
230     Double_t fEffAreaEcalEB;
231     Double_t fEffAreaHcalEB;
232     Double_t fEffAreaTrackEB;
233 loizides 1.10 const PhotonCol *fPhotons; //!photon branch
234 fabstoec 1.15 const TrackCol *fTracks; //!track branch
235     const BeamSpotCol *fBeamspots; //!beamspot branch
236 bendavid 1.16 const PileupEnergyDensityCol *fPileUpDen; //!rho branch
237     const DecayParticleCol *fConversions; //!conversion branch
238     const ElectronCol *fElectrons; //!electron branch
239 bendavid 1.22 const ElectronCol *fGoodElectrons; //!electron branch
240 bendavid 1.21 const VertexCol* fPV; //!
241     const MCParticleCol *fMCParticles;//!
242     const PileupInfoCol *fPileUp; //!
243 bendavid 1.26 const PFCandidateCol *fPFCands; //!
244 fabstoec 1.23
245     Double_t fbdtCutBarrel;
246     Double_t fbdtCutEndcap;
247     int fVariableType;
248     TString fEndcapWeights;
249     TString fBarrelWeights;
250     MVATools fTool;
251    
252 fabstoec 1.25 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 bendavid 1.22 Bool_t fGoodElectronsFromBranch;
266 bendavid 1.21 Bool_t fIsData;
267    
268    
269 fabstoec 1.27 // 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 fabstoec 1.15
305 loizides 1.6 ClassDef(PhotonIDMod, 1) // Photon identification module
306 sixie 1.1 };
307     }
308     #endif