483 std::vector<Jet>
jets = *cJets;
503 float RA1_jet1_pt = 0;
504 float RA1_jet1_eta = 0;
505 float RA1_jet1_emf = 0;
506 float RA1_jet2_pt = 0;
507 float RA1_jet2_eta = 0;
508 float RA1_jet2_emf = 0;
509 float RA1_jet3_pt = 0;
513 for (
typename std::vector<Jet>::const_iterator jet_i = jets.begin(); jet_i != jets.end(); ++jet_i) {
515 RA1_leading = jet_i->p4();
516 RA1_jet1_pt = jet_i->pt();
517 RA1_jet1_eta = jet_i->eta();
518 RA1_jet1_emf = jet_i->emEnergyFraction();
521 RA1_second = jet_i->p4();
522 RA1_jet2_pt = jet_i->pt();
523 RA1_jet2_eta = jet_i->eta();
524 RA1_jet2_emf = jet_i->emEnergyFraction();
527 RA1_jet3_pt = jet_i->pt();
530 RA1_HT += jet_i->pt();
531 RA1_vMHT -= jet_i->p4();
535 float RA1_MHT = RA1_vMHT.pt();
538 float RA1_minDeltaPhi = 9999.;
539 for (
typename std::vector<Jet>::const_iterator jet_i = jets.begin(); jet_i != jets.end(); ++jet_i) {
541 double deltaPhi_tmp = fabs(
deltaPhi(jet_i->phi(), RA1_vMHT.phi()));
542 if (deltaPhi_tmp < RA1_minDeltaPhi)
543 RA1_minDeltaPhi = deltaPhi_tmp;
548 float RA1_alphat = 0;
549 float RA1_DeltaPhiJets = 9999.;
550 if (RA1_nJets >= 2) {
551 RA1_DeltaPhiJets = fabs(
deltaPhi(RA1_leading.phi(), RA1_second.phi()));
555 RA1_alphat = RA1_second.Et() /
sqrt(2* RA1_leading .Et() * RA1_second.Et() * (1 -
cos(RA1_DeltaPhiJets)));
561 float RA2_jet1_pt = 0;
562 float RA2_jet1_eta = 0;
563 float RA2_jet1_emf = 0;
564 float RA2_jet2_pt = 0;
565 float RA2_jet2_eta = 0;
566 float RA2_jet2_emf = 0;
567 float RA2_jet3_pt = 0;
568 float RA2_jet3_eta = 0;
569 float RA2_jet3_emf = 0;
571 for (
typename std::vector<Jet>::const_iterator jet_i = jets.begin(); jet_i != jets.end(); ++jet_i) {
573 RA2_jet1_pt = jet_i->pt();
574 RA2_jet1_eta = jet_i->eta();
575 RA2_jet1_emf = jet_i->emEnergyFraction();
578 RA2_jet2_pt = jet_i->pt();
579 RA2_jet2_eta = jet_i->eta();
580 RA2_jet2_emf = jet_i->emEnergyFraction();
583 RA2_jet3_pt = jet_i->pt();
584 RA2_jet3_eta = jet_i->eta();
585 RA2_jet3_emf = jet_i->emEnergyFraction();
589 RA2_HT += jet_i->pt();
590 RA2_vMHT -= jet_i->p4();
594 float RA2_MHT = RA2_vMHT.pt();
597 float RA2_minDeltaPhi = 9999.;
598 for (
typename std::vector<Jet>::const_iterator jet_i = jets.begin(); jet_i != jets.end(); ++jet_i) {
602 double deltaPhi_tmp = fabs(
deltaPhi(jet_i->phi(), RA1_vMHT.phi()));
603 if (deltaPhi_tmp < RA2_minDeltaPhi)
604 RA2_minDeltaPhi = deltaPhi_tmp;
609 for (reco::VertexCollection::const_iterator vertex = vertices->begin(); vertex != vertices->end(); ++vertex) {
610 bs = vertex->position();
614 float leadingMuPt = 0;
615 for (
typename std::vector<Mu>::const_iterator mu_i = muons->begin(); mu_i != muons->end(); ++mu_i) {
620 float muIso = (Iso_muon.
emEt + Iso_muon.
hadEt + Iso_muon.
sumPt) / mu_i->pt();
623 if (mu_i->pt() > leadingMuPt)
624 leadingMuPt = mu_i->pt();
628 float leadingElecPt = 0;
629 for (
typename std::vector<Ele>::const_iterator ele_i = elecs->begin(); ele_i != elecs->end(); ++ele_i) {
633 float elecIso = (ele_i->dr03TkSumPt() + ele_i->dr03EcalRecHitSumEt() + ele_i->dr03HcalTowerSumEt()) / ele_i->pt();
636 if (ele_i->pt() > leadingElecPt)
637 leadingElecPt = ele_i->pt();
656 for (
int i = 0;
i < 12; ++
i) {
724 for (
int i = 0;
i < 12; ++
i) {
781 for (
typename std::vector<Jet>::const_iterator jet_i = jets.begin(); jet_i != jets.end(); ++jet_i) {
792 sumPt += jet_i->pt();
798 for (
typename std::vector<Met>::const_iterator met_i = mets->begin(); met_i != mets->end(); ++met_i) {
813 for (
typename std::vector<Mu>::const_iterator mu_i = muons->begin(); mu_i != muons->end(); ++mu_i) {
826 for (
typename std::vector<Mu>::const_iterator mu_j = muons->begin(); mu_j != muons->end(); ++mu_j) {
832 inv = (mu_i->p4() + mu_j->p4()).M();
833 if (mu_i->charge() * mu_j->charge() > 0) {
837 if (mu_i->charge() * mu_j->charge() < 0) {
844 for (
typename std::vector<Ele>::const_iterator ele_j = elecs->begin(); ele_j != elecs->end(); ++ele_j) {
847 inv = (mu_i->p4() + ele_j->p4()).M();
848 dR =
deltaR(*mu_i, *ele_j);
850 if (mu_i->charge() * ele_j->charge() > 0) {
854 if (mu_i->charge() * ele_j->charge() < 0) {
866 for (
typename std::vector<Ele>::const_iterator ele_i = elecs->begin(); ele_i != elecs->end(); ++ele_i) {
875 hRAL_Iso_elecs_->
Fill(ele_i->dr03TkSumPt() + ele_i->dr03EcalRecHitSumEt() + ele_i->dr03HcalTowerSumEt());
878 for (
typename std::vector<Ele>::const_iterator ele_j = elecs->begin(); ele_j != elecs->end(); ++ele_j) {
884 inv = (ele_i->p4() + ele_j->p4()).M();
885 if (ele_i->charge() * ele_j->charge() > 0) {
889 if (ele_i->charge() * ele_j->charge() < 0) {
902 if (nElectrons >= 1) {
MonitorElement * hRA1_pt_jet1_nm1_
MonitorElement * hRAL_Iso_elecs_
MonitorElement * hRAL_OSemu_monitor_
MonitorElement * hRA2_pt_jet3_nm1_
MonitorElement * hRA1_eta_jet2_nm1_
MonitorElement * hRA1_eta_jet2_
MonitorElement * hRA2_pt_elecs_
MonitorElement * hRAL_N_muons_
MonitorElement * hRA2_deltaPhi_mht_jets_nm1_
float sumPt
sum-pt of tracks
MonitorElement * hRA2_N_jets_nm1_
MonitorElement * hRA2_pt_jet1_nm1_
MonitorElement * hRA1_pt_elecs_
double deltaPhi(float phi1, float phi2)
MonitorElement * hRAL_Iso_muons_
MonitorElement * hRA2_eta_jet2_
virtual bool goodSusyMuon(const Mu *)
MonitorElement * hRAL_SSmumu_monitor_
MonitorElement * hRAL_Sum_pt_jets_
MonitorElement * hRA1_pt_jet3_nm1_
MonitorElement * hRA1_deltaPhi_mht_jets_
MonitorElement * hRA2_ht_
MonitorElement * hRA2_pt_jet3_
double RAL_jet_sum_pt_cut_
MonitorElement * hRA1_pt_muons_
MonitorElement * hRAL_eta_elecs_
MonitorElement * hRA1_mht_
MonitorElement * hRA1_pt_muons_nm1_
double RAL_jet_max_emf_cut_
MonitorElement * hRA2_N_jets_
MonitorElement * hRA2_eta_jet3_nm1_
MonitorElement * hRAL_mass_SS_mumu_
MonitorElement * hRA2_pt_muons_
MonitorElement * hRA1_pt_elecs_nm1_
MonitorElement * hRA2_pt_jet2_
double RAL_jet_min_emf_cut_
double RA12_elec_iso_cut_
MonitorElement * hRA1_alphat_
MonitorElement * hRAL_mass_SS_emu_
MonitorElement * hRA2_pt_jet2_nm1_
MonitorElement * hRA2_pt_muons_nm1_
MonitorElement * hRA1_pt_jet1_
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
MonitorElement * hRA1_ht_
MonitorElement * hRAL_OSee_monitor_
double RA2_jet_max_emf_cut_
MonitorElement * hRA1_mht_nm1_
double RA2_jet_min_emf_cut_
double RA12_muon_iso_cut_
MonitorElement * hRA1_deltaPhi_jets_nm1_
MonitorElement * hRAL_mass_OS_emu_
MonitorElement * hRAL_mass_OS_mumu_
Cos< T >::type cos(const T &t)
MonitorElement * hRAL_OSmumu_monitor_
virtual bool goodSusyElectron(const Ele *)
MonitorElement * hRA1_alphat_nm1_
MonitorElement * hRAL_Muon_monitor_
MonitorElement * hRAL_pt_elecs_
MonitorElement * hRA2_mht_nm1_
MonitorElement * hRA2_ht_nm1_
MonitorElement * hRAL_eta_muons_
MonitorElement * hRAL_phi_elecs_
MonitorElement * hRA1_pt_jet2_
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
double RA1_deltaPhiJets_cut_
MonitorElement * hRA1_deltaPhi_mht_jets_nm1_
MonitorElement * hRAL_pt_muons_
MonitorElement * hRAL_N_elecs_
double deltaR(double eta1, double eta2, double phi1, double phi2)
MonitorElement * hRA2_eta_jet2_nm1_
MonitorElement * hRAL_Electron_monitor_
MonitorElement * hRA2_deltaPhi_mht_jets_
double RA1_jet_min_emf_cut_
MonitorElement * hRA1_eta_jet1_nm1_
MonitorElement * hRA1_ht_nm1_
MonitorElement * hRAL_Met_
MonitorElement * hRA1_deltaPhi_jets_
MonitorElement * hRA1_pt_jet2_nm1_
MonitorElement * hRA2_eta_jet1_
double RA1_jet_max_emf_cut_
MonitorElement * hRAL_SSee_monitor_
MonitorElement * hRAL_SSemu_monitor_
MonitorElement * hRA2_eta_jet1_nm1_
MonitorElement * hRAL_mass_OS_ee_
MonitorElement * hRA2_pt_jet1_
MonitorElement * hRAL_TriMuon_monitor_
MonitorElement * hRA1_eta_jet1_
MonitorElement * hRA2_pt_elecs_nm1_
MonitorElement * hRAL_phi_muons_
MonitorElement * hRA1_pt_jet3_
MonitorElement * hRAL_mass_SS_ee_
MonitorElement * hRA2_eta_jet3_
MonitorElement * hRAL_dR_emu_
MonitorElement * hRA2_mht_