file analyses/Analysis_CMS_13TeV_2SSLEP_Stop_137invfb.cpp
[No description available]
Namespaces
| Name |
|---|
| Gambit TODO: see if we can use this one: |
| Gambit::ColliderBit |
Classes
| Name | |
|---|---|
| class | Gambit::ColliderBit::Analysis_CMS_13TeV_2SSLEP_Stop_137invfb |
| struct | Gambit::ColliderBit::Analysis_CMS_13TeV_2SSLEP_Stop_137invfb::ptComparison |
Source code
///
/// \author Yang Zhang
/// \date 2019 June
/// *********************************************
// Based on:
// - http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-19-008/index.html
// - http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-19-008/index.html
// Search for physics beyond the standard model in events with two same-sign leptons or at least three leptons and jets in proton-proton collisions at 13 TeV
// Note:
// 1. Not fully validated.
#include <vector>
#include <cmath>
#include <memory>
#include <iomanip>
#include "gambit/ColliderBit/analyses/Analysis.hpp"
#include "gambit/ColliderBit/mt2_bisect.h"
#include "gambit/ColliderBit/CMSEfficiencies.hpp"
#include "gambit/ColliderBit/analyses/Cutflow.hpp"
// #define CHECK_CUTFLOW
using namespace std;
namespace Gambit {
namespace ColliderBit {
class Analysis_CMS_13TeV_2SSLEP_Stop_137invfb : public Analysis {
protected:
// Counters for the number of accepted events for each signal region
std::map<string, EventCounter> _counters = {
// HH
{"ISR1", EventCounter("ISR1")},
{"ISR2", EventCounter("ISR2")},
{"ISR3", EventCounter("ISR3")},
{"ISR4", EventCounter("ISR4")},
{"ISR5", EventCounter("ISR5")},
{"ISR6", EventCounter("ISR6")},
{"ISR7", EventCounter("ISR7")},
{"ISR8", EventCounter("ISR8")},
{"ISR9", EventCounter("ISR9")},
{"ISR10", EventCounter("ISR10")},
{"ISR11", EventCounter("ISR11")},
{"ISR12", EventCounter("ISR12")},
{"ISR13", EventCounter("ISR13")},
{"ISR14", EventCounter("ISR14")},
{"ISR15", EventCounter("ISR15")},
{"ISR16", EventCounter("ISR16")},
{"ISR17", EventCounter("ISR17")},
// {"SRHH-0", EventCounter("SRHH-0")},
// {"SRHH-1", EventCounter("SRHH-1")},
// {"SRHH-2", EventCounter("SRHH-2")},
// {"SRHH-3", EventCounter("SRHH-3")},
// {"SRHH-4", EventCounter("SRHH-4")},
// {"SRHH-5", EventCounter("SRHH-5")},
// {"SRHH-6", EventCounter("SRHH-6")},
// {"SRHH-7", EventCounter("SRHH-7")},
// {"SRHH-8", EventCounter("SRHH-8")},
// {"SRHH-9", EventCounter("SRHH-9")},
// {"SRHH-10", EventCounter("SRHH-10")},
// {"SRHH-11", EventCounter("SRHH-11")},
// {"SRHH-12", EventCounter("SRHH-12")},
// {"SRHH-13", EventCounter("SRHH-13")},
// {"SRHH-14", EventCounter("SRHH-14")},
// {"SRHH-15", EventCounter("SRHH-15")},
// {"SRHH-16", EventCounter("SRHH-16")},
// {"SRHH-17", EventCounter("SRHH-17")},
// {"SRHH-18", EventCounter("SRHH-18")},
// {"SRHH-19", EventCounter("SRHH-19")},
// {"SRHH-20", EventCounter("SRHH-20")},
// {"SRHH-21", EventCounter("SRHH-21")},
// {"SRHH-22", EventCounter("SRHH-22")},
// {"SRHH-23", EventCounter("SRHH-23")},
// {"SRHH-24", EventCounter("SRHH-24")},
// {"SRHH-25", EventCounter("SRHH-25")},
// {"SRHH-26", EventCounter("SRHH-26")},
// {"SRHH-27", EventCounter("SRHH-27")},
// {"SRHH-28", EventCounter("SRHH-28")},
// {"SRHH-29", EventCounter("SRHH-29")},
// {"SRHH-30", EventCounter("SRHH-30")},
// {"SRHH-31", EventCounter("SRHH-31")},
// {"SRHH-32", EventCounter("SRHH-32")},
// {"SRHH-33", EventCounter("SRHH-33")},
// {"SRHH-34", EventCounter("SRHH-34")},
// {"SRHH-35", EventCounter("SRHH-35")},
// {"SRHH-36", EventCounter("SRHH-36")},
// {"SRHH-37", EventCounter("SRHH-37")},
// {"SRHH-38", EventCounter("SRHH-38")},
// {"SRHH-39", EventCounter("SRHH-39")},
// {"SRHH-40", EventCounter("SRHH-40")},
// {"SRHH-41", EventCounter("SRHH-41")},
// {"SRHH-42", EventCounter("SRHH-42")},
// {"SRHH-43", EventCounter("SRHH-43")},
// {"SRHH-44", EventCounter("SRHH-44")},
// {"SRHH-45", EventCounter("SRHH-45")},
// {"SRHH-46", EventCounter("SRHH-46")},
// {"SRHH-47", EventCounter("SRHH-47")},
// {"SRHH-48", EventCounter("SRHH-48")},
// {"SRHH-49", EventCounter("SRHH-49")},
// {"SRHH-50", EventCounter("SRHH-50")},
// {"SRHH-51", EventCounter("SRHH-51")},
// {"SRHH-52", EventCounter("SRHH-52")},
// {"SRHH-53", EventCounter("SRHH-53")},
// {"SRHH-54", EventCounter("SRHH-54")},
// {"SRHH-55", EventCounter("SRHH-55")},
// {"SRHH-56", EventCounter("SRHH-56")},
// {"SRHH-57", EventCounter("SRHH-57")},
// {"SRHH-58", EventCounter("SRHH-58")},
// {"SRHH-59", EventCounter("SRHH-59")},
// {"SRHH-60", EventCounter("SRHH-60")},
// {"SRHH-61", EventCounter("SRHH-61")},
// // HL
// {"SRHL-0", EventCounter("SRHL-0")},
// {"SRHL-1", EventCounter("SRHL-1")},
// {"SRHL-2", EventCounter("SRHL-2")},
// {"SRHL-3", EventCounter("SRHL-3")},
// {"SRHL-4", EventCounter("SRHL-4")},
// {"SRHL-5", EventCounter("SRHL-5")},
// {"SRHL-6", EventCounter("SRHL-6")},
// {"SRHL-7", EventCounter("SRHL-7")},
// {"SRHL-8", EventCounter("SRHL-8")},
// {"SRHL-9", EventCounter("SRHL-9")},
// {"SRHL-10", EventCounter("SRHL-10")},
// {"SRHL-11", EventCounter("SRHL-11")},
// {"SRHL-12", EventCounter("SRHL-12")},
// {"SRHL-13", EventCounter("SRHL-13")},
// {"SRHL-14", EventCounter("SRHL-14")},
// {"SRHL-15", EventCounter("SRHL-15")},
// {"SRHL-16", EventCounter("SRHL-16")},
// {"SRHL-17", EventCounter("SRHL-17")},
// {"SRHL-18", EventCounter("SRHL-18")},
// {"SRHL-19", EventCounter("SRHL-19")},
// {"SRHL-20", EventCounter("SRHL-20")},
// {"SRHL-21", EventCounter("SRHL-21")},
// {"SRHL-22", EventCounter("SRHL-22")},
// {"SRHL-23", EventCounter("SRHL-23")},
// {"SRHL-24", EventCounter("SRHL-24")},
// {"SRHL-25", EventCounter("SRHL-25")},
// {"SRHL-26", EventCounter("SRHL-26")},
// {"SRHL-27", EventCounter("SRHL-27")},
// {"SRHL-28", EventCounter("SRHL-28")},
// {"SRHL-29", EventCounter("SRHL-29")},
// {"SRHL-30", EventCounter("SRHL-30")},
// {"SRHL-31", EventCounter("SRHL-31")},
// {"SRHL-32", EventCounter("SRHL-32")},
// {"SRHL-33", EventCounter("SRHL-33")},
// {"SRHL-34", EventCounter("SRHL-34")},
// {"SRHL-35", EventCounter("SRHL-35")},
// {"SRHL-36", EventCounter("SRHL-36")},
// {"SRHL-37", EventCounter("SRHL-37")},
// {"SRHL-38", EventCounter("SRHL-38")},
// {"SRHL-39", EventCounter("SRHL-39")},
// {"SRHL-40", EventCounter("SRHL-40")},
// {"SRHL-41", EventCounter("SRHL-41")},
// {"SRHL-42", EventCounter("SRHL-42")},
// // LL
// {"SRLL-0", EventCounter("SRLL-0")},
// {"SRLL-1", EventCounter("SRLL-1")},
// {"SRLL-2", EventCounter("SRLL-2")},
// {"SRLL-3", EventCounter("SRLL-3")},
// {"SRLL-4", EventCounter("SRLL-4")},
// {"SRLL-5", EventCounter("SRLL-5")},
// {"SRLL-6", EventCounter("SRLL-6")},
// {"SRLL-7", EventCounter("SRLL-7")},
// // LM
// {"SRLM-0", EventCounter("SRLM-0")},
// {"SRLM-1", EventCounter("SRLM-1")},
// {"SRLM-2", EventCounter("SRLM-2")},
// {"SRLM-3", EventCounter("SRLM-3")},
// {"SRLM-4", EventCounter("SRLM-4")},
// {"SRLM-5", EventCounter("SRLM-5")},
// {"SRLM-6", EventCounter("SRLM-6")},
// {"SRLM-7", EventCounter("SRLM-7")},
// {"SRLM-8", EventCounter("SRLM-8")},
// {"SRLM-9", EventCounter("SRLM-9")},
// {"SRLM-10", EventCounter("SRLM-10")},
// // ML
// {"SRML-0", EventCounter("SRML-0")},
// {"SRML-1", EventCounter("SRML-1")},
// {"SRML-2", EventCounter("SRML-2")},
// {"SRML-3", EventCounter("SRML-3")},
// {"SRML-4", EventCounter("SRML-4")},
// {"SRML-5", EventCounter("SRML-5")},
// {"SRML-6", EventCounter("SRML-6")},
// {"SRML-7", EventCounter("SRML-7")},
// {"SRML-8", EventCounter("SRML-8")},
// {"SRML-9", EventCounter("SRML-9")},
// {"SRML-10", EventCounter("SRML-10")},
// {"SRML-11", EventCounter("SRML-11")},
// {"SRML-12", EventCounter("SRML-12")},
// {"SRML-13", EventCounter("SRML-13")},
// {"SRML-14", EventCounter("SRML-14")},
// {"SRML-15", EventCounter("SRML-15")},
// {"SRML-16", EventCounter("SRML-16")},
// {"SRML-17", EventCounter("SRML-17")},
// {"SRML-18", EventCounter("SRML-18")},
// {"SRML-19", EventCounter("SRML-19")},
// {"SRML-20", EventCounter("SRML-20")},
// {"SRML-21", EventCounter("SRML-21")},
// {"SRML-22", EventCounter("SRML-22")},
// {"SRML-23", EventCounter("SRML-23")},
// {"SRML-24", EventCounter("SRML-24")},
// {"SRML-25", EventCounter("SRML-25")},
// {"SRML-26", EventCounter("SRML-26")},
// {"SRML-27", EventCounter("SRML-27")},
// {"SRML-28", EventCounter("SRML-28")},
// {"SRML-29", EventCounter("SRML-29")},
// {"SRML-30", EventCounter("SRML-30")},
// {"SRML-31", EventCounter("SRML-31")},
// {"SRML-32", EventCounter("SRML-32")},
// {"SRML-33", EventCounter("SRML-33")},
// {"SRML-34", EventCounter("SRML-34")},
// {"SRML-35", EventCounter("SRML-35")},
// {"SRML-36", EventCounter("SRML-36")},
// {"SRML-37", EventCounter("SRML-37")},
// {"SRML-38", EventCounter("SRML-38")},
// {"SRML-39", EventCounter("SRML-39")},
// {"SRML-40", EventCounter("SRML-40")},
// {"SRML-41", EventCounter("SRML-41")},
// {"SRML-42", EventCounter("SRML-42")},
// {"SRML-43", EventCounter("SRML-43")},
};
Cutflow _cutflow;
public:
// Required detector sim
static constexpr const char* detector = "CMS";
Analysis_CMS_13TeV_2SSLEP_Stop_137invfb():
_cutflow("CMS_13TeV_2SSLEP_Stop_137invfb", {"Trigger_and_2leptons", "At_least_one_SS_lepton_pair", "Baseline"})
{
set_analysis_name("CMS_13TeV_2SSLEP_Stop_137invfb");
set_luminosity(137);
}
struct ptComparison {
bool operator() (const HEPUtils::Particle* i,const HEPUtils::Particle* j) {return (i->pT()>j->pT());}
} comparePt;
void run(const HEPUtils::Event* event) {
_cutflow.fillinit();
// Missing energy
double met = event->met();
HEPUtils::P4 ptot = event->missingmom();
// Electrons
//@note Numbers digitized from https://twiki.cern.ch/twiki/pub/CMSPublic/SUSMoriond2017ObjectsEfficiency/2d_full_pteta_el_035_ttbar.pdf
const vector<double> aEl={0., 0.8, 1.442, 1.556, 2., 2.5, DBL_MAX}; // Bin edges in eta
const vector<double> bEl={0., 15., 20., 25., 30., 40., 50, DBL_MAX}; // Bin edges in pT. Assume flat efficiency above 200, where the CMS map stops.
const vector<double> cEl={
// pT: (0,15), (15,20), (20,25), (25,30), (30,40), (40,50), (50,inf)
0.0, 0.398, 0.501, 0.556, 0.619, 0.669, 0.720,// eta: (0, 0.8)
0.0, 0.344, 0.433, 0.498, 0.579, 0.600, 0.671,// eta: (0.8, 1.4429)
0.0, 0.201, 0.156, 0.206, 0.222, 0.255, 0.307,// eta: (1.442, 1.556)
0.0, 0.210, 0.302, 0.338, 0.428, 0.484, 0.561,// eta: (1.556, 2)
0.0, 0.162, 0.172, 0.250, 0.339, 0.396, 0.444,// eta: (2, 2.5)
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0// eta > 2.5
};
HEPUtils::BinnedFn2D<double> _eff2dEl(aEl,bEl,cEl);
vector<const HEPUtils::Particle*> electrons;
for (const HEPUtils::Particle* electron : event->electrons()) {
bool isEl=has_tag(_eff2dEl, fabs(electron->eta()), electron->pT());
if (electron->pT() > 15. && fabs(electron->eta()) < 2.5 && isEl)
electrons.push_back(electron);
}
// Muons
//@note Numbers digitized from https://twiki.cern.ch/twiki/pub/CMSPublic/SUSMoriond2017ObjectsEfficiency/2d_full_pteta_mu_035_ttbar.pdf
const vector<double> aMu={0., 0.9, 1.2, 2.1, 2.4, DBL_MAX}; // Bin edges in eta
const vector<double> bMu={0., 10, 15., 20., 25., 30, 40, 50, DBL_MAX}; // Bin edges in pT. Assume flat efficiency above 200, where the CMS map stops.
const vector<double> cMu={
// pT: (0,10), (10,15), (15,20), (20,25), (25,30), (30,40), (40,50), (50,inf)
0.0, 0.564, 0.645, 0.739, 0.803, 0.860, 0.894, 0.907, // eta: (0, 0.9)
0.0, 0.525, 0.616, 0.700, 0.773, 0.825, 0.891, 0.898, // eta: (0.9, 1.2)
0.0, 0.514, 0.572, 0.697, 0.748, 0.789, 0.837, 0.870, // eta: (1.2, 2.1)
0.0, 0.440, 0.575, 0.604, 0.663, 0.696, 0.784, 0.794,// eta: (2.1, 2.4)
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0// eta > 2.4
};
HEPUtils::BinnedFn2D<double> _eff2dMu(aMu,bMu,cMu);
vector<const HEPUtils::Particle*> muons;
for (const HEPUtils::Particle* muon : event->muons()) {
bool isMu=has_tag(_eff2dMu, fabs(muon->eta()), muon->pT());
if (muon->pT() > 10.&& fabs(muon->eta()) < 2.4 && isMu)
muons.push_back(muon);
}
// Jets
double HT = 0.;
vector<const HEPUtils::Jet*> candJets;
for (const HEPUtils::Jet* jet : event->jets("antikt_R04")) {
if (jet->pT() > 25. && fabs(jet->eta()) < 2.4){
HT += jet->pT();
candJets.push_back(jet);
}
}
// Jets
vector<const HEPUtils::Jet*> bJets;
vector<const HEPUtils::Jet*> nonbJets;
// Find b-jets
// Copied from ATLAS_13TeV_3b_24invfb
double btag = 0.85; double cmisstag = 1/12.; double misstag = 1./381.;
for (const HEPUtils::Jet* jet : candJets) {
// Tag
if( jet->btag() && random_bool(btag) ) bJets.push_back(jet);
// Misstag c-jet
else if( jet->ctag() && random_bool(cmisstag) ) bJets.push_back(jet);
// Misstag light jet
else if( random_bool(misstag) ) bJets.push_back(jet);
// Non b-jet
else if( jet->pT() > 40. ) nonbJets.push_back(jet);
}
// // Overlap removal
// JetLeptonOverlapRemoval(candJets,electrons,0.2);
// LeptonJetOverlapRemoval(electrons,candJets);
// JetLeptonOverlapRemoval(candJets,muons,0.4);
// LeptonJetOverlapRemoval(muons,candJets);
size_t Nb=bJets.size();
size_t Nj=nonbJets.size();
// Leptons = electrons + muons
vector<const HEPUtils::Particle*> leptons;
leptons=electrons;
leptons.insert(leptons.end(),muons.begin(),muons.end());
sort(leptons.begin(),leptons.end(),comparePt);
// At least two light leptons
if (leptons.size()<2) return;
// Find pair same sign (SS) leptons
vector<size_t> SS_1,SS_2;
for (size_t i=0; i<leptons.size(); ++i) {
for (size_t j=i+1; j<leptons.size(); ++j) {
if (leptons[i]->pid()*leptons[j]->pid()>0){
SS_1.push_back(i);
SS_2.push_back(j);
}
// mll>12 for an opposite-sign same flavor pair lepton
if (leptons[i]->pid()+leptons[j]->pid()==0 and (leptons[i]->mom()+leptons[j]->mom()).m()<12) return;
// mll>8 GeV for any pair of leptons
if ((leptons[i]->mom()+leptons[j]->mom()).m()<8) return;
}
}
_cutflow.fill(1);
// One SS lepton pair
if (SS_1.size()==0) return;
_cutflow.fill(2);
// At least two jets and MET>50
if (Nj<2 or met<50) return;
_cutflow.fill(3);
// Find the only SS lepton pair
size_t SS1 = SS_1[0];
size_t SS2 = SS_2[0];
bool find_one_muon = false;
for (size_t i=1; i<SS_1.size(); ++i) {
// SS_1 and SS_2 are already order by lepton PT sum
if (fabs(leptons[SS_1[i]]->pid())==13 and fabs(leptons[SS_1[i]]->pid())==13) {
// both of the leptons are muon
SS1 = SS_1[i];
SS2 = SS_2[i];
break;
}
if ( (not find_one_muon) and (fabs(leptons[SS_1[i]]->pid())==13 or fabs(leptons[SS_1[i]]->pid())==13)){
// one of the leptons is muon
SS1 = SS_1[i];
SS2 = SS_2[i];
find_one_muon = true;
}
}
// M_T^{miss}
double MTmin = sqrt(2.*leptons[SS1]->pT()*met*(1. - cos(leptons[SS1]->mom().deltaPhi(ptot))));
if (MTmin>sqrt(2.*leptons[SS2]->pT()*met*(1. - cos(leptons[SS2]->mom().deltaPhi(ptot))))) {
MTmin = sqrt(2.*leptons[SS2]->pT()*met*(1. - cos(leptons[SS2]->mom().deltaPhi(ptot))));
}
// HH: exactly 2 leptons, both with PT>25 GeV, and MET>50 GeV
if (leptons.size()==2 and leptons[1]->pT() > 25.) {
if (Nb==0 and HT>1000 and met>250 )
_counters.at("ISR1").add_event(event);
if (Nb>=2 and HT>1100 )
_counters.at("ISR2").add_event(event);
if (Nb==0 and met>500 )
_counters.at("ISR3").add_event(event);
if (Nb>=2 and met>300 )
_counters.at("ISR4").add_event(event);
if (Nb==0 and met>250 and MTmin>120 )
_counters.at("ISR5").add_event(event);
if (Nb>=2 and met>200 and MTmin>120 )
_counters.at("ISR6").add_event(event);
if (Nj>=8)
_counters.at("ISR7").add_event(event);
if (Nj>=6 and MTmin>120)
_counters.at("ISR8").add_event(event);
if (Nb>=3 and HT>800)
_counters.at("ISR9").add_event(event);
}
// LL: exactly 2 leptons, both with PT<25 GeV, and MET>50 GeV
if (leptons.size()==2 and leptons[0]->pT() < 25. and leptons[1]->pT() < 25.) {
if (HT>700)
_counters.at("ISR10").add_event(event);
if (met>200)
_counters.at("ISR11").add_event(event);
if (Nj>=6)
_counters.at("ISR12").add_event(event);
if (Nb>=3)
_counters.at("ISR13").add_event(event);
}
// LM: exactly 2 leptons, both with PT>25 GeV, and MET<50 GeV
if (leptons.size()==2 and leptons[0]->pT() > 25. and leptons[1]->pT() > 25.) {
if (Nb==0 and HT>1200 and met<50 )
_counters.at("ISR14").add_event(event);
if (Nb>=2 and HT>1000 and met<50 )
_counters.at("ISR15").add_event(event);
}
// ML: >=3 leptons, at least one with PT>25 GeV, and MET>50 GeV
if (leptons.size()>=3 and leptons[0]->pT() > 25.) {
if (Nb==0 and HT>1000 and met>300 )
_counters.at("ISR16").add_event(event);
if (Nb>=2 and HT>1000)
_counters.at("ISR17").add_event(event);
}
// bool pp = leptons[SS1]->pid()>0;
// bool met_50_200 = met>50 and met<200;
// bool met_200_300 = met>200 and met<300;
// bool met_300_500 = met>300 and met<500;
// bool met_500 = met>500;
// bool MTmin_l_120 = MTmin<120;
// bool MTmin_g_120 = MTmin>120;
// bool Nj_2_4 = Nj>=2 and Nj<=4;
// bool Nj_5 = Nj>=5;
// bool HT_300 = HT<300;
// bool HT_300_1125 = HT>300 and HT<1125;
// bool HT_1125_1300 = HT>1125 and HT<1300;
// bool HT_1300_1600 = HT>1300 and HT<1600;
// bool HT_1300 = HT>1300;
// bool HT_1600 = HT>1600;
// bool SSHH_combine = ((MTmin_l_120&&((met_50_200&&Nj_5)||met_200_300)) or (MTmin_g_120&&met<300) )&&HT_300 ;
// // SSHH: exactly 2 leptons, both with PT>25 GeV, and MET>50 GeV
// if (leptons.size()==2 and leptons[1]->pT() > 25.) {
// if (Nb==0) {
// if (MTmin_l_120 and met_50_200 and Nj_2_4 and HT_300) _counters.at("SRHH-0").add_event(event);
// if (MTmin_l_120 and met_50_200 and Nj_2_4 and HT_300_1125) _counters.at("SRHH-1").add_event(event);
// if (SSHH_combine) _counters.at("SRHH-2").add_event(event);
// if (MTmin_l_120 and met_50_200 and Nj_5 and HT_300_1125) _counters.at("SRHH-3").add_event(event);
// if (MTmin_l_120 and met_200_300 and Nj_2_4 and HT_300_1125 and pp) _counters.at("SRHH-4").add_event(event);
// if (MTmin_l_120 and met_200_300 and Nj_2_4 and HT_300_1125 and !pp) _counters.at("SRHH-5").add_event(event);
// if (MTmin_l_120 and met_200_300 and Nj_5 and HT_300_1125) _counters.at("SRHH-6").add_event(event);
// if (MTmin_g_120 and met_50_200 and Nj_2_4 and HT_300_1125 and pp) _counters.at("SRHH-7").add_event(event);
// if (MTmin_g_120 and met_50_200 and Nj_2_4 and HT_300_1125 and !pp) _counters.at("SRHH-8").add_event(event);
// if (MTmin_g_120 and((met_50_200&&Nj_5)||met_200_300)and HT_300_1125) _counters.at("SRHH-9").add_event(event);
// } else if (Nb==1) {
// if (MTmin_l_120 and met_50_200 and Nj_2_4 and HT_300) _counters.at("SRHH-10").add_event(event);
// if (MTmin_l_120 and met_50_200 and Nj_2_4 and HT_300_1125) _counters.at("SRHH-11").add_event(event);
// if (SSHH_combine and pp) _counters.at("SRHH-12").add_event(event);
// if (SSHH_combine and !pp) _counters.at("SRHH-13").add_event(event);
// if (MTmin_l_120 and met_50_200 and Nj_5 and HT_300_1125 and pp) _counters.at("SRHH-14").add_event(event);
// if (MTmin_l_120 and met_50_200 and Nj_5 and HT_300_1125 and !pp) _counters.at("SRHH-15").add_event(event);
// if (MTmin_l_120 and met_200_300 and Nj_2_4 and HT_300_1125 and pp) _counters.at("SRHH-16").add_event(event);
// if (MTmin_l_120 and met_200_300 and Nj_2_4 and HT_300_1125 and !pp) _counters.at("SRHH-17").add_event(event);
// if (MTmin_l_120 and met_200_300 and Nj_5 and HT_300_1125) _counters.at("SRHH-18").add_event(event);
// if (MTmin_g_120 and met_50_200 and Nj_2_4 and HT_300_1125 and pp) _counters.at("SRHH-19").add_event(event);
// if (MTmin_g_120 and met_50_200 and Nj_2_4 and HT_300_1125 and !pp) _counters.at("SRHH-20").add_event(event);
// if (MTmin_g_120 and((met_50_200&&Nj_5)||met_200_300)and HT_300_1125) _counters.at("SRHH-21").add_event(event);
// } else if (Nb==2){
// if (MTmin_l_120 and met_50_200 and Nj_2_4 and HT_300) _counters.at("SRHH-22").add_event(event);
// if (MTmin_l_120 and met_50_200 and Nj_2_4 and HT_300_1125) _counters.at("SRHH-23").add_event(event);
// if (SSHH_combine and pp) _counters.at("SRHH-24").add_event(event);
// if (SSHH_combine and !pp) _counters.at("SRHH-25").add_event(event);
// if (MTmin_l_120 and met_50_200 and Nj_5 and HT_300_1125 and pp) _counters.at("SRHH-26").add_event(event);
// if (MTmin_l_120 and met_50_200 and Nj_5 and HT_300_1125 and !pp) _counters.at("SRHH-27").add_event(event);
// if (MTmin_l_120 and met_200_300 and Nj_2_4 and HT_300_1125 and pp) _counters.at("SRHH-28").add_event(event);
// if (MTmin_l_120 and met_200_300 and Nj_2_4 and HT_300_1125 and !pp) _counters.at("SRHH-29").add_event(event);
// if (MTmin_l_120 and met_200_300 and Nj_5 and HT_300_1125) _counters.at("SRHH-30").add_event(event);
// if (MTmin_g_120 and met_50_200 and Nj_2_4 and HT_300_1125 and pp) _counters.at("SRHH-31").add_event(event);
// if (MTmin_g_120 and met_50_200 and Nj_2_4 and HT_300_1125 and !pp) _counters.at("SRHH-32").add_event(event);
// if (MTmin_g_120 and((met_50_200&&Nj_5)||met_200_300)and HT_300_1125) _counters.at("SRHH-33").add_event(event);
// } else if (Nb>=3){
// if (MTmin_l_120 and met<300 and HT_300 and pp) _counters.at("SRHH-34").add_event(event);
// if (MTmin_l_120 and met<300 and HT_300 and !pp) _counters.at("SRHH-35").add_event(event);
// if (MTmin_l_120 and met<300 and Nj_2_4 and HT_300_1125 and pp) _counters.at("SRHH-36").add_event(event);
// if (MTmin_l_120 and met<300 and Nj_2_4 and HT_300_1125 and !pp) _counters.at("SRHH-37").add_event(event);
// if (MTmin_l_120 and met<300 and Nj_5 and HT_300_1125 and pp) _counters.at("SRHH-38").add_event(event);
// if (MTmin_l_120 and met<300 and Nj_5 and HT_300_1125 and !pp) _counters.at("SRHH-39").add_event(event);
// if (MTmin_g_120 and HT_300) _counters.at("SRHH-40").add_event(event);
// if (MTmin_g_120 and met<300 and Nj_2_4 and HT_300_1125 and pp) _counters.at("SRHH-41").add_event(event);
// if (MTmin_g_120 and met<300 and Nj_2_4 and HT_300_1125 and !pp) _counters.at("SRHH-42").add_event(event);
// if (MTmin_g_120 and met<300 and Nj_5 and HT_300_1125 and pp) _counters.at("SRHH-43").add_event(event);
// if (MTmin_g_120 and met<300 and Nj_5 and HT_300_1125 and !pp) _counters.at("SRHH-44").add_event(event);
// }
// if (met_300_500 and Nj_2_4 and HT>300 and pp) _counters.at("SRHH-45").add_event(event);
// if (met_300_500 and Nj_2_4 and HT>300 and !pp) _counters.at("SRHH-46").add_event(event);
// if (met_500 and Nj_2_4 and HT>300 and pp) _counters.at("SRHH-47").add_event(event);
// if (met_500 and Nj_2_4 and HT>300 and !pp) _counters.at("SRHH-48").add_event(event);
// if (met_300_500 and Nj_5 and HT>300 and pp) _counters.at("SRHH-49").add_event(event);
// if (met_300_500 and Nj_5 and HT>300 and !pp) _counters.at("SRHH-50").add_event(event);
// if (met_500 and Nj_5 and HT>300 and pp) _counters.at("SRHH-51").add_event(event);
// if (met_500 and Nj_5 and HT>300 and !pp) _counters.at("SRHH-52").add_event(event);
// if (HT_1125_1300 and met<300 and Nj<5) _counters.at("SRHH-53").add_event(event);
// if (HT_1300_1600 and met<300 and Nj<5) _counters.at("SRHH-54").add_event(event);
// if (HT_1600 and met<300 and Nj<5) _counters.at("SRHH-55").add_event(event);
// if (HT_1125_1300 and met<300 and Nj==5 and Nj==6) _counters.at("SRHH-56").add_event(event);
// if (HT_1300_1600 and met<300 and Nj==5 and Nj==6) _counters.at("SRHH-57").add_event(event);
// if (HT_1600 and met<300 and Nj==5 and Nj==6) _counters.at("SRHH-58").add_event(event);
// if (HT_1125_1300 and met<300 and Nj<5) _counters.at("SRHH-59").add_event(event);
// if (HT_1300_1600 and met<300 and Nj<5) _counters.at("SRHH-60").add_event(event);
// if (HT_1600 and met<300 and Nj<5) _counters.at("SRHH-61").add_event(event);
// }
// bool SSHL_combine = MTmin_l_120&&( (met_50_200&&Nj_5) or met_200_300 )&&HT_300 ;
// // SSHL: exactly 2 leptons, one with PT>25 GeV, one with PT<25 GeV, and MET>50 GeV
// if (leptons.size()==2 and leptons[0]->pT() > 25. and leptons[1]->pT() < 25.) {
// if (Nb==0 and MTmin_l_120) {
// if ( met_50_200 and Nj_2_4 and HT_300) _counters.at("SRHL-0").add_event(event);
// if ( met_50_200 and Nj_2_4 and HT_300_1125) _counters.at("SRHL-1").add_event(event);
// if ( SSHL_combine) _counters.at("SRHL-2").add_event(event);
// if ( met_50_200 and Nj_5 and HT_300_1125) _counters.at("SRHL-3").add_event(event);
// if ( met_200_300 and Nj_2_4 and HT_300_1125 and pp) _counters.at("SRHL-4").add_event(event);
// if ( met_200_300 and Nj_2_4 and HT_300_1125 and !pp) _counters.at("SRHL-5").add_event(event);
// if ( met_200_300 and Nj_5 and HT_300_1125) _counters.at("SRHL-6").add_event(event);
// } else if(Nb==1 and MTmin_l_120) {
// if ( met_50_200 and Nj_2_4 and HT_300) _counters.at("SRHL-7").add_event(event);
// if ( met_50_200 and Nj_2_4 and HT_300_1125) _counters.at("SRHL-8").add_event(event);
// if ( SSHL_combine and pp) _counters.at("SRHL-9").add_event(event);
// if ( SSHL_combine and !pp) _counters.at("SRHL-10").add_event(event);
// if ( met_50_200 and Nj_5 and HT_300_1125 and pp) _counters.at("SRHL-11").add_event(event);
// if ( met_50_200 and Nj_5 and HT_300_1125 and !pp) _counters.at("SRHL-12").add_event(event);
// if ( met_200_300 and Nj_2_4 and HT_300_1125) _counters.at("SRHL-13").add_event(event);
// if ( met_200_300 and Nj_5 and HT_300_1125 and pp) _counters.at("SRHL-14").add_event(event);
// if ( met_200_300 and Nj_5 and HT_300_1125 and !pp) _counters.at("SRHL-15").add_event(event);
// } else if(Nb==2 and MTmin_l_120) {
// if ( met_50_200 and Nj_2_4 and HT_300) _counters.at("SRHL-16").add_event(event);
// if ( met_50_200 and Nj_2_4 and HT_300_1125) _counters.at("SRHL-17").add_event(event);
// if ( SSHL_combine and pp) _counters.at("SRHL-18").add_event(event);
// if ( SSHL_combine and !pp) _counters.at("SRHL-19").add_event(event);
// if ( met_50_200 and Nj_5 and HT_300_1125 and pp) _counters.at("SRHL-20").add_event(event);
// if ( met_50_200 and Nj_5 and HT_300_1125 and !pp) _counters.at("SRHL-21").add_event(event);
// if ( met_200_300 and Nj_2_4 and HT_300_1125 and pp) _counters.at("SRHL-22").add_event(event);
// if ( met_200_300 and Nj_2_4 and HT_300_1125 and !pp) _counters.at("SRHL-23").add_event(event);
// if ( met_200_300 and Nj_5 and HT_300_1125) _counters.at("SRHL-24").add_event(event);
// }else if(Nb==3 and MTmin_l_120) {
// if ( met_50_200 and HT_300 and pp) _counters.at("SRHL-25").add_event(event);
// if ( met_50_200 and HT_300 and !pp) _counters.at("SRHL-26").add_event(event);
// if ( met_50_200 and HT_300_1125 and pp) _counters.at("SRHL-27").add_event(event);
// if ( met_50_200 and HT_300_1125 and !pp) _counters.at("SRHL-28").add_event(event);
// if ( met_200_300 and HT_300_1125) _counters.at("SRHL-29").add_event(event);
// }
// if (MTmin_g_120 and met<300 and HT_300) _counters.at("SRHL-30").add_event(event);
// if (MTmin_g_120 and met<300 and HT_300_1125) _counters.at("SRHL-31").add_event(event);
// if (met_300_500 and Nj_2_4 and HT>300 and pp) _counters.at("SRHL-32").add_event(event);
// if (met_300_500 and Nj_2_4 and HT>300 and !pp) _counters.at("SRHL-33").add_event(event);
// if (met_500 and Nj_2_4 and HT>300 and pp) _counters.at("SRHL-34").add_event(event);
// if (met_500 and Nj_2_4 and HT>300 and !pp) _counters.at("SRHL-35").add_event(event);
// if (met_300_500 and Nj_5 and HT>300 and pp) _counters.at("SRHL-36").add_event(event);
// if (met_300_500 and Nj_5 and HT>300 and !pp) _counters.at("SRHL-37").add_event(event);
// if (met_500 and Nj_5 and HT>300) _counters.at("SRHL-38").add_event(event);
// if (HT_1125_1300 and pp) _counters.at("SRHL-39").add_event(event);
// if (HT_1125_1300 and !pp) _counters.at("SRHL-40").add_event(event);
// if (HT_1300 and pp) _counters.at("SRHL-41").add_event(event);
// if (HT_1300 and !pp) _counters.at("SRHL-42").add_event(event);
// }
// // SSLL: exactly 2 leptons, both with PT<25 GeV, and MET>50 GeV
// if (leptons.size()==2 and leptons[0]->pT() < 25. and leptons[1]->pT() < 25.) {
// if (HT>400) {
// if (MTmin_l_120) {
// if (Nb==0) {
// if (met_50_200) _counters.at("SRLL-0").add_event(event);
// else _counters.at("SRLL-1").add_event(event);
// } else if (Nb==1) {
// if (met_50_200) _counters.at("SRLL-2").add_event(event);
// else _counters.at("SRLL-3").add_event(event);
// } else if (Nb==2) {
// if (met_50_200) _counters.at("SRLL-4").add_event(event);
// else _counters.at("SRLL-5").add_event(event);
// } else if (Nb>=3) _counters.at("SRLL-6").add_event(event);
// } else _counters.at("SRLL-7").add_event(event);
// }
// }
// // LM: exactly 2 leptons, both with PT>25 GeV, and MET<50 GeV
// if (leptons.size()==2 and leptons[0]->pT() > 25. and leptons[1]->pT() > 25.) {
// if (HT_300_1125) {
// if (Nb==0) {
// if (Nj_2_4) _counters.at("SRLM-0").add_event(event);
// if (Nj_5) _counters.at("SRLM-1").add_event(event);
// } else if (Nb==1){
// if (Nj_2_4) _counters.at("SRLM-2").add_event(event);
// if (Nj_5) _counters.at("SRLM-3").add_event(event);
// } else if (Nb==2){
// if (Nj_2_4) _counters.at("SRLM-4").add_event(event);
// if (Nj_5) _counters.at("SRLM-5").add_event(event);
// } else _counters.at("SRLM-6").add_event(event);
// } else if (HT_1125_1300){
// if (Nj_2_4) {
// _counters.at("SRLM-7").add_event(event);
// } else if (Nj_5) {
// _counters.at("SRLM-8").add_event(event);
// }
// } else if (HT_1300){
// if (Nj_2_4) {
// _counters.at("SRLM-9").add_event(event);
// } else if (Nj_5) {
// _counters.at("SRLM-1").add_event(event);
// }
// }
// }
// // ML: >=3 leptons, at least one with PT>25 GeV, and MET>50 GeV
// if (leptons.size()>=3 and leptons[0]->pT() > 25.) {
// bool on_Z=false;
// size_t Z_l1;
// size_t Z_l2;
// for (size_t i=0; i<leptons.size(); ++i) {
// for (size_t j=i+1; j<leptons.size(); ++j) {
// if (leptons[i]->pid()+leptons[j]->pid()==0 and fabs((leptons[i]->mom()+leptons[j]->mom()).m()-91)<15) {
// on_Z = true;
// Z_l1=i;
// Z_l2=j;
// break;
// }
// }
// if (on_Z) break;
// }
// if (not on_Z){ // off-Z
// if (Nb==0) {
// if (HT<400) {
// if (met<150) {
// if (MTmin_l_120) _counters.at("SRML-0").add_event(event);
// else _counters.at("SRML-1").add_event(event);
// } else if (met<300) {
// if (MTmin_l_120) _counters.at("SRML-2").add_event(event);
// else _counters.at("SRML-3").add_event(event);
// }
// }
// else if (HT<600) {
// if (met<150) _counters.at("SRML-4").add_event(event);
// else if (met<300) _counters.at("SRML-5").add_event(event);
// }
// } else if (Nb==1) {
// if (HT<400) {
// if (met<150) _counters.at("SRML-6").add_event(event);
// else if (met<300) _counters.at("SRML-7").add_event(event);
// }
// else if (HT<600) {
// if (met<150) _counters.at("SRML-8").add_event(event);
// else if (met<300) _counters.at("SRML-9").add_event(event);
// }
// } else if (Nb==2) {
// if (HT<400) {
// if (met<150) _counters.at("SRML-10").add_event(event);
// else if (met<300) _counters.at("SRML-11").add_event(event);
// }
// else if (HT<600) {
// if (met<150) _counters.at("SRML-12").add_event(event);
// else if (met<300) _counters.at("SRML-13").add_event(event);
// }
// } else if (Nb>=3) {
// if (HT<600) _counters.at("SRML-14").add_event(event);
// }
// if (HT>=600) {
// if (met<150) {
// if (MTmin_l_120) _counters.at("SRML-15").add_event(event);
// else _counters.at("SRML-16").add_event(event);
// } else if (met<300) {
// if (MTmin_l_120) _counters.at("SRML-17").add_event(event);
// else _counters.at("SRML-18").add_event(event);
// }
// }
// if (met>300) {
// if (MTmin_l_120) _counters.at("SRML-19").add_event(event);
// else _counters.at("SRML-20").add_event(event);
// }
// } else { // on-Z
// // M_T^{miss} for on-Z
// double MTmin_Z = 1000;
// for (size_t i=0; i<leptons.size(); ++i) {
// if ( (i != Z_l1) and (i != Z_l2)) {
// double MTmin_try = sqrt(2.*leptons[i]->pT()*met*(1. - cos(leptons[i]->mom().deltaPhi(ptot))));
// if (MTmin_try<MTmin_Z) MTmin_Z=MTmin_try;
// }
// }
// bool MTmin_Z_l_120 = MTmin_Z<120;
// if (Nb==0) {
// if (HT<400) {
// if (met<150) {
// if (MTmin_Z_l_120) _counters.at("SRML-21").add_event(event);
// else _counters.at("SRML-22").add_event(event);
// } else if (met<300) {
// if (MTmin_Z_l_120) _counters.at("SRML-23").add_event(event);
// else _counters.at("SRML-24").add_event(event);
// }
// }
// else if (HT<600) {
// if (met<150) {
// if (MTmin_Z_l_120) _counters.at("SRML-25").add_event(event);
// else _counters.at("SRML-26").add_event(event);
// } else if (met<300) {
// if (MTmin_Z_l_120) _counters.at("SRML-27").add_event(event);
// else _counters.at("SRML-28").add_event(event);
// }
// }
// } else if (Nb==1) {
// if (HT<400) {
// if (met<150) _counters.at("SRML-29").add_event(event);
// else if (met<300) _counters.at("SRML-30").add_event(event);
// }
// else if (HT<600) {
// if (met<150) _counters.at("SRML-31").add_event(event);
// else if (met<300) _counters.at("SRML-32").add_event(event);
// }
// } else if (Nb==2) {
// if (HT<400) {
// if (met<150) _counters.at("SRML-33").add_event(event);
// else if (met<300) _counters.at("SRML-34").add_event(event);
// }
// else if (HT<600) {
// if (met<150) _counters.at("SRML-35").add_event(event);
// else if (met<300) _counters.at("SRML-36").add_event(event);
// }
// } else if (Nb>=3) {
// if (HT<600) _counters.at("SRML-37").add_event(event);
// }
// if (HT>=600) {
// if (met<150) {
// if (MTmin_Z_l_120) _counters.at("SRML-38").add_event(event);
// else _counters.at("SRML-39").add_event(event);
// } else if (met<300) {
// if (MTmin_Z_l_120) _counters.at("SRML-40").add_event(event);
// else _counters.at("SRML-41").add_event(event);
// }
// }
// if (met>300) {
// if (MTmin_Z_l_120) _counters.at("SRML-42").add_event(event);
// else _counters.at("SRML-43").add_event(event);
// }
// }
// }
return;
}
/// Combine the variables of another copy of this analysis (typically on another thread) into this one.
void combine(const Analysis* other)
{
const Analysis_CMS_13TeV_2SSLEP_Stop_137invfb* specificOther
= dynamic_cast<const Analysis_CMS_13TeV_2SSLEP_Stop_137invfb*>(other);
for (auto& pair : _counters) { pair.second += specificOther->_counters.at(pair.first); }
}
void collect_results() {
add_result(SignalRegionData(_counters.at("ISR1"), 16, {12.7, 7.4}));
add_result(SignalRegionData(_counters.at("ISR2"), 14, {11.0, 3.8}));
add_result(SignalRegionData(_counters.at("ISR3"), 13, {10.4, 9.7}));
add_result(SignalRegionData(_counters.at("ISR4"), 17, {11.4, 3.8}));
add_result(SignalRegionData(_counters.at("ISR5"), 10, {6.6, 5.7}));
add_result(SignalRegionData(_counters.at("ISR6"), 8, {6.3, 1.3}));
add_result(SignalRegionData(_counters.at("ISR7"), 12, {7.0, 2.8}));
add_result(SignalRegionData(_counters.at("ISR8"), 10, {6.2, 1.4}));
add_result(SignalRegionData(_counters.at("ISR9"), 8, {7.8, 3.5}));
add_result(SignalRegionData(_counters.at("ISR10"), 12, {10.4, 9.0}));
add_result(SignalRegionData(_counters.at("ISR11"), 13, {12.1, 5.6}));
add_result(SignalRegionData(_counters.at("ISR12"), 7, {7.1, 4.3}));
add_result(SignalRegionData(_counters.at("ISR13"), 3, {1.61, 0.39}));
add_result(SignalRegionData(_counters.at("ISR14"), 3, {3.6, 3.6}));
add_result(SignalRegionData(_counters.at("ISR15"), 4, {2.34, 0.51}));
add_result(SignalRegionData(_counters.at("ISR16"), 7, {5.6, 1.6}));
add_result(SignalRegionData(_counters.at("ISR17"), 7, {5.7, 1.9}));
// HH
// add_result(SignalRegionData(_counters.at("SRHH-0"), 1609, {1510, 310}));
// add_result(SignalRegionData(_counters.at("SRHH-1"), 647, {590, 90}));
// add_result(SignalRegionData(_counters.at("SRHH-2"), 132, {103, 22}));
// add_result(SignalRegionData(_counters.at("SRHH-3"), 51, {38, 7}));
// add_result(SignalRegionData(_counters.at("SRHH-4"), 49, {57, 10}));
// add_result(SignalRegionData(_counters.at("SRHH-5"), 23, {32, 9}));
// add_result(SignalRegionData(_counters.at("SRHH-6"), 7, {5.5, 1.7}));
// add_result(SignalRegionData(_counters.at("SRHH-7"), 31, {25, 6}));
// add_result(SignalRegionData(_counters.at("SRHH-8"), 20, {21, 5}));
// add_result(SignalRegionData(_counters.at("SRHH-9"), 11, {9.4, 1.9}));
// add_result(SignalRegionData(_counters.at("SRHH-10"), 1068, {930, 230}));
// add_result(SignalRegionData(_counters.at("SRHH-11"), 370, {330, 70}));
// add_result(SignalRegionData(_counters.at("SRHH-12"), 38, {36, 7}));
// add_result(SignalRegionData(_counters.at("SRHH-13"), 31, {25, 5}));
// add_result(SignalRegionData(_counters.at("SRHH-14"), 63, {44, 7}));
// add_result(SignalRegionData(_counters.at("SRHH-15"), 38, {39, 8}));
// add_result(SignalRegionData(_counters.at("SRHH-16"), 30, {27, 5}));
// add_result(SignalRegionData(_counters.at("SRHH-17"), 15, {14.8, 3.2}));
// add_result(SignalRegionData(_counters.at("SRHH-18"), 12, {11.5, 3.0}));
// add_result(SignalRegionData(_counters.at("SRHH-19"), 14, {11.8, 2.6}));
// add_result(SignalRegionData(_counters.at("SRHH-20"), 16, {9.6, 2.1}));
// add_result(SignalRegionData(_counters.at("SRHH-21"), 15, {10.0, 1.6}));
// add_result(SignalRegionData(_counters.at("SRHH-22"), 345, {270, 40}));
// add_result(SignalRegionData(_counters.at("SRHH-23"), 169, {143, 20}));
// add_result(SignalRegionData(_counters.at("SRHH-24"), 11, {15.2, 2.4}));
// add_result(SignalRegionData(_counters.at("SRHH-25"), 18, {13.8, 3.4}));
// add_result(SignalRegionData(_counters.at("SRHH-26"), 43, {33, 5}));
// add_result(SignalRegionData(_counters.at("SRHH-27"), 38, {29, 4}));
// add_result(SignalRegionData(_counters.at("SRHH-28"), 9, {11.5, 2.5}));
// add_result(SignalRegionData(_counters.at("SRHH-29"), 5, {6.7, 1.2}));
// add_result(SignalRegionData(_counters.at("SRHH-30"), 6, {7.5, 1.8}));
// add_result(SignalRegionData(_counters.at("SRHH-31"), 14, {5.9, 1.0}));
// add_result(SignalRegionData(_counters.at("SRHH-32"), 7, {6.5, 1.9}));
// add_result(SignalRegionData(_counters.at("SRHH-33"), 11, {6.7, 1.2}));
// add_result(SignalRegionData(_counters.at("SRHH-34"), 17, {10.3, 1.9}));
// add_result(SignalRegionData(_counters.at("SRHH-35"), 11, {8.6, 1.7}));
// add_result(SignalRegionData(_counters.at("SRHH-36"), 6, {10.6, 2.0}));
// add_result(SignalRegionData(_counters.at("SRHH-37"), 5, {7.3, 1.3}));
// add_result(SignalRegionData(_counters.at("SRHH-38"), 8, {9.6, 2.2}));
// add_result(SignalRegionData(_counters.at("SRHH-39"), 11, {9.2, 1.9}));
// add_result(SignalRegionData(_counters.at("SRHH-40"), 2, {1.3, 0.6}));
// add_result(SignalRegionData(_counters.at("SRHH-41"), 1, {0.6, 0.4}));
// add_result(SignalRegionData(_counters.at("SRHH-42"), 0, {0.8, 0.4}));
// add_result(SignalRegionData(_counters.at("SRHH-43"), 1, {0.7, 0.4}));
// add_result(SignalRegionData(_counters.at("SRHH-44"), 1, {0.7, 0.5}));
// add_result(SignalRegionData(_counters.at("SRHH-45"), 59, {42, 7}));
// add_result(SignalRegionData(_counters.at("SRHH-46"), 23, {18, 4}));
// add_result(SignalRegionData(_counters.at("SRHH-47"), 10, {13, 9}));
// add_result(SignalRegionData(_counters.at("SRHH-48"), 4, {2.0, 0.5}));
// add_result(SignalRegionData(_counters.at("SRHH-49"), 13, {6.3, 1.0}));
// add_result(SignalRegionData(_counters.at("SRHH-50"), 4, {3.7, 0.7}));
// add_result(SignalRegionData(_counters.at("SRHH-51"), 4, {1.26, 0.33}));
// add_result(SignalRegionData(_counters.at("SRHH-52"), 2, {0.4, 0.4}));
// add_result(SignalRegionData(_counters.at("SRHH-53"), 24, {10.1, 1.5}));
// add_result(SignalRegionData(_counters.at("SRHH-54"), 4, {7.0, 1.1}));
// add_result(SignalRegionData(_counters.at("SRHH-55"), 5, {4.3, 0.9}));
// add_result(SignalRegionData(_counters.at("SRHH-56"), 7, {5.3, 0.8}));
// add_result(SignalRegionData(_counters.at("SRHH-57"), 6, {6, 6}));
// add_result(SignalRegionData(_counters.at("SRHH-58"), 3, {2.2, 0.4}));
// add_result(SignalRegionData(_counters.at("SRHH-59"), 5, {1.8, 0.5}));
// add_result(SignalRegionData(_counters.at("SRHH-60"), 4, {1.9, 0.4}));
// add_result(SignalRegionData(_counters.at("SRHH-61"), 0, {1.3, 0.9}));
// // HL
// add_result(SignalRegionData(_counters.at("SRHL-0"), 1504, {1300, 310}));
// add_result(SignalRegionData(_counters.at("SRHL-1"), 319, {310, 70}));
// add_result(SignalRegionData(_counters.at("SRHL-2"), 32, {25, 6}));
// add_result(SignalRegionData(_counters.at("SRHL-3"), 32, {32, 8}));
// add_result(SignalRegionData(_counters.at("SRHL-4"), 32, {29, 6}));
// add_result(SignalRegionData(_counters.at("SRHL-5"), 11, {17, 6}));
// add_result(SignalRegionData(_counters.at("SRHL-6"), 6, {4.5, 2.5}));
// add_result(SignalRegionData(_counters.at("SRHL-7"), 1223, {1010, 250}));
// add_result(SignalRegionData(_counters.at("SRHL-8"), 307, {270, 60}));
// add_result(SignalRegionData(_counters.at("SRHL-9"), 5, {7.1, 1.7}));
// add_result(SignalRegionData(_counters.at("SRHL-10"), 7, {6.5, 1.6}));
// add_result(SignalRegionData(_counters.at("SRHL-11"), 42, {39, 9}));
// add_result(SignalRegionData(_counters.at("SRHL-12"), 37, {31, 8}));
// add_result(SignalRegionData(_counters.at("SRHL-13"), 27, {23, 5}));
// add_result(SignalRegionData(_counters.at("SRHL-14"), 7, {2.1, 1.1}));
// add_result(SignalRegionData(_counters.at("SRHL-15"), 2, {1.7, 0.9}));
// add_result(SignalRegionData(_counters.at("SRHL-16"), 256, {210, 40}));
// add_result(SignalRegionData(_counters.at("SRHL-17"), 104, {85, 14}));
// add_result(SignalRegionData(_counters.at("SRHL-18"), 4, {2.5, 1.2}));
// add_result(SignalRegionData(_counters.at("SRHL-19"), 3, {3.0, 1.5}));
// add_result(SignalRegionData(_counters.at("SRHL-20"), 27, {18.9, 3.5}));
// add_result(SignalRegionData(_counters.at("SRHL-21"), 18, {15.9, 2.8}));
// add_result(SignalRegionData(_counters.at("SRHL-22"), 2, {3.3, 0.6}));
// add_result(SignalRegionData(_counters.at("SRHL-23"), 2, {4.4, 1.6}));
// add_result(SignalRegionData(_counters.at("SRHL-24"), 5, {4.5, 1.7}));
// add_result(SignalRegionData(_counters.at("SRHL-25"), 8, {8.2, 2.2}));
// add_result(SignalRegionData(_counters.at("SRHL-26"), 6, {8.1, 2.2}));
// add_result(SignalRegionData(_counters.at("SRHL-27"), 12, {9.7, 2.1}));
// add_result(SignalRegionData(_counters.at("SRHL-28"), 7, {10.8, 2.8}));
// add_result(SignalRegionData(_counters.at("SRHL-29"), 3, {1.1, 0.4}));
// add_result(SignalRegionData(_counters.at("SRHL-30"), 5, {2.2, 0.5}));
// add_result(SignalRegionData(_counters.at("SRHL-31"), 3, {2.6, 0.5}));
// add_result(SignalRegionData(_counters.at("SRHL-32"), 23, {22, 6}));
// add_result(SignalRegionData(_counters.at("SRHL-33"), 8, {7.2, 1.4}));
// add_result(SignalRegionData(_counters.at("SRHL-34"), 4, {2.3, 0.5}));
// add_result(SignalRegionData(_counters.at("SRHL-35"), 1, {0.42, 0.33}));
// add_result(SignalRegionData(_counters.at("SRHL-36"), 3, {3.2, 1.5}));
// add_result(SignalRegionData(_counters.at("SRHL-37"), 0, {1.4, 0.6}));
// add_result(SignalRegionData(_counters.at("SRHL-38"), 0, {0.41, 0.25}));
// add_result(SignalRegionData(_counters.at("SRHL-39"), 7, {3.1, 0.7}));
// add_result(SignalRegionData(_counters.at("SRHL-40"), 0, {4, 4}));
// add_result(SignalRegionData(_counters.at("SRHL-41"), 8, {4.7, 0.9}));
// add_result(SignalRegionData(_counters.at("SRHL-42"), 6, {1.71, 0.35}));
// // LL
// add_result(SignalRegionData(_counters.at("SRLL-0"), 25, {20, 6}));
// add_result(SignalRegionData(_counters.at("SRLL-1"), 6, {4.9, 1.5}));
// add_result(SignalRegionData(_counters.at("SRLL-2"), 29, {23, 6}));
// add_result(SignalRegionData(_counters.at("SRLL-3"), 8, {4, 4}));
// add_result(SignalRegionData(_counters.at("SRLL-4"), 13, {7.7, 2.1}));
// add_result(SignalRegionData(_counters.at("SRLL-5"), 0, {1.5, 0.8}));
// add_result(SignalRegionData(_counters.at("SRLL-6"), 0, {2.5, 1.2}));
// add_result(SignalRegionData(_counters.at("SRLL-7"), 0, {0.07, 0.07}));
// // LM
// add_result(SignalRegionData(_counters.at("SRLM-0"), 314, {240, 50}));
// add_result(SignalRegionData(_counters.at("SRLM-1"), 22, {20, 5}));
// add_result(SignalRegionData(_counters.at("SRLM-2"), 159, {140, 31}));
// add_result(SignalRegionData(_counters.at("SRLM-3"), 37, {32, 7}));
// add_result(SignalRegionData(_counters.at("SRLM-4"), 66, {54, 8}));
// add_result(SignalRegionData(_counters.at("SRLM-5"), 33, {22, 4}));
// add_result(SignalRegionData(_counters.at("SRLM-6"), 23, {9.7, 2.1}));
// add_result(SignalRegionData(_counters.at("SRLM-7"), 3, {1.5, 0.5}));
// add_result(SignalRegionData(_counters.at("SRLM-8"), 1, {1.6, 0.4}));
// add_result(SignalRegionData(_counters.at("SRLM-9"), 1, {2.9, 2.9}));
// add_result(SignalRegionData(_counters.at("SRLM-10"), 3, {1.9, 1.4}));
// // ML
// add_result(SignalRegionData(_counters.at("SRML-0"), 263, {220, 40}));
// add_result(SignalRegionData(_counters.at("SRML-1"), 1, {2.5, 2.5}));
// add_result(SignalRegionData(_counters.at("SRML-2"), 34, {32, 6}));
// add_result(SignalRegionData(_counters.at("SRML-3"), 1, {0.9, 0.5}));
// add_result(SignalRegionData(_counters.at("SRML-4"), 28, {22, 4}));
// add_result(SignalRegionData(_counters.at("SRML-5"), 8, {9.5, 1.8}));
// add_result(SignalRegionData(_counters.at("SRML-6"), 265, {210, 40}));
// add_result(SignalRegionData(_counters.at("SRML-7"), 47, {36, 6}));
// add_result(SignalRegionData(_counters.at("SRML-8"), 20, {21.6, 3.2}));
// add_result(SignalRegionData(_counters.at("SRML-9"), 16, {11.6, 1.9}));
// add_result(SignalRegionData(_counters.at("SRML-10"), 105, {84, 11}));
// add_result(SignalRegionData(_counters.at("SRML-11"), 17, {15.5, 2.1}));
// add_result(SignalRegionData(_counters.at("SRML-12"), 21, {15.7, 2.2}));
// add_result(SignalRegionData(_counters.at("SRML-13"), 8, {5.3, 0.8}));
// add_result(SignalRegionData(_counters.at("SRML-14"), 12, {10.2, 2.1}));
// add_result(SignalRegionData(_counters.at("SRML-15"), 40, {27, 4}));
// add_result(SignalRegionData(_counters.at("SRML-16"), 2, {0.8, 0.5}));
// add_result(SignalRegionData(_counters.at("SRML-17"), 24, {17.8, 2.4}));
// add_result(SignalRegionData(_counters.at("SRML-18"), 0, {1.0, 0.4}));
// add_result(SignalRegionData(_counters.at("SRML-19"), 30, {17.8, 3.0}));
// add_result(SignalRegionData(_counters.at("SRML-20"), 2, {1.26, 0.33}));
// add_result(SignalRegionData(_counters.at("SRML-21"), 955, {830, 180}));
// add_result(SignalRegionData(_counters.at("SRML-22"), 136, {108, 22}));
// add_result(SignalRegionData(_counters.at("SRML-23"), 139, {117, 26}));
// add_result(SignalRegionData(_counters.at("SRML-24"), 8, {11.1, 2.3}));
// add_result(SignalRegionData(_counters.at("SRML-25"), 128, {111, 24}));
// add_result(SignalRegionData(_counters.at("SRML-26"), 20, {21, 5}));
// add_result(SignalRegionData(_counters.at("SRML-27"), 45, {42, 10}));
// add_result(SignalRegionData(_counters.at("SRML-28"), 3, {3.4, 0.9}));
// add_result(SignalRegionData(_counters.at("SRML-29"), 408, {320, 50}));
// add_result(SignalRegionData(_counters.at("SRML-30"), 50, {47, 8}));
// add_result(SignalRegionData(_counters.at("SRML-31"), 62, {51, 9}));
// add_result(SignalRegionData(_counters.at("SRML-32"), 24, {15.1, 2.6}));
// add_result(SignalRegionData(_counters.at("SRML-33"), 157, {131, 24}));
// add_result(SignalRegionData(_counters.at("SRML-34"), 24, {20, 4}));
// add_result(SignalRegionData(_counters.at("SRML-35"), 36, {27, 5}));
// add_result(SignalRegionData(_counters.at("SRML-36"), 11, {7.8, 1.5}));
// add_result(SignalRegionData(_counters.at("SRML-37"), 18, {12.9, 2.6}));
// add_result(SignalRegionData(_counters.at("SRML-38"), 117, {82, 14}));
// add_result(SignalRegionData(_counters.at("SRML-39"), 26, {18, 4}));
// add_result(SignalRegionData(_counters.at("SRML-40"), 29, {39, 8}));
// add_result(SignalRegionData(_counters.at("SRML-41"), 7, {4.9, 0.9}));
// add_result(SignalRegionData(_counters.at("SRML-42"), 44, {46, 10}));
// add_result(SignalRegionData(_counters.at("SRML-43"), 11, {5.7, 1.2}));
// static const vector< vector<double> > BKGCOV = {
// {},
// {}
// };
// set_covariance(BKGCOV);
return;
}
protected:
void analysis_specific_reset() {
for (auto& pair : _counters) { pair.second.reset(); }
}
};
DEFINE_ANALYSIS_FACTORY(CMS_13TeV_2SSLEP_Stop_137invfb)
}
}
Updated on 2024-07-18 at 13:53:35 +0000