304 if (!isFound)
return;
308 if (!isFound)
return;
313 if (!isFound)
return;
314 std::vector<Jet>
jets = *cJets;
315 std::sort(jets.begin(), jets.end(),
PtGreater());
319 if (!isFound)
return;
323 if (!isFound)
return;
330 for (
typename std::vector<Jet>::const_iterator jet_i = jets.begin();
331 jet_i != jets.end(); ++jet_i) {
337 sumPt += jet_i->pt();
343 for (
typename std::vector<Met>::const_iterator met_i = mets->begin();
344 met_i != mets->end(); ++met_i) {
359 for (
typename std::vector<Mu>::const_iterator mu_i = muons->begin();
360 mu_i != muons->end(); ++mu_i) {
372 for (
typename std::vector<Mu>::const_iterator mu_j = muons->begin();
373 mu_j != muons->end(); ++mu_j) {
374 if (mu_i >= mu_j)
continue;
377 inv = (mu_i->p4() + mu_j->p4()).M();
378 if (mu_i->charge() * mu_j->charge() > 0) {
382 if (mu_i->charge() * mu_j->charge() < 0) {
389 for (
typename std::vector<Ele>::const_iterator ele_j = elecs->begin();
390 ele_j != elecs->end(); ++ele_j) {
393 inv = (mu_i->p4() + ele_j->p4()).M();
394 dR =
deltaR(*mu_i, *ele_j);
396 if (mu_i->charge() * ele_j->charge() > 0) {
400 if (mu_i->charge() * ele_j->charge() < 0) {
412 for (
typename std::vector<Ele>::const_iterator ele_i = elecs->begin();
413 ele_i != elecs->end(); ++ele_i) {
423 ele_i->dr03HcalTowerSumEt());
426 for (
typename std::vector<Ele>::const_iterator ele_j = elecs->begin();
427 ele_j != elecs->end(); ++ele_j) {
428 if (ele_i >= ele_j)
continue;
432 inv = (ele_i->p4() + ele_j->p4()).M();
433 if (ele_i->charge() * ele_j->charge() > 0) {
437 if (ele_i->charge() * ele_j->charge() < 0) {
450 if (nElectrons >= 1) {
MonitorElement * hRAL_Iso_elecs_
MonitorElement * hRAL_OSemu_monitor_
MonitorElement * hRAL_N_muons_
edm::EDGetTokenT< std::vector< reco::GsfElectron > > electrons_
float sumPt
sum-pt of tracks
bool getByToken(EDGetToken token, Handle< PROD > &result) const
MonitorElement * hRAL_Iso_muons_
virtual bool goodSusyMuon(const Mu *)
MonitorElement * hRAL_SSmumu_monitor_
MonitorElement * hRAL_Sum_pt_jets_
edm::EDGetTokenT< reco::VertexCollection > vertex_
double RAL_jet_sum_pt_cut_
MonitorElement * hRAL_eta_elecs_
edm::EDGetTokenT< std::vector< reco::Muon > > muons_
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_
edm::EDGetTokenT< std::vector< reco::CaloJet > > jets_
edm::EDGetTokenT< std::vector< reco::CaloMET > > met_
MonitorElement * hRAL_mass_OS_emu_
MonitorElement * hRAL_mass_OS_mumu_
MonitorElement * hRAL_OSmumu_monitor_
virtual bool goodSusyElectron(const Ele *)
MonitorElement * hRAL_Muon_monitor_
MonitorElement * hRAL_pt_elecs_
MonitorElement * hRAL_eta_muons_
MonitorElement * hRAL_phi_elecs_
MonitorElement * hRAL_pt_muons_
MonitorElement * hRAL_N_elecs_
double deltaR(double eta1, double eta2, double phi1, double phi2)
MonitorElement * hRAL_Electron_monitor_
MonitorElement * hRAL_Met_
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 * hRAL_dR_emu_