ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CITCommon/CommonData/interface/ElectronTree.h
(Generate patch)

Comparing UserCode/CITCommon/CommonData/interface/ElectronTree.h (file contents):
Revision 1.2 by sixie, Sat Jul 21 06:57:23 2012 UTC vs.
Revision 1.6 by sixie, Wed Sep 19 09:18:21 2012 UTC

# Line 24 | Line 24 | namespace citana
24                                   kEleTrigger_Ele_CaloIdT_TrkIdVL_CaloIsoVL_TrkIsoVL = 0x000004
25        };
26  
27 +      //*******************************************
28 +      //=== tree versions ===
29 +      //*******************************************
30 +      enum ElectronTreeVersion { kTreeV1, kTreeV2 };
31 +
32        /// variables
33        Float_t                 fWeight;
34        UInt_t                  fRunNumber;
# Line 40 | Line 45 | namespace citana
45        Bool_t                  fEleIsEcalDriven;
46        UInt_t                  fEleTriggerBit;
47        Float_t                 fRho;
48 <      UInt_t                  fNVertices;
48 >      Float_t                 fNVertices;
49  
50        // Conversion and IP
51        Float_t                 fEleD0;
# Line 50 | Line 55 | namespace citana
55        Bool_t                  fEleMatchedConversion;
56        Float_t                 fEleConvDCot;
57        Float_t                 fEleConvDist;
58 <      UInt_t                  fEleNMissHits;
58 >      Float_t                 fEleNMissHits;
59  
60        // E/P variables
61        Float_t                 fEleNBrem;
# Line 80 | Line 85 | namespace citana
85        //track quality
86        Float_t                 fEleGsfTrackChi2OverNdof;    
87        Float_t                 fEleKFTrackChi2OverNDoF;
88 <      Int_t                   fEleKFTrackNHits;
89 <      Int_t                   fEleKFTrackNLayersWithMeasurement;
88 >      Float_t                 fEleKFTrackNHits;
89 >      Float_t                 fEleKFTrackNLayersWithMeasurement;
90        Float_t                 fEleOneMinusSeedE1x5OverE5x5;
91  
92        //Isolation Variables
# Line 111 | Line 116 | namespace citana
116  
117        Bool_t                  fElePassTriggerDenominator;
118  
119 +      //Regression Variables
120 +      Bool_t                  fIsEB;          
121 +      Bool_t                  fIsEE;          
122 +      Float_t                 fSCRawEnergy;
123 +      Float_t                 fNClusters;
124 +      Float_t                 fEtaSeed;
125 +      Float_t                 fPhiSeed;
126 +      Float_t                 fESeed;
127 +      Float_t                 fE3x3Seed;
128 +      Float_t                 fE5x5Seed;
129 +      Float_t                 fEMaxSeed;
130 +      Float_t                 fE2ndSeed;
131 +      Float_t                 fETopSeed;
132 +      Float_t                 fEBottomSeed;
133 +      Float_t                 fELeftSeed;
134 +      Float_t                 fERightSeed;
135 +      Float_t                 fE2x5MaxSeed;
136 +      Float_t                 fE2x5TopSeed;
137 +      Float_t                 fE2x5BottomSeed;
138 +      Float_t                 fE2x5LeftSeed;
139 +      Float_t                 fE2x5RightSeed;
140 +      Float_t                 fIEtaSeed;
141 +      Float_t                 fIPhiSeed;
142 +      Float_t                 fEtaCrySeed;
143 +      Float_t                 fPhiCrySeed;
144 +      Float_t                 fEcalEnergyError;
145 +      Float_t                 fGsfTrackPIn;
146 +      Float_t                 fTrackMomentumError;
147 +      Float_t                 fpmeangsf;
148 +      Float_t                 fpmeankf;
149 +      Float_t                 fCharge;
150 +      Float_t                 fGeneratedEnergy;
151 +      Float_t                 fGeneratedEnergyStatus1;
152 +      Float_t                 fGeneratedEnergyStatus3;
153 +
154      public:
155        /// this is the main element
156        TTree *tree_;
# Line 134 | Line 174 | namespace citana
174          fEventNumber                   = 0.0;
175          fEleEventNumberParity          = 0.0;
176          fElePt                         = 0.0;
177 <        fEleEta                                = 0.0;
178 <        fElePhi                                = 0.0;
179 <        fEleSCEt                               = 0.0;
177 >        fEleEta                        = 0.0;
178 >        fElePhi                        = 0.0;
179 >        fEleSCEt                       = 0.0;
180          fEleSCEta                      = 0.0;
181          fEleSCPhi                      = 0.0;
182          fEleEcalEnergy                 = 0.0;
183 <        fEleIsEcalDriven                       = 0.0;
183 >        fEleIsEcalDriven               = 0.0;
184          fEleTriggerBit                 = 0.0;
185          fRho                           = 0.0;
186          fNVertices                     = 0.0;
187          fEleD0                         = 0.0;
188          fEleDZ                         = 0.0;
189 <        fEleIP3d                               = 0.0;
189 >        fEleIP3d                       = 0.0;
190          fEleIP3dSig                    = 0.0;
191          fEleMatchedConversion          = 0.0;
192          fEleConvDCot                   = 0.0;
# Line 155 | Line 195 | namespace citana
195          fEleNBrem                      = 0.0;
196          fEleFBrem                      = 0.0;
197          fEleEOverP                     = 0.0;
198 <        fEleESeedClusterOverPIn                = 0.0;
199 <        fEleESeedClusterOverPout               = 0.0;
200 <        fEleEEleClusterOverPout                = 0.0;
201 <        fEleOneOverEMinusOneOverP        = 0.0;
198 >        fEleESeedClusterOverPIn        = 0.0;
199 >        fEleESeedClusterOverPout       = 0.0;
200 >        fEleEEleClusterOverPout        = 0.0;
201 >        fEleOneOverEMinusOneOverP      = 0.0;
202          fEleDEtaIn                     = 0.0;
203          fEleDPhiIn                     = 0.0;
204          fEledEtaCalo                   = 0.0;
205          fEledPhiCalo                   = 0.0;
206          fEleSigmaIEtaIEta              = 0.0;
207          fEleSigmaIPhiIPhi              = 0.0;
208 <        fEleSigmaIEtaIPhi                      = 0.0;
208 >        fEleSigmaIEtaIPhi              = 0.0;
209          fEleSCEtaWidth                 = 0.0;
210          fEleSCPhiWidth                 = 0.0;
211          fEleR9                         = 0.0;
212          fElePreShowerOverRaw           = 0.0;
213          fEleHoverE                     = 0.0;
214 <        fEleGsfTrackChi2OverNdof         = 0.0;
214 >        fEleGsfTrackChi2OverNdof       = 0.0;
215          fEleKFTrackChi2OverNDoF        = 0.0;
216 <        fEleKFTrackNHits                       = 0.0;
216 >        fEleKFTrackNHits               = 0.0;
217          fEleKFTrackNLayersWithMeasurement = 0.0;
218 <        fEleOneMinusSeedE1x5OverE5x5     = 0.0;
219 <        fElePFMVA                              = 0.0;
218 >        fEleOneMinusSeedE1x5OverE5x5   = 0.0;
219 >        fElePFMVA                      = 0.0;
220          fEleTrkIso03                   = 0.0;
221          fEleEMIso03                    = 0.0;
222          fEleHadIso03                   = 0.0;
223          fEleTrkIso04                   = 0.0;
224          fEleEMIso04                    = 0.0;
225 <        fEleHadIso04                           = 0.0;
225 >        fEleHadIso04                   = 0.0;
226          fElePFIso04                    = 0.0;
227          fChargedIso_DR0p0To0p1         = 0.0;
228          fChargedIso_DR0p1To0p2         = 0.0;
# Line 194 | Line 234 | namespace citana
234          fGammaIso_DR0p2To0p3           = 0.0;
235          fGammaIso_DR0p3To0p4           = 0.0;
236          fGammaIso_DR0p4To0p5           = 0.0;
237 <        fNeutralHadronIso_DR0p0To0p1     = 0.0;
238 <        fNeutralHadronIso_DR0p1To0p2     = 0.0;
239 <        fNeutralHadronIso_DR0p2To0p3     = 0.0;
240 <        fNeutralHadronIso_DR0p3To0p4     = 0.0;
241 <        fNeutralHadronIso_DR0p4To0p5     = 0.0;
242 <        fElePassTriggerDenominator       = 0.0;
237 >        fNeutralHadronIso_DR0p0To0p1   = 0.0;
238 >        fNeutralHadronIso_DR0p1To0p2   = 0.0;
239 >        fNeutralHadronIso_DR0p2To0p3   = 0.0;
240 >        fNeutralHadronIso_DR0p3To0p4   = 0.0;
241 >        fNeutralHadronIso_DR0p4To0p5   = 0.0;
242 >        fElePassTriggerDenominator     = 0.0;
243 >        fIsEB                          = 0.0;          
244 >        fIsEE                          = 0.0;          
245 >        fSCRawEnergy                   = 0.0;
246 >        fNClusters                     = 0.0;
247 >        fEtaSeed                       = 0.0;
248 >        fPhiSeed                       = 0.0;
249 >        fESeed                         = 0.0;
250 >        fE3x3Seed                      = 0.0;
251 >        fE5x5Seed                      = 0.0;
252 >        fEMaxSeed                      = 0.0;
253 >        fE2ndSeed                      = 0.0;
254 >        fETopSeed                      = 0.0;
255 >        fEBottomSeed                   = 0.0;
256 >        fELeftSeed                     = 0.0;
257 >        fERightSeed                    = 0.0;
258 >        fE2x5MaxSeed                   = 0.0;
259 >        fE2x5TopSeed                   = 0.0;
260 >        fE2x5BottomSeed                = 0.0;
261 >        fE2x5LeftSeed                  = 0.0;
262 >        fE2x5RightSeed                 = 0.0;
263 >        fIEtaSeed                      = 0.0;
264 >        fIPhiSeed                      = 0.0;
265 >        fEtaCrySeed                    = 0.0;
266 >        fPhiCrySeed                    = 0.0;
267 >        fEcalEnergyError               = 0.0;
268 >        fGsfTrackPIn                   = 0.0;
269 >        fTrackMomentumError            = 0.0;
270 >        fpmeangsf                      = 0.0;
271 >        fpmeankf                       = 0.0;
272 >        fCharge                        = 0.0;
273 >        fGeneratedEnergy               = 0.0;
274 >        fGeneratedEnergyStatus1        = 0.0;
275 >        fGeneratedEnergyStatus3        = 0.0;
276 >
277 >
278        }
279      
280        /// load a ElectronTree
281 <      void LoadTree(const char* file){
281 >      void LoadTree(const char* file, int version=kTreeV1){
282          f_ = TFile::Open(file);
283          assert(f_);
284 <        tree_ = dynamic_cast<TTree*>(f_->Get("Electrons"));
284 >        if(version==kTreeV1) tree_ = dynamic_cast<TTree*>(f_->Get("Electrons"));
285 >        else if(version==kTreeV2) tree_ = dynamic_cast<TTree*>(f_->Get("eleIDdir/T1"));
286          assert(tree_);
287        }
288      
# Line 231 | Line 307 | namespace citana
307          tree_->Branch("ecaldriven",&fEleIsEcalDriven,"ecaldriven/O");
308          tree_->Branch("triggerBit",&fEleTriggerBit,"triggerBit/i");
309          tree_->Branch("rho",&fRho,"rho/F");
310 <        tree_->Branch("vertices",&fNVertices,"vertices/i");
310 >        tree_->Branch("vertices",&fNVertices,"vertices/F");
311          tree_->Branch("d0",&fEleD0,"d0/F");
312          tree_->Branch("dz",&fEleDZ,"dz/F");
313          tree_->Branch("ip3d",&fEleIP3d,"ip3d/F");
# Line 239 | Line 315 | namespace citana
315          tree_->Branch("matchConv",&fEleMatchedConversion,"matchConv/O");
316          tree_->Branch("dcot",&fEleConvDCot,"dcot/F");
317          tree_->Branch("dist",&fEleConvDist,"dist/F");
318 <        tree_->Branch("missHits",&fEleNMissHits,"NMissHits/i");
318 >        tree_->Branch("missHits",&fEleNMissHits,"NMissHits/F");
319          tree_->Branch("nbrems",&fEleNBrem,"nbrems/F");
320          tree_->Branch("fbrem",&fEleFBrem,"fbrem/F");
321          tree_->Branch("EoP",&fEleEOverP,"EoP/F");
322          tree_->Branch("EoPin",&fEleESeedClusterOverPIn,"EoPin/F");
323          tree_->Branch("ESeedoPout",&fEleESeedClusterOverPout,"ESeedoPout/F");
324          tree_->Branch("EEleoPout",&fEleEEleClusterOverPout,"EEleoPout/F");
325 < //       tree_->Branch("eleEoPout",&fEleEEleClusterOverPout,"eleEoPout/F");
326 <        tree_->Branch("IoEmIoP",&fEleOneOverEMinusOneOverP,"IoEmIoP/F");
325 > //       tree_->Branch("eleEoPout",&fEleEEleClusterOverPout,"eleEoPout/F"); //old name of the branch
326 >        tree_->Branch("IoEmIoP",&fEleOneOverEMinusOneOverP,"IoEmIoP/F");  
327          tree_->Branch("deta",&fEleDEtaIn,"deta/F");
328          tree_->Branch("dphi",&fEleDPhiIn,"dphi/F");
329          tree_->Branch("detacalo",&fEledEtaCalo,"detacalo/F");
# Line 262 | Line 338 | namespace citana
338          tree_->Branch("HoE",&fEleHoverE,"HoE/F");
339          tree_->Branch("gsfchi2",&fEleGsfTrackChi2OverNdof,"gsfchi2/F");
340          tree_->Branch("kfchi2",&fEleKFTrackChi2OverNDoF,"kfchi2/F");
341 <        tree_->Branch("kfhits",&fEleKFTrackNHits,"kfhits/I");
342 <        tree_->Branch("kflayers",&fEleKFTrackNLayersWithMeasurement,"kflayers/I");
341 >        tree_->Branch("kfhits",&fEleKFTrackNHits,"kfhits/F");
342 >        tree_->Branch("kflayers",&fEleKFTrackNLayersWithMeasurement,"kflayers/F");
343          tree_->Branch("OneMinusSeedE1x5OverE5x5",&fEleOneMinusSeedE1x5OverE5x5,"OneMinusSeedE1x5OverE5x5/F");
344          tree_->Branch("PFMVA",&fElePFMVA,"PFMVA/F");
345          tree_->Branch("trkIso03",&fEleTrkIso03,"trkIso03/F");
# Line 289 | Line 365 | namespace citana
365          tree_->Branch("NeutralHadronIso_DR0p3To0p4",&fNeutralHadronIso_DR0p3To0p4,"NeutralHadronIso_DR0p3To0p4/F");
366          tree_->Branch("NeutralHadronIso_DR0p4To0p5",&fNeutralHadronIso_DR0p4To0p5,"NeutralHadronIso_DR0p4To0p5/F");
367          tree_->Branch("PassTriggerDenominator",&fElePassTriggerDenominator,"PassTriggerDenominator/O");
368 +        tree_->Branch("IsEB",&fIsEB,"IsEB/O");
369 +        tree_->Branch("IsEE",&fIsEE,"IsEE/O");
370 +        tree_->Branch("SCRawEnergy",&fSCRawEnergy,"SCRawEnergy/F");
371 +        tree_->Branch("NClusters",&fNClusters,"NClusters/F");
372 +        tree_->Branch("EtaSeed",&fEtaSeed,"EtaSeed/F");
373 +        tree_->Branch("PhiSeed",&fPhiSeed,"PhiSeed/F");
374 +        tree_->Branch("ESeed",&fESeed,"ESeed/F");
375 +        tree_->Branch("E3x3Seed",&fE3x3Seed,"E3x3Seed/F");
376 +        tree_->Branch("E5x5Seed",&fE5x5Seed,"E5x5Seed/F");
377 +        tree_->Branch("EMaxSeed",&fEMaxSeed,"EMaxSeed/F");
378 +        tree_->Branch("E2ndSeed",&fE2ndSeed,"E2ndSeed/F");
379 +        tree_->Branch("ETopSeed",&fETopSeed,"ETopSeed/F");
380 +        tree_->Branch("EBottomSeed",&fEBottomSeed,"EBottomSeed/F");
381 +        tree_->Branch("ELeftSeed",&fELeftSeed,"ELeftSeed/F");
382 +        tree_->Branch("ERightSeed",&fERightSeed,"ERightSeed/F");
383 +        tree_->Branch("E2x5MaxSeed",&fE2x5MaxSeed,"E2x5MaxSeed/F");
384 +        tree_->Branch("E2x5TopSeed",&fE2x5TopSeed,"E2x5TopSeed/F");
385 +        tree_->Branch("E2x5BottomSeed",&fE2x5BottomSeed,"E2x5BottomSeed/F");
386 +        tree_->Branch("E2x5LeftSeed",&fE2x5LeftSeed,"E2x5LeftSeed/F");
387 +        tree_->Branch("E2x5RightSeed",&fE2x5RightSeed,"E2x5RightSeed/F");
388 +        tree_->Branch("IEtaSeed",&fIEtaSeed,"IEtaSeed/F");
389 +        tree_->Branch("IPhiSeed",&fIPhiSeed,"IPhiSeed/F");
390 +        tree_->Branch("EtaCrySeed",&fEtaCrySeed,"EtaCrySeed/F");
391 +        tree_->Branch("PhiCrySeed",&fPhiCrySeed,"PhiCrySeed/F");
392 +        tree_->Branch("ecalenergyerror",&fEcalEnergyError,"ecalenergyerror/F");
393 +        tree_->Branch("pmodegsf",&fGsfTrackPIn,"pmodegsf/F");
394 +        tree_->Branch("perror",&fTrackMomentumError,"perror/F");
395 +        tree_->Branch("pmeangsf",&fpmeangsf,"pmodegsf/F");
396 +        tree_->Branch("pmeankf",&fpmeankf,"pmeankf/F");
397 +        tree_->Branch("Charge",&fCharge,"Charge/F");
398 +        tree_->Branch("GeneratedEnergy",&fGeneratedEnergy,"GeneratedEnergy/F");
399 +        tree_->Branch("GeneratedEnergyStatus1",&fGeneratedEnergyStatus1,"GeneratedEnergyStatus1/F");
400 +        tree_->Branch("GeneratedEnergyStatus3",&fGeneratedEnergyStatus3,"GeneratedEnergyStatus3/F");
401  
402  
403        }
# Line 374 | Line 483 | namespace citana
483          tree_->SetBranchAddress("NeutralHadronIso_DR0p3To0p4",&fNeutralHadronIso_DR0p3To0p4);
484          tree_->SetBranchAddress("NeutralHadronIso_DR0p4To0p5",&fNeutralHadronIso_DR0p4To0p5);
485          tree_->SetBranchAddress("PassTriggerDenominator",&fElePassTriggerDenominator);
486 +        tree_->SetBranchAddress("IsEB",&fIsEB);
487 +        tree_->SetBranchAddress("IsEE",&fIsEE);
488 +        tree_->SetBranchAddress("SCRawEnergy",&fSCRawEnergy);
489 +        tree_->SetBranchAddress("NClusters",&fNClusters);
490 +        tree_->SetBranchAddress("EtaSeed",&fEtaSeed);
491 +        tree_->SetBranchAddress("PhiSeed",&fPhiSeed);
492 +        tree_->SetBranchAddress("ESeed",&fESeed);
493 +        tree_->SetBranchAddress("E3x3Seed",&fE3x3Seed);
494 +        tree_->SetBranchAddress("E5x5Seed",&fE5x5Seed);
495 +        tree_->SetBranchAddress("EMaxSeed",&fEMaxSeed);
496 +        tree_->SetBranchAddress("E2ndSeed",&fE2ndSeed);
497 +        tree_->SetBranchAddress("ETopSeed",&fETopSeed);
498 +        tree_->SetBranchAddress("EBottomSeed",&fEBottomSeed);
499 +        tree_->SetBranchAddress("ELeftSeed",&fELeftSeed);
500 +        tree_->SetBranchAddress("ERightSeed",&fERightSeed);
501 +        tree_->SetBranchAddress("E2x5MaxSeed",&fE2x5MaxSeed);
502 +        tree_->SetBranchAddress("E2x5TopSeed",&fE2x5TopSeed);
503 +        tree_->SetBranchAddress("E2x5BottomSeed",&fE2x5BottomSeed);
504 +        tree_->SetBranchAddress("E2x5LeftSeed",&fE2x5LeftSeed);
505 +        tree_->SetBranchAddress("E2x5RightSeed",&fE2x5RightSeed);
506 +        tree_->SetBranchAddress("IEtaSeed",&fIEtaSeed);
507 +        tree_->SetBranchAddress("IPhiSeed",&fIPhiSeed);
508 +        tree_->SetBranchAddress("EtaCrySeed",&fEtaCrySeed);
509 +        tree_->SetBranchAddress("PhiCrySeed",&fPhiCrySeed);
510 +        tree_->SetBranchAddress("ecalenergyerror",&fEcalEnergyError);
511 +        tree_->SetBranchAddress("pmodegsf",&fGsfTrackPIn);
512 +        tree_->SetBranchAddress("perror",&fTrackMomentumError);
513 +        tree_->SetBranchAddress("pmeangsf",&fpmeankf);
514 +        tree_->SetBranchAddress("pmeankf",&fpmeankf);
515 +        tree_->SetBranchAddress("Charge",&fCharge);
516 +        tree_->SetBranchAddress("GeneratedEnergy",&fGeneratedEnergy);
517 +        tree_->SetBranchAddress("GeneratedEnergyStatus1",&fGeneratedEnergyStatus1);
518 +        tree_->SetBranchAddress("GeneratedEnergyStatus3",&fGeneratedEnergyStatus3);
519  
520          gErrorIgnoreLevel = currentState;
521        }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines