46 return (i.pt() > j.pt());
50 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
144 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
158 muon_eta_cut_ = pset.
getParameter<
double> (
"muon_eta_cut");
159 muon_nHits_cut_ = pset.
getParameter<
double> (
"muon_nHits_cut");
160 muon_nChi2_cut_ = pset.
getParameter<
double> (
"muon_nChi2_cut");
161 muon_d0_cut_ = pset.
getParameter<
double> (
"muon_d0_cut");
163 elec_eta_cut_ = pset.
getParameter<
double> (
"elec_eta_cut");
164 elec_mva_cut_ = pset.
getParameter<
double> (
"elec_mva_cut");
165 elec_d0_cut_ = pset.
getParameter<
double> (
"elec_d0_cut");
167 RAL_muon_pt_cut_ = pset.
getParameter<
double> (
"RAL_muon_pt_cut");
168 RAL_muon_iso_cut_ = pset.
getParameter<
double> (
"RAL_muon_iso_cut");
170 RAL_elec_pt_cut_ = pset.
getParameter<
double> (
"RAL_elec_pt_cut");
171 RAL_elec_iso_cut_ = pset.
getParameter<
double> (
"RAL_elec_iso_cut");
173 RAL_jet_pt_cut_ = pset.
getParameter<
double> (
"RAL_jet_pt_cut");
174 RAL_jet_sum_pt_cut_ = pset.
getParameter<
double> (
"RAL_jet_sum_pt_cut");
175 RAL_jet_eta_cut_ = pset.
getParameter<
double> (
"RAL_jet_eta_cut");
176 RAL_jet_min_emf_cut_ = pset.
getParameter<
double> (
"RAL_jet_min_emf_cut");
177 RAL_jet_max_emf_cut_ = pset.
getParameter<
double> (
"RAL_jet_max_emf_cut");
179 RAL_met_cut_ = pset.
getParameter<
double> (
"RAL_met_cut");
193 hRAL_Sum_pt_jets_ = 0;
196 hRAL_mass_OS_mumu_ = 0;
197 hRAL_mass_OS_ee_ = 0;
198 hRAL_mass_OS_emu_ = 0;
199 hRAL_mass_SS_mumu_ = 0;
200 hRAL_mass_SS_ee_ = 0;
201 hRAL_mass_SS_emu_ = 0;
202 hRAL_Muon_monitor_ = 0;
203 hRAL_Electron_monitor_ = 0;
204 hRAL_OSee_monitor_ = 0;
205 hRAL_OSemu_monitor_ = 0;
206 hRAL_OSmumu_monitor_ = 0;
207 hRAL_SSee_monitor_ = 0;
208 hRAL_SSemu_monitor_ = 0;
209 hRAL_SSmumu_monitor_ = 0;
210 hRAL_TriMuon_monitor_ = 0;
214 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
219 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
224 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
229 hRAL_N_muons_ =
dbe_->
book1D(
"RAL_N_muons",
"RAL_N_muons", 10, 0., 10.);
230 hRAL_pt_muons_ =
dbe_->
book1D(
"RAL_pt_muons",
"RAL_pt_muons", 50, 0., 300.);
231 hRAL_eta_muons_ =
dbe_->
book1D(
"RAL_eta_muons",
"RAL_eta_muons", 50, -2.5, 2.5);
232 hRAL_phi_muons_ =
dbe_->
book1D(
"RAL_phi_muons",
"RAL_phi_muons", 50, -4., 4.);
233 hRAL_Iso_muons_ =
dbe_->
book1D(
"RAL_Iso_muons",
"RAL_Iso_muons", 50, 0., 25.);
235 hRAL_N_elecs_ =
dbe_->
book1D(
"RAL_N_elecs",
"RAL_N_elecs", 10, 0., 10.);
236 hRAL_pt_elecs_ =
dbe_->
book1D(
"RAL_pt_elecs",
"RAL_pt_elecs", 50, 0., 300.);
237 hRAL_eta_elecs_ =
dbe_->
book1D(
"RAL_eta_elecs",
"RAL_eta_elecs", 50, -2.5, 2.5);
238 hRAL_phi_elecs_ =
dbe_->
book1D(
"RAL_phi_elecs",
"RAL_phi_elecs", 50, -4., 4.);
239 hRAL_Iso_elecs_ =
dbe_->
book1D(
"RAL_Iso_elecs",
"RAL_Iso_elecs", 50, 0., 25.);
241 hRAL_Sum_pt_jets_ =
dbe_->
book1D(
"RAL_Sum_pt_jets",
"RAL_Sum_pt_jets", 50, 0., 2000.);
242 hRAL_Met_ =
dbe_->
book1D(
"RAL_Met",
"RAL_Met", 50, 0., 1000.);
244 hRAL_dR_emu_ =
dbe_->
book1D(
"RAL_deltaR_emu",
"RAL_deltaR_emu", 50, 0., 10.);
246 hRAL_mass_OS_mumu_ =
dbe_->
book1D(
"RAL_mass_OS_mumu",
"RAL_mass_OS_mumu", 50, 0., 300.);
247 hRAL_mass_OS_ee_ =
dbe_->
book1D(
"RAL_mass_OS_ee",
"RAL_mass_OS_ee", 50, 0., 300.);
248 hRAL_mass_OS_emu_ =
dbe_->
book1D(
"RAL_mass_OS_emu",
"RAL_mass_OS_emu", 50, 0., 300.);
249 hRAL_mass_SS_mumu_ =
dbe_->
book1D(
"RAL_mass_SS_mumu",
"RAL_mass_SS_mumu", 50, 0., 300.);
250 hRAL_mass_SS_ee_ =
dbe_->
book1D(
"RAL_mass_SS_ee",
"RAL_mass_SS_ee", 50, 0., 300.);
251 hRAL_mass_SS_emu_ =
dbe_->
book1D(
"RAL_mass_SS_emu",
"RAL_mass_SS_emu", 50, 0., 300.);
253 hRAL_Muon_monitor_ =
dbe_->
book2D(
"RAL_Single_Muon_Selection",
"RAL_Single_Muon_Selection", 50, 0., 1000., 50, 0.,
255 hRAL_Electron_monitor_ =
dbe_->
book2D(
"RAL_Single_Electron_Selection",
"RAL_Single_Electron_Selection", 50, 0.,
256 1000., 50, 0., 1000.);
257 hRAL_OSee_monitor_ =
dbe_->
book2D(
"RAL_OS_Electron_Selection",
"RAL_OS_Electron_Selection", 50, 0., 1000., 50, 0.,
259 hRAL_OSemu_monitor_ =
dbe_->
book2D(
"RAL_OS_ElectronMuon_Selection",
"RAL_OS_ElectronMuon_Selection", 50, 0., 1000.,
261 hRAL_OSmumu_monitor_ =
dbe_->
book2D(
"RAL_OS_Muon_Selection",
"RAL_OS_Muon_Selection", 50, 0., 1000., 50, 0., 1000.);
262 hRAL_SSee_monitor_ =
dbe_->
book2D(
"RAL_SS_Electron_Selection",
"RAL_SS_Electron_Selection", 50, 0., 1000., 50, 0.,
264 hRAL_SSemu_monitor_ =
dbe_->
book2D(
"RAL_SS_ElectronMuon_Selection",
"RAL_SS_ElectronMuon_Selection", 50, 0., 1000.,
266 hRAL_SSmumu_monitor_ =
dbe_->
book2D(
"RAL_SS_Muon_Selection",
"RAL_SS_Muon_Selection", 50, 0., 1000., 50, 0., 1000.);
267 hRAL_TriMuon_monitor_ =
dbe_->
book2D(
"RAL_Tri_Muon_Selection",
"RAL_Tri_Muon_Selection", 50, 0., 1000., 50, 0.,
272 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
277 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
281 if (fabs(ele->eta()) > elec_eta_cut_)
285 if (fabs(ele->gsfTrack()->dxy(bs)) > elec_d0_cut_)
290 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
294 if (fabs(mu->eta()) > muon_eta_cut_)
296 if (!mu->isGlobalMuon())
298 if (mu->innerTrack()->numberOfValidHits() < muon_nHits_cut_)
300 if (mu->globalTrack()->normalizedChi2() > muon_nChi2_cut_)
302 if (fabs(mu->innerTrack()->dxy(bs)) > muon_d0_cut_)
307 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
328 std::vector<Jet>
jets = *cJets;
346 for (
typename std::vector<Jet>::const_iterator jet_i = jets.begin(); jet_i != jets.end(); ++jet_i) {
347 if (jet_i->pt() < RAL_jet_pt_cut_)
349 if (fabs(jet_i->eta()) > RAL_jet_eta_cut_)
351 if (fabs(jet_i->eta()) > RAL_jet_eta_cut_)
353 if (jet_i->emEnergyFraction() < RAL_jet_min_emf_cut_)
355 if (jet_i->emEnergyFraction() > RAL_jet_max_emf_cut_)
357 sumPt += jet_i->pt();
360 hRAL_Sum_pt_jets_->Fill(sumPt);
363 for (
typename std::vector<Met>::const_iterator met_i = mets->begin(); met_i != mets->end(); ++met_i) {
368 hRAL_Met_->Fill(MET);
378 for (
typename std::vector<Mu>::const_iterator mu_i = muons->begin(); mu_i != muons->end(); ++mu_i) {
379 if (!(goodSusyMuon(&(*mu_i)) && mu_i->pt() > RAL_muon_pt_cut_))
383 hRAL_pt_muons_->Fill(mu_i->pt());
384 hRAL_eta_muons_->Fill(mu_i->eta());
385 hRAL_phi_muons_->Fill(mu_i->phi());
391 for (
typename std::vector<Mu>::const_iterator mu_j = muons->begin(); mu_j != muons->end(); ++mu_j) {
394 if (!(goodSusyMuon(&(*mu_j)) && mu_j->pt() > RAL_muon_pt_cut_))
397 inv = (mu_i->p4() + mu_j->p4()).M();
398 if (mu_i->charge() * mu_j->charge() > 0) {
400 hRAL_mass_SS_mumu_->Fill(inv);
402 if (mu_i->charge() * mu_j->charge() < 0) {
404 hRAL_mass_OS_mumu_->Fill(inv);
409 for (
typename std::vector<Ele>::const_iterator ele_j = elecs->begin(); ele_j != elecs->end(); ++ele_j) {
410 if (!(goodSusyElectron(&(*ele_j)) && ele_j->pt() > RAL_elec_pt_cut_))
412 inv = (mu_i->p4() + ele_j->p4()).M();
413 dR =
deltaR(*mu_i, *ele_j);
414 hRAL_dR_emu_->Fill(dR);
415 if (mu_i->charge() * ele_j->charge() > 0) {
417 hRAL_mass_SS_emu_->Fill(inv);
419 if (mu_i->charge() * ele_j->charge() < 0) {
421 hRAL_mass_OS_emu_->Fill(inv);
426 hRAL_N_muons_->Fill(nMuons);
431 for (
typename std::vector<Ele>::const_iterator ele_i = elecs->begin(); ele_i != elecs->end(); ++ele_i) {
432 if (!(goodSusyElectron(&(*ele_i)) && ele_i->pt() > RAL_elec_pt_cut_))
436 hRAL_pt_elecs_->Fill(ele_i->pt());
437 hRAL_eta_elecs_->Fill(ele_i->eta());
438 hRAL_phi_elecs_->Fill(ele_i->phi());
440 hRAL_Iso_elecs_->Fill(ele_i->dr03TkSumPt() + ele_i->dr03EcalRecHitSumEt() + ele_i->dr03HcalTowerSumEt());
443 for (
typename std::vector<Ele>::const_iterator ele_j = elecs->begin(); ele_j != elecs->end(); ++ele_j) {
446 if (!(goodSusyElectron(&(*ele_j)) && ele_j->pt() > RAL_elec_pt_cut_))
449 inv = (ele_i->p4() + ele_j->p4()).M();
450 if (ele_i->charge() * ele_j->charge() > 0) {
452 hRAL_mass_SS_ee_->Fill(inv);
454 if (ele_i->charge() * ele_j->charge() < 0) {
456 hRAL_mass_OS_ee_->Fill(inv);
461 hRAL_N_elecs_->Fill(nElectrons);
463 if (MET > RAL_met_cut_ && sumPt > RAL_jet_sum_pt_cut_) {
465 hRAL_Muon_monitor_->Fill(sumPt, MET);
467 if (nElectrons >= 1) {
468 hRAL_Electron_monitor_->Fill(sumPt, MET);
471 hRAL_OSee_monitor_->Fill(sumPt, MET);
474 hRAL_OSemu_monitor_->Fill(sumPt, MET);
477 hRAL_OSmumu_monitor_->Fill(sumPt, MET);
480 hRAL_SSee_monitor_->Fill(sumPt, MET);
483 hRAL_SSemu_monitor_->Fill(sumPt, MET);
486 hRAL_SSmumu_monitor_->Fill(sumPt, MET);
490 hRAL_TriMuon_monitor_->Fill(sumPt, MET);
495 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
500 template<
typename Mu,
typename Ele,
typename Jet,
typename Met>
MonitorElement * hRAL_Iso_elecs_
MonitorElement * hRAL_OSemu_monitor_
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
MonitorElement * hRAL_N_muons_
float sumPt
sum-pt of tracks
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
SusyDQM(const edm::ParameterSet &)
MonitorElement * hRAL_Iso_muons_
virtual bool goodSusyMuon(const Mu *)
MonitorElement * hRAL_SSmumu_monitor_
MonitorElement * hRAL_Sum_pt_jets_
double RAL_jet_sum_pt_cut_
MonitorElement * hRAL_eta_elecs_
double RAL_jet_max_emf_cut_
MonitorElement * hRAL_mass_SS_mumu_
double RAL_jet_min_emf_cut_
MonitorElement * hRAL_mass_SS_emu_
MonitorElement * hRAL_OSee_monitor_
MonitorElement * hRAL_mass_OS_emu_
MonitorElement * hRAL_mass_OS_mumu_
MonitorElement * hRAL_OSmumu_monitor_
virtual bool goodSusyElectron(const Ele *)
void endRun(const edm::Run &)
MonitorElement * hRAL_Muon_monitor_
MonitorElement * hRAL_pt_elecs_
MonitorElement * hRAL_eta_muons_
MonitorElement * hRAL_phi_elecs_
edm::ParameterSet parameters_
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
virtual void analyze(const edm::Event &, const edm::EventSetup &)
MonitorElement * hRAL_pt_muons_
MonitorElement * hRAL_N_elecs_
double deltaR(double eta1, double eta2, double phi1, double phi2)
MonitorElement * hRAL_Electron_monitor_
void beginRun(const edm::Run &)
MonitorElement * hRAL_Met_
XYZPointD XYZPoint
point in space with cartesian internal representation
MonitorElement * hRAL_SSee_monitor_
MonitorElement * hRAL_SSemu_monitor_
MonitorElement * hRAL_mass_OS_ee_
MonitorElement * hRAL_TriMuon_monitor_
MonitorElement * hRAL_phi_muons_
MonitorElement * hRAL_mass_SS_ee_
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
bool operator()(const T &i, const T &j)
void setCurrentFolder(const std::string &fullpath)
MonitorElement * hRAL_dR_emu_
SusyDQM< reco::Muon, reco::GsfElectron, reco::CaloJet, reco::CaloMET > RecoSusyDQM