CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
JetAnalyzer_HeavyIons.cc
Go to the documentation of this file.
1 //
2 // Jet Analyzer class for heavy ion jets. for DQM jet analysis monitoring
3 // For CMSSW_7_4_X, especially reading background subtracted jets
4 // author: Raghav Kunnawalkam Elayavalli, Mohammed Zakaria (co Author)
5 // Jan 12th 2015
6 // Rutgers University, email: raghav.k.e at CERN dot CH
7 // UIC, email: mzakaria @ CERN dot CH
8 
9 
10 #include "DQMOffline/JetMET/interface/JetAnalyzer_HeavyIons.h"
11 
12 using namespace edm;
13 using namespace reco;
14 using namespace std;
15 
16 // declare the constructors:
17 
18 JetAnalyzer_HeavyIons::JetAnalyzer_HeavyIons(const edm::ParameterSet& iConfig) :
19  mInputCollection (iConfig.getParameter<edm::InputTag> ("src")),
20  mInputPFCandCollection (iConfig.getParameter<edm::InputTag> ("PFcands")),
21  mOutputFile (iConfig.getUntrackedParameter<std::string>("OutputFile","")),
22  JetType (iConfig.getUntrackedParameter<std::string>("JetType")),
23  UEAlgo (iConfig.getUntrackedParameter<std::string>("UEAlgo")),
24  Background (iConfig.getParameter<edm::InputTag> ("Background")),
25  mRecoJetPtThreshold (iConfig.getParameter<double> ("recoJetPtThreshold")),
26  mReverseEnergyFractionThreshold(iConfig.getParameter<double> ("reverseEnergyFractionThreshold")),
27  mRThreshold (iConfig.getParameter<double> ("RThreshold")),
28  JetCorrectionService (iConfig.getParameter<std::string> ("JetCorrections"))
29 {
30  std::string inputCollectionLabel(mInputCollection.label());
31 
32  isCaloJet = (std::string("calo")==JetType);
33  isJPTJet = (std::string("jpt") ==JetType);
34  isPFJet = (std::string("pf") ==JetType);
35 
36  //consumes
37  pvToken_ = consumes<std::vector<reco::Vertex> >(edm::InputTag("offlinePrimaryVertices"));
38  caloTowersToken_ = consumes<CaloTowerCollection>(edm::InputTag("towerMaker"));
39  if (isCaloJet) caloJetsToken_ = consumes<reco::CaloJetCollection>(mInputCollection);
40  if (isJPTJet) jptJetsToken_ = consumes<reco::JPTJetCollection>(mInputCollection);
41  if (isPFJet) {
42  if(std::string("Pu")==UEAlgo) basicJetsToken_ = consumes<reco::BasicJetCollection>(mInputCollection);
43  if(std::string("Vs")==UEAlgo) pfJetsToken_ = consumes<reco::PFJetCollection>(mInputCollection);
44  }
45 
46  pfCandToken_ = consumes<reco::PFCandidateCollection>(mInputPFCandCollection);
47  pfCandViewToken_ = consumes<reco::CandidateView>(mInputPFCandCollection);
48  caloCandViewToken_ = consumes<reco::CandidateView>(edm::InputTag("towerMaker"));
49  backgrounds_ = consumes<edm::ValueMap<reco::VoronoiBackground> >(Background);
50  backgrounds_value_ = consumes<std::vector<float> >(Background);
51 
52  centralityTag_ = iConfig.getParameter<InputTag>("centralitycollection");
53  centralityToken = consumes<reco::Centrality>(centralityTag_);
54 
55  centralityBinTag_ = (iConfig.getParameter<edm::InputTag> ("centralitybincollection"));
56  centralityBinToken = consumes<int>(centralityBinTag_);
57 
58  hiVertexToken_ = consumes<std::vector<reco::Vertex> >(edm::InputTag("hiSelectedVertex"));
59 
60  // need to initialize the PF cand histograms : which are also event variables
61  if(isPFJet){
62  mNPFpart = 0;
63  mPFPt = 0;
64  mPFEta = 0;
65  mPFPhi = 0;
66  mPFVsPt = 0;
67  mPFVsPtInitial = 0;
68  mPFVsPtInitial_HF = 0;
69  mPFVsPtInitial_Barrel_Centrality_100_95 = 0;
70  mPFVsPtInitial_Barrel_Centrality_95_90 = 0;
71  mPFVsPtInitial_Barrel_Centrality_90_85 = 0;
72  mPFVsPtInitial_Barrel_Centrality_85_80 = 0;
73  mPFVsPtInitial_Barrel_Centrality_80_75 = 0;
74  mPFVsPtInitial_Barrel_Centrality_75_70 = 0;
75  mPFVsPtInitial_Barrel_Centrality_70_65 = 0;
76  mPFVsPtInitial_Barrel_Centrality_65_60 = 0;
77  mPFVsPtInitial_Barrel_Centrality_60_55 = 0;
78  mPFVsPtInitial_Barrel_Centrality_55_50 = 0;
79  mPFVsPtInitial_Barrel_Centrality_50_45 = 0;
80  mPFVsPtInitial_Barrel_Centrality_45_40 = 0;
81  mPFVsPtInitial_Barrel_Centrality_40_35 = 0;
82  mPFVsPtInitial_Barrel_Centrality_35_30 = 0;
83  mPFVsPtInitial_Barrel_Centrality_30_25 = 0;
84  mPFVsPtInitial_Barrel_Centrality_25_20 = 0;
85  mPFVsPtInitial_Barrel_Centrality_20_15 = 0;
86  mPFVsPtInitial_Barrel_Centrality_15_10 = 0;
87  mPFVsPtInitial_Barrel_Centrality_10_5 = 0;
88  mPFVsPtInitial_Barrel_Centrality_5_0 = 0;
89 
90  mPFVsPtInitial_EndCap_Centrality_100_95 = 0;
91  mPFVsPtInitial_EndCap_Centrality_95_90 = 0;
92  mPFVsPtInitial_EndCap_Centrality_90_85 = 0;
93  mPFVsPtInitial_EndCap_Centrality_85_80 = 0;
94  mPFVsPtInitial_EndCap_Centrality_80_75 = 0;
95  mPFVsPtInitial_EndCap_Centrality_75_70 = 0;
96  mPFVsPtInitial_EndCap_Centrality_70_65 = 0;
97  mPFVsPtInitial_EndCap_Centrality_65_60 = 0;
98  mPFVsPtInitial_EndCap_Centrality_60_55 = 0;
99  mPFVsPtInitial_EndCap_Centrality_55_50 = 0;
100  mPFVsPtInitial_EndCap_Centrality_50_45 = 0;
101  mPFVsPtInitial_EndCap_Centrality_45_40 = 0;
102  mPFVsPtInitial_EndCap_Centrality_40_35 = 0;
103  mPFVsPtInitial_EndCap_Centrality_35_30 = 0;
104  mPFVsPtInitial_EndCap_Centrality_30_25 = 0;
105  mPFVsPtInitial_EndCap_Centrality_25_20 = 0;
106  mPFVsPtInitial_EndCap_Centrality_20_15 = 0;
107  mPFVsPtInitial_EndCap_Centrality_15_10 = 0;
108  mPFVsPtInitial_EndCap_Centrality_10_5 = 0;
109  mPFVsPtInitial_EndCap_Centrality_5_0 = 0;
110 
111  mPFVsPtInitial_HF_Centrality_100_95 = 0;
112  mPFVsPtInitial_HF_Centrality_95_90 = 0;
113  mPFVsPtInitial_HF_Centrality_90_85 = 0;
114  mPFVsPtInitial_HF_Centrality_85_80 = 0;
115  mPFVsPtInitial_HF_Centrality_80_75 = 0;
116  mPFVsPtInitial_HF_Centrality_75_70 = 0;
117  mPFVsPtInitial_HF_Centrality_70_65 = 0;
118  mPFVsPtInitial_HF_Centrality_65_60 = 0;
119  mPFVsPtInitial_HF_Centrality_60_55 = 0;
120  mPFVsPtInitial_HF_Centrality_55_50 = 0;
121  mPFVsPtInitial_HF_Centrality_50_45 = 0;
122  mPFVsPtInitial_HF_Centrality_45_40 = 0;
123  mPFVsPtInitial_HF_Centrality_40_35 = 0;
124  mPFVsPtInitial_HF_Centrality_35_30 = 0;
125  mPFVsPtInitial_HF_Centrality_30_25 = 0;
126  mPFVsPtInitial_HF_Centrality_25_20 = 0;
127  mPFVsPtInitial_HF_Centrality_20_15 = 0;
128  mPFVsPtInitial_HF_Centrality_15_10 = 0;
129  mPFVsPtInitial_HF_Centrality_10_5 = 0;
130  mPFVsPtInitial_HF_Centrality_5_0 = 0;
131 
132  mPFArea = 0;
133  mPFDeltaR = 0;
134  mPFDeltaR_Scaled_R = 0;
135 
136  mPFDeltaR_pTCorrected = 0;
137  mPFDeltaR_pTCorrected_PFpT_20To30 = 0;
138  mPFDeltaR_pTCorrected_PFpT_30To50 = 0;
139  mPFDeltaR_pTCorrected_PFpT_50To80 = 0;
140  mPFDeltaR_pTCorrected_PFpT_80To120 = 0;
141  mPFDeltaR_pTCorrected_PFpT_120To180 = 0;
142  mPFDeltaR_pTCorrected_PFpT_180To300 = 0;
143  mPFDeltaR_pTCorrected_PFpT_300ToInf = 0;
144 
145  mPFDeltaR_pTCorrected_PFVsInitialpT_20To30 = 0;
146  mPFDeltaR_pTCorrected_PFVsInitialpT_30To50 = 0;
147  mPFDeltaR_pTCorrected_PFVsInitialpT_50To80 = 0;
148  mPFDeltaR_pTCorrected_PFVsInitialpT_80To120 = 0;
149  mPFDeltaR_pTCorrected_PFVsInitialpT_120To180 = 0;
150  mPFDeltaR_pTCorrected_PFVsInitialpT_180To300 = 0;
151  mPFDeltaR_pTCorrected_PFVsInitialpT_300ToInf = 0;
152 
153 
154  mPFVsPtInitialDeltaR_pTCorrected = 0;
155  mPFVsPtDeltaR_pTCorrected = 0;
156 
157  mSumPFVsPt = 0;
158  mSumPFVsPtInitial = 0;
159  mSumPFPt = 0;
160  mSumPFVsPtInitial_eta = 0;
161  mSumPFVsPt_eta = 0;
162  mSumPFPt_eta = 0;
163  mSumSquaredPFVsPt = 0;
164  mSumSquaredPFVsPtInitial = 0;
165  mSumSquaredPFPt = 0;
166  mSumSquaredPFVsPtInitial_eta = 0;
167  mSumSquaredPFVsPt_eta = 0;
168  mSumSquaredPFPt_eta = 0;
169  mSumPFVsPtInitial_HF = 0;
170  mSumPFVsPt_HF = 0;
171  mSumPFPt_HF = 0;
172 
173  mSumPFVsPtInitial_n5p191_n2p650 = 0;
174  mSumPFVsPtInitial_n2p650_n2p043 = 0;
175  mSumPFVsPtInitial_n2p043_n1p740 = 0;
176  mSumPFVsPtInitial_n1p740_n1p479 = 0;
177  mSumPFVsPtInitial_n1p479_n1p131 = 0;
178  mSumPFVsPtInitial_n1p131_n0p783 = 0;
179  mSumPFVsPtInitial_n0p783_n0p522 = 0;
180  mSumPFVsPtInitial_n0p522_0p522 = 0;
181  mSumPFVsPtInitial_0p522_0p783 = 0;
182  mSumPFVsPtInitial_0p783_1p131 = 0;
183  mSumPFVsPtInitial_1p131_1p479 = 0;
184  mSumPFVsPtInitial_1p479_1p740 = 0;
185  mSumPFVsPtInitial_1p740_2p043 = 0;
186  mSumPFVsPtInitial_2p043_2p650 = 0;
187  mSumPFVsPtInitial_2p650_5p191 = 0;
188 
189  mSumPFVsPt_n5p191_n2p650 = 0;
190  mSumPFVsPt_n2p650_n2p043 = 0;
191  mSumPFVsPt_n2p043_n1p740 = 0;
192  mSumPFVsPt_n1p740_n1p479 = 0;
193  mSumPFVsPt_n1p479_n1p131 = 0;
194  mSumPFVsPt_n1p131_n0p783 = 0;
195  mSumPFVsPt_n0p783_n0p522 = 0;
196  mSumPFVsPt_n0p522_0p522 = 0;
197  mSumPFVsPt_0p522_0p783 = 0;
198  mSumPFVsPt_0p783_1p131 = 0;
199  mSumPFVsPt_1p131_1p479 = 0;
200  mSumPFVsPt_1p479_1p740 = 0;
201  mSumPFVsPt_1p740_2p043 = 0;
202  mSumPFVsPt_2p043_2p650 = 0;
203  mSumPFVsPt_2p650_5p191 = 0;
204 
205  mSumPFPt_n5p191_n2p650 = 0;
206  mSumPFPt_n2p650_n2p043 = 0;
207  mSumPFPt_n2p043_n1p740 = 0;
208  mSumPFPt_n1p740_n1p479 = 0;
209  mSumPFPt_n1p479_n1p131 = 0;
210  mSumPFPt_n1p131_n0p783 = 0;
211  mSumPFPt_n0p783_n0p522 = 0;
212  mSumPFPt_n0p522_0p522 = 0;
213  mSumPFPt_0p522_0p783 = 0;
214  mSumPFPt_0p783_1p131 = 0;
215  mSumPFPt_1p131_1p479 = 0;
216  mSumPFPt_1p479_1p740 = 0;
217  mSumPFPt_1p740_2p043 = 0;
218  mSumPFPt_2p043_2p650 = 0;
219  mSumPFPt_2p650_5p191 = 0;
220 
221  mPFCandpT_vs_eta_Unknown = 0; // pf id 0
222  mPFCandpT_vs_eta_ChargedHadron = 0; // pf id - 1
223  mPFCandpT_vs_eta_electron = 0; // pf id - 2
224  mPFCandpT_vs_eta_muon = 0; // pf id - 3
225  mPFCandpT_vs_eta_photon = 0; // pf id - 4
226  mPFCandpT_vs_eta_NeutralHadron = 0; // pf id - 5
227  mPFCandpT_vs_eta_HadE_inHF = 0; // pf id - 6
228  mPFCandpT_vs_eta_EME_inHF = 0; // pf id - 7
229 
230  mPFCandpT_Barrel_Unknown = 0; // pf id 0
231  mPFCandpT_Barrel_ChargedHadron = 0; // pf id - 1
232  mPFCandpT_Barrel_electron = 0; // pf id - 2
233  mPFCandpT_Barrel_muon = 0; // pf id - 3
234  mPFCandpT_Barrel_photon = 0; // pf id - 4
235  mPFCandpT_Barrel_NeutralHadron = 0; // pf id - 5
236  mPFCandpT_Barrel_HadE_inHF = 0; // pf id - 6
237  mPFCandpT_Barrel_EME_inHF = 0; // pf id - 7
238 
239  mPFCandpT_Endcap_Unknown = 0; // pf id 0
240  mPFCandpT_Endcap_ChargedHadron = 0; // pf id - 1
241  mPFCandpT_Endcap_electron = 0; // pf id - 2
242  mPFCandpT_Endcap_muon = 0; // pf id - 3
243  mPFCandpT_Endcap_photon = 0; // pf id - 4
244  mPFCandpT_Endcap_NeutralHadron = 0; // pf id - 5
245  mPFCandpT_Endcap_HadE_inHF = 0; // pf id - 6
246  mPFCandpT_Endcap_EME_inHF = 0; // pf id - 7
247 
248  mPFCandpT_Forward_Unknown = 0; // pf id 0
249  mPFCandpT_Forward_ChargedHadron = 0; // pf id - 1
250  mPFCandpT_Forward_electron = 0; // pf id - 2
251  mPFCandpT_Forward_muon = 0; // pf id - 3
252  mPFCandpT_Forward_photon = 0; // pf id - 4
253  mPFCandpT_Forward_NeutralHadron = 0; // pf id - 5
254  mPFCandpT_Forward_HadE_inHF = 0; // pf id - 6
255  mPFCandpT_Forward_EME_inHF = 0; // pf id - 7
256 
257  }
258  if(isCaloJet){
259  mNCalopart = 0;
260  mCaloPt = 0;
261  mCaloEta = 0;
262  mCaloPhi = 0;
263  mCaloVsPt = 0;
264  mCaloVsPtInitial = 0;
265  mCaloArea = 0;
266 
267  mSumCaloVsPt = 0;
268  mSumCaloVsPtInitial = 0;
269  mSumCaloPt = 0;
270  mSumCaloVsPtInitial_eta = 0;
271  mSumCaloVsPt_eta = 0;
272  mSumCaloPt_eta = 0;
273  mSumSquaredCaloVsPt = 0;
274  mSumSquaredCaloVsPtInitial = 0;
275  mSumSquaredCaloPt = 0;
276  mSumSquaredCaloVsPtInitial_eta = 0;
277  mSumSquaredCaloVsPt_eta = 0;
278  mSumSquaredCaloPt_eta = 0;
279  mSumCaloVsPtInitial_HF = 0;
280  mSumCaloVsPt_HF = 0;
281  mSumCaloPt_HF = 0;
282 
283  mSumCaloVsPtInitial_n5p191_n2p650 = 0;
284  mSumCaloVsPtInitial_n2p650_n2p043 = 0;
285  mSumCaloVsPtInitial_n2p043_n1p740 = 0;
286  mSumCaloVsPtInitial_n1p740_n1p479 = 0;
287  mSumCaloVsPtInitial_n1p479_n1p131 = 0;
288  mSumCaloVsPtInitial_n1p131_n0p783 = 0;
289  mSumCaloVsPtInitial_n0p783_n0p522 = 0;
290  mSumCaloVsPtInitial_n0p522_0p522 = 0;
291  mSumCaloVsPtInitial_0p522_0p783 = 0;
292  mSumCaloVsPtInitial_0p783_1p131 = 0;
293  mSumCaloVsPtInitial_1p131_1p479 = 0;
294  mSumCaloVsPtInitial_1p479_1p740 = 0;
295  mSumCaloVsPtInitial_1p740_2p043 = 0;
296  mSumCaloVsPtInitial_2p043_2p650 = 0;
297  mSumCaloVsPtInitial_2p650_5p191 = 0;
298 
299  mSumCaloVsPt_n5p191_n2p650 = 0;
300  mSumCaloVsPt_n2p650_n2p043 = 0;
301  mSumCaloVsPt_n2p043_n1p740 = 0;
302  mSumCaloVsPt_n1p740_n1p479 = 0;
303  mSumCaloVsPt_n1p479_n1p131 = 0;
304  mSumCaloVsPt_n1p131_n0p783 = 0;
305  mSumCaloVsPt_n0p783_n0p522 = 0;
306  mSumCaloVsPt_n0p522_0p522 = 0;
307  mSumCaloVsPt_0p522_0p783 = 0;
308  mSumCaloVsPt_0p783_1p131 = 0;
309  mSumCaloVsPt_1p131_1p479 = 0;
310  mSumCaloVsPt_1p479_1p740 = 0;
311  mSumCaloVsPt_1p740_2p043 = 0;
312  mSumCaloVsPt_2p043_2p650 = 0;
313  mSumCaloVsPt_2p650_5p191 = 0;
314 
315  mSumCaloPt_n5p191_n2p650 = 0;
316  mSumCaloPt_n2p650_n2p043 = 0;
317  mSumCaloPt_n2p043_n1p740 = 0;
318  mSumCaloPt_n1p740_n1p479 = 0;
319  mSumCaloPt_n1p479_n1p131 = 0;
320  mSumCaloPt_n1p131_n0p783 = 0;
321  mSumCaloPt_n0p783_n0p522 = 0;
322  mSumCaloPt_n0p522_0p522 = 0;
323  mSumCaloPt_0p522_0p783 = 0;
324  mSumCaloPt_0p783_1p131 = 0;
325  mSumCaloPt_1p131_1p479 = 0;
326  mSumCaloPt_1p479_1p740 = 0;
327  mSumCaloPt_1p740_2p043 = 0;
328  mSumCaloPt_2p043_2p650 = 0;
329  mSumCaloPt_2p650_5p191 = 0;
330 
331  }
332 
333  mSumpt = 0;
334 
335  // Events variables
336  mNvtx = 0;
337  mHF = 0;
338 
339  // added Jan 12th 2015
340 
341  // Jet parameters
342  mEta = 0;
343  mPhi = 0;
344  mEnergy = 0;
345  mP = 0;
346  mPt = 0;
347  mMass = 0;
348  mConstituents = 0;
349  mJetArea = 0;
350  mjetpileup = 0;
351  mNJets_40 = 0;
352  mNJets = 0;
353 
354 
355 }
356 
357 void JetAnalyzer_HeavyIons::bookHistograms(DQMStore::IBooker & ibooker, edm::Run const & iRun,edm::EventSetup const &)
358 {
359 
360  ibooker.setCurrentFolder("JetMET/HIJetValidation/"+mInputCollection.label());
361 
362  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 };
363 
364  TH2F *h2D_etabins_vs_pt2 = new TH2F("h2D_etabins_vs_pt2",";#eta;sum p_{T}^{2}",etaBins_,edge_pseudorapidity,10000,0,10000);
365  TH2F *h2D_etabins_vs_pt = new TH2F("h2D_etabins_vs_pt",";#eta;sum p_{T}",etaBins_,edge_pseudorapidity,10000,-1000,1000);
366  TH2F *h2D_pfcand_etabins_vs_pt = new TH2F("h2D_etabins_vs_pt",";#eta;sum p_{T}",etaBins_,edge_pseudorapidity,300, 0, 300);
367 
368  if(isPFJet){
369 
370  mNPFpart = ibooker.book1D("NPFpart","No of particle flow candidates",1000,0,1000);
371  mPFPt = ibooker.book1D("PFPt","PF candidate p_{T}",10000,-500,500);
372  mPFEta = ibooker.book1D("PFEta","PF candidate #eta",120,-6,6);
373  mPFPhi = ibooker.book1D("PFPhi","PF candidate #phi",70,-3.5,3.5);
374  mPFVsPt = ibooker.book1D("PFVsPt","Vs PF candidate p_{T}",1000,-500,500);
375  mPFVsPtInitial = ibooker.book1D("PFVsPtInitial","Vs background subtracted PF candidate p_{T}",1000,-500,500);
376  mPFVsPtInitial_HF = ibooker.book2D("PFVsPtInitial_HF","Vs background subtracted PF candidate p_{T} versus HF", 600,0,6000,1000,-500,500);
377 
378  mPFVsPtInitial_Barrel_Centrality_100_95 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_100_95","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW100And95 Barrel",1000,-100,100);
379  mPFVsPtInitial_Barrel_Centrality_95_90 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_95_90","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW95And90 Barrel",1000,-100,100);
380  mPFVsPtInitial_Barrel_Centrality_90_85 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_90_85","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW90And85 Barrel",1000,-100,100);
381  mPFVsPtInitial_Barrel_Centrality_85_80 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_85_80","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW85And80 Barrel",1000,-100,100);
382  mPFVsPtInitial_Barrel_Centrality_80_75 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_80_75","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW80And75 Barrel",1000,-100,100);
383  mPFVsPtInitial_Barrel_Centrality_75_70 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_75_70","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW75And70 Barrel",1000,-100,100);
384  mPFVsPtInitial_Barrel_Centrality_70_65 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_70_65","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW70And65 Barrel",1000,-100,100);
385  mPFVsPtInitial_Barrel_Centrality_65_60 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_65_60","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW65And60 Barrel",1000,-100,100);
386  mPFVsPtInitial_Barrel_Centrality_60_55 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_60_55","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW60And55 Barrel",1000,-100,100);
387  mPFVsPtInitial_Barrel_Centrality_55_50 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_55_50","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW55And50 Barrel",1000,-100,100);
388  mPFVsPtInitial_Barrel_Centrality_50_45 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_50_45","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW50And45 Barrel",1000,-100,100);
389  mPFVsPtInitial_Barrel_Centrality_45_40 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_45_40","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW45And40 Barrel",1000,-100,100);
390  mPFVsPtInitial_Barrel_Centrality_40_35 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_40_35","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW40And35 Barrel",1000,-100,100);
391  mPFVsPtInitial_Barrel_Centrality_35_30 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_35_30","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW35And30 Barrel",1000,-100,100);
392  mPFVsPtInitial_Barrel_Centrality_30_25 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_30_25","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW30And25 Barrel",1000,-100,100);
393  mPFVsPtInitial_Barrel_Centrality_25_20 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_25_20","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW25And20 Barrel",1000,-100,100);
394  mPFVsPtInitial_Barrel_Centrality_20_15 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_20_15","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW20And15 Barrel",1000,-100,100);
395  mPFVsPtInitial_Barrel_Centrality_15_10 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_15_10","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW15And10 Barrel",1000,-100,100);
396  mPFVsPtInitial_Barrel_Centrality_10_5 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_10_5","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW10And05 Barrel",1000,-100,100);
397  mPFVsPtInitial_Barrel_Centrality_5_0 = ibooker.book1D("mPFVsPtInitial_Barrel_Centrality_5_0","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW05And00 Barrel",1000,-100,100);
398 
399  mPFVsPtInitial_EndCap_Centrality_100_95 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_100_95","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW100And95 EndCap",1000,-100,100);
400  mPFVsPtInitial_EndCap_Centrality_95_90 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_95_90","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW95And90 EndCap",1000,-100,100);
401  mPFVsPtInitial_EndCap_Centrality_90_85 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_90_85","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW90And85 EndCap",1000,-100,100);
402  mPFVsPtInitial_EndCap_Centrality_85_80 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_85_80","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW85And80 EndCap",1000,-100,100);
403  mPFVsPtInitial_EndCap_Centrality_80_75 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_80_75","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW80And75 EndCap",1000,-100,100);
404  mPFVsPtInitial_EndCap_Centrality_75_70 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_75_70","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW75And70 EndCap",1000,-100,100);
405  mPFVsPtInitial_EndCap_Centrality_70_65 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_70_65","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW70And65 EndCap",1000,-100,100);
406  mPFVsPtInitial_EndCap_Centrality_65_60 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_65_60","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW65And60 EndCap",1000,-100,100);
407  mPFVsPtInitial_EndCap_Centrality_60_55 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_60_55","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW60And55 EndCap",1000,-100,100);
408  mPFVsPtInitial_EndCap_Centrality_55_50 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_55_50","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW55And50 EndCap",1000,-100,100);
409  mPFVsPtInitial_EndCap_Centrality_50_45 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_50_45","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW50And45 EndCap",1000,-100,100);
410  mPFVsPtInitial_EndCap_Centrality_45_40 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_45_40","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW45And40 EndCap",1000,-100,100);
411  mPFVsPtInitial_EndCap_Centrality_40_35 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_40_35","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW40And35 EndCap",1000,-100,100);
412  mPFVsPtInitial_EndCap_Centrality_35_30 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_35_30","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW35And30 EndCap",1000,-100,100);
413  mPFVsPtInitial_EndCap_Centrality_30_25 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_30_25","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW30And25 EndCap",1000,-100,100);
414  mPFVsPtInitial_EndCap_Centrality_25_20 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_25_20","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW25And20 EndCap",1000,-100,100);
415  mPFVsPtInitial_EndCap_Centrality_20_15 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_20_15","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW20And15 EndCap",1000,-100,100);
416  mPFVsPtInitial_EndCap_Centrality_15_10 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_15_10","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW15And10 EndCap",1000,-100,100);
417  mPFVsPtInitial_EndCap_Centrality_10_5 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_10_5","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW10And05 EndCap",1000,-100,100);
418  mPFVsPtInitial_EndCap_Centrality_5_0 = ibooker.book1D("mPFVsPtInitial_EndCap_Centrality_5_0","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW05And00 EndCap",1000,-100,100);
419 
420  mPFVsPtInitial_HF_Centrality_100_95 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_100_95","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW100And95 HF",1000,-100,100);
421  mPFVsPtInitial_HF_Centrality_95_90 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_95_90","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW95And90 HF",1000,-100,100);
422  mPFVsPtInitial_HF_Centrality_90_85 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_90_85","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW90And85 HF",1000,-100,100);
423  mPFVsPtInitial_HF_Centrality_85_80 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_85_80","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW85And80 HF",1000,-100,100);
424  mPFVsPtInitial_HF_Centrality_80_75 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_80_75","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW80And75 HF",1000,-100,100);
425  mPFVsPtInitial_HF_Centrality_75_70 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_75_70","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW75And70 HF",1000,-100,100);
426  mPFVsPtInitial_HF_Centrality_70_65 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_70_65","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW70And65 HF",1000,-100,100);
427  mPFVsPtInitial_HF_Centrality_65_60 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_65_60","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW65And60 HF",1000,-100,100);
428  mPFVsPtInitial_HF_Centrality_60_55 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_60_55","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW60And55 HF",1000,-100,100);
429  mPFVsPtInitial_HF_Centrality_55_50 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_55_50","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW55And50 HF",1000,-100,100);
430  mPFVsPtInitial_HF_Centrality_50_45 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_50_45","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW50And45 HF",1000,-100,100);
431  mPFVsPtInitial_HF_Centrality_45_40 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_45_40","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW45And40 HF",1000,-100,100);
432  mPFVsPtInitial_HF_Centrality_40_35 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_40_35","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW40And35 HF",1000,-100,100);
433  mPFVsPtInitial_HF_Centrality_35_30 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_35_30","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW35And30 HF",1000,-100,100);
434  mPFVsPtInitial_HF_Centrality_30_25 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_30_25","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW30And25 HF",1000,-100,100);
435  mPFVsPtInitial_HF_Centrality_25_20 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_25_20","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW25And20 HF",1000,-100,100);
436  mPFVsPtInitial_HF_Centrality_20_15 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_20_15","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW20And15 HF",1000,-100,100);
437  mPFVsPtInitial_HF_Centrality_15_10 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_15_10","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW15And10 HF",1000,-100,100);
438  mPFVsPtInitial_HF_Centrality_10_5 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_10_5","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW10And05 HF",1000,-100,100);
439  mPFVsPtInitial_HF_Centrality_5_0 = ibooker.book1D("mPFVsPtInitial_HF_Centrality_5_0","Vs background subtracted PF candidate p_{T} versus HF Energy Dist. BTW05And00 HF",1000,-100,100);
440 
441  mPFArea = ibooker.book1D("PFArea","VS PF candidate area",100,0,4);
442  mPFDeltaR = ibooker.book1D("PFDeltaR","VS PF candidate DeltaR",100,0,4); //MZ
443  mPFDeltaR_Scaled_R = ibooker.book1D("PFDeltaR_Scaled_R","VS PF candidate DeltaR Divided by DeltaR square",100,0,4); //MZ
444  mPFDeltaR_pTCorrected = ibooker.book1D("PFDeltaR_pTCorrected","VS PF candidate DeltaR time pFpT over jet pT",100,0,4); //MZ
445  mPFDeltaR_pTCorrected_PFpT_20To30 = ibooker.book1D("PFDeltaR_pTCorrected_PFpT_20_30","VS PF candidate DeltaR time pFpT over jet pT for PFpT 20 to 30 GeV",100,0,4); //MZ
446  mPFDeltaR_pTCorrected_PFpT_30To50 = ibooker.book1D("PFDeltaR_pTCorrected_PFpT_30_50","VS PF candidate DeltaR time pFpT over jet pT for PFpT 30 to 50 GeV",100,0,4); //MZ
447  mPFDeltaR_pTCorrected_PFpT_50To80 = ibooker.book1D("PFDeltaR_pTCorrected_PFpT_50_80","VS PF candidate DeltaR time pFpT over jet pT for PFpT 50 to 80 GeV",100,0,4); //MZ
448  mPFDeltaR_pTCorrected_PFpT_80To120 = ibooker.book1D("PFDeltaR_pTCorrected_PFpT_80_120","VS PF candidate DeltaR time pFpT over jet pT for PFpT 80 to 120 GeV",100,0,4); //MZ
449  mPFDeltaR_pTCorrected_PFpT_120To180 = ibooker.book1D("PFDeltaR_pTCorrected_PFpT_120_180","VS PF candidate DeltaR time pFpT over jet pT for PFpT 120 to 180 GeV",100,0,4); //MZ
450  mPFDeltaR_pTCorrected_PFpT_180To300 = ibooker.book1D("PFDeltaR_pTCorrected_PFpT_180_300","VS PF candidate DeltaR time pFpT over jet pT for PFpT 180 to 300 GeV",100,0,4); //MZ
451  mPFDeltaR_pTCorrected_PFpT_300ToInf = ibooker.book1D("PFDeltaR_pTCorrected_PFpT_300_Inf","VS PF candidate DeltaR time pFpT over jet pT for PFpT 300 to Inf GeV",100,0,4); //MZ
452 
453  mPFDeltaR_pTCorrected_PFVsInitialpT_20To30 = ibooker.book1D("PFVsInitialDeltaR_pTCorrected_PFpT_20_30","VS PF candidate DeltaR time pFVsInitiaal pT over jet pT for PFpT 20 to 30 GeV",100,0,4); //MZ
454  mPFDeltaR_pTCorrected_PFVsInitialpT_30To50 = ibooker.book1D("PFVsInitialDeltaR_pTCorrected_PFpT_30_50","VS PF candidate DeltaR time pFVsInitiaal pT over jet pT for PFpT 30 to 50 GeV",100,0,4); //MZ
455  mPFDeltaR_pTCorrected_PFVsInitialpT_50To80 = ibooker.book1D("PFVsInitialDeltaR_pTCorrected_PFpT_50_80","VS PF candidate DeltaR time pFVsInitiaal pT over jet pT for PFpT 50 to 80 GeV",100,0,4); //MZ
456  mPFDeltaR_pTCorrected_PFVsInitialpT_80To120 = ibooker.book1D("PFVsInitialDeltaR_pTCorrected_PFpT_80_120","VS PF candidate DeltaR time pFVsInitiaal pT over jet pT for PFpT 80 to 120 GeV",100,0,4); //MZ
457  mPFDeltaR_pTCorrected_PFVsInitialpT_120To180 = ibooker.book1D("PFVsInitialDeltaR_pTCorrected_PFpT_120_180","VS PF candidate DeltaR time pFVsInitiaal pT over jet pT for PFpT 120 to 180 GeV",100,0,4); //MZ
458  mPFDeltaR_pTCorrected_PFVsInitialpT_180To300 = ibooker.book1D("PFVsInitialDeltaR_pTCorrected_PFpT_180_300","VS PF candidate DeltaR time pFVsInitiaal pT over jet pT for PFpT 180 to 300 GeV",100,0,4); //MZ
459  mPFDeltaR_pTCorrected_PFVsInitialpT_300ToInf = ibooker.book1D("PFVsInitialDeltaR_pTCorrected_PFpT_300_Inf","VS PF candidate DeltaR time pFVsInitiaal pT over jet pT for PFpT 300 to Inf GeV",100,0,4); //MZ
460 
461 
462  mPFVsPtInitialDeltaR_pTCorrected = ibooker.book1D("PFVsPtInitialDeltaR_pTCorrected","VS PF Initial candidate DeltaR time pFpTInitial over jet pT",100,0,4); //MZ
463  mPFVsPtDeltaR_pTCorrected = ibooker.book1D("PFVsPtDeltaR_pTCorrected","VS PF candidate DeltaR time pFpT over jet pT",100,0,4); //MZ
464 
465  mSumPFVsPt = ibooker.book1D("SumPFVsPt","Sum of final PF VS p_{T}",1000,-10000,10000);
466  mSumPFVsPtInitial= ibooker.book1D("SumPFVsPtInitial","Sum PF VS p_{T} after subtraction",1000,-10000,10000);
467  mSumPFPt = ibooker.book1D("SumPFPt","Sum of initial PF p_{T}",1000,-10000,10000);
468  mSumPFVsPt_eta = ibooker.book2D("SumPFVsPt_etaBins",h2D_etabins_vs_pt);
469  mSumPFVsPtInitial_eta = ibooker.book2D("SumPFVsPtInitial_etaBins",h2D_etabins_vs_pt);
470  mSumPFPt_eta = ibooker.book2D("SumPFPt_etaBins",h2D_etabins_vs_pt);
471 
472  mSumSquaredPFVsPt = ibooker.book1D("SumSquaredPFVsPt","Sum PF Vs p_{T} square",10000,0,10000);
473  mSumSquaredPFVsPtInitial= ibooker.book1D("SumSquaredPFVsPtInitial","Sum PF Vs p_{T} square after subtraction ",10000,0,10000);
474  mSumSquaredPFPt = ibooker.book1D("SumSquaredPFPt","Sum of initial PF p_{T} squared",10000,0,10000);
475  mSumSquaredPFVsPt_eta = ibooker.book2D("SumSquaredPFVsPt_etaBins",h2D_etabins_vs_pt2);
476  mSumSquaredPFVsPtInitial_eta = ibooker.book2D("SumSquaredPFVsPtInitial_etaBins",h2D_etabins_vs_pt2);
477  mSumSquaredPFPt_eta = ibooker.book2D("SumSquaredPFPt_etaBins",h2D_etabins_vs_pt2);
478 
479  mSumPFVsPtInitial_HF = ibooker.book2D("SumPFVsPtInitial_HF","HF Energy (y axis) vs Sum PF Vs p_{T} before subtraction (x axis)",1000,-1000,1000,1000,0,10000);
480  mSumPFVsPt_HF = ibooker.book2D("SumPFVsPt_HF","HF energy (y axis) vs Sum PF Vs p_{T} final (x axis)",1000,-1000,1000,1000,0,10000);
481  mSumPFPt_HF = ibooker.book2D("SumPFPt_HF","HF energy (y axis) vs Sum initial PF p_{T} (x axis)",1000,-1000,1000,1000,0,10000);
482 
483  mSumPFVsPtInitial_n5p191_n2p650 = ibooker.book1D("mSumPFVsPtInitial_n5p191_n2p650","Sum PFVsPt Initial variable in the eta range -5.191 to -2.650",1000,-5000,5000);
484  mSumPFVsPtInitial_n2p650_n2p043 = ibooker.book1D("mSumPFVsPtInitial_n2p650_n2p043","Sum PFVsPt Initial variable in the eta range -2.650 to -2.043 ",1000,-5000,5000);
485  mSumPFVsPtInitial_n2p043_n1p740 = ibooker.book1D("mSumPFVsPtInitial_n2p043_n1p740","Sum PFVsPt Initial variable in the eta range -2.043 to -1.740",1000,-1000,1000);
486  mSumPFVsPtInitial_n1p740_n1p479 = ibooker.book1D("mSumPFVsPtInitial_n1p740_n1p479","Sum PFVsPt Initial variable in the eta range -1.740 to -1.479",1000,-1000,1000);
487  mSumPFVsPtInitial_n1p479_n1p131 = ibooker.book1D("mSumPFVsPtInitial_n1p479_n1p131","Sum PFVsPt Initial variable in the eta range -1.479 to -1.131",1000,-1000,1000);
488  mSumPFVsPtInitial_n1p131_n0p783 = ibooker.book1D("mSumPFVsPtInitial_n1p131_n0p783","Sum PFVsPt Initial variable in the eta range -1.131 to -0.783",1000,-1000,1000);
489  mSumPFVsPtInitial_n0p783_n0p522 = ibooker.book1D("mSumPFVsPtInitial_n0p783_n0p522","Sum PFVsPt Initial variable in the eta range -0.783 to -0.522",1000,-1000,1000);
490  mSumPFVsPtInitial_n0p522_0p522 = ibooker.book1D("mSumPFVsPtInitial_n0p522_0p522","Sum PFVsPt Initial variable in the eta range -0.522 to 0.522",1000,-1000,1000);
491  mSumPFVsPtInitial_0p522_0p783 = ibooker.book1D("mSumPFVsPtInitial_0p522_0p783","Sum PFVsPt Initial variable in the eta range 0.522 to 0.783",1000,-1000,1000);
492  mSumPFVsPtInitial_0p783_1p131 = ibooker.book1D("mSumPFVsPtInitial_0p783_1p131","Sum PFVsPt Initial variable in the eta range 0.783 to 1.131",1000,-1000,1000);
493  mSumPFVsPtInitial_1p131_1p479 = ibooker.book1D("mSumPFVsPtInitial_1p131_1p479","Sum PFVsPt Initial variable in the eta range 1.131 to 1.479",1000,-1000,1000);
494  mSumPFVsPtInitial_1p479_1p740 = ibooker.book1D("mSumPFVsPtInitial_1p479_1p740","Sum PFVsPt Initial variable in the eta range 1.479 to 1.740",1000,-1000,1000);
495  mSumPFVsPtInitial_1p740_2p043 = ibooker.book1D("mSumPFVsPtInitial_1p740_2p043","Sum PFVsPt Initial variable in the eta range 1.740 to 2.043",1000,-1000,1000);
496  mSumPFVsPtInitial_2p043_2p650 = ibooker.book1D("mSumPFVsPtInitial_2p043_2p650","Sum PFVsPt Initial variable in the eta range 2.043 to 2.650",1000,-5000,5000);
497  mSumPFVsPtInitial_2p650_5p191 = ibooker.book1D("mSumPFVsPtInitial_2p650_5p191","Sum PFVsPt Initial variable in the eta range 2.650 to 5.191",1000,-5000,5000);
498 
499  mSumPFVsPt_n5p191_n2p650 = ibooker.book1D("mSumPFVsPt_n5p191_n2p650","Sum PFVsPt variable in the eta range -5.191 to -2.650",1000,-5000,5000);
500  mSumPFVsPt_n2p650_n2p043 = ibooker.book1D("mSumPFVsPt_n2p650_n2p043","Sum PFVsPt variable in the eta range -2.650 to -2.043 ",1000,-5000,5000);
501  mSumPFVsPt_n2p043_n1p740 = ibooker.book1D("mSumPFVsPt_n2p043_n1p740","Sum PFVsPt variable in the eta range -2.043 to -1.740",1000,-1000,1000);
502  mSumPFVsPt_n1p740_n1p479 = ibooker.book1D("mSumPFVsPt_n1p740_n1p479","Sum PFVsPt variable in the eta range -1.740 to -1.479",1000,-1000,1000);
503  mSumPFVsPt_n1p479_n1p131 = ibooker.book1D("mSumPFVsPt_n1p479_n1p131","Sum PFVsPt variable in the eta range -1.479 to -1.131",1000,-1000,1000);
504  mSumPFVsPt_n1p131_n0p783 = ibooker.book1D("mSumPFVsPt_n1p131_n0p783","Sum PFVsPt variable in the eta range -1.131 to -0.783",1000,-1000,1000);
505  mSumPFVsPt_n0p783_n0p522 = ibooker.book1D("mSumPFVsPt_n0p783_n0p522","Sum PFVsPt variable in the eta range -0.783 to -0.522",1000,-1000,1000);
506  mSumPFVsPt_n0p522_0p522 = ibooker.book1D("mSumPFVsPt_n0p522_0p522","Sum PFVsPt variable in the eta range -0.522 to 0.522",1000,-1000,1000);
507  mSumPFVsPt_0p522_0p783 = ibooker.book1D("mSumPFVsPt_0p522_0p783","Sum PFVsPt variable in the eta range 0.522 to 0.783",1000,-1000,1000);
508  mSumPFVsPt_0p783_1p131 = ibooker.book1D("mSumPFVsPt_0p783_1p131","Sum PFVsPt variable in the eta range 0.783 to 1.131",1000,-1000,1000);
509  mSumPFVsPt_1p131_1p479 = ibooker.book1D("mSumPFVsPt_1p131_1p479","Sum PFVsPt variable in the eta range 1.131 to 1.479",1000,-1000,1000);
510  mSumPFVsPt_1p479_1p740 = ibooker.book1D("mSumPFVsPt_1p479_1p740","Sum PFVsPt variable in the eta range 1.479 to 1.740",1000,-1000,1000);
511  mSumPFVsPt_1p740_2p043 = ibooker.book1D("mSumPFVsPt_1p740_2p043","Sum PFVsPt variable in the eta range 1.740 to 2.043",1000,-1000,1000);
512  mSumPFVsPt_2p043_2p650 = ibooker.book1D("mSumPFVsPt_2p043_2p650","Sum PFVsPt variable in the eta range 2.043 to 2.650",1000,-5000,5000);
513  mSumPFVsPt_2p650_5p191 = ibooker.book1D("mSumPFVsPt_2p650_5p191","Sum PFVsPt variable in the eta range 2.650 to 5.191",1000,-5000,5000);
514 
515  mSumPFPt_n5p191_n2p650 = ibooker.book1D("mSumPFPt_n5p191_n2p650","Sum PFPt in the eta range -5.191 to -2.650",1000,-5000,5000);
516  mSumPFPt_n2p650_n2p043 = ibooker.book1D("mSumPFPt_n2p650_n2p043","Sum PFPt in the eta range -2.650 to -2.043 ",1000,-5000,5000);
517  mSumPFPt_n2p043_n1p740 = ibooker.book1D("mSumPFPt_n2p043_n1p740","Sum PFPt in the eta range -2.043 to -1.740",1000,-1000,1000);
518  mSumPFPt_n1p740_n1p479 = ibooker.book1D("mSumPFPt_n1p740_n1p479","Sum PFPt in the eta range -1.740 to -1.479",1000,-1000,1000);
519  mSumPFPt_n1p479_n1p131 = ibooker.book1D("mSumPFPt_n1p479_n1p131","Sum PFPt in the eta range -1.479 to -1.131",1000,-1000,1000);
520  mSumPFPt_n1p131_n0p783 = ibooker.book1D("mSumPFPt_n1p131_n0p783","Sum PFPt in the eta range -1.131 to -0.783",1000,-1000,1000);
521  mSumPFPt_n0p783_n0p522 = ibooker.book1D("mSumPFPt_n0p783_n0p522","Sum PFPt in the eta range -0.783 to -0.522",1000,-1000,1000);
522  mSumPFPt_n0p522_0p522 = ibooker.book1D("mSumPFPt_n0p522_0p522","Sum PFPt in the eta range -0.522 to 0.522",1000,-1000,1000);
523  mSumPFPt_0p522_0p783 = ibooker.book1D("mSumPFPt_0p522_0p783","Sum PFPt in the eta range 0.522 to 0.783",1000,-1000,1000);
524  mSumPFPt_0p783_1p131 = ibooker.book1D("mSumPFPt_0p783_1p131","Sum PFPt in the eta range 0.783 to 1.131",1000,-1000,1000);
525  mSumPFPt_1p131_1p479 = ibooker.book1D("mSumPFPt_1p131_1p479","Sum PFPt in the eta range 1.131 to 1.479",1000,-1000,1000);
526  mSumPFPt_1p479_1p740 = ibooker.book1D("mSumPFPt_1p479_1p740","Sum PFPt in the eta range 1.479 to 1.740",1000,-1000,1000);
527  mSumPFPt_1p740_2p043 = ibooker.book1D("mSumPFPt_1p740_2p043","Sum PFPt in the eta range 1.740 to 2.043",1000,-1000,1000);
528  mSumPFPt_2p043_2p650 = ibooker.book1D("mSumPFPt_2p043_2p650","Sum PFPt in the eta range 2.043 to 2.650",1000,-5000,5000);
529  mSumPFPt_2p650_5p191 = ibooker.book1D("mSumPFPt_2p650_5p191","Sum PFPt in the eta range 2.650 to 5.191",1000,-5000,5000);
530 
531  mPFCandpT_vs_eta_Unknown = ibooker.book2D("PF_cand_X_unknown",h2D_pfcand_etabins_vs_pt); // pf id 0
532  mPFCandpT_vs_eta_ChargedHadron = ibooker.book2D("PF_cand_chargedHad",h2D_pfcand_etabins_vs_pt); // pf id - 1
533  mPFCandpT_vs_eta_electron = ibooker.book2D("PF_cand_electron",h2D_pfcand_etabins_vs_pt); // pf id - 2
534  mPFCandpT_vs_eta_muon = ibooker.book2D("PF_cand_muon",h2D_pfcand_etabins_vs_pt); // pf id - 3
535  mPFCandpT_vs_eta_photon = ibooker.book2D("PF_cand_photon",h2D_pfcand_etabins_vs_pt); // pf id - 4
536  mPFCandpT_vs_eta_NeutralHadron = ibooker.book2D("PF_cand_neutralHad",h2D_pfcand_etabins_vs_pt); // pf id - 5
537  mPFCandpT_vs_eta_HadE_inHF = ibooker.book2D("PF_cand_HadEner_inHF",h2D_pfcand_etabins_vs_pt); // pf id - 6
538  mPFCandpT_vs_eta_EME_inHF = ibooker.book2D("PF_cand_EMEner_inHF",h2D_pfcand_etabins_vs_pt); // pf id - 7
539 
540  mPFCandpT_Barrel_Unknown = ibooker.book1D("mPFCandpT_Barrel_Unknown",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 0
541  mPFCandpT_Barrel_ChargedHadron = ibooker.book1D("mPFCandpT_Barrel_ChargedHadron",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 1
542  mPFCandpT_Barrel_electron = ibooker.book1D("mPFCandpT_Barrel_electron",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 2
543  mPFCandpT_Barrel_muon = ibooker.book1D("mPFCandpT_Barrel_muon",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 3
544  mPFCandpT_Barrel_photon = ibooker.book1D("mPFCandpT_Barrel_photon",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 4
545  mPFCandpT_Barrel_NeutralHadron = ibooker.book1D("mPFCandpT_Barrel_NeutralHadron",Form(";PF candidate p_{T}, |#eta|<%2.2f; counts", BarrelEta),300, 0, 300); // pf id - 5
546  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
547  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
548 
549  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
550  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
551  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
552  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
553  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
554  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
555  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
556  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
557 
558  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
559  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
560  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
561  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
562  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
563  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
564  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
565  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
566 
567  }
568 
569  if(isCaloJet){
570 
571  mNCalopart = ibooker.book1D("NCalopart","No of particle flow candidates",1000,0,10000);
572  mCaloPt = ibooker.book1D("CaloPt","Calo candidate p_{T}",1000,-5000,5000);
573  mCaloEta = ibooker.book1D("CaloEta","Calo candidate #eta",120,-6,6);
574  mCaloPhi = ibooker.book1D("CaloPhi","Calo candidate #phi",70,-3.5,3.5);
575  mCaloVsPt = ibooker.book1D("CaloVsPt","Vs Calo candidate p_{T}",1000,-5000,5000);
576  mCaloVsPtInitial = ibooker.book1D("CaloVsPtInitial","Vs background subtracted Calo candidate p_{T}",1000,-500,500);
577  mCaloArea = ibooker.book1D("CaloArea","VS Calo candidate area",100,0,4);
578 
579  mSumCaloVsPt = ibooker.book1D("SumCaloVsPt","Sum of final Calo VS p_{T} ",1000,-10000,10000);
580  mSumCaloVsPtInitial= ibooker.book1D("SumCaloVsPtInitial","Sum Calo VS p_{T} after subtraction",1000,-10000,10000);
581  mSumCaloPt = ibooker.book1D("SumCaloPt","Sum Calo p_{T}",1000,-10000,10000);
582  mSumCaloVsPt_eta = ibooker.book2D("SumCaloVsPt_etaBins",h2D_etabins_vs_pt);
583  mSumCaloVsPtInitial_eta = ibooker.book2D("SumCaloVsPtInitial_etaBins",h2D_etabins_vs_pt);
584  mSumCaloPt_eta = ibooker.book2D("SumCaloPt_etaBins",h2D_etabins_vs_pt);
585 
586  mSumSquaredCaloVsPt = ibooker.book1D("SumSquaredCaloVsPt","Sum of final Calo VS p_{T} squared",10000,0,10000);
587  mSumSquaredCaloVsPtInitial= ibooker.book1D("SumSquaredCaloVsPtInitial","Sum of subtracted Calo VS p_{T} squared",10000,0,10000);
588  mSumSquaredCaloPt = ibooker.book1D("SumSquaredCaloPt","Sum of initial Calo tower p_{T} squared",10000,0,10000);
589  mSumSquaredCaloVsPt_eta = ibooker.book2D("SumSquaredCaloVsPt_etaBins",h2D_etabins_vs_pt2);
590  mSumSquaredCaloVsPtInitial_eta = ibooker.book2D("SumSquaredCaloVsPtInitial_etaBins",h2D_etabins_vs_pt2);
591  mSumSquaredCaloPt_eta = ibooker.book2D("SumSquaredCaloPt_etaBins",h2D_etabins_vs_pt2);
592 
593  mSumCaloVsPtInitial_HF = ibooker.book2D("SumCaloVsPtInitial_HF","HF Energy (y axis) vs Sum Calo Vs p_{T} before subtraction (x axis)",1000,-1000,1000,1000,0,10000);
594  mSumCaloVsPt_HF = ibooker.book2D("SumCaloVsPt_HF","HF Energy (y axis) vs Sum Calo Vs p_{T} (x axis)",1000,-1000,1000,1000,0,10000);
595  mSumCaloPt_HF = ibooker.book2D("SumCaloPt_HF","HF Energy (y axis) vs Sum Calo tower p_{T}",1000,-1000,1000,1000,0,10000);
596 
597  mSumCaloVsPtInitial_n5p191_n2p650 = ibooker.book1D("mSumCaloVsPtInitial_n5p191_n2p650","Sum CaloVsPt Initial variable in the eta range -5.191 to -2.650",1000,-5000,5000);
598  mSumCaloVsPtInitial_n2p650_n2p043 = ibooker.book1D("mSumCaloVsPtInitial_n2p650_n2p043","Sum CaloVsPt Initial variable in the eta range -2.650 to -2.043 ",1000,-5000,5000);
599  mSumCaloVsPtInitial_n2p043_n1p740 = ibooker.book1D("mSumCaloVsPtInitial_n2p043_n1p740","Sum CaloVsPt Initial variable in the eta range -2.043 to -1.740",1000,-1000,1000);
600  mSumCaloVsPtInitial_n1p740_n1p479 = ibooker.book1D("mSumCaloVsPtInitial_n1p740_n1p479","Sum CaloVsPt Initial variable in the eta range -1.740 to -1.479",1000,-1000,1000);
601  mSumCaloVsPtInitial_n1p479_n1p131 = ibooker.book1D("mSumCaloVsPtInitial_n1p479_n1p131","Sum CaloVsPt Initial variable in the eta range -1.479 to -1.131",1000,-1000,1000);
602  mSumCaloVsPtInitial_n1p131_n0p783 = ibooker.book1D("mSumCaloVsPtInitial_n1p131_n0p783","Sum CaloVsPt Initial variable in the eta range -1.131 to -0.783",1000,-1000,1000);
603  mSumCaloVsPtInitial_n0p783_n0p522 = ibooker.book1D("mSumCaloVsPtInitial_n0p783_n0p522","Sum CaloVsPt Initial variable in the eta range -0.783 to -0.522",1000,-1000,1000);
604  mSumCaloVsPtInitial_n0p522_0p522 = ibooker.book1D("mSumCaloVsPtInitial_n0p522_0p522","Sum CaloVsPt Initial variable in the eta range -0.522 to 0.522",1000,-1000,1000);
605  mSumCaloVsPtInitial_0p522_0p783 = ibooker.book1D("mSumCaloVsPtInitial_0p522_0p783","Sum CaloVsPt Initial variable in the eta range 0.522 to 0.783",1000,-1000,1000);
606  mSumCaloVsPtInitial_0p783_1p131 = ibooker.book1D("mSumCaloVsPtInitial_0p783_1p131","Sum CaloVsPt Initial variable in the eta range 0.783 to 1.131",1000,-1000,1000);
607  mSumCaloVsPtInitial_1p131_1p479 = ibooker.book1D("mSumCaloVsPtInitial_1p131_1p479","Sum CaloVsPt Initial variable in the eta range 1.131 to 1.479",1000,-1000,1000);
608  mSumCaloVsPtInitial_1p479_1p740 = ibooker.book1D("mSumCaloVsPtInitial_1p479_1p740","Sum CaloVsPt Initial variable in the eta range 1.479 to 1.740",1000,-1000,1000);
609  mSumCaloVsPtInitial_1p740_2p043 = ibooker.book1D("mSumCaloVsPtInitial_1p740_2p043","Sum CaloVsPt Initial variable in the eta range 1.740 to 2.043",1000,-1000,1000);
610  mSumCaloVsPtInitial_2p043_2p650 = ibooker.book1D("mSumCaloVsPtInitial_2p043_2p650","Sum CaloVsPt Initial variable in the eta range 2.043 to 2.650",1000,-5000,5000);
611  mSumCaloVsPtInitial_2p650_5p191 = ibooker.book1D("mSumCaloVsPtInitial_2p650_5p191","Sum CaloVsPt Initial variable in the eta range 2.650 to 5.191",1000,-5000,5000);
612 
613  mSumCaloVsPt_n5p191_n2p650 = ibooker.book1D("mSumCaloVsPt_n5p191_n2p650","Sum CaloVsPt variable in the eta range -5.191 to -2.650",1000,-5000,5000);
614  mSumCaloVsPt_n2p650_n2p043 = ibooker.book1D("mSumCaloVsPt_n2p650_n2p043","Sum CaloVsPt variable in the eta range -2.650 to -2.043",1000,-5000,5000);
615  mSumCaloVsPt_n2p043_n1p740 = ibooker.book1D("mSumCaloVsPt_n2p043_n1p740","Sum CaloVsPt variable in the eta range -2.043 to -1.740",1000,-1000,1000);
616  mSumCaloVsPt_n1p740_n1p479 = ibooker.book1D("mSumCaloVsPt_n1p740_n1p479","Sum CaloVsPt variable in the eta range -1.740 to -1.479",1000,-1000,1000);
617  mSumCaloVsPt_n1p479_n1p131 = ibooker.book1D("mSumCaloVsPt_n1p479_n1p131","Sum CaloVsPt variable in the eta range -1.479 to -1.131",1000,-1000,1000);
618  mSumCaloVsPt_n1p131_n0p783 = ibooker.book1D("mSumCaloVsPt_n1p131_n0p783","Sum CaloVsPt variable in the eta range -1.131 to -0.783",1000,-1000,1000);
619  mSumCaloVsPt_n0p783_n0p522 = ibooker.book1D("mSumCaloVsPt_n0p783_n0p522","Sum CaloVsPt variable in the eta range -0.783 to -0.522",1000,-1000,1000);
620  mSumCaloVsPt_n0p522_0p522 = ibooker.book1D("mSumCaloVsPt_n0p522_0p522","Sum CaloVsPt variable in the eta range -0.522 to 0.522",1000,-1000,1000);
621  mSumCaloVsPt_0p522_0p783 = ibooker.book1D("mSumCaloVsPt_0p522_0p783","Sum CaloVsPt variable in the eta range 0.522 to 0.783",1000,-1000,1000);
622  mSumCaloVsPt_0p783_1p131 = ibooker.book1D("mSumCaloVsPt_0p783_1p131","Sum CaloVsPt variable in the eta range 0.783 to 1.131",1000,-1000,1000);
623  mSumCaloVsPt_1p131_1p479 = ibooker.book1D("mSumCaloVsPt_1p131_1p479","Sum CaloVsPt variable in the eta range 1.131 to 1.479",1000,-1000,1000);
624  mSumCaloVsPt_1p479_1p740 = ibooker.book1D("mSumCaloVsPt_1p479_1p740","Sum CaloVsPt variable in the eta range 1.479 to 1.740",1000,-1000,1000);
625  mSumCaloVsPt_1p740_2p043 = ibooker.book1D("mSumCaloVsPt_1p740_2p043","Sum CaloVsPt variable in the eta range 1.740 to 2.043",1000,-1000,1000);
626  mSumCaloVsPt_2p043_2p650 = ibooker.book1D("mSumCaloVsPt_2p043_2p650","Sum CaloVsPt variable in the eta range 2.043 to 2.650",1000,-5000,5000);
627  mSumCaloVsPt_2p650_5p191 = ibooker.book1D("mSumCaloVsPt_2p650_5p191","Sum CaloVsPt variable in the eta range 2.650 to 5.191",1000,-5000,5000);
628 
629  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);
630  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);
631  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);
632  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);
633  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);
634  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);
635  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);
636  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);
637  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);
638  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);
639  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);
640  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);
641  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);
642  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);
643  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);
644 
645 
646  }
647 
648  // particle flow variables histograms
649  mSumpt = ibooker.book1D("SumpT","Sum p_{T} of all the PF candidates per event",1000,0,10000);
650 
651  // Event variables
652  mNvtx = ibooker.book1D("Nvtx", "number of vertices", 60, 0, 60);
653  mHF = ibooker.book1D("HF", "HF energy distribution",1000,0,10000);
654 
655  // Jet parameters
656  mEta = ibooker.book1D("Eta", "Eta", 120, -6, 6);
657  mPhi = ibooker.book1D("Phi", "Phi", 70, -3.5, 3.5);
658  mPt = ibooker.book1D("Pt", "Pt", 100, 0, 1000);
659  mP = ibooker.book1D("P", "P", 100, 0, 1000);
660  mEnergy = ibooker.book1D("Energy", "Energy", 100, 0, 1000);
661  mMass = ibooker.book1D("Mass", "Mass", 100, 0, 200);
662  mConstituents = ibooker.book1D("Constituents", "Constituents", 100, 0, 100);
663  mJetArea = ibooker.book1D("JetArea", "JetArea", 100, 0, 4);
664  mjetpileup = ibooker.book1D("jetPileUp","jetPileUp",100,0,150);
665  mNJets_40 = ibooker.book1D("NJets_pt_greater_40", "NJets pT > 40 GeV", 100, 0, 100);
666  mNJets = ibooker.book1D("NJets", "NJets", 100, 0, 50);
667 
668  if (mOutputFile.empty ())
669  LogInfo("OutputInfo") << " Histograms will NOT be saved";
670  else
671  LogInfo("OutputInfo") << " Histograms will be saved to file:" << mOutputFile;
672 
673  delete h2D_etabins_vs_pt2;
674  delete h2D_etabins_vs_pt;
675  delete h2D_pfcand_etabins_vs_pt;
676 
677 }
678 
679 
680 
681 //------------------------------------------------------------------------------
682 // ~JetAnalyzer_HeavyIons
683 //------------------------------------------------------------------------------
684 JetAnalyzer_HeavyIons::~JetAnalyzer_HeavyIons() {}
685 
686 
687 //------------------------------------------------------------------------------
688 // beginJob
689 //------------------------------------------------------------------------------
690 //void JetAnalyzer_HeavyIons::beginJob() {
691 // std:://cout<<"inside the begin job function"<<endl;
692 //}
693 
694 
695 //------------------------------------------------------------------------------
696 // endJob
697 //------------------------------------------------------------------------------
698 //void JetAnalyzer_HeavyIons::endJob()
699 //{
700 // if (!mOutputFile.empty() && &*edm::Service<DQMStore>())
701 // {
702 // edm::Service<DQMStore>()->save(mOutputFile);
703 // }
704 //}
705 
706 
707 //------------------------------------------------------------------------------
708 // analyze
709 //------------------------------------------------------------------------------
710 void JetAnalyzer_HeavyIons::analyze(const edm::Event& mEvent, const edm::EventSetup& mSetup)
711 {
712 
713  // switch(mEvent.id().event() == 15296770)
714  // case 1:
715  // break;
716 
717  // Get the primary vertices
718  edm::Handle<vector<reco::Vertex> > pvHandle;
719  mEvent.getByToken(pvToken_, pvHandle);
720  reco::Vertex::Point vtx(0,0,0);
721  edm::Handle<reco::VertexCollection> vtxs;
722 
723  mEvent.getByToken(hiVertexToken_, vtxs);
724  int greatestvtx = 0;
725  int nVertex = vtxs->size();
726 
727  for (unsigned int i = 0 ; i< vtxs->size(); ++i){
728  unsigned int daughter = (*vtxs)[i].tracksSize();
729  if( daughter > (*vtxs)[greatestvtx].tracksSize()) greatestvtx = i;
730  }
731 
732  if(nVertex<=0){
733  vtx = reco::Vertex::Point(0,0,0);
734  }
735  vtx = (*vtxs)[greatestvtx].position();
736 
737  int nGoodVertices = 0;
738 
739  if (pvHandle.isValid())
740  {
741  for (unsigned i=0; i<pvHandle->size(); i++)
742  {
743  if ((*pvHandle)[i].ndof() > 4 &&
744  (fabs((*pvHandle)[i].z()) <= 24) &&
745  (fabs((*pvHandle)[i].position().rho()) <= 2))
746  nGoodVertices++;
747  }
748  }
749 
750  mNvtx->Fill(nGoodVertices);
751 
752  // Get the Jet collection
753  //----------------------------------------------------------------------------
754 
755  std::vector<Jet> recoJets;
756  recoJets.clear();
757 
758  edm::Handle<CaloJetCollection> caloJets;
759  edm::Handle<JPTJetCollection> jptJets;
760  edm::Handle<PFJetCollection> pfJets;
761  edm::Handle<BasicJetCollection> basicJets;
762 
763  // Get the Particle flow candidates and the Voronoi variables
764  edm::Handle<reco::PFCandidateCollection> pfCandidates;
765  edm::Handle<CaloTowerCollection> caloCandidates;
766  edm::Handle<reco::CandidateView> pfcandidates_;
767  edm::Handle<reco::CandidateView> calocandidates_;
768 
769  edm::Handle<edm::ValueMap<VoronoiBackground>> VsBackgrounds;
770  edm::Handle<std::vector<float>> vn_;
771 
772  if (isCaloJet) mEvent.getByToken(caloJetsToken_, caloJets);
773  if (isJPTJet) mEvent.getByToken(jptJetsToken_, jptJets);
774  if (isPFJet) {
775  if(std::string("Pu")==UEAlgo) mEvent.getByToken(basicJetsToken_, basicJets);
776  if(std::string("Vs")==UEAlgo) mEvent.getByToken(pfJetsToken_, pfJets);
777  }
778 
779  mEvent.getByToken(pfCandToken_, pfCandidates);
780  mEvent.getByToken(pfCandViewToken_, pfcandidates_);
781 
782  mEvent.getByToken(caloTowersToken_, caloCandidates);
783  mEvent.getByToken(caloCandViewToken_, calocandidates_);
784 
785  mEvent.getByToken(backgrounds_, VsBackgrounds);
786  mEvent.getByToken(backgrounds_value_, vn_);
787 
788  // get the centrality
789  edm::Handle<reco::Centrality> cent;
790  mEvent.getByToken(centralityToken, cent); //_centralitytag comes from the cfg
791 
792  mHF->Fill(cent->EtHFtowerSum());
793  Float_t HF_energy = cent->EtHFtowerSum();
794 
795  edm::Handle<int> cbin;
796  mEvent.getByToken(centralityBinToken, cbin);
797 
798  if (!cent.isValid()) return;
799  int hibin = -999;
800  if(cent.isValid())
801  hibin = *cbin;
802 
803 
804  const reco::PFCandidateCollection *pfCandidateColl = pfCandidates.product();
805 
806  Float_t vsPt=0;
807  Float_t vsPtInitial = 0;
808  Float_t vsArea = 0;
809  Int_t NPFpart = 0;
810  Int_t NCaloTower = 0;
811  Float_t pfPt = 0;
812  Float_t pfEta = 0;
813  Float_t pfPhi = 0;
814  Int_t pfID = 0;
815  Float_t pfDeltaR = 0;
816  Float_t caloPt = 0;
817  Float_t caloEta = 0;
818  Float_t caloPhi = 0;
819  Float_t SumPt_value = 0;
820 
821  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 };
822 
823 
824  vector < vector <float> > numbers;
825  vector <float> tempVector;
826  numbers.clear();
827  tempVector.clear();
828 
829  if(isCaloJet){
830 
831  Float_t SumCaloVsPtInitial[etaBins_];
832  Float_t SumCaloVsPt[etaBins_];
833  Float_t SumCaloPt[etaBins_];
834 
835  Float_t SumSquaredCaloVsPtInitial[etaBins_];
836  Float_t SumSquaredCaloVsPt[etaBins_];
837  Float_t SumSquaredCaloPt[etaBins_];
838 
839  // Need to set up histograms to get the RMS values for each pT bin
840  TH1F *hSumCaloVsPtInitial[nedge_pseudorapidity-1], *hSumCaloVsPt[nedge_pseudorapidity-1], *hSumCaloPt[nedge_pseudorapidity-1];
841 
842  for(int i = 0;i<etaBins_;++i)
843  {
844 
845  SumCaloVsPtInitial[i] = 0;
846  SumCaloVsPt[i] = 0;
847  SumCaloPt[i] = 0;
848  SumSquaredCaloVsPtInitial[i] = 0;
849  SumSquaredCaloVsPt[i] = 0;
850  SumSquaredCaloPt[i] = 0;
851 
852  hSumCaloVsPtInitial[i] = new TH1F(Form("hSumCaloVsPtInitial_%d",i),"",10000,-10000,10000);
853  hSumCaloVsPt[i] = new TH1F(Form("hSumCaloVsPt_%d",i),"",10000,-10000,10000);
854  hSumCaloPt[i] = new TH1F(Form("hSumCaloPt_%d",i),"",10000,-10000,10000);
855 
856  }
857 
858  for(unsigned icand = 0;icand<caloCandidates->size(); icand++){
859 
860  const CaloTower & tower = (*caloCandidates)[icand];
861  reco::CandidateViewRef ref(calocandidates_,icand);
862  //10 is tower pT min
863  if(tower.p4(vtx).Et() < 0.1) continue;
864 
865  vsPt = 0;
866  vsPtInitial = 0;
867  vsArea = 0;
868 
869  if(std::string("Vs")==UEAlgo) {
870  const reco::VoronoiBackground& voronoi = (*VsBackgrounds)[ref];
871  vsPt = voronoi.pt();
872  vsPtInitial = voronoi.pt_subtracted();
873  vsArea = voronoi.area();
874  }
875 
876  NCaloTower++;
877 
878  caloPt = tower.p4(vtx).Et();
879  caloEta = tower.p4(vtx).Eta();
880  caloPhi = tower.p4(vtx).Phi();
881 
882 
883  for(size_t k = 0;k<nedge_pseudorapidity-1; k++){
884  if(caloEta >= edge_pseudorapidity[k] && caloEta < edge_pseudorapidity[k+1]){
885  SumCaloVsPtInitial[k] = SumCaloVsPtInitial[k] + vsPtInitial;
886  SumCaloVsPt[k] = SumCaloVsPt[k] + vsPt;
887  SumCaloPt[k] = SumCaloPt[k] + caloPt;
888  SumSquaredCaloVsPtInitial[k] = SumSquaredCaloVsPtInitial[k] + vsPtInitial*vsPtInitial;
889  SumSquaredCaloVsPt[k] = SumSquaredCaloVsPt[k] + vsPt*vsPt;
890  SumSquaredCaloPt[k] = SumSquaredCaloPt[k] + caloPt*caloPt;
891  break;
892  }// eta selection statement
893 
894  }// eta bin loop
895 
896  SumPt_value = SumPt_value + caloPt;
897 
898  mCaloPt->Fill(caloPt);
899  mCaloEta->Fill(caloEta);
900  mCaloPhi->Fill(caloPhi);
901  mCaloVsPt->Fill(vsPt);
902  mCaloVsPtInitial->Fill(vsPtInitial);
903  mCaloArea->Fill(vsArea);
904 
905  }// calo tower candidate loop
906 
907  for(int k = 0;k<nedge_pseudorapidity-1;k++){
908 
909  hSumCaloVsPtInitial[k]->Fill(SumCaloVsPtInitial[k]);
910  hSumCaloVsPt[k]->Fill(SumCaloVsPt[k]);
911  hSumCaloPt[k]->Fill(SumCaloPt[k]);
912 
913  }// eta bin loop
914 
915  Float_t Evt_SumCaloVsPt = 0;
916  Float_t Evt_SumCaloVsPtInitial = 0;
917  Float_t Evt_SumCaloPt = 0;
918 
919  Float_t Evt_SumSquaredCaloVsPt = 0;
920  Float_t Evt_SumSquaredCaloVsPtInitial = 0;
921  Float_t Evt_SumSquaredCaloPt = 0;
922 
923 
924  mSumCaloVsPtInitial_n5p191_n2p650->Fill(SumCaloVsPtInitial[0]);
925  mSumCaloVsPtInitial_n2p650_n2p043->Fill(SumCaloVsPtInitial[1]);
926  mSumCaloVsPtInitial_n2p043_n1p740->Fill(SumCaloVsPtInitial[2]);
927  mSumCaloVsPtInitial_n1p740_n1p479->Fill(SumCaloVsPtInitial[3]);
928  mSumCaloVsPtInitial_n1p479_n1p131->Fill(SumCaloVsPtInitial[4]);
929  mSumCaloVsPtInitial_n1p131_n0p783->Fill(SumCaloVsPtInitial[5]);
930  mSumCaloVsPtInitial_n0p783_n0p522->Fill(SumCaloVsPtInitial[6]);
931  mSumCaloVsPtInitial_n0p522_0p522->Fill(SumCaloVsPtInitial[7]);
932  mSumCaloVsPtInitial_0p522_0p783->Fill(SumCaloVsPtInitial[8]);
933  mSumCaloVsPtInitial_0p783_1p131->Fill(SumCaloVsPtInitial[9]);
934  mSumCaloVsPtInitial_1p131_1p479->Fill(SumCaloVsPtInitial[10]);
935  mSumCaloVsPtInitial_1p479_1p740->Fill(SumCaloVsPtInitial[11]);
936  mSumCaloVsPtInitial_1p740_2p043->Fill(SumCaloVsPtInitial[12]);
937  mSumCaloVsPtInitial_2p043_2p650->Fill(SumCaloVsPtInitial[13]);
938  mSumCaloVsPtInitial_2p650_5p191->Fill(SumCaloVsPtInitial[14]);
939 
940  mSumCaloVsPt_n5p191_n2p650->Fill(SumCaloVsPt[0]);
941  mSumCaloVsPt_n2p650_n2p043->Fill(SumCaloVsPt[1]);
942  mSumCaloVsPt_n2p043_n1p740->Fill(SumCaloVsPt[2]);
943  mSumCaloVsPt_n1p740_n1p479->Fill(SumCaloVsPt[3]);
944  mSumCaloVsPt_n1p479_n1p131->Fill(SumCaloVsPt[4]);
945  mSumCaloVsPt_n1p131_n0p783->Fill(SumCaloVsPt[5]);
946  mSumCaloVsPt_n0p783_n0p522->Fill(SumCaloVsPt[6]);
947  mSumCaloVsPt_n0p522_0p522->Fill(SumCaloVsPt[7]);
948  mSumCaloVsPt_0p522_0p783->Fill(SumCaloVsPt[8]);
949  mSumCaloVsPt_0p783_1p131->Fill(SumCaloVsPt[9]);
950  mSumCaloVsPt_1p131_1p479->Fill(SumCaloVsPt[10]);
951  mSumCaloVsPt_1p479_1p740->Fill(SumCaloVsPt[11]);
952  mSumCaloVsPt_1p740_2p043->Fill(SumCaloVsPt[12]);
953  mSumCaloVsPt_2p043_2p650->Fill(SumCaloVsPt[13]);
954  mSumCaloVsPt_2p650_5p191->Fill(SumCaloVsPt[14]);
955 
956  mSumCaloPt_n5p191_n2p650->Fill(SumCaloPt[0]);
957  mSumCaloPt_n2p650_n2p043->Fill(SumCaloPt[1]);
958  mSumCaloPt_n2p043_n1p740->Fill(SumCaloPt[2]);
959  mSumCaloPt_n1p740_n1p479->Fill(SumCaloPt[3]);
960  mSumCaloPt_n1p479_n1p131->Fill(SumCaloPt[4]);
961  mSumCaloPt_n1p131_n0p783->Fill(SumCaloPt[5]);
962  mSumCaloPt_n0p783_n0p522->Fill(SumCaloPt[6]);
963  mSumCaloPt_n0p522_0p522->Fill(SumCaloPt[7]);
964  mSumCaloPt_0p522_0p783->Fill(SumCaloPt[8]);
965  mSumCaloPt_0p783_1p131->Fill(SumCaloPt[9]);
966  mSumCaloPt_1p131_1p479->Fill(SumCaloPt[10]);
967  mSumCaloPt_1p479_1p740->Fill(SumCaloPt[11]);
968  mSumCaloPt_1p740_2p043->Fill(SumCaloPt[12]);
969  mSumCaloPt_2p043_2p650->Fill(SumCaloPt[13]);
970  mSumCaloPt_2p650_5p191->Fill(SumCaloPt[14]);
971 
972 
973  for(size_t k = 0;k<nedge_pseudorapidity-1;k++){
974 
975  Evt_SumCaloVsPtInitial = Evt_SumCaloVsPtInitial + SumCaloVsPtInitial[k];
976  Evt_SumCaloVsPt = Evt_SumCaloVsPt + SumCaloVsPt[k];
977  Evt_SumCaloPt = Evt_SumCaloPt + SumCaloPt[k];
978 
979  mSumCaloVsPtInitial_eta->Fill(edge_pseudorapidity[k],SumCaloVsPtInitial[k]);
980  mSumCaloVsPt_eta->Fill(edge_pseudorapidity[k],SumCaloVsPt[k]);
981  mSumCaloPt_eta->Fill(edge_pseudorapidity[k],SumCaloPt[k]);
982 
983  Evt_SumSquaredCaloVsPtInitial = Evt_SumSquaredCaloVsPtInitial + SumSquaredCaloVsPtInitial[k];
984  Evt_SumSquaredCaloVsPt = Evt_SumSquaredCaloVsPt + SumSquaredCaloVsPt[k];
985  Evt_SumSquaredCaloPt = Evt_SumSquaredCaloPt + SumSquaredCaloPt[k];
986 
987  mSumSquaredCaloVsPtInitial_eta->Fill(edge_pseudorapidity[k],hSumCaloVsPtInitial[k]->GetRMS(1));
988  mSumSquaredCaloVsPt_eta->Fill(edge_pseudorapidity[k],hSumCaloVsPt[k]->GetRMS(1));
989  mSumSquaredCaloPt_eta->Fill(edge_pseudorapidity[k],hSumCaloPt[k]->GetRMS(1));
990 
991  delete hSumCaloVsPtInitial[k];
992  delete hSumCaloVsPt[k];
993  delete hSumCaloPt[k];
994 
995  }// eta bin loop
996 
997  mSumCaloVsPtInitial->Fill(Evt_SumCaloVsPtInitial);
998  mSumCaloVsPt->Fill(Evt_SumCaloVsPt);
999  mSumCaloPt->Fill(Evt_SumCaloPt);
1000  mSumCaloVsPtInitial_HF->Fill(Evt_SumCaloVsPtInitial,HF_energy);
1001  mSumCaloVsPt_HF->Fill(Evt_SumCaloVsPt,HF_energy);
1002  mSumCaloPt_HF->Fill(Evt_SumCaloPt,HF_energy);
1003 
1004  mSumSquaredCaloVsPtInitial->Fill(Evt_SumSquaredCaloVsPtInitial);
1005  mSumSquaredCaloVsPt->Fill(Evt_SumSquaredCaloVsPt);
1006  mSumSquaredCaloPt->Fill(Evt_SumSquaredCaloPt);
1007 
1008  mNCalopart->Fill(NCaloTower);
1009  mSumpt->Fill(SumPt_value);
1010 
1011  }// is calo jet
1012 
1013  if(isPFJet){
1014 
1015  Float_t SumPFVsPtInitial[etaBins_];
1016  Float_t SumPFVsPt[etaBins_];
1017  Float_t SumPFPt[etaBins_];
1018 
1019  Float_t SumSquaredPFVsPtInitial[etaBins_];
1020  Float_t SumSquaredPFVsPt[etaBins_];
1021  Float_t SumSquaredPFPt[etaBins_];
1022 
1023  // Need to set up histograms to get the RMS values for each pT bin
1024  TH1F *hSumPFVsPtInitial[nedge_pseudorapidity-1], *hSumPFVsPt[nedge_pseudorapidity-1], *hSumPFPt[nedge_pseudorapidity-1];
1025 
1026  for(int i = 0;i<etaBins_;i++){
1027 
1028  SumPFVsPtInitial[i] = 0;
1029  SumPFVsPt[i] = 0;
1030  SumPFPt[i] = 0;
1031  SumSquaredPFVsPtInitial[i] = 0;
1032  SumSquaredPFVsPt[i] = 0;
1033  SumSquaredPFPt[i] = 0;
1034 
1035  hSumPFVsPtInitial[i] = new TH1F(Form("hSumPFVsPtInitial_%d",i),"",10000,-10000,10000);
1036  hSumPFVsPt[i] = new TH1F(Form("hSumPFVsPt_%d",i),"",10000,-10000,10000);
1037  hSumPFPt[i] = new TH1F(Form("hSumPFPt_%d",i),"",10000,-10000,10000);
1038 
1039  }
1040 
1041  vector<vector <float> > PF_Space(1,vector<float>(3));
1042 
1043  for(unsigned icand=0;icand<pfCandidateColl->size(); icand++)
1044  {
1045  const reco::PFCandidate pfCandidate = pfCandidateColl->at(icand);
1046  reco::CandidateViewRef ref(pfcandidates_,icand);
1047  if(pfCandidate.pt() < 5) continue;
1048 
1049  if(std::string("Vs")==UEAlgo)
1050  {
1051 
1052  const reco::VoronoiBackground& voronoi = (*VsBackgrounds)[ref];
1053  vsPt = voronoi.pt();
1054  vsPtInitial = voronoi.pt_subtracted();
1055  vsArea = voronoi.area();
1056 
1057  }
1058 
1059  NPFpart++;
1060  pfPt = pfCandidate.pt();
1061  pfEta = pfCandidate.eta();
1062  pfPhi = pfCandidate.phi();
1063  pfID = pfCandidate.particleId();
1064 
1065  bool isBarrel = false;
1066  bool isEndcap = false;
1067  bool isForward = false;
1068 
1069  if(fabs(pfEta)<BarrelEta) isBarrel = true;
1070  if(fabs(pfEta)>=BarrelEta && fabs(pfEta)<EndcapEta) isEndcap = true;
1071  if(fabs(pfEta)>=EndcapEta && fabs(pfEta)<ForwardEta) isForward = true;
1072 
1073  switch(pfID){
1074  case 0 :
1075  mPFCandpT_vs_eta_Unknown->Fill(pfPt, pfEta);
1076  if(isBarrel) mPFCandpT_Barrel_Unknown->Fill(pfPt);
1077  if(isEndcap) mPFCandpT_Endcap_Unknown->Fill(pfPt);
1078  if(isForward) mPFCandpT_Forward_Unknown->Fill(pfPt);
1079  case 1 :
1080  mPFCandpT_vs_eta_ChargedHadron->Fill(pfPt, pfEta);
1081  if(isBarrel) mPFCandpT_Barrel_ChargedHadron->Fill(pfPt);
1082  if(isEndcap) mPFCandpT_Endcap_ChargedHadron->Fill(pfPt);
1083  if(isForward) mPFCandpT_Forward_ChargedHadron->Fill(pfPt);
1084  case 2 :
1085  mPFCandpT_vs_eta_electron->Fill(pfPt, pfEta);
1086  if(isBarrel) mPFCandpT_Barrel_electron->Fill(pfPt);
1087  if(isEndcap) mPFCandpT_Endcap_electron->Fill(pfPt);
1088  if(isForward) mPFCandpT_Forward_electron->Fill(pfPt);
1089  case 3 :
1090  mPFCandpT_vs_eta_muon->Fill(pfPt, pfEta);
1091  if(isBarrel) mPFCandpT_Barrel_muon->Fill(pfPt);
1092  if(isEndcap) mPFCandpT_Endcap_muon->Fill(pfPt);
1093  if(isForward) mPFCandpT_Forward_muon->Fill(pfPt);
1094  case 4 :
1095  mPFCandpT_vs_eta_photon->Fill(pfPt, pfEta);
1096  if(isBarrel) mPFCandpT_Barrel_photon->Fill(pfPt);
1097  if(isEndcap) mPFCandpT_Endcap_photon->Fill(pfPt);
1098  if(isForward) mPFCandpT_Forward_photon->Fill(pfPt);
1099  case 5 :
1100  mPFCandpT_vs_eta_NeutralHadron->Fill(pfPt, pfEta);
1101  if(isBarrel) mPFCandpT_Barrel_NeutralHadron->Fill(pfPt);
1102  if(isEndcap) mPFCandpT_Endcap_NeutralHadron->Fill(pfPt);
1103  if(isForward) mPFCandpT_Forward_NeutralHadron->Fill(pfPt);
1104  case 6 :
1105  mPFCandpT_vs_eta_HadE_inHF->Fill(pfPt, pfEta);
1106  if(isBarrel) mPFCandpT_Barrel_HadE_inHF->Fill(pfPt);
1107  if(isEndcap) mPFCandpT_Endcap_HadE_inHF->Fill(pfPt);
1108  if(isForward) mPFCandpT_Forward_HadE_inHF->Fill(pfPt);
1109  case 7 :
1110  mPFCandpT_vs_eta_EME_inHF->Fill(pfPt, pfEta);
1111  if(isBarrel) mPFCandpT_Barrel_EME_inHF->Fill(pfPt);
1112  if(isEndcap) mPFCandpT_Endcap_EME_inHF->Fill(pfPt);
1113  if(isForward) mPFCandpT_Forward_EME_inHF->Fill(pfPt);
1114  }
1115 
1116 
1117  //Fill 2d vector matrix
1118  tempVector.push_back(pfPt);
1119  tempVector.push_back(pfEta);
1120  tempVector.push_back(pfPhi);
1121  tempVector.push_back(vsPtInitial);
1122  tempVector.push_back(vsPt);
1123 
1124 
1125  numbers.push_back(tempVector);
1126  tempVector.clear();
1127 
1128  for(size_t k = 0;k<nedge_pseudorapidity-1; k++)
1129  {
1130  if(pfEta >= edge_pseudorapidity[k] && pfEta < edge_pseudorapidity[k+1]){
1131  SumPFVsPtInitial[k] = SumPFVsPtInitial[k] + vsPtInitial;
1132  SumPFVsPt[k] = SumPFVsPt[k] + vsPt;
1133  SumPFPt[k] = SumPFPt[k] + pfPt;
1134 
1135  SumSquaredPFVsPtInitial[k] = SumSquaredPFVsPtInitial[k] + vsPtInitial*vsPtInitial;
1136  SumSquaredPFVsPt[k] = SumSquaredPFVsPt[k] + vsPt*vsPt;
1137  SumSquaredPFPt[k] = SumSquaredPFPt[k] + pfPt*pfPt;
1138  break;
1139  }// eta selection statement
1140 
1141  }// eta bin loop
1142 
1143  SumPt_value = SumPt_value + pfPt;
1144 
1145  mPFPt->Fill(pfPt);
1146  mPFEta->Fill(pfEta);
1147  mPFPhi->Fill(pfPhi);
1148  mPFVsPt->Fill(vsPt);
1149  mPFVsPtInitial->Fill(vsPtInitial); //filled with zeros (need a fix)
1150  mPFVsPtInitial_HF->Fill(cent->EtHFtowerSum(), vsPtInitial,1); //filled with zeros (need a fix)
1151 
1152  if(isBarrel)
1153  {
1154  if(hibin >= 190 && hibin < 200)
1155  mPFVsPtInitial_Barrel_Centrality_100_95->Fill(vsPtInitial);
1156 
1157  if(hibin >= 180 && hibin < 190)
1158  mPFVsPtInitial_Barrel_Centrality_95_90->Fill(vsPtInitial);
1159 
1160  if(hibin >= 170 && hibin < 180)
1161  mPFVsPtInitial_Barrel_Centrality_90_85->Fill(vsPtInitial);
1162 
1163  if(hibin >= 160 && hibin < 170)
1164  mPFVsPtInitial_Barrel_Centrality_85_80->Fill(vsPtInitial);
1165 
1166  if(hibin >= 150 && hibin < 160)
1167  mPFVsPtInitial_Barrel_Centrality_80_75->Fill(vsPtInitial);
1168 
1169  if(hibin >= 140 && hibin < 150)
1170  mPFVsPtInitial_Barrel_Centrality_75_70->Fill(vsPtInitial);
1171 
1172  if(hibin >= 130 && hibin < 140)
1173  mPFVsPtInitial_Barrel_Centrality_70_65->Fill(vsPtInitial);
1174 
1175  if(hibin >= 120 && hibin < 130)
1176  mPFVsPtInitial_Barrel_Centrality_65_60->Fill(vsPtInitial);
1177 
1178  if(hibin >= 110 && hibin < 120)
1179  mPFVsPtInitial_Barrel_Centrality_60_55->Fill(vsPtInitial);
1180 
1181  if(hibin >= 100 && hibin < 110)
1182  mPFVsPtInitial_Barrel_Centrality_55_50->Fill(vsPtInitial);
1183 
1184  if(hibin >= 90 && hibin < 100)
1185  mPFVsPtInitial_Barrel_Centrality_50_45->Fill(vsPtInitial);
1186 
1187  if(hibin >= 80 && hibin < 90)
1188  mPFVsPtInitial_Barrel_Centrality_45_40->Fill(vsPtInitial);
1189 
1190  if(hibin >= 70 && hibin < 80)
1191  mPFVsPtInitial_Barrel_Centrality_40_35->Fill(vsPtInitial);
1192 
1193  if(hibin >= 60 && hibin < 70)
1194  mPFVsPtInitial_Barrel_Centrality_35_30->Fill(vsPtInitial);
1195 
1196  if(hibin >= 50 && hibin < 60)
1197  mPFVsPtInitial_Barrel_Centrality_30_25->Fill(vsPtInitial);
1198 
1199  if(hibin >= 40 && hibin < 50)
1200  mPFVsPtInitial_Barrel_Centrality_25_20->Fill(vsPtInitial);
1201 
1202  if(hibin >= 30 && hibin < 40)
1203  mPFVsPtInitial_Barrel_Centrality_20_15->Fill(vsPtInitial);
1204 
1205  if(hibin >= 20 && hibin < 30)
1206  mPFVsPtInitial_Barrel_Centrality_15_10->Fill(vsPtInitial);
1207 
1208  if(hibin >= 10 && hibin < 20)
1209  mPFVsPtInitial_Barrel_Centrality_10_5->Fill(vsPtInitial);
1210 
1211  if(hibin >= 0 && hibin < 10)
1212  mPFVsPtInitial_Barrel_Centrality_5_0->Fill(vsPtInitial);
1213  }
1214 
1215 
1216  if(isEndcap)
1217  {
1218 
1219  if(hibin >= 190 && hibin < 200)
1220  mPFVsPtInitial_EndCap_Centrality_100_95->Fill(vsPtInitial);
1221 
1222  if(hibin >= 180 && hibin < 190)
1223  mPFVsPtInitial_EndCap_Centrality_95_90->Fill(vsPtInitial);
1224 
1225  if(hibin >= 170 && hibin < 180)
1226  mPFVsPtInitial_EndCap_Centrality_90_85->Fill(vsPtInitial);
1227 
1228  if(hibin >= 160 && hibin < 170)
1229  mPFVsPtInitial_EndCap_Centrality_85_80->Fill(vsPtInitial);
1230 
1231  if(hibin >= 150 && hibin < 160)
1232  mPFVsPtInitial_EndCap_Centrality_80_75->Fill(vsPtInitial);
1233 
1234  if(hibin >= 140 && hibin < 150)
1235  mPFVsPtInitial_EndCap_Centrality_75_70->Fill(vsPtInitial);
1236 
1237  if(hibin >= 130 && hibin < 140)
1238  mPFVsPtInitial_EndCap_Centrality_70_65->Fill(vsPtInitial);
1239 
1240  if(hibin >= 120 && hibin < 130)
1241  mPFVsPtInitial_EndCap_Centrality_65_60->Fill(vsPtInitial);
1242 
1243  if(hibin >= 110 && hibin < 120)
1244  mPFVsPtInitial_EndCap_Centrality_60_55->Fill(vsPtInitial);
1245 
1246  if(hibin >= 100 && hibin < 110)
1247  mPFVsPtInitial_EndCap_Centrality_55_50->Fill(vsPtInitial);
1248 
1249  if(hibin >= 90 && hibin < 100)
1250  mPFVsPtInitial_EndCap_Centrality_50_45->Fill(vsPtInitial);
1251 
1252  if(hibin >= 80 && hibin < 90)
1253  mPFVsPtInitial_EndCap_Centrality_45_40->Fill(vsPtInitial);
1254 
1255  if(hibin >= 70 && hibin < 80)
1256  mPFVsPtInitial_EndCap_Centrality_40_35->Fill(vsPtInitial);
1257 
1258  if(hibin >= 60 && hibin < 70)
1259  mPFVsPtInitial_EndCap_Centrality_35_30->Fill(vsPtInitial);
1260 
1261  if(hibin >= 50 && hibin < 60)
1262  mPFVsPtInitial_EndCap_Centrality_30_25->Fill(vsPtInitial);
1263 
1264  if(hibin >= 40 && hibin < 50)
1265  mPFVsPtInitial_EndCap_Centrality_25_20->Fill(vsPtInitial);
1266 
1267  if(hibin >= 30 && hibin < 40)
1268  mPFVsPtInitial_EndCap_Centrality_20_15->Fill(vsPtInitial);
1269 
1270  if(hibin >= 20 && hibin < 30)
1271  mPFVsPtInitial_EndCap_Centrality_15_10->Fill(vsPtInitial);
1272 
1273  if(hibin >= 10 && hibin < 20)
1274  mPFVsPtInitial_EndCap_Centrality_10_5->Fill(vsPtInitial);
1275 
1276  if(hibin >= 0 && hibin < 10)
1277  mPFVsPtInitial_EndCap_Centrality_5_0->Fill(vsPtInitial);
1278 
1279  }
1280 
1281 
1282  if(isForward)
1283  {
1284 
1285 
1286  if(hibin >= 190 && hibin < 200)
1287  mPFVsPtInitial_HF_Centrality_100_95->Fill(vsPtInitial);
1288 
1289  if(hibin >= 180 && hibin < 190)
1290  mPFVsPtInitial_HF_Centrality_95_90->Fill(vsPtInitial);
1291 
1292  if(hibin >= 170 && hibin < 180)
1293  mPFVsPtInitial_HF_Centrality_90_85->Fill(vsPtInitial);
1294 
1295  if(hibin >= 160 && hibin < 170)
1296  mPFVsPtInitial_HF_Centrality_85_80->Fill(vsPtInitial);
1297 
1298  if(hibin >= 150 && hibin < 160)
1299  mPFVsPtInitial_HF_Centrality_80_75->Fill(vsPtInitial);
1300 
1301  if(hibin >= 140 && hibin < 150)
1302  mPFVsPtInitial_HF_Centrality_75_70->Fill(vsPtInitial);
1303 
1304  if(hibin >= 130 && hibin < 140)
1305  mPFVsPtInitial_HF_Centrality_70_65->Fill(vsPtInitial);
1306 
1307  if(hibin >= 120 && hibin < 130)
1308  mPFVsPtInitial_HF_Centrality_65_60->Fill(vsPtInitial);
1309 
1310  if(hibin >= 110 && hibin < 120)
1311  mPFVsPtInitial_HF_Centrality_60_55->Fill(vsPtInitial);
1312 
1313  if(hibin >= 100 && hibin < 110)
1314  mPFVsPtInitial_HF_Centrality_55_50->Fill(vsPtInitial);
1315 
1316  if(hibin >= 90 && hibin < 100)
1317  mPFVsPtInitial_HF_Centrality_50_45->Fill(vsPtInitial);
1318 
1319  if(hibin >= 80 && hibin < 90)
1320  mPFVsPtInitial_HF_Centrality_45_40->Fill(vsPtInitial);
1321 
1322  if(hibin >= 70 && hibin < 80)
1323  mPFVsPtInitial_HF_Centrality_40_35->Fill(vsPtInitial);
1324 
1325  if(hibin >= 60 && hibin < 70)
1326  mPFVsPtInitial_HF_Centrality_35_30->Fill(vsPtInitial);
1327 
1328  if(hibin >= 50 && hibin < 60)
1329  mPFVsPtInitial_HF_Centrality_30_25->Fill(vsPtInitial);
1330 
1331  if(hibin >= 40 && hibin < 50)
1332  mPFVsPtInitial_HF_Centrality_25_20->Fill(vsPtInitial);
1333 
1334  if(hibin >= 30 && hibin < 40)
1335  mPFVsPtInitial_HF_Centrality_20_15->Fill(vsPtInitial);
1336 
1337  if(hibin >= 20 && hibin < 30)
1338  mPFVsPtInitial_HF_Centrality_15_10->Fill(vsPtInitial);
1339 
1340  if(hibin >= 10 && hibin < 20)
1341  mPFVsPtInitial_HF_Centrality_10_5->Fill(vsPtInitial);
1342 
1343  if(hibin >= 0 && hibin < 10)
1344  mPFVsPtInitial_HF_Centrality_5_0->Fill(vsPtInitial);
1345 
1346 
1347  }
1348 
1349 
1350 
1351 
1352  mPFArea->Fill(vsArea);
1353 
1354 
1355 
1356  }// pf candidate loop
1357 
1358  for(int k = 0;k<nedge_pseudorapidity-1;k++){
1359 
1360  hSumPFVsPtInitial[k]->Fill(SumPFVsPtInitial[k]);
1361  hSumPFVsPt[k]->Fill(SumPFVsPt[k]);
1362  hSumPFPt[k]->Fill(SumPFPt[k]);
1363 
1364  }// eta bin loop
1365 
1366  Float_t Evt_SumPFVsPt = 0;
1367  Float_t Evt_SumPFVsPtInitial = 0;
1368  Float_t Evt_SumPFPt = 0;
1369  Float_t Evt_SumSquaredPFVsPt = 0;
1370  Float_t Evt_SumSquaredPFVsPtInitial = 0;
1371  Float_t Evt_SumSquaredPFPt = 0;
1372 
1373  mSumPFVsPtInitial_n5p191_n2p650->Fill(SumPFVsPtInitial[0]);
1374  mSumPFVsPtInitial_n2p650_n2p043->Fill(SumPFVsPtInitial[1]);
1375  mSumPFVsPtInitial_n2p043_n1p740->Fill(SumPFVsPtInitial[2]);
1376  mSumPFVsPtInitial_n1p740_n1p479->Fill(SumPFVsPtInitial[3]);
1377  mSumPFVsPtInitial_n1p479_n1p131->Fill(SumPFVsPtInitial[4]);
1378  mSumPFVsPtInitial_n1p131_n0p783->Fill(SumPFVsPtInitial[5]);
1379  mSumPFVsPtInitial_n0p783_n0p522->Fill(SumPFVsPtInitial[6]);
1380  mSumPFVsPtInitial_n0p522_0p522->Fill(SumPFVsPtInitial[7]);
1381  mSumPFVsPtInitial_0p522_0p783->Fill(SumPFVsPtInitial[8]);
1382  mSumPFVsPtInitial_0p783_1p131->Fill(SumPFVsPtInitial[9]);
1383  mSumPFVsPtInitial_1p131_1p479->Fill(SumPFVsPtInitial[10]);
1384  mSumPFVsPtInitial_1p479_1p740->Fill(SumPFVsPtInitial[11]);
1385  mSumPFVsPtInitial_1p740_2p043->Fill(SumPFVsPtInitial[12]);
1386  mSumPFVsPtInitial_2p043_2p650->Fill(SumPFVsPtInitial[13]);
1387  mSumPFVsPtInitial_2p650_5p191->Fill(SumPFVsPtInitial[14]);
1388 
1389  mSumPFVsPt_n5p191_n2p650->Fill(SumPFVsPt[0]);
1390  mSumPFVsPt_n2p650_n2p043->Fill(SumPFVsPt[1]);
1391  mSumPFVsPt_n2p043_n1p740->Fill(SumPFVsPt[2]);
1392  mSumPFVsPt_n1p740_n1p479->Fill(SumPFVsPt[3]);
1393  mSumPFVsPt_n1p479_n1p131->Fill(SumPFVsPt[4]);
1394  mSumPFVsPt_n1p131_n0p783->Fill(SumPFVsPt[5]);
1395  mSumPFVsPt_n0p783_n0p522->Fill(SumPFVsPt[6]);
1396  mSumPFVsPt_n0p522_0p522->Fill(SumPFVsPt[7]);
1397  mSumPFVsPt_0p522_0p783->Fill(SumPFVsPt[8]);
1398  mSumPFVsPt_0p783_1p131->Fill(SumPFVsPt[9]);
1399  mSumPFVsPt_1p131_1p479->Fill(SumPFVsPt[10]);
1400  mSumPFVsPt_1p479_1p740->Fill(SumPFVsPt[11]);
1401  mSumPFVsPt_1p740_2p043->Fill(SumPFVsPt[12]);
1402  mSumPFVsPt_2p043_2p650->Fill(SumPFVsPt[13]);
1403  mSumPFVsPt_2p650_5p191->Fill(SumPFVsPt[14]);
1404 
1405  mSumPFPt_n5p191_n2p650->Fill(SumPFPt[0]);
1406  mSumPFPt_n2p650_n2p043->Fill(SumPFPt[1]);
1407  mSumPFPt_n2p043_n1p740->Fill(SumPFPt[2]);
1408  mSumPFPt_n1p740_n1p479->Fill(SumPFPt[3]);
1409  mSumPFPt_n1p479_n1p131->Fill(SumPFPt[4]);
1410  mSumPFPt_n1p131_n0p783->Fill(SumPFPt[5]);
1411  mSumPFPt_n0p783_n0p522->Fill(SumPFPt[6]);
1412  mSumPFPt_n0p522_0p522->Fill(SumPFPt[7]);
1413  mSumPFPt_0p522_0p783->Fill(SumPFPt[8]);
1414  mSumPFPt_0p783_1p131->Fill(SumPFPt[9]);
1415  mSumPFPt_1p131_1p479->Fill(SumPFPt[10]);
1416  mSumPFPt_1p479_1p740->Fill(SumPFPt[11]);
1417  mSumPFPt_1p740_2p043->Fill(SumPFPt[12]);
1418  mSumPFPt_2p043_2p650->Fill(SumPFPt[13]);
1419  mSumPFPt_2p650_5p191->Fill(SumPFPt[14]);
1420 
1421  for(size_t k = 0;k<nedge_pseudorapidity-1;k++){
1422 
1423  Evt_SumPFVsPtInitial = Evt_SumPFVsPtInitial + SumPFVsPtInitial[k];
1424  Evt_SumPFVsPt = Evt_SumPFVsPt + SumPFVsPt[k];
1425  Evt_SumPFPt = Evt_SumPFPt + SumPFPt[k];
1426 
1427  mSumPFVsPtInitial_eta->Fill(edge_pseudorapidity[k],SumPFVsPtInitial[k]);
1428  mSumPFVsPt_eta->Fill(edge_pseudorapidity[k],SumPFVsPt[k]);
1429  mSumPFPt_eta->Fill(edge_pseudorapidity[k],SumPFPt[k]);
1430 
1431  Evt_SumSquaredPFVsPtInitial = Evt_SumSquaredPFVsPtInitial + SumSquaredPFVsPtInitial[k];
1432  Evt_SumSquaredPFVsPt = Evt_SumSquaredPFVsPt + SumSquaredPFVsPt[k];
1433  Evt_SumSquaredPFPt = Evt_SumSquaredPFPt + SumSquaredPFPt[k];
1434 
1435  mSumSquaredPFVsPtInitial_eta->Fill(edge_pseudorapidity[k],hSumPFVsPtInitial[k]->GetRMS(1));
1436  mSumSquaredPFVsPt_eta->Fill(edge_pseudorapidity[k],hSumPFVsPt[k]->GetRMS(1));
1437  mSumSquaredPFPt_eta->Fill(edge_pseudorapidity[k],hSumPFPt[k]->GetRMS(1));
1438 
1439  delete hSumPFVsPtInitial[k];
1440  delete hSumPFVsPt[k];
1441  delete hSumPFPt[k];
1442 
1443  }// eta bin loop
1444 
1445  mSumPFVsPtInitial->Fill(Evt_SumPFVsPtInitial);
1446  mSumPFVsPt->Fill(Evt_SumPFVsPt);
1447  mSumPFPt->Fill(Evt_SumPFPt);
1448  mSumPFVsPtInitial_HF->Fill(Evt_SumPFVsPtInitial,HF_energy);
1449  mSumPFVsPt_HF->Fill(Evt_SumPFVsPt,HF_energy);
1450  mSumPFPt_HF->Fill(Evt_SumPFPt,HF_energy);
1451 
1452  mSumSquaredPFVsPtInitial->Fill(Evt_SumSquaredPFVsPtInitial);
1453  mSumSquaredPFVsPt->Fill(Evt_SumSquaredPFVsPt);
1454  mSumSquaredPFPt->Fill(Evt_SumSquaredPFPt);
1455 
1456  mNPFpart->Fill(NPFpart);
1457  mSumpt->Fill(SumPt_value);
1458 
1459  }
1460 
1461  if (isCaloJet)
1462  {
1463  for (unsigned ijet=0; ijet<caloJets->size(); ijet++)
1464  {
1465  recoJets.push_back((*caloJets)[ijet]);
1466  }
1467  }
1468 
1469  if (isJPTJet)
1470  {
1471  for (unsigned ijet=0; ijet<jptJets->size(); ijet++)
1472  recoJets.push_back((*jptJets)[ijet]);
1473  }
1474 
1475  if (isPFJet)
1476  {
1477  if(std::string("Pu")==UEAlgo)
1478  {
1479  for (unsigned ijet=0; ijet<basicJets->size();ijet++)
1480  {
1481  recoJets.push_back((*basicJets)[ijet]);
1482  }
1483  }
1484  if(std::string("Vs")==UEAlgo){
1485  for (unsigned ijet=0; ijet<pfJets->size(); ijet++)
1486  {
1487  recoJets.push_back((*pfJets)[ijet]);
1488  }
1489  }
1490  }
1491 
1492 
1493  if (isCaloJet && !caloJets.isValid())
1494  {
1495  return;
1496  }
1497  if (isJPTJet && !jptJets.isValid())
1498  {
1499  return;
1500  }
1501  if (isPFJet){
1502  if(std::string("Pu")==UEAlgo){if(!basicJets.isValid()) return;}
1503  if(std::string("Vs")==UEAlgo){if(!pfJets.isValid()) return;}
1504  }
1505 
1506 
1507  int nJet_40 = 0;
1508 
1509  mNJets->Fill(recoJets.size());
1510 
1511  for (unsigned ijet=0; ijet<recoJets.size(); ijet++) {
1512  if (recoJets[ijet].pt() > mRecoJetPtThreshold) {
1513  //counting forward and barrel jets
1514  // get an idea of no of jets with pT>40 GeV
1515  if(recoJets[ijet].pt() > 40)
1516  nJet_40++;
1517  if (mEta) mEta->Fill(recoJets[ijet].eta());
1518  if (mjetpileup) mjetpileup->Fill(recoJets[ijet].pileup());
1519  if (mJetArea) mJetArea ->Fill(recoJets[ijet].jetArea());
1520  if (mPhi) mPhi ->Fill(recoJets[ijet].phi());
1521  if (mEnergy) mEnergy ->Fill(recoJets[ijet].energy());
1522  if (mP) mP ->Fill(recoJets[ijet].p());
1523  if (mPt) mPt ->Fill(recoJets[ijet].pt());
1524  if (mMass) mMass ->Fill(recoJets[ijet].mass());
1525  if (mConstituents) mConstituents->Fill(recoJets[ijet].nConstituents());
1526 
1527  for(size_t iii = 0 ; iii < numbers.size() ; iii++)
1528  {
1529  pfDeltaR = sqrt((numbers[iii][2]-recoJets[ijet].phi())*(numbers[iii][2]-recoJets[ijet].phi()) + (numbers[iii][1]-recoJets[ijet].eta())*(numbers[iii][1]-recoJets[ijet].eta())); //MZ
1530 
1531  mPFVsPtInitialDeltaR_pTCorrected->Fill(pfDeltaR,numbers[iii][3]/recoJets[ijet].pt()); //MZ
1532 
1533 
1534  mPFVsPtDeltaR_pTCorrected->Fill(pfDeltaR,numbers[iii][4]/recoJets[ijet].pt()); //MZ
1535 
1536  //if(pfDeltaR < 1.0)
1537  // mPFDeltaR ->Fill(pfDeltaR,numbers[iii][0]/recoJets[ijet].pt()); //MZ
1538  mPFDeltaR ->Fill(pfDeltaR); //MZ
1539  mPFDeltaR_Scaled_R->Fill(pfDeltaR,1. / pow(pfDeltaR,2)); //MZ
1540  //mPFDeltaR_pTCorrected->Fill(pfDeltaR,numbers[iii][0]/recoJets[ijet].pt()); //MZ
1541  mPFDeltaR_pTCorrected->Fill(pfDeltaR,numbers[iii][0]/recoJets[ijet].pt()); //MZ
1542 
1543  //mPFVsPtInitialDeltaR_pTCorrected->Fill(pfDeltaR,numbers[iii][3]/recoJets[ijet].pt()); //MZ
1544 
1545  if(recoJets[ijet].pt() > 20 && recoJets[ijet].pt() < 30)
1546  {
1547  mPFDeltaR_pTCorrected_PFpT_20To30->Fill(pfDeltaR,numbers[iii][0]/recoJets[ijet].pt()); //MZ
1548  mPFDeltaR_pTCorrected_PFVsInitialpT_20To30->Fill(pfDeltaR,numbers[iii][3]/recoJets[ijet].pt()); //MZ
1549  }
1550 
1551  if(recoJets[ijet].pt() > 30 && recoJets[ijet].pt() < 50)
1552  {
1553  mPFDeltaR_pTCorrected_PFpT_30To50->Fill(pfDeltaR,numbers[iii][0]/recoJets[ijet].pt()); //MZ
1554  mPFDeltaR_pTCorrected_PFVsInitialpT_30To50->Fill(pfDeltaR,numbers[iii][3]/recoJets[ijet].pt()); //MZ
1555  }
1556 
1557  if(recoJets[ijet].pt() > 50 && recoJets[ijet].pt() < 80)
1558  {
1559  mPFDeltaR_pTCorrected_PFpT_50To80->Fill(pfDeltaR,numbers[iii][0]/recoJets[ijet].pt()); //MZ
1560  mPFDeltaR_pTCorrected_PFVsInitialpT_50To80->Fill(pfDeltaR,numbers[iii][3]/recoJets[ijet].pt()); //MZ
1561  }
1562 
1563 
1564  if(recoJets[ijet].pt() > 80 && recoJets[ijet].pt() < 120)
1565  {
1566  mPFDeltaR_pTCorrected_PFpT_80To120->Fill(pfDeltaR,numbers[iii][0]/recoJets[ijet].pt()); //MZ
1567  mPFDeltaR_pTCorrected_PFVsInitialpT_80To120->Fill(pfDeltaR,numbers[iii][3]/recoJets[ijet].pt()); //MZ
1568  }
1569 
1570 
1571  if(recoJets[ijet].pt() > 120 && recoJets[ijet].pt() < 180)
1572  {
1573  mPFDeltaR_pTCorrected_PFpT_120To180->Fill(pfDeltaR,numbers[iii][0]/recoJets[ijet].pt()); //MZ
1574  mPFDeltaR_pTCorrected_PFVsInitialpT_120To180->Fill(pfDeltaR,numbers[iii][3]/recoJets[ijet].pt()); //MZ
1575  }
1576 
1577 
1578  if(recoJets[ijet].pt() > 180 && recoJets[ijet].pt() < 300)
1579  {
1580  mPFDeltaR_pTCorrected_PFpT_180To300->Fill(pfDeltaR,numbers[iii][0]/recoJets[ijet].pt()); //MZ
1581  mPFDeltaR_pTCorrected_PFVsInitialpT_180To300->Fill(pfDeltaR,numbers[iii][3]/recoJets[ijet].pt()); //MZ
1582  }
1583 
1584  if(recoJets[ijet].pt() > 300)
1585  {
1586  mPFDeltaR_pTCorrected_PFpT_300ToInf->Fill(pfDeltaR,numbers[iii][0]/recoJets[ijet].pt()); //MZ
1587  mPFDeltaR_pTCorrected_PFVsInitialpT_300ToInf->Fill(pfDeltaR,numbers[iii][3]/recoJets[ijet].pt()); //MZ
1588  }
1589 
1590  }
1591 
1592  }
1593  //cout << "End RecoJets loop" << endl;
1594  if (mNJets_40)
1595  mNJets_40->Fill(nJet_40);
1596 
1597  } // end recojet loop 222
1598 
1599  numbers.clear();
1600 
1601 }
JetAnalyzer_HeavyIons::etaBins_
static const Int_t etaBins_
Definition: JetAnalyzer_HeavyIons.h:429
JetAnalyzer_HeavyIons::mPFCandpT_Forward_photon
MonitorElement * mPFCandpT_Forward_photon
Definition: JetAnalyzer_HeavyIons.h:417
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFVsInitialpT_30To50
MonitorElement * mPFDeltaR_pTCorrected_PFVsInitialpT_30To50
Definition: JetAnalyzer_HeavyIons.h:214
JetAnalyzer_HeavyIons::caloJetsToken_
edm::EDGetTokenT< reco::CaloJetCollection > caloJetsToken_
Definition: JetAnalyzer_HeavyIons.h:118
JetAnalyzer_HeavyIons::pfCandToken_
edm::EDGetTokenT< reco::PFCandidateCollection > pfCandToken_
Definition: JetAnalyzer_HeavyIons.h:122
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_25_20
MonitorElement * mPFVsPtInitial_Barrel_Centrality_25_20
Definition: JetAnalyzer_HeavyIons.h:153
JetAnalyzer_HeavyIons::mCaloEta
MonitorElement * mCaloEta
Definition: JetAnalyzer_HeavyIons.h:225
JetAnalyzer_HeavyIons::mSumPFVsPt_2p650_5p191
MonitorElement * mSumPFVsPt_2p650_5p191
Definition: JetAnalyzer_HeavyIons.h:306
i
int i
Definition: DBlmapReader.cc:9
JetAnalyzer_HeavyIons::mSumCaloVsPt_HF
MonitorElement * mSumCaloVsPt_HF
Definition: JetAnalyzer_HeavyIons.h:272
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_n2p043_n1p740
MonitorElement * mSumCaloVsPtInitial_n2p043_n1p740
Definition: JetAnalyzer_HeavyIons.h:327
JetAnalyzer_HeavyIons::mPFEta
MonitorElement * mPFEta
Definition: JetAnalyzer_HeavyIons.h:132
JetAnalyzer_HeavyIons::mSumCaloVsPt_n5p191_n2p650
MonitorElement * mSumCaloVsPt_n5p191_n2p650
Definition: JetAnalyzer_HeavyIons.h:341
JetAnalyzer_HeavyIons::mSumCaloVsPt_n1p131_n0p783
MonitorElement * mSumCaloVsPt_n1p131_n0p783
Definition: JetAnalyzer_HeavyIons.h:346
JetAnalyzer_HeavyIons::mSumSquaredPFPt
MonitorElement * mSumSquaredPFPt
Definition: JetAnalyzer_HeavyIons.h:249
JetAnalyzer_HeavyIons::mPFCandpT_vs_eta_electron
MonitorElement * mPFCandpT_vs_eta_electron
Definition: JetAnalyzer_HeavyIons.h:388
JetAnalyzer_HeavyIons::mSumPFVsPtInitial
MonitorElement * mSumPFVsPtInitial
Definition: JetAnalyzer_HeavyIons.h:232
JetAnalyzer_HeavyIons::caloCandViewToken_
edm::EDGetTokenT< reco::CandidateView > caloCandViewToken_
Definition: JetAnalyzer_HeavyIons.h:124
JetAnalyzer_HeavyIons::mSumCaloPt_n2p043_n1p740
MonitorElement * mSumCaloPt_n2p043_n1p740
Definition: JetAnalyzer_HeavyIons.h:359
JetAnalyzer_HeavyIons::mCaloPhi
MonitorElement * mCaloPhi
Definition: JetAnalyzer_HeavyIons.h:226
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_10_5
MonitorElement * mPFVsPtInitial_Barrel_Centrality_10_5
Definition: JetAnalyzer_HeavyIons.h:156
JetAnalyzer_HeavyIons::mPFCandpT_Forward_EME_inHF
MonitorElement * mPFCandpT_Forward_EME_inHF
Definition: JetAnalyzer_HeavyIons.h:420
JetAnalyzer_HeavyIons::Background
edm::InputTag Background
Definition: JetAnalyzer_HeavyIons.h:101
EndcapEta
const Double_t EndcapEta
Definition: JetAnalyzer_HeavyIons.h:78
JetAnalyzer_HeavyIons::jptJetsToken_
edm::EDGetTokenT< reco::JPTJetCollection > jptJetsToken_
Definition: JetAnalyzer_HeavyIons.h:121
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_HF
MonitorElement * mSumCaloVsPtInitial_HF
Definition: JetAnalyzer_HeavyIons.h:273
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_1p131_1p479
MonitorElement * mSumPFVsPtInitial_1p131_1p479
Definition: JetAnalyzer_HeavyIons.h:286
JetAnalyzer_HeavyIons::mSumSquaredCaloVsPt_eta
MonitorElement * mSumSquaredCaloVsPt_eta
Definition: JetAnalyzer_HeavyIons.h:260
JetAnalyzer_HeavyIons::mSumPFVsPt_n5p191_n2p650
MonitorElement * mSumPFVsPt_n5p191_n2p650
Definition: JetAnalyzer_HeavyIons.h:292
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_60_55
MonitorElement * mPFVsPtInitial_EndCap_Centrality_60_55
Definition: JetAnalyzer_HeavyIons.h:167
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFVsInitialpT_80To120
MonitorElement * mPFDeltaR_pTCorrected_PFVsInitialpT_80To120
Definition: JetAnalyzer_HeavyIons.h:216
JetAnalyzer_HeavyIons::centralityTag_
edm::InputTag centralityTag_
Definition: JetAnalyzer_HeavyIons.h:107
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_n2p650_n2p043
MonitorElement * mSumCaloVsPtInitial_n2p650_n2p043
Definition: JetAnalyzer_HeavyIons.h:326
JetAnalyzer_HeavyIons::mCaloVsPtInitial
MonitorElement * mCaloVsPtInitial
Definition: JetAnalyzer_HeavyIons.h:228
CaloTower::p4
math::PtEtaPhiMLorentzVector p4(double vtxZ) const
Definition: CaloTower.cc:129
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFpT_180To300
MonitorElement * mPFDeltaR_pTCorrected_PFpT_180To300
Definition: JetAnalyzer_HeavyIons.h:210
JetAnalyzer_HeavyIons::mSumCaloVsPt_1p479_1p740
MonitorElement * mSumCaloVsPt_1p479_1p740
Definition: JetAnalyzer_HeavyIons.h:352
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_55_50
MonitorElement * mPFVsPtInitial_Barrel_Centrality_55_50
Definition: JetAnalyzer_HeavyIons.h:147
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_90_85
MonitorElement * mPFVsPtInitial_Barrel_Centrality_90_85
Definition: JetAnalyzer_HeavyIons.h:140
JetAnalyzer_HeavyIons::mSumCaloVsPt_n2p650_n2p043
MonitorElement * mSumCaloVsPt_n2p650_n2p043
Definition: JetAnalyzer_HeavyIons.h:342
JetAnalyzer_HeavyIons::mSumCaloVsPt_n0p783_n0p522
MonitorElement * mSumCaloVsPt_n0p783_n0p522
Definition: JetAnalyzer_HeavyIons.h:347
JetAnalyzer_HeavyIons::mSumSquaredCaloVsPtInitial_eta
MonitorElement * mSumSquaredCaloVsPtInitial_eta
Definition: JetAnalyzer_HeavyIons.h:259
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_15_10
MonitorElement * mPFVsPtInitial_EndCap_Centrality_15_10
Definition: JetAnalyzer_HeavyIons.h:176
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_1p479_1p740
MonitorElement * mSumCaloVsPtInitial_1p479_1p740
Definition: JetAnalyzer_HeavyIons.h:336
JetAnalyzer_HeavyIons::mPFCandpT_Forward_muon
MonitorElement * mPFCandpT_Forward_muon
Definition: JetAnalyzer_HeavyIons.h:416
JetAnalyzer_HeavyIons.h
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_75_70
MonitorElement * mPFVsPtInitial_HF_Centrality_75_70
Definition: JetAnalyzer_HeavyIons.h:185
JetAnalyzer_HeavyIons::mSumPFVsPt_n2p043_n1p740
MonitorElement * mSumPFVsPt_n2p043_n1p740
Definition: JetAnalyzer_HeavyIons.h:294
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_55_50
MonitorElement * mPFVsPtInitial_EndCap_Centrality_55_50
Definition: JetAnalyzer_HeavyIons.h:168
JetAnalyzer_HeavyIons::isCaloJet
bool isCaloJet
Definition: JetAnalyzer_HeavyIons.h:424
JetAnalyzer_HeavyIons::pfJetsToken_
edm::EDGetTokenT< reco::PFJetCollection > pfJetsToken_
Definition: JetAnalyzer_HeavyIons.h:119
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_1p131_1p479
MonitorElement * mSumCaloVsPtInitial_1p131_1p479
Definition: JetAnalyzer_HeavyIons.h:335
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_20_15
MonitorElement * mPFVsPtInitial_HF_Centrality_20_15
Definition: JetAnalyzer_HeavyIons.h:196
JetAnalyzer_HeavyIons::mSumCaloPt_0p783_1p131
MonitorElement * mSumCaloPt_0p783_1p131
Definition: JetAnalyzer_HeavyIons.h:366
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_25_20
MonitorElement * mPFVsPtInitial_EndCap_Centrality_25_20
Definition: JetAnalyzer_HeavyIons.h:174
GeomDetEnumerators::isBarrel
bool isBarrel(GeomDetEnumerators::SubDetector m)
Definition: GeomDetEnumerators.cc:32
JetAnalyzer_HeavyIons::mPFCandpT_Endcap_NeutralHadron
MonitorElement * mPFCandpT_Endcap_NeutralHadron
Definition: JetAnalyzer_HeavyIons.h:409
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF
MonitorElement * mPFVsPtInitial_HF
Definition: JetAnalyzer_HeavyIons.h:136
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_25_20
MonitorElement * mPFVsPtInitial_HF_Centrality_25_20
Definition: JetAnalyzer_HeavyIons.h:195
JetAnalyzer_HeavyIons::mCaloArea
MonitorElement * mCaloArea
Definition: JetAnalyzer_HeavyIons.h:229
JetAnalyzer_HeavyIons::backgrounds_value_
edm::EDGetTokenT< std::vector< float > > backgrounds_value_
Definition: JetAnalyzer_HeavyIons.h:127
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_n1p479_n1p131
MonitorElement * mSumCaloVsPtInitial_n1p479_n1p131
Definition: JetAnalyzer_HeavyIons.h:329
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_5_0
MonitorElement * mPFVsPtInitial_Barrel_Centrality_5_0
Definition: JetAnalyzer_HeavyIons.h:157
JetAnalyzer_HeavyIons::mPFCandpT_vs_eta_EME_inHF
MonitorElement * mPFCandpT_vs_eta_EME_inHF
Definition: JetAnalyzer_HeavyIons.h:393
JetAnalyzer_HeavyIons::mSumPFPt_n0p522_0p522
MonitorElement * mSumPFPt_n0p522_0p522
Definition: JetAnalyzer_HeavyIons.h:315
JetAnalyzer_HeavyIons::mSumSquaredCaloPt
MonitorElement * mSumSquaredCaloPt
Definition: JetAnalyzer_HeavyIons.h:257
ecaldqm::isForward
bool isForward(DetId const &)
Definition: EcalDQMCommonUtils.cc:276
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_95_90
MonitorElement * mPFVsPtInitial_Barrel_Centrality_95_90
Definition: JetAnalyzer_HeavyIons.h:139
JetAnalyzer_HeavyIons::mPhi
MonitorElement * mPhi
Definition: JetAnalyzer_HeavyIons.h:375
JetAnalyzer_HeavyIons::mSumPFPt_2p043_2p650
MonitorElement * mSumPFPt_2p043_2p650
Definition: JetAnalyzer_HeavyIons.h:321
JetAnalyzer_HeavyIons::mSumPFVsPt_0p522_0p783
MonitorElement * mSumPFVsPt_0p522_0p783
Definition: JetAnalyzer_HeavyIons.h:300
JetAnalyzer_HeavyIons::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: JetAnalyzer_HeavyIons.cc:710
JetAnalyzer_HeavyIons::mSumSquaredCaloVsPtInitial
MonitorElement * mSumSquaredCaloVsPtInitial
Definition: JetAnalyzer_HeavyIons.h:256
JetAnalyzer_HeavyIons::mSumCaloPt_1p131_1p479
MonitorElement * mSumCaloPt_1p131_1p479
Definition: JetAnalyzer_HeavyIons.h:367
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_30_25
MonitorElement * mPFVsPtInitial_EndCap_Centrality_30_25
Definition: JetAnalyzer_HeavyIons.h:173
JetAnalyzer_HeavyIons::mSumSquaredPFVsPt_eta
MonitorElement * mSumSquaredPFVsPt_eta
Definition: JetAnalyzer_HeavyIons.h:252
JetAnalyzer_HeavyIons::mSumPFVsPt_eta
MonitorElement * mSumPFVsPt_eta
Definition: JetAnalyzer_HeavyIons.h:236
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_2p650_5p191
MonitorElement * mSumPFVsPtInitial_2p650_5p191
Definition: JetAnalyzer_HeavyIons.h:290
JetAnalyzer_HeavyIons::mSumCaloVsPt_1p131_1p479
MonitorElement * mSumCaloVsPt_1p131_1p479
Definition: JetAnalyzer_HeavyIons.h:351
JetAnalyzer_HeavyIons::mSumPFPt_n1p740_n1p479
MonitorElement * mSumPFPt_n1p740_n1p479
Definition: JetAnalyzer_HeavyIons.h:311
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_85_80
MonitorElement * mPFVsPtInitial_HF_Centrality_85_80
Definition: JetAnalyzer_HeavyIons.h:183
JetAnalyzer_HeavyIons::mEta
MonitorElement * mEta
Definition: JetAnalyzer_HeavyIons.h:374
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_n2p650_n2p043
MonitorElement * mSumPFVsPtInitial_n2p650_n2p043
Definition: JetAnalyzer_HeavyIons.h:277
JetAnalyzer_HeavyIons::mSumPFPt_2p650_5p191
MonitorElement * mSumPFPt_2p650_5p191
Definition: JetAnalyzer_HeavyIons.h:322
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_0p783_1p131
MonitorElement * mSumPFVsPtInitial_0p783_1p131
Definition: JetAnalyzer_HeavyIons.h:285
JetAnalyzer_HeavyIons::mPFCandpT_Forward_electron
MonitorElement * mPFCandpT_Forward_electron
Definition: JetAnalyzer_HeavyIons.h:415
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:106
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFpT_300ToInf
MonitorElement * mPFDeltaR_pTCorrected_PFpT_300ToInf
Definition: JetAnalyzer_HeavyIons.h:211
JetAnalyzer_HeavyIons::pvToken_
edm::EDGetTokenT< std::vector< reco::Vertex > > pvToken_
Definition: JetAnalyzer_HeavyIons.h:116
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_n0p783_n0p522
MonitorElement * mSumPFVsPtInitial_n0p783_n0p522
Definition: JetAnalyzer_HeavyIons.h:282
HLT_25ns14e33_v1_cff.InputTag
tuple InputTag
Definition: HLT_25ns14e33_v1_cff.py:45333
reco::LeafCandidate::eta
virtual double eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:137
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.h:131
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_10_5
MonitorElement * mPFVsPtInitial_EndCap_Centrality_10_5
Definition: JetAnalyzer_HeavyIons.h:177
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
JetAnalyzer_HeavyIons::mSumCaloPt_eta
MonitorElement * mSumCaloPt_eta
Definition: JetAnalyzer_HeavyIons.h:245
JetAnalyzer_HeavyIons::mSumPFPt_0p522_0p783
MonitorElement * mSumPFPt_0p522_0p783
Definition: JetAnalyzer_HeavyIons.h:316
JetAnalyzer_HeavyIons::mPFArea
MonitorElement * mPFArea
Definition: JetAnalyzer_HeavyIons.h:201
JetAnalyzer_HeavyIons::mSumPFPt_n2p650_n2p043
MonitorElement * mSumPFPt_n2p650_n2p043
Definition: JetAnalyzer_HeavyIons.h:309
JetAnalyzer_HeavyIons::mPFCandpT_vs_eta_muon
MonitorElement * mPFCandpT_vs_eta_muon
Definition: JetAnalyzer_HeavyIons.h:389
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_80_75
MonitorElement * mPFVsPtInitial_EndCap_Centrality_80_75
Definition: JetAnalyzer_HeavyIons.h:163
JetAnalyzer_HeavyIons::mSumCaloPt_n1p740_n1p479
MonitorElement * mSumCaloPt_n1p740_n1p479
Definition: JetAnalyzer_HeavyIons.h:360
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_15_10
MonitorElement * mPFVsPtInitial_Barrel_Centrality_15_10
Definition: JetAnalyzer_HeavyIons.h:155
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFVsInitialpT_50To80
MonitorElement * mPFDeltaR_pTCorrected_PFVsInitialpT_50To80
Definition: JetAnalyzer_HeavyIons.h:215
reco::tau::pfCandidates
std::vector< PFCandidatePtr > pfCandidates(const PFJet &jet, int particleId, bool sort=true)
Definition: RecoTauCommonUtilities.cc:34
JetAnalyzer_HeavyIons::mPt
MonitorElement * mPt
Definition: JetAnalyzer_HeavyIons.h:376
JetAnalyzer_HeavyIons::mSumPFVsPt_1p740_2p043
MonitorElement * mSumPFVsPt_1p740_2p043
Definition: JetAnalyzer_HeavyIons.h:304
JetAnalyzer_HeavyIons::mPFCandpT_Endcap_HadE_inHF
MonitorElement * mPFCandpT_Endcap_HadE_inHF
Definition: JetAnalyzer_HeavyIons.h:410
JetAnalyzer_HeavyIons::centralityBinToken
edm::EDGetTokenT< int > centralityBinToken
Definition: JetAnalyzer_HeavyIons.h:112
JetAnalyzer_HeavyIons::mSumCaloPt_0p522_0p783
MonitorElement * mSumCaloPt_0p522_0p783
Definition: JetAnalyzer_HeavyIons.h:365
JetAnalyzer_HeavyIons::mSumPFPt_HF
MonitorElement * mSumPFPt_HF
Definition: JetAnalyzer_HeavyIons.h:270
JetAnalyzer_HeavyIons::mSumCaloVsPt_n1p479_n1p131
MonitorElement * mSumCaloVsPt_n1p479_n1p131
Definition: JetAnalyzer_HeavyIons.h:345
JetAnalyzer_HeavyIons::mPFCandpT_Endcap_ChargedHadron
MonitorElement * mPFCandpT_Endcap_ChargedHadron
Definition: JetAnalyzer_HeavyIons.h:405
JetAnalyzer_HeavyIons::mSumCaloPt_n2p650_n2p043
MonitorElement * mSumCaloPt_n2p650_n2p043
Definition: JetAnalyzer_HeavyIons.h:358
JetAnalyzer_HeavyIons::mSumPFPt_0p783_1p131
MonitorElement * mSumPFPt_0p783_1p131
Definition: JetAnalyzer_HeavyIons.h:317
JetAnalyzer_HeavyIons::mSumCaloPt_1p479_1p740
MonitorElement * mSumCaloPt_1p479_1p740
Definition: JetAnalyzer_HeavyIons.h:368
estimatePileup_makeJSON.pileup
list pileup
Definition: estimatePileup_makeJSON.py:193
JetAnalyzer_HeavyIons::mSumPFVsPt_1p479_1p740
MonitorElement * mSumPFVsPt_1p479_1p740
Definition: JetAnalyzer_HeavyIons.h:303
JetAnalyzer_HeavyIons::mPFCandpT_vs_eta_NeutralHadron
MonitorElement * mPFCandpT_vs_eta_NeutralHadron
Definition: JetAnalyzer_HeavyIons.h:391
JetAnalyzer_HeavyIons::mSumPFVsPt_HF
MonitorElement * mSumPFVsPt_HF
Definition: JetAnalyzer_HeavyIons.h:268
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_45_40
MonitorElement * mPFVsPtInitial_HF_Centrality_45_40
Definition: JetAnalyzer_HeavyIons.h:191
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_0p522_0p783
MonitorElement * mSumCaloVsPtInitial_0p522_0p783
Definition: JetAnalyzer_HeavyIons.h:333
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_100_95
MonitorElement * mPFVsPtInitial_EndCap_Centrality_100_95
Definition: JetAnalyzer_HeavyIons.h:159
JetAnalyzer_HeavyIons::mSumCaloPt_HF
MonitorElement * mSumCaloPt_HF
Definition: JetAnalyzer_HeavyIons.h:274
JetAnalyzer_HeavyIons::mPFVsPtInitial
MonitorElement * mPFVsPtInitial
Definition: JetAnalyzer_HeavyIons.h:135
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFVsInitialpT_20To30
MonitorElement * mPFDeltaR_pTCorrected_PFVsInitialpT_20To30
Definition: JetAnalyzer_HeavyIons.h:213
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_70_65
MonitorElement * mPFVsPtInitial_EndCap_Centrality_70_65
Definition: JetAnalyzer_HeavyIons.h:165
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_n0p522_0p522
MonitorElement * mSumCaloVsPtInitial_n0p522_0p522
Definition: JetAnalyzer_HeavyIons.h:332
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_n1p479_n1p131
MonitorElement * mSumPFVsPtInitial_n1p479_n1p131
Definition: JetAnalyzer_HeavyIons.h:280
JetAnalyzer_HeavyIons::mSumpt
MonitorElement * mSumpt
Definition: JetAnalyzer_HeavyIons.h:230
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_100_95
MonitorElement * mPFVsPtInitial_Barrel_Centrality_100_95
Definition: JetAnalyzer_HeavyIons.h:138
JetAnalyzer_HeavyIons::mPFCandpT_vs_eta_photon
MonitorElement * mPFCandpT_vs_eta_photon
Definition: JetAnalyzer_HeavyIons.h:390
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_20_15
MonitorElement * mPFVsPtInitial_Barrel_Centrality_20_15
Definition: JetAnalyzer_HeavyIons.h:154
JetAnalyzer_HeavyIons::mSumPFVsPt_n2p650_n2p043
MonitorElement * mSumPFVsPt_n2p650_n2p043
Definition: JetAnalyzer_HeavyIons.h:293
JetAnalyzer_HeavyIons::mSumPFVsPt_n0p522_0p522
MonitorElement * mSumPFVsPt_n0p522_0p522
Definition: JetAnalyzer_HeavyIons.h:299
JetAnalyzer_HeavyIons::centralityToken
edm::EDGetTokenT< reco::Centrality > centralityToken
Definition: JetAnalyzer_HeavyIons.h:108
JetAnalyzer_HeavyIons::mSumSquaredPFVsPtInitial
MonitorElement * mSumSquaredPFVsPtInitial
Definition: JetAnalyzer_HeavyIons.h:248
JetAnalyzer_HeavyIons::mPFCandpT_Forward_ChargedHadron
MonitorElement * mPFCandpT_Forward_ChargedHadron
Definition: JetAnalyzer_HeavyIons.h:414
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_n1p131_n0p783
MonitorElement * mSumCaloVsPtInitial_n1p131_n0p783
Definition: JetAnalyzer_HeavyIons.h:330
JetAnalyzer_HeavyIons::mSumSquaredCaloVsPt
MonitorElement * mSumSquaredCaloVsPt
Definition: JetAnalyzer_HeavyIons.h:255
JetAnalyzer_HeavyIons::mSumPFVsPt_1p131_1p479
MonitorElement * mSumPFVsPt_1p131_1p479
Definition: JetAnalyzer_HeavyIons.h:302
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_75_70
MonitorElement * mPFVsPtInitial_Barrel_Centrality_75_70
Definition: JetAnalyzer_HeavyIons.h:143
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
JetAnalyzer_HeavyIons::mSumPFPt_1p131_1p479
MonitorElement * mSumPFPt_1p131_1p479
Definition: JetAnalyzer_HeavyIons.h:318
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_1p479_1p740
MonitorElement * mSumPFVsPtInitial_1p479_1p740
Definition: JetAnalyzer_HeavyIons.h:287
JetAnalyzer_HeavyIons::mOutputFile
std::string mOutputFile
Definition: JetAnalyzer_HeavyIons.h:98
JetAnalyzer_HeavyIons::mSumPFVsPt_n0p783_n0p522
MonitorElement * mSumPFVsPt_n0p783_n0p522
Definition: JetAnalyzer_HeavyIons.h:298
JetAnalyzer_HeavyIons::mPFCandpT_Endcap_muon
MonitorElement * mPFCandpT_Endcap_muon
Definition: JetAnalyzer_HeavyIons.h:407
JetAnalyzer_HeavyIons::mPFCandpT_Endcap_photon
MonitorElement * mPFCandpT_Endcap_photon
Definition: JetAnalyzer_HeavyIons.h:408
JetAnalyzer_HeavyIons::mPFDeltaR_Scaled_R
MonitorElement * mPFDeltaR_Scaled_R
Definition: JetAnalyzer_HeavyIons.h:203
JetAnalyzer_HeavyIons::mSumPFPt
MonitorElement * mSumPFPt
Definition: JetAnalyzer_HeavyIons.h:233
JetAnalyzer_HeavyIons::mSumPFPt_n1p479_n1p131
MonitorElement * mSumPFPt_n1p479_n1p131
Definition: JetAnalyzer_HeavyIons.h:312
JetAnalyzer_HeavyIons::hiVertexToken_
edm::EDGetTokenT< std::vector< reco::Vertex > > hiVertexToken_
Definition: JetAnalyzer_HeavyIons.h:128
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_65_60
MonitorElement * mPFVsPtInitial_EndCap_Centrality_65_60
Definition: JetAnalyzer_HeavyIons.h:166
JetAnalyzer_HeavyIons::mConstituents
MonitorElement * mConstituents
Definition: JetAnalyzer_HeavyIons.h:380
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFpT_80To120
MonitorElement * mPFDeltaR_pTCorrected_PFpT_80To120
Definition: JetAnalyzer_HeavyIons.h:208
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_50_45
MonitorElement * mPFVsPtInitial_EndCap_Centrality_50_45
Definition: JetAnalyzer_HeavyIons.h:169
JetAnalyzer_HeavyIons::isJPTJet
bool isJPTJet
Definition: JetAnalyzer_HeavyIons.h:425
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_35_30
MonitorElement * mPFVsPtInitial_EndCap_Centrality_35_30
Definition: JetAnalyzer_HeavyIons.h:172
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFpT_50To80
MonitorElement * mPFDeltaR_pTCorrected_PFpT_50To80
Definition: JetAnalyzer_HeavyIons.h:207
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_95_90
MonitorElement * mPFVsPtInitial_EndCap_Centrality_95_90
Definition: JetAnalyzer_HeavyIons.h:160
JetAnalyzer_HeavyIons::mPFCandpT_Barrel_photon
MonitorElement * mPFCandpT_Barrel_photon
Definition: JetAnalyzer_HeavyIons.h:399
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_40_35
MonitorElement * mPFVsPtInitial_HF_Centrality_40_35
Definition: JetAnalyzer_HeavyIons.h:192
JetAnalyzer_HeavyIons::mSumPFVsPt_n1p479_n1p131
MonitorElement * mSumPFVsPt_n1p479_n1p131
Definition: JetAnalyzer_HeavyIons.h:296
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_75_70
MonitorElement * mPFVsPtInitial_EndCap_Centrality_75_70
Definition: JetAnalyzer_HeavyIons.h:164
JetAnalyzer_HeavyIons::centralityBinTag_
edm::InputTag centralityBinTag_
Definition: JetAnalyzer_HeavyIons.h:111
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_50_45
MonitorElement * mPFVsPtInitial_HF_Centrality_50_45
Definition: JetAnalyzer_HeavyIons.h:190
JetAnalyzer_HeavyIons::mP
MonitorElement * mP
Definition: JetAnalyzer_HeavyIons.h:377
JetAnalyzer_HeavyIons::mCaloPt
MonitorElement * mCaloPt
Definition: JetAnalyzer_HeavyIons.h:224
JetAnalyzer_HeavyIons::mNvtx
MonitorElement * mNvtx
Definition: JetAnalyzer_HeavyIons.h:264
JetAnalyzer_HeavyIons::mSumSquaredPFVsPt
MonitorElement * mSumSquaredPFVsPt
Definition: JetAnalyzer_HeavyIons.h:247
JetAnalyzer_HeavyIons::mPFVsPtDeltaR_pTCorrected
MonitorElement * mPFVsPtDeltaR_pTCorrected
Definition: JetAnalyzer_HeavyIons.h:222
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected
MonitorElement * mPFDeltaR_pTCorrected
Definition: JetAnalyzer_HeavyIons.h:204
reco::Vertex::Point
math::XYZPoint Point
point in the space
Definition: Vertex.h:39
JetAnalyzer_HeavyIons::mSumCaloVsPt_n1p740_n1p479
MonitorElement * mSumCaloVsPt_n1p740_n1p479
Definition: JetAnalyzer_HeavyIons.h:344
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_1p740_2p043
MonitorElement * mSumCaloVsPtInitial_1p740_2p043
Definition: JetAnalyzer_HeavyIons.h:337
JetAnalyzer_HeavyIons::mSumCaloPt_n0p783_n0p522
MonitorElement * mSumCaloPt_n0p783_n0p522
Definition: JetAnalyzer_HeavyIons.h:363
JetAnalyzer_HeavyIons::mInputCollection
edm::InputTag mInputCollection
Definition: JetAnalyzer_HeavyIons.h:95
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_90_85
MonitorElement * mPFVsPtInitial_HF_Centrality_90_85
Definition: JetAnalyzer_HeavyIons.h:182
JetAnalyzer_HeavyIons::mSumSquaredCaloPt_eta
MonitorElement * mSumSquaredCaloPt_eta
Definition: JetAnalyzer_HeavyIons.h:261
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFpT_20To30
MonitorElement * mPFDeltaR_pTCorrected_PFpT_20To30
Definition: JetAnalyzer_HeavyIons.h:205
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_2p043_2p650
MonitorElement * mSumCaloVsPtInitial_2p043_2p650
Definition: JetAnalyzer_HeavyIons.h:338
JetAnalyzer_HeavyIons::mPFCandpT_Barrel_Unknown
MonitorElement * mPFCandpT_Barrel_Unknown
Definition: JetAnalyzer_HeavyIons.h:395
JetAnalyzer_HeavyIons::mMass
MonitorElement * mMass
Definition: JetAnalyzer_HeavyIons.h:379
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
JetAnalyzer_HeavyIons::mSumCaloPt_n0p522_0p522
MonitorElement * mSumCaloPt_n0p522_0p522
Definition: JetAnalyzer_HeavyIons.h:364
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFVsInitialpT_120To180
MonitorElement * mPFDeltaR_pTCorrected_PFVsInitialpT_120To180
Definition: JetAnalyzer_HeavyIons.h:217
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_n2p043_n1p740
MonitorElement * mSumPFVsPtInitial_n2p043_n1p740
Definition: JetAnalyzer_HeavyIons.h:278
JetAnalyzer_HeavyIons::mSumCaloVsPt_n2p043_n1p740
MonitorElement * mSumCaloVsPt_n2p043_n1p740
Definition: JetAnalyzer_HeavyIons.h:343
JetAnalyzer_HeavyIons::mPFCandpT_Barrel_NeutralHadron
MonitorElement * mPFCandpT_Barrel_NeutralHadron
Definition: JetAnalyzer_HeavyIons.h:400
JetAnalyzer_HeavyIons::isPFJet
bool isPFJet
Definition: JetAnalyzer_HeavyIons.h:426
GeomDetEnumerators::isEndcap
bool isEndcap(GeomDetEnumerators::SubDetector m)
Definition: GeomDetEnumerators.cc:37
JetAnalyzer_HeavyIons::mPFCandpT_Endcap_Unknown
MonitorElement * mPFCandpT_Endcap_Unknown
Definition: JetAnalyzer_HeavyIons.h:404
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFpT_30To50
MonitorElement * mPFDeltaR_pTCorrected_PFpT_30To50
Definition: JetAnalyzer_HeavyIons.h:206
JetAnalyzer_HeavyIons::mSumCaloVsPt
MonitorElement * mSumCaloVsPt
Definition: JetAnalyzer_HeavyIons.h:239
JetAnalyzer_HeavyIons::mSumPFPt_n5p191_n2p650
MonitorElement * mSumPFPt_n5p191_n2p650
Definition: JetAnalyzer_HeavyIons.h:308
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFVsInitialpT_300ToInf
MonitorElement * mPFDeltaR_pTCorrected_PFVsInitialpT_300ToInf
Definition: JetAnalyzer_HeavyIons.h:219
JetAnalyzer_HeavyIons::mSumPFVsPt_2p043_2p650
MonitorElement * mSumPFVsPt_2p043_2p650
Definition: JetAnalyzer_HeavyIons.h:305
JetAnalyzer_HeavyIons::mSumCaloVsPt_1p740_2p043
MonitorElement * mSumCaloVsPt_1p740_2p043
Definition: JetAnalyzer_HeavyIons.h:353
ForwardEta
const Double_t ForwardEta
Definition: JetAnalyzer_HeavyIons.h:79
JetAnalyzer_HeavyIons::mNJets
MonitorElement * mNJets
Definition: JetAnalyzer_HeavyIons.h:383
JetAnalyzer_HeavyIons::mPFVsPt
MonitorElement * mPFVsPt
Definition: JetAnalyzer_HeavyIons.h:134
JetAnalyzer_HeavyIons::mSumCaloVsPt_n0p522_0p522
MonitorElement * mSumCaloVsPt_n0p522_0p522
Definition: JetAnalyzer_HeavyIons.h:348
JetAnalyzer_HeavyIons::mSumCaloVsPt_2p650_5p191
MonitorElement * mSumCaloVsPt_2p650_5p191
Definition: JetAnalyzer_HeavyIons.h:355
JetAnalyzer_HeavyIons::UEAlgo
std::string UEAlgo
Definition: JetAnalyzer_HeavyIons.h:100
JetAnalyzer_HeavyIons::mPFCandpT_vs_eta_ChargedHadron
MonitorElement * mPFCandpT_vs_eta_ChargedHadron
Definition: JetAnalyzer_HeavyIons.h:387
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_n5p191_n2p650
MonitorElement * mSumPFVsPtInitial_n5p191_n2p650
Definition: JetAnalyzer_HeavyIons.h:276
JetAnalyzer_HeavyIons::mPFCandpT_Forward_HadE_inHF
MonitorElement * mPFCandpT_Forward_HadE_inHF
Definition: JetAnalyzer_HeavyIons.h:419
JetAnalyzer_HeavyIons::basicJetsToken_
edm::EDGetTokenT< reco::BasicJetCollection > basicJetsToken_
Definition: JetAnalyzer_HeavyIons.h:120
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_HF
MonitorElement * mSumPFVsPtInitial_HF
Definition: JetAnalyzer_HeavyIons.h:269
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_n1p740_n1p479
MonitorElement * mSumCaloVsPtInitial_n1p740_n1p479
Definition: JetAnalyzer_HeavyIons.h:328
JetAnalyzer_HeavyIons::mPFCandpT_Barrel_ChargedHadron
MonitorElement * mPFCandpT_Barrel_ChargedHadron
Definition: JetAnalyzer_HeavyIons.h:396
JetAnalyzer_HeavyIons::mPFCandpT_Barrel_HadE_inHF
MonitorElement * mPFCandpT_Barrel_HadE_inHF
Definition: JetAnalyzer_HeavyIons.h:401
JetAnalyzer_HeavyIons::mPFDeltaR
MonitorElement * mPFDeltaR
Definition: JetAnalyzer_HeavyIons.h:202
reco::PFCandidateCollection
std::vector< reco::PFCandidate > PFCandidateCollection
collection of PFCandidates
Definition: PFCandidateFwd.h:12
JetAnalyzer_HeavyIons::mPFPhi
MonitorElement * mPFPhi
Definition: JetAnalyzer_HeavyIons.h:133
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_15_10
MonitorElement * mPFVsPtInitial_HF_Centrality_15_10
Definition: JetAnalyzer_HeavyIons.h:197
JetAnalyzer_HeavyIons::mSumPFPt_1p479_1p740
MonitorElement * mSumPFPt_1p479_1p740
Definition: JetAnalyzer_HeavyIons.h:319
JetAnalyzer_HeavyIons::mSumCaloPt_2p650_5p191
MonitorElement * mSumCaloPt_2p650_5p191
Definition: JetAnalyzer_HeavyIons.h:371
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_60_55
MonitorElement * mPFVsPtInitial_Barrel_Centrality_60_55
Definition: JetAnalyzer_HeavyIons.h:146
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_30_25
MonitorElement * mPFVsPtInitial_Barrel_Centrality_30_25
Definition: JetAnalyzer_HeavyIons.h:152
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_85_80
MonitorElement * mPFVsPtInitial_Barrel_Centrality_85_80
Definition: JetAnalyzer_HeavyIons.h:141
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_80_75
MonitorElement * mPFVsPtInitial_HF_Centrality_80_75
Definition: JetAnalyzer_HeavyIons.h:184
JetAnalyzer_HeavyIons::mInputPFCandCollection
edm::InputTag mInputPFCandCollection
Definition: JetAnalyzer_HeavyIons.h:96
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_n1p131_n0p783
MonitorElement * mSumPFVsPtInitial_n1p131_n0p783
Definition: JetAnalyzer_HeavyIons.h:281
JetAnalyzer_HeavyIons::mCaloVsPt
MonitorElement * mCaloVsPt
Definition: JetAnalyzer_HeavyIons.h:227
edm::Handle::product
T const * product() const
Definition: Handle.h:81
JetAnalyzer_HeavyIons::backgrounds_
edm::EDGetTokenT< edm::ValueMap< reco::VoronoiBackground > > backgrounds_
Definition: JetAnalyzer_HeavyIons.h:126
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
JetAnalyzer_HeavyIons::mSumCaloPt
MonitorElement * mSumCaloPt
Definition: JetAnalyzer_HeavyIons.h:241
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_95_90
MonitorElement * mPFVsPtInitial_HF_Centrality_95_90
Definition: JetAnalyzer_HeavyIons.h:181
JetAnalyzer_HeavyIons::mSumSquaredPFPt_eta
MonitorElement * mSumSquaredPFPt_eta
Definition: JetAnalyzer_HeavyIons.h:253
JetAnalyzer_HeavyIons::mSumPFPt_n1p131_n0p783
MonitorElement * mSumPFPt_n1p131_n0p783
Definition: JetAnalyzer_HeavyIons.h:313
JetAnalyzer_HeavyIons::mPFCandpT_Endcap_EME_inHF
MonitorElement * mPFCandpT_Endcap_EME_inHF
Definition: JetAnalyzer_HeavyIons.h:411
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_n5p191_n2p650
MonitorElement * mSumCaloVsPtInitial_n5p191_n2p650
Definition: JetAnalyzer_HeavyIons.h:325
JetAnalyzer_HeavyIons::mEnergy
MonitorElement * mEnergy
Definition: JetAnalyzer_HeavyIons.h:378
JetAnalyzer_HeavyIons::mPFCandpT_vs_eta_HadE_inHF
MonitorElement * mPFCandpT_vs_eta_HadE_inHF
Definition: JetAnalyzer_HeavyIons.h:392
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_65_60
MonitorElement * mPFVsPtInitial_HF_Centrality_65_60
Definition: JetAnalyzer_HeavyIons.h:187
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_100_95
MonitorElement * mPFVsPtInitial_HF_Centrality_100_95
Definition: JetAnalyzer_HeavyIons.h:180
RecoJets_cff.recoJets
tuple recoJets
Definition: RecoJets_cff.py:56
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial
MonitorElement * mSumCaloVsPtInitial
Definition: JetAnalyzer_HeavyIons.h:240
JetAnalyzer_HeavyIons::pfCandViewToken_
edm::EDGetTokenT< reco::CandidateView > pfCandViewToken_
Definition: JetAnalyzer_HeavyIons.h:123
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_eta
MonitorElement * mSumCaloVsPtInitial_eta
Definition: JetAnalyzer_HeavyIons.h:243
JetAnalyzer_HeavyIons::caloTowersToken_
edm::EDGetTokenT< CaloTowerCollection > caloTowersToken_
Definition: JetAnalyzer_HeavyIons.h:117
JetAnalyzer_HeavyIons::mPFVsPtInitialDeltaR_pTCorrected
MonitorElement * mPFVsPtInitialDeltaR_pTCorrected
Definition: JetAnalyzer_HeavyIons.h:221
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_0p783_1p131
MonitorElement * mSumCaloVsPtInitial_0p783_1p131
Definition: JetAnalyzer_HeavyIons.h:334
JetAnalyzer_HeavyIons::mSumPFPt_1p740_2p043
MonitorElement * mSumPFPt_1p740_2p043
Definition: JetAnalyzer_HeavyIons.h:320
JetAnalyzer_HeavyIons::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: JetAnalyzer_HeavyIons.cc:357
JetAnalyzer_HeavyIons::mPFCandpT_Forward_Unknown
MonitorElement * mPFCandpT_Forward_Unknown
Definition: JetAnalyzer_HeavyIons.h:413
JetAnalyzer_HeavyIons::mSumCaloPt_n5p191_n2p650
MonitorElement * mSumCaloPt_n5p191_n2p650
Definition: JetAnalyzer_HeavyIons.h:357
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_60_55
MonitorElement * mPFVsPtInitial_HF_Centrality_60_55
Definition: JetAnalyzer_HeavyIons.h:188
JetAnalyzer_HeavyIons::mPFPt
MonitorElement * mPFPt
Definition: JetAnalyzer_HeavyIons.h:131
JetAnalyzer_HeavyIons::mSumPFPt_n0p783_n0p522
MonitorElement * mSumPFPt_n0p783_n0p522
Definition: JetAnalyzer_HeavyIons.h:314
JetAnalyzer_HeavyIons::mSumCaloPt_2p043_2p650
MonitorElement * mSumCaloPt_2p043_2p650
Definition: JetAnalyzer_HeavyIons.h:370
JetAnalyzer_HeavyIons::mHF
MonitorElement * mHF
Definition: JetAnalyzer_HeavyIons.h:265
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:43
JetAnalyzer_HeavyIons::mSumCaloPt_n1p131_n0p783
MonitorElement * mSumCaloPt_n1p131_n0p783
Definition: JetAnalyzer_HeavyIons.h:362
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_45_40
MonitorElement * mPFVsPtInitial_Barrel_Centrality_45_40
Definition: JetAnalyzer_HeavyIons.h:149
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_50_45
MonitorElement * mPFVsPtInitial_Barrel_Centrality_50_45
Definition: JetAnalyzer_HeavyIons.h:148
JetAnalyzer_HeavyIons::mSumPFPt_n2p043_n1p740
MonitorElement * mSumPFPt_n2p043_n1p740
Definition: JetAnalyzer_HeavyIons.h:310
reco::VoronoiBackground::pt_subtracted
double pt_subtracted() const
Definition: VoronoiBackground.h:18
reco::PFCandidate
Particle reconstructed by the particle flow algorithm.
Definition: PFCandidate.h:39
JetAnalyzer_HeavyIons::mSumPFVsPt
MonitorElement * mSumPFVsPt
Definition: JetAnalyzer_HeavyIons.h:231
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_5_0
MonitorElement * mPFVsPtInitial_EndCap_Centrality_5_0
Definition: JetAnalyzer_HeavyIons.h:178
JetAnalyzer_HeavyIons::mPFCandpT_Barrel_electron
MonitorElement * mPFCandpT_Barrel_electron
Definition: JetAnalyzer_HeavyIons.h:397
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFVsInitialpT_180To300
MonitorElement * mPFDeltaR_pTCorrected_PFVsInitialpT_180To300
Definition: JetAnalyzer_HeavyIons.h:218
position
static int position[264][3]
Definition: ReadPGInfo.cc:509
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_45_40
MonitorElement * mPFVsPtInitial_EndCap_Centrality_45_40
Definition: JetAnalyzer_HeavyIons.h:170
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_n0p783_n0p522
MonitorElement * mSumCaloVsPtInitial_n0p783_n0p522
Definition: JetAnalyzer_HeavyIons.h:331
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_5_0
MonitorElement * mPFVsPtInitial_HF_Centrality_5_0
Definition: JetAnalyzer_HeavyIons.h:199
JetAnalyzer_HeavyIons::mPFCandpT_Endcap_electron
MonitorElement * mPFCandpT_Endcap_electron
Definition: JetAnalyzer_HeavyIons.h:406
JetAnalyzer_HeavyIons::mNPFpart
MonitorElement * mNPFpart
Definition: JetAnalyzer_HeavyIons.h:130
JetAnalyzer_HeavyIons::nedge_pseudorapidity
static const Int_t nedge_pseudorapidity
Definition: JetAnalyzer_HeavyIons.h:431
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_20_15
MonitorElement * mPFVsPtInitial_EndCap_Centrality_20_15
Definition: JetAnalyzer_HeavyIons.h:175
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_40_35
MonitorElement * mPFVsPtInitial_Barrel_Centrality_40_35
Definition: JetAnalyzer_HeavyIons.h:150
JetAnalyzer_HeavyIons::mSumPFPt_eta
MonitorElement * mSumPFPt_eta
Definition: JetAnalyzer_HeavyIons.h:237
JetAnalyzer_HeavyIons::mSumCaloVsPtInitial_2p650_5p191
MonitorElement * mSumCaloVsPtInitial_2p650_5p191
Definition: JetAnalyzer_HeavyIons.h:339
JetAnalyzer_HeavyIons::mNJets_40
MonitorElement * mNJets_40
Definition: JetAnalyzer_HeavyIons.h:384
relval_parameters_module.energy
string energy
Definition: relval_parameters_module.py:29
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_90_85
MonitorElement * mPFVsPtInitial_EndCap_Centrality_90_85
Definition: JetAnalyzer_HeavyIons.h:161
BarrelEta
const Double_t BarrelEta
Definition: JetAnalyzer_HeavyIons.h:77
JetAnalyzer_HeavyIons::mRecoJetPtThreshold
double mRecoJetPtThreshold
Definition: JetAnalyzer_HeavyIons.h:102
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_0p522_0p783
MonitorElement * mSumPFVsPtInitial_0p522_0p783
Definition: JetAnalyzer_HeavyIons.h:284
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_35_30
MonitorElement * mPFVsPtInitial_Barrel_Centrality_35_30
Definition: JetAnalyzer_HeavyIons.h:151
JetAnalyzer_HeavyIons::JetType
std::string JetType
Definition: JetAnalyzer_HeavyIons.h:99
JetAnalyzer_HeavyIons::mSumPFVsPt_0p783_1p131
MonitorElement * mSumPFVsPt_0p783_1p131
Definition: JetAnalyzer_HeavyIons.h:301
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_10_5
MonitorElement * mPFVsPtInitial_HF_Centrality_10_5
Definition: JetAnalyzer_HeavyIons.h:198
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_70_65
MonitorElement * mPFVsPtInitial_Barrel_Centrality_70_65
Definition: JetAnalyzer_HeavyIons.h:144
JetAnalyzer_HeavyIons::mJetArea
MonitorElement * mJetArea
Definition: JetAnalyzer_HeavyIons.h:381
reco::PFCandidate::particleId
virtual ParticleType particleId() const
Definition: PFCandidate.h:373
JetAnalyzer_HeavyIons::mSumCaloPt_n1p479_n1p131
MonitorElement * mSumCaloPt_n1p479_n1p131
Definition: JetAnalyzer_HeavyIons.h:361
JetAnalyzer_HeavyIons::mPFDeltaR_pTCorrected_PFpT_120To180
MonitorElement * mPFDeltaR_pTCorrected_PFpT_120To180
Definition: JetAnalyzer_HeavyIons.h:209
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_n1p740_n1p479
MonitorElement * mSumPFVsPtInitial_n1p740_n1p479
Definition: JetAnalyzer_HeavyIons.h:279
JetAnalyzer_HeavyIons::JetAnalyzer_HeavyIons
JetAnalyzer_HeavyIons(const edm::ParameterSet &)
Definition: JetAnalyzer_HeavyIons.cc:18
relval_steps.k
list k
Definition: relval_steps.py:1462
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_1p740_2p043
MonitorElement * mSumPFVsPtInitial_1p740_2p043
Definition: JetAnalyzer_HeavyIons.h:288
tempVector
vector< int > tempVector
Definition: DBlmapReader.cc:10
JetAnalyzer_HeavyIons::mPFCandpT_Barrel_EME_inHF
MonitorElement * mPFCandpT_Barrel_EME_inHF
Definition: JetAnalyzer_HeavyIons.h:402
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_40_35
MonitorElement * mPFVsPtInitial_EndCap_Centrality_40_35
Definition: JetAnalyzer_HeavyIons.h:171
pfJets_cff.pfJets
tuple pfJets
Definition: pfJets_cff.py:8
reco::LeafCandidate::phi
virtual double phi() const
momentum azimuthal angle
Definition: LeafCandidate.h:133
JetAnalyzer_HeavyIons::mjetpileup
MonitorElement * mjetpileup
Definition: JetAnalyzer_HeavyIons.h:382
JetAnalyzer_HeavyIons::mSumPFVsPt_n1p740_n1p479
MonitorElement * mSumPFVsPt_n1p740_n1p479
Definition: JetAnalyzer_HeavyIons.h:295
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_35_30
MonitorElement * mPFVsPtInitial_HF_Centrality_35_30
Definition: JetAnalyzer_HeavyIons.h:193
JetAnalyzer_HeavyIons::mSumSquaredPFVsPtInitial_eta
MonitorElement * mSumSquaredPFVsPtInitial_eta
Definition: JetAnalyzer_HeavyIons.h:251
JetAnalyzer_HeavyIons::mSumPFVsPt_n1p131_n0p783
MonitorElement * mSumPFVsPt_n1p131_n0p783
Definition: JetAnalyzer_HeavyIons.h:297
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_2p043_2p650
MonitorElement * mSumPFVsPtInitial_2p043_2p650
Definition: JetAnalyzer_HeavyIons.h:289
JetAnalyzer_HeavyIons::mSumCaloVsPt_eta
MonitorElement * mSumCaloVsPt_eta
Definition: JetAnalyzer_HeavyIons.h:244
JetAnalyzer_HeavyIons::mSumCaloVsPt_0p522_0p783
MonitorElement * mSumCaloVsPt_0p522_0p783
Definition: JetAnalyzer_HeavyIons.h:349
JetAnalyzer_HeavyIons::mPFCandpT_vs_eta_Unknown
MonitorElement * mPFCandpT_vs_eta_Unknown
Definition: JetAnalyzer_HeavyIons.h:386
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_65_60
MonitorElement * mPFVsPtInitial_Barrel_Centrality_65_60
Definition: JetAnalyzer_HeavyIons.h:145
JetAnalyzer_HeavyIons::mPFCandpT_Forward_NeutralHadron
MonitorElement * mPFCandpT_Forward_NeutralHadron
Definition: JetAnalyzer_HeavyIons.h:418
JetAnalyzer_HeavyIons::mPFVsPtInitial_Barrel_Centrality_80_75
MonitorElement * mPFVsPtInitial_Barrel_Centrality_80_75
Definition: JetAnalyzer_HeavyIons.h:142
JetAnalyzer_HeavyIons::~JetAnalyzer_HeavyIons
virtual ~JetAnalyzer_HeavyIons()
Definition: JetAnalyzer_HeavyIons.cc:684
JetAnalyzer_HeavyIons::mNCalopart
MonitorElement * mNCalopart
Definition: JetAnalyzer_HeavyIons.h:223
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_eta
MonitorElement * mSumPFVsPtInitial_eta
Definition: JetAnalyzer_HeavyIons.h:235
edm::Run
Definition: Run.h:41
JetAnalyzer_HeavyIons::mPFVsPtInitial_EndCap_Centrality_85_80
MonitorElement * mPFVsPtInitial_EndCap_Centrality_85_80
Definition: JetAnalyzer_HeavyIons.h:162
JetAnalyzer_HeavyIons::mSumCaloVsPt_0p783_1p131
MonitorElement * mSumCaloVsPt_0p783_1p131
Definition: JetAnalyzer_HeavyIons.h:350
JetAnalyzer_HeavyIons::mSumCaloPt_1p740_2p043
MonitorElement * mSumCaloPt_1p740_2p043
Definition: JetAnalyzer_HeavyIons.h:369
JetAnalyzer_HeavyIons::mSumPFVsPtInitial_n0p522_0p522
MonitorElement * mSumPFVsPtInitial_n0p522_0p522
Definition: JetAnalyzer_HeavyIons.h:283
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_70_65
MonitorElement * mPFVsPtInitial_HF_Centrality_70_65
Definition: JetAnalyzer_HeavyIons.h:186
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_30_25
MonitorElement * mPFVsPtInitial_HF_Centrality_30_25
Definition: JetAnalyzer_HeavyIons.h:194
JetAnalyzer_HeavyIons::mPFVsPtInitial_HF_Centrality_55_50
MonitorElement * mPFVsPtInitial_HF_Centrality_55_50
Definition: JetAnalyzer_HeavyIons.h:189
JetAnalyzer_HeavyIons::mSumCaloVsPt_2p043_2p650
MonitorElement * mSumCaloVsPt_2p043_2p650
Definition: JetAnalyzer_HeavyIons.h:354
JetAnalyzer_HeavyIons::mPFCandpT_Barrel_muon
MonitorElement * mPFCandpT_Barrel_muon
Definition: JetAnalyzer_HeavyIons.h:398
Generated for CMSSW Reference Manual by doxygen 1.8.5