ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/interface/HiInclusiveJetAnalyzer.h
(Generate patch)

Comparing UserCode/CmsHi/JetAnalysis/interface/HiInclusiveJetAnalyzer.h (file contents):
Revision 1.2 by yilmaz, Tue Sep 20 20:06:06 2011 UTC vs.
Revision 1.9 by yilmaz, Thu May 10 16:08:29 2012 UTC

# Line 10 | Line 10
10   #include "FWCore/Framework/interface/Frameworkfwd.h"
11   #include "FWCore/Framework/interface/EDAnalyzer.h"
12   #include "FWCore/ParameterSet/interface/ParameterSet.h"
13 –
13   #include "FWCore/Framework/interface/Event.h"
14   #include "FWCore/Framework/interface/MakerMacros.h"
15  
17 –
18 –
16   #include "TFile.h"
17   #include "TTree.h"
18   #include "TH1.h"
# Line 24 | Line 21
21  
22   #include "DataFormats/HeavyIonEvent/interface/CentralityProvider.h"
23  
24 + #include "DataFormats/PatCandidates/interface/Jet.h"
25 + #include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
26   //
27  
28   /**\class HiInclusiveJetAnalyzer
# Line 62 | Line 61 | class HiInclusiveJetAnalyzer : public ed
61  
62   private:
63    
64 +  int getPFJetMuon(const pat::Jet& pfJet, const reco::PFCandidateCollection *pfCandidateColl);
65  
66 +  double getPtRel(const reco::PFCandidate lep, const pat::Jet& jet );
67  
68 <  edm::InputTag   jetTag_, vtxTag_, genjetTag_, eventInfoTag_, L1gtReadout_;
68 >  void saveDaughters( const reco::GenParticle & gen);
69 >  void saveDaughters( const reco::Candidate & gen);
70  
71 +  edm::InputTag   jetTag_, vtxTag_, genjetTag_, eventInfoTag_, L1gtReadout_, pfCandidateLabel_;
72 +
73 +  std::vector<float> usedStringPts;
74  
75    /// verbose ?
76    bool   verbose_;
# Line 75 | Line 80 | class HiInclusiveJetAnalyzer : public ed
80    bool useJEC_;
81    bool usePat_;
82    bool isMC_;
83 +  bool doTrigger_;
84  
85 +  bool doSubEvent_;
86 +  double genPtMin_;
87 +  bool doLifeTimeTagging_;
88 +  bool doLifeTimeTaggingExtras_;
89 +  bool saveBfragments_;
90  
91    TTree *t;
92    edm::Service<TFileService> fs1;
# Line 95 | Line 106 | class HiInclusiveJetAnalyzer : public ed
106  
107  
108  
109 <  static const int MAXJETS = 50000;
110 <  static const int MAXHLTBITS = 500000;
111 <
109 >  static const int MAXJETS = 500;
110 >  static const int MAXTRACKS = 5000;
111 >  static const int MAXHLTBITS = 5000;
112 >  static const int MAXBFRAG = 500;
113  
114    struct JRA{
115      
# Line 117 | Line 129 | class HiInclusiveJetAnalyzer : public ed
129      float jty[MAXJETS];
130      float jtpu[MAXJETS];
131  
132 +    float discr_csvMva[MAXJETS];
133 +    float discr_csvSimple[MAXJETS];
134 +    float discr_muByIp3[MAXJETS];
135 +    float discr_muByPt[MAXJETS];
136 +    float discr_prob[MAXJETS];
137 +    float discr_probb[MAXJETS];
138 +    float discr_tcHighEff[MAXJETS];    
139 +    float discr_tcHighPur[MAXJETS];
140 +
141 +    int nsvtx[MAXJETS];    
142 +    int svtxntrk[MAXJETS];
143 +    float svtxdl[MAXJETS];  
144 +    float svtxdls[MAXJETS];  
145 +    float svtxm[MAXJETS];    
146 +    float svtxpt[MAXJETS];  
147 +
148 +    int nIPtrk[MAXJETS];
149 +    int nselIPtrk[MAXJETS];
150 +    
151 +    int nIP;
152 +    int ipJetIndex[MAXTRACKS];
153 +    float ipPt[MAXTRACKS];
154 +    float ipProb0[MAXTRACKS];
155 +    float ipProb1[MAXTRACKS];
156 +    float ip2d[MAXTRACKS];
157 +    float ip2dSig[MAXTRACKS];
158 +    float ip3d[MAXTRACKS];
159 +    float ip3dSig[MAXTRACKS];
160 +    float ipDist2Jet[MAXTRACKS];
161 +    float ipDist2JetSig[MAXTRACKS];
162 +    float ipClosest2Jet[MAXTRACKS];
163 +
164 +    float mue[MAXJETS];    
165 +    float mupt[MAXJETS];    
166 +    float mueta[MAXJETS];  
167 +    float muphi[MAXJETS];  
168 +    float mudr[MAXJETS];    
169 +    float muptrel[MAXJETS];
170 +    int muchg[MAXJETS];  
171 +    
172      float refpt[MAXJETS];
173      float refeta[MAXJETS];
174      float refphi[MAXJETS];
# Line 124 | Line 176 | class HiInclusiveJetAnalyzer : public ed
176      float refdphijt[MAXJETS];
177      float refdrjt[MAXJETS];
178      float refparton_pt[MAXJETS];
179 <    float refparton_flavor[MAXJETS];
179 >    int refparton_flavor[MAXJETS];
180 >    int refparton_flavorForB[MAXJETS];
181  
182      float pthat;
183 +    int beamId1, beamId2;
184      int ngen;
185      int genmatchindex[MAXJETS];
186      float genpt[MAXJETS];
# Line 135 | Line 189 | class HiInclusiveJetAnalyzer : public ed
189      float geny[MAXJETS];
190      float gendphijt[MAXJETS];
191      float gendrjt[MAXJETS];
192 +    int gensubid[MAXJETS];
193  
194      // hlt
195      int nHLTBit;
# Line 146 | Line 201 | class HiInclusiveJetAnalyzer : public ed
201      int nL1ABit;
202      bool l1ABit[MAXHLTBITS];
203  
204 +    int bMult;
205 +    int bJetIndex[MAXBFRAG];
206 +    int bStatus[MAXBFRAG];
207 +    int bPdg[MAXBFRAG];
208 +    int bChg[MAXBFRAG];
209 +    float bVx[MAXBFRAG];
210 +    float bVy[MAXBFRAG];
211 +    float bVz[MAXBFRAG];
212 +    float bPt[MAXBFRAG];
213 +    float bEta[MAXBFRAG];
214 +    float bPhi[MAXBFRAG];
215 +
216 +
217    };
218  
219    JRA jets_;

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines