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)
MonitorElement * m_HT_finalSel
MonitorElement * m_MjjWide_num_NOdeta
MonitorElement * m_MjjWide_num_detaL4
MonitorElement * m_MjjWide_den_detaL3
edm::InputTag metCleanCollectionTag_
MonitorElement * m_HT_inclusive
MonitorElement * m_MjjWide_den_NOdeta
bool configurationUpdated() const
MonitorElement * m_MjjWide_num_detaL2
unsigned int numwidejets_
MonitorElement * m_selJets_towersArea
edm::InputTag metCollectionTag_
MonitorElement * m_metCaseNoMetClean
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
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::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
MonitorElement * m_selJets_hadEnergyFraction
MonitorElement * m_metCases