320 if (!isFound)
return;
324 if (!isFound)
return;
329 if (!isFound)
return;
330 std::vector<Jet>
jets = *cJets;
335 if (!isFound)
return;
339 if (!isFound)
return;
346 for (
typename std::vector<Jet>::const_iterator jet_i = jets.begin();
347 jet_i != jets.end(); ++jet_i) {
353 sumPt += jet_i->pt();
359 for (
typename std::vector<Met>::const_iterator met_i = mets->begin();
360 met_i != mets->end(); ++met_i) {
375 for (
typename std::vector<Mu>::const_iterator mu_i = muons->begin();
376 mu_i != muons->end(); ++mu_i) {
388 for (
typename std::vector<Mu>::const_iterator mu_j = muons->begin();
389 mu_j != muons->end(); ++mu_j) {
390 if (mu_i >= mu_j)
continue;
393 inv = (mu_i->p4() + mu_j->p4()).M();
394 if (mu_i->charge() * mu_j->charge() > 0) {
398 if (mu_i->charge() * mu_j->charge() < 0) {
405 for (
typename std::vector<Ele>::const_iterator ele_j = elecs->begin();
406 ele_j != elecs->end(); ++ele_j) {
409 inv = (mu_i->p4() + ele_j->p4()).M();
410 dR =
deltaR(*mu_i, *ele_j);
412 if (mu_i->charge() * ele_j->charge() > 0) {
416 if (mu_i->charge() * ele_j->charge() < 0) {
428 for (
typename std::vector<Ele>::const_iterator ele_i = elecs->begin();
429 ele_i != elecs->end(); ++ele_i) {
439 ele_i->dr03HcalTowerSumEt());
442 for (
typename std::vector<Ele>::const_iterator ele_j = elecs->begin();
443 ele_j != elecs->end(); ++ele_j) {
444 if (ele_i >= ele_j)
continue;
448 inv = (ele_i->p4() + ele_j->p4()).M();
449 if (ele_i->charge() * ele_j->charge() > 0) {
453 if (ele_i->charge() * ele_j->charge() < 0) {
466 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_