76 for(reco::CaloJetCollection::const_iterator it = calojets_handle->begin(); it != calojets_handle->end(); ++it)
98 TLorentzVector wdijet;
101 double DeltaEtaJJWide = -1;
102 double DeltaPhiJJWide = -1;
104 if( widejets_handle->size() >= 2 )
106 wj1.SetPtEtaPhiM(widejets_handle->at(0).pt(),
107 widejets_handle->at(0).eta(),
108 widejets_handle->at(0).phi(),
109 widejets_handle->at(0).mass()
111 wj2.SetPtEtaPhiM(widejets_handle->at(1).pt(),
112 widejets_handle->at(1).eta(),
113 widejets_handle->at(1).phi(),
114 widejets_handle->at(1).mass()
119 MJJWide = wdijet.M();
120 DeltaEtaJJWide = fabs(wj1.Eta()-wj2.Eta());
121 DeltaPhiJJWide = fabs(wj1.DeltaPhi(wj2));
140 if( calomet_handle.isValid() && calometClean_handle.isValid() )
148 m_metDiff->
Fill( (calomet_handle->front()).
pt() - (calometClean_handle->front()).
pt() );
151 else if( calomet_handle.isValid() && !calometClean_handle.isValid() )
156 else if( !calomet_handle.isValid() && calometClean_handle.isValid() )
162 bool pass_nocut=
false;
163 bool pass_twowidejets=
false;
164 bool pass_etaptcuts=
false;
165 bool pass_deta=
false;
166 bool pass_JetIDtwojets=
true;
167 bool pass_dphi=
false;
168 bool pass_metFilter=
true;
170 bool pass_deta_L4=
false;
171 bool pass_deta_L3=
false;
172 bool pass_deta_L2=
false;
174 bool pass_fullsel_NOdeta=
false;
175 bool pass_fullsel_detaL4=
false;
176 bool pass_fullsel_detaL3=
false;
177 bool pass_fullsel_detaL2=
false;
178 bool pass_fullsel=
false;
187 pass_twowidejets=
true;
207 if( DeltaEtaJJWide < 4.0 )
210 if( DeltaEtaJJWide < 3.0 )
213 if( DeltaEtaJJWide < 2.0 )
223 reco::CaloJetCollection::const_iterator thisJet = calojets_handle->begin();
226 pass_JetIDtwojets=
false;
232 pass_JetIDtwojets=
false;
236 if( calomet_handle.isValid() && calometClean_handle.isValid() )
238 if( fabs ( (calomet_handle->front()).
pt() - (calometClean_handle->front()).
pt() ) > 0.1 )
239 pass_metFilter=
false;
243 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_JetIDtwojets && pass_dphi && pass_metFilter )
244 pass_fullsel_NOdeta=
true;
247 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_JetIDtwojets && pass_dphi && pass_metFilter && pass_deta_L4 )
248 pass_fullsel_detaL4=
true;
249 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_JetIDtwojets && pass_dphi && pass_metFilter && pass_deta_L3 )
250 pass_fullsel_detaL3=
true;
251 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_JetIDtwojets && pass_dphi && pass_metFilter && pass_deta_L2 )
252 pass_fullsel_detaL2=
true;
255 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_deta && pass_JetIDtwojets && pass_dphi && pass_metFilter )
263 if( pass_nocut && pass_twowidejets )
265 if( pass_nocut && pass_twowidejets && pass_etaptcuts )
267 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_deta )
269 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_deta && pass_JetIDtwojets )
271 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_deta && pass_JetIDtwojets && pass_dphi )
288 if( pass_nocut && pass_twowidejets && pass_etaptcuts && pass_deta )
295 if( pass_fullsel_NOdeta )
300 if( DeltaEtaJJWide >= 0.0 && DeltaEtaJJWide < 0.5 )
302 if( DeltaEtaJJWide >= 0.5 && DeltaEtaJJWide < 1.0 )
304 if( DeltaEtaJJWide >= 1.0 && DeltaEtaJJWide < 1.5 )
306 if( DeltaEtaJJWide >= 1.5 && DeltaEtaJJWide < 2.0 )
308 if( DeltaEtaJJWide >= 2.0 && DeltaEtaJJWide < 2.5 )
310 if( DeltaEtaJJWide >= 2.5 && DeltaEtaJJWide < 3.0 )
312 if( DeltaEtaJJWide >= 3.0 )
329 int HLTpathMain_fired = -1;
330 int HLTpathMonitor_fired = -1;
357 if( pass_fullsel_NOdeta && HLTpathMonitor_fired == 1 )
362 if( HLTpathMain_fired == 1)
369 if( pass_fullsel_detaL4 && HLTpathMonitor_fired == 1 )
374 if( HLTpathMain_fired == 1)
381 if( pass_fullsel_detaL3 && HLTpathMonitor_fired == 1 )
386 if( HLTpathMain_fired == 1)
393 if( pass_fullsel_detaL2 && HLTpathMonitor_fired == 1 )
398 if( HLTpathMain_fired == 1)
405 if( pass_fullsel && HLTpathMonitor_fired == 1 )
410 if( HLTpathMain_fired == 1)
MonitorElement * m_HT_finalSel
MonitorElement * m_MjjWide_num_NOdeta
MonitorElement * m_MjjWide_num_detaL4
MonitorElement * m_MjjWide_den_detaL3
MonitorElement * m_HT_inclusive
MonitorElement * m_MjjWide_den_NOdeta
bool configurationUpdated() const
MonitorElement * m_MjjWide_num_detaL2
unsigned int numwidejets_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
MonitorElement * m_selJets_towersArea
MonitorElement * m_metCaseNoMetClean
edm::EDGetTokenT< reco::CaloJetCollection > jetCollectionTagToken_
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
MonitorElement * m_selJets_phi
MonitorElement * m_MjjWide_num
MonitorElement * m_selJets_pt
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::EDGetTokenT< reco::CaloMETCollection > metCollectionTagToken_
triggerExpression::Evaluator * HLTpathMain_
MonitorElement * m_DphijjWide_finalSel
MonitorElement * m_MjjWide_deta_0p5_1p0
edm::EDGetTokenT< reco::CaloMETCollection > metCleanCollectionTagToken_
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
edm::EDGetTokenT< std::vector< math::PtEtaPhiMLorentzVector > > widejetsCollectionTagToken_
MonitorElement * m_DetajjWide
MonitorElement * m_selJets_hadEnergyFraction
MonitorElement * m_metCases