104 outputMEsInRootFile = pSet.
getParameter<
bool>(
"OutputMEsInRootFile");
112 fill_jet_high_level_histo=pSet.
getParameter<
bool>(
"filljetHighLevel"),
118 if (isCaloJet_) caloJetsToken_ = consumes<reco::CaloJetCollection>(mInputCollection_);
120 if (isPFJet_) pfJetsToken_ = consumes<reco::PFJetCollection>(mInputCollection_);
122 JetIDQuality_ = pSet.
getParameter<
string>(
"JetIDQuality");
123 JetIDVersion_ = pSet.
getParameter<
string>(
"JetIDVersion");
128 jetID_ValueMapToken_= consumes< edm::ValueMap<reco::JetID> >(inputJetIDValueMap);
129 if(JetIDVersion_==
"PURE09"){
131 }
else if (JetIDVersion_==
"DQM09"){
133 }
else if (JetIDVersion_==
"CRAFT08"){
136 if (verbose_)
std::cout<<
"no Valid JetID version given"<<std::endl;
138 if(JetIDQuality_==
"MINIMAL"){
140 }
else if (JetIDQuality_==
"LOOSE_AOD"){
142 }
else if (JetIDQuality_==
"LOOSE"){
144 }
else if (JetIDQuality_==
"TIGHT"){
147 if (verbose_)
std::cout<<
"no Valid JetID quality given"<<std::endl;
155 if(JetIDVersion_==
"FIRSTDATA"){
158 if (verbose_)
std::cout<<
"no valid PF JetID version given"<<std::endl;
160 if (JetIDQuality_==
"LOOSE"){
162 }
else if (JetIDQuality_==
"TIGHT"){
165 if (verbose_)
std::cout<<
"no Valid PFJetID quality given"<<std::endl;
174 ptThresholdUnc_ = 20.;
175 asymmetryThirdJetCut_ = 5.;
176 balanceThirdJetCut_ = 0.2;
185 jetCleaningFlag_ =
false;
194 DCSFilterForDCSMonitoring_ =
new JetMETDQMDCSFilter(
"ecal:hbhe:hf:ho:pixel:sistrip:es:muon", iC);
203 highPtJetExpr_ = highptjetparms.
getParameter<std::vector<std::string> >(
"hltPaths");
204 lowPtJetExpr_ = lowptjetparms .getParameter<std::vector<std::string> >(
"hltPaths");
211 bypassAllPVChecks_= cleaningParameters_.
getParameter<
bool>(
"bypassAllPVChecks");
212 vertexLabel_ = cleaningParameters_.getParameter<
edm::InputTag>(
"vertexCollection");
213 vertexToken_ = consumes<std::vector<reco::Vertex> >(
edm::InputTag(vertexLabel_));
215 gtLabel_ = cleaningParameters_.getParameter<
edm::InputTag>(
"gtLabel");
216 gtToken_ = consumes<L1GlobalTriggerReadoutRecord>(
edm::InputTag(gtLabel_));
218 std::string inputCollectionLabel(mInputCollection_.label());
219 verbose_= parameters_.getParameter<
int>(
"verbose");
221 etaBin_ = parameters_.getParameter<
int>(
"etaBin");
222 etaMin_ = parameters_.getParameter<
double>(
"etaMin");
223 etaMax_ = parameters_.getParameter<
double>(
"etaMax");
225 phiBin_ = parameters_.getParameter<
int>(
"phiBin");
226 phiMin_ = parameters_.getParameter<
double>(
"phiMin");
227 phiMax_ = parameters_.getParameter<
double>(
"phiMax");
229 ptBin_ = parameters_.getParameter<
int>(
"ptBin");
230 ptMin_ = parameters_.getParameter<
double>(
"ptMin");
231 ptMax_ = parameters_.getParameter<
double>(
"ptMax");
233 eBin_ = parameters_.getParameter<
int>(
"eBin");
234 eMin_ = parameters_.getParameter<
double>(
"eMin");
235 eMax_ = parameters_.getParameter<
double>(
"eMax");
237 pBin_ = parameters_.getParameter<
int>(
"pBin");
238 pMin_ = parameters_.getParameter<
double>(
"pMin");
239 pMax_ = parameters_.getParameter<
double>(
"pMax");
241 nbinsPV_ = parameters_.getParameter<
int>(
"pVBin");
242 nPVlow_ = parameters_.getParameter<
double>(
"pVMin");
243 nPVhigh_ = parameters_.getParameter<
double>(
"pVMax");
245 ptThreshold_ = parameters_.getParameter<
double>(
"ptThreshold");
246 ptThresholdUnc_=parameters_.getParameter<
double>(
"ptThresholdUnc");
247 asymmetryThirdJetCut_ = parameters_.getParameter<
double>(
"asymmetryThirdJetCut");
248 balanceThirdJetCut_ = parameters_.getParameter<
double>(
"balanceThirdJetCut");
255 delete highPtJetEventFlag_;
256 delete lowPtJetEventFlag_;
258 delete DCSFilterForDCSMonitoring_;
259 delete DCSFilterForJetMonitoring_;
269 if(jetCleaningFlag_){
271 DirName =
"JetMET/Jet/Cleaned"+mInputCollection_.label();
274 DirName =
"JetMET/Jet/Uncleaned"+mInputCollection_.label();
277 jetME = ibooker.
book1D(
"jetReco",
"jetReco", 3, 1, 4);
279 jetME->setBinLabel(2,
"PFJets",1);
280 jetME->setBinLabel(3,
"JPTJets",1);
282 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"jetReco",jetME));
285 mPt = ibooker.
book1D(
"Pt",
"pt", ptBin_, ptMin_, ptMax_);
286 mEta = ibooker.
book1D(
"Eta",
"eta", etaBin_, etaMin_, etaMax_);
287 mPhi = ibooker.
book1D(
"Phi",
"phi", phiBin_, phiMin_, phiMax_);
288 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt" ,mPt));
289 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Eta",mEta));
290 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi",mPhi));
293 mConstituents = ibooker.
book1D(
"Constituents",
"# of constituents", 50, 0, 100);
294 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Constituents",mConstituents));
296 mJetEnergyCorr= ibooker.
book1D(
"JetEnergyCorr",
"jet energy correction factor", 50, 0.0,3.0);
297 mJetEnergyCorrVSEta= ibooker.
bookProfile(
"JetEnergyCorrVSEta",
"jet energy correction factor VS eta", etaBin_, etaMin_,etaMax_, 0.0,3.0);
298 mJetEnergyCorrVSPt= ibooker.
bookProfile(
"JetEnergyCorrVSPt",
"jet energy correction factor VS pt", ptBin_, ptMin_,ptMax_, 0.0,3.0);
299 mHFrac = ibooker.
book1D(
"HFrac",
"HFrac", 140, -0.2, 1.2);
300 mEFrac = ibooker.
book1D(
"EFrac",
"EFrac", 52, -0.02, 1.02);
301 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"JetEnergyCorr" ,mJetEnergyCorr));
302 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"JetEnergyCorrVSEta" ,mJetEnergyCorrVSEta));
303 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"JetEnergyCorrVSPt" ,mJetEnergyCorrVSPt));
304 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFrac" ,mHFrac));
305 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"EFrac" ,mEFrac));
307 mPt_uncor = ibooker.
book1D(
"Pt_uncor",
"pt for uncorrected jets", ptBin_, ptThresholdUnc_, ptMax_);
308 mEta_uncor = ibooker.
book1D(
"Eta_uncor",
"eta for uncorrected jets", etaBin_, etaMin_, etaMax_);
309 mPhi_uncor = ibooker.
book1D(
"Phi_uncor",
"phi for uncorrected jets", phiBin_, phiMin_, phiMax_);
310 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_uncor" ,mPt_uncor));
311 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Eta_uncor",mEta_uncor));
312 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_uncor",mPhi_uncor));
314 mConstituents_uncor = ibooker.
book1D(
"Constituents_uncor",
"# of constituents for uncorrected jets", 50, 0, 100);
315 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Constituents_uncor",mConstituents_uncor));
318 mDPhi = ibooker.
book1D(
"DPhi",
"dPhi btw the two leading jets", 100, 0., acos(-1.));
319 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"DPhi",mDPhi));
323 mPt_profile = ibooker.
bookProfile(
"Pt_profile",
"pt", nbinsPV_, nPVlow_, nPVhigh_, ptBin_, ptMin_, ptMax_);
324 mEta_profile = ibooker.
bookProfile(
"Eta_profile",
"eta", nbinsPV_, nPVlow_, nPVhigh_, etaBin_, etaMin_, etaMax_);
325 mPhi_profile = ibooker.
bookProfile(
"Phi_profile",
"phi", nbinsPV_, nPVlow_, nPVhigh_, phiBin_, phiMin_, phiMax_);
327 mConstituents_profile = ibooker.
bookProfile(
"Constituents_profile",
"# of constituents", nbinsPV_, nPVlow_, nPVhigh_, 50, 0, 100);
329 mHFrac_profile = ibooker.
bookProfile(
"HFrac_profile",
"HFrac", nbinsPV_, nPVlow_, nPVhigh_, 140, -0.2, 1.2);
330 mEFrac_profile = ibooker.
bookProfile(
"EFrac_profile",
"EFrac", nbinsPV_, nPVlow_, nPVhigh_, 52, -0.02, 1.02);
331 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_profile" ,mPt_profile));
332 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Eta_profile",mEta_profile));
333 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_profile",mPhi_profile));
334 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFrac_profile",mHFrac_profile));
335 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"EFrac_profile",mEFrac_profile));
338 if(!jetCleaningFlag_){
339 mLooseJIDPassFractionVSeta = ibooker.
bookProfile(
"JetIDPassFractionVSeta",
"JetIDPassFractionVSeta",etaBin_, etaMin_, etaMax_,0.,1.2);
340 mLooseJIDPassFractionVSpt = ibooker.
bookProfile(
"JetIDPassFractionVSpt",
"JetIDPassFractionVSpt",ptBin_, ptMin_, ptMax_,0.,1.2);
341 mLooseJIDPassFractionVSptNoHF = ibooker.
bookProfile(
"JetIDPassFractionVSptNoHF",
"JetIDPassFractionVSptNoHF",ptBin_, ptMin_, ptMax_,0.,1.2);
342 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"JetIDPassFractionVSeta" ,mLooseJIDPassFractionVSeta));
343 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"JetIDPassFractionVSpt" ,mLooseJIDPassFractionVSpt));
344 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"JetIDPassFractionVSptNoHF",mLooseJIDPassFractionVSptNoHF));
347 mNJets_profile = ibooker.
bookProfile(
"NJets_profile",
"number of jets", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
353 mEta_profile ->setAxisTitle(
"nvtx",1);
354 mPhi_profile ->setAxisTitle(
"nvtx",1);
356 mConstituents_profile->setAxisTitle(
"nvtx",1);
357 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Constituents_profile",mConstituents_profile));
359 mHFrac_profile ->setAxisTitle(
"nvtx",1);
360 mEFrac_profile ->setAxisTitle(
"nvtx",1);
362 mNJets_profile->setAxisTitle(
"nvtx",1);
364 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_profile" ,mPt_profile));
365 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Eta_profile",mEta_profile));
366 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_profile",mPhi_profile));
367 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFrac_profile" ,mHFrac_profile));
368 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"EFrac_profile" ,mEFrac_profile));
369 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NJets_profile" ,mNJets_profile));
372 mPhiVSEta = ibooker.
book2D(
"PhiVSEta",
"PhiVSEta", 50, etaMin_, etaMax_, 24, phiMin_, phiMax_);
373 mPhiVSEta->
getTH2F()->SetOption(
"colz");
374 mPhiVSEta->setAxisTitle(
"#eta",1);
375 mPhiVSEta->setAxisTitle(
"#phi",2);
376 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhiVSEta" ,mPhiVSEta));
378 mPt_1 = ibooker.
book1D(
"Pt_1",
"Pt spectrum of jets - range 1", 20, 0, 100);
379 mPt_2 = ibooker.
book1D(
"Pt_2",
"Pt spectrum of jets - range 2", 60, 0, 300);
380 mPt_3 = ibooker.
book1D(
"Pt_3",
"Pt spectrum of jets - range 3", 100, 0, 5000);
382 mPt_Lo = ibooker.
book1D(
"Pt_Lo",
"Pt (Pass Low Pt Jet Trigger)", 20, 0, 100);
384 mPhi_Lo = ibooker.
book1D(
"Phi_Lo",
"Phi (Pass Low Pt Jet Trigger)", phiBin_, phiMin_, phiMax_);
386 mPt_Hi = ibooker.
book1D(
"Pt_Hi",
"Pt (Pass Hi Pt Jet Trigger)", 60, 0, 300);
387 mEta_Hi = ibooker.
book1D(
"Eta_Hi",
"Eta (Pass Hi Pt Jet Trigger)", etaBin_, etaMin_, etaMax_);
388 mPhi_Hi = ibooker.
book1D(
"Phi_Hi",
"Phi (Pass Hi Pt Jet Trigger)", phiBin_, phiMin_, phiMax_);
389 mNJets = ibooker.
book1D(
"NJets",
"number of jets", 100, 0, 100);
391 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_1" ,mPt_1));
392 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_2" ,mPt_2));
393 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_3" ,mPt_3));
394 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_Lo" ,mPt_Lo));
395 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_Lo" ,mPhi_Lo));
396 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_Hi" ,mPt_Hi));
397 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Eta_Hi" ,mEta_Hi));
398 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_Hi" ,mPhi_Hi));
399 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NJets" ,mNJets));
404 mConstituents_Barrel = ibooker.
book1D(
"Constituents_Barrel",
"Constituents Barrel", 50, 0, 100);
405 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Constituents_Barrel",mConstituents_Barrel));
407 mHFrac_Barrel = ibooker.
book1D(
"HFrac_Barrel",
"HFrac Barrel", 100, 0, 1);
408 mEFrac_Barrel = ibooker.
book1D(
"EFrac_Barrel",
"EFrac Barrel", 52, -0.02, 1.02);
409 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFrac_Barrel" ,mHFrac_Barrel));
410 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"EFrac_Barrel" ,mEFrac_Barrel));
415 mConstituents_EndCap = ibooker.
book1D(
"Constituents_EndCap",
"Constituents EndCap", 50, 0, 100);
416 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Constituents_EndCap",mConstituents_EndCap));
418 mHFrac_EndCap = ibooker.
book1D(
"HFrac_EndCap",
"HFrac EndCap", 100, 0, 1);
419 mEFrac_EndCap = ibooker.
book1D(
"EFrac_EndCap",
"EFrac EndCap", 52, -0.02, 1.02);
420 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFrac_EndCap" ,mHFrac_EndCap));
421 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"EFrac_EndCap" ,mEFrac_EndCap));
426 mConstituents_Forward = ibooker.
book1D(
"Constituents_Forward",
"Constituents Forward", 50, 0, 100);
427 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Constituents_Forward",mConstituents_Forward));
429 mHFrac_Forward = ibooker.
book1D(
"HFrac_Forward",
"HFrac Forward", 140, -0.2, 1.2);
430 mEFrac_Forward = ibooker.
book1D(
"EFrac_Forward",
"EFrac Forward", 52, -0.02, 1.02);
431 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFrac_Forward" ,mHFrac_Forward));
432 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"EFrac_Forward" ,mEFrac_Forward));
435 mPt_Barrel_Hi = ibooker.
book1D(
"Pt_Barrel_Hi",
"Pt Barrel (Pass Hi Pt Jet Trigger)", 60, 0, 300);
436 mPhi_Barrel_Hi = ibooker.
book1D(
"Phi_Barrel_Hi",
"Phi Barrel (Pass Hi Pt Jet Trigger)", phiBin_, phiMin_, phiMax_);
438 mPt_EndCap_Hi = ibooker.
book1D(
"Pt_EndCap_Hi",
"Pt EndCap (Pass Hi Pt Jet Trigger)", 60, 0, 300);
439 mPhi_EndCap_Hi = ibooker.
book1D(
"Phi_EndCap_Hi",
"Phi EndCap (Pass Hi Pt Jet Trigger)", phiBin_, phiMin_, phiMax_);
441 mPt_Forward_Hi = ibooker.
book1D(
"Pt_Forward_Hi",
"Pt Forward (Pass Hi Pt Jet Trigger)", 60, 0, 300);
442 mPhi_Forward_Hi = ibooker.
book1D(
"Phi_Forward_Hi",
"Phi Forward (Pass Hi Pt Jet Trigger)", phiBin_, phiMin_, phiMax_);
444 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_Barrel_Hi" ,mPt_Barrel_Hi));
445 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_Barrel_Hi",mPhi_Barrel_Hi));
446 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_EndCap_Hi" ,mPt_EndCap_Hi));
447 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_EndCap_Hi",mPhi_EndCap_Hi));
448 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_Forward_Hi" ,mPt_Forward_Hi));
449 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_Forward_Hi",mPhi_Forward_Hi));
451 mPhi_Barrel = ibooker.
book1D(
"Phi_Barrel",
"Phi_Barrel", phiBin_, phiMin_, phiMax_);
452 mPt_Barrel = ibooker.
book1D(
"Pt_Barrel",
"Pt_Barrel", ptBin_, ptMin_, ptMax_);
454 mPhi_EndCap = ibooker.
book1D(
"Phi_EndCap",
"Phi_EndCap", phiBin_, phiMin_, phiMax_);
455 mPt_EndCap = ibooker.
book1D(
"Pt_EndCap",
"Pt_EndCap", ptBin_, ptMin_, ptMax_);
457 mPhi_Forward = ibooker.
book1D(
"Phi_Forward",
"Phi_Forward", phiBin_, phiMin_, phiMax_);
458 mPt_Forward = ibooker.
book1D(
"Pt_Forward",
"Pt_Forward", ptBin_, ptMin_, ptMax_);
460 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_Barrel" ,mPt_Barrel));
461 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_Barrel",mPhi_Barrel));
462 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_EndCap" ,mPt_EndCap));
463 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_EndCap",mPhi_EndCap));
464 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Pt_Forward" ,mPt_Forward));
465 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"Phi_Forward",mPhi_Forward));
468 mEtaFirst = ibooker.
book1D(
"EtaFirst",
"EtaFirst", 100, -5, 5);
469 mPhiFirst = ibooker.
book1D(
"PhiFirst",
"PhiFirst", 70, -3.5, 3.5);
470 mPtFirst = ibooker.
book1D(
"PtFirst",
"PtFirst", ptBin_, ptMin_, ptMax_);
472 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"EtaFirst" ,mEtaFirst));
473 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PtFirst" ,mPtFirst));
474 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhiFirst" ,mPhiFirst));
480 mMaxEInEmTowers = ibooker.
book1D(
"MaxEInEmTowers",
"MaxEInEmTowers", 150, 0, 150);
481 mMaxEInHadTowers = ibooker.
book1D(
"MaxEInHadTowers",
"MaxEInHadTowers", 150, 0, 150);
483 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MaxEInEmTowers" ,mMaxEInEmTowers));
484 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MaxEInHadTowers" ,mMaxEInHadTowers));
486 mHadEnergyInHO = ibooker.
book1D(
"HadEnergyInHO",
"HadEnergyInHO", 100, 0, 20);
487 mHadEnergyInHB = ibooker.
book1D(
"HadEnergyInHB",
"HadEnergyInHB", 100, 0, 100);
488 mHadEnergyInHF = ibooker.
book1D(
"HadEnergyInHF",
"HadEnergyInHF", 100, 0, 100);
489 mHadEnergyInHE = ibooker.
book1D(
"HadEnergyInHE",
"HadEnergyInHE", 100, 0, 200);
490 mEmEnergyInEB = ibooker.
book1D(
"EmEnergyInEB",
"EmEnergyInEB", 100, 0, 100);
491 mEmEnergyInEE = ibooker.
book1D(
"EmEnergyInEE",
"EmEnergyInEE", 100, 0, 100);
492 mEmEnergyInHF = ibooker.
book1D(
"EmEnergyInHF",
"EmEnergyInHF", 120, -20, 200);
493 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HadEnergyInHO" ,mHadEnergyInHO));
494 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HadEnergyInHB" ,mHadEnergyInHB));
495 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HadEnergyInHF" ,mHadEnergyInHF));
496 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HadEnergyInHE" ,mHadEnergyInHE));
497 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"EmEnergyInEB" ,mEmEnergyInEB));
498 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"EmEnergyInEE" ,mEmEnergyInEE));
499 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"EmEnergyInHF" ,mEmEnergyInHF));
502 mresEMF = ibooker.
book1D(
"resEMF",
"resEMF", 50, 0., 1.);
503 mN90Hits = ibooker.
book1D(
"N90Hits",
"N90Hits", 100, 0., 100);
504 mfHPD = ibooker.
book1D(
"fHPD",
"fHPD", 50, 0., 1.);
505 mfRBX = ibooker.
book1D(
"fRBX",
"fRBX", 50, 0., 1.);
506 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"resEMF" ,mresEMF));
507 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"N90Hits" ,mN90Hits));
508 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"fHPD" ,mfHPD));
509 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"fRBX" ,mfRBX));
698 mCHFracVSeta_lowPt= ibooker.
bookProfile(
"CHFracVSeta_lowPt",
"CHFracVSeta_lowPt",etaBin_, etaMin_, etaMax_,0.,1.2);
699 mNHFracVSeta_lowPt= ibooker.
bookProfile(
"NHFracVSeta_lowPt",
"NHFracVSeta_lowPt",etaBin_, etaMin_, etaMax_,0.,1.2);
700 mPhFracVSeta_lowPt= ibooker.
bookProfile(
"PhFracVSeta_lowPt",
"PhFracVSeta_lowPt",etaBin_, etaMin_, etaMax_,0.,1.2);
701 mElFracVSeta_lowPt= ibooker.
bookProfile(
"ElFracVSeta_lowPt",
"ElFracVSeta_lowPt",etaBin_, etaMin_, etaMax_,0.,1.2);
702 mMuFracVSeta_lowPt= ibooker.
bookProfile(
"MuFracVSeta_lowPt",
"MuFracVSeta_lowPt",etaBin_, etaMin_, etaMax_,0.,1.2);
703 mCHFracVSeta_mediumPt= ibooker.
bookProfile(
"CHFracVSeta_mediumPt",
"CHFracVSeta_mediumPt",etaBin_, etaMin_, etaMax_,0.,1.2);
704 mNHFracVSeta_mediumPt= ibooker.
bookProfile(
"NHFracVSeta_mediumPt",
"NHFracVSeta_mediumPt",etaBin_, etaMin_, etaMax_,0.,1.2);
705 mPhFracVSeta_mediumPt= ibooker.
bookProfile(
"PhFracVSeta_mediumPt",
"PhFracVSeta_mediumPt",etaBin_, etaMin_, etaMax_,0.,1.2);
706 mElFracVSeta_mediumPt= ibooker.
bookProfile(
"ElFracVSeta_mediumPt",
"ElFracVSeta_mediumPt",etaBin_, etaMin_, etaMax_,0.,1.2);
707 mMuFracVSeta_mediumPt= ibooker.
bookProfile(
"MuFracVSeta_mediumPt",
"MuFracVSeta_mediumPt",etaBin_, etaMin_, etaMax_,0.,1.2);
708 mCHFracVSeta_highPt= ibooker.
bookProfile(
"CHFracVSeta_highPt",
"CHFracVSeta_highPt",etaBin_, etaMin_, etaMax_,0.,1.2);
709 mNHFracVSeta_highPt= ibooker.
bookProfile(
"NHFracVSeta_highPt",
"NHFracVSeta_highPt",etaBin_, etaMin_, etaMax_,0.,1.2);
710 mPhFracVSeta_highPt= ibooker.
bookProfile(
"PhFracVSeta_highPt",
"PhFracVSeta_highPt",etaBin_, etaMin_, etaMax_,0.,1.2);
711 mElFracVSeta_highPt= ibooker.
bookProfile(
"ElFracVSeta_highPt",
"ElFracVSeta_highPt",etaBin_, etaMin_, etaMax_,0.,1.2);
712 mMuFracVSeta_highPt= ibooker.
bookProfile(
"MuFracVSeta_highPt",
"MuFracVSeta_highPt",etaBin_, etaMin_, etaMax_,0.,1.2);
714 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFracVSeta_lowPt" ,mCHFracVSeta_lowPt));
715 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFracVSeta_lowPt" ,mNHFracVSeta_lowPt));
716 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFracVSeta_lowPt" ,mPhFracVSeta_lowPt));
717 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFracVSeta_lowPt" ,mElFracVSeta_lowPt));
718 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuFracVSeta_lowPt" ,mMuFracVSeta_lowPt));
719 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFracVSeta_mediumPt" ,mCHFracVSeta_mediumPt));
720 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFracVSeta_mediumPt" ,mNHFracVSeta_mediumPt));
721 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFracVSeta_mediumPt" ,mPhFracVSeta_mediumPt));
722 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFracVSeta_mediumPt" ,mElFracVSeta_mediumPt));
723 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuFracVSeta_mediumPt" ,mMuFracVSeta_mediumPt));
724 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFracVSeta_highPt" ,mCHFracVSeta_highPt));
725 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFracVSeta_highPt" ,mNHFracVSeta_highPt));
726 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFracVSeta_highPt" ,mPhFracVSeta_highPt));
727 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFracVSeta_highPt" ,mElFracVSeta_highPt));
728 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuFracVSeta_highPt" ,mMuFracVSeta_highPt));
732 mCHFrac_lowPt_Barrel = ibooker.
book1D(
"CHFrac_lowPt_Barrel",
"CHFrac_lowPt_Barrel", 120, -0.1, 1.1);
733 mNHFrac_lowPt_Barrel = ibooker.
book1D(
"NHFrac_lowPt_Barrel",
"NHFrac_lowPt_Barrel", 120, -0.1, 1.1);
734 mPhFrac_lowPt_Barrel = ibooker.
book1D(
"PhFrac_lowPt_Barrel",
"PhFrac_lowPt_Barrel", 120, -0.1, 1.1);
735 mElFrac_lowPt_Barrel = ibooker.
book1D(
"ElFrac_lowPt_Barrel",
"ElFrac_lowPt_Barrel", 120, -0.1, 1.1);
736 mMuFrac_lowPt_Barrel = ibooker.
book1D(
"MuFrac_lowPt_Barrel",
"MuFrac_lowPt_Barrel", 120, -0.1, 1.1);
737 mCHFrac_mediumPt_Barrel = ibooker.
book1D(
"CHFrac_mediumPt_Barrel",
"CHFrac_mediumPt_Barrel", 120, -0.1, 1.1);
738 mNHFrac_mediumPt_Barrel = ibooker.
book1D(
"NHFrac_mediumPt_Barrel",
"NHFrac_mediumPt_Barrel", 120, -0.1, 1.1);
739 mPhFrac_mediumPt_Barrel = ibooker.
book1D(
"PhFrac_mediumPt_Barrel",
"PhFrac_mediumPt_Barrel", 120, -0.1, 1.1);
740 mElFrac_mediumPt_Barrel = ibooker.
book1D(
"ElFrac_mediumPt_Barrel",
"ElFrac_mediumPt_Barrel", 120, -0.1, 1.1);
741 mMuFrac_mediumPt_Barrel = ibooker.
book1D(
"MuFrac_mediumPt_Barrel",
"MuFrac_mediumPt_Barrel", 120, -0.1, 1.1);
742 mCHFrac_highPt_Barrel = ibooker.
book1D(
"CHFrac_highPt_Barrel",
"CHFrac_highPt_Barrel", 120, -0.1, 1.1);
743 mNHFrac_highPt_Barrel = ibooker.
book1D(
"NHFrac_highPt_Barrel",
"NHFrac_highPt_Barrel", 120, -0.1, 1.1);
744 mPhFrac_highPt_Barrel = ibooker.
book1D(
"PhFrac_highPt_Barrel",
"PhFrac_highPt_Barrel", 120, -0.1, 1.1);
745 mElFrac_highPt_Barrel = ibooker.
book1D(
"ElFrac_highPt_Barrel",
"ElFrac_highPt_Barrel", 120, -0.1, 1.1);
746 mMuFrac_highPt_Barrel = ibooker.
book1D(
"MuFrac_highPt_Barrel",
"MuFrac_highPt_Barrel", 120, -0.1, 1.1);
748 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFrac_lowPt_Barrel" ,mCHFrac_lowPt_Barrel));
749 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFrac_lowPt_Barrel" ,mNHFrac_lowPt_Barrel));
750 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFrac_lowPt_Barrel" ,mPhFrac_lowPt_Barrel));
751 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFrac_lowPt_Barrel" ,mElFrac_lowPt_Barrel));
752 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuFrac_lowPt_Barrel" ,mMuFrac_lowPt_Barrel));
753 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFrac_mediumPt_Barrel" ,mCHFrac_mediumPt_Barrel));
754 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFrac_mediumPt_Barrel" ,mNHFrac_mediumPt_Barrel));
755 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFrac_mediumPt_Barrel" ,mPhFrac_mediumPt_Barrel));
756 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFrac_mediumPt_Barrel" ,mElFrac_mediumPt_Barrel));
757 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuFrac_mediumPt_Barrel" ,mMuFrac_mediumPt_Barrel));
758 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFrac_highPt_Barrel" ,mCHFrac_highPt_Barrel));
759 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFrac_highPt_Barrel" ,mNHFrac_highPt_Barrel));
760 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFrac_highPt_Barrel" ,mPhFrac_highPt_Barrel));
761 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFrac_highPt_Barrel" ,mElFrac_highPt_Barrel));
762 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuFrac_highPt_Barrel" ,mMuFrac_highPt_Barrel));
765 mCHEn_lowPt_Barrel = ibooker.
book1D(
"CHEn_lowPt_Barrel",
"CHEn_lowPt_Barrel", ptBin_, 0., ptMax_);
766 mNHEn_lowPt_Barrel = ibooker.
book1D(
"NHEn_lowPt_Barrel",
"NHEn_lowPt_Barrel", ptBin_, 0., ptMax_);
767 mPhEn_lowPt_Barrel = ibooker.
book1D(
"PhEn_lowPt_Barrel",
"PhEn_lowPt_Barrel", ptBin_, 0., ptMax_);
768 mElEn_lowPt_Barrel = ibooker.
book1D(
"ElEn_lowPt_Barrel",
"ElEn_lowPt_Barrel", ptBin_, 0., ptMax_);
769 mMuEn_lowPt_Barrel = ibooker.
book1D(
"MuEn_lowPt_Barrel",
"MuEn_lowPt_Barrel", ptBin_, 0., ptMax_);
770 mCHEn_mediumPt_Barrel = ibooker.
book1D(
"CHEn_mediumPt_Barrel",
"CHEn_mediumPt_Barrel", ptBin_, 0., ptMax_);
771 mNHEn_mediumPt_Barrel = ibooker.
book1D(
"NHEn_mediumPt_Barrel",
"NHEn_mediumPt_Barrel", ptBin_, 0., ptMax_);
772 mPhEn_mediumPt_Barrel = ibooker.
book1D(
"PhEn_mediumPt_Barrel",
"PhEn_mediumPt_Barrel", ptBin_, 0., ptMax_);
773 mElEn_mediumPt_Barrel = ibooker.
book1D(
"ElEn_mediumPt_Barrel",
"ElEn_mediumPt_Barrel", ptBin_, 0., ptMax_);
774 mMuEn_mediumPt_Barrel = ibooker.
book1D(
"MuEn_mediumPt_Barrel",
"MuEn_mediumPt_Barrel", ptBin_, 0., ptMax_);
775 mCHEn_highPt_Barrel = ibooker.
book1D(
"CHEn_highPt_Barrel",
"CHEn_highPt_Barrel", ptBin_, 0., 1.1*ptMax_);
776 mNHEn_highPt_Barrel = ibooker.
book1D(
"NHEn_highPt_Barrel",
"NHEn_highPt_Barrel", ptBin_, 0., ptMax_);
777 mPhEn_highPt_Barrel = ibooker.
book1D(
"PhEn_highPt_Barrel",
"PhEn_highPt_Barrel", ptBin_, 0., ptMax_);
778 mElEn_highPt_Barrel = ibooker.
book1D(
"ElEn_highPt_Barrel",
"ElEn_highPt_Barrel", ptBin_, 0., ptMax_);
779 mMuEn_highPt_Barrel = ibooker.
book1D(
"MuEn_highPt_Barrel",
"MuEn_highPt_Barrel", ptBin_, 0., ptMax_);
781 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHEn_lowPt_Barrel" ,mCHEn_lowPt_Barrel));
782 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHEn_lowPt_Barrel" ,mNHEn_lowPt_Barrel));
783 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhEn_lowPt_Barrel" ,mPhEn_lowPt_Barrel));
784 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElEn_lowPt_Barrel" ,mElEn_lowPt_Barrel));
785 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuEn_lowPt_Barrel" ,mMuEn_lowPt_Barrel));
786 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHEn_mediumPt_Barrel" ,mCHEn_mediumPt_Barrel));
787 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHEn_mediumPt_Barrel" ,mNHEn_mediumPt_Barrel));
788 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhEn_mediumPt_Barrel" ,mPhEn_mediumPt_Barrel));
789 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElEn_mediumPt_Barrel" ,mElEn_mediumPt_Barrel));
790 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuEn_mediumPt_Barrel" ,mMuEn_mediumPt_Barrel));
791 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHEn_highPt_Barrel" ,mCHEn_highPt_Barrel));
792 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHEn_highPt_Barrel" ,mNHEn_highPt_Barrel));
793 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhEn_highPt_Barrel" ,mPhEn_highPt_Barrel));
794 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElEn_highPt_Barrel" ,mElEn_highPt_Barrel));
795 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuEn_highPt_Barrel" ,mMuEn_highPt_Barrel));
798 mChMultiplicity_lowPt_Barrel = ibooker.
book1D(
"ChMultiplicity_lowPt_Barrel",
"ChMultiplicity_lowPt_Barrel", 60,0,60);
799 mNeutMultiplicity_lowPt_Barrel = ibooker.
book1D(
"NeutMultiplicity_lowPt_Barrel",
"NeutMultiplicity_lowPt_Barrel", 60,0,60);
800 mMuMultiplicity_lowPt_Barrel = ibooker.
book1D(
"MuMultiplicity_lowPt_Barrel",
"MuMultiplicity_lowPt_Barrel", 10,0,10);
801 mChMultiplicity_mediumPt_Barrel = ibooker.
book1D(
"ChMultiplicity_mediumPt_Barrel",
"ChMultiplicity_mediumPt_Barrel", 60,0,60);
802 mNeutMultiplicity_mediumPt_Barrel = ibooker.
book1D(
"NeutMultiplicity_mediumPt_Barrel",
"NeutMultiplicity_mediumPt_Barrel", 60,0,60);
803 mMuMultiplicity_mediumPt_Barrel = ibooker.
book1D(
"MuMultiplicity_mediumPt_Barrel",
"MuMultiplicity_mediumPt_Barrel", 10,0,10);
804 mChMultiplicity_highPt_Barrel = ibooker.
book1D(
"ChMultiplicity_highPt_Barrel",
"ChMultiplicity_highPt_Barrel", 60,0,60);
805 mNeutMultiplicity_highPt_Barrel = ibooker.
book1D(
"NeutMultiplicity_highPt_Barrel",
"NeutMultiplicity_highPt_Barrel", 60,0,60);
806 mMuMultiplicity_highPt_Barrel = ibooker.
book1D(
"MuMultiplicity_highPt_Barrel",
"MuMultiplicity_highPt_Barrel", 10,0,10);
808 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChMultiplicity_lowPt_Barrel" ,mChMultiplicity_lowPt_Barrel));
809 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutMultiplicity_lowPt_Barrel" ,mNeutMultiplicity_lowPt_Barrel));
810 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuMultiplicity_lowPt_Barrel" ,mMuMultiplicity_lowPt_Barrel));
811 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChMultiplicity_mediumPt_Barrel" ,mChMultiplicity_mediumPt_Barrel));
812 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutMultiplicity_mediumPt_Barrel" ,mNeutMultiplicity_mediumPt_Barrel));
813 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuMultiplicity_mediumPt_Barrel" ,mMuMultiplicity_mediumPt_Barrel));
814 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChMultiplicity_highPt_Barrel" ,mChMultiplicity_highPt_Barrel));
815 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutMultiplicity_highPt_Barrel" ,mNeutMultiplicity_highPt_Barrel));
816 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuMultiplicity_highPt_Barrel" ,mMuMultiplicity_highPt_Barrel));
819 mCHFracVSpT_Barrel= ibooker.
bookProfile(
"CHFracVSpT_Barrel",
"CHFracVSpT_Barrel",ptBin_, ptMin_, ptMax_,0.,1.2);
820 mNHFracVSpT_Barrel= ibooker.
bookProfile(
"NHFracVSpT_Barrel",
"NHFracVSpT_Barrel",ptBin_, ptMin_, ptMax_,0.,1.2);
821 mPhFracVSpT_Barrel= ibooker.
bookProfile(
"PhFracVSpT_Barrel",
"PhFracVSpT_Barrel",ptBin_, ptMin_, ptMax_,0.,1.2);
822 mElFracVSpT_Barrel= ibooker.
bookProfile(
"ElFracVSpT_Barrel",
"ElFracVSpT_Barrel",ptBin_, ptMin_, ptMax_,0.,1.2);
823 mMuFracVSpT_Barrel= ibooker.
bookProfile(
"MuFracVSpT_Barrel",
"MuFracVSpT_Barrel",ptBin_, ptMin_, ptMax_,0.,1.2);
824 mCHFracVSpT_EndCap= ibooker.
bookProfile(
"CHFracVSpT_EndCap",
"CHFracVSpT_EndCap",ptBin_, ptMin_, ptMax_,0.,1.2);
825 mNHFracVSpT_EndCap= ibooker.
bookProfile(
"NHFracVSpT_EndCap",
"NHFracVSpT_EndCap",ptBin_, ptMin_, ptMax_,0.,1.2);
826 mPhFracVSpT_EndCap= ibooker.
bookProfile(
"PhFracVSpT_EndCap",
"PhFracVSpT_EndCap",ptBin_, ptMin_, ptMax_,0.,1.2);
827 mElFracVSpT_EndCap= ibooker.
bookProfile(
"ElFracVSpT_EndCap",
"ElFracVSpT_EndCap",ptBin_, ptMin_, ptMax_,0.,1.2);
828 mMuFracVSpT_EndCap= ibooker.
bookProfile(
"MuFracVSpT_EndCap",
"MuFracVSpT_EndCap",ptBin_, ptMin_, ptMax_,0.,1.2);
829 mHFHFracVSpT_Forward= ibooker.
bookProfile(
"HFHFracVSpT_Forward",
"HFHFracVSpT_Forward",ptBin_, ptMin_, ptMax_,-0.2,1.2);
830 mHFEFracVSpT_Forward= ibooker.
bookProfile(
"HFEFracVSpT_Forward",
"HFEFracVSpT_Forward",ptBin_, ptMin_, ptMax_,-0.2,1.2);
832 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFracVSpT_Barrel" ,mCHFracVSpT_Barrel));
833 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFracVSpT_Barrel" ,mNHFracVSpT_Barrel));
834 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFracVSpT_Barrel" ,mPhFracVSpT_Barrel));
835 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFracVSpT_Barrel" ,mElFracVSpT_Barrel));
836 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuFracVSpT_Barrel" ,mMuFracVSpT_Barrel));
837 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFracVSpT_EndCap" ,mCHFracVSpT_EndCap));
838 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFracVSpT_EndCap" ,mNHFracVSpT_EndCap));
839 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFracVSpT_EndCap" ,mPhFracVSpT_EndCap));
840 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFracVSpT_EndCap" ,mElFracVSpT_EndCap));
841 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFHFracVSpT_Forward" ,mHFHFracVSpT_Forward));
842 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFEFracVSpT_Forward" ,mHFEFracVSpT_Forward));
846 mCHFrac_lowPt_EndCap = ibooker.
book1D(
"CHFrac_lowPt_EndCap",
"CHFrac_lowPt_EndCap", 120, -0.1, 1.1);
847 mNHFrac_lowPt_EndCap = ibooker.
book1D(
"NHFrac_lowPt_EndCap",
"NHFrac_lowPt_EndCap", 120, -0.1, 1.1);
848 mPhFrac_lowPt_EndCap = ibooker.
book1D(
"PhFrac_lowPt_EndCap",
"PhFrac_lowPt_EndCap", 120, -0.1, 1.1);
849 mElFrac_lowPt_EndCap = ibooker.
book1D(
"ElFrac_lowPt_EndCap",
"ElFrac_lowPt_EndCap", 120, -0.1, 1.1);
850 mMuFrac_lowPt_EndCap = ibooker.
book1D(
"MuFrac_lowPt_EndCap",
"MuFrac_lowPt_EndCap", 120, -0.1, 1.1);
851 mCHFrac_mediumPt_EndCap = ibooker.
book1D(
"CHFrac_mediumPt_EndCap",
"CHFrac_mediumPt_EndCap", 120, -0.1, 1.1);
852 mNHFrac_mediumPt_EndCap = ibooker.
book1D(
"NHFrac_mediumPt_EndCap",
"NHFrac_mediumPt_EndCap", 120, -0.1, 1.1);
853 mPhFrac_mediumPt_EndCap = ibooker.
book1D(
"PhFrac_mediumPt_EndCap",
"PhFrac_mediumPt_EndCap", 120, -0.1, 1.1);
854 mElFrac_mediumPt_EndCap = ibooker.
book1D(
"ElFrac_mediumPt_EndCap",
"ElFrac_mediumPt_EndCap", 120, -0.1, 1.1);
855 mMuFrac_mediumPt_EndCap = ibooker.
book1D(
"MuFrac_mediumPt_EndCap",
"MuFrac_mediumPt_EndCap", 120, -0.1, 1.1);
856 mCHFrac_highPt_EndCap = ibooker.
book1D(
"CHFrac_highPt_EndCap",
"CHFrac_highPt_EndCap", 120, -0.1, 1.1);
857 mNHFrac_highPt_EndCap = ibooker.
book1D(
"NHFrac_highPt_EndCap",
"NHFrac_highPt_EndCap", 120, -0.1, 1.1);
858 mPhFrac_highPt_EndCap = ibooker.
book1D(
"PhFrac_highPt_EndCap",
"PhFrac_highPt_EndCap", 120, -0.1, 1.1);
859 mElFrac_highPt_EndCap = ibooker.
book1D(
"ElFrac_highPt_EndCap",
"ElFrac_highPt_EndCap", 120, -0.1, 1.1);
860 mMuFrac_highPt_EndCap = ibooker.
book1D(
"MuFrac_highPt_EndCap",
"MuFrac_highPt_EndCap", 120, -0.1, 1.1);
862 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFrac_lowPt_EndCap" ,mCHFrac_lowPt_EndCap));
863 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFrac_lowPt_EndCap" ,mNHFrac_lowPt_EndCap));
864 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFrac_lowPt_EndCap" ,mPhFrac_lowPt_EndCap));
865 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFrac_lowPt_EndCap" ,mElFrac_lowPt_EndCap));
866 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuFrac_lowPt_EndCap" ,mMuFrac_lowPt_EndCap));
867 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFrac_mediumPt_EndCap" ,mCHFrac_mediumPt_EndCap));
868 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFrac_mediumPt_EndCap" ,mNHFrac_mediumPt_EndCap));
869 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFrac_mediumPt_EndCap" ,mPhFrac_mediumPt_EndCap));
870 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFrac_mediumPt_EndCap" ,mElFrac_mediumPt_EndCap));
871 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuFrac_mediumPt_EndCap" ,mMuFrac_mediumPt_EndCap));
872 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHFrac_highPt_EndCap" ,mCHFrac_highPt_EndCap));
873 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHFrac_highPt_EndCap" ,mNHFrac_highPt_EndCap));
874 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhFrac_highPt_EndCap" ,mPhFrac_highPt_EndCap));
875 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElFrac_highPt_EndCap" ,mElFrac_highPt_EndCap));
876 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuFrac_highPt_EndCap" ,mMuFrac_highPt_EndCap));
879 mCHEn_lowPt_EndCap = ibooker.
book1D(
"CHEn_lowPt_EndCap",
"CHEn_lowPt_EndCap", ptBin_, 0., ptMax_);
880 mNHEn_lowPt_EndCap = ibooker.
book1D(
"NHEn_lowPt_EndCap",
"NHEn_lowPt_EndCap", ptBin_, 0., ptMax_);
881 mPhEn_lowPt_EndCap = ibooker.
book1D(
"PhEn_lowPt_EndCap",
"PhEn_lowPt_EndCap", ptBin_, 0., ptMax_);
882 mElEn_lowPt_EndCap = ibooker.
book1D(
"ElEn_lowPt_EndCap",
"ElEn_lowPt_EndCap", ptBin_, 0., ptMax_);
883 mMuEn_lowPt_EndCap = ibooker.
book1D(
"MuEn_lowPt_EndCap",
"MuEn_lowPt_EndCap", ptBin_, 0., ptMax_);
884 mCHEn_mediumPt_EndCap = ibooker.
book1D(
"CHEn_mediumPt_EndCap",
"CHEn_mediumPt_EndCap", ptBin_, 0., ptMax_);
885 mNHEn_mediumPt_EndCap = ibooker.
book1D(
"NHEn_mediumPt_EndCap",
"NHEn_mediumPt_EndCap", ptBin_, 0., ptMax_);
886 mPhEn_mediumPt_EndCap = ibooker.
book1D(
"PhEn_mediumPt_EndCap",
"PhEn_mediumPt_EndCap", ptBin_, 0., ptMax_);
887 mElEn_mediumPt_EndCap = ibooker.
book1D(
"ElEn_mediumPt_EndCap",
"ElEn_mediumPt_EndCap", ptBin_, 0., ptMax_);
888 mMuEn_mediumPt_EndCap = ibooker.
book1D(
"MuEn_mediumPt_EndCap",
"MuEn_mediumPt_EndCap", ptBin_, 0., ptMax_);
889 mCHEn_highPt_EndCap = ibooker.
book1D(
"CHEn_highPt_EndCap",
"CHEn_highPt_EndCap", ptBin_, 0., 1.5*ptMax_);
890 mNHEn_highPt_EndCap = ibooker.
book1D(
"NHEn_highPt_EndCap",
"NHEn_highPt_EndCap", ptBin_, 0., 1.5*ptMax_);
891 mPhEn_highPt_EndCap = ibooker.
book1D(
"PhEn_highPt_EndCap",
"PhEn_highPt_EndCap", ptBin_, 0., 1.5*ptMax_);
892 mElEn_highPt_EndCap = ibooker.
book1D(
"ElEn_highPt_EndCap",
"ElEn_highPt_EndCap", ptBin_, 0., ptMax_);
893 mMuEn_highPt_EndCap = ibooker.
book1D(
"MuEn_highPt_EndCap",
"MuEn_highPt_EndCap", ptBin_, 0., ptMax_);
895 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHEn_lowPt_EndCap" ,mCHEn_lowPt_EndCap));
896 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHEn_lowPt_EndCap" ,mNHEn_lowPt_EndCap));
897 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhEn_lowPt_EndCap" ,mPhEn_lowPt_EndCap));
898 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElEn_lowPt_EndCap" ,mElEn_lowPt_EndCap));
899 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuEn_lowPt_EndCap" ,mMuEn_lowPt_EndCap));
900 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHEn_mediumPt_EndCap" ,mCHEn_mediumPt_EndCap));
901 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHEn_mediumPt_EndCap" ,mNHEn_mediumPt_EndCap));
902 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhEn_mediumPt_EndCap" ,mPhEn_mediumPt_EndCap));
903 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElEn_mediumPt_EndCap" ,mElEn_mediumPt_EndCap));
904 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuEn_mediumPt_EndCap" ,mMuEn_mediumPt_EndCap));
905 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"CHEn_highPt_EndCap" ,mCHEn_highPt_EndCap));
906 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NHEn_highPt_EndCap" ,mNHEn_highPt_EndCap));
907 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"PhEn_highPt_EndCap" ,mPhEn_highPt_EndCap));
908 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ElEn_highPt_EndCap" ,mElEn_highPt_EndCap));
909 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuEn_highPt_EndCap" ,mMuEn_highPt_EndCap));
912 mChMultiplicity_lowPt_EndCap = ibooker.
book1D(
"ChMultiplicity_lowPt_EndCap",
"ChMultiplicity_lowPt_EndCap", 60,0,60);
913 mNeutMultiplicity_lowPt_EndCap = ibooker.
book1D(
"NeutMultiplicity_lowPt_EndCap",
"NeutMultiplicity_lowPt_EndCap", 60,0,60);
914 mMuMultiplicity_lowPt_EndCap = ibooker.
book1D(
"MuMultiplicity_lowPt_EndCap",
"MuMultiplicity_lowPt_EndCap", 10,0,10);
915 mChMultiplicity_mediumPt_EndCap = ibooker.
book1D(
"ChMultiplicity_mediumPt_EndCap",
"ChMultiplicity_mediumPt_EndCap", 60,0,60);
916 mNeutMultiplicity_mediumPt_EndCap = ibooker.
book1D(
"NeutMultiplicity_mediumPt_EndCap",
"NeutMultiplicity_mediumPt_EndCap", 60,0,60);
917 mMuMultiplicity_mediumPt_EndCap = ibooker.
book1D(
"MuMultiplicity_mediumPt_EndCap",
"MuMultiplicity_mediumPt_EndCap", 10,0,10);
918 mChMultiplicity_highPt_EndCap = ibooker.
book1D(
"ChMultiplicity_highPt_EndCap",
"ChMultiplicity_highPt_EndCap", 60,0,60);
919 mNeutMultiplicity_highPt_EndCap = ibooker.
book1D(
"NeutMultiplicity_highPt_EndCap",
"NeutMultiplicity_highPt_EndCap", 60,0,60);
920 mMuMultiplicity_highPt_EndCap = ibooker.
book1D(
"MuMultiplicity_highPt_EndCap",
"MuMultiplicity_highPt_EndCap", 10,0,10);
922 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChMultiplicity_lowPt_EndCap" ,mChMultiplicity_lowPt_EndCap));
923 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutMultiplicity_lowPt_EndCap" ,mNeutMultiplicity_lowPt_EndCap));
924 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuMultiplicity_lowPt_EndCap" ,mMuMultiplicity_lowPt_EndCap));
925 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChMultiplicity_mediumPt_EndCap" ,mChMultiplicity_mediumPt_EndCap));
926 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutMultiplicity_mediumPt_EndCap" ,mNeutMultiplicity_mediumPt_EndCap));
927 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuMultiplicity_mediumPt_EndCap" ,mMuMultiplicity_mediumPt_EndCap));
928 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChMultiplicity_highPt_EndCap" ,mChMultiplicity_highPt_EndCap));
929 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutMultiplicity_highPt_EndCap" ,mNeutMultiplicity_highPt_EndCap));
930 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuMultiplicity_highPt_EndCap" ,mMuMultiplicity_highPt_EndCap));
934 mHFEFrac_lowPt_Forward = ibooker.
book1D(
"HFEFrac_lowPt_Forward",
"HFEFrac_lowPt_Forward", 140, -0.2, 1.2);
935 mHFHFrac_lowPt_Forward = ibooker.
book1D(
"HFHFrac_lowPt_Forward",
"HFHFrac_lowPt_Forward", 140, -0.2, 1.2);
936 mHFEFrac_mediumPt_Forward = ibooker.
book1D(
"HFEFrac_mediumPt_Forward",
"HFEFrac_mediumPt_Forward", 140, -0.2, 1.2);
937 mHFHFrac_mediumPt_Forward = ibooker.
book1D(
"HFHFrac_mediumPt_Forward",
"HFHFrac_mediumPt_Forward", 140, -0.2, 1.2);
938 mHFEFrac_highPt_Forward = ibooker.
book1D(
"HFEFrac_highPt_Forward",
"HFEFrac_highPt_Forward", 140, -0.2, 1.2);
939 mHFHFrac_highPt_Forward = ibooker.
book1D(
"HFHFrac_highPt_Forward",
"HFHFrac_highPt_Forward", 140, -0.2, 1.2);
940 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFHFrac_lowPt_Forward" ,mHFHFrac_lowPt_Forward));
941 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFEFrac_lowPt_Forward" ,mHFEFrac_lowPt_Forward));
942 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFHFrac_mediumPt_Forward" ,mHFHFrac_mediumPt_Forward));
943 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFEFrac_mediumPt_Forward" ,mHFEFrac_mediumPt_Forward));
944 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFHFrac_highPt_Forward" ,mHFHFrac_highPt_Forward));
945 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFEFrac_highPt_Forward" ,mHFEFrac_highPt_Forward));
948 mHFEEn_lowPt_Forward = ibooker.
book1D(
"HFEEn_lowPt_Forward",
"HFEEn_lowPt_Forward", ptBin_, 0., ptMax_);
949 mHFHEn_lowPt_Forward = ibooker.
book1D(
"HFHEn_lowPt_Forward",
"HFHEn_lowPt_Forward", ptBin_, 0., 2.0*ptMax_);
950 mHFEEn_mediumPt_Forward = ibooker.
book1D(
"HFEEn_mediumPt_Forward",
"HFEEn_mediumPt_Forward", ptBin_, 0., 1.5*ptMax_);
951 mHFHEn_mediumPt_Forward = ibooker.
book1D(
"HFHEn_mediumPt_Forward",
"HFHEn_mediumPt_Forward", ptBin_, 0., 2.5*ptMax_);
952 mHFEEn_highPt_Forward = ibooker.
book1D(
"HFEEn_highPt_Forward",
"HFEEn_highPt_Forward", ptBin_, 0., 1.5*ptMax_);
953 mHFHEn_highPt_Forward = ibooker.
book1D(
"HFHEn_highPt_Forward",
"HFHEn_highPt_Forward", ptBin_, 0., 5.0*ptMax_);
954 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFHEn_lowPt_Forward" ,mHFHEn_lowPt_Forward));
955 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFEEn_lowPt_Forward" ,mHFEEn_lowPt_Forward));
956 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFHEn_mediumPt_Forward" ,mHFHEn_mediumPt_Forward));
957 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFEEn_mediumPt_Forward" ,mHFEEn_mediumPt_Forward));
958 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFHEn_highPt_Forward" ,mHFHEn_highPt_Forward));
959 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"HFEEn_highPt_Forward" ,mHFEEn_highPt_Forward));
961 mChMultiplicity_lowPt_Forward = ibooker.
book1D(
"ChMultiplicity_lowPt_Forward",
"ChMultiplicity_lowPt_Forward", 60,0,60);
962 mNeutMultiplicity_lowPt_Forward = ibooker.
book1D(
"NeutMultiplicity_lowPt_Forward",
"NeutMultiplicity_lowPt_Forward", 60,0,60);
963 mChMultiplicity_mediumPt_Forward = ibooker.
book1D(
"ChMultiplicity_mediumPt_Forward",
"ChMultiplicity_mediumPt_Forward", 60,0,60);
964 mNeutMultiplicity_mediumPt_Forward = ibooker.
book1D(
"NeutMultiplicity_mediumPt_Forward",
"NeutMultiplicity_mediumPt_Forward", 60,0,60);
965 mChMultiplicity_highPt_Forward = ibooker.
book1D(
"ChMultiplicity_highPt_Forward",
"ChMultiplicity_highPt_Forward", 60,0,60);
966 mNeutMultiplicity_highPt_Forward = ibooker.
book1D(
"NeutMultiplicity_highPt_Forward",
"NeutMultiplicity_highPt_Forward", 60,0,60);
968 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChMultiplicity_lowPt_Forward" ,mChMultiplicity_lowPt_Forward));
969 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutMultiplicity_lowPt_Forward" ,mNeutMultiplicity_lowPt_Forward));
970 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChMultiplicity_mediumPt_Forward" ,mChMultiplicity_mediumPt_Forward));
971 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutMultiplicity_mediumPt_Forward" ,mNeutMultiplicity_mediumPt_Forward));
972 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChMultiplicity_highPt_Forward" ,mChMultiplicity_highPt_Forward));
973 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutMultiplicity_highPt_Forward" ,mNeutMultiplicity_highPt_Forward));
975 mChargedHadronEnergy = ibooker.
book1D(
"ChargedHadronEnergy",
"charged HAD energy", 100, 0, 100);
976 mNeutralHadronEnergy = ibooker.
book1D(
"NeutralHadronEnergy",
"neutral HAD energy", 100, 0, 100);
977 mChargedEmEnergy = ibooker.
book1D(
"ChargedEmEnergy",
"charged EM energy ", 100, 0, 100);
978 mChargedMuEnergy = ibooker.
book1D(
"ChargedMuEnergy",
"charged Mu energy", 100, 0, 100);
979 mNeutralEmEnergy = ibooker.
book1D(
"NeutralEmEnergy",
"neutral EM energy", 100, 0, 100);
980 mChargedMultiplicity = ibooker.
book1D(
"ChargedMultiplicity",
"charged multiplicity ", 100, 0, 100);
981 mNeutralMultiplicity = ibooker.
book1D(
"NeutralMultiplicity",
"neutral multiplicity", 100, 0, 100);
982 mMuonMultiplicity = ibooker.
book1D(
"MuonMultiplicity",
"muon multiplicity", 100, 0, 100);
983 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChargedHadronEnergy" ,mChargedHadronEnergy));
984 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutralHadronEnergy" ,mNeutralHadronEnergy));
985 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChargedEmEnergy" ,mChargedEmEnergy));
986 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChargedMuEnergy" ,mChargedMuEnergy));
987 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutralEmEnergy" ,mNeutralEmEnergy));
988 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChargedMultiplicity" ,mChargedMultiplicity));
989 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutralMultiplicity" ,mNeutralMultiplicity));
990 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuonMultiplicity" ,mMuonMultiplicity));
994 mChargedHadronEnergy_profile = ibooker.
bookProfile(
"ChargedHadronEnergy_profile",
"charged HAD energy", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
995 mNeutralHadronEnergy_profile = ibooker.
bookProfile(
"NeutralHadronEnergy_profile",
"neutral HAD energy", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
996 mChargedEmEnergy_profile = ibooker.
bookProfile(
"ChargedEmEnergy_profile",
"charged EM energy", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
997 mChargedMuEnergy_profile = ibooker.
bookProfile(
"ChargedMuEnergy_profile",
"charged Mu energy", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
998 mNeutralEmEnergy_profile = ibooker.
bookProfile(
"NeutralEmEnergy_profile",
"neutral EM energy", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
999 mChargedMultiplicity_profile = ibooker.
bookProfile(
"ChargedMultiplicity_profile",
"charged multiplicity", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
1000 mNeutralMultiplicity_profile = ibooker.
bookProfile(
"NeutralMultiplicity_profile",
"neutral multiplicity", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
1001 mMuonMultiplicity_profile = ibooker.
bookProfile(
"MuonMultiplicity_profile",
"muon multiplicity", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
1006 mNeutralHadronEnergy_profile->setAxisTitle(
"nvtx",1);
1007 mChargedEmEnergy_profile ->setAxisTitle(
"nvtx",1);
1008 mChargedMuEnergy_profile ->setAxisTitle(
"nvtx",1);
1009 mNeutralEmEnergy_profile ->setAxisTitle(
"nvtx",1);
1010 mChargedMultiplicity_profile->setAxisTitle(
"nvtx",1);
1011 mNeutralMultiplicity_profile->setAxisTitle(
"nvtx",1);
1012 mMuonMultiplicity_profile ->setAxisTitle(
"nvtx",1);
1014 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChargedHadronEnergy_profile" ,mChargedHadronEnergy_profile));
1015 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutralHadronEnergy_profile" ,mNeutralHadronEnergy_profile));
1016 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChargedEmEnergy_profile" ,mChargedEmEnergy_profile));
1017 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChargedMuEnergy_profile" ,mChargedMuEnergy_profile));
1018 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutralEmEnergy_profile" ,mNeutralEmEnergy_profile));
1019 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"ChargedMultiplicity_profile" ,mChargedMultiplicity_profile));
1020 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutralMultiplicity_profile" ,mNeutralMultiplicity_profile));
1021 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"MuonMultiplicity_profile" ,mMuonMultiplicity_profile));
1023 mNeutralFraction = ibooker.
book1D(
"NeutralConstituentsFraction",
"Neutral Constituents Fraction",100,0,1);
1024 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
DirName+
"/"+
"NeutralConstituentsFraction" ,mNeutralFraction));
1028 if(jetCleaningFlag_){
1030 folderNames_.push_back(
"DiJet");
1032 for (std::vector<std::string>::const_iterator ic = folderNames_.begin();
1033 ic != folderNames_.end(); ic++){
1034 bookMESetSelection(
DirName+
"/"+*ic, ibooker);
1039 cleanupME = ibooker.
book1D(
"cleanup",
"cleanup", 10, 0., 10.);
1041 cleanupME->setBinLabel(2,
"DCS::Pixel");
1042 cleanupME->setBinLabel(3,
"DCS::SiStrip");
1043 cleanupME->setBinLabel(4,
"DCS::ECAL");
1044 cleanupME->setBinLabel(5,
"DCS::ES");
1045 cleanupME->setBinLabel(6,
"DCS::HBHE");
1046 cleanupME->setBinLabel(7,
"DCS::HF");
1047 cleanupME->setBinLabel(8,
"DCS::HO");
1048 cleanupME->setBinLabel(9,
"DCS::Muon");
1049 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
"JetMET/cleanup" ,cleanupME));
1051 verticesME = ibooker.
book1D(
"vertices",
"vertices", 100, 0, 100);
1052 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(
"JetMET/vertices" ,verticesME));
1062 mPt = ibooker.
book1D(
"Pt",
"pt", ptBin_, ptMin_, ptMax_);
1063 mEta = ibooker.
book1D(
"Eta",
"eta", etaBin_, etaMin_, etaMax_);
1064 mPhi = ibooker.
book1D(
"Phi",
"phi", phiBin_, phiMin_, phiMax_);
1065 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"Pt" ,mPt));
1066 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"Eta" ,mEta));
1067 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"Phi" ,mPhi));
1069 mConstituents = ibooker.
book1D(
"Constituents",
"# of constituents", 50, 0, 100);
1070 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"Constituents" ,mConstituents));
1072 mJetEnergyCorr= ibooker.
book1D(
"JetEnergyCorr",
"jet energy correction factor", 50, 0.0,3.0);
1073 mJetEnergyCorrVSEta= ibooker.
bookProfile(
"JetEnergyCorrVSEta",
"jet energy correction factor VS eta", etaBin_, etaMin_,etaMax_, 0.0,3.0);
1074 mJetEnergyCorrVSPt= ibooker.
bookProfile(
"JetEnergyCorrVSPt",
"jet energy correction factor VS pt", ptBin_, ptMin_,ptMax_, 0.0,3.0);
1075 mHFrac = ibooker.
book1D(
"HFrac",
"HFrac", 140, -0.2, 1.2);
1076 mEFrac = ibooker.
book1D(
"EFrac",
"EFrac", 140, -0.2, 1.2);
1077 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"JetEnergyCorr" ,mJetEnergyCorr));
1078 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"JetEnergyCorrVSEta" ,mJetEnergyCorrVSEta));
1079 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"JetEnergyCorrVSPt" ,mJetEnergyCorrVSPt));
1080 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"HFrac" ,mHFrac));
1081 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"EFrac" ,mEFrac));
1083 mDPhi = ibooker.
book1D(
"DPhi",
"dPhi btw the two leading jets", 100, 0., acos(-1.));
1084 mDijetAsymmetry = ibooker.
book1D(
"DijetAsymmetry",
"DijetAsymmetry", 100, -1., 1.);
1085 mDijetBalance = ibooker.
book1D(
"DijetBalance",
"DijetBalance", 100, -2., 2.);
1086 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"DPhi" ,mDPhi));
1087 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"DijetAsymmetry" ,mDijetAsymmetry));
1088 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"DijetBalance" ,mDijetBalance));
1092 mPt_profile = ibooker.
bookProfile(
"Pt_profile",
"pt", nbinsPV_, nPVlow_, nPVhigh_, ptBin_, ptMin_, ptMax_);
1093 mEta_profile = ibooker.
bookProfile(
"Eta_profile",
"eta", nbinsPV_, nPVlow_, nPVhigh_, etaBin_, etaMin_, etaMax_);
1094 mPhi_profile = ibooker.
bookProfile(
"Phi_profile",
"phi", nbinsPV_, nPVlow_, nPVhigh_, phiBin_, phiMin_, phiMax_);
1096 mConstituents_profile = ibooker.
bookProfile(
"Constituents_profile",
"# of constituents", nbinsPV_, nPVlow_, nPVhigh_, 50, 0, 100);
1098 mHFrac_profile = ibooker.
bookProfile(
"HFrac_profile",
"HFrac", nbinsPV_, nPVlow_, nPVhigh_, 140, -0.2, 1.2);
1099 mEFrac_profile = ibooker.
bookProfile(
"EFrac_profile",
"EFrac", nbinsPV_, nPVlow_, nPVhigh_, 140, -0.2, 1.2);
1103 mEta_profile ->setAxisTitle(
"nvtx",1);
1104 mPhi_profile ->setAxisTitle(
"nvtx",1);
1106 mConstituents_profile->setAxisTitle(
"nvtx",1);
1107 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"Constituents_profile",mConstituents_profile));
1109 mHFrac_profile ->setAxisTitle(
"nvtx",1);
1110 mEFrac_profile ->setAxisTitle(
"nvtx",1);
1112 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"Pt_profile" ,mPt_profile));
1113 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"Eta_profile",mEta_profile));
1114 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"Phi_profile",mPhi_profile));
1115 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"HFrac_profile",mHFrac_profile));
1116 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"EFrac_profile",mEFrac_profile));
1122 mMaxEInEmTowers = ibooker.
book1D(
"MaxEInEmTowers",
"MaxEInEmTowers", 150, 0, 150);
1123 mMaxEInHadTowers = ibooker.
book1D(
"MaxEInHadTowers",
"MaxEInHadTowers", 150, 0, 150);
1125 mresEMF = ibooker.
book1D(
"resEMF",
"resEMF", 50, 0., 1.);
1126 mN90Hits = ibooker.
book1D(
"N90Hits",
"N90Hits", 100, 0., 100);
1127 mfHPD = ibooker.
book1D(
"fHPD",
"fHPD", 50, 0., 1.);
1128 mfRBX = ibooker.
book1D(
"fRBX",
"fRBX", 50, 0., 1.);
1130 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MaxEInEmTowers" ,mMaxEInEmTowers));
1131 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MaxEInHadTowers" ,mMaxEInHadTowers));
1132 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"resEMF" ,mresEMF));
1133 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"N90Hits",mN90Hits));
1134 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"fHPD" ,mfHPD));
1135 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"fRBX" ,mfRBX));
1142 mCHFrac = ibooker.
book1D(
"CHFrac",
"CHFrac", 120, -0.1, 1.1);
1143 mNHFrac = ibooker.
book1D(
"NHFrac",
"NHFrac", 120, -0.1, 1.1);
1144 mPhFrac = ibooker.
book1D(
"PhFrac",
"PhFrac", 120, -0.1, 1.1);
1145 mElFrac = ibooker.
book1D(
"ElFrac",
"ElFrac", 120, -0.1, 1.1);
1146 mMuFrac = ibooker.
book1D(
"MuFrac",
"MuFrac", 120, -0.1, 1.1);
1147 mHFEMFrac = ibooker.
book1D(
"HFEMFrac",
"HFEMFrac", 120, -0.1, 1.1);
1148 mHFHFrac = ibooker.
book1D(
"HFHFrac",
"HFHFrac", 120, -0.1, 1.1);
1150 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CHFrac" ,mCHFrac));
1151 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"NHFrac" ,mNHFrac));
1152 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PhFrac" ,mPhFrac));
1153 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"ElFrac" ,mElFrac));
1154 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MuFrac" ,mMuFrac));
1155 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"HFEMFrac",mHFEMFrac));
1156 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"HFHFrac" ,mHFHFrac));
1158 mChargedMultiplicity = ibooker.
book1D(
"ChargedMultiplicity",
"charged multiplicity ", 100, 0, 100);
1159 mNeutralMultiplicity = ibooker.
book1D(
"NeutralMultiplicity",
"neutral multiplicity", 100, 0, 100);
1160 mMuonMultiplicity = ibooker.
book1D(
"MuonMultiplicity",
"muon multiplicity", 100, 0, 100);
1162 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"ChargedMultiplicity" ,mChargedMultiplicity));
1163 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"NeutralMultiplicity" ,mNeutralMultiplicity));
1164 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MuonMultiplicity" ,mMuonMultiplicity));
1168 mCHFrac_profile = ibooker.
bookProfile(
"CHFrac_profile",
"charged HAD fraction profile", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 1);
1169 mNHFrac_profile = ibooker.
bookProfile(
"NHFrac_profile",
"neutral HAD fraction profile", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 1);
1170 mElFrac_profile = ibooker.
bookProfile(
"ElFrac_profile",
"Electron Fraction Profile", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 1);
1171 mMuFrac_profile = ibooker.
bookProfile(
"MuFrac_profile",
"Muon Fraction Profile", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 1);
1172 mPhFrac_profile = ibooker.
bookProfile(
"PhFrac_profile",
"Photon Fraction Profile", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 1);
1173 mHFEMFrac_profile = ibooker.
bookProfile(
"HFEMFrac_profile",
"HF electomagnetic fraction Profile", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 1);
1174 mHFHFrac_profile = ibooker.
bookProfile(
"HFHFrac_profile",
"HF hadronic fraction profile", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 1);
1175 mChargedMultiplicity_profile = ibooker.
bookProfile(
"ChargedMultiplicity_profile",
"charged multiplicity", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
1176 mNeutralMultiplicity_profile = ibooker.
bookProfile(
"NeutralMultiplicity_profile",
"neutral multiplicity", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
1177 mMuonMultiplicity_profile = ibooker.
bookProfile(
"MuonMultiplicity_profile",
"muon multiplicity", nbinsPV_, nPVlow_, nPVhigh_, 100, 0, 100);
1182 mNHFrac_profile ->setAxisTitle(
"nvtx",1);
1183 mElFrac_profile ->setAxisTitle(
"nvtx",1);
1184 mMuFrac_profile ->setAxisTitle(
"nvtx",1);
1185 mPhFrac_profile ->setAxisTitle(
"nvtx",1);
1186 mHFEMFrac_profile ->setAxisTitle(
"nvtx",1);
1187 mHFHFrac_profile ->setAxisTitle(
"nvtx",1);
1188 mChargedMultiplicity_profile->setAxisTitle(
"nvtx",1);
1189 mNeutralMultiplicity_profile->setAxisTitle(
"nvtx",1);
1190 mMuonMultiplicity_profile ->setAxisTitle(
"nvtx",1);
1192 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"CHFrac_profile" ,mCHFrac_profile));
1193 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"NHFrac_profile" ,mNHFrac_profile));
1194 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"PhFrac_profile" ,mPhFrac_profile));
1195 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"ElFrac_profile" ,mElFrac_profile));
1196 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MuFrac_profile" ,mMuFrac_profile));
1197 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"HFEMFrac_profile",mHFEMFrac_profile));
1198 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"HFHFrac_profile" ,mHFHFrac_profile));
1199 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"ChargedMultiplicity_profile" ,mChargedMultiplicity_profile));
1200 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"NeutralMultiplicity_profile" ,mNeutralMultiplicity_profile));
1201 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"MuonMultiplicity_profile" ,mMuonMultiplicity_profile));
1203 mNeutralFraction = ibooker.
book1D(
"NeutralConstituentsFraction",
"Neutral Constituents Fraction",100,0,1);
1204 map_of_MEs.insert(std::pair<std::string,MonitorElement*>(DirName+
"/"+
"NeutralConstituentsFraction" ,mNeutralFraction));
1215 if ( highPtJetEventFlag_->on() ) highPtJetEventFlag_->initRun( iRun, iSetup );
1216 if ( lowPtJetEventFlag_ ->on() ) lowPtJetEventFlag_ ->initRun( iRun, iSetup );
1218 if (highPtJetEventFlag_->on() && highPtJetEventFlag_->expressionsFromDB(highPtJetEventFlag_->hltDBKey(), iSetup)[0] !=
"CONFIG_ERROR")
1219 highPtJetExpr_ = highPtJetEventFlag_->expressionsFromDB(highPtJetEventFlag_->hltDBKey(), iSetup);
1220 if (lowPtJetEventFlag_->on() && lowPtJetEventFlag_->expressionsFromDB(lowPtJetEventFlag_->hltDBKey(), iSetup)[0] !=
"CONFIG_ERROR")
1221 lowPtJetExpr_ = lowPtJetEventFlag_->expressionsFromDB(lowPtJetEventFlag_->hltDBKey(), iSetup);
1228 hltInitialized_ = hltConfig_.init(iRun,iSetup,processname_,changed);
1229 if (!hltInitialized_) {
1231 processname_ =
"FU";
1232 hltInitialized_ = hltConfig_.init(iRun,iSetup,processname_,changed);
1233 if(!hltInitialized_){
1235 LogDebug(
"JetAnalyzer") <<
"HLTConfigProvider failed to initialize.";
1251 if(jetCleaningFlag_){
1253 DirName =
"JetMET/Jet/Cleaned"+mInputCollection_.label();
1256 DirName =
"JetMET/Jet/Uncleaned"+mInputCollection_.label();
1264 Int_t JetLoPass = 0;
1265 Int_t JetHiPass = 0;
1267 if (triggerResults.
isValid()){
1270 const unsigned int nTrig(triggerNames.
size());
1271 for (
unsigned int i=0;
i<nTrig;++
i)
1273 if (triggerNames.
triggerName(
i).find(highPtJetExpr_[0].substr(0,highPtJetExpr_[0].rfind(
"_v")+2))!=std::string::npos && triggerResults->accept(
i))
1275 else if (triggerNames.
triggerName(
i).find(lowPtJetExpr_[0].substr(0,lowPtJetExpr_[0].rfind(
"_v")+2))!=std::string::npos && triggerResults->accept(
i))
1281 if (verbose_)
std::cout <<
"trigger label " << theTriggerResultsLabel_ << std::endl;
1285 std::cout <<
">>> Trigger Lo = " << JetLoPass
1286 <<
" Hi = " << JetHiPass
1293 iEvent.
getByToken(vertexToken_, vertexHandle);
1295 if (!vertexHandle.
isValid()) {
1296 LogDebug(
"") <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1297 if (verbose_)
std::cout <<
"CaloMETAnalyzer: Could not find vertex collection" << std::endl;
1300 if ( vertexHandle.
isValid() ){
1302 numPV = vertexCollection.size();
1304 bool bPrimaryVertex = (bypassAllPVChecks_ || (numPV>0));
1305 if(fill_jet_high_level_histo){
1306 verticesME = map_of_MEs[
"JetMET/vertices"];
if(verticesME && verticesME->getRootObject())verticesME->Fill(numPV);
1311 iEvent.
getByToken(gtToken_, gtReadoutRecord);
1314 if (!gtReadoutRecord.
isValid()) {
1315 LogInfo(
"JetAnalyzer") <<
"JetAnalyzer: Could not find GT readout record" << std::endl;
1316 if (verbose_)
std::cout <<
"JetAnalyzer: Could not find GT readout record product" << std::endl;
1318 DCSFilterForDCSMonitoring_->filter(iEvent, iSetup);
1319 if(fill_jet_high_level_histo){
1320 cleanupME = map_of_MEs[
"JetMET/cleanup"];
if(cleanupME && cleanupME->getRootObject()){
1321 if (bPrimaryVertex) cleanupME->Fill(0.5);
1322 if ( DCSFilterForDCSMonitoring_->passPIX ) cleanupME->Fill(1.5);
1323 if ( DCSFilterForDCSMonitoring_->passSiStrip ) cleanupME->Fill(2.5);
1324 if ( DCSFilterForDCSMonitoring_->passECAL ) cleanupME->Fill(3.5);
1325 if ( DCSFilterForDCSMonitoring_->passES ) cleanupME->Fill(4.5);
1326 if ( DCSFilterForDCSMonitoring_->passHBHE ) cleanupME->Fill(5.5);
1327 if ( DCSFilterForDCSMonitoring_->passHF ) cleanupME->Fill(6.5);
1328 if ( DCSFilterForDCSMonitoring_->passHO ) cleanupME->Fill(7.5);
1329 if ( DCSFilterForDCSMonitoring_->passMuon ) cleanupME->Fill(8.5);
1336 if (isCaloJet_) iEvent.
getByToken(caloJetsToken_, caloJets);
1338 if (isPFJet_) iEvent.
getByToken(pfJetsToken_, pfJets);
1343 iEvent.
getByToken(jetID_ValueMapToken_,jetID_ValueMap_Handle);
1348 bool dcsDecision = DCSFilterForJetMonitoring_->filter(iEvent, iSetup);
1349 bool jetCollectionIsValid =
false;
1350 if (isCaloJet_) jetCollectionIsValid = caloJets.
isValid();
1352 if (isPFJet_) jetCollectionIsValid = pfJets.
isValid();
1355 if (jetCleaningFlag_ && (!jetCollectionIsValid || !bPrimaryVertex || !dcsDecision))
return;
1357 unsigned int collSize=-1;
1358 if (isCaloJet_) collSize = caloJets->size();
1367 if (isPFJet_) collSize=pfJets->size();
1371 std::vector<Jet> corJets;
1375 unsigned int ind1=-1;
1377 unsigned int ind2=-1;
1379 unsigned int ind3=-1;
1382 bool cleaned_first_jet=
false;
1383 bool cleaned_second_jet=
false;
1384 bool cleaned_third_jet=
false;
1393 for (
unsigned int ijet=0; ijet<collSize; ijet++) {
1396 bool pass_uncorrected=
false;
1397 bool pass_corrected=
false;
1399 correctedJet=(*caloJets)[ijet];
1405 correctedJet=(*pfJets)[ijet];
1407 if(correctedJet.
pt()>ptThresholdUnc_){
1408 pass_uncorrected=
true;
1410 if (!jetCorrectionService_.empty()) {
1415 scale = corrector->
correction((*caloJets)[ijet], iEvent, iSetup);
1421 scale = corrector->
correction((*pfJets)[ijet], iEvent, iSetup);
1426 if(correctedJet.
pt()> ptThreshold_){
1427 pass_corrected=
true;
1430 if (!pass_corrected && !pass_uncorrected)
continue;
1433 recoJets.push_back(correctedJet);
1436 bool jetpassid=
true;
1437 bool Thiscleaned=
true;
1443 jetpassid = jetIDFunctor((*caloJets)[ijet], jetID);
1444 if(jetCleaningFlag_){
1445 Thiscleaned=jetpassid;
1447 if(Thiscleaned && pass_corrected){
1448 mN90Hits = map_of_MEs[
DirName+
"/"+
"N90Hits"];
if (mN90Hits && mN90Hits->getRootObject()) mN90Hits->Fill (jetID.
n90Hits);
1449 mfHPD = map_of_MEs[
DirName+
"/"+
"fHPD"];
if (mfHPD && mfHPD->getRootObject()) mfHPD->Fill (jetID.
fHPD);
1450 mresEMF = map_of_MEs[
DirName+
"/"+
"resEMF"];
if (mresEMF && mresEMF->getRootObject()) mresEMF->Fill (jetID.
restrictedEMF);
1451 mfRBX = map_of_MEs[
DirName+
"/"+
"fRBX"];
if (mfRBX && mfRBX->getRootObject()) mfRBX->Fill (jetID.
fRBX);
1454 if(jetCleaningFlag_){
1455 Thiscleaned=jetpassid;
1457 if(Thiscleaned && pass_uncorrected){
1458 mPt_uncor = map_of_MEs[
DirName+
"/"+
"Pt_uncor"];
if (mPt_uncor && mPt_uncor->getRootObject()) mPt_uncor->Fill ((*caloJets)[ijet].pt());
1459 mEta_uncor = map_of_MEs[
DirName+
"/"+
"Eta_uncor"];
if (mEta_uncor && mEta_uncor->getRootObject()) mEta_uncor->Fill ((*caloJets)[ijet].eta());
1460 mPhi_uncor = map_of_MEs[
DirName+
"/"+
"Phi_uncor"];
if (mPhi_uncor && mPhi_uncor->getRootObject()) mPhi_uncor->Fill ((*caloJets)[ijet].phi());
1462 mConstituents_uncor = map_of_MEs[
DirName+
"/"+
"Constituents_uncor"];
if (mConstituents_uncor && mConstituents_uncor->getRootObject()) mConstituents_uncor->Fill ((*caloJets)[ijet].nConstituents());
1466 if(Thiscleaned && pass_corrected){
1467 mHFrac = map_of_MEs[
DirName+
"/"+
"HFrac"];
if (mHFrac && mHFrac->getRootObject()) mHFrac->Fill ((*caloJets)[ijet].energyFractionHadronic());
1468 mEFrac = map_of_MEs[
DirName+
"/"+
"EFrac"];
if (mEFrac && mHFrac->getRootObject()) mEFrac->Fill ((*caloJets)[ijet].emEnergyFraction());
1469 mHFrac_profile = map_of_MEs[
DirName+
"/"+
"HFrac_profile"];
if (mHFrac_profile && mHFrac_profile->getRootObject()) mHFrac_profile ->
Fill(numPV, (*caloJets)[ijet].energyFractionHadronic());
1470 mEFrac_profile = map_of_MEs[
DirName+
"/"+
"EFrac_profile"];
if (mEFrac_profile && mEFrac_profile->getRootObject()) mEFrac_profile ->
Fill(numPV, (*caloJets)[ijet].emEnergyFraction());
1471 if (fabs((*caloJets)[ijet].
eta()) <= 1.3) {
1472 mHFrac_Barrel = map_of_MEs[
DirName+
"/"+
"HFrac_Barrel"];
if (mHFrac_Barrel && mHFrac_Barrel->getRootObject()) mHFrac_Barrel->Fill((*caloJets)[ijet].energyFractionHadronic());
1473 mEFrac_Barrel = map_of_MEs[
DirName+
"/"+
"EFrac_Barrel"];
if (mEFrac_Barrel && mEFrac_Barrel->getRootObject()) mEFrac_Barrel->Fill((*caloJets)[ijet].emEnergyFraction());
1474 }
else if(fabs((*caloJets)[ijet].
eta()) <3.0){
1475 mHFrac_EndCap = map_of_MEs[
DirName+
"/"+
"HFrac_EndCap"];
if (mHFrac_EndCap && mHFrac_EndCap->getRootObject()) mHFrac_EndCap->Fill((*caloJets)[ijet].energyFractionHadronic());
1476 mEFrac_EndCap = map_of_MEs[
DirName+
"/"+
"EFrac_EndCap"];
if (mEFrac_EndCap && mEFrac_EndCap->getRootObject()) mEFrac_EndCap->Fill((*caloJets)[ijet].emEnergyFraction());
1478 mHFrac_Forward = map_of_MEs[
DirName+
"/"+
"HFrac_Forward"];
if (mHFrac_Forward && mHFrac_Forward->getRootObject()) mHFrac_Forward->Fill((*caloJets)[ijet].energyFractionHadronic());
1479 mEFrac_Forward = map_of_MEs[
DirName+
"/"+
"EFrac_Forward"];
if (mEFrac_Forward && mEFrac_Forward->getRootObject()) mEFrac_Forward->Fill((*caloJets)[ijet].emEnergyFraction());
1481 mMaxEInEmTowers = map_of_MEs[
DirName+
"/"+
"MaxEInEmTowers"];
if (mMaxEInEmTowers && mMaxEInEmTowers->getRootObject()) mMaxEInEmTowers->Fill ((*caloJets)[ijet].maxEInEmTowers());
1482 mMaxEInHadTowers = map_of_MEs[
DirName+
"/"+
"MaxEInHadTowers"];
if (mMaxEInHadTowers && mMaxEInHadTowers->getRootObject()) mMaxEInHadTowers->Fill ((*caloJets)[ijet].maxEInHadTowers());
1484 mHadEnergyInHO = map_of_MEs[
DirName+
"/"+
"HadEnergyInHO"];
if (mHadEnergyInHO && mHadEnergyInHO->getRootObject()) mHadEnergyInHO->Fill ((*caloJets)[ijet].hadEnergyInHO());
1485 mHadEnergyInHB = map_of_MEs[
DirName+
"/"+
"HadEnergyInHB"];
if (mHadEnergyInHB && mHadEnergyInHB->getRootObject()) mHadEnergyInHB->Fill ((*caloJets)[ijet].hadEnergyInHB());
1486 mHadEnergyInHF = map_of_MEs[
DirName+
"/"+
"HadEnergyInHF"];
if (mHadEnergyInHF && mHadEnergyInHF->getRootObject()) mHadEnergyInHF->Fill ((*caloJets)[ijet].hadEnergyInHF());
1487 mHadEnergyInHE = map_of_MEs[
DirName+
"/"+
"HadEnergyInHE"];
if (mHadEnergyInHE && mHadEnergyInHE->getRootObject()) mHadEnergyInHE->Fill ((*caloJets)[ijet].hadEnergyInHE());
1488 mEmEnergyInEB = map_of_MEs[
DirName+
"/"+
"EmEnergyInEB"];
if (mEmEnergyInEB && mEmEnergyInEB->getRootObject()) mEmEnergyInEB->Fill ((*caloJets)[ijet].emEnergyInEB());
1489 mEmEnergyInEE = map_of_MEs[
DirName+
"/"+
"EmEnergyInEE"];
if (mEmEnergyInEE && mEmEnergyInEE->getRootObject()) mEmEnergyInEE->Fill ((*caloJets)[ijet].emEnergyInEE());
1490 mEmEnergyInHF = map_of_MEs[
DirName+
"/"+
"EmEnergyInHF"];
if (mEmEnergyInHF && mEmEnergyInHF->getRootObject()) mEmEnergyInHF->Fill ((*caloJets)[ijet].emEnergyInHF());
1718 jetpassid = pfjetIDFunctor((*pfJets)[ijet]);
1719 if(jetCleaningFlag_){
1720 Thiscleaned = jetpassid;
1722 if(Thiscleaned && pass_uncorrected){
1723 mPt_uncor = map_of_MEs[
DirName+
"/"+
"Pt_uncor"];
if (mPt_uncor && mPt_uncor->getRootObject())
if (mPt_uncor) mPt_uncor->Fill ((*pfJets)[ijet].
pt());
1724 mEta_uncor = map_of_MEs[
DirName+
"/"+
"Eta_uncor"];
if (mEta_uncor && mEta_uncor->getRootObject())
if (mEta_uncor) mEta_uncor->Fill ((*pfJets)[ijet].
eta());
1725 mPhi_uncor = map_of_MEs[
DirName+
"/"+
"Phi_uncor"];
if (mPhi_uncor && mPhi_uncor->getRootObject())
if (mPhi_uncor) mPhi_uncor->Fill ((*pfJets)[ijet].
phi());
1727 mConstituents_uncor = map_of_MEs[
DirName+
"/"+
"Constituents_uncor"];
if (mConstituents_uncor && mConstituents_uncor->getRootObject())
if (mConstituents_uncor) mConstituents_uncor->Fill ((*pfJets)[ijet].nConstituents());
1730 if(Thiscleaned && pass_corrected){
1731 mHFrac = map_of_MEs[
DirName+
"/"+
"HFrac"];
if (mHFrac && mHFrac->getRootObject()) mHFrac->Fill ((*pfJets)[ijet].chargedHadronEnergyFraction()+(*pfJets)[ijet].neutralHadronEnergyFraction()+(*pfJets)[ijet].HFHadronEnergyFraction ());
1732 mEFrac = map_of_MEs[
DirName+
"/"+
"EFrac"];
if (mEFrac && mHFrac->getRootObject()) mEFrac->Fill ((*pfJets)[ijet].chargedEmEnergyFraction() +(*pfJets)[ijet].neutralEmEnergyFraction()+(*pfJets)[ijet].HFEMEnergyFraction ());
1733 if ((*pfJets)[ijet].pt()<= 50) {
1734 mCHFracVSeta_lowPt = map_of_MEs[
DirName+
"/"+
"CHFracVSeta_lowPt"];
if (mCHFracVSeta_lowPt && mCHFracVSeta_lowPt->getRootObject()) mCHFracVSeta_lowPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].chargedHadronEnergyFraction());
1735 mNHFracVSeta_lowPt = map_of_MEs[
DirName+
"/"+
"NHFracVSeta_lowPt"];
if (mNHFracVSeta_lowPt && mNHFracVSeta_lowPt->getRootObject()) mNHFracVSeta_lowPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].neutralHadronEnergyFraction());
1736 mPhFracVSeta_lowPt = map_of_MEs[
DirName+
"/"+
"PhFracVSeta_lowPt"];
if (mPhFracVSeta_lowPt && mPhFracVSeta_lowPt->getRootObject()) mPhFracVSeta_lowPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].neutralEmEnergyFraction());
1737 mElFracVSeta_lowPt = map_of_MEs[
DirName+
"/"+
"ElFracVSeta_lowPt"];
if (mElFracVSeta_lowPt && mElFracVSeta_lowPt->getRootObject()) mElFracVSeta_lowPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].chargedEmEnergyFraction());
1738 mMuFracVSeta_lowPt = map_of_MEs[
DirName+
"/"+
"MuFracVSeta_lowPt"];
if (mMuFracVSeta_lowPt && mMuFracVSeta_lowPt->getRootObject()) mMuFracVSeta_lowPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].chargedMuEnergyFraction());
1740 if ((*pfJets)[ijet].pt()>50. && (*pfJets)[ijet].pt()<=140.) {
1741 mCHFracVSeta_mediumPt = map_of_MEs[
DirName+
"/"+
"CHFracVSeta_mediumPt"];
if (mCHFracVSeta_mediumPt && mCHFracVSeta_mediumPt->getRootObject()) mCHFracVSeta_mediumPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].chargedHadronEnergyFraction());
1742 mNHFracVSeta_mediumPt = map_of_MEs[
DirName+
"/"+
"NHFracVSeta_mediumPt"];
if (mNHFracVSeta_mediumPt && mNHFracVSeta_mediumPt->getRootObject()) mNHFracVSeta_mediumPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].neutralHadronEnergyFraction());
1743 mPhFracVSeta_mediumPt = map_of_MEs[
DirName+
"/"+
"PhFracVSeta_mediumPt"];
if (mPhFracVSeta_mediumPt && mPhFracVSeta_mediumPt->getRootObject()) mPhFracVSeta_mediumPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].neutralEmEnergyFraction());
1744 mElFracVSeta_mediumPt = map_of_MEs[
DirName+
"/"+
"ElFracVSeta_mediumPt"];
if (mElFracVSeta_mediumPt && mElFracVSeta_mediumPt->getRootObject()) mElFracVSeta_mediumPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].chargedEmEnergyFraction());
1745 mMuFracVSeta_mediumPt = map_of_MEs[
DirName+
"/"+
"MuFracVSeta_mediumPt"];
if (mMuFracVSeta_mediumPt && mMuFracVSeta_mediumPt->getRootObject()) mMuFracVSeta_mediumPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].chargedMuEnergyFraction());
1747 if ((*pfJets)[ijet].pt()>140.) {
1748 mCHFracVSeta_highPt = map_of_MEs[
DirName+
"/"+
"CHFracVSeta_highPt"];
if (mCHFracVSeta_highPt && mCHFracVSeta_highPt->getRootObject()) mCHFracVSeta_highPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].chargedHadronEnergyFraction());
1749 mNHFracVSeta_highPt = map_of_MEs[
DirName+
"/"+
"NHFracVSeta_highPt"];
if (mNHFracVSeta_highPt && mNHFracVSeta_highPt->getRootObject()) mNHFracVSeta_highPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].neutralHadronEnergyFraction());
1750 mPhFracVSeta_highPt = map_of_MEs[
DirName+
"/"+
"PhFracVSeta_highPt"];
if (mPhFracVSeta_highPt && mPhFracVSeta_highPt->getRootObject()) mPhFracVSeta_highPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].neutralEmEnergyFraction());
1751 mElFracVSeta_highPt = map_of_MEs[
DirName+
"/"+
"ElFracVSeta_highPt"];
if (mElFracVSeta_highPt && mElFracVSeta_highPt->getRootObject()) mElFracVSeta_highPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].chargedEmEnergyFraction());
1752 mMuFracVSeta_highPt = map_of_MEs[
DirName+
"/"+
"MuFracVSeta_highPt"];
if (mMuFracVSeta_highPt && mMuFracVSeta_highPt->getRootObject()) mMuFracVSeta_highPt->Fill((*pfJets)[ijet].eta(),(*pfJets)[ijet].chargedMuEnergyFraction());
1754 if (fabs((*pfJets)[ijet].
eta()) <= 1.3) {
1755 mHFrac_Barrel = map_of_MEs[
DirName+
"/"+
"HFrac_Barrel"];
if (mHFrac_Barrel && mHFrac_Barrel->getRootObject()) mHFrac_Barrel->Fill((*pfJets)[ijet].chargedHadronEnergyFraction() + (*pfJets)[ijet].neutralHadronEnergyFraction() );
1756 mEFrac_Barrel = map_of_MEs[
DirName+
"/"+
"EFrac_Barrel"];
if (mEFrac_Barrel && mEFrac_Barrel->getRootObject()) mEFrac_Barrel->Fill ((*pfJets)[ijet].chargedEmEnergyFraction() + (*pfJets)[ijet].neutralEmEnergyFraction());
1758 if ((*pfJets)[ijet].pt()<=50.) {
1759 mCHFrac_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"CHFrac_lowPt_Barrel"];
if (mCHFrac_lowPt_Barrel && mCHFrac_lowPt_Barrel->getRootObject()) mCHFrac_lowPt_Barrel->Fill((*pfJets)[ijet].chargedHadronEnergyFraction());
1760 mNHFrac_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"NHFrac_lowPt_Barrel"];
if (mNHFrac_lowPt_Barrel && mNHFrac_lowPt_Barrel->getRootObject()) mNHFrac_lowPt_Barrel->Fill((*pfJets)[ijet].neutralHadronEnergyFraction());
1761 mPhFrac_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"PhFrac_lowPt_Barrel"];
if (mPhFrac_lowPt_Barrel && mPhFrac_lowPt_Barrel->getRootObject()) mPhFrac_lowPt_Barrel->Fill((*pfJets)[ijet].neutralEmEnergyFraction());
1762 mElFrac_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"ElFrac_lowPt_Barrel"];
if (mElFrac_lowPt_Barrel && mElFrac_lowPt_Barrel->getRootObject()) mElFrac_lowPt_Barrel->Fill((*pfJets)[ijet].chargedEmEnergyFraction());
1763 mMuFrac_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"MuFrac_lowPt_Barrel"];
if (mMuFrac_lowPt_Barrel && mMuFrac_lowPt_Barrel->getRootObject()) mMuFrac_lowPt_Barrel->Fill((*pfJets)[ijet].chargedMuEnergyFraction());
1764 mCHEn_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"CHEn_lowPt_Barrel"];
if (mCHEn_lowPt_Barrel && mCHEn_lowPt_Barrel->getRootObject()) mCHEn_lowPt_Barrel->Fill((*pfJets)[ijet].chargedHadronEnergy());
1765 mNHEn_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"NHEn_lowPt_Barrel"];
if (mNHEn_lowPt_Barrel && mNHEn_lowPt_Barrel->getRootObject()) mNHEn_lowPt_Barrel->Fill((*pfJets)[ijet].neutralHadronEnergy());
1766 mPhEn_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"PhEn_lowPt_Barrel"];
if (mPhEn_lowPt_Barrel && mPhEn_lowPt_Barrel->getRootObject()) mPhEn_lowPt_Barrel->Fill((*pfJets)[ijet].neutralEmEnergy());
1767 mElEn_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"ElEn_lowPt_Barrel"];
if (mElEn_lowPt_Barrel && mElEn_lowPt_Barrel->getRootObject()) mElEn_lowPt_Barrel->Fill((*pfJets)[ijet].chargedEmEnergy());
1768 mMuEn_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"MuEn_lowPt_Barrel"];
if (mMuEn_lowPt_Barrel && mMuEn_lowPt_Barrel->getRootObject()) mMuEn_lowPt_Barrel->Fill((*pfJets)[ijet].chargedMuEnergy());
1769 mChMultiplicity_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"ChMultiplicity_lowPt_Barrel"];
if(mChMultiplicity_lowPt_Barrel && mChMultiplicity_lowPt_Barrel->getRootObject()) mChMultiplicity_lowPt_Barrel->Fill((*pfJets)[ijet].chargedMultiplicity());
1770 mNeutMultiplicity_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"NeutMultiplicity_lowPt_Barrel"];
if(mNeutMultiplicity_lowPt_Barrel && mNeutMultiplicity_lowPt_Barrel->getRootObject()) mNeutMultiplicity_lowPt_Barrel->Fill((*pfJets)[ijet].neutralMultiplicity());
1771 mMuMultiplicity_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"MuMultiplicity_lowPt_Barrel"];
if(mMuMultiplicity_lowPt_Barrel && mMuMultiplicity_lowPt_Barrel->getRootObject()) mMuMultiplicity_lowPt_Barrel->Fill((*pfJets)[ijet].muonMultiplicity());
1773 if ((*pfJets)[ijet].pt()>50. && (*pfJets)[ijet].pt()<=140.) {
1774 mCHFrac_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"CHFrac_mediumPt_Barrel"];
if (mCHFrac_mediumPt_Barrel && mCHFrac_mediumPt_Barrel->getRootObject()) mCHFrac_mediumPt_Barrel->Fill((*pfJets)[ijet].chargedHadronEnergyFraction());
1775 mNHFrac_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"NHFrac_mediumPt_Barrel"];
if (mNHFrac_mediumPt_Barrel && mNHFrac_mediumPt_Barrel->getRootObject()) mNHFrac_mediumPt_Barrel->Fill((*pfJets)[ijet].neutralHadronEnergyFraction());
1776 mPhFrac_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"PhFrac_mediumPt_Barrel"];
if (mPhFrac_mediumPt_Barrel && mPhFrac_mediumPt_Barrel->getRootObject()) mPhFrac_mediumPt_Barrel->Fill((*pfJets)[ijet].neutralEmEnergyFraction());
1777 mElFrac_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"ElFrac_mediumPt_Barrel"];
if (mElFrac_mediumPt_Barrel && mElFrac_mediumPt_Barrel->getRootObject()) mElFrac_mediumPt_Barrel->Fill((*pfJets)[ijet].chargedEmEnergyFraction());
1778 mMuFrac_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"MuFrac_mediumPt_Barrel"];
if (mMuFrac_mediumPt_Barrel && mMuFrac_mediumPt_Barrel->getRootObject()) mMuFrac_mediumPt_Barrel->Fill((*pfJets)[ijet].chargedMuEnergyFraction());
1779 mCHEn_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"CHEn_mediumPt_Barrel"];
if (mCHEn_mediumPt_Barrel && mCHEn_mediumPt_Barrel->getRootObject()) mCHEn_mediumPt_Barrel->Fill((*pfJets)[ijet].chargedHadronEnergy());
1780 mNHEn_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"NHEn_mediumPt_Barrel"];
if (mNHEn_mediumPt_Barrel && mNHEn_mediumPt_Barrel->getRootObject()) mNHEn_mediumPt_Barrel->Fill((*pfJets)[ijet].neutralHadronEnergy());
1781 mPhEn_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"PhEn_mediumPt_Barrel"];
if (mPhEn_mediumPt_Barrel && mPhEn_mediumPt_Barrel->getRootObject()) mPhEn_mediumPt_Barrel->Fill((*pfJets)[ijet].neutralEmEnergy());
1782 mElEn_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"ElEn_mediumPt_Barrel"];
if (mElEn_mediumPt_Barrel && mElEn_mediumPt_Barrel->getRootObject()) mElEn_mediumPt_Barrel->Fill((*pfJets)[ijet].chargedEmEnergy());
1783 mMuEn_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"MuEn_mediumPt_Barrel"];
if (mMuEn_mediumPt_Barrel && mMuEn_mediumPt_Barrel->getRootObject()) mMuEn_mediumPt_Barrel->Fill((*pfJets)[ijet].chargedMuEnergy());
1784 mChMultiplicity_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"ChMultiplicity_mediumPt_Barrel"];
if(mChMultiplicity_mediumPt_Barrel && mChMultiplicity_mediumPt_Barrel->getRootObject()) mChMultiplicity_mediumPt_Barrel->Fill((*pfJets)[ijet].chargedMultiplicity());
1785 mNeutMultiplicity_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"NeutMultiplicity_mediumPt_Barrel"];
if(mNeutMultiplicity_mediumPt_Barrel && mNeutMultiplicity_mediumPt_Barrel->getRootObject()) mNeutMultiplicity_mediumPt_Barrel->Fill((*pfJets)[ijet].neutralMultiplicity());
1786 mMuMultiplicity_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"MuMultiplicity_mediumPt_Barrel"];
if(mMuMultiplicity_mediumPt_Barrel && mMuMultiplicity_mediumPt_Barrel->getRootObject()) mMuMultiplicity_mediumPt_Barrel->Fill((*pfJets)[ijet].muonMultiplicity());
1788 if ((*pfJets)[ijet].pt()>140.) {
1789 mCHFrac_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"CHFrac_highPt_Barrel"];
if (mCHFrac_highPt_Barrel && mCHFrac_highPt_Barrel->getRootObject()) mCHFrac_highPt_Barrel->Fill((*pfJets)[ijet].chargedHadronEnergyFraction());
1790 mNHFrac_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"NHFrac_highPt_Barrel"];
if (mNHFrac_highPt_Barrel && mNHFrac_highPt_Barrel->getRootObject()) mNHFrac_highPt_Barrel->Fill((*pfJets)[ijet].neutralHadronEnergyFraction());
1791 mPhFrac_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"PhFrac_highPt_Barrel"];
if (mPhFrac_highPt_Barrel && mPhFrac_highPt_Barrel->getRootObject()) mPhFrac_highPt_Barrel->Fill((*pfJets)[ijet].neutralEmEnergyFraction());
1792 mElFrac_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"ElFrac_highPt_Barrel"];
if (mElFrac_highPt_Barrel && mElFrac_highPt_Barrel->getRootObject()) mElFrac_highPt_Barrel->Fill((*pfJets)[ijet].chargedEmEnergyFraction());
1793 mMuFrac_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"MuFrac_highPt_Barrel"];
if (mMuFrac_highPt_Barrel && mMuFrac_highPt_Barrel->getRootObject()) mMuFrac_highPt_Barrel->Fill((*pfJets)[ijet].chargedMuEnergyFraction());
1794 mCHEn_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"CHEn_highPt_Barrel"];
if (mCHEn_highPt_Barrel && mCHEn_highPt_Barrel->getRootObject()) mCHEn_highPt_Barrel->Fill((*pfJets)[ijet].chargedHadronEnergy());
1795 mNHEn_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"NHEn_highPt_Barrel"];
if (mNHEn_highPt_Barrel && mNHEn_highPt_Barrel->getRootObject()) mNHEn_highPt_Barrel->Fill((*pfJets)[ijet].neutralHadronEnergy());
1796 mPhEn_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"PhEn_highPt_Barrel"];
if (mPhEn_highPt_Barrel && mPhEn_highPt_Barrel->getRootObject()) mPhEn_highPt_Barrel->Fill((*pfJets)[ijet].neutralEmEnergy());
1797 mElEn_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"ElEn_highPt_Barrel"];
if (mElEn_highPt_Barrel && mElEn_highPt_Barrel->getRootObject()) mElEn_highPt_Barrel->Fill((*pfJets)[ijet].chargedEmEnergy());
1798 mMuEn_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"MuEn_highPt_Barrel"];
if (mMuEn_highPt_Barrel && mMuEn_highPt_Barrel->getRootObject()) mMuEn_highPt_Barrel->Fill((*pfJets)[ijet].chargedMuEnergy());
1799 mChMultiplicity_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"ChMultiplicity_highPt_Barrel"];
if(mChMultiplicity_highPt_Barrel && mChMultiplicity_highPt_Barrel->getRootObject()) mChMultiplicity_highPt_Barrel->Fill((*pfJets)[ijet].chargedMultiplicity());
1800 mNeutMultiplicity_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"NeutMultiplicity_highPt_Barrel"];
if(mNeutMultiplicity_highPt_Barrel && mNeutMultiplicity_highPt_Barrel->getRootObject()) mNeutMultiplicity_highPt_Barrel->Fill((*pfJets)[ijet].neutralMultiplicity());
1801 mMuMultiplicity_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"MuMultiplicity_highPt_Barrel"];
if(mMuMultiplicity_highPt_Barrel && mMuMultiplicity_highPt_Barrel->getRootObject()) mMuMultiplicity_highPt_Barrel->Fill((*pfJets)[ijet].muonMultiplicity());
1803 mCHFracVSpT_Barrel = map_of_MEs[
DirName+
"/"+
"CHFracVSpT_Barrel"];
if(mCHFracVSpT_Barrel && mCHFracVSpT_Barrel->getRootObject()) mCHFracVSpT_Barrel->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].chargedHadronEnergyFraction());
1804 mNHFracVSpT_Barrel = map_of_MEs[
DirName+
"/"+
"NHFracVSpT_Barrel"];
if (mNHFracVSpT_Barrel && mNHFracVSpT_Barrel->getRootObject()) mNHFracVSpT_Barrel->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].neutralHadronEnergyFraction());
1805 mPhFracVSpT_Barrel = map_of_MEs[
DirName+
"/"+
"PhFracVSpT_Barrel"];
if (mPhFracVSpT_Barrel && mPhFracVSpT_Barrel->getRootObject()) mPhFracVSpT_Barrel->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].neutralEmEnergyFraction());
1806 mElFracVSpT_Barrel = map_of_MEs[
DirName+
"/"+
"ElFracVSpT_Barrel"];
if (mElFracVSpT_Barrel && mElFracVSpT_Barrel->getRootObject()) mElFracVSpT_Barrel->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].chargedEmEnergyFraction());
1807 mMuFracVSpT_Barrel = map_of_MEs[
DirName+
"/"+
"MuFracVSpT_Barrel"];
if (mMuFracVSpT_Barrel && mMuFracVSpT_Barrel->getRootObject()) mMuFracVSpT_Barrel->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].chargedMuEnergyFraction());
1808 }
else if(fabs((*pfJets)[ijet].
eta()) <= 3) {
1809 mHFrac_EndCap = map_of_MEs[
DirName+
"/"+
"HFrac_EndCap"];
if (mHFrac_EndCap && mHFrac_EndCap->getRootObject()) mHFrac_EndCap->Fill((*pfJets)[ijet].chargedHadronEnergyFraction() + (*pfJets)[ijet].neutralHadronEnergyFraction()+(*pfJets)[ijet].HFHadronEnergyFraction ());
1810 mEFrac_EndCap = map_of_MEs[
DirName+
"/"+
"EFrac_EndCap"];
if (mEFrac_EndCap && mEFrac_EndCap->getRootObject()) mEFrac_EndCap->Fill ((*pfJets)[ijet].chargedEmEnergyFraction() + (*pfJets)[ijet].neutralEmEnergyFraction()+(*pfJets)[ijet].HFEMEnergyFraction ());
1812 if ((*pfJets)[ijet].pt()<=50.) {
1813 mCHFrac_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"CHFrac_lowPt_EndCap"];
if (mCHFrac_lowPt_EndCap && mCHFrac_lowPt_EndCap->getRootObject()) mCHFrac_lowPt_EndCap->Fill((*pfJets)[ijet].chargedHadronEnergyFraction());
1814 mNHFrac_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"NHFrac_lowPt_EndCap"];
if (mNHFrac_lowPt_EndCap && mNHFrac_lowPt_EndCap->getRootObject()) mNHFrac_lowPt_EndCap->Fill((*pfJets)[ijet].neutralHadronEnergyFraction());
1815 mPhFrac_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"PhFrac_lowPt_EndCap"];
if (mPhFrac_lowPt_EndCap && mPhFrac_lowPt_EndCap->getRootObject()) mPhFrac_lowPt_EndCap->Fill((*pfJets)[ijet].neutralEmEnergyFraction());
1816 mElFrac_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"ElFrac_lowPt_EndCap"];
if (mElFrac_lowPt_EndCap && mElFrac_lowPt_EndCap->getRootObject()) mElFrac_lowPt_EndCap->Fill((*pfJets)[ijet].chargedEmEnergyFraction());
1817 mMuFrac_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"MuFrac_lowPt_EndCap"];
if (mMuFrac_lowPt_EndCap && mMuFrac_lowPt_EndCap->getRootObject()) mMuFrac_lowPt_EndCap->Fill((*pfJets)[ijet].chargedMuEnergyFraction());
1818 mCHEn_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"CHEn_lowPt_EndCap"];
if (mCHEn_lowPt_EndCap && mCHEn_lowPt_EndCap->getRootObject()) mCHEn_lowPt_EndCap->Fill((*pfJets)[ijet].chargedHadronEnergy());
1819 mNHEn_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"NHEn_lowPt_EndCap"];
if (mNHEn_lowPt_EndCap && mNHEn_lowPt_EndCap->getRootObject()) mNHEn_lowPt_EndCap->Fill((*pfJets)[ijet].neutralHadronEnergy());
1820 mPhEn_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"PhEn_lowPt_EndCap"];
if (mPhEn_lowPt_EndCap && mPhEn_lowPt_EndCap->getRootObject()) mPhEn_lowPt_EndCap->Fill((*pfJets)[ijet].neutralEmEnergy());
1821 mElEn_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"ElEn_lowPt_EndCap"];
if (mElEn_lowPt_EndCap && mElEn_lowPt_EndCap->getRootObject()) mElEn_lowPt_EndCap->Fill((*pfJets)[ijet].chargedEmEnergy());
1822 mMuEn_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"MuEn_lowPt_EndCap"];
if (mMuEn_lowPt_EndCap && mMuEn_lowPt_EndCap->getRootObject()) mMuEn_lowPt_EndCap->Fill((*pfJets)[ijet].chargedMuEnergy());
1823 mChMultiplicity_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"ChMultiplicity_lowPt_EndCap"];
if(mChMultiplicity_lowPt_EndCap && mChMultiplicity_lowPt_EndCap->getRootObject()) mChMultiplicity_lowPt_EndCap->Fill((*pfJets)[ijet].chargedMultiplicity());
1824 mNeutMultiplicity_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"NeutMultiplicity_lowPt_EndCap"];
if(mNeutMultiplicity_lowPt_EndCap && mNeutMultiplicity_lowPt_EndCap->getRootObject()) mNeutMultiplicity_lowPt_EndCap->Fill((*pfJets)[ijet].neutralMultiplicity());
1825 mMuMultiplicity_lowPt_EndCap = map_of_MEs[
DirName+
"/"+
"MuMultiplicity_lowPt_EndCap"];
if(mMuMultiplicity_lowPt_EndCap && mMuMultiplicity_lowPt_EndCap->getRootObject()) mMuMultiplicity_lowPt_EndCap->Fill((*pfJets)[ijet].muonMultiplicity());
1827 if ((*pfJets)[ijet].pt()>50. && (*pfJets)[ijet].pt()<=140.) {
1828 mCHFrac_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"CHFrac_mediumPt_EndCap"];
if (mCHFrac_mediumPt_EndCap && mCHFrac_mediumPt_EndCap->getRootObject()) mCHFrac_mediumPt_EndCap->Fill((*pfJets)[ijet].chargedHadronEnergyFraction());
1829 mNHFrac_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"NHFrac_mediumPt_EndCap"];
if (mNHFrac_mediumPt_EndCap && mNHFrac_mediumPt_EndCap->getRootObject()) mNHFrac_mediumPt_EndCap->Fill((*pfJets)[ijet].neutralHadronEnergyFraction());
1830 mPhFrac_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"PhFrac_mediumPt_EndCap"];
if (mPhFrac_mediumPt_EndCap && mPhFrac_mediumPt_EndCap->getRootObject()) mPhFrac_mediumPt_EndCap->Fill((*pfJets)[ijet].neutralEmEnergyFraction());
1831 mElFrac_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"ElFrac_mediumPt_EndCap"];
if (mElFrac_mediumPt_EndCap && mElFrac_mediumPt_EndCap->getRootObject()) mElFrac_mediumPt_EndCap->Fill((*pfJets)[ijet].chargedEmEnergyFraction());
1832 mMuFrac_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"MuFrac_mediumPt_EndCap"];
if (mMuFrac_mediumPt_EndCap && mMuFrac_mediumPt_EndCap->getRootObject()) mMuFrac_mediumPt_EndCap->Fill((*pfJets)[ijet].chargedMuEnergyFraction());
1833 mCHEn_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"CHEn_mediumPt_EndCap"];
if (mCHEn_mediumPt_EndCap && mCHEn_mediumPt_EndCap->getRootObject()) mCHEn_mediumPt_EndCap->Fill((*pfJets)[ijet].chargedHadronEnergy());
1834 mNHEn_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"NHEn_mediumPt_EndCap"];
if (mNHEn_mediumPt_EndCap && mNHEn_mediumPt_EndCap->getRootObject()) mNHEn_mediumPt_EndCap->Fill((*pfJets)[ijet].neutralHadronEnergy());
1835 mPhEn_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"PhEn_mediumPt_EndCap"];
if (mPhEn_mediumPt_EndCap && mPhEn_mediumPt_EndCap->getRootObject()) mPhEn_mediumPt_EndCap->Fill((*pfJets)[ijet].neutralEmEnergy());
1836 mElEn_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"ElEn_mediumPt_EndCap"];
if (mElEn_mediumPt_EndCap && mElEn_mediumPt_EndCap->getRootObject()) mElEn_mediumPt_EndCap->Fill((*pfJets)[ijet].chargedEmEnergy());
1837 mMuEn_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"MuEn_mediumPt_EndCap"];
if (mMuEn_mediumPt_EndCap && mMuEn_mediumPt_EndCap->getRootObject()) mMuEn_mediumPt_EndCap->Fill((*pfJets)[ijet].chargedMuEnergy());
1838 mChMultiplicity_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"ChMultiplicity_mediumPt_EndCap"];
if(mChMultiplicity_mediumPt_EndCap && mChMultiplicity_mediumPt_EndCap->getRootObject()) mChMultiplicity_mediumPt_EndCap->Fill((*pfJets)[ijet].chargedMultiplicity());
1839 mNeutMultiplicity_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"NeutMultiplicity_mediumPt_EndCap"];
if(mNeutMultiplicity_mediumPt_EndCap && mNeutMultiplicity_mediumPt_EndCap->getRootObject()) mNeutMultiplicity_mediumPt_EndCap->Fill((*pfJets)[ijet].neutralMultiplicity());
1840 mMuMultiplicity_mediumPt_EndCap = map_of_MEs[
DirName+
"/"+
"MuMultiplicity_mediumPt_EndCap"];
if(mMuMultiplicity_mediumPt_EndCap && mMuMultiplicity_mediumPt_EndCap->getRootObject()) mMuMultiplicity_mediumPt_EndCap->Fill((*pfJets)[ijet].muonMultiplicity());
1842 if ((*pfJets)[ijet].pt()>140.) {
1843 mCHFrac_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"CHFrac_highPt_EndCap"];
if (mCHFrac_highPt_EndCap && mCHFrac_highPt_EndCap->getRootObject()) mCHFrac_highPt_EndCap->Fill((*pfJets)[ijet].chargedHadronEnergyFraction());
1844 mNHFrac_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"NHFrac_highPt_EndCap"];
if (mNHFrac_highPt_EndCap && mNHFrac_highPt_EndCap->getRootObject()) mNHFrac_highPt_EndCap->Fill((*pfJets)[ijet].neutralHadronEnergyFraction());
1845 mPhFrac_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"PhFrac_highPt_EndCap"];
if (mPhFrac_highPt_EndCap && mPhFrac_highPt_EndCap->getRootObject()) mPhFrac_highPt_EndCap->Fill((*pfJets)[ijet].neutralEmEnergyFraction());
1846 mElFrac_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"ElFrac_highPt_EndCap"];
if (mElFrac_highPt_EndCap && mElFrac_highPt_EndCap->getRootObject()) mElFrac_highPt_EndCap->Fill((*pfJets)[ijet].chargedEmEnergyFraction());
1847 mMuFrac_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"MuFrac_highPt_EndCap"];
if (mMuFrac_highPt_EndCap && mMuFrac_highPt_EndCap->getRootObject()) mMuFrac_highPt_EndCap->Fill((*pfJets)[ijet].chargedMuEnergyFraction());
1848 mCHEn_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"CHEn_highPt_EndCap"];
if (mCHEn_highPt_EndCap && mCHEn_highPt_EndCap->getRootObject()) mCHEn_highPt_EndCap->Fill((*pfJets)[ijet].chargedHadronEnergy());
1849 mNHEn_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"NHEn_highPt_EndCap"];
if (mNHEn_highPt_EndCap && mNHEn_highPt_EndCap->getRootObject()) mNHEn_highPt_EndCap->Fill((*pfJets)[ijet].neutralHadronEnergy());
1850 mPhEn_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"PhEn_highPt_EndCap"];
if (mPhEn_highPt_EndCap && mPhEn_highPt_EndCap->getRootObject()) mPhEn_highPt_EndCap->Fill((*pfJets)[ijet].neutralEmEnergy());
1851 mElEn_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"ElEn_highPt_EndCap"];
if (mElEn_highPt_EndCap && mElEn_highPt_EndCap->getRootObject()) mElEn_highPt_EndCap->Fill((*pfJets)[ijet].chargedEmEnergy());
1852 mMuEn_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"MuEn_highPt_EndCap"];
if (mMuEn_highPt_EndCap && mMuEn_highPt_EndCap->getRootObject()) mMuEn_highPt_EndCap->Fill((*pfJets)[ijet].chargedMuEnergy());
1853 mChMultiplicity_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"ChMultiplicity_highPt_EndCap"];
if(mChMultiplicity_highPt_EndCap && mChMultiplicity_highPt_EndCap->getRootObject()) mChMultiplicity_highPt_EndCap->Fill((*pfJets)[ijet].chargedMultiplicity());
1854 mNeutMultiplicity_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"NeutMultiplicity_highPt_EndCap"];
if(mNeutMultiplicity_highPt_EndCap && mNeutMultiplicity_highPt_EndCap->getRootObject()) mNeutMultiplicity_highPt_EndCap->Fill((*pfJets)[ijet].neutralMultiplicity());
1855 mMuMultiplicity_highPt_EndCap = map_of_MEs[
DirName+
"/"+
"MuMultiplicity_highPt_EndCap"];
if(mMuMultiplicity_highPt_EndCap && mMuMultiplicity_highPt_EndCap->getRootObject()) mMuMultiplicity_highPt_EndCap->Fill((*pfJets)[ijet].muonMultiplicity());
1857 mCHFracVSpT_EndCap = map_of_MEs[
DirName+
"/"+
"CHFracVSpT_EndCap"];
if(mCHFracVSpT_EndCap && mCHFracVSpT_EndCap->getRootObject()) mCHFracVSpT_EndCap->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].chargedHadronEnergyFraction());
1858 mNHFracVSpT_EndCap = map_of_MEs[
DirName+
"/"+
"NHFracVSpT_EndCap"];
if (mNHFracVSpT_EndCap && mNHFracVSpT_EndCap->getRootObject()) mNHFracVSpT_EndCap->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].neutralHadronEnergyFraction());
1859 mPhFracVSpT_EndCap = map_of_MEs[
DirName+
"/"+
"PhFracVSpT_EndCap"];
if (mPhFracVSpT_EndCap && mPhFracVSpT_EndCap->getRootObject()) mPhFracVSpT_EndCap->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].neutralEmEnergyFraction());
1860 mElFracVSpT_EndCap = map_of_MEs[
DirName+
"/"+
"ElFracVSpT_EndCap"];
if (mElFracVSpT_EndCap && mElFracVSpT_EndCap->getRootObject()) mElFracVSpT_EndCap->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].chargedEmEnergyFraction());
1861 mMuFracVSpT_EndCap = map_of_MEs[
DirName+
"/"+
"MuFracVSpT_EndCap"];
if (mMuFracVSpT_EndCap && mMuFracVSpT_EndCap->getRootObject()) mMuFracVSpT_EndCap->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].chargedMuEnergyFraction());
1863 mHFrac_Forward = map_of_MEs[
DirName+
"/"+
"HFrac_Forward"];
if (mHFrac_Forward && mHFrac_Forward->getRootObject()) mHFrac_Forward->Fill((*pfJets)[ijet].chargedHadronEnergyFraction() + (*pfJets)[ijet].neutralHadronEnergyFraction()+(*pfJets)[ijet].HFHadronEnergyFraction ());
1864 mEFrac_Forward = map_of_MEs[
DirName+
"/"+
"EFrac_Forward"];
if (mEFrac_Forward && mEFrac_Forward->getRootObject()) mEFrac_Forward->Fill ((*pfJets)[ijet].chargedEmEnergyFraction() + (*pfJets)[ijet].neutralEmEnergyFraction()+(*pfJets)[ijet].HFEMEnergyFraction ());
1865 mHFHFracVSpT_Forward = map_of_MEs[
DirName+
"/"+
"HFHFracVSpT_Forward"];
if (mHFHFracVSpT_Forward && mHFHFracVSpT_Forward->getRootObject()) mHFHFracVSpT_Forward->Fill((*pfJets)[ijet].pt(),(*pfJets)[ijet].chargedHadronEnergyFraction() + (*pfJets)[ijet].neutralHadronEnergyFraction()+(*pfJets)[ijet].HFHadronEnergyFraction ());
1866 mHFEFracVSpT_Forward = map_of_MEs[
DirName+
"/"+
"HFEFracVSpT_Forward"];
if (mHFEFracVSpT_Forward && mHFEFracVSpT_Forward->getRootObject()) mHFEFracVSpT_Forward->Fill ((*pfJets)[ijet].pt(),(*pfJets)[ijet].chargedEmEnergyFraction() + (*pfJets)[ijet].neutralEmEnergyFraction()+(*pfJets)[ijet].HFEMEnergyFraction ());
1868 if ((*pfJets)[ijet].pt()<=50.) {
1869 mHFEFrac_lowPt_Forward = map_of_MEs[
DirName+
"/"+
"HFEFrac_lowPt_Forward"];
if(mHFEFrac_lowPt_Forward && mHFEFrac_lowPt_Forward->getRootObject()) mHFEFrac_lowPt_Forward->Fill((*pfJets)[ijet].HFEMEnergyFraction());
1870 mHFEFrac_lowPt_Forward = map_of_MEs[
DirName+
"/"+
"HFEFrac_lowPt_Forward"];
if(mHFHFrac_lowPt_Forward && mHFHFrac_lowPt_Forward->getRootObject()) mHFHFrac_lowPt_Forward->Fill((*pfJets)[ijet].HFHadronEnergyFraction());
1871 mHFEEn_lowPt_Forward = map_of_MEs[
DirName+
"/"+
"HFEEn_lowPt_Forward"];
if(mHFEEn_lowPt_Forward && mHFEEn_lowPt_Forward->getRootObject()) mHFEEn_lowPt_Forward->Fill((*pfJets)[ijet].HFEMEnergy());
1872 mHFHEn_lowPt_Forward = map_of_MEs[
DirName+
"/"+
"HFHEn_lowPt_Forward"];
if(mHFHEn_lowPt_Forward && mHFHEn_lowPt_Forward->getRootObject()) mHFHEn_lowPt_Forward->Fill((*pfJets)[ijet].HFHadronEnergy());
1873 mChMultiplicity_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"ChMultiplicity_lowPt_Barrel"];
if(mChMultiplicity_lowPt_Forward && mChMultiplicity_lowPt_Forward->getRootObject()) mChMultiplicity_lowPt_Forward->Fill((*pfJets)[ijet].chargedMultiplicity());
1874 mNeutMultiplicity_lowPt_Barrel = map_of_MEs[
DirName+
"/"+
"NeutMultiplicity_lowPt_Barrel"];
if(mNeutMultiplicity_lowPt_Forward && mNeutMultiplicity_lowPt_Forward->getRootObject()) mNeutMultiplicity_lowPt_Forward->Fill((*pfJets)[ijet].neutralMultiplicity());
1876 if ((*pfJets)[ijet].pt()>50. && (*pfJets)[ijet].pt()<=140.) {
1877 mHFEFrac_mediumPt_Forward = map_of_MEs[
DirName+
"/"+
"HFEFrac_mediumPt_Forward"];
if(mHFEFrac_mediumPt_Forward && mHFEFrac_mediumPt_Forward->getRootObject()) mHFEFrac_mediumPt_Forward->Fill((*pfJets)[ijet].HFEMEnergyFraction());
1878 mHFEFrac_mediumPt_Forward = map_of_MEs[
DirName+
"/"+
"HFEFrac_mediumPt_Forward"];
if(mHFHFrac_mediumPt_Forward && mHFHFrac_mediumPt_Forward->getRootObject()) mHFHFrac_mediumPt_Forward->Fill((*pfJets)[ijet].HFHadronEnergyFraction());
1879 mHFEEn_mediumPt_Forward = map_of_MEs[
DirName+
"/"+
"HFEEn_mediumPt_Forward"];
if(mHFEEn_mediumPt_Forward && mHFEEn_mediumPt_Forward->getRootObject()) mHFEEn_mediumPt_Forward->Fill((*pfJets)[ijet].HFEMEnergy());
1880 mHFHEn_mediumPt_Forward = map_of_MEs[
DirName+
"/"+
"HFHEn_mediumPt_Forward"];
if(mHFHEn_mediumPt_Forward && mHFHEn_mediumPt_Forward->getRootObject()) mHFHEn_mediumPt_Forward->Fill((*pfJets)[ijet].HFHadronEnergy());
1881 mChMultiplicity_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"ChMultiplicity_mediumPt_Barrel"];
if(mChMultiplicity_mediumPt_Forward && mChMultiplicity_mediumPt_Forward->getRootObject()) mChMultiplicity_mediumPt_Forward->Fill((*pfJets)[ijet].chargedMultiplicity());
1882 mNeutMultiplicity_mediumPt_Barrel = map_of_MEs[
DirName+
"/"+
"NeutMultiplicity_mediumPt_Barrel"];
if(mNeutMultiplicity_mediumPt_Forward && mNeutMultiplicity_mediumPt_Forward->getRootObject()) mNeutMultiplicity_mediumPt_Forward->Fill((*pfJets)[ijet].neutralMultiplicity());
1884 if ((*pfJets)[ijet].pt()>140.) {
1885 mHFEFrac_highPt_Forward = map_of_MEs[
DirName+
"/"+
"HFEFrac_highPt_Forward"];
if(mHFEFrac_highPt_Forward && mHFEFrac_highPt_Forward->getRootObject()) mHFEFrac_highPt_Forward->Fill((*pfJets)[ijet].HFEMEnergyFraction());
1886 mHFEFrac_highPt_Forward = map_of_MEs[
DirName+
"/"+
"HFEFrac_highPt_Forward"];
if(mHFHFrac_highPt_Forward && mHFHFrac_highPt_Forward->getRootObject()) mHFHFrac_highPt_Forward->Fill((*pfJets)[ijet].HFHadronEnergyFraction());
1887 mHFEEn_highPt_Forward = map_of_MEs[
DirName+
"/"+
"HFEEn_highPt_Forward"];
if(mHFEEn_highPt_Forward && mHFEEn_highPt_Forward->getRootObject()) mHFEEn_highPt_Forward->Fill((*pfJets)[ijet].HFEMEnergy());
1888 mHFHEn_highPt_Forward = map_of_MEs[
DirName+
"/"+
"HFHEn_highPt_Forward"];
if(mHFHEn_highPt_Forward && mHFHEn_highPt_Forward->getRootObject()) mHFHEn_highPt_Forward->Fill((*pfJets)[ijet].HFHadronEnergy());
1889 mChMultiplicity_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"ChMultiplicity_highPt_Barrel"];
if(mChMultiplicity_highPt_Forward && mChMultiplicity_highPt_Forward->getRootObject()) mChMultiplicity_highPt_Forward->Fill((*pfJets)[ijet].chargedMultiplicity());
1890 mNeutMultiplicity_highPt_Barrel = map_of_MEs[
DirName+
"/"+
"NeutMultiplicity_highPt_Barrel"];
if(mNeutMultiplicity_highPt_Forward && mNeutMultiplicity_highPt_Forward->getRootObject()) mNeutMultiplicity_highPt_Forward->Fill((*pfJets)[ijet].neutralMultiplicity());
1893 mChargedHadronEnergy = map_of_MEs[
DirName+
"/"+
"ChargedHadronEnergy"];
if (mChargedHadronEnergy && mChargedHadronEnergy->getRootObject()) mChargedHadronEnergy->Fill ((*pfJets)[ijet].chargedHadronEnergy());
1894 mNeutralHadronEnergy = map_of_MEs[
DirName+
"/"+
"NeutralHadronEnergy"];
if (mNeutralHadronEnergy && mNeutralHadronEnergy->getRootObject()) mNeutralHadronEnergy->Fill ((*pfJets)[ijet].neutralHadronEnergy());
1895 mChargedEmEnergy = map_of_MEs[
DirName+
"/"+
"ChargedEmEnergy"];
if (mChargedEmEnergy && mChargedEmEnergy->getRootObject()) mChargedEmEnergy->Fill((*pfJets)[ijet].chargedEmEnergy());
1896 mChargedMuEnergy = map_of_MEs[
DirName+
"/"+
"ChargedMuEnergy"];
if (mChargedMuEnergy && mChargedMuEnergy->getRootObject()) mChargedMuEnergy->Fill ((*pfJets)[ijet].chargedMuEnergy ());
1897 mNeutralEmEnergy = map_of_MEs[
DirName+
"/"+
"NeutralEmEnergy"];
if (mNeutralEmEnergy && mNeutralEmEnergy->getRootObject()) mNeutralEmEnergy->Fill((*pfJets)[ijet].neutralEmEnergy());
1898 mChargedMultiplicity = map_of_MEs[
DirName+
"/"+
"ChargedMultiplicity"];
if (mChargedMultiplicity && mChargedMultiplicity->getRootObject()) mChargedMultiplicity->Fill((*pfJets)[ijet].chargedMultiplicity());
1899 mNeutralMultiplicity = map_of_MEs[
DirName+
"/"+
"NeutralMultiplicity"];
if (mNeutralMultiplicity && mNeutralMultiplicity->getRootObject()) mNeutralMultiplicity->Fill((*pfJets)[ijet].neutralMultiplicity());
1900 mMuonMultiplicity = map_of_MEs[
DirName+
"/"+
"MuonMultiplicity"];
if (mMuonMultiplicity && mMuonMultiplicity->getRootObject()) mMuonMultiplicity->Fill ((*pfJets)[ijet].
muonMultiplicity());
1902 mNeutralFraction = map_of_MEs[
DirName+
"/"+
"NeutralConstituentsFraction"];
if (mNeutralFraction && mNeutralFraction->getRootObject()) mNeutralFraction->Fill ((
double)(*pfJets)[ijet].neutralMultiplicity()/(double)(*pfJets)[ijet].nConstituents());
1904 mHFrac_profile = map_of_MEs[
DirName+
"/"+
"HFrac_profile"];
if (mHFrac_profile && mHFrac_profile->getRootObject()) mHFrac_profile ->
Fill(numPV, (*pfJets)[ijet].chargedHadronEnergyFraction() + (*pfJets)[ijet].neutralHadronEnergyFraction());
1905 mEFrac_profile = map_of_MEs[
DirName+
"/"+
"EFrac_profile"];
if (mEFrac_profile && mEFrac_profile->getRootObject()) mEFrac_profile ->
Fill(numPV, (*pfJets)[ijet].chargedEmEnergyFraction() + (*pfJets)[ijet].neutralEmEnergyFraction());
1907 mChargedHadronEnergy_profile = map_of_MEs[
DirName+
"/"+
"ChargedHadronEnergy_profile"];
if (mChargedHadronEnergy_profile && mChargedHadronEnergy_profile->getRootObject()) mChargedHadronEnergy_profile->Fill(numPV, (*pfJets)[ijet].chargedHadronEnergy());
1908 mNeutralHadronEnergy_profile = map_of_MEs[
DirName+
"/"+
"NeutralHadronEnergy_profile"];
if (mNeutralHadronEnergy_profile && mNeutralHadronEnergy_profile->getRootObject()) mNeutralHadronEnergy_profile->Fill(numPV, (*pfJets)[ijet].neutralHadronEnergy());
1909 mChargedEmEnergy_profile = map_of_MEs[
DirName+
"/"+
"ChargedEmEnergy_profile"];
if (mChargedEmEnergy_profile && mChargedEmEnergy_profile->getRootObject()) mChargedEmEnergy_profile ->
Fill(numPV, (*pfJets)[ijet].chargedEmEnergy());
1910 mChargedMuEnergy_profile = map_of_MEs[
DirName+
"/"+
"ChargedMuEnergy_profile"];
if (mChargedMuEnergy_profile && mChargedMuEnergy_profile->getRootObject()) mChargedMuEnergy_profile ->
Fill(numPV, (*pfJets)[ijet].chargedMuEnergy ());
1911 mNeutralEmEnergy_profile = map_of_MEs[
DirName+
"/"+
"NeutralEmEnergy_profile"];
if (mNeutralEmEnergy_profile && mNeutralEmEnergy_profile->getRootObject()) mNeutralEmEnergy_profile ->
Fill(numPV, (*pfJets)[ijet].neutralEmEnergy());
1912 mChargedMultiplicity_profile = map_of_MEs[
DirName+
"/"+
"ChargedMultiplicity_profile"];
if (mChargedMultiplicity_profile && mChargedMultiplicity_profile->getRootObject()) mChargedMultiplicity_profile->Fill(numPV, (*pfJets)[ijet].chargedMultiplicity());
1913 mNeutralMultiplicity_profile = map_of_MEs[
DirName+
"/"+
"NeutralMultiplicity_profile"];
if (mNeutralMultiplicity_profile && mNeutralMultiplicity_profile->getRootObject()) mNeutralMultiplicity_profile->Fill(numPV, (*pfJets)[ijet].neutralMultiplicity());
1914 mMuonMultiplicity_profile = map_of_MEs[
DirName+
"/"+
"MuonMultiplicity_profile"];
if (mMuonMultiplicity_profile && mMuonMultiplicity_profile->getRootObject()) mMuonMultiplicity_profile ->
Fill(numPV, (*pfJets)[ijet].muonMultiplicity());
1921 if(!pass_corrected){
1924 if(correctedJet.
pt()>pt1){
1927 cleaned_third_jet=cleaned_second_jet;
1930 cleaned_second_jet=cleaned_first_jet;
1931 pt1=correctedJet.
pt();
1933 cleaned_first_jet=jetpassid;
1934 }
else if(correctedJet.
pt()>pt2){
1937 cleaned_third_jet=cleaned_second_jet;
1938 pt2=correctedJet.
pt();
1940 cleaned_second_jet=jetpassid;
1941 }
else if(correctedJet.
pt()>pt3){
1942 pt3=correctedJet.
pt();
1944 cleaned_third_jet=jetpassid;
1946 if(cleaned_third_jet){
1950 if(!jetCleaningFlag_){
1952 mLooseJIDPassFractionVSeta = map_of_MEs[
DirName+
"/"+
"JetIDPassFractionVSeta"];
if (mLooseJIDPassFractionVSeta && mLooseJIDPassFractionVSeta->getRootObject()) mLooseJIDPassFractionVSeta->Fill(correctedJet.
eta(),1.);
1953 mLooseJIDPassFractionVSpt = map_of_MEs[
DirName+
"/"+
"JetIDPassFractionVSpt"];
if (mLooseJIDPassFractionVSpt && mLooseJIDPassFractionVSpt->getRootObject()) mLooseJIDPassFractionVSpt->Fill(correctedJet.
pt(),1.);
1954 if(correctedJet.
eta()<3.0){
1955 mLooseJIDPassFractionVSptNoHF= map_of_MEs[
DirName+
"/"+
"JetIDPassFractionVSptNoHF"];
if (mLooseJIDPassFractionVSptNoHF && mLooseJIDPassFractionVSptNoHF->getRootObject()) mLooseJIDPassFractionVSptNoHF->Fill(correctedJet.
pt(),1.);
1958 mLooseJIDPassFractionVSeta = map_of_MEs[
DirName+
"/"+
"JetIDPassFractionVSeta"];
if (mLooseJIDPassFractionVSeta && mLooseJIDPassFractionVSeta->getRootObject()) mLooseJIDPassFractionVSeta->Fill(correctedJet.
eta(),0.);
1959 mLooseJIDPassFractionVSpt = map_of_MEs[
DirName+
"/"+
"JetIDPassFractionVSpt"];
if (mLooseJIDPassFractionVSpt && mLooseJIDPassFractionVSpt->getRootObject()) mLooseJIDPassFractionVSpt->Fill(correctedJet.
pt(),0.);
1960 if(correctedJet.
eta()<3.0){
1961 mLooseJIDPassFractionVSptNoHF= map_of_MEs[
DirName+
"/"+
"JetIDPassFractionVSptNoHF"];
if (mLooseJIDPassFractionVSptNoHF && mLooseJIDPassFractionVSptNoHF->getRootObject()) mLooseJIDPassFractionVSptNoHF->Fill(correctedJet.
pt(),0.);
1966 if(Thiscleaned && pass_corrected){
1970 jetME = map_of_MEs[
DirName+
"/"+
"jetReco"];
if(jetME && jetME->getRootObject()) jetME->Fill(1);
1971 mJetEnergyCorr = map_of_MEs[
DirName+
"/"+
"JetEnergyCorr"];
if(mJetEnergyCorr && mJetEnergyCorr->getRootObject()) mJetEnergyCorr->Fill(correctedJet.
pt()/(*caloJets)[ijet].pt());
1972 mJetEnergyCorrVSEta = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSEta"];
if(mJetEnergyCorrVSEta && mJetEnergyCorrVSEta->getRootObject()) mJetEnergyCorrVSEta->Fill(correctedJet.
eta(),correctedJet.
pt()/(*caloJets)[ijet].pt());
1973 mJetEnergyCorrVSPt = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSPt"];
if(mJetEnergyCorrVSPt && mJetEnergyCorrVSPt->getRootObject()) mJetEnergyCorrVSPt->Fill(correctedJet.
pt(),correctedJet.
pt()/(*caloJets)[ijet].pt());
1976 jetME = map_of_MEs[
DirName+
"/"+
"jetReco"];
if(jetME && jetME->getRootObject()) jetME->Fill(2);
1977 mJetEnergyCorr = map_of_MEs[
DirName+
"/"+
"JetEnergyCorr"];
if(mJetEnergyCorr && mJetEnergyCorr->getRootObject()) mJetEnergyCorr->Fill(correctedJet.
pt()/(*pfJets)[ijet].pt());
1978 mJetEnergyCorrVSEta = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSEta"];
if(mJetEnergyCorrVSEta && mJetEnergyCorrVSEta->getRootObject())mJetEnergyCorrVSEta->Fill(correctedJet.
eta(),correctedJet.
pt()/(*pfJets)[ijet].pt());
1979 mJetEnergyCorrVSPt = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSPt"];
if(mJetEnergyCorrVSPt && mJetEnergyCorrVSPt->getRootObject()) mJetEnergyCorrVSPt->Fill(correctedJet.
pt(),correctedJet.
pt()/(*pfJets)[ijet].pt());
1988 if (jetLoPass_ == 1) {
1989 mPhi_Lo = map_of_MEs[
DirName+
"/"+
"Phi_Lo"];
if (mPhi_Lo && mPhi_Lo->getRootObject()) mPhi_Lo->Fill (correctedJet.
phi());
1990 mPt_Lo = map_of_MEs[
DirName+
"/"+
"Pt_Lo"];
if (mPt_Lo && mPt_Lo->getRootObject()) mPt_Lo->Fill (correctedJet.
pt());
1994 if (jetHiPass_ == 1&& correctedJet.
pt()>100. ) {
1995 mEta_Hi = map_of_MEs[
DirName+
"/"+
"Eta_Hi"];
if (mEta_Hi && mEta_Hi->getRootObject()) mEta_Hi->Fill (correctedJet.
eta());
1996 mPhi_Hi = map_of_MEs[
DirName+
"/"+
"Phi_Hi"];
if (mPhi_Hi && mPhi_Hi->getRootObject()) mPhi_Hi->Fill (correctedJet.
phi());
1997 mPt_Hi = map_of_MEs[
DirName+
"/"+
"Pt_Hi"];
if (mPt_Hi && mPt_Hi->getRootObject()) mPt_Hi->Fill (correctedJet.
pt());
1999 mPt = map_of_MEs[
DirName+
"/"+
"Pt"];
if (mPt && mPt->getRootObject()) mPt->Fill (correctedJet.
pt());
2000 mPt_1 = map_of_MEs[
DirName+
"/"+
"Pt_1"];
if (mPt_1 && mPt_1->getRootObject()) mPt_1->Fill (correctedJet.
pt());
2001 mPt_2 = map_of_MEs[
DirName+
"/"+
"Pt_2"];
if (mPt_2 && mPt_2->getRootObject()) mPt_2->Fill (correctedJet.
pt());
2002 mPt_3 = map_of_MEs[
DirName+
"/"+
"Pt_3"];
if (mPt_3 && mPt_3->getRootObject()) mPt_3->Fill (correctedJet.
pt());
2003 mEta = map_of_MEs[
DirName+
"/"+
"Eta"];
if (mEta && mEta->getRootObject()) mEta->Fill (correctedJet.
eta());
2004 mPhi = map_of_MEs[
DirName+
"/"+
"Phi"];
if (mPhi && mPhi->getRootObject()) mPhi->Fill (correctedJet.
phi());
2006 mPhiVSEta = map_of_MEs[
DirName+
"/"+
"PhiVSEta"];
if (mPhiVSEta && mPhiVSEta->getRootObject()) mPhiVSEta->Fill(correctedJet.
eta(),correctedJet.
phi());
2008 mConstituents = map_of_MEs[
DirName+
"/"+
"Constituents"];
if (mConstituents && mConstituents->getRootObject()) mConstituents->Fill (correctedJet.
nConstituents());
2012 mPt_profile = map_of_MEs[
DirName+
"/"+
"Pt_profile"];
if (mPt_profile && mPt_profile->getRootObject()) mPt_profile ->
Fill(numPV, correctedJet.
pt());
2013 mEta_profile = map_of_MEs[
DirName+
"/"+
"Eta_profile"];
if (mEta_profile && mEta_profile->getRootObject()) mEta_profile ->
Fill(numPV, correctedJet.
eta());
2014 mPhi_profile = map_of_MEs[
DirName+
"/"+
"Phi_profile"];
if (mPhi_profile && mPhi_profile->getRootObject()) mPhi_profile ->
Fill(numPV, correctedJet.
phi());
2016 mConstituents_profile = map_of_MEs[
DirName+
"/"+
"Constituents_profile"];
if (mConstituents_profile && mConstituents_profile->getRootObject()) mConstituents_profile->Fill(numPV, correctedJet.
nConstituents());
2018 if (fabs(correctedJet.
eta()) <= 1.3) {
2019 mPt_Barrel = map_of_MEs[
DirName+
"/"+
"Pt_Barrel"];
if (mPt_Barrel && mPt_Barrel->getRootObject()) mPt_Barrel->Fill (correctedJet.
pt());
2020 mPhi_Barrel = map_of_MEs[
DirName+
"/"+
"Phi_Barrel"];
if (mPhi_Barrel && mPhi_Barrel->getRootObject()) mPhi_Barrel->Fill (correctedJet.
phi());
2023 mConstituents_Barrel = map_of_MEs[
DirName+
"/"+
"Constituents_Barrel"];
if (mConstituents_Barrel && mConstituents_Barrel->getRootObject()) mConstituents_Barrel->Fill(correctedJet.
nConstituents());
2025 }
else if (fabs(correctedJet.
eta()) <= 3) {
2026 mPt_EndCap = map_of_MEs[
DirName+
"/"+
"Pt_EndCap"];
if (mPt_EndCap && mPt_EndCap->getRootObject()) mPt_EndCap->Fill (correctedJet.
pt());
2027 mPhi_EndCap = map_of_MEs[
DirName+
"/"+
"Phi_EndCap"];
if (mPhi_EndCap && mPhi_EndCap->getRootObject()) mPhi_EndCap->Fill (correctedJet.
phi());
2030 mConstituents_EndCap = map_of_MEs[
DirName+
"/"+
"Constituents_EndCap"];
if (mConstituents_EndCap && mConstituents_EndCap->getRootObject()) mConstituents_EndCap->Fill(correctedJet.
nConstituents());
2033 mPt_Forward = map_of_MEs[
DirName+
"/"+
"Pt_Forward"];
if (mPt_Forward && mPt_Forward->getRootObject()) mPt_Forward->Fill (correctedJet.
pt());
2034 mPhi_Forward = map_of_MEs[
DirName+
"/"+
"Phi_Forward"];
if (mPhi_Forward && mPhi_Forward->getRootObject()) mPhi_Forward->Fill (correctedJet.
phi());
2037 mConstituents_Forward = map_of_MEs[
DirName+
"/"+
"Constituents_Forward"];
if (mConstituents_Forward && mConstituents_Forward->getRootObject()) mConstituents_Forward->Fill(correctedJet.
nConstituents());
2044 mNJets = map_of_MEs[
DirName+
"/"+
"NJets"];
if (mNJets && mNJets->getRootObject()) mNJets->Fill (numofjets);
2045 mNJets_profile = map_of_MEs[
DirName+
"/"+
"NJets_profile"];
if (mNJets_profile && mNJets_profile->getRootObject()) mNJets_profile->Fill(numPV, numofjets);
2047 sort(recoJets.begin(),recoJets.end(),jetSortingRule);
2053 if(cleaned_first_jet){
2054 mEtaFirst = map_of_MEs[
DirName+
"/"+
"EtaFirst"];
if (mEtaFirst && mEtaFirst->getRootObject()) mEtaFirst->Fill ((recoJets)[0].eta());
2055 mPhiFirst = map_of_MEs[
DirName+
"/"+
"PhiFirst"];
if (mPhiFirst && mPhiFirst->getRootObject()) mPhiFirst->Fill ((recoJets)[0].phi());
2056 mPtFirst = map_of_MEs[
DirName+
"/"+
"PtFirst"];
if (mPtFirst && mPtFirst->getRootObject()) mPtFirst->Fill ((recoJets)[0].pt());
2058 if(numofjets>1 && cleaned_second_jet) {
2059 double dphi=fabs((recoJets)[0].
phi()-(recoJets)[1].
phi());
2061 dphi=2*acos(-1.)-dphi;
2063 mDPhi = map_of_MEs[
DirName+
"/"+
"DPhi"];
if (mDPhi && mDPhi->getRootObject()) mDPhi->Fill (dphi);
2074 if(jetCleaningFlag_ && recoJets.size()>1 && cleaned_first_jet && cleaned_second_jet ){
2076 if(jetCleaningFlag_){
2077 DirName =
"JetMET/Jet/Cleaned"+mInputCollection_.label()+
"/DiJet";
2082 double dphi=fabs((recoJets)[0].
phi()-(recoJets)[1].
phi());
2084 dphi=2*acos(-1.)-dphi;
2086 mDPhi = map_of_MEs[
DirName+
"/"+
"DPhi"];
if (mDPhi && mDPhi->getRootObject())
if (mDPhi) mDPhi->Fill (dphi);
2092 reco::JetID jetID1 = (*jetID_ValueMap_Handle)[calojetref1];
2094 reco::JetID jetID2 = (*jetID_ValueMap_Handle)[calojetref2];
2095 mN90Hits = map_of_MEs[
DirName+
"/"+
"N90Hits"];
if (mN90Hits && mN90Hits->getRootObject()) mN90Hits->Fill (jetID1.
n90Hits);
2096 mfHPD = map_of_MEs[
DirName+
"/"+
"fHPD"];
if (mfHPD && mfHPD->getRootObject()) mfHPD->Fill (jetID1.
fHPD);
2097 mresEMF = map_of_MEs[
DirName+
"/"+
"resEMF"];
if (mresEMF && mresEMF->getRootObject()) mresEMF->Fill (jetID1.
restrictedEMF);
2098 mfRBX = map_of_MEs[
DirName+
"/"+
"fRBX"];
if (mfRBX && mfRBX->getRootObject()) mfRBX->Fill (jetID1.
fRBX);
2099 mN90Hits = map_of_MEs[
DirName+
"/"+
"N90Hits"];
if (mN90Hits && mN90Hits->getRootObject()) mN90Hits->Fill (jetID2.
n90Hits);
2100 mfHPD = map_of_MEs[
DirName+
"/"+
"fHPD"];
if (mfHPD && mfHPD->getRootObject()) mfHPD->Fill (jetID2.
fHPD);
2101 mresEMF = map_of_MEs[
DirName+
"/"+
"resEMF"];
if (mresEMF && mresEMF->getRootObject()) mresEMF->Fill (jetID2.
restrictedEMF);
2102 mfRBX = map_of_MEs[
DirName+
"/"+
"fRBX"];
if (mfRBX && mfRBX->getRootObject()) mfRBX->Fill (jetID2.
fRBX);
2104 mHFrac = map_of_MEs[
DirName+
"/"+
"HFrac"];
if (mHFrac && mHFrac->getRootObject()) mHFrac->Fill ((*caloJets)[ind1].energyFractionHadronic());
2105 mEFrac = map_of_MEs[
DirName+
"/"+
"EFrac"];
if (mEFrac && mHFrac->getRootObject()) mEFrac->Fill ((*caloJets)[ind1].emEnergyFraction());
2106 mMaxEInEmTowers = map_of_MEs[
DirName+
"/"+
"MaxEInEmTowers"];
if (mMaxEInEmTowers && mMaxEInEmTowers->getRootObject()) mMaxEInEmTowers->Fill ((*caloJets)[ind1].maxEInEmTowers());
2107 mMaxEInHadTowers = map_of_MEs[
DirName+
"/"+
"MaxEInHadTowers"];
if (mMaxEInHadTowers && mMaxEInHadTowers->getRootObject()) mMaxEInHadTowers->Fill ((*caloJets)[ind1].maxEInHadTowers());
2108 mHFrac = map_of_MEs[
DirName+
"/"+
"HFrac"];
if (mHFrac && mHFrac->getRootObject()) mHFrac->Fill ((*caloJets)[ind2].energyFractionHadronic());
2109 mEFrac = map_of_MEs[
DirName+
"/"+
"EFrac"];
if (mEFrac && mHFrac->getRootObject()) mEFrac->Fill ((*caloJets)[ind2].emEnergyFraction());
2110 mMaxEInEmTowers = map_of_MEs[
DirName+
"/"+
"MaxEInEmTowers"];
if (mMaxEInEmTowers && mMaxEInEmTowers->getRootObject()) mMaxEInEmTowers->Fill ((*caloJets)[ind2].maxEInEmTowers());
2111 mMaxEInHadTowers = map_of_MEs[
DirName+
"/"+
"MaxEInHadTowers"];
if (mMaxEInHadTowers && mMaxEInHadTowers->getRootObject()) mMaxEInHadTowers->Fill ((*caloJets)[ind2].maxEInHadTowers());
2112 mHFrac_profile = map_of_MEs[
DirName+
"/"+
"HFrac_profile"];
if (mHFrac_profile && mHFrac_profile->getRootObject()) mHFrac_profile ->
Fill(numPV, (*caloJets)[ind1].energyFractionHadronic());
2113 mEFrac_profile = map_of_MEs[
DirName+
"/"+
"EFrac_profile"];
if (mEFrac_profile && mEFrac_profile->getRootObject()) mEFrac_profile ->
Fill(numPV, (*caloJets)[ind1].emEnergyFraction());
2114 mHFrac_profile = map_of_MEs[
DirName+
"/"+
"HFrac_profile"];
if (mHFrac_profile && mHFrac_profile->getRootObject()) mHFrac_profile ->
Fill(numPV, (*caloJets)[ind2].energyFractionHadronic());
2115 mEFrac_profile = map_of_MEs[
DirName+
"/"+
"EFrac_profile"];
if (mEFrac_profile && mEFrac_profile->getRootObject()) mEFrac_profile ->
Fill(numPV, (*caloJets)[ind2].emEnergyFraction());
2117 mJetEnergyCorr = map_of_MEs[
DirName+
"/"+
"JetEnergyCorr"];
if(mJetEnergyCorr && mJetEnergyCorr->getRootObject()) mJetEnergyCorr->Fill(recoJets[0].pt()/(*caloJets)[ind1].pt());
2118 mJetEnergyCorrVSEta = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSEta"];
if(mJetEnergyCorrVSEta && mJetEnergyCorrVSEta->getRootObject()) mJetEnergyCorrVSEta->Fill(recoJets[0].eta(),recoJets[0].pt()/(*caloJets)[ind1].pt());
2119 mJetEnergyCorrVSPt = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSPt"];
if(mJetEnergyCorrVSPt && mJetEnergyCorrVSPt->getRootObject()) mJetEnergyCorrVSPt->Fill(recoJets[0].pt(),recoJets[0].pt()/(*caloJets)[ind1].pt());
2120 mJetEnergyCorr = map_of_MEs[
DirName+
"/"+
"JetEnergyCorr"];
if(mJetEnergyCorr && mJetEnergyCorr->getRootObject()) mJetEnergyCorr->Fill(recoJets[1].pt()/(*caloJets)[ind2].pt());
2121 mJetEnergyCorrVSEta = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSEta"];
if(mJetEnergyCorrVSEta && mJetEnergyCorrVSEta->getRootObject()) mJetEnergyCorrVSEta->Fill(recoJets[1].eta(),recoJets[1].pt()/(*caloJets)[ind2].pt());
2122 mJetEnergyCorrVSPt = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSPt"];
if(mJetEnergyCorrVSPt && mJetEnergyCorrVSPt->getRootObject()) mJetEnergyCorrVSPt->Fill(recoJets[1].pt(),recoJets[1].pt()/(*caloJets)[ind2].pt());
2142 mHFrac = map_of_MEs[
DirName+
"/"+
"HFrac"];
if (mHFrac && mHFrac->getRootObject()) mHFrac->Fill ((*pfJets)[ind1].chargedHadronEnergyFraction()+(*pfJets)[ind1].neutralHadronEnergyFraction()+(*pfJets)[ind1].HFHadronEnergyFraction());
2143 mEFrac = map_of_MEs[
DirName+
"/"+
"EFrac"];
if (mEFrac && mHFrac->getRootObject()) mEFrac->Fill ((*pfJets)[ind1].chargedEmEnergyFraction() +(*pfJets)[ind1].neutralEmEnergyFraction()+(*pfJets)[ind1].HFEMEnergyFraction());
2145 mCHFrac = map_of_MEs[
DirName+
"/"+
"CHFrac"];
if (mCHFrac && mCHFrac->getRootObject()) mCHFrac ->
Fill((*pfJets)[ind1].chargedHadronEnergyFraction());
2146 mNHFrac = map_of_MEs[
DirName+
"/"+
"NHFrac"];
if (mNHFrac && mNHFrac->getRootObject()) mNHFrac ->
Fill((*pfJets)[ind1].neutralHadronEnergyFraction());
2147 mPhFrac = map_of_MEs[
DirName+
"/"+
"PhFrac"];
if (mPhFrac && mPhFrac->getRootObject()) mPhFrac ->
Fill((*pfJets)[ind1].neutralEmEnergyFraction());
2148 mElFrac = map_of_MEs[
DirName+
"/"+
"ElFrac"];
if (mElFrac && mElFrac->getRootObject()) mElFrac ->
Fill((*pfJets)[ind1].chargedEmEnergyFraction());
2149 mMuFrac = map_of_MEs[
DirName+
"/"+
"MuFrac"];
if (mMuFrac && mMuFrac->getRootObject()) mMuFrac ->
Fill((*pfJets)[ind1].chargedMuEnergyFraction());
2150 mHFEMFrac = map_of_MEs[
DirName+
"/"+
"HFEMFrac"];
if (mHFEMFrac && mHFEMFrac->getRootObject()) mHFEMFrac ->
Fill((*pfJets)[ind1].HFEMEnergyFraction());
2151 mHFHFrac = map_of_MEs[
DirName+
"/"+
"HFHFrac"];
if (mHFHFrac && mHFHFrac->getRootObject()) mHFHFrac ->
Fill((*pfJets)[ind1].HFHadronEnergyFraction());
2153 mJetEnergyCorr = map_of_MEs[
DirName+
"/"+
"JetEnergyCorr"];
if(mJetEnergyCorr && mJetEnergyCorr->getRootObject()) mJetEnergyCorr->Fill(recoJets[0].pt()/(*pfJets)[ind1].pt());
2154 mJetEnergyCorrVSEta = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSEta"];
if(mJetEnergyCorrVSEta && mJetEnergyCorrVSEta->getRootObject()) mJetEnergyCorrVSEta->Fill(recoJets[0].eta(),recoJets[0].pt()/(*pfJets)[ind1].pt());
2155 mJetEnergyCorrVSPt = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSPt"];
if(mJetEnergyCorrVSPt && mJetEnergyCorrVSPt->getRootObject()) mJetEnergyCorrVSPt->Fill(recoJets[0].pt(),recoJets[0].pt()/(*pfJets)[ind1].pt());
2156 mJetEnergyCorr = map_of_MEs[
DirName+
"/"+
"JetEnergyCorr"];
if(mJetEnergyCorr && mJetEnergyCorr->getRootObject()) mJetEnergyCorr->Fill(recoJets[1].pt()/(*pfJets)[ind2].pt());
2157 mJetEnergyCorrVSEta = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSEta"];
if(mJetEnergyCorrVSEta && mJetEnergyCorrVSEta->getRootObject()) mJetEnergyCorrVSEta->Fill(recoJets[1].eta(),recoJets[1].pt()/(*pfJets)[ind2].pt());
2158 mJetEnergyCorrVSPt = map_of_MEs[
DirName+
"/"+
"JetEnergyCorrVSPt"];
if(mJetEnergyCorrVSPt && mJetEnergyCorrVSPt->getRootObject()) mJetEnergyCorrVSPt->Fill(recoJets[1].pt(),recoJets[1].pt()/(*pfJets)[ind2].pt());
2160 mChargedMultiplicity = map_of_MEs[
DirName+
"/"+
"ChargedMultiplicity"];
if(mChargedMultiplicity && mChargedMultiplicity->getRootObject()) mChargedMultiplicity->Fill((*pfJets)[ind1].chargedMultiplicity());
2161 mNeutralMultiplicity = map_of_MEs[
DirName+
"/"+
"NeutralMultiplicity"];
if(mNeutralMultiplicity && mNeutralMultiplicity->getRootObject()) mNeutralMultiplicity->Fill((*pfJets)[ind1].neutralMultiplicity());
2162 mMuonMultiplicity = map_of_MEs[
DirName+
"/"+
"MuonMultiplicity"];
if(mMuonMultiplicity && mMuonMultiplicity->getRootObject()) mMuonMultiplicity->Fill((*pfJets)[ind1].muonMultiplicity());
2164 mHFrac = map_of_MEs[
DirName+
"/"+
"HFrac"];
if (mHFrac && mHFrac->getRootObject()) mHFrac->Fill ((*pfJets)[ind2].chargedHadronEnergyFraction()+(*pfJets)[ind2].neutralHadronEnergyFraction()+(*pfJets)[ind2].HFHadronEnergyFraction());
2165 mEFrac = map_of_MEs[
DirName+
"/"+
"EFrac"];
if (mEFrac && mHFrac->getRootObject()) mEFrac->Fill ((*pfJets)[ind2].chargedEmEnergyFraction() +(*pfJets)[ind2].neutralEmEnergyFraction()+(*pfJets)[ind2].HFEMEnergyFraction());
2167 mCHFrac = map_of_MEs[
DirName+
"/"+
"CHFrac"];
if (mCHFrac && mCHFrac->getRootObject()) mCHFrac ->
Fill((*pfJets)[ind2].chargedHadronEnergyFraction());
2168 mNHFrac = map_of_MEs[
DirName+
"/"+
"NHFrac"];
if (mNHFrac && mNHFrac->getRootObject()) mNHFrac ->
Fill((*pfJets)[ind2].neutralHadronEnergyFraction());
2169 mPhFrac = map_of_MEs[
DirName+
"/"+
"PhFrac"];
if (mPhFrac && mPhFrac->getRootObject()) mPhFrac ->
Fill((*pfJets)[ind2].neutralEmEnergyFraction());
2170 mElFrac = map_of_MEs[
DirName+
"/"+
"ElFrac"];
if (mElFrac && mElFrac->getRootObject()) mElFrac ->
Fill((*pfJets)[ind2].chargedEmEnergyFraction());
2171 mMuFrac = map_of_MEs[
DirName+
"/"+
"MuFrac"];
if (mMuFrac && mMuFrac->getRootObject()) mMuFrac ->
Fill((*pfJets)[ind2].chargedMuEnergyFraction());
2172 mHFEMFrac = map_of_MEs[
DirName+
"/"+
"HFEMFrac"];
if (mHFEMFrac && mHFEMFrac->getRootObject()) mHFEMFrac ->
Fill((*pfJets)[ind2].HFEMEnergyFraction());
2173 mHFHFrac = map_of_MEs[
DirName+
"/"+
"HFHFrac"];
if (mHFHFrac && mHFHFrac->getRootObject()) mHFHFrac ->
Fill((*pfJets)[ind2].HFHadronEnergyFraction());
2175 mNeutralFraction = map_of_MEs[
DirName+
"/"+
"NeutralConstituentsFraction"];
if (mNeutralFraction && mNeutralFraction->getRootObject()) mNeutralFraction->Fill ((
double)(*pfJets)[ind1].neutralMultiplicity()/(double)(*pfJets)[ind1].nConstituents());
2177 mChargedMultiplicity = map_of_MEs[
DirName+
"/"+
"ChargedMultiplicity"];
if(mChargedMultiplicity && mChargedMultiplicity->getRootObject()) mChargedMultiplicity->Fill((*pfJets)[ind2].chargedMultiplicity());
2178 mNeutralMultiplicity = map_of_MEs[
DirName+
"/"+
"NeutralMultiplicity"];
if(mNeutralMultiplicity && mNeutralMultiplicity->getRootObject()) mNeutralMultiplicity->Fill((*pfJets)[ind2].neutralMultiplicity());
2179 mMuonMultiplicity = map_of_MEs[
DirName+
"/"+
"MuonMultiplicity"];
if(mMuonMultiplicity && mMuonMultiplicity->getRootObject()) mMuonMultiplicity->Fill((*pfJets)[ind2].muonMultiplicity());
2182 mHFrac_profile = map_of_MEs[
DirName+
"/"+
"HFrac_profile"];
if (mHFrac_profile && mHFrac_profile->getRootObject()) mHFrac_profile ->
Fill(numPV, (*pfJets)[ind1].chargedHadronEnergyFraction() + (*pfJets)[ind1].neutralHadronEnergyFraction());
2183 mEFrac_profile = map_of_MEs[
DirName+
"/"+
"EFrac_profile"];
if (mEFrac_profile && mEFrac_profile->getRootObject()) mEFrac_profile ->
Fill(numPV, (*pfJets)[ind1].chargedEmEnergyFraction() + (*pfJets)[ind1].neutralEmEnergyFraction());
2184 mCHFrac_profile = map_of_MEs[
DirName+
"/"+
"CHFrac_profile"];
if (mCHFrac_profile && mCHFrac_profile->getRootObject()) mCHFrac_profile ->
Fill(numPV, (*pfJets)[ind1].chargedHadronEnergyFraction());
2185 mNHFrac_profile = map_of_MEs[
DirName+
"/"+
"NHFrac_profile"];
if (mNHFrac_profile && mNHFrac_profile->getRootObject()) mNHFrac_profile ->
Fill(numPV, (*pfJets)[ind1].neutralHadronEnergyFraction());
2186 mPhFrac_profile = map_of_MEs[
DirName+
"/"+
"PhFrac_profile"];
if (mPhFrac_profile && mPhFrac_profile->getRootObject()) mPhFrac_profile ->
Fill(numPV, (*pfJets)[ind1].neutralEmEnergyFraction());
2187 mElFrac_profile = map_of_MEs[
DirName+
"/"+
"ElFrac_profile"];
if (mElFrac_profile && mElFrac_profile->getRootObject()) mElFrac_profile ->
Fill(numPV, (*pfJets)[ind1].chargedEmEnergyFraction());
2188 mMuFrac_profile = map_of_MEs[
DirName+
"/"+
"MuFrac_profile"];
if (mMuFrac_profile && mMuFrac_profile->getRootObject()) mMuFrac_profile ->
Fill(numPV, (*pfJets)[ind1].chargedMuEnergyFraction());
2189 mHFEMFrac_profile = map_of_MEs[
DirName+
"/"+
"HFEMFrac_profile"];
if (mHFEMFrac_profile && mHFEMFrac_profile->getRootObject()) mHFEMFrac_profile ->
Fill(numPV, (*pfJets)[ind1].HFEMEnergyFraction());
2190 mHFHFrac_profile = map_of_MEs[
DirName+
"/"+
"HFHFrac_profile"];
if (mHFHFrac_profile && mHFHFrac_profile->getRootObject()) mHFHFrac_profile ->
Fill(numPV, (*pfJets)[ind1].HFHadronEnergyFraction());
2192 mNeutralFraction = map_of_MEs[
DirName+
"/"+
"NeutralConstituentsFraction"];
if (mNeutralFraction && mNeutralFraction->getRootObject()) mNeutralFraction->Fill ((
double)(*pfJets)[ind2].neutralMultiplicity()/(double)(*pfJets)[ind2].nConstituents());
2194 mChargedMultiplicity_profile = map_of_MEs[
DirName+
"/"+
"ChargedMultiplicity_profile"];
if(mChargedMultiplicity_profile && mChargedMultiplicity_profile->getRootObject()) mChargedMultiplicity_profile->Fill(numPV, (*pfJets)[ind1].chargedMultiplicity());
2195 mNeutralMultiplicity_profile = map_of_MEs[
DirName+
"/"+
"NeutralMultiplicity_profile"];
if(mNeutralMultiplicity_profile && mNeutralMultiplicity_profile->getRootObject()) mNeutralMultiplicity_profile->Fill(numPV, (*pfJets)[ind1].neutralMultiplicity());
2196 mMuonMultiplicity_profile = map_of_MEs[
DirName+
"/"+
"MuonMultiplicity_profile"];
if(mMuonMultiplicity_profile && mMuonMultiplicity_profile->getRootObject()) mMuonMultiplicity->Fill(numPV, (*pfJets)[ind1].muonMultiplicity());
2198 mHFrac_profile = map_of_MEs[
DirName+
"/"+
"HFrac_profile"];
if (mHFrac_profile && mHFrac_profile->getRootObject()) mHFrac_profile ->
Fill(numPV, (*pfJets)[ind2].chargedHadronEnergyFraction() + (*pfJets)[ind1].neutralHadronEnergyFraction());
2199 mEFrac_profile = map_of_MEs[
DirName+
"/"+
"EFrac_profile"];
if (mEFrac_profile && mEFrac_profile->getRootObject()) mEFrac_profile ->
Fill(numPV, (*pfJets)[ind2].chargedEmEnergyFraction() + (*pfJets)[ind1].neutralEmEnergyFraction());
2200 mCHFrac_profile = map_of_MEs[
DirName+
"/"+
"CHFrac_profile"];
if (mCHFrac_profile && mCHFrac_profile->getRootObject()) mCHFrac_profile ->
Fill(numPV, (*pfJets)[ind2].chargedHadronEnergyFraction());
2201 mNHFrac_profile = map_of_MEs[
DirName+
"/"+
"NHFrac_profile"];
if (mNHFrac_profile && mNHFrac_profile->getRootObject()) mNHFrac_profile ->
Fill(numPV, (*pfJets)[ind2].neutralHadronEnergyFraction());
2202 mPhFrac_profile = map_of_MEs[
DirName+
"/"+
"PhFrac_profile"];
if (mPhFrac_profile && mPhFrac_profile->getRootObject()) mPhFrac_profile ->
Fill(numPV, (*pfJets)[ind2].neutralEmEnergyFraction());
2203 mElFrac_profile = map_of_MEs[
DirName+
"/"+
"ElFrac_profile"];
if (mElFrac_profile && mElFrac_profile->getRootObject()) mElFrac_profile ->
Fill(numPV, (*pfJets)[ind2].chargedEmEnergyFraction());
2204 mMuFrac_profile = map_of_MEs[
DirName+
"/"+
"MuFrac_profile"];
if (mMuFrac_profile && mMuFrac_profile->getRootObject()) mMuFrac_profile ->
Fill(numPV, (*pfJets)[ind2].chargedMuEnergyFraction());
2205 mHFEMFrac_profile = map_of_MEs[
DirName+
"/"+
"HFEMFrac_profile"];
if (mHFEMFrac_profile && mHFEMFrac_profile->getRootObject()) mHFEMFrac_profile ->
Fill(numPV, (*pfJets)[ind2].HFEMEnergyFraction());
2206 mHFHFrac_profile = map_of_MEs[
DirName+
"/"+
"HFHFrac_profile"];
if (mHFHFrac_profile && mHFHFrac_profile->getRootObject()) mHFHFrac_profile ->
Fill(numPV, (*pfJets)[ind2].HFHadronEnergyFraction());
2208 mChargedMultiplicity_profile = map_of_MEs[
DirName+
"/"+
"ChargedMultiplicity_profile"];
if(mChargedMultiplicity_profile && mChargedMultiplicity_profile->getRootObject()) mChargedMultiplicity_profile->Fill(numPV, (*pfJets)[ind2].chargedMultiplicity());
2209 mNeutralMultiplicity_profile = map_of_MEs[
DirName+
"/"+
"NeutralMultiplicity_profile"];
if(mNeutralMultiplicity_profile && mNeutralMultiplicity_profile->getRootObject()) mNeutralMultiplicity_profile->Fill(numPV, (*pfJets)[ind2].neutralMultiplicity());
2210 mMuonMultiplicity_profile = map_of_MEs[
DirName+
"/"+
"MuonMultiplicity_profile"];
if(mMuonMultiplicity_profile && mMuonMultiplicity_profile->getRootObject()) mMuonMultiplicity_profile->Fill(numPV, (*pfJets)[ind2].muonMultiplicity());
2213 mPt = map_of_MEs[
DirName+
"/"+
"Pt"];
if (mPt && mPt->getRootObject()) mPt->Fill (recoJets[0].pt());
2214 mEta = map_of_MEs[
DirName+
"/"+
"Eta"];
if (mEta && mEta->getRootObject()) mEta->Fill (recoJets[0].eta());
2215 mPhi = map_of_MEs[
DirName+
"/"+
"Phi"];
if (mPhi && mPhi->getRootObject()) mPhi->Fill (recoJets[0].phi());
2216 mPhiVSEta = map_of_MEs[
DirName+
"/"+
"PhiVSEta"];
if (mPhiVSEta && mPhiVSEta->getRootObject()) mPhiVSEta->Fill(recoJets[0].eta(),recoJets[0].phi());
2218 mConstituents = map_of_MEs[
DirName+
"/"+
"Constituents"];
if (mConstituents && mConstituents->getRootObject()) mConstituents->Fill (recoJets[0].nConstituents());
2220 mPt = map_of_MEs[
DirName+
"/"+
"Pt"];
if (mPt && mPt->getRootObject()) mPt->Fill (recoJets[1].pt());
2221 mEta = map_of_MEs[
DirName+
"/"+
"Eta"];
if (mEta && mEta->getRootObject()) mEta->Fill (recoJets[1].eta());
2222 mPhi = map_of_MEs[
DirName+
"/"+
"Phi"];
if (mPhi && mPhi->getRootObject()) mPhi->Fill (recoJets[1].phi());
2223 mPhiVSEta = map_of_MEs[
DirName+
"/"+
"PhiVSEta"];
if (mPhiVSEta && mPhiVSEta->getRootObject()) mPhiVSEta->Fill(recoJets[1].eta(),recoJets[1].phi());
2225 mConstituents = map_of_MEs[
DirName+
"/"+
"Constituents"];
if (mConstituents && mConstituents->getRootObject()) mConstituents->Fill (recoJets[1].nConstituents());
2228 mPt_profile = map_of_MEs[
DirName+
"/"+
"Pt_profile"];
if (mPt_profile && mPt_profile->getRootObject()) mPt_profile ->
Fill(numPV, recoJets[0].
pt());
2229 mEta_profile = map_of_MEs[
DirName+
"/"+
"Eta_profile"];
if (mEta_profile && mEta_profile->getRootObject()) mEta_profile ->
Fill(numPV, recoJets[0].
eta());
2230 mPhi_profile = map_of_MEs[
DirName+
"/"+
"Phi_profile"];
if (mPhi_profile && mPhi_profile->getRootObject()) mPhi_profile ->
Fill(numPV, recoJets[0].
phi());
2232 mConstituents_profile = map_of_MEs[
DirName+
"/"+
"Constituents_profile"];
if (mConstituents_profile && mConstituents_profile->getRootObject()) mConstituents_profile->Fill(numPV, recoJets[0].nConstituents());
2234 mPt_profile = map_of_MEs[
DirName+
"/"+
"Pt_profile"];
if (mPt_profile && mPt_profile->getRootObject()) mPt_profile ->
Fill(numPV, recoJets[1].
pt());
2235 mEta_profile = map_of_MEs[
DirName+
"/"+
"Eta_profile"];
if (mEta_profile && mEta_profile->getRootObject()) mEta_profile ->
Fill(numPV, recoJets[1].
eta());
2236 mPhi_profile = map_of_MEs[
DirName+
"/"+
"Phi_profile"];
if (mPhi_profile && mPhi_profile->getRootObject()) mPhi_profile ->
Fill(numPV, recoJets[1].
phi());
2238 mConstituents_profile = map_of_MEs[
DirName+
"/"+
"Constituents_profile"];
if (mConstituents_profile && mConstituents_profile->getRootObject()) mConstituents_profile->Fill(numPV, recoJets[1].nConstituents());
2240 if (fabs(recoJets[0].
eta() < 1.4)) {
2241 double pt_dijet = (recoJets[0].pt() + recoJets[1].pt())/2;
2247 int randJet =iEvent.
id().
event()%2;
2248 if (fabs(recoJets[1].
eta() < 1.4)) {
2264 bool thirdJetCut =
true;
2272 if (pt3 > asymmetryThirdJetCut_){
2277 double dijetAsymmetry =(recoJets[jet1].pt() - recoJets[jet2].pt()) / (recoJets[jet1].
pt() + recoJets[jet2].pt());
2278 mDijetAsymmetry = map_of_MEs[
DirName+
"/"+
"DijetAsymmetry"];
if (mDijetAsymmetry && mDijetAsymmetry->getRootObject()) mDijetAsymmetry->Fill(dijetAsymmetry);
2287 pt_barrel = recoJets[jet1].pt();
2288 pt_probe = recoJets[jet2].pt();
2298 bool thirdJetCut =
true;
2300 if (pt3/pt_dijet > balanceThirdJetCut_){
2301 thirdJetCut =
false;
2305 double dijetBalance = (pt_probe - pt_barrel) / pt_dijet;
2306 mDijetBalance = map_of_MEs[
DirName+
"/"+
"DijetBalance"];
if (mDijetBalance && mDijetBalance->getRootObject()) mDijetBalance->Fill(dijetBalance);
T getParameter(std::string const &) const
EventNumber_t event() const
T getUntrackedParameter(std::string const &, T const &) const
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
virtual void scaleEnergy(double fScale)
scale energy of the jet
virtual float pt() const
transverse momentum
JetAnalyzer(const edm::ParameterSet &)
Constructor.
MonitorElement * bookProfile(Args &&...args)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
const std::string metname
virtual float phi() const
momentum azimuthal angle
Base class for all types of Jets.
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
std::vector< Vertex > VertexCollection
collection of Vertex objects
double correction(const LorentzVector &fJet) const
get correction using Jet information only
Provides a code based selection for trigger and DCS information in order to have no failing filters i...
Strings::size_type size() const
triggerResultsToken_(consumes< edm::TriggerResults >(edm::InputTag("TriggerResults")))
void endRun(const edm::Run &, const edm::EventSetup &)
Finish up a run.
virtual ~JetAnalyzer()
Destructor.
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
virtual float eta() const
momentum pseudorapidity
MonitorElement * book1D(Args &&...args)
PF Jet selector for pat::Jets.
Jet selector for pat::Jets and for CaloJets.
void setCurrentFolder(const std::string &fullpath)
void dqmBeginRun(const edm::Run &, const edm::EventSetup &)
Initialize run-based parameters.
T const * product() const
MonitorElement * book2D(Args &&...args)
std::string const & triggerName(unsigned int index) const
static std::string const triggerResults("TriggerResults")
static const JetCorrector * getJetCorrector(const std::string &fName, const edm::EventSetup &fSetup)
retrieve corrector from the event setup. troughs exception if something is missing ...
void analyze(const edm::Event &, const edm::EventSetup &)
Get the analysis.
virtual int nConstituents() const
of constituents
void bookMESetSelection(std::string, DQMStore::IBooker &)
TH2F * getTH2F(void) const
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Inizialize parameters for histo binning.
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)