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

Diff Legend

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