CMS 3D CMS Logo
JetTester_HeavyIons.cc
Go to the documentation of this file.
1 // Producer for validation histograms for Calo and PF background subtracted objects
2 // Modified by Raghav Kunnawalkam Elayavalli, Aug 18th 2014
3 // , Oct 22nd 2014 to run in 73X
4 
5 #include "JetTester_HeavyIons.h"
6 
7 using namespace edm;
8 using namespace reco;
9 using namespace std;
10 
11 JetTester_HeavyIons::JetTester_HeavyIons(const edm::ParameterSet& iConfig) :
12  mInputCollection (iConfig.getParameter<edm::InputTag> ("src")),
13  mInputGenCollection (iConfig.getParameter<edm::InputTag> ("srcGen")),
14  mInputPFCandCollection (iConfig.getParameter<edm::InputTag> ("PFcands")),
15  // mInputCandCollection (iConfig.getParameter<edm::InputTag> ("Cands")),
16  // rhoTag (iConfig.getParameter<edm::InputTag> ("srcRho")),
17  mOutputFile (iConfig.getUntrackedParameter<std::string>("OutputFile","")),
18  JetType (iConfig.getUntrackedParameter<std::string>("JetType")),
19  UEAlgo (iConfig.getUntrackedParameter<std::string>("UEAlgo")),
20  Background (iConfig.getParameter<edm::InputTag> ("Background")),
21  mRecoJetPtThreshold (iConfig.getParameter<double> ("recoJetPtThreshold")),
22  mMatchGenPtThreshold (iConfig.getParameter<double> ("matchGenPtThreshold")),
23  mGenEnergyFractionThreshold (iConfig.getParameter<double> ("genEnergyFractionThreshold")),
24  mReverseEnergyFractionThreshold(iConfig.getParameter<double> ("reverseEnergyFractionThreshold")),
25  mRThreshold (iConfig.getParameter<double> ("RThreshold")),
26  JetCorrectionService (iConfig.getParameter<std::string> ("JetCorrections"))
27 {
28  std::string inputCollectionLabel(mInputCollection.label());
29 
30  isCaloJet = (std::string("calo")==JetType);
31  isJPTJet = (std::string("jpt") ==JetType);
32  isPFJet = (std::string("pf") ==JetType);
33 
34  //consumes
35  pvToken_ = consumes<std::vector<reco::Vertex> >(edm::InputTag("offlinePrimaryVertices"));
36  caloTowersToken_ = consumes<CaloTowerCollection>(edm::InputTag("towerMaker"));
37  if (isCaloJet) caloJetsToken_ = consumes<reco::CaloJetCollection>(mInputCollection);
38  if (isJPTJet) jptJetsToken_ = consumes<reco::JPTJetCollection>(mInputCollection);
39  if (isPFJet) {
40  if(std::string("Pu")==UEAlgo) basicJetsToken_ = consumes<reco::BasicJetCollection>(mInputCollection);
41  }
42 
43  genJetsToken_ = consumes<reco::GenJetCollection>(edm::InputTag(mInputGenCollection));
44  evtToken_ = consumes<GenEventInfoProduct>(edm::InputTag("generator"));
45  pfCandToken_ = consumes<reco::PFCandidateCollection>(mInputPFCandCollection);
46  pfCandViewToken_ = consumes<reco::CandidateView>(mInputPFCandCollection);
47  caloCandViewToken_ = consumes<reco::CandidateView>(edm::InputTag("towerMaker"));
48  backgrounds_ = consumes<edm::ValueMap<reco::VoronoiBackground> >(Background);
49  backgrounds_value_ = consumes<std::vector<float> >(Background);
50  centralityTag_ = iConfig.getParameter<InputTag>("centralitycollection");
51  centralityToken = consumes<reco::Centrality>(centralityTag_);
52 
53  centralityBinTag_ = (iConfig.getParameter<edm::InputTag> ("centralitybincollection"));
54  centralityBinToken = consumes<int>(centralityBinTag_);
55  hiVertexToken_ = consumes<std::vector<reco::Vertex> >(edm::InputTag("hiSelectedVertex"));
56 
57  // need to initialize the PF cand histograms : which are also event variables
58  if(isPFJet){
59  mNPFpart = 0;
60  mPFPt = 0;
61  mPFEta = 0;
62  mPFPhi = 0;
63  mPFArea = 0;
64  mSumPFPt = 0;
65  mSumSquaredPFPt = 0;
66  mSumPFPt_HF = 0;
67 
68  mSumPFPt_n5p191_n2p650 = 0;
69  mSumPFPt_n2p650_n2p043 = 0;
70  mSumPFPt_n2p043_n1p740 = 0;
71  mSumPFPt_n1p740_n1p479 = 0;
72  mSumPFPt_n1p479_n1p131 = 0;
73  mSumPFPt_n1p131_n0p783 = 0;
74  mSumPFPt_n0p783_n0p522 = 0;
75  mSumPFPt_n0p522_0p522 = 0;
76  mSumPFPt_0p522_0p783 = 0;
77  mSumPFPt_0p783_1p131 = 0;
78  mSumPFPt_1p131_1p479 = 0;
79  mSumPFPt_1p479_1p740 = 0;
80  mSumPFPt_1p740_2p043 = 0;
81  mSumPFPt_2p043_2p650 = 0;
82  mSumPFPt_2p650_5p191 = 0;
83 
84  mPFCandpT_vs_eta_Unknown = 0; // pf id 0
85  mPFCandpT_vs_eta_ChargedHadron = 0; // pf id - 1
86  mPFCandpT_vs_eta_electron = 0; // pf id - 2
87  mPFCandpT_vs_eta_muon = 0; // pf id - 3
88  mPFCandpT_vs_eta_photon = 0; // pf id - 4
89  mPFCandpT_vs_eta_NeutralHadron = 0; // pf id - 5
90  mPFCandpT_vs_eta_HadE_inHF = 0; // pf id - 6
91  mPFCandpT_vs_eta_EME_inHF = 0; // pf id - 7
92 
93  mPFCandpT_Barrel_Unknown = 0; // pf id 0
94  mPFCandpT_Barrel_ChargedHadron = 0; // pf id - 1
95  mPFCandpT_Barrel_electron = 0; // pf id - 2
96  mPFCandpT_Barrel_muon = 0; // pf id - 3
97  mPFCandpT_Barrel_photon = 0; // pf id - 4
98  mPFCandpT_Barrel_NeutralHadron = 0; // pf id - 5
99  mPFCandpT_Barrel_HadE_inHF = 0; // pf id - 6
100  mPFCandpT_Barrel_EME_inHF = 0; // pf id - 7
101 
102  mPFCandpT_Endcap_Unknown = 0; // pf id 0
103  mPFCandpT_Endcap_ChargedHadron = 0; // pf id - 1
104  mPFCandpT_Endcap_electron = 0; // pf id - 2
105  mPFCandpT_Endcap_muon = 0; // pf id - 3
106  mPFCandpT_Endcap_photon = 0; // pf id - 4
107  mPFCandpT_Endcap_NeutralHadron = 0; // pf id - 5
108  mPFCandpT_Endcap_HadE_inHF = 0; // pf id - 6
109  mPFCandpT_Endcap_EME_inHF = 0; // pf id - 7
110 
111  mPFCandpT_Forward_Unknown = 0; // pf id 0
112  mPFCandpT_Forward_ChargedHadron = 0; // pf id - 1
113  mPFCandpT_Forward_electron = 0; // pf id - 2
114  mPFCandpT_Forward_muon = 0; // pf id - 3
115  mPFCandpT_Forward_photon = 0; // pf id - 4
116  mPFCandpT_Forward_NeutralHadron = 0; // pf id - 5
117  mPFCandpT_Forward_HadE_inHF = 0; // pf id - 6
118  mPFCandpT_Forward_EME_inHF = 0; // pf id - 7
119 
120  }
121  if(isCaloJet){
122  mNCalopart = 0;
123  mCaloPt = 0;
124  mCaloEta = 0;
125  mCaloPhi = 0;
126  mCaloArea = 0;
127 
128  mSumCaloPt = 0;
129  mSumSquaredCaloPt = 0;
130  mSumCaloPt_HF = 0;
131 
132  mSumCaloPt_n5p191_n2p650 = 0;
133  mSumCaloPt_n2p650_n2p043 = 0;
134  mSumCaloPt_n2p043_n1p740 = 0;
135  mSumCaloPt_n1p740_n1p479 = 0;
136  mSumCaloPt_n1p479_n1p131 = 0;
137  mSumCaloPt_n1p131_n0p783 = 0;
138  mSumCaloPt_n0p783_n0p522 = 0;
139  mSumCaloPt_n0p522_0p522 = 0;
140  mSumCaloPt_0p522_0p783 = 0;
141  mSumCaloPt_0p783_1p131 = 0;
142  mSumCaloPt_1p131_1p479 = 0;
143  mSumCaloPt_1p479_1p740 = 0;
144  mSumCaloPt_1p740_2p043 = 0;
145  mSumCaloPt_2p043_2p650 = 0;
146  mSumCaloPt_2p650_5p191 = 0;
147 
148  }
149  mSumpt = 0;
150 
151 
152 
153  // Events variables
154  mNvtx = 0;
155  mHF = 0;
156 
157 
158  // Jet parameters
159  mEta = 0;
160  mPhi = 0;
161  mEnergy = 0;
162  mP = 0;
163  mPt = 0;
164  mMass = 0;
165  mConstituents = 0;
166  mJetArea = 0;
167  mjetpileup = 0;
168  mNJets_40 = 0;
169  mNJets = 0;
170 
171  mGenEta = 0;
172  mGenPhi = 0;
173  mGenPt = 0;
174  mPtHat = 0;
175 
176  mPtRecoOverGen_B_20_30_Cent_0_10 = 0;
177  mPtRecoOverGen_E_20_30_Cent_0_10 = 0;
178  mPtRecoOverGen_F_20_30_Cent_0_10 = 0;
179  mPtRecoOverGen_B_30_50_Cent_0_10 = 0;
180  mPtRecoOverGen_E_30_50_Cent_0_10 = 0;
181  mPtRecoOverGen_F_30_50_Cent_0_10 = 0;
182  mPtRecoOverGen_B_50_80_Cent_0_10 = 0;
183  mPtRecoOverGen_E_50_80_Cent_0_10 = 0;
184  mPtRecoOverGen_F_50_80_Cent_0_10 = 0;
185  mPtRecoOverGen_B_80_120_Cent_0_10 = 0;
186  mPtRecoOverGen_E_80_120_Cent_0_10 = 0;
187  mPtRecoOverGen_F_80_120_Cent_0_10 = 0;
188  mPtRecoOverGen_B_120_180_Cent_0_10 = 0;
189  mPtRecoOverGen_E_120_180_Cent_0_10 = 0;
190  mPtRecoOverGen_F_120_180_Cent_0_10 = 0;
191  mPtRecoOverGen_B_180_300_Cent_0_10 = 0;
192  mPtRecoOverGen_E_180_300_Cent_0_10 = 0;
193  mPtRecoOverGen_F_180_300_Cent_0_10 = 0;
194  mPtRecoOverGen_B_300_Inf_Cent_0_10 = 0;
195  mPtRecoOverGen_E_300_Inf_Cent_0_10 = 0;
196  mPtRecoOverGen_F_300_Inf_Cent_0_10 = 0;
197 
198  mPtRecoOverGen_B_20_30_Cent_10_30 = 0;
199  mPtRecoOverGen_E_20_30_Cent_10_30 = 0;
200  mPtRecoOverGen_F_20_30_Cent_10_30 = 0;
201  mPtRecoOverGen_B_30_50_Cent_10_30 = 0;
202  mPtRecoOverGen_E_30_50_Cent_10_30 = 0;
203  mPtRecoOverGen_F_30_50_Cent_10_30 = 0;
204  mPtRecoOverGen_B_50_80_Cent_10_30 = 0;
205  mPtRecoOverGen_E_50_80_Cent_10_30 = 0;
206  mPtRecoOverGen_F_50_80_Cent_10_30 = 0;
207  mPtRecoOverGen_B_80_120_Cent_10_30 = 0;
208  mPtRecoOverGen_E_80_120_Cent_10_30 = 0;
209  mPtRecoOverGen_F_80_120_Cent_10_30 = 0;
210  mPtRecoOverGen_B_120_180_Cent_10_30 = 0;
211  mPtRecoOverGen_E_120_180_Cent_10_30 = 0;
212  mPtRecoOverGen_F_120_180_Cent_10_30 = 0;
213  mPtRecoOverGen_B_180_300_Cent_10_30 = 0;
214  mPtRecoOverGen_E_180_300_Cent_10_30 = 0;
215  mPtRecoOverGen_F_180_300_Cent_10_30 = 0;
216  mPtRecoOverGen_B_300_Inf_Cent_10_30 = 0;
217  mPtRecoOverGen_E_300_Inf_Cent_10_30 = 0;
218  mPtRecoOverGen_F_300_Inf_Cent_10_30 = 0;
219 
220  mPtRecoOverGen_B_20_30_Cent_30_50 = 0;
221  mPtRecoOverGen_E_20_30_Cent_30_50 = 0;
222  mPtRecoOverGen_F_20_30_Cent_30_50 = 0;
223  mPtRecoOverGen_B_30_50_Cent_30_50 = 0;
224  mPtRecoOverGen_E_30_50_Cent_30_50 = 0;
225  mPtRecoOverGen_F_30_50_Cent_30_50 = 0;
226  mPtRecoOverGen_B_50_80_Cent_30_50 = 0;
227  mPtRecoOverGen_E_50_80_Cent_30_50 = 0;
228  mPtRecoOverGen_F_50_80_Cent_30_50 = 0;
229  mPtRecoOverGen_B_80_120_Cent_30_50 = 0;
230  mPtRecoOverGen_E_80_120_Cent_30_50 = 0;
231  mPtRecoOverGen_F_80_120_Cent_30_50 = 0;
232  mPtRecoOverGen_B_120_180_Cent_30_50 = 0;
233  mPtRecoOverGen_E_120_180_Cent_30_50 = 0;
234  mPtRecoOverGen_F_120_180_Cent_30_50 = 0;
235  mPtRecoOverGen_B_180_300_Cent_30_50 = 0;
236  mPtRecoOverGen_E_180_300_Cent_30_50 = 0;
237  mPtRecoOverGen_F_180_300_Cent_30_50 = 0;
238  mPtRecoOverGen_B_300_Inf_Cent_30_50 = 0;
239  mPtRecoOverGen_E_300_Inf_Cent_30_50 = 0;
240  mPtRecoOverGen_F_300_Inf_Cent_30_50 = 0;
241 
242  mPtRecoOverGen_B_20_30_Cent_50_80 = 0;
243  mPtRecoOverGen_E_20_30_Cent_50_80 = 0;
244  mPtRecoOverGen_F_20_30_Cent_50_80 = 0;
245  mPtRecoOverGen_B_30_50_Cent_50_80 = 0;
246  mPtRecoOverGen_E_30_50_Cent_50_80 = 0;
247  mPtRecoOverGen_F_30_50_Cent_50_80 = 0;
248  mPtRecoOverGen_B_50_80_Cent_50_80 = 0;
249  mPtRecoOverGen_E_50_80_Cent_50_80 = 0;
250  mPtRecoOverGen_F_50_80_Cent_50_80 = 0;
251  mPtRecoOverGen_B_80_120_Cent_50_80 = 0;
252  mPtRecoOverGen_E_80_120_Cent_50_80 = 0;
253  mPtRecoOverGen_F_80_120_Cent_50_80 = 0;
254  mPtRecoOverGen_B_120_180_Cent_50_80 = 0;
255  mPtRecoOverGen_E_120_180_Cent_50_80 = 0;
256  mPtRecoOverGen_F_120_180_Cent_50_80 = 0;
257  mPtRecoOverGen_B_180_300_Cent_50_80 = 0;
258  mPtRecoOverGen_E_180_300_Cent_50_80 = 0;
259  mPtRecoOverGen_F_180_300_Cent_50_80 = 0;
260  mPtRecoOverGen_B_300_Inf_Cent_50_80 = 0;
261  mPtRecoOverGen_E_300_Inf_Cent_50_80 = 0;
262  mPtRecoOverGen_F_300_Inf_Cent_50_80 = 0;
263 
264  mPtRecoOverGen_GenPt_B_Cent_0_10 = 0;
265  mPtRecoOverGen_GenPt_E_Cent_0_10 = 0;
266  mPtRecoOverGen_GenPt_F_Cent_0_10 = 0;
267  mPtRecoOverGen_GenPt_B_Cent_10_30 = 0;
268  mPtRecoOverGen_GenPt_E_Cent_10_30 = 0;
269  mPtRecoOverGen_GenPt_F_Cent_10_30 = 0;
270  mPtRecoOverGen_GenPt_B_Cent_30_50 = 0;
271  mPtRecoOverGen_GenPt_E_Cent_30_50 = 0;
272  mPtRecoOverGen_GenPt_F_Cent_30_50 = 0;
273  mPtRecoOverGen_GenPt_B_Cent_50_80 = 0;
274  mPtRecoOverGen_GenPt_E_Cent_50_80 = 0;
275  mPtRecoOverGen_GenPt_F_Cent_50_80 = 0;
276 
277  mPtRecoOverGen_GenEta_20_30_Cent_0_10 = 0;
278  mPtRecoOverGen_GenEta_30_50_Cent_0_10 = 0;
279  mPtRecoOverGen_GenEta_50_80_Cent_0_10 = 0;
280  mPtRecoOverGen_GenEta_80_120_Cent_0_10 = 0;
281  mPtRecoOverGen_GenEta_120_180_Cent_0_10 = 0;
282  mPtRecoOverGen_GenEta_180_300_Cent_0_10 = 0;
283  mPtRecoOverGen_GenEta_300_Inf_Cent_0_10 = 0;
284 
285  mPtRecoOverGen_GenEta_20_30_Cent_10_30 = 0;
286  mPtRecoOverGen_GenEta_30_50_Cent_10_30 = 0;
287  mPtRecoOverGen_GenEta_50_80_Cent_10_30 = 0;
288  mPtRecoOverGen_GenEta_80_120_Cent_10_30 = 0;
289  mPtRecoOverGen_GenEta_120_180_Cent_10_30 = 0;
290  mPtRecoOverGen_GenEta_180_300_Cent_10_30 = 0;
291  mPtRecoOverGen_GenEta_300_Inf_Cent_10_30 = 0;
292 
293  mPtRecoOverGen_GenEta_20_30_Cent_30_50 = 0;
294  mPtRecoOverGen_GenEta_30_50_Cent_30_50 = 0;
295  mPtRecoOverGen_GenEta_50_80_Cent_30_50 = 0;
296  mPtRecoOverGen_GenEta_80_120_Cent_30_50 = 0;
297  mPtRecoOverGen_GenEta_120_180_Cent_30_50 = 0;
298  mPtRecoOverGen_GenEta_180_300_Cent_30_50 = 0;
299  mPtRecoOverGen_GenEta_300_Inf_Cent_30_50 = 0;
300 
301  mPtRecoOverGen_GenEta_20_30_Cent_50_80 = 0;
302  mPtRecoOverGen_GenEta_30_50_Cent_50_80 = 0;
303  mPtRecoOverGen_GenEta_50_80_Cent_50_80 = 0;
304  mPtRecoOverGen_GenEta_80_120_Cent_50_80 = 0;
305  mPtRecoOverGen_GenEta_120_180_Cent_50_80 = 0;
306  mPtRecoOverGen_GenEta_180_300_Cent_50_80 = 0;
307  mPtRecoOverGen_GenEta_300_Inf_Cent_50_80 = 0;
308 
309 
310 }
311 
312 void JetTester_HeavyIons::bookHistograms(DQMStore::IBooker & ibooker, edm::Run const & iRun,edm::EventSetup const &)
313 {
314 
315  ibooker.setCurrentFolder("JetMET/JetValidation/"+mInputCollection.label());
316 
317  double log10PtMin = 0.50;
318  double log10PtMax = 3.75;
319  int log10PtBins = 26;
320 
321  static const size_t ncms_hcal_edge_pseudorapidity = 82 + 1;
322  static const double cms_hcal_edge_pseudorapidity[ncms_hcal_edge_pseudorapidity] = {
323  -5.191, -4.889, -4.716, -4.538, -4.363, -4.191, -4.013,
324  -3.839, -3.664, -3.489, -3.314, -3.139, -2.964, -2.853,
325  -2.650, -2.500, -2.322, -2.172, -2.043, -1.930, -1.830,
326  -1.740, -1.653, -1.566, -1.479, -1.392, -1.305, -1.218,
327  -1.131, -1.044, -0.957, -0.879, -0.783, -0.696, -0.609,
328  -0.522, -0.435, -0.348, -0.261, -0.174, -0.087,
329  0.000,
330  0.087, 0.174, 0.261, 0.348, 0.435, 0.522, 0.609,
331  0.696, 0.783, 0.879, 0.957, 1.044, 1.131, 1.218,
332  1.305, 1.392, 1.479, 1.566, 1.653, 1.740, 1.830,
333  1.930, 2.043, 2.172, 2.322, 2.500, 2.650, 2.853,
334  2.964, 3.139, 3.314, 3.489, 3.664, 3.839, 4.013,
335  4.191, 4.363, 4.538, 4.716, 4.889, 5.191
336  };
337 
338  double etaRange[91] = {-6.0, -5.8, -5.6, -5.4, -5.2, -5.0, -4.8, -4.6, -4.4, -4.2,
339  -4.0, -3.8, -3.6, -3.4, -3.2, -3.0, -2.9, -2.8, -2.7, -2.6,
340  -2.5, -2.4, -2.3, -2.2, -2.1, -2.0, -1.9, -1.8, -1.7, -1.6,
341  -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.9, -0.8, -0.7, -0.6,
342  -0.5, -0.4, -0.3, -0.2, -0.1,
343  0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9,
344  1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9,
345  2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9,
346  3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8,
347  5.0, 5.2, 5.4, 5.6, 5.8, 6.0};
348 
349  double edge_pseudorapidity[etaBins_ +1] = {-5.191, -2.650, -2.043, -1.740, -1.479, -1.131, -0.783, -0.522, 0.522, 0.783, 1.131, 1.479, 1.740, 2.043, 2.650, 5.191 };
350 
351  TH2F *h2D_etabins_vs_pt2 = new TH2F("h2D_etabins_vs_pt2","etaBins (x axis), sum pt^{2} (y axis)",etaBins_,edge_pseudorapidity,10000,0,10000);
352  TH2F *h2D_etabins_vs_pt = new TH2F("h2D_etabins_vs_pt","etaBins (x axis), sum pt (y axis)",etaBins_,edge_pseudorapidity,10000,-1000,1000);
353  TH2F *h2D_etabins_vs_phi = new TH2F("h2D_etabins_vs_phi","candidate map, eta(x axis), phi (y axis), pt (z axis)",ncms_hcal_edge_pseudorapidity-1,cms_hcal_edge_pseudorapidity,36,-TMath::Pi(),TMath::Pi());
354  TH2F *h2D_pfcand_etabins_vs_pt = new TH2F("h2D_etabins_vs_pt",";#eta;sum p_{T}",etaBins_,edge_pseudorapidity,300, 0, 300);
355 
356  if(isPFJet){
357 
358  mNPFpart = ibooker.book1D("NPFpart","No of particle flow candidates",1000,0,10000);
359  mPFPt = ibooker.book1D("PFPt","PF candidate p_{T}",1000,-5000,5000);
360  mPFEta = ibooker.book1D("PFEta","PF candidate #eta",120,-6,6);
361  mPFPhi = ibooker.book1D("PFPhi","PF candidate #phi",70,-3.5,3.5);
362  mPFArea = ibooker.book1D("PFArea","VS PF candidate area",100,0,4);
363 
364  mSumPFPt = ibooker.book1D("SumPFPt","Sum of initial PF p_{T}",1000,-10000,10000);
365 
366  mSumSquaredPFPt = ibooker.book1D("SumSquaredPFPt","Sum of initial PF p_{T} squared",10000,0,10000);
367 
368  mSumPFPt_HF = ibooker.book2D("SumPFPt_HF","HF energy (y axis) vs Sum initial PF p_{T} (x axis)",1000,-1000,1000,1000,0,10000);
369 
370 
371  mSumPFPt_n5p191_n2p650 = ibooker.book1D("mSumPFPt_n5p191_n2p650","Sum PFPt in the eta range -5.191 to -2.650",1000,-5000,5000);
372  mSumPFPt_n2p650_n2p043 = ibooker.book1D("mSumPFPt_n2p650_n2p043","Sum PFPt in the eta range -2.650 to -2.043 ",1000,-5000,5000);
373  mSumPFPt_n2p043_n1p740 = ibooker.book1D("mSumPFPt_n2p043_n1p740","Sum PFPt in the eta range -2.043 to -1.740",1000,-1000,1000);
374  mSumPFPt_n1p740_n1p479 = ibooker.book1D("mSumPFPt_n1p740_n1p479","Sum PFPt in the eta range -1.740 to -1.479",1000,-1000,1000);
375  mSumPFPt_n1p479_n1p131 = ibooker.book1D("mSumPFPt_n1p479_n1p131","Sum PFPt in the eta range -1.479 to -1.131",1000,-1000,1000);
376  mSumPFPt_n1p131_n0p783 = ibooker.book1D("mSumPFPt_n1p131_n0p783","Sum PFPt in the eta range -1.131 to -0.783",1000,-1000,1000);
377  mSumPFPt_n0p783_n0p522 = ibooker.book1D("mSumPFPt_n0p783_n0p522","Sum PFPt in the eta range -0.783 to -0.522",1000,-1000,1000);
378  mSumPFPt_n0p522_0p522 = ibooker.book1D("mSumPFPt_n0p522_0p522","Sum PFPt in the eta range -0.522 to 0.522",1000,-1000,1000);
379  mSumPFPt_0p522_0p783 = ibooker.book1D("mSumPFPt_0p522_0p783","Sum PFPt in the eta range 0.522 to 0.783",1000,-1000,1000);
380  mSumPFPt_0p783_1p131 = ibooker.book1D("mSumPFPt_0p783_1p131","Sum PFPt in the eta range 0.783 to 1.131",1000,-1000,1000);
381  mSumPFPt_1p131_1p479 = ibooker.book1D("mSumPFPt_1p131_1p479","Sum PFPt in the eta range 1.131 to 1.479",1000,-1000,1000);
382  mSumPFPt_1p479_1p740 = ibooker.book1D("mSumPFPt_1p479_1p740","Sum PFPt in the eta range 1.479 to 1.740",1000,-1000,1000);
383  mSumPFPt_1p740_2p043 = ibooker.book1D("mSumPFPt_1p740_2p043","Sum PFPt in the eta range 1.740 to 2.043",1000,-1000,1000);
384  mSumPFPt_2p043_2p650 = ibooker.book1D("mSumPFPt_2p043_2p650","Sum PFPt in the eta range 2.043 to 2.650",1000,-5000,5000);
385  mSumPFPt_2p650_5p191 = ibooker.book1D("mSumPFPt_2p650_5p191","Sum PFPt in the eta range 2.650 to 5.191",1000,-5000,5000);
386 
387  mPFCandpT_vs_eta_Unknown = ibooker.book2D("PF_cand_X_unknown",h2D_pfcand_etabins_vs_pt); // pf id 0
388  mPFCandpT_vs_eta_ChargedHadron = ibooker.book2D("PF_cand_chargedHad",h2D_pfcand_etabins_vs_pt); // pf id - 1
389  mPFCandpT_vs_eta_electron = ibooker.book2D("PF_cand_electron",h2D_pfcand_etabins_vs_pt); // pf id - 2
390  mPFCandpT_vs_eta_muon = ibooker.book2D("PF_cand_muon",h2D_pfcand_etabins_vs_pt); // pf id - 3
391  mPFCandpT_vs_eta_photon = ibooker.book2D("PF_cand_photon",h2D_pfcand_etabins_vs_pt); // pf id - 4
392  mPFCandpT_vs_eta_NeutralHadron = ibooker.book2D("PF_cand_neutralHad",h2D_pfcand_etabins_vs_pt); // pf id - 5
393  mPFCandpT_vs_eta_HadE_inHF = ibooker.book2D("PF_cand_HadEner_inHF",h2D_pfcand_etabins_vs_pt); // pf id - 6
394  mPFCandpT_vs_eta_EME_inHF = ibooker.book2D("PF_cand_EMEner_inHF",h2D_pfcand_etabins_vs_pt); // pf id - 7
395 
396  mPFCandpT_Barrel_Unknown = ibooker.book1D("mPFCandpT_Barrel_Unknown",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 0
397  mPFCandpT_Barrel_ChargedHadron = ibooker.book1D("mPFCandpT_Barrel_ChargedHadron",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 1
398  mPFCandpT_Barrel_electron = ibooker.book1D("mPFCandpT_Barrel_electron",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 2
399  mPFCandpT_Barrel_muon = ibooker.book1D("mPFCandpT_Barrel_muon",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 3
400  mPFCandpT_Barrel_photon = ibooker.book1D("mPFCandpT_Barrel_photon",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 4
401  mPFCandpT_Barrel_NeutralHadron = ibooker.book1D("mPFCandpT_Barrel_NeutralHadron",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 5
402  mPFCandpT_Barrel_HadE_inHF = ibooker.book1D("mPFCandpT_Barrel_HadE_inHF",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 6
403  mPFCandpT_Barrel_EME_inHF = ibooker.book1D("mPFCandpT_Barrel_EME_inHF",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 7
404 
405  mPFCandpT_Endcap_Unknown = ibooker.book1D("mPFCandpT_Endcap_Unknown",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", BarrelEta, EndcapEta),300, 0, 300); // pf id - 0
406  mPFCandpT_Endcap_ChargedHadron = ibooker.book1D("mPFCandpT_Endcap_ChargedHadron",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", BarrelEta, EndcapEta),300, 0, 300); // pf id - 1
407  mPFCandpT_Endcap_electron = ibooker.book1D("mPFCandpT_Endcap_electron",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", BarrelEta, EndcapEta),300, 0, 300); // pf id - 2
408  mPFCandpT_Endcap_muon = ibooker.book1D("mPFCandpT_Endcap_muon",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", BarrelEta, EndcapEta),300, 0, 300); // pf id - 3
409  mPFCandpT_Endcap_photon = ibooker.book1D("mPFCandpT_Endcap_photon",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", BarrelEta, EndcapEta),300, 0, 300); // pf id - 4
410  mPFCandpT_Endcap_NeutralHadron = ibooker.book1D("mPFCandpT_Endcap_NeutralHadron",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", BarrelEta, EndcapEta),300, 0, 300); // pf id - 5
411  mPFCandpT_Endcap_HadE_inHF = ibooker.book1D("mPFCandpT_Endcap_HadE_inHF",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", BarrelEta, EndcapEta),300, 0, 300); // pf id - 6
412  mPFCandpT_Endcap_EME_inHF = ibooker.book1D("mPFCandpT_Endcap_EME_inHF",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", BarrelEta, EndcapEta),300, 0, 300); // pf id - 7
413 
414  mPFCandpT_Forward_Unknown = ibooker.book1D("mPFCandpT_Forward_Unknown",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", EndcapEta, ForwardEta),300, 0, 300); // pf id - 0
415  mPFCandpT_Forward_ChargedHadron = ibooker.book1D("mPFCandpT_Forward_ChargedHadron",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", EndcapEta, ForwardEta),300, 0, 300); // pf id - 1
416  mPFCandpT_Forward_electron = ibooker.book1D("mPFCandpT_Forward_electron",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", EndcapEta, ForwardEta),300, 0, 300); // pf id - 2
417  mPFCandpT_Forward_muon = ibooker.book1D("mPFCandpT_Forward_muon",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", EndcapEta, ForwardEta),300, 0, 300); // pf id - 3
418  mPFCandpT_Forward_photon = ibooker.book1D("mPFCandpT_Forward_photon",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", EndcapEta, ForwardEta),300, 0, 300); // pf id - 4
419  mPFCandpT_Forward_NeutralHadron = ibooker.book1D("mPFCandpT_Forward_NeutralHadron",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", EndcapEta, ForwardEta),300, 0, 300); // pf id - 5
420  mPFCandpT_Forward_HadE_inHF = ibooker.book1D("mPFCandpT_Forward_HadE_inHF",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", EndcapEta, ForwardEta),300, 0, 300); // pf id - 6
421  mPFCandpT_Forward_EME_inHF = ibooker.book1D("mPFCandpT_Forward_EME_inHF",Form(";PF candidate p_{T}, %2.2f<|#eta|<%2.2f; counts", EndcapEta, ForwardEta),300, 0, 300); // pf id - 7
422 
423 
424  }
425 
426  if(isCaloJet){
427 
428  mNCalopart = ibooker.book1D("NCalopart","No of particle flow candidates",1000,0,10000);
429  mCaloPt = ibooker.book1D("CaloPt","Calo candidate p_{T}",1000,-5000,5000);
430  mCaloEta = ibooker.book1D("CaloEta","Calo candidate #eta",120,-6,6);
431  mCaloPhi = ibooker.book1D("CaloPhi","Calo candidate #phi",70,-3.5,3.5);
432  mCaloArea = ibooker.book1D("CaloArea","VS Calo candidate area",100,0,4);
433 
434  mSumCaloPt = ibooker.book1D("SumCaloPt","Sum Calo p_{T}",1000,-10000,10000);
435 
436  mSumSquaredCaloPt = ibooker.book1D("SumSquaredCaloPt","Sum of initial Calo tower p_{T} squared",10000,0,10000);
437 
438  mSumCaloPt_HF = ibooker.book2D("SumCaloPt_HF","HF Energy (y axis) vs Sum Calo tower p_{T}",1000,-1000,1000,1000,0,10000);
439 
440 
441  mSumCaloPt_n5p191_n2p650 = ibooker.book1D("mSumCaloPt_n5p191_n2p650","Sum Calo tower pT variable in the eta range -5.191 to -2.650",1000,-5000,5000);
442  mSumCaloPt_n2p650_n2p043 = ibooker.book1D("mSumCaloPt_n2p650_n2p043","Sum Calo tower pT variable in the eta range -2.650 to -2.043",1000,-5000,5000);
443  mSumCaloPt_n2p043_n1p740 = ibooker.book1D("mSumCaloPt_n2p043_n1p740","Sum Calo tower pT variable in the eta range -2.043 to -1.740",1000,-1000,1000);
444  mSumCaloPt_n1p740_n1p479 = ibooker.book1D("mSumCaloPt_n1p740_n1p479","Sum Calo tower pT variable in the eta range -1.740 to -1.479",1000,-1000,1000);
445  mSumCaloPt_n1p479_n1p131 = ibooker.book1D("mSumCaloPt_n1p479_n1p131","Sum Calo tower pT variable in the eta range -1.479 to -1.131",1000,-1000,1000);
446  mSumCaloPt_n1p131_n0p783 = ibooker.book1D("mSumCaloPt_n1p131_n0p783","Sum Calo tower pT variable in the eta range -1.131 to -0.783",1000,-1000,1000);
447  mSumCaloPt_n0p783_n0p522 = ibooker.book1D("mSumCaloPt_n0p783_n0p522","Sum Calo tower pT variable in the eta range -0.783 to -0.522",1000,-1000,1000);
448  mSumCaloPt_n0p522_0p522 = ibooker.book1D("mSumCaloPt_n0p522_0p522","Sum Calo tower pT variable in the eta range -0.522 to 0.522",1000,-1000,1000);
449  mSumCaloPt_0p522_0p783 = ibooker.book1D("mSumCaloPt_0p522_0p783","Sum Calo tower pT variable in the eta range 0.522 to 0.783",1000,-1000,1000);
450  mSumCaloPt_0p783_1p131 = ibooker.book1D("mSumCaloPt_0p783_1p131","Sum Calo tower pT variable in the eta range 0.783 to 1.131",1000,-1000,1000);
451  mSumCaloPt_1p131_1p479 = ibooker.book1D("mSumCaloPt_1p131_1p479","Sum Calo tower pT variable in the eta range 1.131 to 1.479",1000,-1000,1000);
452  mSumCaloPt_1p479_1p740 = ibooker.book1D("mSumCaloPt_1p479_1p740","Sum Calo tower pT variable in the eta range 1.479 to 1.740",1000,-1000,1000);
453  mSumCaloPt_1p740_2p043 = ibooker.book1D("mSumCaloPt_1p740_2p043","Sum Calo tower pT variable in the eta range 1.740 to 2.043",1000,-1000,1000);
454  mSumCaloPt_2p043_2p650 = ibooker.book1D("mSumCaloPt_2p043_2p650","Sum Calo tower pT variable in the eta range 2.043 to 2.650",1000,-5000,5000);
455  mSumCaloPt_2p650_5p191 = ibooker.book1D("mSumCaloPt_2p650_5p191","Sum Calo tower pT variable in the eta range 2.650 to 5.191",1000,-5000,5000);
456 
457  }
458 
459  // particle flow variables histograms
460  mSumpt = ibooker.book1D("SumpT","Sum p_{T} of all the PF candidates per event",1000,0,10000);
461 
462  // Event variables
463  mNvtx = ibooker.book1D("Nvtx", "number of vertices", 60, 0, 60);
464  mHF = ibooker.book1D("HF", "HF energy distribution",1000,0,10000);
465 
466 
467  // Jet parameters
468  mEta = ibooker.book1D("Eta", "Eta", 120, -6, 6);
469  mPhi = ibooker.book1D("Phi", "Phi", 70, -3.5, 3.5);
470  mPt = ibooker.book1D("Pt", "Pt", 100, 0, 1000);
471  mP = ibooker.book1D("P", "P", 100, 0, 1000);
472  mEnergy = ibooker.book1D("Energy", "Energy", 100, 0, 1000);
473  mMass = ibooker.book1D("Mass", "Mass", 100, 0, 200);
474  mConstituents = ibooker.book1D("Constituents", "Constituents", 100, 0, 100);
475  mJetArea = ibooker.book1D("JetArea", "JetArea", 100, 0, 4);
476  mjetpileup = ibooker.book1D("jetPileUp","jetPileUp",100,0,150);
477  mNJets_40 = ibooker.book1D("NJets_pt_greater_40", "NJets pT > 40 GeV", 50, 0, 100);
478  mNJets = ibooker.book1D("NJets", "NJets", 50, 0, 100);
479 
480  mGenEta = ibooker.book1D("Gen Eta",";gen jet #eta;counts",120,-6,6);
481  mGenPhi = ibooker.book1D("Gen Phi",";gen jet #phi;counts",70,-3.5,3.5);
482  mGenPt = ibooker.book1D("Gen pT","gen jet p_{T}",250,0,1000);
483  mPtHat = ibooker.book1D("pThat","#hat{p_{T}}",250,0,1000);
484 
485  mPtRecoOverGen_B_20_30_Cent_0_10 = ibooker.book1D("PtRecoOverGen_B_20_30_Cent_0_10", "20<genpt<30; recopt/genpt (0-10%) (Barrel);counts", 90, 0, 2);
486  mPtRecoOverGen_E_20_30_Cent_0_10 = ibooker.book1D("PtRecoOverGen_E_20_30_Cent_0_10", "20<genpt<30; recopt/genpt (0-10%) (EndCap);counts", 90, 0, 2);
487  mPtRecoOverGen_F_20_30_Cent_0_10 = ibooker.book1D("PtRecoOverGen_F_20_30_Cent_0_10", "20<genpt<30; recopt/genpt (0-10%) (Forward);counts", 90, 0, 2);
488  mPtRecoOverGen_B_30_50_Cent_0_10 = ibooker.book1D("PtRecoOverGen_B_30_50_Cent_0_10", "30<genpt<50; recopt/genpt (0-10%) (Barrel);counts", 90, 0, 2);
489  mPtRecoOverGen_E_30_50_Cent_0_10 = ibooker.book1D("PtRecoOverGen_E_30_50_Cent_0_10", "30<genpt<50; recopt/genpt (0-10%) (EndCap);counts", 90, 0, 2);
490  mPtRecoOverGen_F_30_50_Cent_0_10 = ibooker.book1D("PtRecoOverGen_F_30_50_Cent_0_10", "30<genpt<50; recopt/genpt (0-10%) (Forward);counts", 90, 0, 2);
491  mPtRecoOverGen_B_50_80_Cent_0_10 = ibooker.book1D("PtRecoOverGen_B_50_80_Cent_0_10", "50<genpt<80; recopt/genpt (0-10%) (Barrel);counts", 90, 0, 2);
492  mPtRecoOverGen_E_50_80_Cent_0_10 = ibooker.book1D("PtRecoOverGen_E_50_80_Cent_0_10", "50<genpt<80; recopt/genpt (0-10%) (EndCap);counts", 90, 0, 2);
493  mPtRecoOverGen_F_50_80_Cent_0_10 = ibooker.book1D("PtRecoOverGen_F_50_80_Cent_0_10", "50<genpt<80; recopt/genpt (0-10%) (Forward);counts", 90, 0, 2);
494  mPtRecoOverGen_B_80_120_Cent_0_10 = ibooker.book1D("PtRecoOverGen_B_80_120_Cent_0_10", "80<genpt<120; recopt/genpt (0-10%) (Barrel);counts", 90, 0, 2);
495  mPtRecoOverGen_E_80_120_Cent_0_10 = ibooker.book1D("PtRecoOverGen_E_80_120_Cent_0_10", "80<genpt<120; recopt/genpt (0-10%) (EndCap);counts", 90, 0, 2);
496  mPtRecoOverGen_F_80_120_Cent_0_10 = ibooker.book1D("PtRecoOverGen_F_80_120_Cent_0_10", "80<genpt<120; recopt/genpt (0-10%) (Forward);counts", 90, 0, 2);
497  mPtRecoOverGen_B_120_180_Cent_0_10 = ibooker.book1D("PtRecoOverGen_B_120_180_Cent_0_10", "120<genpt<180; recopt/genpt (0-10%) (Barrel);counts", 90, 0, 2);
498  mPtRecoOverGen_E_120_180_Cent_0_10 = ibooker.book1D("PtRecoOverGen_E_120_180_Cent_0_10", "120<genpt<180; recopt/genpt (0-10%) (EndCap);counts", 90, 0, 2);
499  mPtRecoOverGen_F_120_180_Cent_0_10 = ibooker.book1D("PtRecoOverGen_F_120_180_Cent_0_10", "120<genpt<180; recopt/genpt (0-10%) (Forward);counts", 90, 0, 2);
500  mPtRecoOverGen_B_180_300_Cent_0_10 = ibooker.book1D("PtRecoOverGen_B_180_300_Cent_0_10", "180<genpt<300; recopt/genpt (0-10%) (Barrel);counts", 90, 0, 2);
501  mPtRecoOverGen_E_180_300_Cent_0_10 = ibooker.book1D("PtRecoOverGen_E_180_300_Cent_0_10", "180<genpt<300; recopt/genpt (0-10%) (EndCap);counts", 90, 0, 2);
502  mPtRecoOverGen_F_180_300_Cent_0_10 = ibooker.book1D("PtRecoOverGen_F_180_300_Cent_0_10", "180<genpt<300; recopt/genpt (0-10%) (Forward);counts", 90, 0, 2);
503  mPtRecoOverGen_B_300_Inf_Cent_0_10 = ibooker.book1D("PtRecoOverGen_B_300_Inf_Cent_0_10", "300<genpt<Inf; recopt/genpt (0-10%) (Barrel);counts", 90, 0, 2);
504  mPtRecoOverGen_E_300_Inf_Cent_0_10 = ibooker.book1D("PtRecoOverGen_E_300_Inf_Cent_0_10", "300<genpt<Inf; recopt/genpt (0-10%) (EndCap);counts", 90, 0, 2);
505  mPtRecoOverGen_F_300_Inf_Cent_0_10 = ibooker.book1D("PtRecoOverGen_F_300_Inf_Cent_0_10", "300<genpt<Inf; recopt/genpt (0-10%) (Forward);counts", 90, 0, 2);
506 
507 
508  mPtRecoOverGen_B_20_30_Cent_10_30 = ibooker.book1D("PtRecoOverGen_B_20_30_Cent_10_30", "20<genpt<30; recopt/genpt (10-30%) (Barrel);counts", 90, 0, 2);
509  mPtRecoOverGen_E_20_30_Cent_10_30 = ibooker.book1D("PtRecoOverGen_E_20_30_Cent_10_30", "20<genpt<30; recopt/genpt (10-30%) (EndCap);counts", 90, 0, 2);
510  mPtRecoOverGen_F_20_30_Cent_10_30 = ibooker.book1D("PtRecoOverGen_F_20_30_Cent_10_30", "20<genpt<30; recopt/genpt (10-30%) (Forward);counts", 90, 0, 2);
511  mPtRecoOverGen_B_30_50_Cent_10_30 = ibooker.book1D("PtRecoOverGen_B_30_50_Cent_10_30", "30<genpt<50; recopt/genpt (10-30%) (Barrel);counts", 90, 0, 2);
512  mPtRecoOverGen_E_30_50_Cent_10_30 = ibooker.book1D("PtRecoOverGen_E_30_50_Cent_10_30", "30<genpt<50; recopt/genpt (10-30%) (EndCap);counts", 90, 0, 2);
513  mPtRecoOverGen_F_30_50_Cent_10_30 = ibooker.book1D("PtRecoOverGen_F_30_50_Cent_10_30", "30<genpt<50; recopt/genpt (10-30%) (Forward);counts", 90, 0, 2);
514  mPtRecoOverGen_B_50_80_Cent_10_30 = ibooker.book1D("PtRecoOverGen_B_50_80_Cent_10_30", "50<genpt<80; recopt/genpt (10-30%) (Barrel);counts", 90, 0, 2);
515  mPtRecoOverGen_E_50_80_Cent_10_30 = ibooker.book1D("PtRecoOverGen_E_50_80_Cent_10_30", "50<genpt<80; recopt/genpt (10-30%) (EndCap);counts", 90, 0, 2);
516  mPtRecoOverGen_F_50_80_Cent_10_30 = ibooker.book1D("PtRecoOverGen_F_50_80_Cent_10_30", "50<genpt<80; recopt/genpt (10-30%) (Forward);counts", 90, 0, 2);
517  mPtRecoOverGen_B_80_120_Cent_10_30 = ibooker.book1D("PtRecoOverGen_B_80_120_Cent_10_30", "80<genpt<120; recopt/genpt (10-30%) (Barrel);counts", 90, 0, 2);
518  mPtRecoOverGen_E_80_120_Cent_10_30 = ibooker.book1D("PtRecoOverGen_E_80_120_Cent_10_30", "80<genpt<120; recopt/genpt (10-30%) (EndCap);counts", 90, 0, 2);
519  mPtRecoOverGen_F_80_120_Cent_10_30 = ibooker.book1D("PtRecoOverGen_F_80_120_Cent_10_30", "80<genpt<120; recopt/genpt (10-30%) (Forward);counts", 90, 0, 2);
520  mPtRecoOverGen_B_120_180_Cent_10_30 = ibooker.book1D("PtRecoOverGen_B_120_180_Cent_10_30", "120<genpt<180; recopt/genpt (10-30%) (Barrel);counts", 90, 0, 2);
521  mPtRecoOverGen_E_120_180_Cent_10_30 = ibooker.book1D("PtRecoOverGen_E_120_180_Cent_10_30", "120<genpt<180; recopt/genpt (10-30%) (EndCap);counts", 90, 0, 2);
522  mPtRecoOverGen_F_120_180_Cent_10_30 = ibooker.book1D("PtRecoOverGen_F_120_180_Cent_10_30", "120<genpt<180; recopt/genpt (10-30%) (Forward);counts", 90, 0, 2);
523  mPtRecoOverGen_B_180_300_Cent_10_30 = ibooker.book1D("PtRecoOverGen_B_180_300_Cent_10_30", "180<genpt<300; recopt/genpt (10-30%) (Barrel);counts", 90, 0, 2);
524  mPtRecoOverGen_E_180_300_Cent_10_30 = ibooker.book1D("PtRecoOverGen_E_180_300_Cent_10_30", "180<genpt<300; recopt/genpt (10-30%) (EndCap);counts", 90, 0, 2);
525  mPtRecoOverGen_F_180_300_Cent_10_30 = ibooker.book1D("PtRecoOverGen_F_180_300_Cent_10_30", "180<genpt<300; recopt/genpt (10-30%) (Forward);counts", 90, 0, 2);
526  mPtRecoOverGen_B_300_Inf_Cent_10_30 = ibooker.book1D("PtRecoOverGen_B_300_Inf_Cent_10_30", "300<genpt<Inf; recopt/genpt (10-30%) (Barrel);counts", 90, 0, 2);
527  mPtRecoOverGen_E_300_Inf_Cent_10_30 = ibooker.book1D("PtRecoOverGen_E_300_Inf_Cent_10_30", "300<genpt<Inf; recopt/genpt (10-30%) (EndCap);counts", 90, 0, 2);
528  mPtRecoOverGen_F_300_Inf_Cent_10_30 = ibooker.book1D("PtRecoOverGen_F_300_Inf_Cent_10_30", "300<genpt<Inf; recopt/genpt (10-30%) (Forward);counts", 90, 0, 2);
529 
530 
531  mPtRecoOverGen_B_20_30_Cent_30_50 = ibooker.book1D("PtRecoOverGen_B_20_30_Cent_30_50", "20<genpt<30; recopt/genpt (30-50%) (Barrel);counts", 90, 0, 2);
532  mPtRecoOverGen_E_20_30_Cent_30_50 = ibooker.book1D("PtRecoOverGen_E_20_30_Cent_30_50", "20<genpt<30; recopt/genpt (30-50%) (EndCap);counts", 90, 0, 2);
533  mPtRecoOverGen_F_20_30_Cent_30_50 = ibooker.book1D("PtRecoOverGen_F_20_30_Cent_30_50", "20<genpt<30; recopt/genpt (30-50%) (Forward);counts", 90, 0, 2);
534  mPtRecoOverGen_B_30_50_Cent_30_50 = ibooker.book1D("PtRecoOverGen_B_30_50_Cent_30_50", "30<genpt<50; recopt/genpt (30-50%) (Barrel);counts", 90, 0, 2);
535  mPtRecoOverGen_E_30_50_Cent_30_50 = ibooker.book1D("PtRecoOverGen_E_30_50_Cent_30_50", "30<genpt<50; recopt/genpt (30-50%) (EndCap);counts", 90, 0, 2);
536  mPtRecoOverGen_F_30_50_Cent_30_50 = ibooker.book1D("PtRecoOverGen_F_30_50_Cent_30_50", "30<genpt<50; recopt/genpt (30-50%) (Forward);counts", 90, 0, 2);
537  mPtRecoOverGen_B_50_80_Cent_30_50 = ibooker.book1D("PtRecoOverGen_B_50_80_Cent_30_50", "50<genpt<80; recopt/genpt (30-50%) (Barrel);counts", 90, 0, 2);
538  mPtRecoOverGen_E_50_80_Cent_30_50 = ibooker.book1D("PtRecoOverGen_E_50_80_Cent_30_50", "50<genpt<80; recopt/genpt (30-50%) (EndCap);counts", 90, 0, 2);
539  mPtRecoOverGen_F_50_80_Cent_30_50 = ibooker.book1D("PtRecoOverGen_F_50_80_Cent_30_50", "50<genpt<80; recopt/genpt (30-50%) (Forward);counts", 90, 0, 2);
540  mPtRecoOverGen_B_80_120_Cent_30_50 = ibooker.book1D("PtRecoOverGen_B_80_120_Cent_30_50", "80<genpt<120; recopt/genpt (30-50%) (Barrel);counts", 90, 0, 2);
541  mPtRecoOverGen_E_80_120_Cent_30_50 = ibooker.book1D("PtRecoOverGen_E_80_120_Cent_30_50", "80<genpt<120; recopt/genpt (30-50%) (EndCap);counts", 90, 0, 2);
542  mPtRecoOverGen_F_80_120_Cent_30_50 = ibooker.book1D("PtRecoOverGen_F_80_120_Cent_30_50", "80<genpt<120; recopt/genpt (30-50%) (Forward);counts", 90, 0, 2);
543  mPtRecoOverGen_B_120_180_Cent_30_50 = ibooker.book1D("PtRecoOverGen_B_120_180_Cent_30_50", "120<genpt<180; recopt/genpt (30-50%) (Barrel);counts", 90, 0, 2);
544  mPtRecoOverGen_E_120_180_Cent_30_50 = ibooker.book1D("PtRecoOverGen_E_120_180_Cent_30_50", "120<genpt<180; recopt/genpt (30-50%) (EndCap);counts", 90, 0, 2);
545  mPtRecoOverGen_F_120_180_Cent_30_50 = ibooker.book1D("PtRecoOverGen_F_120_180_Cent_30_50", "120<genpt<180; recopt/genpt (30-50%) (Forward);counts", 90, 0, 2);
546  mPtRecoOverGen_B_180_300_Cent_30_50 = ibooker.book1D("PtRecoOverGen_B_180_300_Cent_30_50", "180<genpt<300; recopt/genpt (30-50%) (Barrel);counts", 90, 0, 2);
547  mPtRecoOverGen_E_180_300_Cent_30_50 = ibooker.book1D("PtRecoOverGen_E_180_300_Cent_30_50", "180<genpt<300; recopt/genpt (30-50%) (EndCap);counts", 90, 0, 2);
548  mPtRecoOverGen_F_180_300_Cent_30_50 = ibooker.book1D("PtRecoOverGen_F_180_300_Cent_30_50", "180<genpt<300; recopt/genpt (30-50%) (Forward);counts", 90, 0, 2);
549  mPtRecoOverGen_B_300_Inf_Cent_30_50 = ibooker.book1D("PtRecoOverGen_B_300_Inf_Cent_30_50", "300<genpt<Inf; recopt/genpt (30-50%) (Barrel);counts", 90, 0, 2);
550  mPtRecoOverGen_E_300_Inf_Cent_30_50 = ibooker.book1D("PtRecoOverGen_E_300_Inf_Cent_30_50", "300<genpt<Inf; recopt/genpt (30-50%) (EndCap);counts", 90, 0, 2);
551  mPtRecoOverGen_F_300_Inf_Cent_30_50 = ibooker.book1D("PtRecoOverGen_F_300_Inf_Cent_30_50", "300<genpt<Inf; recopt/genpt (30-50%) (Forward);counts", 90, 0, 2);
552 
553  mPtRecoOverGen_B_20_30_Cent_50_80 = ibooker.book1D("PtRecoOverGen_B_20_30_Cent_50_80", "20<genpt<30; recopt/genpt (50-80%) (Barrel);counts", 90, 0, 2);
554  mPtRecoOverGen_E_20_30_Cent_50_80 = ibooker.book1D("PtRecoOverGen_E_20_30_Cent_50_80", "20<genpt<30; recopt/genpt (50-80%) (EndCap);counts", 90, 0, 2);
555  mPtRecoOverGen_F_20_30_Cent_50_80 = ibooker.book1D("PtRecoOverGen_F_20_30_Cent_50_80", "20<genpt<30; recopt/genpt (50-80%) (Forward);counts", 90, 0, 2);
556  mPtRecoOverGen_B_30_50_Cent_50_80 = ibooker.book1D("PtRecoOverGen_B_30_50_Cent_50_80", "30<genpt<50; recopt/genpt (50-80%) (Barrel);counts", 90, 0, 2);
557  mPtRecoOverGen_E_30_50_Cent_50_80 = ibooker.book1D("PtRecoOverGen_E_30_50_Cent_50_80", "30<genpt<50; recopt/genpt (50-80%) (EndCap);counts", 90, 0, 2);
558  mPtRecoOverGen_F_30_50_Cent_50_80 = ibooker.book1D("PtRecoOverGen_F_30_50_Cent_50_80", "30<genpt<50; recopt/genpt (50-80%) (Forward);counts", 90, 0, 2);
559  mPtRecoOverGen_B_50_80_Cent_50_80 = ibooker.book1D("PtRecoOverGen_B_50_80_Cent_50_80", "50<genpt<80; recopt/genpt (50-80%) (Barrel);counts", 90, 0, 2);
560  mPtRecoOverGen_E_50_80_Cent_50_80 = ibooker.book1D("PtRecoOverGen_E_50_80_Cent_50_80", "50<genpt<80; recopt/genpt (50-80%) (EndCap);counts", 90, 0, 2);
561  mPtRecoOverGen_F_50_80_Cent_50_80 = ibooker.book1D("PtRecoOverGen_F_50_80_Cent_50_80", "50<genpt<80; recopt/genpt (50-80%) (Forward);counts", 90, 0, 2);
562  mPtRecoOverGen_B_80_120_Cent_50_80 = ibooker.book1D("PtRecoOverGen_B_80_120_Cent_50_80", "80<genpt<120; recopt/genpt (50-80%) (Barrel);counts", 90, 0, 2);
563  mPtRecoOverGen_E_80_120_Cent_50_80 = ibooker.book1D("PtRecoOverGen_E_80_120_Cent_50_80", "80<genpt<120; recopt/genpt (50-80%) (EndCap);counts", 90, 0, 2);
564  mPtRecoOverGen_F_80_120_Cent_50_80 = ibooker.book1D("PtRecoOverGen_F_80_120_Cent_50_80", "80<genpt<120; recopt/genpt (50-80%) (Forward);counts", 90, 0, 2);
565  mPtRecoOverGen_B_120_180_Cent_50_80 = ibooker.book1D("PtRecoOverGen_B_120_180_Cent_50_80", "120<genpt<180; recopt/genpt (50-80%) (Barrel);counts", 90, 0, 2);
566  mPtRecoOverGen_E_120_180_Cent_50_80 = ibooker.book1D("PtRecoOverGen_E_120_180_Cent_50_80", "120<genpt<180; recopt/genpt (50-80%) (EndCap);counts", 90, 0, 2);
567  mPtRecoOverGen_F_120_180_Cent_50_80 = ibooker.book1D("PtRecoOverGen_F_120_180_Cent_50_80", "120<genpt<180; recopt/genpt (50-80%) (Forward);counts", 90, 0, 2);
568  mPtRecoOverGen_B_180_300_Cent_50_80 = ibooker.book1D("PtRecoOverGen_B_180_300_Cent_50_80", "180<genpt<300; recopt/genpt (50-80%) (Barrel);counts", 90, 0, 2);
569  mPtRecoOverGen_E_180_300_Cent_50_80 = ibooker.book1D("PtRecoOverGen_E_180_300_Cent_50_80", "180<genpt<300; recopt/genpt (50-80%) (EndCap);counts", 90, 0, 2);
570  mPtRecoOverGen_F_180_300_Cent_50_80 = ibooker.book1D("PtRecoOverGen_F_180_300_Cent_50_80", "180<genpt<300; recopt/genpt (50-80%) (Forward);counts", 90, 0, 2);
571  mPtRecoOverGen_B_300_Inf_Cent_50_80 = ibooker.book1D("PtRecoOverGen_B_300_Inf_Cent_50_80", "300<genpt<Inf; recopt/genpt (50-80%) (Barrel);counts", 90, 0, 2);
572  mPtRecoOverGen_E_300_Inf_Cent_50_80 = ibooker.book1D("PtRecoOverGen_E_300_Inf_Cent_50_80", "300<genpt<Inf; recopt/genpt (50-80%) (EndCap);counts", 90, 0, 2);
573  mPtRecoOverGen_F_300_Inf_Cent_50_80 = ibooker.book1D("PtRecoOverGen_F_300_Inf_Cent_50_80", "300<genpt<Inf; recopt/genpt (50-80%) (Forward);counts", 90, 0, 2);
574 
575 
576  mPtRecoOverGen_GenPt_B_Cent_0_10 = ibooker.bookProfile("PtRecoOverGen_GenPt_B_Cent_0_10", Form("|#eta|<%2.2f, (0-10cent);genpt;recopt/genpt", BarrelEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
577  mPtRecoOverGen_GenPt_E_Cent_0_10 = ibooker.bookProfile("PtRecoOverGen_GenPt_E_Cent_0_10", Form("%2.2f<|#eta|<%2.2f, (0-10cent);genpt;recopt/genpt", BarrelEta, EndcapEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
578  mPtRecoOverGen_GenPt_F_Cent_0_10 = ibooker.bookProfile("PtRecoOverGen_GenPt_F_Cent_0_10", Form("%2.2f<|#eta|<%2.2f, (0-10cent);genpt;recopt/genpt", EndcapEta, ForwardEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
579  mPtRecoOverGen_GenPt_B_Cent_10_30 = ibooker.bookProfile("PtRecoOverGen_GenPt_B_Cent_10_30", Form("|#eta|<%2.2f, (10-30cent);genpt;recopt/genpt", BarrelEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
580  mPtRecoOverGen_GenPt_E_Cent_10_30 = ibooker.bookProfile("PtRecoOverGen_GenPt_E_Cent_10_30", Form("%2.2f<|#eta|<%2.2f, (10-30cent);genpt;recopt/genpt", BarrelEta, EndcapEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
581  mPtRecoOverGen_GenPt_F_Cent_10_30 = ibooker.bookProfile("PtRecoOverGen_GenPt_F_Cent_10_30", Form("%2.2f<|#eta|<%2.2f, (10-30cent);genpt;recopt/genpt", EndcapEta, ForwardEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
582  mPtRecoOverGen_GenPt_B_Cent_30_50 = ibooker.bookProfile("PtRecoOverGen_GenPt_B_Cent_30_50", Form("|#eta|<%2.2f, (30-50cent);genpt;recopt/genpt", BarrelEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
583  mPtRecoOverGen_GenPt_E_Cent_30_50 = ibooker.bookProfile("PtRecoOverGen_GenPt_E_Cent_30_50", Form("%2.2f<|#eta|<%2.2f, (30-50cent);genpt;recopt/genpt", BarrelEta, EndcapEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
584  mPtRecoOverGen_GenPt_F_Cent_30_50 = ibooker.bookProfile("PtRecoOverGen_GenPt_F_Cent_30_50", Form("%2.2f<|#eta|<%2.2f, (30-50cent);genpt;recopt/genpt", EndcapEta, ForwardEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
585  mPtRecoOverGen_GenPt_B_Cent_50_80 = ibooker.bookProfile("PtRecoOverGen_GenPt_B_Cent_50_80", Form("|#eta|<%2.2f, (50-80cent);genpt;recopt/genpt", BarrelEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
586  mPtRecoOverGen_GenPt_E_Cent_50_80 = ibooker.bookProfile("PtRecoOverGen_GenPt_E_Cent_50_80", Form("%2.2f<|#eta|<%2.2f, (50-80cent);genpt;recopt/genpt", BarrelEta, EndcapEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
587  mPtRecoOverGen_GenPt_F_Cent_50_80 = ibooker.bookProfile("PtRecoOverGen_GenPt_F_Cent_50_80", Form("%2.2f<|#eta|<%2.2f, (50-80cent);genpt;recopt/genpt", EndcapEta, ForwardEta), log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
588 
589  mPtRecoOverGen_GenEta_20_30_Cent_0_10 = ibooker.bookProfile("PtRecoOverGen_GenEta_20_30_Cent_0_10","20<genpt<30 (0-10%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
590  mPtRecoOverGen_GenEta_30_50_Cent_0_10 = ibooker.bookProfile("PtRecoOverGen_GenEta_30_50_Cent_0_10","30<genpt<50 (0-10%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
591  mPtRecoOverGen_GenEta_50_80_Cent_0_10 = ibooker.bookProfile("PtRecoOverGen_GenEta_50_80_Cent_0_10","50<genpt<80 (0-10%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
592  mPtRecoOverGen_GenEta_80_120_Cent_0_10 = ibooker.bookProfile("PtRecoOverGen_GenEta_80_120_Cent_0_10","80<genpt<120 (0-10%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
593  mPtRecoOverGen_GenEta_120_180_Cent_0_10 = ibooker.bookProfile("PtRecoOverGen_GenEta_120_180_Cent_0_10","120<genpt<180 (0-10%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
594  mPtRecoOverGen_GenEta_180_300_Cent_0_10 = ibooker.bookProfile("PtRecoOverGen_GenEta_180_300_Cent_0_10","180<genpt<300 (0-10%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
595  mPtRecoOverGen_GenEta_300_Inf_Cent_0_10 = ibooker.bookProfile("PtRecoOverGen_GenEta_300_Inf_Cent_0_10","300<genpt<Inf (0-10%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
596 
597 
598  mPtRecoOverGen_GenEta_20_30_Cent_10_30 = ibooker.bookProfile("PtRecoOverGen_GenEta_20_30_Cent_10_30","20<genpt<30 (10-30%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
599  mPtRecoOverGen_GenEta_30_50_Cent_10_30 = ibooker.bookProfile("PtRecoOverGen_GenEta_30_50_Cent_10_30","30<genpt<50 (10-30%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
600  mPtRecoOverGen_GenEta_50_80_Cent_10_30 = ibooker.bookProfile("PtRecoOverGen_GenEta_50_80_Cent_10_30","50<genpt<80 (10-30%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
601  mPtRecoOverGen_GenEta_80_120_Cent_10_30 = ibooker.bookProfile("PtRecoOverGen_GenEta_80_120_Cent_10_30","80<genpt<120 (10-30%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
602  mPtRecoOverGen_GenEta_120_180_Cent_10_30 = ibooker.bookProfile("PtRecoOverGen_GenEta_120_180_Cent_10_30","120<genpt<180 (10-30%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
603  mPtRecoOverGen_GenEta_180_300_Cent_10_30 = ibooker.bookProfile("PtRecoOverGen_GenEta_180_300_Cent_10_30","180<genpt<300 (10-30%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
604  mPtRecoOverGen_GenEta_300_Inf_Cent_10_30 = ibooker.bookProfile("PtRecoOverGen_GenEta_300_Inf_Cent_10_30","300<genpt<Inf (10-30%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
605 
606  mPtRecoOverGen_GenEta_20_30_Cent_30_50 = ibooker.bookProfile("PtRecoOverGen_GenEta_20_30_Cent_30_50","20<genpt<30 (30-50%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
607  mPtRecoOverGen_GenEta_30_50_Cent_30_50 = ibooker.bookProfile("PtRecoOverGen_GenEta_30_50_Cent_30_50","30<genpt<50 (30-50%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
608  mPtRecoOverGen_GenEta_50_80_Cent_30_50 = ibooker.bookProfile("PtRecoOverGen_GenEta_50_80_Cent_30_50","50<genpt<80 (30-50%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
609  mPtRecoOverGen_GenEta_80_120_Cent_30_50 = ibooker.bookProfile("PtRecoOverGen_GenEta_80_120_Cent_30_50","80<genpt<120 (30-50%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
610  mPtRecoOverGen_GenEta_120_180_Cent_30_50 = ibooker.bookProfile("PtRecoOverGen_GenEta_120_180_Cent_30_50","120<genpt<180 (30-50%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
611  mPtRecoOverGen_GenEta_180_300_Cent_30_50 = ibooker.bookProfile("PtRecoOverGen_GenEta_180_300_Cent_30_50","180<genpt<300 (30-50%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
612  mPtRecoOverGen_GenEta_300_Inf_Cent_30_50 = ibooker.bookProfile("PtRecoOverGen_GenEta_300_Inf_Cent_30_50","300<genpt<Inf (30-50%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
613 
614  mPtRecoOverGen_GenEta_20_30_Cent_50_80 = ibooker.bookProfile("PtRecoOverGen_GenEta_20_30_Cent_50_80","20<genpt<30 (50-80%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
615  mPtRecoOverGen_GenEta_30_50_Cent_50_80 = ibooker.bookProfile("PtRecoOverGen_GenEta_30_50_Cent_50_80","30<genpt<50 (50-80%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
616  mPtRecoOverGen_GenEta_50_80_Cent_50_80 = ibooker.bookProfile("PtRecoOverGen_GenEta_50_80_Cent_50_80","50<genpt<80 (50-80%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
617  mPtRecoOverGen_GenEta_80_120_Cent_50_80 = ibooker.bookProfile("PtRecoOverGen_GenEta_80_120_Cent_50_80","80<genpt<120 (50-80%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
618  mPtRecoOverGen_GenEta_120_180_Cent_50_80 = ibooker.bookProfile("PtRecoOverGen_GenEta_120_180_Cent_50_80","120<genpt<180 (50-80%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
619  mPtRecoOverGen_GenEta_180_300_Cent_50_80 = ibooker.bookProfile("PtRecoOverGen_GenEta_180_300_Cent_50_80","180<genpt<300 (50-80%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
620  mPtRecoOverGen_GenEta_300_Inf_Cent_50_80 = ibooker.bookProfile("PtRecoOverGen_GenEta_300_Inf_Cent_50_80","300<genpt<Inf (50-80%);geneta;recopt/genpt",90, etaRange, 0, 2, " ");
621 
622 
623  if (mOutputFile.empty ())
624  LogInfo("OutputInfo") << " Histograms will NOT be saved";
625  else
626  LogInfo("OutputInfo") << " Histograms will be saved to file:" << mOutputFile;
627 
628  delete h2D_etabins_vs_pt2;
629  delete h2D_etabins_vs_pt;
630  delete h2D_etabins_vs_phi;
631  delete h2D_pfcand_etabins_vs_pt;
632 }
633 
634 
635 
636 //------------------------------------------------------------------------------
637 // ~JetTester_HeavyIons
638 //------------------------------------------------------------------------------
639 JetTester_HeavyIons::~JetTester_HeavyIons() {}
640 
641 
642 //------------------------------------------------------------------------------
643 // beginJob
644 //------------------------------------------------------------------------------
645 //void JetTester_HeavyIons::beginJob() {
646 // //std::cout<<"inside the begin job function"<<endl;
647 //}
648 
649 
650 //------------------------------------------------------------------------------
651 // endJob
652 //------------------------------------------------------------------------------
653 //void JetTester_HeavyIons::endJob()
654 //{
655 // if (!mOutputFile.empty() && &*edm::Service<DQMStore>())
656 // {
657 // edm::Service<DQMStore>()->save(mOutputFile);
658 // }
659 //}
660 
661 
662 //------------------------------------------------------------------------------
663 // analyze
664 //------------------------------------------------------------------------------
665 void JetTester_HeavyIons::analyze(const edm::Event& mEvent, const edm::EventSetup& mSetup)
666 {
667  // Get the primary vertices
668  //----------------------------------------------------------------------------
669  edm::Handle<vector<reco::Vertex> > pvHandle;
670  mEvent.getByToken(pvToken_, pvHandle);
671  reco::Vertex::Point vtx(0,0,0);
672  edm::Handle<reco::VertexCollection> vtxs;
673  //vtx = getVtx(mEvent);
674 
675  mEvent.getByToken(hiVertexToken_, vtxs);
676  int greatestvtx = 0;
677  int nVertex = vtxs->size();
678 
679  for (unsigned int i = 0 ; i< vtxs->size(); ++i){
680  unsigned int daughter = (*vtxs)[i].tracksSize();
681  if( daughter > (*vtxs)[greatestvtx].tracksSize()) greatestvtx = i;
682  }
683 
684  if(nVertex<=0){
685  vtx = reco::Vertex::Point(0,0,0);
686  }
687  vtx = (*vtxs)[greatestvtx].position();
688 
689  int nGoodVertices = 0;
690 
691  if (pvHandle.isValid())
692  {
693  for (unsigned i=0; i<pvHandle->size(); i++)
694  {
695  if ((*pvHandle)[i].ndof() > 4 &&
696  (fabs((*pvHandle)[i].z()) <= 24) &&
697  (fabs((*pvHandle)[i].position().rho()) <= 2))
698  nGoodVertices++;
699  }
700  }
701 
702  mNvtx->Fill(nGoodVertices);
703 
704  // Get the Jet collection
705  //----------------------------------------------------------------------------
706  //math::XYZTLorentzVector p4tmp[2];
707 
708  std::vector<Jet> recoJets;
709  recoJets.clear();
710 
711  edm::Handle<CaloJetCollection> caloJets;
712  edm::Handle<JPTJetCollection> jptJets;
713  edm::Handle<PFJetCollection> pfJets;
714  edm::Handle<BasicJetCollection> basicJets;
715 
716  // Get the Particle flow candidates and the Voronoi variables
717  edm::Handle<reco::PFCandidateCollection> pfCandidates;
718  edm::Handle<CaloTowerCollection> caloCandidates;
719  edm::Handle<reco::CandidateView> pfcandidates_;
720  edm::Handle<reco::CandidateView> calocandidates_;
721 
722 
723  // get the centrality
724  edm::Handle<reco::Centrality> cent;
725  mEvent.getByToken(centralityToken, cent); //_centralitytag comes from the cfg
726 
727  mHF->Fill(cent->EtHFtowerSum());
728  Float_t HF_energy = cent->EtHFtowerSum();
729 
730  edm::Handle<int> cbin;
731  mEvent.getByToken(centralityBinToken, cbin);
732  if (!cent.isValid()) return;
733 
734  int hibin = -999;
735  if(cbin.isValid()){
736  hibin = *cbin;
737  }
738  //else edm::LogWarning("JetTester_HeavyIons") << "invalid collection: centralityBin " << std::endl;
739 
740  bool isCentral = false;
741  bool ismidCentral = false;
742  bool ismidPeripheral = false;
743  bool isPeripheral = false;
744 
745  if(hibin < 20) isCentral = true;
746  if(hibin >= 20 && hibin < 60) ismidCentral = true;
747  if(hibin >= 60 && hibin < 100) ismidPeripheral = true;
748  if(hibin >= 100 && hibin < 160) isPeripheral = true;
749 
750 
751  if (isCaloJet) mEvent.getByToken(caloJetsToken_, caloJets);
752  if (isJPTJet) mEvent.getByToken(jptJetsToken_, jptJets);
753  if (isPFJet) {
754  if(std::string("Pu")==UEAlgo) mEvent.getByToken(basicJetsToken_, basicJets);
755  }
756 
757  mEvent.getByToken(pfCandToken_, pfCandidates);
758  mEvent.getByToken(pfCandViewToken_, pfcandidates_);
759 
760  mEvent.getByToken(caloTowersToken_, caloCandidates);
761  mEvent.getByToken(caloCandViewToken_, calocandidates_);
762 
763 
764 
765  const reco::PFCandidateCollection *pfCandidateColl = pfCandidates.product();
766 
767  Int_t NPFpart = 0;
768  Int_t NCaloTower = 0;
769  Float_t pfPt = 0;
770  Float_t pfEta = 0;
771  Int_t pfID = 0;
772  Float_t pfPhi = 0;
773  Float_t caloPt = 0;
774  Float_t caloEta = 0;
775  Float_t caloPhi = 0;
776  Float_t SumPt_value = 0;
777 
778  double edge_pseudorapidity[etaBins_ +1] = {-5.191, -2.650, -2.043, -1.740, -1.479, -1.131, -0.783, -0.522, 0.522, 0.783, 1.131, 1.479, 1.740, 2.043, 2.650, 5.191 };
779 
780 
781  if(isCaloJet){
782 
783  Float_t SumCaloPt[etaBins_];
784  Float_t SumSquaredCaloPt[etaBins_];
785 
786  for(int i = 0;i<etaBins_;i++){
787 
788  SumCaloPt[i] = 0;
789  SumSquaredCaloPt[i] = 0;
790 
791  }
792 
793  for(unsigned icand = 0;icand<caloCandidates->size(); icand++){
794 
795  const CaloTower & tower = (*caloCandidates)[icand];
796  reco::CandidateViewRef ref(calocandidates_,icand);
797  if(tower.p4(vtx).Et() < 0.1) continue;
798 
799 
800  NCaloTower++;
801 
802  caloPt = tower.p4(vtx).Et();
803  caloEta = tower.p4(vtx).Eta();
804  caloPhi = tower.p4(vtx).Phi();
805 
806 
807 
808  for(size_t k = 0;k<nedge_pseudorapidity-1; k++){
809  if(caloEta >= edge_pseudorapidity[k] && caloEta < edge_pseudorapidity[k+1]){
810  SumCaloPt[k] = SumCaloPt[k] + caloPt;
811  SumSquaredCaloPt[k] = SumSquaredCaloPt[k] + caloPt*caloPt;
812  }// eta selection statement
813 
814  }// eta bin loop
815 
816  SumPt_value = SumPt_value + caloPt;
817 
818  mCaloPt->Fill(caloPt);
819  mCaloEta->Fill(caloEta);
820  mCaloPhi->Fill(caloPhi);
821 
822  }// calo tower candidate loop
823 
824  Float_t Evt_SumCaloPt = 0;
825 
826  Float_t Evt_SumSquaredCaloPt = 0;
827 
828 
829  mSumCaloPt_n5p191_n2p650->Fill(SumCaloPt[0]);
830  mSumCaloPt_n2p650_n2p043->Fill(SumCaloPt[1]);
831  mSumCaloPt_n2p043_n1p740->Fill(SumCaloPt[2]);
832  mSumCaloPt_n1p740_n1p479->Fill(SumCaloPt[3]);
833  mSumCaloPt_n1p479_n1p131->Fill(SumCaloPt[4]);
834  mSumCaloPt_n1p131_n0p783->Fill(SumCaloPt[5]);
835  mSumCaloPt_n0p783_n0p522->Fill(SumCaloPt[6]);
836  mSumCaloPt_n0p522_0p522->Fill(SumCaloPt[7]);
837  mSumCaloPt_0p522_0p783->Fill(SumCaloPt[8]);
838  mSumCaloPt_0p783_1p131->Fill(SumCaloPt[9]);
839  mSumCaloPt_1p131_1p479->Fill(SumCaloPt[10]);
840  mSumCaloPt_1p479_1p740->Fill(SumCaloPt[11]);
841  mSumCaloPt_1p740_2p043->Fill(SumCaloPt[12]);
842  mSumCaloPt_2p043_2p650->Fill(SumCaloPt[13]);
843  mSumCaloPt_2p650_5p191->Fill(SumCaloPt[14]);
844 
845  for(size_t k = 0;k<nedge_pseudorapidity-1;k++){
846 
847  Evt_SumCaloPt = Evt_SumCaloPt + SumCaloPt[k];
848 
849  Evt_SumSquaredCaloPt = Evt_SumSquaredCaloPt + SumSquaredCaloPt[k];
850 
851  }// eta bin loop
852 
853 
854  mSumCaloPt->Fill(Evt_SumCaloPt);
855 
856  mSumSquaredCaloPt->Fill(Evt_SumSquaredCaloPt);
857 
858  mSumCaloPt_HF->Fill(Evt_SumCaloPt,HF_energy);
859 
860  mNCalopart->Fill(NCaloTower);
861  mSumpt->Fill(SumPt_value);
862 
863  }// is calo jet
864 
865  if(isPFJet){
866 
867  Float_t SumPFPt[etaBins_];
868 
869  Float_t SumSquaredPFPt[etaBins_];
870 
871  for(int i = 0;i<etaBins_;i++){
872 
873  SumPFPt[i] = 0;
874  SumSquaredPFPt[i] = 0;
875 
876  }
877 
878  for(unsigned icand=0;icand<pfCandidateColl->size(); icand++){
879 
880  const reco::PFCandidate pfCandidate = pfCandidateColl->at(icand);
881  reco::CandidateViewRef ref(pfcandidates_,icand);
882 
883  if(pfCandidate.pt() < 0.5) continue;
884 
885 
886  NPFpart++;
887  pfPt = pfCandidate.pt();
888  pfEta = pfCandidate.eta();
889  pfPhi = pfCandidate.phi();
890  pfID = pfCandidate.particleId();
891 
892  bool inBarrel = false;
893  bool inEndcap = false;
894  bool inForward = false;
895 
896  if(fabs(pfEta)<BarrelEta) inBarrel = true;
897  if(fabs(pfEta)>=BarrelEta && fabs(pfEta)<EndcapEta) inEndcap = true;
898  if(fabs(pfEta)>=EndcapEta && fabs(pfEta)<ForwardEta) inForward = true;
899 
900  switch(pfID){
901  case 0 :
902  mPFCandpT_vs_eta_Unknown->Fill(pfPt, pfEta);
903  if(inBarrel) mPFCandpT_Barrel_Unknown->Fill(pfPt);
904  if(inEndcap) mPFCandpT_Endcap_Unknown->Fill(pfPt);
905  if(inForward) mPFCandpT_Forward_Unknown->Fill(pfPt);
906  case 1 :
907  mPFCandpT_vs_eta_ChargedHadron->Fill(pfPt, pfEta);
908  if(inBarrel) mPFCandpT_Barrel_ChargedHadron->Fill(pfPt);
909  if(inEndcap) mPFCandpT_Endcap_ChargedHadron->Fill(pfPt);
910  if(inForward) mPFCandpT_Forward_ChargedHadron->Fill(pfPt);
911  case 2 :
912  mPFCandpT_vs_eta_electron->Fill(pfPt, pfEta);
913  if(inBarrel) mPFCandpT_Barrel_electron->Fill(pfPt);
914  if(inEndcap) mPFCandpT_Endcap_electron->Fill(pfPt);
915  if(inForward) mPFCandpT_Forward_electron->Fill(pfPt);
916  case 3 :
917  mPFCandpT_vs_eta_muon->Fill(pfPt, pfEta);
918  if(inBarrel) mPFCandpT_Barrel_muon->Fill(pfPt);
919  if(inEndcap) mPFCandpT_Endcap_muon->Fill(pfPt);
920  if(inForward) mPFCandpT_Forward_muon->Fill(pfPt);
921  case 4 :
922  mPFCandpT_vs_eta_photon->Fill(pfPt, pfEta);
923  if(inBarrel) mPFCandpT_Barrel_photon->Fill(pfPt);
924  if(inEndcap) mPFCandpT_Endcap_photon->Fill(pfPt);
925  if(inForward) mPFCandpT_Forward_photon->Fill(pfPt);
926  case 5 :
927  mPFCandpT_vs_eta_NeutralHadron->Fill(pfPt, pfEta);
928  if(inBarrel) mPFCandpT_Barrel_NeutralHadron->Fill(pfPt);
929  if(inEndcap) mPFCandpT_Endcap_NeutralHadron->Fill(pfPt);
930  if(inForward) mPFCandpT_Forward_NeutralHadron->Fill(pfPt);
931  case 6 :
932  mPFCandpT_vs_eta_HadE_inHF->Fill(pfPt, pfEta);
933  if(inBarrel) mPFCandpT_Barrel_HadE_inHF->Fill(pfPt);
934  if(inEndcap) mPFCandpT_Endcap_HadE_inHF->Fill(pfPt);
935  if(inForward) mPFCandpT_Forward_HadE_inHF->Fill(pfPt);
936  case 7 :
937  mPFCandpT_vs_eta_EME_inHF->Fill(pfPt, pfEta);
938  if(inBarrel) mPFCandpT_Barrel_EME_inHF->Fill(pfPt);
939  if(inEndcap) mPFCandpT_Endcap_EME_inHF->Fill(pfPt);
940  if(inForward) mPFCandpT_Forward_EME_inHF->Fill(pfPt);
941  }
942 
943 
944 
945 
946  for(size_t k = 0;k<nedge_pseudorapidity-1; k++){
947  if(pfEta >= edge_pseudorapidity[k] && pfEta < edge_pseudorapidity[k+1]){
948  SumPFPt[k] = SumPFPt[k] + pfPt;
949 
950  SumSquaredPFPt[k] = SumSquaredPFPt[k] + pfPt*pfPt;
951 
952  }// eta selection statement
953 
954  }// eta bin loop
955 
956  SumPt_value = SumPt_value + pfPt;
957 
958  mPFPt->Fill(pfPt);
959  mPFEta->Fill(pfEta);
960  mPFPhi->Fill(pfPhi);
961 
962 
963  }// pf candidate loop
964 
965  Float_t Evt_SumPFPt = 0;
966 
967  Float_t Evt_SumSquaredPFPt = 0;
968 
969 
970  mSumPFPt_n5p191_n2p650->Fill(SumPFPt[0]);
971  mSumPFPt_n2p650_n2p043->Fill(SumPFPt[1]);
972  mSumPFPt_n2p043_n1p740->Fill(SumPFPt[2]);
973  mSumPFPt_n1p740_n1p479->Fill(SumPFPt[3]);
974  mSumPFPt_n1p479_n1p131->Fill(SumPFPt[4]);
975  mSumPFPt_n1p131_n0p783->Fill(SumPFPt[5]);
976  mSumPFPt_n0p783_n0p522->Fill(SumPFPt[6]);
977  mSumPFPt_n0p522_0p522->Fill(SumPFPt[7]);
978  mSumPFPt_0p522_0p783->Fill(SumPFPt[8]);
979  mSumPFPt_0p783_1p131->Fill(SumPFPt[9]);
980  mSumPFPt_1p131_1p479->Fill(SumPFPt[10]);
981  mSumPFPt_1p479_1p740->Fill(SumPFPt[11]);
982  mSumPFPt_1p740_2p043->Fill(SumPFPt[12]);
983  mSumPFPt_2p043_2p650->Fill(SumPFPt[13]);
984  mSumPFPt_2p650_5p191->Fill(SumPFPt[14]);
985 
986  for(size_t k = 0;k<nedge_pseudorapidity-1;k++){
987 
988  Evt_SumPFPt = Evt_SumPFPt + SumPFPt[k];
989 
990  Evt_SumSquaredPFPt = Evt_SumSquaredPFPt + SumSquaredPFPt[k];
991 
992  }// eta bin loop
993 
994  mSumPFPt->Fill(Evt_SumPFPt);
995 
996  mSumSquaredPFPt->Fill(Evt_SumSquaredPFPt);
997 
998  mSumPFPt_HF->Fill(Evt_SumPFPt,HF_energy);
999 
1000  mNPFpart->Fill(NPFpart);
1001  mSumpt->Fill(SumPt_value);
1002 
1003  }
1004 
1005  if (isCaloJet)
1006  {
1007  for (unsigned ijet=0; ijet<caloJets->size(); ijet++) recoJets.push_back((*caloJets)[ijet]);
1008  }
1009 
1010  if (isJPTJet)
1011  {
1012  for (unsigned ijet=0; ijet<jptJets->size(); ijet++) recoJets.push_back((*jptJets)[ijet]);
1013  }
1014 
1015  if (isPFJet) {
1016  if(std::string("Pu")==UEAlgo){
1017  for (unsigned ijet=0; ijet<basicJets->size();ijet++) recoJets.push_back((*basicJets)[ijet]);
1018  }
1019  }
1020 
1021  if (isCaloJet && !caloJets.isValid()) return;
1022  if (isJPTJet && !jptJets.isValid()) return;
1023  if (isPFJet){
1024  if(std::string("Pu")==UEAlgo){if(!basicJets.isValid()) return;}
1025  }
1026 
1027  int nJet_40 = 0;
1028 
1029  mNJets->Fill(recoJets.size());
1030 
1031  for (unsigned ijet=0; ijet<recoJets.size(); ijet++) {
1032 
1033  if (recoJets[ijet].pt() > mRecoJetPtThreshold) {
1034  //counting forward and barrel jets
1035  // get an idea of no of jets with pT>40 GeV
1036  if(recoJets[ijet].pt() > 40)
1037  nJet_40++;
1038 
1039  if (mEta) mEta->Fill(recoJets[ijet].eta());
1040  if (mjetpileup) mjetpileup->Fill(recoJets[ijet].pileup());
1041  if (mJetArea) mJetArea ->Fill(recoJets[ijet].jetArea());
1042  if (mPhi) mPhi ->Fill(recoJets[ijet].phi());
1043  if (mEnergy) mEnergy ->Fill(recoJets[ijet].energy());
1044  if (mP) mP ->Fill(recoJets[ijet].p());
1045  if (mPt) mPt ->Fill(recoJets[ijet].pt());
1046  if (mMass) mMass ->Fill(recoJets[ijet].mass());
1047  if (mConstituents) mConstituents->Fill(recoJets[ijet].nConstituents());
1048 
1049  }
1050  }
1051 
1052  if (mNJets_40) mNJets_40->Fill(nJet_40);
1053 
1054 
1055  // Gen level information:
1056  if (!mEvent.isRealData()){
1057 
1058  // Get ptHat
1059  //------------------------------------------------------------------------
1060  edm::Handle<GenEventInfoProduct> myGenEvt;
1061  mEvent.getByToken(evtToken_, myGenEvt);
1062 
1063  if (myGenEvt.isValid()) {
1064  if(myGenEvt->hasBinningValues()){
1065  double ptHat = myGenEvt->binningValues()[0];
1066  if (mPtHat) mPtHat->Fill(ptHat);
1067  }
1068  }
1069  // Gen jets
1070  //------------------------------------------------------------------------
1071  edm::Handle<GenJetCollection> genJets;
1072  mEvent.getByToken(genJetsToken_, genJets);
1073 
1074  if (!genJets.isValid()) return;
1075 
1076  for (GenJetCollection::const_iterator gjet=genJets->begin(); gjet!=genJets->end(); gjet++) {
1077  if(gjet->pt() > mMatchGenPtThreshold){
1078  if (mGenEta) mGenEta->Fill(gjet->eta());
1079  if (mGenPhi) mGenPhi->Fill(gjet->phi());
1080  if (mGenPt) mGenPt ->Fill(gjet->pt());
1081  }
1082  }
1083 
1084  if (!(mInputGenCollection.label().empty())) {
1085  for (GenJetCollection::const_iterator gjet=genJets->begin(); gjet!=genJets->end(); gjet++) {
1086  if (fabs(gjet->eta()) > 6.) continue; // Out of the detector
1087  if (gjet->pt() < mMatchGenPtThreshold) continue;
1088  if (recoJets.size() <= 0) continue;
1089 
1090  bool inBarrel = false;
1091  bool inEndcap = false;
1092  bool inForward = false;
1093 
1094  if(fabs(gjet->eta())<BarrelEta) inBarrel = true;
1095  if(fabs(gjet->eta())>=BarrelEta && fabs(gjet->eta())<EndcapEta) inEndcap = true;
1096  if(fabs(gjet->eta())>=EndcapEta && fabs(gjet->eta())<ForwardEta) inForward = true;
1097 
1098  // pt response
1099  //------------------------------------------------------------
1100  int iMatch = -1;
1101  double deltaRBest = 999;
1102  double JetPtBest = 0;
1103  for (unsigned ijet=0; ijet<recoJets.size(); ++ijet) {
1104  double recoPt = recoJets[ijet].pt();
1105  if (recoPt > 10) {
1106  double delR = deltaR(gjet->eta(), gjet->phi(), recoJets[ijet].eta(), recoJets[ijet].phi());
1107  if (delR < deltaRBest) {
1108  deltaRBest = delR;
1109  JetPtBest = recoPt;
1110  iMatch = ijet;
1111  }
1112  }
1113  }
1114  if (iMatch<0) continue;
1115 
1116  //fillMatchHists(gjet->eta(), gjet->phi(), gjet->pt(), recoJets[iMatch].eta(), recoJets[iMatch].phi(), recoJets[iMatch].pt(), hibin);
1117  if(deltaRBest < mRThreshold){
1118  double genpt = gjet->pt();
1119  double geneta = gjet->eta();
1120  double response = JetPtBest / genpt;
1121  // Fill all the response histograms here: for each pT bin, eta region, centrality bin
1122 
1123  if(inBarrel) {
1124  if(isCentral)
1125  mPtRecoOverGen_GenPt_B_Cent_0_10->Fill(log10(genpt), response);
1126  if(ismidCentral)
1127  mPtRecoOverGen_GenPt_B_Cent_10_30->Fill(log10(genpt), response);
1128  if(ismidPeripheral)
1129  mPtRecoOverGen_GenPt_B_Cent_30_50->Fill(log10(genpt), response);
1130  if(isPeripheral)
1131  mPtRecoOverGen_GenPt_B_Cent_50_80->Fill(log10(genpt), response);
1132  }
1133  if(inEndcap) {
1134  if(isCentral)
1135  mPtRecoOverGen_GenPt_E_Cent_0_10->Fill(log10(genpt), response);
1136  if(ismidCentral)
1137  mPtRecoOverGen_GenPt_E_Cent_10_30->Fill(log10(genpt), response);
1138  if(ismidPeripheral)
1139  mPtRecoOverGen_GenPt_E_Cent_30_50->Fill(log10(genpt), response);
1140  if(isPeripheral)
1141  mPtRecoOverGen_GenPt_E_Cent_50_80->Fill(log10(genpt), response);
1142  }
1143  if(inForward) {
1144  if(isCentral)
1145  mPtRecoOverGen_GenPt_F_Cent_0_10->Fill(log10(genpt), response);
1146  if(ismidCentral)
1147  mPtRecoOverGen_GenPt_F_Cent_10_30->Fill(log10(genpt), response);
1148  if(ismidPeripheral)
1149  mPtRecoOverGen_GenPt_F_Cent_30_50->Fill(log10(genpt), response);
1150  if(isPeripheral)
1151  mPtRecoOverGen_GenPt_F_Cent_50_80->Fill(log10(genpt), response);
1152  }
1153 
1154  if( gjet->pt() >= 20 && gjet->pt() < 30 ){
1155  if(isCentral){
1156  mPtRecoOverGen_GenEta_20_30_Cent_0_10->Fill(geneta, response);
1157  if(inBarrel)
1158  mPtRecoOverGen_B_20_30_Cent_0_10->Fill(response);
1159  if(inEndcap)
1160  mPtRecoOverGen_E_20_30_Cent_0_10->Fill(response);
1161  if(inForward)
1162  mPtRecoOverGen_F_20_30_Cent_0_10->Fill(response);
1163  }//
1164  if(ismidCentral){
1165  mPtRecoOverGen_GenEta_20_30_Cent_10_30->Fill(geneta, response);
1166  if(inBarrel)
1167  mPtRecoOverGen_B_20_30_Cent_10_30->Fill(response);
1168  if(inEndcap)
1169  mPtRecoOverGen_E_20_30_Cent_10_30->Fill(response);
1170  if(inForward)
1171  mPtRecoOverGen_F_20_30_Cent_10_30->Fill(response);
1172  }//
1173  if(ismidPeripheral){
1174  mPtRecoOverGen_GenEta_20_30_Cent_30_50->Fill(geneta, response);
1175  if(inBarrel)
1176  mPtRecoOverGen_B_20_30_Cent_30_50->Fill(response);
1177  if(inEndcap)
1178  mPtRecoOverGen_E_20_30_Cent_30_50->Fill(response);
1179  if(inForward)
1180  mPtRecoOverGen_F_20_30_Cent_30_50->Fill(response);
1181  }//
1182  if(isPeripheral){
1183  mPtRecoOverGen_GenEta_20_30_Cent_50_80->Fill(geneta, response);
1184  if(inBarrel)
1185  mPtRecoOverGen_B_20_30_Cent_50_80->Fill(response);
1186  if(inEndcap)
1187  mPtRecoOverGen_E_20_30_Cent_50_80->Fill(response);
1188  if(inForward)
1189  mPtRecoOverGen_F_20_30_Cent_50_80->Fill(response);
1190  }//
1191  }// pt bin 20-30
1192 
1193  if( gjet->pt() >= 30 && gjet->pt() < 50 ){
1194  if(isCentral){
1195  mPtRecoOverGen_GenEta_30_50_Cent_0_10->Fill(geneta, response);
1196  if(inBarrel)
1197  mPtRecoOverGen_B_30_50_Cent_0_10->Fill(response);
1198  if(inEndcap)
1199  mPtRecoOverGen_E_30_50_Cent_0_10->Fill(response);
1200  if(inForward)
1201  mPtRecoOverGen_F_30_50_Cent_0_10->Fill(response);
1202  }//
1203  if(ismidCentral){
1204  mPtRecoOverGen_GenEta_30_50_Cent_10_30->Fill(geneta, response);
1205  if(inBarrel)
1206  mPtRecoOverGen_B_30_50_Cent_10_30->Fill(response);
1207  if(inEndcap)
1208  mPtRecoOverGen_E_30_50_Cent_10_30->Fill(response);
1209  if(inForward)
1210  mPtRecoOverGen_F_30_50_Cent_10_30->Fill(response);
1211  }//
1212  if(ismidPeripheral){
1213  mPtRecoOverGen_GenEta_30_50_Cent_30_50->Fill(geneta, response);
1214  if(inBarrel)
1215  mPtRecoOverGen_B_30_50_Cent_30_50->Fill(response);
1216  if(inEndcap)
1217  mPtRecoOverGen_E_30_50_Cent_30_50->Fill(response);
1218  if(inForward)
1219  mPtRecoOverGen_F_30_50_Cent_30_50->Fill(response);
1220  }//
1221  if(isPeripheral){
1222  mPtRecoOverGen_GenEta_30_50_Cent_50_80->Fill(geneta, response);
1223  if(inBarrel)
1224  mPtRecoOverGen_B_30_50_Cent_50_80->Fill(response);
1225  if(inEndcap)
1226  mPtRecoOverGen_E_30_50_Cent_50_80->Fill(response);
1227  if(inForward)
1228  mPtRecoOverGen_F_30_50_Cent_50_80->Fill(response);
1229  }//
1230  }// pt bin 30-50
1231 
1232  if( gjet->pt() >= 50 && gjet->pt() < 80 ){
1233  if(isCentral){
1234  mPtRecoOverGen_GenEta_50_80_Cent_0_10->Fill(geneta, response);
1235  if(inBarrel)
1236  mPtRecoOverGen_B_50_80_Cent_0_10->Fill(response);
1237  if(inEndcap)
1238  mPtRecoOverGen_E_50_80_Cent_0_10->Fill(response);
1239  if(inForward)
1240  mPtRecoOverGen_F_50_80_Cent_0_10->Fill(response);
1241  }//
1242  if(ismidCentral){
1243  mPtRecoOverGen_GenEta_50_80_Cent_10_30->Fill(geneta, response);
1244  if(inBarrel)
1245  mPtRecoOverGen_B_50_80_Cent_10_30->Fill(response);
1246  if(inEndcap)
1247  mPtRecoOverGen_E_50_80_Cent_10_30->Fill(response);
1248  if(inForward)
1249  mPtRecoOverGen_F_50_80_Cent_10_30->Fill(response);
1250  }//
1251  if(ismidPeripheral){
1252  mPtRecoOverGen_GenEta_50_80_Cent_30_50->Fill(geneta, response);
1253  if(inBarrel)
1254  mPtRecoOverGen_B_50_80_Cent_30_50->Fill(response);
1255  if(inEndcap)
1256  mPtRecoOverGen_E_50_80_Cent_30_50->Fill(response);
1257  if(inForward)
1258  mPtRecoOverGen_F_50_80_Cent_30_50->Fill(response);
1259  }//
1260  if(isPeripheral){
1261  mPtRecoOverGen_GenEta_50_80_Cent_50_80->Fill(geneta, response);
1262  if(inBarrel)
1263  mPtRecoOverGen_B_50_80_Cent_50_80->Fill(response);
1264  if(inEndcap)
1265  mPtRecoOverGen_E_50_80_Cent_50_80->Fill(response);
1266  if(inForward)
1267  mPtRecoOverGen_F_50_80_Cent_50_80->Fill(response);
1268  }//
1269  }// pt bin 50-80
1270 
1271  if( gjet->pt() >= 80 && gjet->pt() < 120 ){
1272  if(isCentral){
1273  mPtRecoOverGen_GenEta_80_120_Cent_0_10->Fill(geneta, response);
1274  if(inBarrel)
1275  mPtRecoOverGen_B_80_120_Cent_0_10->Fill(response);
1276  if(inEndcap)
1277  mPtRecoOverGen_E_80_120_Cent_0_10->Fill(response);
1278  if(inForward)
1279  mPtRecoOverGen_F_80_120_Cent_0_10->Fill(response);
1280  }//
1281  if(ismidCentral){
1282  mPtRecoOverGen_GenEta_80_120_Cent_10_30->Fill(geneta, response);
1283  if(inBarrel)
1284  mPtRecoOverGen_B_80_120_Cent_10_30->Fill(response);
1285  if(inEndcap)
1286  mPtRecoOverGen_E_80_120_Cent_10_30->Fill(response);
1287  if(inForward)
1288  mPtRecoOverGen_F_80_120_Cent_10_30->Fill(response);
1289  }//
1290  if(ismidPeripheral){
1291  mPtRecoOverGen_GenEta_80_120_Cent_30_50->Fill(geneta, response);
1292  if(inBarrel)
1293  mPtRecoOverGen_B_80_120_Cent_30_50->Fill(response);
1294  if(inEndcap)
1295  mPtRecoOverGen_E_80_120_Cent_30_50->Fill(response);
1296  if(inForward)
1297  mPtRecoOverGen_F_80_120_Cent_30_50->Fill(response);
1298  }//
1299  if(isPeripheral){
1300  mPtRecoOverGen_GenEta_80_120_Cent_50_80->Fill(geneta, response);
1301  if(inBarrel)
1302  mPtRecoOverGen_B_80_120_Cent_50_80->Fill(response);
1303  if(inEndcap)
1304  mPtRecoOverGen_E_80_120_Cent_50_80->Fill(response);
1305  if(inForward)
1306  mPtRecoOverGen_F_80_120_Cent_50_80->Fill(response);
1307  }//
1308  }// pt bin 80-120
1309 
1310  if( gjet->pt() >= 120 && gjet->pt() < 180 ){
1311  if(isCentral){
1312  mPtRecoOverGen_GenEta_120_180_Cent_0_10->Fill(geneta, response);
1313  if(inBarrel)
1314  mPtRecoOverGen_B_120_180_Cent_0_10->Fill(response);
1315  if(inEndcap)
1316  mPtRecoOverGen_E_120_180_Cent_0_10->Fill(response);
1317  if(inForward)
1318  mPtRecoOverGen_F_120_180_Cent_0_10->Fill(response);
1319  }//
1320  if(ismidCentral){
1321  mPtRecoOverGen_GenEta_120_180_Cent_10_30->Fill(geneta, response);
1322  if(inBarrel)
1323  mPtRecoOverGen_B_120_180_Cent_10_30->Fill(response);
1324  if(inEndcap)
1325  mPtRecoOverGen_E_120_180_Cent_10_30->Fill(response);
1326  if(inForward)
1327  mPtRecoOverGen_F_120_180_Cent_10_30->Fill(response);
1328  }//
1329  if(ismidPeripheral){
1330  mPtRecoOverGen_GenEta_120_180_Cent_30_50->Fill(geneta, response);
1331  if(inBarrel)
1332  mPtRecoOverGen_B_120_180_Cent_30_50->Fill(response);
1333  if(inEndcap)
1334  mPtRecoOverGen_E_120_180_Cent_30_50->Fill(response);
1335  if(inForward)
1336  mPtRecoOverGen_F_120_180_Cent_30_50->Fill(response);
1337  }//
1338  if(isPeripheral){
1339  mPtRecoOverGen_GenEta_120_180_Cent_50_80->Fill(geneta, response);
1340  if(inBarrel)
1341  mPtRecoOverGen_B_120_180_Cent_50_80->Fill(response);
1342  if(inEndcap)
1343  mPtRecoOverGen_E_120_180_Cent_50_80->Fill(response);
1344  if(inForward)
1345  mPtRecoOverGen_F_120_180_Cent_50_80->Fill(response);
1346  }//
1347  }// pt bin 120-180
1348 
1349  if( gjet->pt() >= 180 && gjet->pt() < 300 ){
1350  if(isCentral){
1351  mPtRecoOverGen_GenEta_180_300_Cent_0_10->Fill(geneta, response);
1352  if(inBarrel)
1353  mPtRecoOverGen_B_180_300_Cent_0_10->Fill(response);
1354  if(inEndcap)
1355  mPtRecoOverGen_E_180_300_Cent_0_10->Fill(response);
1356  if(inForward)
1357  mPtRecoOverGen_F_180_300_Cent_0_10->Fill(response);
1358  }//
1359  if(ismidCentral){
1360  mPtRecoOverGen_GenEta_180_300_Cent_10_30->Fill(geneta, response);
1361  if(inBarrel)
1362  mPtRecoOverGen_B_180_300_Cent_10_30->Fill(response);
1363  if(inEndcap)
1364  mPtRecoOverGen_E_180_300_Cent_10_30->Fill(response);
1365  if(inForward)
1366  mPtRecoOverGen_F_180_300_Cent_10_30->Fill(response);
1367  }//
1368  if(ismidPeripheral){
1369  mPtRecoOverGen_GenEta_180_300_Cent_30_50->Fill(geneta, response);
1370  if(inBarrel)
1371  mPtRecoOverGen_B_180_300_Cent_30_50->Fill(response);
1372  if(inEndcap)
1373  mPtRecoOverGen_E_180_300_Cent_30_50->Fill(response);
1374  if(inForward)
1375  mPtRecoOverGen_F_180_300_Cent_30_50->Fill(response);
1376  }//
1377  if(isPeripheral){
1378  mPtRecoOverGen_GenEta_180_300_Cent_50_80->Fill(geneta, response);
1379  if(inBarrel)
1380  mPtRecoOverGen_B_180_300_Cent_50_80->Fill(response);
1381  if(inEndcap)
1382  mPtRecoOverGen_E_180_300_Cent_50_80->Fill(response);
1383  if(inForward)
1384  mPtRecoOverGen_F_180_300_Cent_50_80->Fill(response);
1385  }//
1386  }// pt bin 180-300
1387 
1388  if(gjet->pt() >= 300){
1389  if(isCentral){
1390  mPtRecoOverGen_GenEta_300_Inf_Cent_0_10->Fill(geneta, response);
1391  if(inBarrel)
1392  mPtRecoOverGen_B_300_Inf_Cent_0_10->Fill(response);
1393  if(inEndcap)
1394  mPtRecoOverGen_E_300_Inf_Cent_0_10->Fill(response);
1395  if(inForward)
1396  mPtRecoOverGen_F_300_Inf_Cent_0_10->Fill(response);
1397  }//
1398  if(ismidCentral){
1399  mPtRecoOverGen_GenEta_300_Inf_Cent_10_30->Fill(geneta, response);
1400  if(inBarrel)
1401  mPtRecoOverGen_B_300_Inf_Cent_10_30->Fill(response);
1402  if(inEndcap)
1403  mPtRecoOverGen_E_300_Inf_Cent_10_30->Fill(response);
1404  if(inForward)
1405  mPtRecoOverGen_F_300_Inf_Cent_10_30->Fill(response);
1406  }//
1407  if(ismidPeripheral){
1408  mPtRecoOverGen_GenEta_300_Inf_Cent_30_50->Fill(geneta, response);
1409  if(inBarrel)
1410  mPtRecoOverGen_B_300_Inf_Cent_30_50->Fill(response);
1411  if(inEndcap)
1412  mPtRecoOverGen_E_300_Inf_Cent_30_50->Fill(response);
1413  if(inForward)
1414  mPtRecoOverGen_F_300_Inf_Cent_30_50->Fill(response);
1415  }//
1416  if(isPeripheral){
1417  mPtRecoOverGen_GenEta_300_Inf_Cent_50_80->Fill(geneta, response);
1418  if(inBarrel)
1419  mPtRecoOverGen_B_300_Inf_Cent_50_80->Fill(response);
1420  if(inEndcap)
1421  mPtRecoOverGen_E_300_Inf_Cent_50_80->Fill(response);
1422  if(inForward)
1423  mPtRecoOverGen_F_300_Inf_Cent_50_80->Fill(response);
1424  }//
1425  }// pt bin 300-Inf
1426 
1427  }// delta R < mRthreshold
1428 
1429  }// gen jet collection loop
1430 
1431 
1432  }// not empty gen collection
1433 
1434  }// is the event real
1435 
1436 }
1437 
JetTester_HeavyIons::mSumCaloPt_0p522_0p783
MonitorElement * mSumCaloPt_0p522_0p783
Definition: JetTester_HeavyIons.h:270
JetTester_HeavyIons::mPtRecoOverGen_B_50_80_Cent_0_10
MonitorElement * mPtRecoOverGen_B_50_80_Cent_0_10
Definition: JetTester_HeavyIons.h:303
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
JetTester_HeavyIons::~JetTester_HeavyIons
virtual ~JetTester_HeavyIons()
Definition: JetTester_HeavyIons.cc:639
JetTester_HeavyIons::mPtRecoOverGen_F_20_30_Cent_30_50
MonitorElement * mPtRecoOverGen_F_20_30_Cent_30_50
Definition: JetTester_HeavyIons.h:343
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
JetTester_HeavyIons::mPtRecoOverGen_F_120_180_Cent_30_50
MonitorElement * mPtRecoOverGen_F_120_180_Cent_30_50
Definition: JetTester_HeavyIons.h:355
JetTester_HeavyIons::mPtRecoOverGen_E_30_50_Cent_50_80
MonitorElement * mPtRecoOverGen_E_30_50_Cent_50_80
Definition: JetTester_HeavyIons.h:367
JetTester_HeavyIons::mHF
MonitorElement * mHF
Definition: JetTester_HeavyIons.h:167
JetTester_HeavyIons::mPtRecoOverGen_F_300_Inf_Cent_0_10
MonitorElement * mPtRecoOverGen_F_300_Inf_Cent_0_10
Definition: JetTester_HeavyIons.h:317
JetTester_HeavyIons::mOutputFile
std::string mOutputFile
Definition: JetTester_HeavyIons.h:104
JetTester_HeavyIons::mSumCaloPt_n5p191_n2p650
MonitorElement * mSumCaloPt_n5p191_n2p650
Definition: JetTester_HeavyIons.h:262
JetTester_HeavyIons::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: JetTester_HeavyIons.cc:312
JetTester_HeavyIons::mPtRecoOverGen_F_30_50_Cent_0_10
MonitorElement * mPtRecoOverGen_F_30_50_Cent_0_10
Definition: JetTester_HeavyIons.h:302
JetTester_HeavyIons::mSumCaloPt_2p043_2p650
MonitorElement * mSumCaloPt_2p043_2p650
Definition: JetTester_HeavyIons.h:275
JetTester_HeavyIons::mSumPFPt_n0p783_n0p522
MonitorElement * mSumPFPt_n0p783_n0p522
Definition: JetTester_HeavyIons.h:219
JetTester_HeavyIons::mPFCandpT_Barrel_photon
MonitorElement * mPFCandpT_Barrel_photon
Definition: JetTester_HeavyIons.h:444
JetTester_HeavyIons::mPtRecoOverGen_GenEta_20_30_Cent_30_50
MonitorElement * mPtRecoOverGen_GenEta_20_30_Cent_30_50
Definition: JetTester_HeavyIons.h:415
mps_fire.i
i
Definition: mps_fire.py:269
JetTester_HeavyIons::isCaloJet
bool isCaloJet
Definition: JetTester_HeavyIons.h:470
reco::LeafCandidate::eta
double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:137
JetTester_HeavyIons::mPtRecoOverGen_F_120_180_Cent_50_80
MonitorElement * mPtRecoOverGen_F_120_180_Cent_50_80
Definition: JetTester_HeavyIons.h:377
GenEventInfoProduct::binningValues
const std::vector< double > & binningValues() const
Definition: GenEventInfoProduct.h:44
JetTester_HeavyIons::mInputCollection
edm::InputTag mInputCollection
Definition: JetTester_HeavyIons.h:91
JetTester_HeavyIons::caloTowersToken_
edm::EDGetTokenT< CaloTowerCollection > caloTowersToken_
Definition: JetTester_HeavyIons.h:117
EndcapEta
const Double_t EndcapEta
Definition: JetAnalyzer_HeavyIons.h:78
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
JetTester_HeavyIons::mPtRecoOverGen_F_30_50_Cent_30_50
MonitorElement * mPtRecoOverGen_F_30_50_Cent_30_50
Definition: JetTester_HeavyIons.h:346
JetTester_HeavyIons::mPtRecoOverGen_GenEta_50_80_Cent_30_50
MonitorElement * mPtRecoOverGen_GenEta_50_80_Cent_30_50
Definition: JetTester_HeavyIons.h:417
JetTester_HeavyIons::mNCalopart
MonitorElement * mNCalopart
Definition: JetTester_HeavyIons.h:140
JetTester_HeavyIons::mCaloArea
MonitorElement * mCaloArea
Definition: JetTester_HeavyIons.h:146
JetTester_HeavyIons::mPtRecoOverGen_GenEta_300_Inf_Cent_50_80
MonitorElement * mPtRecoOverGen_GenEta_300_Inf_Cent_50_80
Definition: JetTester_HeavyIons.h:429
CaloTower::p4
math::PtEtaPhiMLorentzVector p4(double vtxZ) const
Definition: CaloTower.cc:129
JetTester_HeavyIons::JetTester_HeavyIons
JetTester_HeavyIons(const edm::ParameterSet &)
Definition: JetTester_HeavyIons.cc:11
JetValidationHeavyIons_cff.UEAlgo
UEAlgo
Definition: JetValidationHeavyIons_cff.py:22
JetTester_HeavyIons::mSumCaloPt_1p479_1p740
MonitorElement * mSumCaloPt_1p479_1p740
Definition: JetTester_HeavyIons.h:273
GenEventInfoProduct::hasBinningValues
bool hasBinningValues() const
Definition: GenEventInfoProduct.h:45
DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
JetTester_HeavyIons::mSumCaloPt_n0p522_0p522
MonitorElement * mSumCaloPt_n0p522_0p522
Definition: JetTester_HeavyIons.h:269
JetTester_HeavyIons::mEnergy
MonitorElement * mEnergy
Definition: JetTester_HeavyIons.h:283
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:508
JetTester_HeavyIons::etaBins_
static const Int_t etaBins_
Definition: JetTester_HeavyIons.h:475
JetTester_HeavyIons::mPtRecoOverGen_F_80_120_Cent_30_50
MonitorElement * mPtRecoOverGen_F_80_120_Cent_30_50
Definition: JetTester_HeavyIons.h:352
JetTester_HeavyIons::mPFPhi
MonitorElement * mPFPhi
Definition: JetTester_HeavyIons.h:136
JetTester_HeavyIons::mPtRecoOverGen_E_50_80_Cent_30_50
MonitorElement * mPtRecoOverGen_E_50_80_Cent_30_50
Definition: JetTester_HeavyIons.h:348
JetTester_HeavyIons::mSumCaloPt
MonitorElement * mSumCaloPt
Definition: JetTester_HeavyIons.h:155
JetTester_HeavyIons::mSumCaloPt_n1p131_n0p783
MonitorElement * mSumCaloPt_n1p131_n0p783
Definition: JetTester_HeavyIons.h:267
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
JetTester_HeavyIons::mPtRecoOverGen_B_120_180_Cent_30_50
MonitorElement * mPtRecoOverGen_B_120_180_Cent_30_50
Definition: JetTester_HeavyIons.h:353
JetTester_HeavyIons::mPtRecoOverGen_GenEta_180_300_Cent_30_50
MonitorElement * mPtRecoOverGen_GenEta_180_300_Cent_30_50
Definition: JetTester_HeavyIons.h:420
JetTester_HeavyIons::pfCandToken_
edm::EDGetTokenT< reco::PFCandidateCollection > pfCandToken_
Definition: JetTester_HeavyIons.h:124
JetTester_HeavyIons::mSumSquaredCaloPt
MonitorElement * mSumSquaredCaloPt
Definition: JetTester_HeavyIons.h:163
JetTester_HeavyIons::mPtRecoOverGen_E_300_Inf_Cent_10_30
MonitorElement * mPtRecoOverGen_E_300_Inf_Cent_10_30
Definition: JetTester_HeavyIons.h:338
JetTester_HeavyIons::mCaloPt
MonitorElement * mCaloPt
Definition: JetTester_HeavyIons.h:141
isolatedTracks_cfi.caloJets
caloJets
Definition: isolatedTracks_cfi.py:33
JetTester_HeavyIons::centralityBinTag_
edm::InputTag centralityBinTag_
Definition: JetTester_HeavyIons.h:100
JetTester_HeavyIons::mPtRecoOverGen_E_120_180_Cent_50_80
MonitorElement * mPtRecoOverGen_E_120_180_Cent_50_80
Definition: JetTester_HeavyIons.h:376
JetTester_HeavyIons::mPFCandpT_vs_eta_EME_inHF
MonitorElement * mPFCandpT_vs_eta_EME_inHF
Definition: JetTester_HeavyIons.h:438
JetTester_HeavyIons::mSumCaloPt_n2p043_n1p740
MonitorElement * mSumCaloPt_n2p043_n1p740
Definition: JetTester_HeavyIons.h:264
JetTester_HeavyIons::mPtRecoOverGen_F_30_50_Cent_10_30
MonitorElement * mPtRecoOverGen_F_30_50_Cent_10_30
Definition: JetTester_HeavyIons.h:324
JetTester_HeavyIons::mPtRecoOverGen_B_300_Inf_Cent_30_50
MonitorElement * mPtRecoOverGen_B_300_Inf_Cent_30_50
Definition: JetTester_HeavyIons.h:359
JetTester_HeavyIons::mNJets
MonitorElement * mNJets
Definition: JetTester_HeavyIons.h:288
JetTester_HeavyIons::mPtRecoOverGen_E_50_80_Cent_10_30
MonitorElement * mPtRecoOverGen_E_50_80_Cent_10_30
Definition: JetTester_HeavyIons.h:326
JetTester_HeavyIons::mPtRecoOverGen_E_30_50_Cent_10_30
MonitorElement * mPtRecoOverGen_E_30_50_Cent_10_30
Definition: JetTester_HeavyIons.h:323
JetTester_HeavyIons::mPtRecoOverGen_GenPt_E_Cent_0_10
MonitorElement * mPtRecoOverGen_GenPt_E_Cent_0_10
Definition: JetTester_HeavyIons.h:387
JetTester_HeavyIons::mPtRecoOverGen_E_20_30_Cent_50_80
MonitorElement * mPtRecoOverGen_E_20_30_Cent_50_80
Definition: JetTester_HeavyIons.h:364
JetTester_HeavyIons::mPtRecoOverGen_GenEta_50_80_Cent_0_10
MonitorElement * mPtRecoOverGen_GenEta_50_80_Cent_0_10
Definition: JetTester_HeavyIons.h:401
JetTester_HeavyIons::mPtRecoOverGen_E_80_120_Cent_50_80
MonitorElement * mPtRecoOverGen_E_80_120_Cent_50_80
Definition: JetTester_HeavyIons.h:373
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:131
JetTester_HeavyIons::backgrounds_value_
edm::EDGetTokenT< std::vector< float > > backgrounds_value_
Definition: JetTester_HeavyIons.h:129
JetTester_HeavyIons::mPtRecoOverGen_B_300_Inf_Cent_0_10
MonitorElement * mPtRecoOverGen_B_300_Inf_Cent_0_10
Definition: JetTester_HeavyIons.h:315
JetTester_HeavyIons::mInputPFCandCollection
edm::InputTag mInputPFCandCollection
Definition: JetTester_HeavyIons.h:93
JetTester_HeavyIons::mSumSquaredPFPt
MonitorElement * mSumSquaredPFPt
Definition: JetTester_HeavyIons.h:159
JetTester_HeavyIons::mSumPFPt_n2p650_n2p043
MonitorElement * mSumPFPt_n2p650_n2p043
Definition: JetTester_HeavyIons.h:214
JetTester_HeavyIons::mSumCaloPt_n2p650_n2p043
MonitorElement * mSumCaloPt_n2p650_n2p043
Definition: JetTester_HeavyIons.h:263
JetTester_HeavyIons::mPtRecoOverGen_GenEta_20_30_Cent_10_30
MonitorElement * mPtRecoOverGen_GenEta_20_30_Cent_10_30
Definition: JetTester_HeavyIons.h:407
JetTester_HeavyIons::mPtRecoOverGen_GenPt_B_Cent_10_30
MonitorElement * mPtRecoOverGen_GenPt_B_Cent_10_30
Definition: JetTester_HeavyIons.h:389
JetValidationHeavyIons_cff.Background
Background
Definition: JetValidationHeavyIons_cff.py:27
JetTester_HeavyIons::mPtRecoOverGen_E_20_30_Cent_0_10
MonitorElement * mPtRecoOverGen_E_20_30_Cent_0_10
Definition: JetTester_HeavyIons.h:298
JetTester_HeavyIons::mGenPt
MonitorElement * mGenPt
Definition: JetTester_HeavyIons.h:294
JetTester_HeavyIons::mPtRecoOverGen_GenPt_F_Cent_50_80
MonitorElement * mPtRecoOverGen_GenPt_F_Cent_50_80
Definition: JetTester_HeavyIons.h:397
JetTester_HeavyIons::mPtRecoOverGen_B_20_30_Cent_0_10
MonitorElement * mPtRecoOverGen_B_20_30_Cent_0_10
Definition: JetTester_HeavyIons.h:297
JetTester_HeavyIons::mSumCaloPt_n1p740_n1p479
MonitorElement * mSumCaloPt_n1p740_n1p479
Definition: JetTester_HeavyIons.h:265
JetTester_HeavyIons::mPhi
MonitorElement * mPhi
Definition: JetTester_HeavyIons.h:280
JetTester_HeavyIons::mPFCandpT_Barrel_ChargedHadron
MonitorElement * mPFCandpT_Barrel_ChargedHadron
Definition: JetTester_HeavyIons.h:441
edm::EventBase::isRealData
bool isRealData() const
Definition: EventBase.h:64
JetTester_HeavyIons::mPtRecoOverGen_B_50_80_Cent_30_50
MonitorElement * mPtRecoOverGen_B_50_80_Cent_30_50
Definition: JetTester_HeavyIons.h:347
JetTester_HeavyIons::mPtRecoOverGen_GenEta_180_300_Cent_0_10
MonitorElement * mPtRecoOverGen_GenEta_180_300_Cent_0_10
Definition: JetTester_HeavyIons.h:404
JetTester_HeavyIons::mPtRecoOverGen_GenEta_30_50_Cent_0_10
MonitorElement * mPtRecoOverGen_GenEta_30_50_Cent_0_10
Definition: JetTester_HeavyIons.h:400
JetTester_HeavyIons::mPFCandpT_vs_eta_muon
MonitorElement * mPFCandpT_vs_eta_muon
Definition: JetTester_HeavyIons.h:434
JetTester_HeavyIons::mPtRecoOverGen_B_180_300_Cent_10_30
MonitorElement * mPtRecoOverGen_B_180_300_Cent_10_30
Definition: JetTester_HeavyIons.h:334
JetTester_HeavyIons::caloCandViewToken_
edm::EDGetTokenT< reco::CandidateView > caloCandViewToken_
Definition: JetTester_HeavyIons.h:126
JetTester_HeavyIons::mSumPFPt_0p783_1p131
MonitorElement * mSumPFPt_0p783_1p131
Definition: JetTester_HeavyIons.h:222
JetTester_HeavyIons::mSumPFPt_1p479_1p740
MonitorElement * mSumPFPt_1p479_1p740
Definition: JetTester_HeavyIons.h:224
JetTester_HeavyIons::mPtRecoOverGen_B_20_30_Cent_30_50
MonitorElement * mPtRecoOverGen_B_20_30_Cent_30_50
Definition: JetTester_HeavyIons.h:341
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
JetTester_HeavyIons::mPtRecoOverGen_B_20_30_Cent_10_30
MonitorElement * mPtRecoOverGen_B_20_30_Cent_10_30
Definition: JetTester_HeavyIons.h:319
JetTester_HeavyIons::mPtRecoOverGen_F_30_50_Cent_50_80
MonitorElement * mPtRecoOverGen_F_30_50_Cent_50_80
Definition: JetTester_HeavyIons.h:368
JetTester_HeavyIons::mPFCandpT_Barrel_HadE_inHF
MonitorElement * mPFCandpT_Barrel_HadE_inHF
Definition: JetTester_HeavyIons.h:446
JetTester_HeavyIons::Background
edm::InputTag Background
Definition: JetTester_HeavyIons.h:107
JetTester_HeavyIons::mRThreshold
double mRThreshold
Definition: JetTester_HeavyIons.h:112
JetTester_HeavyIons::mPtHat
MonitorElement * mPtHat
Definition: JetTester_HeavyIons.h:295
JetTester_HeavyIons::mPtRecoOverGen_GenPt_B_Cent_50_80
MonitorElement * mPtRecoOverGen_GenPt_B_Cent_50_80
Definition: JetTester_HeavyIons.h:395
JetTester_HeavyIons::mPFCandpT_vs_eta_NeutralHadron
MonitorElement * mPFCandpT_vs_eta_NeutralHadron
Definition: JetTester_HeavyIons.h:436
JetTester_HeavyIons::mEta
MonitorElement * mEta
Definition: JetTester_HeavyIons.h:279
JetTester_HeavyIons::mSumPFPt_n1p131_n0p783
MonitorElement * mSumPFPt_n1p131_n0p783
Definition: JetTester_HeavyIons.h:218
estimatePileup_makeJSON.pileup
pileup
Definition: estimatePileup_makeJSON.py:193
JetTester_HeavyIons::mPtRecoOverGen_E_20_30_Cent_30_50
MonitorElement * mPtRecoOverGen_E_20_30_Cent_30_50
Definition: JetTester_HeavyIons.h:342
JetTester_HeavyIons::mPFCandpT_Barrel_Unknown
MonitorElement * mPFCandpT_Barrel_Unknown
Definition: JetTester_HeavyIons.h:440
JetTester_HeavyIons::mPFCandpT_Barrel_EME_inHF
MonitorElement * mPFCandpT_Barrel_EME_inHF
Definition: JetTester_HeavyIons.h:447
JetTester_HeavyIons::mPFCandpT_Forward_Unknown
MonitorElement * mPFCandpT_Forward_Unknown
Definition: JetTester_HeavyIons.h:458
JetTester_HeavyIons::mPtRecoOverGen_F_180_300_Cent_10_30
MonitorElement * mPtRecoOverGen_F_180_300_Cent_10_30
Definition: JetTester_HeavyIons.h:336
slimmedMuons_cfi.pfCandidates
pfCandidates
Definition: slimmedMuons_cfi.py:6
JetTester_HeavyIons::mSumCaloPt_1p131_1p479
MonitorElement * mSumCaloPt_1p131_1p479
Definition: JetTester_HeavyIons.h:272
JetTester_HeavyIons::mNJets_40
MonitorElement * mNJets_40
Definition: JetTester_HeavyIons.h:289
JetTester_HeavyIons::mPtRecoOverGen_F_300_Inf_Cent_50_80
MonitorElement * mPtRecoOverGen_F_300_Inf_Cent_50_80
Definition: JetTester_HeavyIons.h:383
JetTester_HeavyIons::centralityToken
edm::EDGetTokenT< reco::Centrality > centralityToken
Definition: JetTester_HeavyIons.h:97
JetTester_HeavyIons::mPtRecoOverGen_F_20_30_Cent_10_30
MonitorElement * mPtRecoOverGen_F_20_30_Cent_10_30
Definition: JetTester_HeavyIons.h:321
JetTester_HeavyIons::mPtRecoOverGen_GenEta_120_180_Cent_30_50
MonitorElement * mPtRecoOverGen_GenEta_120_180_Cent_30_50
Definition: JetTester_HeavyIons.h:419
JetTester_HeavyIons::basicJetsToken_
edm::EDGetTokenT< reco::BasicJetCollection > basicJetsToken_
Definition: JetTester_HeavyIons.h:120
JetTester_HeavyIons::mPtRecoOverGen_E_80_120_Cent_10_30
MonitorElement * mPtRecoOverGen_E_80_120_Cent_10_30
Definition: JetTester_HeavyIons.h:329
JetTester_HeavyIons::mPtRecoOverGen_B_80_120_Cent_30_50
MonitorElement * mPtRecoOverGen_B_80_120_Cent_30_50
Definition: JetTester_HeavyIons.h:350
JetTester_HeavyIons::mPtRecoOverGen_F_50_80_Cent_10_30
MonitorElement * mPtRecoOverGen_F_50_80_Cent_10_30
Definition: JetTester_HeavyIons.h:327
JetTester_HeavyIons::mPFCandpT_Endcap_NeutralHadron
MonitorElement * mPFCandpT_Endcap_NeutralHadron
Definition: JetTester_HeavyIons.h:454
JetTester_HeavyIons::mPtRecoOverGen_E_80_120_Cent_0_10
MonitorElement * mPtRecoOverGen_E_80_120_Cent_0_10
Definition: JetTester_HeavyIons.h:307
MuonErrorMatrixValues_cff.etaRange
etaRange
Definition: MuonErrorMatrixValues_cff.py:4
JetTester_HeavyIons::mPtRecoOverGen_E_30_50_Cent_0_10
MonitorElement * mPtRecoOverGen_E_30_50_Cent_0_10
Definition: JetTester_HeavyIons.h:301
JetTester_HeavyIons::mSumPFPt_n0p522_0p522
MonitorElement * mSumPFPt_n0p522_0p522
Definition: JetTester_HeavyIons.h:220
JetTester_HeavyIons::mSumPFPt_n5p191_n2p650
MonitorElement * mSumPFPt_n5p191_n2p650
Definition: JetTester_HeavyIons.h:213
JetTester_HeavyIons::mPtRecoOverGen_GenPt_E_Cent_10_30
MonitorElement * mPtRecoOverGen_GenPt_E_Cent_10_30
Definition: JetTester_HeavyIons.h:390
JetTester_HeavyIons::hiVertexToken_
edm::EDGetTokenT< std::vector< reco::Vertex > > hiVertexToken_
Definition: JetTester_HeavyIons.h:130
JetTester_HeavyIons::mPtRecoOverGen_E_30_50_Cent_30_50
MonitorElement * mPtRecoOverGen_E_30_50_Cent_30_50
Definition: JetTester_HeavyIons.h:345
JetTester_HeavyIons::mJetArea
MonitorElement * mJetArea
Definition: JetTester_HeavyIons.h:286
JetTester_HeavyIons::mPFCandpT_Forward_NeutralHadron
MonitorElement * mPFCandpT_Forward_NeutralHadron
Definition: JetTester_HeavyIons.h:463
JetTester_HeavyIons::mPtRecoOverGen_GenEta_30_50_Cent_50_80
MonitorElement * mPtRecoOverGen_GenEta_30_50_Cent_50_80
Definition: JetTester_HeavyIons.h:424
JetTester_HeavyIons::mPFCandpT_vs_eta_photon
MonitorElement * mPFCandpT_vs_eta_photon
Definition: JetTester_HeavyIons.h:435
JetTester_HeavyIons::mPFCandpT_Forward_photon
MonitorElement * mPFCandpT_Forward_photon
Definition: JetTester_HeavyIons.h:462
JetTester_HeavyIons::mPFCandpT_Endcap_photon
MonitorElement * mPFCandpT_Endcap_photon
Definition: JetTester_HeavyIons.h:453
JetTester_HeavyIons::mPtRecoOverGen_B_80_120_Cent_10_30
MonitorElement * mPtRecoOverGen_B_80_120_Cent_10_30
Definition: JetTester_HeavyIons.h:328
JetTester_HeavyIons::mSumPFPt_n2p043_n1p740
MonitorElement * mSumPFPt_n2p043_n1p740
Definition: JetTester_HeavyIons.h:215
JetTester_HeavyIons::mSumCaloPt_2p650_5p191
MonitorElement * mSumCaloPt_2p650_5p191
Definition: JetTester_HeavyIons.h:276
JetTester_HeavyIons::isJPTJet
bool isJPTJet
Definition: JetTester_HeavyIons.h:471
JetTester_HeavyIons::mPtRecoOverGen_F_20_30_Cent_0_10
MonitorElement * mPtRecoOverGen_F_20_30_Cent_0_10
Definition: JetTester_HeavyIons.h:299
JetTester_HeavyIons::mPtRecoOverGen_GenEta_20_30_Cent_0_10
MonitorElement * mPtRecoOverGen_GenEta_20_30_Cent_0_10
Definition: JetTester_HeavyIons.h:399
JetTester_HeavyIons::mSumPFPt
MonitorElement * mSumPFPt
Definition: JetTester_HeavyIons.h:151
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
JetTester_HeavyIons::mPtRecoOverGen_B_30_50_Cent_0_10
MonitorElement * mPtRecoOverGen_B_30_50_Cent_0_10
Definition: JetTester_HeavyIons.h:300
JetTester_HeavyIons::mPtRecoOverGen_GenEta_50_80_Cent_10_30
MonitorElement * mPtRecoOverGen_GenEta_50_80_Cent_10_30
Definition: JetTester_HeavyIons.h:409
JetTester_HeavyIons::mPtRecoOverGen_E_180_300_Cent_10_30
MonitorElement * mPtRecoOverGen_E_180_300_Cent_10_30
Definition: JetTester_HeavyIons.h:335
JetTester_HeavyIons::mPtRecoOverGen_F_50_80_Cent_50_80
MonitorElement * mPtRecoOverGen_F_50_80_Cent_50_80
Definition: JetTester_HeavyIons.h:371
JetTester_HeavyIons::mPtRecoOverGen_F_180_300_Cent_30_50
MonitorElement * mPtRecoOverGen_F_180_300_Cent_30_50
Definition: JetTester_HeavyIons.h:358
JetTester_HeavyIons::mPtRecoOverGen_F_300_Inf_Cent_30_50
MonitorElement * mPtRecoOverGen_F_300_Inf_Cent_30_50
Definition: JetTester_HeavyIons.h:361
JetTester_HeavyIons::mPtRecoOverGen_GenPt_E_Cent_30_50
MonitorElement * mPtRecoOverGen_GenPt_E_Cent_30_50
Definition: JetTester_HeavyIons.h:393
JetTester_HeavyIons::mPtRecoOverGen_B_30_50_Cent_10_30
MonitorElement * mPtRecoOverGen_B_30_50_Cent_10_30
Definition: JetTester_HeavyIons.h:322
reco::Vertex::Point
math::XYZPoint Point
point in the space
Definition: Vertex.h:39
JetTester_HeavyIons::mNPFpart
MonitorElement * mNPFpart
Definition: JetTester_HeavyIons.h:133
JetTester_HeavyIons::mPtRecoOverGen_GenPt_F_Cent_10_30
MonitorElement * mPtRecoOverGen_GenPt_F_Cent_10_30
Definition: JetTester_HeavyIons.h:391
JetTester_HeavyIons::mPFPt
MonitorElement * mPFPt
Definition: JetTester_HeavyIons.h:134
JetTester_HeavyIons::mPtRecoOverGen_F_80_120_Cent_50_80
MonitorElement * mPtRecoOverGen_F_80_120_Cent_50_80
Definition: JetTester_HeavyIons.h:374
JetTester_HeavyIons::evtToken_
edm::EDGetTokenT< GenEventInfoProduct > evtToken_
Definition: JetTester_HeavyIons.h:123
JetTester_HeavyIons::mPtRecoOverGen_F_180_300_Cent_0_10
MonitorElement * mPtRecoOverGen_F_180_300_Cent_0_10
Definition: JetTester_HeavyIons.h:314
JetTester_HeavyIons::mPtRecoOverGen_F_50_80_Cent_0_10
MonitorElement * mPtRecoOverGen_F_50_80_Cent_0_10
Definition: JetTester_HeavyIons.h:305
JetTester_HeavyIons::mPtRecoOverGen_B_120_180_Cent_10_30
MonitorElement * mPtRecoOverGen_B_120_180_Cent_10_30
Definition: JetTester_HeavyIons.h:331
RecoJets_cff.recoJets
recoJets
Definition: RecoJets_cff.py:16
JetTester_HeavyIons::jptJetsToken_
edm::EDGetTokenT< reco::JPTJetCollection > jptJetsToken_
Definition: JetTester_HeavyIons.h:121
JetTester_HeavyIons::mPFCandpT_Endcap_EME_inHF
MonitorElement * mPFCandpT_Endcap_EME_inHF
Definition: JetTester_HeavyIons.h:456
JetTester_HeavyIons::mPFCandpT_Barrel_muon
MonitorElement * mPFCandpT_Barrel_muon
Definition: JetTester_HeavyIons.h:443
JetTester_HeavyIons::mPtRecoOverGen_GenEta_120_180_Cent_0_10
MonitorElement * mPtRecoOverGen_GenEta_120_180_Cent_0_10
Definition: JetTester_HeavyIons.h:403
JetTester_HeavyIons::mSumPFPt_0p522_0p783
MonitorElement * mSumPFPt_0p522_0p783
Definition: JetTester_HeavyIons.h:221
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
JetTester_HeavyIons::mPtRecoOverGen_GenEta_120_180_Cent_50_80
MonitorElement * mPtRecoOverGen_GenEta_120_180_Cent_50_80
Definition: JetTester_HeavyIons.h:427
JetTester_HeavyIons::mPtRecoOverGen_E_300_Inf_Cent_30_50
MonitorElement * mPtRecoOverGen_E_300_Inf_Cent_30_50
Definition: JetTester_HeavyIons.h:360
JetTester_HeavyIons::nedge_pseudorapidity
static const size_t nedge_pseudorapidity
Definition: JetTester_HeavyIons.h:477
JetTester_HeavyIons::mPtRecoOverGen_B_300_Inf_Cent_10_30
MonitorElement * mPtRecoOverGen_B_300_Inf_Cent_10_30
Definition: JetTester_HeavyIons.h:337
JetTester_HeavyIons.h
JetTester_HeavyIons::mPtRecoOverGen_F_120_180_Cent_10_30
MonitorElement * mPtRecoOverGen_F_120_180_Cent_10_30
Definition: JetTester_HeavyIons.h:333
JetTester_HeavyIons::mPtRecoOverGen_B_50_80_Cent_50_80
MonitorElement * mPtRecoOverGen_B_50_80_Cent_50_80
Definition: JetTester_HeavyIons.h:369
JetTester_HeavyIons::mPtRecoOverGen_E_20_30_Cent_10_30
MonitorElement * mPtRecoOverGen_E_20_30_Cent_10_30
Definition: JetTester_HeavyIons.h:320
JetTester_HeavyIons::mSumPFPt_1p131_1p479
MonitorElement * mSumPFPt_1p131_1p479
Definition: JetTester_HeavyIons.h:223
JetTester_HeavyIons::mPtRecoOverGen_F_80_120_Cent_10_30
MonitorElement * mPtRecoOverGen_F_80_120_Cent_10_30
Definition: JetTester_HeavyIons.h:330
JetTester_HeavyIons::centralityTag_
edm::InputTag centralityTag_
Definition: JetTester_HeavyIons.h:96
JetTester_HeavyIons::mPFCandpT_Forward_ChargedHadron
MonitorElement * mPFCandpT_Forward_ChargedHadron
Definition: JetTester_HeavyIons.h:459
gen::k
int k[5][pyjets_maxn]
Definition: Cascade2Hadronizer.cc:79
ForwardEta
const Double_t ForwardEta
Definition: JetAnalyzer_HeavyIons.h:79
jets_cff.nConstituents
nConstituents
Definition: jets_cff.py:201
JetTester_HeavyIons::mPFCandpT_Endcap_HadE_inHF
MonitorElement * mPFCandpT_Endcap_HadE_inHF
Definition: JetTester_HeavyIons.h:455
deltaR
double deltaR(double eta1, double eta2, double phi1, double phi2)
Definition: TreeUtility.cc:17
JetTester_HeavyIons::mPFCandpT_Forward_electron
MonitorElement * mPFCandpT_Forward_electron
Definition: JetTester_HeavyIons.h:460
JetTester_HeavyIons::mPtRecoOverGen_F_120_180_Cent_0_10
MonitorElement * mPtRecoOverGen_F_120_180_Cent_0_10
Definition: JetTester_HeavyIons.h:311
JetTester_HeavyIons::mConstituents
MonitorElement * mConstituents
Definition: JetTester_HeavyIons.h:285
JetTester_HeavyIons::mPFCandpT_vs_eta_ChargedHadron
MonitorElement * mPFCandpT_vs_eta_ChargedHadron
Definition: JetTester_HeavyIons.h:432
JetTester_HeavyIons::mSumPFPt_1p740_2p043
MonitorElement * mSumPFPt_1p740_2p043
Definition: JetTester_HeavyIons.h:225
JetTester_HeavyIons::JetType
std::string JetType
Definition: JetTester_HeavyIons.h:105
JetTester_HeavyIons::pfCandViewToken_
edm::EDGetTokenT< reco::CandidateView > pfCandViewToken_
Definition: JetTester_HeavyIons.h:125
reco::PFCandidateCollection
std::vector< reco::PFCandidate > PFCandidateCollection
collection of PFCandidates
Definition: PFCandidateFwd.h:12
JetTester_HeavyIons::mPtRecoOverGen_GenEta_80_120_Cent_50_80
MonitorElement * mPtRecoOverGen_GenEta_80_120_Cent_50_80
Definition: JetTester_HeavyIons.h:426
JetTester_HeavyIons::mPtRecoOverGen_B_30_50_Cent_30_50
MonitorElement * mPtRecoOverGen_B_30_50_Cent_30_50
Definition: JetTester_HeavyIons.h:344
JetTester_HeavyIons::mPtRecoOverGen_B_80_120_Cent_50_80
MonitorElement * mPtRecoOverGen_B_80_120_Cent_50_80
Definition: JetTester_HeavyIons.h:372
JetTester_HeavyIons::genJetsToken_
edm::EDGetTokenT< reco::GenJetCollection > genJetsToken_
Definition: JetTester_HeavyIons.h:122
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
JetTester_HeavyIons::isPFJet
bool isPFJet
Definition: JetTester_HeavyIons.h:472
reco::Centrality::EtHFtowerSum
double EtHFtowerSum() const
Definition: Centrality.h:23
JetTester_HeavyIons::pvToken_
edm::EDGetTokenT< std::vector< reco::Vertex > > pvToken_
Definition: JetTester_HeavyIons.h:116
JetTester_HeavyIons::mPFCandpT_Forward_muon
MonitorElement * mPFCandpT_Forward_muon
Definition: JetTester_HeavyIons.h:461
edm::Handle::product
T const * product() const
Definition: Handle.h:81
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
JetTester_HeavyIons::mPtRecoOverGen_B_80_120_Cent_0_10
MonitorElement * mPtRecoOverGen_B_80_120_Cent_0_10
Definition: JetTester_HeavyIons.h:306
JetTester_HeavyIons::mPtRecoOverGen_GenEta_300_Inf_Cent_30_50
MonitorElement * mPtRecoOverGen_GenEta_300_Inf_Cent_30_50
Definition: JetTester_HeavyIons.h:421
JetTester_HeavyIons::mSumPFPt_n1p479_n1p131
MonitorElement * mSumPFPt_n1p479_n1p131
Definition: JetTester_HeavyIons.h:217
JetTester_HeavyIons::mPtRecoOverGen_B_120_180_Cent_0_10
MonitorElement * mPtRecoOverGen_B_120_180_Cent_0_10
Definition: JetTester_HeavyIons.h:309
JetTester_HeavyIons::mSumPFPt_2p043_2p650
MonitorElement * mSumPFPt_2p043_2p650
Definition: JetTester_HeavyIons.h:226
JetTester_HeavyIons::mPtRecoOverGen_GenPt_F_Cent_30_50
MonitorElement * mPtRecoOverGen_GenPt_F_Cent_30_50
Definition: JetTester_HeavyIons.h:394
JetTester_HeavyIons::UEAlgo
std::string UEAlgo
Definition: JetTester_HeavyIons.h:106
JetTester_HeavyIons::mPFCandpT_vs_eta_HadE_inHF
MonitorElement * mPFCandpT_vs_eta_HadE_inHF
Definition: JetTester_HeavyIons.h:437
JetTester_HeavyIons::mCaloEta
MonitorElement * mCaloEta
Definition: JetTester_HeavyIons.h:142
JetTester_HeavyIons::mPtRecoOverGen_E_120_180_Cent_10_30
MonitorElement * mPtRecoOverGen_E_120_180_Cent_10_30
Definition: JetTester_HeavyIons.h:332
JetTester_HeavyIons::mPtRecoOverGen_GenEta_30_50_Cent_30_50
MonitorElement * mPtRecoOverGen_GenEta_30_50_Cent_30_50
Definition: JetTester_HeavyIons.h:416
JetTester_HeavyIons::mPFCandpT_Endcap_Unknown
MonitorElement * mPFCandpT_Endcap_Unknown
Definition: JetTester_HeavyIons.h:449
JetTester_HeavyIons::mSumCaloPt_n0p783_n0p522
MonitorElement * mSumCaloPt_n0p783_n0p522
Definition: JetTester_HeavyIons.h:268
JetTester_HeavyIons::mPFEta
MonitorElement * mPFEta
Definition: JetTester_HeavyIons.h:135
JetTester_HeavyIons::mGenPhi
MonitorElement * mGenPhi
Definition: JetTester_HeavyIons.h:293
JetTester_HeavyIons::mPtRecoOverGen_B_120_180_Cent_50_80
MonitorElement * mPtRecoOverGen_B_120_180_Cent_50_80
Definition: JetTester_HeavyIons.h:375
JetTester_HeavyIons::mPtRecoOverGen_F_80_120_Cent_0_10
MonitorElement * mPtRecoOverGen_F_80_120_Cent_0_10
Definition: JetTester_HeavyIons.h:308
JetTester_HeavyIons::mPFCandpT_vs_eta_electron
MonitorElement * mPFCandpT_vs_eta_electron
Definition: JetTester_HeavyIons.h:433
JetTester_HeavyIons::mPFCandpT_vs_eta_Unknown
MonitorElement * mPFCandpT_vs_eta_Unknown
Definition: JetTester_HeavyIons.h:431
JetTester_HeavyIons::mCaloPhi
MonitorElement * mCaloPhi
Definition: JetTester_HeavyIons.h:143
JetTester_HeavyIons::mPFCandpT_Endcap_ChargedHadron
MonitorElement * mPFCandpT_Endcap_ChargedHadron
Definition: JetTester_HeavyIons.h:450
JetTester_HeavyIons::mP
MonitorElement * mP
Definition: JetTester_HeavyIons.h:282
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
JetTester_HeavyIons::mPtRecoOverGen_F_50_80_Cent_30_50
MonitorElement * mPtRecoOverGen_F_50_80_Cent_30_50
Definition: JetTester_HeavyIons.h:349
JetTester_HeavyIons::mPFCandpT_Forward_EME_inHF
MonitorElement * mPFCandpT_Forward_EME_inHF
Definition: JetTester_HeavyIons.h:465
JetTester_HeavyIons::mSumCaloPt_0p783_1p131
MonitorElement * mSumCaloPt_0p783_1p131
Definition: JetTester_HeavyIons.h:271
JetTester_HeavyIons::mPtRecoOverGen_E_80_120_Cent_30_50
MonitorElement * mPtRecoOverGen_E_80_120_Cent_30_50
Definition: JetTester_HeavyIons.h:351
reco::PFCandidate
Particle reconstructed by the particle flow algorithm.
Definition: PFCandidate.h:39
JetTester_HeavyIons::mSumpt
MonitorElement * mSumpt
Definition: JetTester_HeavyIons.h:147
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:43
JetTester_HeavyIons::mPtRecoOverGen_GenEta_80_120_Cent_0_10
MonitorElement * mPtRecoOverGen_GenEta_80_120_Cent_0_10
Definition: JetTester_HeavyIons.h:402
JetTester_HeavyIons::mPtRecoOverGen_E_180_300_Cent_30_50
MonitorElement * mPtRecoOverGen_E_180_300_Cent_30_50
Definition: JetTester_HeavyIons.h:357
edm
HLT enums.
Definition: AlignableModifier.h:17
JetTester_HeavyIons::mPtRecoOverGen_E_50_80_Cent_50_80
MonitorElement * mPtRecoOverGen_E_50_80_Cent_50_80
Definition: JetTester_HeavyIons.h:370
JetTester_HeavyIons::centralityBinToken
edm::EDGetTokenT< int > centralityBinToken
Definition: JetTester_HeavyIons.h:101
JetTester_HeavyIons::mSumPFPt_HF
MonitorElement * mSumPFPt_HF
Definition: JetTester_HeavyIons.h:172
JetTester_HeavyIons::mPtRecoOverGen_GenPt_F_Cent_0_10
MonitorElement * mPtRecoOverGen_GenPt_F_Cent_0_10
Definition: JetTester_HeavyIons.h:388
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:236
JetTester_HeavyIons::mPtRecoOverGen_GenPt_B_Cent_0_10
MonitorElement * mPtRecoOverGen_GenPt_B_Cent_0_10
Definition: JetTester_HeavyIons.h:386
position
static int position[264][3]
Definition: ReadPGInfo.cc:509
JetTester_HeavyIons::mPtRecoOverGen_B_180_300_Cent_30_50
MonitorElement * mPtRecoOverGen_B_180_300_Cent_30_50
Definition: JetTester_HeavyIons.h:356
JetTester_HeavyIons::mPtRecoOverGen_F_180_300_Cent_50_80
MonitorElement * mPtRecoOverGen_F_180_300_Cent_50_80
Definition: JetTester_HeavyIons.h:380
JetTester_HeavyIons::mPtRecoOverGen_GenEta_300_Inf_Cent_10_30
MonitorElement * mPtRecoOverGen_GenEta_300_Inf_Cent_10_30
Definition: JetTester_HeavyIons.h:413
JetTester_HeavyIons::mSumPFPt_2p650_5p191
MonitorElement * mSumPFPt_2p650_5p191
Definition: JetTester_HeavyIons.h:227
JetTester_HeavyIons::mGenEta
MonitorElement * mGenEta
Definition: JetTester_HeavyIons.h:292
JetTester_HeavyIons::mPtRecoOverGen_E_300_Inf_Cent_50_80
MonitorElement * mPtRecoOverGen_E_300_Inf_Cent_50_80
Definition: JetTester_HeavyIons.h:382
JetTester_HeavyIons::mPtRecoOverGen_GenPt_B_Cent_30_50
MonitorElement * mPtRecoOverGen_GenPt_B_Cent_30_50
Definition: JetTester_HeavyIons.h:392
JetTester_HeavyIons::mPtRecoOverGen_GenEta_80_120_Cent_10_30
MonitorElement * mPtRecoOverGen_GenEta_80_120_Cent_10_30
Definition: JetTester_HeavyIons.h:410
JetTester_HeavyIons::mPtRecoOverGen_E_120_180_Cent_30_50
MonitorElement * mPtRecoOverGen_E_120_180_Cent_30_50
Definition: JetTester_HeavyIons.h:354
BarrelEta
const Double_t BarrelEta
Definition: JetAnalyzer_HeavyIons.h:77
JetTester_HeavyIons::mPtRecoOverGen_B_180_300_Cent_0_10
MonitorElement * mPtRecoOverGen_B_180_300_Cent_0_10
Definition: JetTester_HeavyIons.h:312
JetTester_HeavyIons::mPtRecoOverGen_B_180_300_Cent_50_80
MonitorElement * mPtRecoOverGen_B_180_300_Cent_50_80
Definition: JetTester_HeavyIons.h:378
JetTester_HeavyIons::mPFCandpT_Barrel_NeutralHadron
MonitorElement * mPFCandpT_Barrel_NeutralHadron
Definition: JetTester_HeavyIons.h:445
JetTester_HeavyIons::mSumCaloPt_HF
MonitorElement * mSumCaloPt_HF
Definition: JetTester_HeavyIons.h:178
JetTester_HeavyIons::mPFCandpT_Forward_HadE_inHF
MonitorElement * mPFCandpT_Forward_HadE_inHF
Definition: JetTester_HeavyIons.h:464
JetTester_HeavyIons::mNvtx
MonitorElement * mNvtx
Definition: JetTester_HeavyIons.h:166
badGlobalMuonTaggersAOD_cff.vtx
vtx
Definition: badGlobalMuonTaggersAOD_cff.py:5
JetTester_HeavyIons::mPtRecoOverGen_GenEta_180_300_Cent_50_80
MonitorElement * mPtRecoOverGen_GenEta_180_300_Cent_50_80
Definition: JetTester_HeavyIons.h:428
JetTester_HeavyIons::mMass
MonitorElement * mMass
Definition: JetTester_HeavyIons.h:284
JetTester_HeavyIons::mMatchGenPtThreshold
double mMatchGenPtThreshold
Definition: JetTester_HeavyIons.h:109
JetTester_HeavyIons::mPtRecoOverGen_GenEta_300_Inf_Cent_0_10
MonitorElement * mPtRecoOverGen_GenEta_300_Inf_Cent_0_10
Definition: JetTester_HeavyIons.h:405
reco::PFCandidate::particleId
virtual ParticleType particleId() const
Definition: PFCandidate.h:373
JetTester_HeavyIons::mInputGenCollection
edm::InputTag mInputGenCollection
Definition: JetTester_HeavyIons.h:92
JetTester_HeavyIons::mPtRecoOverGen_GenEta_50_80_Cent_50_80
MonitorElement * mPtRecoOverGen_GenEta_50_80_Cent_50_80
Definition: JetTester_HeavyIons.h:425
JetTester_HeavyIons::mSumPFPt_n1p740_n1p479
MonitorElement * mSumPFPt_n1p740_n1p479
Definition: JetTester_HeavyIons.h:216
JetTester_HeavyIons::mPFArea
MonitorElement * mPFArea
Definition: JetTester_HeavyIons.h:139
JetTester_HeavyIons::mPtRecoOverGen_GenEta_120_180_Cent_10_30
MonitorElement * mPtRecoOverGen_GenEta_120_180_Cent_10_30
Definition: JetTester_HeavyIons.h:411
JetTester_HeavyIons::mPtRecoOverGen_GenEta_20_30_Cent_50_80
MonitorElement * mPtRecoOverGen_GenEta_20_30_Cent_50_80
Definition: JetTester_HeavyIons.h:423
JetTester_HeavyIons::mPFCandpT_Endcap_electron
MonitorElement * mPFCandpT_Endcap_electron
Definition: JetTester_HeavyIons.h:451
reco::LeafCandidate::phi
double phi() const final
momentum azimuthal angle
Definition: LeafCandidate.h:133
JetTester_HeavyIons::mPtRecoOverGen_B_20_30_Cent_50_80
MonitorElement * mPtRecoOverGen_B_20_30_Cent_50_80
Definition: JetTester_HeavyIons.h:363
pfJetBenchmark_cfi.pfJets
pfJets
Definition: pfJetBenchmark_cfi.py:4
JetTester_HeavyIons::mPt
MonitorElement * mPt
Definition: JetTester_HeavyIons.h:281
JetTester_HeavyIons::mPFCandpT_Barrel_electron
MonitorElement * mPFCandpT_Barrel_electron
Definition: JetTester_HeavyIons.h:442
JetTester_HeavyIons::mPtRecoOverGen_GenEta_180_300_Cent_10_30
MonitorElement * mPtRecoOverGen_GenEta_180_300_Cent_10_30
Definition: JetTester_HeavyIons.h:412
JetTester_HeavyIons::backgrounds_
edm::EDGetTokenT< edm::ValueMap< reco::VoronoiBackground > > backgrounds_
Definition: JetTester_HeavyIons.h:128
JetTester_HeavyIons::mPtRecoOverGen_E_180_300_Cent_50_80
MonitorElement * mPtRecoOverGen_E_180_300_Cent_50_80
Definition: JetTester_HeavyIons.h:379
JetTester_HeavyIons::mPtRecoOverGen_B_50_80_Cent_10_30
MonitorElement * mPtRecoOverGen_B_50_80_Cent_10_30
Definition: JetTester_HeavyIons.h:325
JetTester_HeavyIons::mPtRecoOverGen_GenEta_30_50_Cent_10_30
MonitorElement * mPtRecoOverGen_GenEta_30_50_Cent_10_30
Definition: JetTester_HeavyIons.h:408
JetTester_HeavyIons::caloJetsToken_
edm::EDGetTokenT< reco::CaloJetCollection > caloJetsToken_
Definition: JetTester_HeavyIons.h:118
JetTester_HeavyIons::mPFCandpT_Endcap_muon
MonitorElement * mPFCandpT_Endcap_muon
Definition: JetTester_HeavyIons.h:452
JetTester_HeavyIons::mPtRecoOverGen_F_300_Inf_Cent_10_30
MonitorElement * mPtRecoOverGen_F_300_Inf_Cent_10_30
Definition: JetTester_HeavyIons.h:339
JetTester_HeavyIons::mPtRecoOverGen_E_180_300_Cent_0_10
MonitorElement * mPtRecoOverGen_E_180_300_Cent_0_10
Definition: JetTester_HeavyIons.h:313
JetTester_HeavyIons::mSumCaloPt_1p740_2p043
MonitorElement * mSumCaloPt_1p740_2p043
Definition: JetTester_HeavyIons.h:274
JetTester_HeavyIons::mPtRecoOverGen_GenEta_80_120_Cent_30_50
MonitorElement * mPtRecoOverGen_GenEta_80_120_Cent_30_50
Definition: JetTester_HeavyIons.h:418
edm::Run
Definition: Run.h:43
JetTester_HeavyIons::mRecoJetPtThreshold
double mRecoJetPtThreshold
Definition: JetTester_HeavyIons.h:108
JetTester_HeavyIons::mPtRecoOverGen_E_50_80_Cent_0_10
MonitorElement * mPtRecoOverGen_E_50_80_Cent_0_10
Definition: JetTester_HeavyIons.h:304
JetTester_HeavyIons::mPtRecoOverGen_GenPt_E_Cent_50_80
MonitorElement * mPtRecoOverGen_GenPt_E_Cent_50_80
Definition: JetTester_HeavyIons.h:396
JetTester_HeavyIons::mPtRecoOverGen_B_30_50_Cent_50_80
MonitorElement * mPtRecoOverGen_B_30_50_Cent_50_80
Definition: JetTester_HeavyIons.h:366
JetTester_HeavyIons::mSumCaloPt_n1p479_n1p131
MonitorElement * mSumCaloPt_n1p479_n1p131
Definition: JetTester_HeavyIons.h:266
JetTester_HeavyIons::mPtRecoOverGen_B_300_Inf_Cent_50_80
MonitorElement * mPtRecoOverGen_B_300_Inf_Cent_50_80
Definition: JetTester_HeavyIons.h:381
JetTester_HeavyIons::mPtRecoOverGen_F_20_30_Cent_50_80
MonitorElement * mPtRecoOverGen_F_20_30_Cent_50_80
Definition: JetTester_HeavyIons.h:365
JetTester_HeavyIons::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: JetTester_HeavyIons.cc:665
ttbarCategorization_cff.genJets
genJets
Definition: ttbarCategorization_cff.py:30
JetTester_HeavyIons::mPtRecoOverGen_E_300_Inf_Cent_0_10
MonitorElement * mPtRecoOverGen_E_300_Inf_Cent_0_10
Definition: JetTester_HeavyIons.h:316
JetTester_HeavyIons::mPtRecoOverGen_E_120_180_Cent_0_10
MonitorElement * mPtRecoOverGen_E_120_180_Cent_0_10
Definition: JetTester_HeavyIons.h:310
JetTester_HeavyIons::mjetpileup
MonitorElement * mjetpileup
Definition: JetTester_HeavyIons.h:287
Generated for CMSSW Reference Manual by doxygen 1.8.11