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.1 by yilmaz, Tue Sep 20 19:45:05 2011 UTC vs.
Revision 1.10 by ylai, Fri May 11 14:39:31 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 73 | Line 78 | class HiInclusiveJetAnalyzer : public ed
78    bool useCentrality_;
79    bool useVtx_;
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 94 | 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 114 | Line 127 | class HiInclusiveJetAnalyzer : public ed
127      float jteta[MAXJETS];
128      float jtphi[MAXJETS];
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 discr_fr01[MAXJETS];
173 +
174      float refpt[MAXJETS];
175      float refeta[MAXJETS];
176      float refphi[MAXJETS];
# Line 121 | Line 178 | class HiInclusiveJetAnalyzer : public ed
178      float refdphijt[MAXJETS];
179      float refdrjt[MAXJETS];
180      float refparton_pt[MAXJETS];
181 <    float refparton_flavor[MAXJETS];
181 >    int refparton_flavor[MAXJETS];
182 >    int refparton_flavorForB[MAXJETS];
183  
184      float pthat;
185 +    int beamId1, beamId2;
186      int ngen;
187      int genmatchindex[MAXJETS];
188      float genpt[MAXJETS];
# Line 132 | Line 191 | class HiInclusiveJetAnalyzer : public ed
191      float geny[MAXJETS];
192      float gendphijt[MAXJETS];
193      float gendrjt[MAXJETS];
194 +    int gensubid[MAXJETS];
195  
196      // hlt
197      int nHLTBit;
# Line 143 | Line 203 | class HiInclusiveJetAnalyzer : public ed
203      int nL1ABit;
204      bool l1ABit[MAXHLTBITS];
205  
206 +    int bMult;
207 +    int bJetIndex[MAXBFRAG];
208 +    int bStatus[MAXBFRAG];
209 +    int bPdg[MAXBFRAG];
210 +    int bChg[MAXBFRAG];
211 +    float bVx[MAXBFRAG];
212 +    float bVy[MAXBFRAG];
213 +    float bVz[MAXBFRAG];
214 +    float bPt[MAXBFRAG];
215 +    float bEta[MAXBFRAG];
216 +    float bPhi[MAXBFRAG];
217 +
218 +
219    };
220  
221    JRA jets_;

Diff Legend

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