27 jetCollectionTag_ (conf.getUntrackedParameter<edm::InputTag>(
"jetCollectionTag")),
29 widejetsCollectionTag_ (conf.getUntrackedParameter<edm::InputTag>(
"widejetsCollectionTag")),
30 metCollectionTag_ (conf.getUntrackedParameter<edm::InputTag>(
"metCollectionTag")),
31 metCleanCollectionTag_ (conf.getUntrackedParameter<edm::InputTag>(
"metCleanCollectionTag")),
32 numwidejets_ (conf.getParameter<unsigned int>(
"numwidejets")),
33 etawidejets_ (conf.getParameter<double>(
"etawidejets")),
34 ptwidejets_ (conf.getParameter<double>(
"ptwidejets")),
35 detawidejets_ (conf.getParameter<double>(
"detawidejets")),
36 dphiwidejets_ (conf.getParameter<double>(
"dphiwidejets")),
37 maxEMfraction_ (conf.getParameter<double>(
"maxEMfraction")),
38 maxHADfraction_ (conf.getParameter<double>(
"maxHADfraction")),
39 HLTpathMain_ (triggerExpression::
parse( conf.getParameter<std::
string>(
"HLTpathMain") )),
40 HLTpathMonitor_ (triggerExpression::
parse( conf.getParameter<std::
string>(
"HLTpathMonitor") )),
70 for(reco::CaloJetCollection::const_iterator it = calojets_handle->begin(); it != calojets_handle->end(); ++it)
92 TLorentzVector wdijet;
95 double DeltaEtaJJWide = -1;
96 double DeltaPhiJJWide = -1;
98 if( widejets_handle->size() >= 2 )
100 wj1.SetPtEtaPhiM(widejets_handle->at(0).pt(),
101 widejets_handle->at(0).eta(),
102 widejets_handle->at(0).phi(),
103 widejets_handle->at(0).mass()
105 wj2.SetPtEtaPhiM(widejets_handle->at(1).pt(),
106 widejets_handle->at(1).eta(),
107 widejets_handle->at(1).phi(),
108 widejets_handle->at(1).mass()
113 MJJWide = wdijet.M();
114 DeltaEtaJJWide = fabs(wj1.Eta()-wj2.Eta());
115 DeltaPhiJJWide = fabs(wj1.DeltaPhi(wj2));
134 if( calomet_handle.isValid() && calometClean_handle.isValid() )
142 m_metDiff->
Fill( (calomet_handle->front()).pt() - (calometClean_handle->front()).pt() );
143 m_metVSmetclean->
Fill( (calomet_handle->front()).pt() , (calometClean_handle->front()).pt() );
145 else if( calomet_handle.isValid() && !calometClean_handle.isValid() )
150 else if( !calomet_handle.isValid() && calometClean_handle.isValid() )
156 bool pass_nocut=
false;
157 bool pass_twowidejets=
false;
158 bool pass_etaptcuts=
false;
159 bool pass_deta=
false;
160 bool pass_JetIDtwojets=
true;
161 bool pass_dphi=
false;
162 bool pass_metFilter=
true;
164 bool pass_deta_L4=
false;
165 bool pass_deta_L3=
false;
166 bool pass_deta_L2=
false;
168 bool pass_fullsel_NOdeta=
false;
169 bool pass_fullsel_detaL4=
false;
170 bool pass_fullsel_detaL3=
false;
171 bool pass_fullsel_detaL2=
false;
172 bool pass_fullsel=
false;
181 pass_twowidejets=
true;
201 if( DeltaEtaJJWide < 4.0 )
204 if( DeltaEtaJJWide < 3.0 )
207 if( DeltaEtaJJWide < 2.0 )
217 reco::CaloJetCollection::const_iterator thisJet = calojets_handle->begin();
220 pass_JetIDtwojets=
false;
226 pass_JetIDtwojets=
false;
230 if( calomet_handle.isValid() && calometClean_handle.isValid() )
232 if( fabs ( (calomet_handle->front()).pt() - (calometClean_handle->front()).pt() ) > 0.1 )
233 pass_metFilter=
false;
237 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_JetIDtwojets && pass_dphi && pass_metFilter )
238 pass_fullsel_NOdeta=
true;
241 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_JetIDtwojets && pass_dphi && pass_metFilter && pass_deta_L4 )
242 pass_fullsel_detaL4=
true;
243 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_JetIDtwojets && pass_dphi && pass_metFilter && pass_deta_L3 )
244 pass_fullsel_detaL3=
true;
245 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_JetIDtwojets && pass_dphi && pass_metFilter && pass_deta_L2 )
246 pass_fullsel_detaL2=
true;
249 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_deta && pass_JetIDtwojets && pass_dphi && pass_metFilter )
257 if( pass_nocut && pass_twowidejets )
259 if( pass_nocut && pass_twowidejets && pass_etaptcuts )
261 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_deta )
263 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_deta && pass_JetIDtwojets )
265 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_deta && pass_JetIDtwojets && pass_dphi )
282 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_deta )
289 if( pass_fullsel_NOdeta )
294 if( DeltaEtaJJWide >= 0.0 && DeltaEtaJJWide < 0.5 )
296 if( DeltaEtaJJWide >= 0.5 && DeltaEtaJJWide < 1.0 )
298 if( DeltaEtaJJWide >= 1.0 && DeltaEtaJJWide < 1.5 )
300 if( DeltaEtaJJWide >= 1.5 && DeltaEtaJJWide < 2.0 )
302 if( DeltaEtaJJWide >= 2.0 && DeltaEtaJJWide < 2.5 )
304 if( DeltaEtaJJWide >= 2.5 && DeltaEtaJJWide < 3.0 )
306 if( DeltaEtaJJWide >= 3.0 )
323 int HLTpathMain_fired = -1;
324 int HLTpathMonitor_fired = -1;
351 if( pass_fullsel_NOdeta && HLTpathMonitor_fired == 1 )
356 if( HLTpathMain_fired == 1)
363 if( pass_fullsel_detaL4 && HLTpathMonitor_fired == 1 )
368 if( HLTpathMain_fired == 1)
375 if( pass_fullsel_detaL3 && HLTpathMonitor_fired == 1 )
380 if( HLTpathMain_fired == 1)
387 if( pass_fullsel_detaL2 && HLTpathMonitor_fired == 1 )
392 if( HLTpathMain_fired == 1)
399 if( pass_fullsel && HLTpathMonitor_fired == 1 )
404 if( HLTpathMain_fired == 1)
421 const int N_mass_bins=83;
422 float massBins[N_mass_bins+1] = {1, 3, 6, 10, 16, 23, 31, 40, 50, 61, 74, 88, 103, 119, 137, 156, 176, 197, 220, 244, 270, 296, 325, 354, 386, 419, 453, 489, 526, 565, 606, 649, 693, 740, 788, 838, 890, 944, 1000, 1058, 1118, 1181, 1246, 1313, 1383, 1455, 1530, 1607, 1687, 1770, 1856, 1945, 2037, 2132, 2231, 2332, 2438, 2546, 2659, 2775, 2895, 3019, 3147, 3279, 3416, 3558, 3704, 3854, 4010, 4171, 4337, 4509, 4686, 4869, 5058, 5253, 5455, 5663, 5877, 6099, 6328, 6564, 6808, 7000};
440 "M_{jj} WideJets (final selection)",
442 "M_{jj} WideJets [GeV]",
447 "M_{jj} WideJets (final selection)",
448 N_mass_bins, massBins,
449 "M_{jj} WideJets [GeV]",
454 "M_{jj} WideJets (final selection, without noise filters)",
456 "M_{jj} WideJets [GeV]",
461 "M_{jj} WideJets (final selection, without noise filters)",
462 N_mass_bins, massBins,
463 "M_{jj} WideJets [GeV]",
468 "M_{jj} WideJets (0.0<=#Delta#eta<0.5)",
470 "M_{jj} WideJets [GeV]",
475 "M_{jj} WideJets (0.5<=#Delta#eta<1.0)",
477 "M_{jj} WideJets [GeV]",
482 "M_{jj} WideJets (1.0<=#Delta#eta<1.5)",
484 "M_{jj} WideJets [GeV]",
489 "M_{jj} WideJets (1.5<=#Delta#eta<2.0)",
491 "M_{jj} WideJets [GeV]",
496 "M_{jj} WideJets (2.0<=#Delta#eta<2.5)",
498 "M_{jj} WideJets [GeV]",
503 "M_{jj} WideJets (2.5<=#Delta#eta<3.0)",
505 "M_{jj} WideJets [GeV]",
510 "M_{jj} WideJets (#Delta#eta>=3.0)",
512 "M_{jj} WideJets [GeV]",
517 "HLT Efficiency Studies (no deta cut)",
519 "M_{jj} WideJets [GeV]",
524 "HLT Efficiency Studies (no deta cut)",
526 "M_{jj} WideJets [GeV]",
531 "HLT Efficiency Studies (deta cut < 4.0)",
533 "M_{jj} WideJets [GeV]",
538 "HLT Efficiency Studies (deta cut < 4.0)",
540 "M_{jj} WideJets [GeV]",
545 "HLT Efficiency Studies (deta cut < 3.0)",
547 "M_{jj} WideJets [GeV]",
552 "HLT Efficiency Studies (deta cut < 3.0)",
554 "M_{jj} WideJets [GeV]",
559 "HLT Efficiency Studies (deta cut < 2.0)",
561 "M_{jj} WideJets [GeV]",
566 "HLT Efficiency Studies (deta cut < 2.0)",
568 "M_{jj} WideJets [GeV]",
573 "HLT Efficiency Studies (default deta cut)",
575 "M_{jj} WideJets [GeV]",
580 "HLT Efficiency Studies (default deta cut)",
582 "M_{jj} WideJets [GeV]",
587 "#Delta#eta_{jj} WideJets (final selection)",
589 "#Delta#eta_{jj} WideJets",
594 "#Delta#eta_{jj} WideJets (final selection except #Delta#eta cut)",
596 "#Delta#eta_{jj} WideJets",
601 "#Delta#phi_{jj} WideJets (final selection)",
603 "#Delta#phi_{jj} WideJets [rad.]",
632 "HAD Energy Fraction",
639 "EM Energy Fraction",
653 "met - metcleaned [GeV]",
683 "HT (final selection)",
692 "#Delta#eta_{jj} vs M_{jj} WideJets",
695 "M_{jj} WideJets [GeV]",
696 "#Delta#eta_{jj} WideJets");
699 "#Delta#eta_{jj} vs M_{jj} WideJets",
702 "M_{jj} WideJets [GeV]",
703 "#Delta#eta_{jj} WideJets");
MonitorElement * m_HT_finalSel
MonitorElement * m_MjjWide_num_NOdeta
MonitorElement * m_MjjWide_num_detaL4
MonitorElement * m_MjjWide_den_detaL3
edm::InputTag metCleanCollectionTag_
virtual void endRun(edm::Run const &, edm::EventSetup const &)
MonitorElement * bookH2withSumw2(const std::string &name, const std::string &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="COLZ")
MonitorElement * m_HT_inclusive
MonitorElement * m_MjjWide_den_NOdeta
bool configurationUpdated() const
MonitorElement * m_MjjWide_num_detaL2
unsigned int numwidejets_
Evaluator * parse(const T &text)
MonitorElement * m_selJets_towersArea
MonitorElement * bookH1withSumw2(const std::string &name, const std::string &title, int nchX, double lowX, double highX, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")
edm::InputTag metCollectionTag_
MonitorElement * m_metCaseNoMetClean
ParameterSet const & getParameterSet(ParameterSetID const &id)
edm::InputTag jetCollectionTag_
MonitorElement * m_MjjWide_deta_0p0_0p5
triggerExpression::Evaluator * HLTpathMonitor_
triggerExpression::Data triggerConfiguration_
MonitorElement * m_metDiff
MonitorElement * m_MjjWide_num_detaL3
MonitorElement * m_MjjWide_deta_3p0_inf
MonitorElement * m_MjjWide_deta_1p0_1p5
MonitorElement * m_selJets_emEnergyFraction
DiJetVarAnalyzer(const edm::ParameterSet &)
MonitorElement * m_selJets_phi
MonitorElement * bookH1withSumw2BinArray(const std::string &name, const std::string &title, int nchX, float *xbinsize, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")
MonitorElement * m_MjjWide_num
virtual void analyze(const edm::Event &, const edm::EventSetup &)
MonitorElement * m_selJets_pt
virtual void beginRun(const edm::Run &, const edm::EventSetup &)
MonitorElement * m_MjjWide_finalSel_varbin
MonitorElement * m_cutFlow
virtual void init(const Data &data)
MonitorElement * m_MjjWide_finalSel_WithoutNoiseFilter
MonitorElement * m_metVSmetclean
MonitorElement * m_MjjWide_den
edm::InputTag widejetsCollectionTag_
triggerExpression::Evaluator * HLTpathMain_
MonitorElement * m_DphijjWide_finalSel
MonitorElement * m_MjjWide_deta_0p5_1p0
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
MonitorElement * m_MjjWide_deta_2p0_2p5
MonitorElement * m_MjjWide_den_detaL2
MonitorElement * m_DetajjVsMjjWide_rebin
MonitorElement * m_selJets_eta
MonitorElement * m_MjjWide_deta_2p5_3p0
MonitorElement * m_DetajjVsMjjWide
bool setEvent(const edm::Event &event, const edm::EventSetup &setup)
MonitorElement * m_MjjWide_deta_1p5_2p0
MonitorElement * m_MjjWide_finalSel_WithoutNoiseFilter_varbin
MonitorElement * m_DetajjWide_finalSel
MonitorElement * m_MjjWide_finalSel
MonitorElement * m_MjjWide_den_detaL4
MonitorElement * m_DetajjWide
virtual ~DiJetVarAnalyzer()
MonitorElement * m_selJets_hadEnergyFraction
MonitorElement * m_metCases