ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Mods/interface/PhotonTreeWriter.h
Revision: 1.15
Committed: Mon May 14 16:38:23 2012 UTC (12 years, 11 months ago) by fabstoec
Content type: text/plain
Branch: MAIN
Changes since 1.14: +7 -4 lines
Log Message: 
improved Btagging

File Contents

# User Rev Content
1 bendavid 1.1 //--------------------------------------------------------------------------------------------------
2 fabstoec 1.15 // $Id: PhotonTreeWriter.h,v 1.14 2012/05/03 15:49:34 fabstoec Exp $
3 bendavid 1.1 //
4     // PhotonTreeWriter
5     //
6     // Authors: J. Bendavid
7     //--------------------------------------------------------------------------------------------------
8    
9     #ifndef MITPHYSICS_MODS_PHOTONTREEWRITER_H
10     #define MITPHYSICS_MODS_PHOTONTREEWRITER_H
11    
12     #include "MitAna/TreeMod/interface/BaseMod.h"
13     #include "MitAna/DataTree/interface/PhotonFwd.h"
14     #include "MitAna/DataTree/interface/TrackCol.h"
15     #include "MitAna/DataTree/interface/VertexCol.h"
16     #include "MitAna/DataTree/interface/BeamSpotCol.h"
17     #include "MitAna/DataTree/interface/PFCandidateCol.h"
18     #include "MitAna/DataTree/interface/PileupEnergyDensityCol.h"
19     #include "MitAna/DataTree/interface/DecayParticleCol.h"
20     #include "MitAna/DataTree/interface/ElectronCol.h"
21 fabstoec 1.10 #include "MitAna/DataTree/interface/MuonCol.h"
22 bendavid 1.1 #include "MitAna/DataTree/interface/DecayParticleCol.h"
23     #include "MitAna/DataTree/interface/PileupInfoCol.h"
24     #include "MitAna/DataTree/interface/MCParticleCol.h"
25 bendavid 1.2 #include "MitAna/DataTree/interface/SuperClusterCol.h"
26     #include "MitAna/DataTree/interface/PFMetCol.h"
27 bendavid 1.5 #include "MitAna/DataTree/interface/JetCol.h"
28 fabstoec 1.13 #include "MitAna/DataTree/interface/PFJetCol.h"
29     #include "MitAna/DataTree/interface/GenJetCol.h"
30 bendavid 1.2 #include "MitPhysics/Utils/interface/PhotonFix.h"
31 bendavid 1.1 #include "MitPhysics/Utils/interface/PhotonTools.h"
32    
33     class TNtuple;
34     class TRandom3;
35    
36     namespace mithep
37     {
38    
39     class PhotonTreeWriterPhoton
40     {
41     public:
42 bendavid 1.8 void SetVars(const Photon *p, const DecayParticle *c, const Electron *ele, const SuperCluster *pfsc, const MCParticle *m, PhotonFix &phfixph, PhotonFix &phfixele, const TrackCol* trackCol,const VertexCol* vtxCol,Double_t _tRho, const ElectronCol* els=0, Bool_t applyElectronVeto=kTRUE);
43 bendavid 1.1 Float_t Ecor() const { return ecor; };
44     Float_t Ecorerr() const { return ecorerr; };
45 bendavid 1.2 Float_t Ecorele() const { return ecorele; };
46     Float_t Ecoreleerr() const { return ecoreleerr; };
47 bendavid 1.1
48     private:
49 bendavid 1.5 UChar_t hasphoton;
50 bendavid 1.1 Float_t e;
51     Float_t pt;
52     Float_t eta;
53     Float_t phi;
54     Float_t r9;
55     Float_t e3x3;
56     Float_t e5x5;
57     Float_t sce;
58     Float_t scrawe;
59     Float_t scpse;
60     Float_t sceta;
61     Float_t scphi;
62 paus 1.9 UInt_t scnclusters;
63     UInt_t scnhits;
64 bendavid 1.2 Float_t scetawidth;
65     Float_t scphiwidth;
66 bendavid 1.1 Float_t hovere;
67     Float_t sigietaieta;
68 paus 1.9 Bool_t isbarrel;
69     Bool_t isr9reco;
70     Bool_t isr9cat;
71     Char_t phcat;
72 bendavid 1.2 Float_t eerr;
73     Float_t eerrsmeared;
74 bendavid 1.3 Float_t esmearing;
75 bendavid 1.5 Float_t idmva;
76     Float_t ecalisodr03;
77     Float_t hcalisodr03;
78     Float_t trkisohollowdr03;
79 bendavid 1.8 Float_t ecalisodr04;
80     Float_t hcalisodr04;
81     Float_t trkisohollowdr04;
82     Float_t trackiso1;
83     Float_t trackiso2;
84     Float_t combiso1;
85     Float_t combiso2;
86 bendavid 1.5
87     //quantities from seed basic cluster
88 bendavid 1.2 Float_t eseed;
89 bendavid 1.4 Float_t etaseed;
90     Float_t phiseed;
91 paus 1.9 Int_t ietaseed;
92     Int_t iphiseed;
93     Int_t ixseed;
94     Int_t iyseed;
95 bendavid 1.4 Float_t etacryseed;
96     Float_t phicryseed;
97     Float_t xcryseed;
98     Float_t ycryseed;
99     Float_t thetaaxisseed;
100     Float_t phiaxisseed;
101 bendavid 1.5 Float_t sigietaietaseed;
102     Float_t sigiphiphiseed;
103     Float_t covietaiphiseed;
104     Float_t e3x3seed;
105     Float_t e5x5seed;
106     Float_t emaxseed;
107     Float_t e2ndseed;
108     Float_t etopseed;
109     Float_t ebottomseed;
110     Float_t eleftseed;
111     Float_t erightseed;
112     Float_t e1x3seed;
113     Float_t e3x1seed;
114     Float_t e1x5seed;
115     Float_t e2x2seed;
116     Float_t e4x4seed;
117     Float_t e2x5maxseed;
118     Float_t e2x5topseed;
119     Float_t e2x5bottomseed;
120     Float_t e2x5leftseed;
121     Float_t e2x5rightseed;
122     Float_t xseedseed;
123     Float_t yseedseed;
124     Float_t zseedseed;
125     UInt_t nhitsseed;
126    
127 bendavid 1.4 //quantities from second basic cluster, if present
128     Float_t ebc2;
129     Float_t etabc2;
130     Float_t phibc2;
131 paus 1.9 Int_t ietabc2;
132     Int_t iphibc2;
133     Int_t ixbc2;
134     Int_t iybc2;
135 bendavid 1.4 Float_t etacrybc2;
136     Float_t phicrybc2;
137     Float_t xcrybc2;
138     Float_t ycrybc2;
139     Float_t thetaaxisbc2;
140 bendavid 1.5 Float_t phiaxisbc2;
141     Float_t sigietaietabc2;
142     Float_t sigiphiphibc2;
143     Float_t covietaiphibc2;
144     Float_t e3x3bc2;
145     Float_t e5x5bc2;
146     Float_t emaxbc2;
147     Float_t e2ndbc2;
148     Float_t etopbc2;
149     Float_t ebottombc2;
150     Float_t eleftbc2;
151     Float_t erightbc2;
152     Float_t e1x3bc2;
153     Float_t e3x1bc2;
154     Float_t e1x5bc2;
155     Float_t e2x2bc2;
156     Float_t e4x4bc2;
157     Float_t e2x5maxbc2;
158     Float_t e2x5topbc2;
159     Float_t e2x5bottombc2;
160     Float_t e2x5leftbc2;
161     Float_t e2x5rightbc2;
162     Float_t xbc2bc2;
163     Float_t ybc2bc2;
164     Float_t zbc2bc2;
165     UInt_t nhitsbc2;
166    
167     //quantities from lowest energy basic cluster if present
168     Float_t ebclast;
169     Float_t etabclast;
170     Float_t phibclast;
171 paus 1.9 Int_t ietabclast;
172     Int_t iphibclast;
173     Int_t ixbclast;
174     Int_t iybclast;
175 bendavid 1.5 Float_t etacrybclast;
176     Float_t phicrybclast;
177     Float_t xcrybclast;
178     Float_t ycrybclast;
179     Float_t thetaaxisbclast;
180     Float_t phiaxisbclast;
181     Float_t sigietaietabclast;
182     Float_t sigiphiphibclast;
183     Float_t covietaiphibclast;
184     Float_t e3x3bclast;
185     Float_t e5x5bclast;
186     UInt_t nhitsbclast;
187    
188     //quantities from second lowest energy basic cluster if present
189     Float_t ebclast2;
190     Float_t etabclast2;
191     Float_t phibclast2;
192 paus 1.9 Int_t ietabclast2;
193     Int_t iphibclast2;
194     Int_t ixbclast2;
195     Int_t iybclast2;
196 bendavid 1.5 Float_t etacrybclast2;
197     Float_t phicrybclast2;
198     Float_t xcrybclast2;
199     Float_t ycrybclast2;
200     Float_t thetaaxisbclast2;
201     Float_t phiaxisbclast2;
202     Float_t sigietaietabclast2;
203     Float_t sigiphiphibclast2;
204     Float_t covietaiphibclast2;
205     Float_t e3x3bclast2;
206     Float_t e5x5bclast2;
207     UInt_t nhitsbclast2;
208 bendavid 1.1
209     //energy correction quantities from PhotonFix
210     Float_t ecor;
211     Float_t ecorerr;
212 bendavid 1.2 Float_t ecorele;
213     Float_t ecoreleerr;
214 bendavid 1.1 Float_t etac;
215     Float_t etas;
216     Float_t etam;
217     Float_t phic;
218     Float_t phis;
219     Float_t phim;
220     Float_t xz;
221     Float_t xc;
222     Float_t xs;
223     Float_t xm;
224     Float_t yz;
225     Float_t yc;
226     Float_t ys;
227     Float_t ym;
228    
229     //conversion quantities
230 bendavid 1.5 UChar_t hasconversion;
231 bendavid 1.1 Float_t convp;
232     Float_t convpt;
233     Float_t conveta;
234     Float_t convphi;
235     Float_t convdeta;
236     Float_t convdphi;
237     Float_t convvtxrho;
238     Float_t convvtxz;
239     Float_t convvtxphi;
240     Float_t convleadpt;
241     Float_t convtrailpt;
242     Float_t convleadtrackpt;
243 paus 1.9 Char_t convleadtrackalgo;
244     Char_t convleadtrackalgos;
245     Char_t convleadtrackcharge;
246 bendavid 1.1 Float_t convtrailtrackpt;
247 paus 1.9 Char_t convtrailtrackalgo;
248     Char_t convtrailtrackalgos;
249     Char_t trailtrackcharge;
250 bendavid 1.1
251 bendavid 1.2 //electron quantities
252 bendavid 1.5 UChar_t haselectron;
253     UChar_t eleisecaldriven;
254     UChar_t eleistrackerdriven;
255 bendavid 1.2 Float_t elee;
256     Float_t elept;
257     Float_t eleeta;
258     Float_t elephi;
259 paus 1.9 Char_t elecharge;
260 bendavid 1.2 Float_t elefbrem;
261     Float_t eledeta;
262     Float_t eledphi;
263     Float_t elep;
264     Float_t elepin;
265     Float_t elepout;
266    
267     //pf supercluster quantities
268 bendavid 1.5 UChar_t haspfsc;
269 bendavid 1.2 Float_t pfsce;
270     Float_t pfscrawe;
271     Float_t pfsceta;
272     Float_t pfscphi;
273 bendavid 1.1
274     //generator level quantities
275 bendavid 1.5 UChar_t ispromptgen;
276 bendavid 1.1 Float_t gene;
277     Float_t genpt;
278     Float_t geneta;
279     Float_t genphi;
280     Float_t genz;
281 bendavid 1.5 Int_t pdgid;
282     Int_t motherpdgid;
283 bendavid 1.1 };
284    
285     class PhotonTreeWriterDiphotonEvent
286     {
287     public:
288 fabstoec 1.12 // ------------ BTAG STUFF -------------------
289 fabstoec 1.15 Float_t btagJet1;
290     Float_t btagJet1Pt;
291     Float_t btagJet1Eta;
292     Float_t btagJet2;
293     Float_t btagJet2Pt;
294     Float_t btagJet2Eta;
295 fabstoec 1.12 // ----------- LEPTON TAG STUFF -------------
296 fabstoec 1.10 Int_t leptonTag;
297 fabstoec 1.11 // ---------- MUON STUFF --------------------
298     Float_t muonPt;
299     Float_t muonEta;
300     Float_t muDR1;
301     Float_t muDR2;
302     Float_t muIso1;
303     Float_t muIso2;
304     Float_t muIso3;
305     Float_t muIso4;
306     Float_t muD0;
307     Float_t muDZ;
308     Int_t muNhits;
309     Float_t muChi2;
310     Int_t muNpixhits;
311     Int_t muNegs;
312     Int_t muNMatch;
313     // ----------- ELECTRON STUFF --------------
314     Float_t elePt;
315     Float_t eleEta;
316     Float_t eleSCEta;
317     Float_t eleIso1;
318     Float_t eleIso2;
319     Float_t eleIso3;
320     Float_t eleIso4;
321     Float_t eleDist;
322     Float_t eleDcot;
323     Float_t eleCoviee;
324     Float_t eleDphiin;
325     Float_t eleDetain;
326     Float_t eleDR1;
327     Float_t eleDR2;
328     Float_t eleMass1;
329     Float_t eleMass2;
330     Int_t eleNinnerHits;
331     // -----------------------------------------
332 bendavid 1.1 Float_t rho;
333     Float_t genHiggspt;
334     Float_t genHiggsZ;
335     Float_t genmass;
336     Float_t gencostheta;
337 bendavid 1.2 Float_t bsX;
338     Float_t bsY;
339     Float_t bsZ;
340 bendavid 1.7 Float_t bsSigmaZ;
341 bendavid 1.2 Float_t vtxX;
342     Float_t vtxY;
343 bendavid 1.1 Float_t vtxZ;
344     Int_t nVtx;
345     Int_t numPU;
346     Int_t numPUminus;
347     Int_t numPUplus;
348     Float_t mass;
349 bendavid 1.2 Float_t masserr;
350     Float_t masserrsmeared;
351 bendavid 1.4 Float_t masserrwrongvtx;
352     Float_t masserrsmearedwrongvtx;
353     Float_t vtxprob;
354 bendavid 1.5 Float_t deltamvtx;
355 bendavid 1.1 Float_t ptgg;
356 bendavid 1.5 Float_t etagg;
357     Float_t phigg;
358 bendavid 1.1 Float_t costheta;
359 bendavid 1.2 Float_t massele;
360     Float_t ptee;
361     Float_t costhetaele;
362     Float_t mt;
363     Float_t cosphimet;
364     Float_t mtele;
365     Float_t cosphimetele;
366 bendavid 1.1 UInt_t evt;
367     UInt_t run;
368     UInt_t lumi;
369     UChar_t evtcat;
370 bendavid 1.2 UInt_t nobj;
371 fabstoec 1.13 Double_t fromZ; //added: Heng 2/14/2012
372     Double_t fromW; //added: Heng
373     Float_t zpt; //added: Heng
374     Float_t allZpt; //added: Heng
375     Float_t zEta; //added: Heng
376     Float_t allZEta; //added: Heng
377     Float_t corrpfmet;
378     Float_t corrpfmetphi;
379     Float_t corrpfmetx;
380     Float_t corrpfmety;
381     Double_t dphiMetgg;
382     Double_t cosdphiMetgg;
383     Double_t dphiPhPh;
384 bendavid 1.2 Float_t pfmet;
385     Float_t pfmetphi;
386     Float_t pfmetx;
387     Float_t pfmety;
388 bendavid 1.5 UChar_t ismc;
389 bendavid 1.1
390     //corrected quantities from PhotonFix corrections
391     Float_t masscor;
392     Float_t masscorerr;
393 bendavid 1.2 Float_t masscorele;
394     Float_t masscoreleerr;
395 bendavid 1.1
396 bendavid 1.5 //jet quantities
397     Float_t jet1pt;
398     Float_t jet1eta;
399     Float_t jet1phi;
400     Float_t jet1mass;
401     Float_t jet2pt;
402     Float_t jet2eta;
403     Float_t jet2phi;
404     Float_t jet2mass;
405     Float_t jetcentralpt;
406     Float_t jetcentraleta;
407     Float_t jetcentralphi;
408     Float_t jetcentralmass;
409     Float_t dijetpt;
410     Float_t dijeteta;
411     Float_t dijetphi;
412     Float_t dijetmass;
413     Float_t jetetaplus;
414     Float_t jetetaminus;
415    
416     Float_t zeppenfeld;
417     Float_t dphidijetgg;
418    
419 bendavid 1.1 PhotonTreeWriterPhoton photons[2];
420     };
421    
422     class PhotonTreeWriter : public BaseMod
423     {
424     public:
425     PhotonTreeWriter(const char *name ="PhotonTreeWriter",
426     const char *title="Selecting PhotonPairs");
427    
428     ~PhotonTreeWriter();
429    
430     // setting all the input Names
431     void SetInputPhotonsName(const char *n){ fPhotonBranchName= n; }
432     void SetPhotonsFromBranch(bool b) { fPhotonsFromBranch = b; }
433     void SetTrackName(const char *n) { fTrackBranchName = n; }
434     void SetElectronName(const char *n) { fElectronName = n; }
435     void SetConversionName(const char *n) { fConversionName = n; }
436     void SetPUDensityName(const char *n) { fPileUpDenName = n; }
437     void SetPVName(const char *n) { fPVName = n; }
438     void SetPVFromBranch(bool b) { fPVFromBranch = b; }
439     void SetMCParticle(const char *n) { fMCParticleName = n; }
440     void SetPUInfoName(const char *n) { fPileUpName = n; }
441     void SetBeamspotName(const char *n) { fBeamspotName = n; }
442     void SetPFCandName(const char *n) { fPFCandName = n; }
443 bendavid 1.2 void SetSuperClusterName(const char *n) { fSuperClusterName = n; }
444 bendavid 1.5 void SetPFJetName(const char *n) { fPFJetName = n; }
445 fabstoec 1.13 void SetGenJetName(const char *n) { fGenJetName = n; }
446     void SetuncorrPFJetName(const char *n) { funcorrPFJetName = n; }
447    
448 bendavid 1.5 void SetPFJetsFromBranch(Bool_t b) { fPFJetsFromBranch = b; }
449     void SetEnableJets(Bool_t b) { fEnableJets = b; }
450 fabstoec 1.10 void SetApplyLeptonTag(Bool_t b) { fApplyLeptonTag = b; }
451 fabstoec 1.12 void SetApplyBTag(Bool_t b) { fApplyBTag = b; }
452 fabstoec 1.14 void SetApplyPFMetCorr(Bool_t b) { fApplyPFMetCorrections = b; }
453 bendavid 1.2 void SetPhFixDataFile(const char *n) { fPhFixDataFile = n; }
454 bendavid 1.1
455     // set basic Cut variables (FOR PRE-SELECTION)
456    
457     // is Data Or Not?
458 paus 1.9 void SetIsData (Bool_t b) { fIsData = b; };
459 bendavid 1.1
460    
461 bendavid 1.6 void SetApplyElectronVeto(Bool_t b) { fApplyElectronVeto = b; }
462 bendavid 1.1
463 paus 1.9 void SetTupleName(const char* c) { fTupleName = c; }
464 bendavid 1.1 void SetGoodElectronsFromBranch(Bool_t b) { fGoodElectronsFromBranch = b; }
465 paus 1.9 void SetGoodElectronName(TString name) { fGoodElectronName = name; }
466     void SetWriteDiphotonTree(Bool_t b) { fWriteDiphotonTree = b; }
467     void SetWriteSingleTree(Bool_t b) { fWriteSingleTree = b; }
468     void SetLoopOnGoodElectrons(Bool_t b) { fLoopOnGoodElectrons = b; }
469     void SetExcludeSinglePrompt(Bool_t b) { fExcludeSinglePrompt = b; }
470     void SetExcludeDoublePrompt(Bool_t b) { fExcludeDoublePrompt = b; }
471 bendavid 1.1
472 fabstoec 1.10 void SetLeptonTagElectronsName(TString name) { fLeptonTagElectronsName = name; }
473     void SetLeptonTagMuonsName (TString name) { fLeptonTagMuonsName = name; }
474    
475 bendavid 1.1 protected:
476     void Process();
477     void SlaveBegin();
478 paus 1.9 // Private auxiliary methods...
479 bendavid 1.1 void FindHiggsPtAndZ(Float_t& pt, Float_t& z, Float_t& mass);
480 paus 1.9 Float_t GetEventCat (PhotonTools::CiCBaseLineCats cat1,
481     PhotonTools::CiCBaseLineCats cat2);
482 bendavid 1.1
483     // Names for the input Collections
484     TString fPhotonBranchName;
485     TString fElectronName;
486     TString fGoodElectronName;
487     TString fConversionName;
488     TString fTrackBranchName;
489     TString fPileUpDenName;
490     TString fPVName;
491     TString fBeamspotName;
492     TString fPFCandName;
493     TString fMCParticleName;
494     TString fPileUpName;
495 bendavid 1.2 TString fSuperClusterName;
496     TString fPFMetName;
497 bendavid 1.5 TString fPFJetName;
498 fabstoec 1.10
499 fabstoec 1.13 TString fGenJetName; //added to do pfmet correction 05/01/2012
500     TString funcorrPFJetName;
501    
502 fabstoec 1.10 TString fLeptonTagElectronsName;
503     TString fLeptonTagMuonsName;
504    
505 bendavid 1.1
506     // is it Data or MC?
507     Bool_t fIsData;
508    
509     // in case there's some PV pre-selection
510     Bool_t fPhotonsFromBranch;
511     Bool_t fPVFromBranch;
512     Bool_t fGoodElectronsFromBranch;
513 bendavid 1.5 Bool_t fPFJetsFromBranch;
514 bendavid 1.1
515 paus 1.9 const PhotonCol *fPhotons;
516     const ElectronCol *fElectrons;
517     const ElectronCol *fGoodElectrons;
518     const DecayParticleCol *fConversions;
519     const TrackCol *fTracks;
520     const PileupEnergyDensityCol *fPileUpDen;
521     const VertexCol *fPV;
522     const BeamSpotCol *fBeamspot;
523     const PFCandidateCol *fPFCands;
524     const MCParticleCol *fMCParticles;
525     const PileupInfoCol *fPileUp;
526     const SuperClusterCol *fSuperClusters;
527     const PFMetCol *fPFMet;
528     const JetCol *fPFJets;
529 fabstoec 1.13 const GenJetCol *fGenJets;
530     const PFJetCol *funcorrPFJets;
531 bendavid 1.2
532 fabstoec 1.10 const ElectronCol *fLeptonTagElectrons;
533     const MuonCol *fLeptonTagMuons;
534    
535 bendavid 1.1 // --------------------------------
536 paus 1.9 Bool_t fLoopOnGoodElectrons; //loop over good elecs instead of photons
537     Bool_t fApplyElectronVeto; //invert elec veto (for cic sel. only atm)
538     Bool_t fWriteDiphotonTree;
539     Bool_t fWriteSingleTree;
540 bendavid 1.1
541 paus 1.9 Bool_t fExcludeSinglePrompt;
542     Bool_t fExcludeDoublePrompt;
543 bendavid 1.2
544 paus 1.9 Bool_t fEnableJets;
545 fabstoec 1.10
546     Bool_t fApplyLeptonTag;
547 fabstoec 1.12 Bool_t fApplyBTag;
548 fabstoec 1.14 Bool_t fApplyPFMetCorrections;
549 fabstoec 1.12
550 paus 1.9 TString fPhFixDataFile;
551     PhotonFix fPhfixph;
552     PhotonFix fPhfixele;
553 bendavid 1.1
554     // --------------------------------
555     // validation Tuple
556 paus 1.9 TString fTupleName;
557 bendavid 1.1 PhotonTreeWriterDiphotonEvent* fDiphotonEvent;
558 paus 1.9 PhotonTreeWriterPhoton* fSinglePhoton;
559     TTree* hCiCTuple;
560     TTree* hCiCTupleSingle;
561 bendavid 1.1
562     ClassDef(PhotonTreeWriter, 1) // Photon identification module
563     };
564     }
565     #endif