ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/CmsHi/JetAnalysis/interface/HiInclusiveJetAnalyzer.h
Revision: 1.3
Committed: Tue Mar 6 15:03:23 2012 UTC (13 years, 2 months ago) by mnguyen
Content type: text/plain
Branch: MAIN
Changes since 1.2: +36 -5 lines
Log Message: 
Added option to do btagging

File Contents

# User Rev Content
1 yilmaz 1.1 #ifndef MNguyen_HiInclusiveJetAnalyzer_inclusiveJetAnalyzer_
2     #define MNguyen_HiInclusiveJetAnalyzer_inclusiveJetAnalyzer_
3    
4     // system include files
5     #include <memory>
6     #include <string>
7     #include <iostream>
8    
9     // user include files
10     #include "FWCore/Framework/interface/Frameworkfwd.h"
11     #include "FWCore/Framework/interface/EDAnalyzer.h"
12     #include "FWCore/ParameterSet/interface/ParameterSet.h"
13     #include "FWCore/Framework/interface/Event.h"
14     #include "FWCore/Framework/interface/MakerMacros.h"
15    
16     #include "TFile.h"
17     #include "TTree.h"
18     #include "TH1.h"
19     #include "CommonTools/UtilAlgos/interface/TFileService.h"
20     #include "FWCore/ServiceRegistry/interface/Service.h"
21    
22     #include "DataFormats/HeavyIonEvent/interface/CentralityProvider.h"
23    
24 mnguyen 1.3 #include "DataFormats/PatCandidates/interface/Jet.h"
25     #include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
26 yilmaz 1.1 //
27    
28     /**\class HiInclusiveJetAnalyzer
29    
30     \author Matt Nguyen
31     \date November 2010
32     */
33    
34    
35    
36    
37     class HiInclusiveJetAnalyzer : public edm::EDAnalyzer {
38     public:
39    
40     explicit HiInclusiveJetAnalyzer(const edm::ParameterSet&);
41    
42     ~HiInclusiveJetAnalyzer();
43    
44     virtual void analyze(const edm::Event&, const edm::EventSetup&);
45    
46     virtual void beginRun(const edm::Run & r, const edm::EventSetup & c);
47    
48     virtual void beginJob();
49    
50     void fillL1Bits(const edm::Event &iEvent);
51    
52     void fillHLTBits(const edm::Event &iEvent);
53    
54     template <typename TYPE>
55     void getProduct(const std::string name, edm::Handle<TYPE> &prod,
56     const edm::Event &event) const;
57     template <typename TYPE>
58     bool getProductSafe(const std::string name, edm::Handle<TYPE> &prod,
59     const edm::Event &event) const;
60    
61    
62     private:
63    
64 mnguyen 1.3 int getPFJetMuon(const pat::Jet& pfJet, const reco::PFCandidateCollection *pfCandidateColl);
65 yilmaz 1.1
66 mnguyen 1.3 double getPtRel(const reco::PFCandidate lep, const pat::Jet& jet );
67 yilmaz 1.1
68 mnguyen 1.3 edm::InputTag jetTag_, vtxTag_, genjetTag_, eventInfoTag_, L1gtReadout_, pfCandidateLabel_;
69 yilmaz 1.1
70    
71     /// verbose ?
72     bool verbose_;
73    
74     bool useCentrality_;
75     bool useVtx_;
76     bool useJEC_;
77 yilmaz 1.2 bool usePat_;
78 yilmaz 1.1 bool isMC_;
79    
80 mnguyen 1.3 bool doLifeTimeTagging_;
81 yilmaz 1.1
82     TTree *t;
83     edm::Service<TFileService> fs1;
84    
85     CentralityProvider * centrality_;
86    
87    
88    
89     std::string hltResName_; //HLT trigger results name
90     std::vector<std::string> hltProcNames_; //HLT process name(s)
91     std::vector<std::string> hltTrgNames_; //HLT trigger name(s)
92    
93     std::vector<int> hltTrgBits_; //HLT trigger bit(s)
94     std::vector<bool> hltTrgDeci_; //HLT trigger descision(s)
95     std::vector<std::string> hltTrgUsedNames_; //HLT used trigger name(s)
96     std::string hltUsedResName_; //used HLT trigger results name
97    
98    
99    
100     static const int MAXJETS = 50000;
101     static const int MAXHLTBITS = 500000;
102    
103    
104     struct JRA{
105    
106     int nref;
107     int run;
108     int evt;
109     int lumi;
110     int bin;
111     float vx, vy, vz;
112     float b;
113     float hf;
114    
115     float rawpt[MAXJETS];
116     float jtpt[MAXJETS];
117     float jteta[MAXJETS];
118     float jtphi[MAXJETS];
119     float jty[MAXJETS];
120 yilmaz 1.2 float jtpu[MAXJETS];
121    
122 mnguyen 1.3 float discr_csvMva[MAXJETS];
123     float discr_csvSimple[MAXJETS];
124     float discr_muByIp3[MAXJETS];
125     float discr_muByPt[MAXJETS];
126     float discr_prob[MAXJETS];
127     float discr_probb[MAXJETS];
128     float discr_tcHighEff[MAXJETS];
129     float discr_tcHighPur[MAXJETS];
130    
131     unsigned nsvtx[MAXJETS];
132     unsigned svtxntrk[MAXJETS];
133     float svtxdl[MAXJETS];
134     float svtxdls[MAXJETS];
135     float svtxm[MAXJETS];
136     float svtxpt[MAXJETS];
137    
138     unsigned nIPtracks[MAXJETS];
139     unsigned nselIPtracks[MAXJETS];
140    
141     float mue[MAXJETS];
142     float mupt[MAXJETS];
143     float mueta[MAXJETS];
144     float muphi[MAXJETS];
145     float mudr[MAXJETS];
146     float muptrel[MAXJETS];
147     int muchg[MAXJETS];
148    
149 yilmaz 1.1 float refpt[MAXJETS];
150     float refeta[MAXJETS];
151     float refphi[MAXJETS];
152     float refy[MAXJETS];
153     float refdphijt[MAXJETS];
154     float refdrjt[MAXJETS];
155     float refparton_pt[MAXJETS];
156 mnguyen 1.3 int refparton_flavor[MAXJETS];
157     int refparton_flavorForB[MAXJETS];
158 yilmaz 1.1
159     float pthat;
160     int ngen;
161     int genmatchindex[MAXJETS];
162     float genpt[MAXJETS];
163     float geneta[MAXJETS];
164     float genphi[MAXJETS];
165     float geny[MAXJETS];
166     float gendphijt[MAXJETS];
167     float gendrjt[MAXJETS];
168    
169     // hlt
170     int nHLTBit;
171     bool hltBit[MAXHLTBITS];
172    
173     // l1
174     int nL1TBit;
175     bool l1TBit[MAXHLTBITS];
176     int nL1ABit;
177     bool l1ABit[MAXHLTBITS];
178    
179 mnguyen 1.3
180 yilmaz 1.1 };
181    
182     JRA jets_;
183    
184     };
185    
186     #endif