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.7 by mnguyen, Thu May 3 21:09:56 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 76 | Line 81 | class HiInclusiveJetAnalyzer : public ed
81    bool usePat_;
82    bool isMC_;
83  
84 +  double genPtMin_;
85 +  bool doLifeTimeTagging_;
86 +  bool doLifeTimeTaggingExtras_;
87 +  bool saveBfragments_;
88  
89    TTree *t;
90    edm::Service<TFileService> fs1;
# Line 95 | Line 104 | class HiInclusiveJetAnalyzer : public ed
104  
105  
106  
107 <  static const int MAXJETS = 50000;
108 <  static const int MAXHLTBITS = 500000;
109 <
107 >  static const int MAXJETS = 500;
108 >  static const int MAXTRACKS = 5000;
109 >  static const int MAXHLTBITS = 5000;
110 >  static const int MAXBFRAG = 500;
111  
112    struct JRA{
113      
# Line 117 | Line 127 | class HiInclusiveJetAnalyzer : public ed
127      float jty[MAXJETS];
128      float jtpu[MAXJETS];
129  
130 +    float discr_csvMva[MAXJETS];
131 +    float discr_csvSimple[MAXJETS];
132 +    float discr_muByIp3[MAXJETS];
133 +    float discr_muByPt[MAXJETS];
134 +    float discr_prob[MAXJETS];
135 +    float discr_probb[MAXJETS];
136 +    float discr_tcHighEff[MAXJETS];    
137 +    float discr_tcHighPur[MAXJETS];
138 +
139 +    int nsvtx[MAXJETS];    
140 +    int svtxntrk[MAXJETS];
141 +    float svtxdl[MAXJETS];  
142 +    float svtxdls[MAXJETS];  
143 +    float svtxm[MAXJETS];    
144 +    float svtxpt[MAXJETS];  
145 +
146 +    int nIPtrk[MAXJETS];
147 +    int nselIPtrk[MAXJETS];
148 +    
149 +    int nIP;
150 +    int ipJetIndex[MAXTRACKS];
151 +    float ipPt[MAXTRACKS];
152 +    float ipProb0[MAXTRACKS];
153 +    float ipProb1[MAXTRACKS];
154 +    float ip2d[MAXTRACKS];
155 +    float ip2dSig[MAXTRACKS];
156 +    float ip3d[MAXTRACKS];
157 +    float ip3dSig[MAXTRACKS];
158 +    float ipDist2Jet[MAXTRACKS];
159 +    float ipDist2JetSig[MAXTRACKS];
160 +    float ipClosest2Jet[MAXTRACKS];
161 +
162 +    float mue[MAXJETS];    
163 +    float mupt[MAXJETS];    
164 +    float mueta[MAXJETS];  
165 +    float muphi[MAXJETS];  
166 +    float mudr[MAXJETS];    
167 +    float muptrel[MAXJETS];
168 +    int muchg[MAXJETS];  
169 +    
170      float refpt[MAXJETS];
171      float refeta[MAXJETS];
172      float refphi[MAXJETS];
# Line 124 | Line 174 | class HiInclusiveJetAnalyzer : public ed
174      float refdphijt[MAXJETS];
175      float refdrjt[MAXJETS];
176      float refparton_pt[MAXJETS];
177 <    float refparton_flavor[MAXJETS];
177 >    int refparton_flavor[MAXJETS];
178 >    int refparton_flavorForB[MAXJETS];
179  
180      float pthat;
181 +    int beamId1, beamId2;
182      int ngen;
183      int genmatchindex[MAXJETS];
184      float genpt[MAXJETS];
# Line 146 | Line 198 | class HiInclusiveJetAnalyzer : public ed
198      int nL1ABit;
199      bool l1ABit[MAXHLTBITS];
200  
201 +    int bMult;
202 +    int bJetIndex[MAXBFRAG];
203 +    int bStatus[MAXBFRAG];
204 +    int bPdg[MAXBFRAG];
205 +    int bChg[MAXBFRAG];
206 +    float bVx[MAXBFRAG];
207 +    float bVy[MAXBFRAG];
208 +    float bVz[MAXBFRAG];
209 +    float bPt[MAXBFRAG];
210 +    float bEta[MAXBFRAG];
211 +    float bPhi[MAXBFRAG];
212 +
213 +
214    };
215  
216    JRA jets_;

Diff Legend

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