file analyses/Analysis_ATLAS_13TeV_MultiLEP_confnote_36invfb.cpp

[No description available]

Namespaces

Name
Gambit
TODO: see if we can use this one:
Gambit::ColliderBit

Classes

Name
classGambit::ColliderBit::Analysis_ATLAS_13TeV_MultiLEP_confnote_36invfb
structGambit::ColliderBit::Analysis_ATLAS_13TeV_MultiLEP_confnote_36invfb::ptComparison
structGambit::ColliderBit::Analysis_ATLAS_13TeV_MultiLEP_confnote_36invfb::ptJetComparison

Source code

///
///  \author Rose Kudzman-Blais
///  \date 2017 May
///
///  *********************************************

// Based on https://cds.cern.ch/record/2267406

#include <vector>
#include <cmath>
#include <memory>
#include <iomanip>
#include <algorithm>
#include <fstream>

#include "gambit/ColliderBit/analyses/Analysis.hpp"
#include "gambit/ColliderBit/ATLASEfficiencies.hpp"
#include "gambit/ColliderBit/mt2_bisect.h"

using namespace std;

namespace Gambit {
  namespace ColliderBit {

    class Analysis_ATLAS_13TeV_MultiLEP_confnote_36invfb : public Analysis {
    private:

      // Numbers passing cuts
      std::map<string, EventCounter> _counters = {
        {"SR2_SF_loose", EventCounter("SR2_SF_loose")},
        {"SR2_SF_tight", EventCounter("SR2_SF_tight")},
        {"SR2_DF_100", EventCounter("SR2_DF_100")},
        {"SR2_DF_150", EventCounter("SR2_DF_150")},
        {"SR2_DF_200", EventCounter("SR2_DF_200")},
        {"SR2_DF_300", EventCounter("SR2_DF_300")},
        {"SR2_int", EventCounter("SR2_int")},
        {"SR2_high", EventCounter("SR2_high")},
        {"SR2_low", EventCounter("SR2_low")},
        {"SR3_slep_a", EventCounter("SR3_slep_a")},
        {"SR3_slep_b", EventCounter("SR3_slep_b")},
        {"SR3_slep_c", EventCounter("SR3_slep_c")},
        {"SR3_slep_d", EventCounter("SR3_slep_d")},
        {"SR3_slep_e", EventCounter("SR3_slep_e")},
        {"SR3_WZ_0Ja", EventCounter("SR3_WZ_0Ja")},
        {"SR3_WZ_0Jb", EventCounter("SR3_WZ_0Jb")},
        {"SR3_WZ_0Jc", EventCounter("SR3_WZ_0Jc")},
        {"SR3_WZ_1Ja", EventCounter("SR3_WZ_1Ja")},
        {"SR3_WZ_1Jb", EventCounter("SR3_WZ_1Jb")},
        {"SR3_WZ_1Jc", EventCounter("SR3_WZ_1Jc")},
      };

      vector<int> cutFlowVector1;
      vector<string> cutFlowVector1_str;
      size_t NCUTS1;
      // vector<double> cutFlowVector1ATLAS_200_100;
      // double xsec1ATLAS_200_100;

      vector<int> cutFlowVector2;
      vector<string> cutFlowVector2_str;
      size_t NCUTS2;
      // vector<double> cutFlowVector2ATLAS_400_200;
      // double xsec2ATLAS_400_200;
      // vector<double> cutFlowVector2ATLAS_500_100;
      // double xsec2ATLAS_500_100;

      vector<int> cutFlowVector3;
      vector<string> cutFlowVector3_str;
      size_t NCUTS3;
      // vector<double> cutFlowVector3ATLAS_200_100;
      // double xsec3ATLAS_200_100;

      vector<int> cutFlowVector4;
      vector<string> cutFlowVector4_str;
      size_t NCUTS4;
      // vector<double> cutFlowVector4ATLAS_800_600;
      // double xsec4ATLAS_800_600;

      vector<int> cutFlowVector5;
      vector<string> cutFlowVector5_str;
      size_t NCUTS5;
      // vector<double> cutFlowVector5ATLAS_401_1;
      // double xsec5ATLAS_401_1;
      // vector<double> cutFlowVector5ATLAS_300_150;
      // double xsec5ATLAS_300_150;

      // ofstream cutflowFile;


    public:

      // Required detector sim
      static constexpr const char* detector = "ATLAS";

      Analysis_ATLAS_13TeV_MultiLEP_confnote_36invfb() {

        set_analysis_name("ATLAS_13TeV_MultiLEP_confnote_36invfb");
        set_luminosity(36.1);

        NCUTS1=22;
        // xsec1ATLAS_200_100=1807.4;
        for (size_t i=0;i<NCUTS1;i++){
          cutFlowVector1.push_back(0);
          // cutFlowVector1ATLAS_200_100.push_back(0);
          cutFlowVector1_str.push_back("");
        }

        NCUTS2=14;
        // xsec2ATLAS_400_200=121.0269;
        // xsec2ATLAS_500_100=46.3576;
        for (size_t i=0;i<NCUTS2;i++){
          cutFlowVector2.push_back(0);
          // cutFlowVector2ATLAS_400_200.push_back(0);
          // cutFlowVector2ATLAS_500_100.push_back(0);
          cutFlowVector2_str.push_back("");
        }

        NCUTS3=24;
        // xsec3ATLAS_200_100=1807.4;
        for (size_t i=0;i<NCUTS3;i++){
          cutFlowVector3.push_back(0);
          // cutFlowVector3ATLAS_200_100.push_back(0);
          cutFlowVector3_str.push_back("");
        }

        NCUTS4=12;
        // xsec4ATLAS_800_600=3.803;
        for (size_t i=0;i<NCUTS4;i++){
          cutFlowVector4.push_back(0);
          // cutFlowVector4ATLAS_800_600.push_back(0);
          cutFlowVector4_str.push_back("");
        }

        NCUTS5=11;
        // xsec5ATLAS_401_1=5.43;
        // xsec5ATLAS_300_150=190.159;
        for (size_t i=0;i<NCUTS5;i++){
          cutFlowVector5.push_back(0);
          // cutFlowVector5ATLAS_401_1.push_back(0);
          // cutFlowVector5ATLAS_300_150.push_back(0);
          cutFlowVector5_str.push_back("");
        }

      }

      struct ptComparison {
        bool operator() (const HEPUtils::Particle* i,const HEPUtils::Particle* j) {return (i->pT()>j->pT());}
      } comparePt;

      struct ptJetComparison {
        bool operator() (const HEPUtils::Jet* i,const HEPUtils::Jet* j) {return (i->pT()>j->pT());}
      } compareJetPt;

      void run(const HEPUtils::Event* event) {

        double met = event->met();

        // Baseline electrons
        vector<const HEPUtils::Particle*> baselineElectrons;
        for (const HEPUtils::Particle* electron : event->electrons()) {
          if (electron->pT()>10. && electron->abseta()<2.47)baselineElectrons.push_back(electron);
        }

        // Apply electron efficiency
        ATLAS::applyElectronEff(baselineElectrons);

        // Apply loose electron selection
        ATLAS::applyLooseIDElectronSelectionR2(baselineElectrons);

        // Baseline muons
        vector<const HEPUtils::Particle*> baselineMuons;
        for (const HEPUtils::Particle* muon : event->muons()) {
          if (muon->pT()>10. && muon->abseta()<2.4)baselineMuons.push_back(muon);
        }

        // Apply muon efficiency
        ATLAS::applyMuonEff(baselineMuons);

        vector<const HEPUtils::Jet*> baselineJets;
        for (const HEPUtils::Jet* jet : event->jets("antikt_R04")) {
          if (jet->pT()>20. && jet->abseta()<4.5)baselineJets.push_back(jet);
        }

        //Overlap Removal + Signal Objects
        vector<const HEPUtils::Particle*> signalElectrons;
        vector<const HEPUtils::Particle*> signalMuons;
        vector<const HEPUtils::Particle*> signalLeptons;
        vector<const HEPUtils::Jet*> signalJets;
        vector<const HEPUtils::Jet*> signalBJets;

        const vector<double> aBJet={0,10.};
        const vector<double> bBJet={0,30., 40., 50., 70., 80., 90., 100.,150., 200., 10000.};
        const vector<double> cBJet={0.63, 0.705, 0.745, 0.76, 0.775, 0.79,0.795, 0.805, 0.795, 0.76};
        HEPUtils::BinnedFn2D<double> _eff2d(aBJet,bBJet,cBJet);

        vector<const HEPUtils::Jet*> overlapJet;
        for (size_t iJet=0;iJet<baselineJets.size();iJet++) {
          vector<const HEPUtils::Particle*> overlapEl;
          bool hasTag=has_tag(_eff2d, baselineJets.at(iJet)->abseta(), baselineJets.at(iJet)->pT());
          for (size_t iEl=0;iEl<baselineElectrons.size();iEl++) {
            if (baselineElectrons.at(iEl)->mom().deltaR_eta(baselineJets.at(iJet)->mom())<0.2)overlapEl.push_back(baselineElectrons.at(iEl));
          }
          if (overlapEl.size()>0 && (baselineJets.at(iJet)->btag() && hasTag)) {
            for (size_t iO=0;iO<overlapEl.size();iO++) {
              baselineElectrons.erase(remove(baselineElectrons.begin(), baselineElectrons.end(), overlapEl.at(iO)), baselineElectrons.end());
            }
          }
          if (overlapEl.size()>0 && !(baselineJets.at(iJet)->btag() && hasTag))overlapJet.push_back(baselineJets.at(iJet));
        }
        for (size_t iO=0;iO<overlapJet.size();iO++) {
          baselineJets.erase(remove(baselineJets.begin(), baselineJets.end(), overlapJet.at(iO)), baselineJets.end());
        }

        for (size_t iEl=0;iEl<baselineElectrons.size();iEl++) {
          bool overlap=false;
          for (size_t iJet=0;iJet<baselineJets.size();iJet++) {
            if (baselineElectrons.at(iEl)->mom().deltaR_eta(baselineJets.at(iJet)->mom())<0.4)overlap=true;
          }
          if (!overlap)signalElectrons.push_back(baselineElectrons.at(iEl));
        }
        ATLAS::applyMediumIDElectronSelectionR2(signalElectrons);

        for (size_t iJet=0;iJet<baselineJets.size();iJet++) {
          bool overlap=false;
          for (size_t iMu=0;iMu<baselineMuons.size();iMu++) {
            if (baselineMuons.at(iMu)->mom().deltaR_eta(baselineJets.at(iJet)->mom())<0.2 && baselineMuons.at(iMu)->pT()>0.7*baselineJets.at(iJet)->pT())overlap=true;
          }
          if (!overlap) {
            bool hasTag=has_tag(_eff2d, baselineJets.at(iJet)->abseta(), baselineJets.at(iJet)->pT());
            signalJets.push_back(baselineJets.at(iJet));
            if (baselineJets.at(iJet)->btag() && hasTag && baselineJets.at(iJet)->abseta()<2.4)signalBJets.push_back(baselineJets.at(iJet));
          }
        }

        for (size_t iMu=0;iMu<baselineMuons.size();iMu++) {
          bool overlap=false;
          for (size_t iJet=0;iJet<signalJets.size();iJet++) {
            if (baselineMuons.at(iMu)->mom().deltaR_eta(signalJets.at(iJet)->mom())<0.4)overlap=true;
          }
          if (!overlap)signalMuons.push_back(baselineMuons.at(iMu));
        }

        signalLeptons=signalElectrons;
        signalLeptons.insert(signalLeptons.end(),signalMuons.begin(),signalMuons.end());
        sort(signalJets.begin(),signalJets.end(),compareJetPt);
        sort(signalLeptons.begin(),signalLeptons.end(),comparePt);
        size_t nBaselineLeptons=baselineElectrons.size()+baselineMuons.size();
        size_t nSignalLeptons=signalLeptons.size();
        size_t nSignalJets=signalJets.size();
        size_t nSignalBJets=signalBJets.size();

        vector<vector<const HEPUtils::Particle*>> SFOSpairs=getSFOSpairs(signalLeptons);
        vector<vector<const HEPUtils::Particle*>> OSpairs=getOSpairs(signalLeptons);

        //Variables
        double pT_l0=0.;
        double pT_l1=0.;
        double pT_l2=0.;
        // double mlll=0.;
        double pTlll=999.;
        double mll=999.;
        double mT2=0;
        double deltaR_ll=999.;

        double pT_j0=0.;
        double pT_j1=0.;
        double pT_j2=0.;
        double mjj=0;
        double deltaR_jj=999.;

        HEPUtils::P4 Z;
        double deltaPhi_met_Z=999.;

        HEPUtils::P4 W;
        vector<HEPUtils::P4> W_ISR;
        double deltaPhi_met_W=0.;
        double deltaPhi_met_ISR=0.;
        double deltaPhi_met_jet0=0.;

        double mTmin=999;
        double mSFOS=999;

        bool central_jet_veto=true;
        bool bjet_veto=false;

        for (size_t iJet=0;iJet<nSignalJets;iJet++) {
          if (signalJets.at(iJet)->pT()>60 && signalJets.at(iJet)->abseta()<2.4)central_jet_veto=false;
        }
        if (nSignalBJets==0)bjet_veto=true;

        if (nSignalLeptons>0)pT_l0=signalLeptons.at(0)->pT();
        if (nSignalLeptons>1) {
          pT_l1=signalLeptons.at(1)->pT();
          mll=(signalLeptons.at(0)->mom()+signalLeptons.at(1)->mom()).m();
          deltaR_ll=signalLeptons.at(0)->mom().deltaR_eta(signalLeptons.at(1)->mom());

          double pLep1[3] = {signalLeptons.at(0)->mass(), signalLeptons.at(0)->mom().px(), signalLeptons.at(0)->mom().py()};
          double pLep2[3] = {signalLeptons.at(1)->mass(), signalLeptons.at(1)->mom().px(), signalLeptons.at(1)->mom().py()};
          double pMiss[3] = {0., event->missingmom().px(), event->missingmom().py() };
          double mn = 0.;

          mt2_bisect::mt2 mt2_calc;
          mt2_calc.set_momenta(pLep1,pLep2,pMiss);
          mt2_calc.set_mn(mn);
          mT2 = mt2_calc.get_mt2();

          Z=signalLeptons.at(0)->mom()+signalLeptons.at(1)->mom();
          deltaPhi_met_Z=Z.deltaPhi(event->missingmom());
          for (size_t iPa=0;iPa<SFOSpairs.size();iPa++) {
            for (size_t iLep=0;iLep<signalLeptons.size();iLep++) {
              if (signalLeptons.at(iLep)!=SFOSpairs.at(iPa).at(0) && signalLeptons.at(iLep)!=SFOSpairs.at(iPa).at(1)) {
                double mT=sqrt(2*signalLeptons.at(iLep)->pT()*met*(1-cos(signalLeptons.at(iLep)->mom().deltaPhi(event->missingmom()))));
                if (mT<mTmin) {
                  mTmin=mT;
                  mSFOS=(SFOSpairs.at(iPa).at(0)->mom()+SFOSpairs.at(iPa).at(1)->mom()).m();
                }
              }
            }
          }
        }

        if (nSignalLeptons>2) {
          pT_l2=signalLeptons.at(2)->pT();
          // mlll=(signalLeptons.at(0)->mom()+signalLeptons.at(1)->mom()+signalLeptons.at(2)->mom()).m();
          pTlll=(signalLeptons.at(0)->mom()+signalLeptons.at(1)->mom()+signalLeptons.at(2)->mom()).pT();
        }

        if (nSignalJets>0) {
          pT_j0=signalJets.at(0)->pT();
          deltaPhi_met_jet0=signalJets.at(0)->mom().deltaPhi(event->missingmom());
        }
        if (nSignalJets>1) {
          pT_j1=signalJets.at(1)->pT();
          if (nSignalJets<3 && bjet_veto) {
            W=signalJets.at(0)->mom()+signalJets.at(1)->mom();
            mjj=W.m();
            deltaR_jj=signalJets.at(0)->mom().deltaR_eta(signalJets.at(1)->mom());
            deltaPhi_met_W=W.deltaPhi(event->missingmom());
          }
          if (nSignalJets>2 && nSignalJets<6 && nSignalLeptons>1 && bjet_veto) {
            W_ISR=get_W_ISR(signalJets,Z,event->missingmom());
            W=W_ISR.at(0);
            mjj=W.m();
            deltaR_jj=W_ISR.at(3).deltaR_eta(W_ISR.at(2));
            deltaPhi_met_W=W.deltaPhi(event->missingmom());
            deltaPhi_met_ISR=W_ISR.at(1).deltaPhi(event->missingmom());
          }
        }
        if (nSignalJets>2)pT_j2=signalJets.at(2)->pT();

        bool preselection=false;
        if ((nSignalLeptons==2 || nSignalLeptons==3) && nBaselineLeptons==nSignalLeptons && pT_l0>25 && pT_l1>20)preselection=true;

        //Signal Regions
        //2lep+0jet
        if (preselection && nSignalLeptons==2 && OSpairs.size()==1 && mll>40 && central_jet_veto && bjet_veto) {
          if (SFOSpairs.size()==1) {
            if (mT2>100 && mll>111) _counters.at("SR2_SF_loose").add_event(event);
            if (mT2>130 && mll>300) _counters.at("SR2_SF_tight").add_event(event);
          }
          if (SFOSpairs.size()==0) {
            if (mT2>100) _counters.at("SR2_DF_100").add_event(event);
            if (mT2>150) _counters.at("SR2_DF_150").add_event(event);
            if (mT2>200) _counters.at("SR2_DF_200").add_event(event);
            if (mT2>300) _counters.at("SR2_DF_300").add_event(event);
          }
        }

        //2lep+jets
        if (preselection && nSignalLeptons==2 && SFOSpairs.size()==1 && bjet_veto && nSignalJets>1 && pT_j0>30 && pT_j1>30 && pT_l1>25) {
          //SR2_int + SR2_high
          if (mll>81. && mll<101. && mjj>70. && mjj<100. && Z.pT()>80. && W.pT()>100. && mT2>100. && deltaR_jj<1.5 && deltaR_ll<1.8 && deltaPhi_met_W>0.5 && deltaPhi_met_W<3.0) {
            if (met>150) _counters.at("SR2_int").add_event(event);
            if (met>250) _counters.at("SR2_high").add_event(event);
          }
          //SR2_low_2J
          if (nSignalJets==2 && mll>81. && mll<101. && mjj>70. && mjj<90. && met>100. && Z.pT()>60. && deltaPhi_met_Z<0.8 && deltaPhi_met_W>1.5 && (met/Z.pT())>0.6 && (met/Z.pT())<1.6 && (met/W.pT())<0.8) _counters.at("SR2_low").add_event(event);
          //SR2_low_3J
          if (nSignalJets>2 && nSignalJets<6 && mll>86 && mll<96 && mjj>70. && mjj<90. && met>100 && Z.pT()>40 && deltaR_jj<2.2 && deltaPhi_met_W<2.2 && deltaPhi_met_ISR>2.4 && deltaPhi_met_jet0>2.6 && (met/W_ISR.at(1).pT())>0.4 && (met/W_ISR.at(1).pT())<0.8 && Z.abseta()<1.6 && pT_j2>30.) _counters.at("SR2_low").add_event(event);
        }

        //3lep
        if (preselection && nSignalLeptons==3 && bjet_veto && SFOSpairs.size()) {
          if (mSFOS<81.2 && met>130. && mTmin>110.) {
            if (pT_l2>20. && pT_l2<30.) _counters.at("SR3_slep_a").add_event(event);
            if (pT_l2>30.) _counters.at("SR3_slep_b").add_event(event);
          }
          if (mSFOS>101.2 && met>130. && mTmin>110.) {
            if (pT_l2>20. && pT_l2<50.) _counters.at("SR3_slep_c").add_event(event);
            if (pT_l2>50. && pT_l2<80.) _counters.at("SR3_slep_d").add_event(event);
            if (pT_l2>80.) _counters.at("SR3_slep_e").add_event(event);
          }
          if (mSFOS>81.2 && mSFOS<101.2 && nSignalJets==0 && mTmin>110.) {
            if (met>60. && met<120.) _counters.at("SR3_WZ_0Ja").add_event(event);
            if (met>120. && met<170.) _counters.at("SR3_WZ_0Jb").add_event(event);
            if (met>170.) _counters.at("SR3_WZ_0Jc").add_event(event);
          }
          if (mSFOS>81.2 && mSFOS<101.2 && nSignalJets>0) {
            if (met>120. && met<200. && mTmin>110. && pTlll<120. && pT_j1>70.) _counters.at("SR3_WZ_1Ja").add_event(event);
            if (met>200. && mTmin>110. && mTmin<160.) _counters.at("SR3_WZ_1Jb").add_event(event);
            if (met>200. && pT_l2>35. && mTmin>160.) _counters.at("SR3_WZ_1Jc").add_event(event);
          }
        }

        // if (analysis_name().find("200_100") != string::npos) {

 //       cutFlowVector1_str[0] = "All events";
 //          cutFlowVector1_str[1] = "$\\geq$ 2 signal leptons \\& SFOS";
 //          cutFlowVector1_str[2] = "3 signal leptons \\& extra lepton veto";
 //          cutFlowVector1_str[3] = "B-jet veto";
 //          cutFlowVector1_str[4] = "$p_{T}^{l0} > 25 GeV$";
 //          cutFlowVector1_str[5] = "$p_{T}^{l2} > 20 GeV$";
 //          cutFlowVector1_str[6] = "$m_{lll} > 20 GeV$";
 //          cutFlowVector1_str[7] = "$| m_{ll} - m_{Z} | < 10 GeV$";
 //          cutFlowVector1_str[8] = "0 jets";
 //          cutFlowVector1_str[9] = "$60 < E^{miss}_{T} < 120 GeV$";
 //          cutFlowVector1_str[10] = "$m_{T}^{min} > 110 GeV$";
 //          cutFlowVector1_str[11] = "$120 < E^{miss}_{T} < 170 GeV$";
 //          cutFlowVector1_str[12] = "$m_{T}^{min} > 110 GeV$";
 //          cutFlowVector1_str[13] = "$E^{miss}_{T} > 170 GeV$";
 //          cutFlowVector1_str[14] = "$m_{T}^{min} > 110 GeV$";
 //          cutFlowVector1_str[15] = "$\\geq$ 1 jet";
 //          cutFlowVector1_str[16] = "$120 < E^{miss}_{T} < 200 GeV$";
 //          cutFlowVector1_str[17] = "$m_{T}^{min} > 110 GeV$";
 //          cutFlowVector1_str[18] = "$p_{T}^{lll} < 120 GeV$";
 //          cutFlowVector1_str[19] = "$p_{T}^{j0} > 70 GeV$";
 //          cutFlowVector1_str[20] = "$E^{miss}_{T} > 200 GeV$";
 //          cutFlowVector1_str[21] = "$110 < m_{T}^{min} < 160 GeV$";

 //          cutFlowVector1ATLAS_200_100[0]=17682.;
        //   cutFlowVector1ATLAS_200_100[1]=425.63;
 //       cutFlowVector1ATLAS_200_100[2]=424.59;
        //   cutFlowVector1ATLAS_200_100[3]=414.43;
        //   cutFlowVector1ATLAS_200_100[4]=413.98;
        //   cutFlowVector1ATLAS_200_100[5]=306.91;
        //   cutFlowVector1ATLAS_200_100[6]=301.70;
        //   cutFlowVector1ATLAS_200_100[7]=227.15;
        //   cutFlowVector1ATLAS_200_100[8]=110.35;
        //   cutFlowVector1ATLAS_200_100[9]=43.24;
        //   cutFlowVector1ATLAS_200_100[10]=8.91;
        //   cutFlowVector1ATLAS_200_100[11]=6.02;
        //   cutFlowVector1ATLAS_200_100[12]=1.1;
        //   cutFlowVector1ATLAS_200_100[13]=3.15;
        //   cutFlowVector1ATLAS_200_100[14]=0.49;
        //   cutFlowVector1ATLAS_200_100[15]=116.81;
        //   cutFlowVector1ATLAS_200_100[16]=18.86;
        //   cutFlowVector1ATLAS_200_100[17]=5.8;
        //   cutFlowVector1ATLAS_200_100[18]=4.63;
        //   cutFlowVector1ATLAS_200_100[19]=3.18;
        //   cutFlowVector1ATLAS_200_100[20]=7.32;
        //   cutFlowVector1ATLAS_200_100[21]=1.85;

   //        for (size_t j=0;j<NCUTS1;j++){
   //          if(
   //            (j==0) ||

          //     (j==1 && nSignalLeptons>=2 && SFOSpairs.size()>0) ||

          //     (j==2 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3) ||

          //     (j==3 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto) ||

          //     (j==4 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25.) ||

          //     (j==5 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20.) ||

          //     (j==6 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20.) ||

          //     (j==7 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10.) ||

          //     (j==8 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets==0) ||

          //     (j==9 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && met>60. && met<120. && nSignalJets==0) ||

          //     (j==10 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && met>60. && met<120.&& nSignalJets==0 && mTmin>110.) ||

          //     (j==11 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets==0 && met>120. && met<170.) ||

          //     (j==12 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets==0 && met>120. && met<170. && mTmin>110.) ||

          //     (j==13 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets==0 && met>170.) ||

          //     (j==14 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets==0 && met>170. && mTmin>110.) ||

          //     (j==15 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets>0) ||

          //     (j==16 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets>0 && met>120. && met<200.) ||

          //     (j==17 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets>0 && met>120. && met<200. && mTmin>110.) ||

          //     (j==18 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets>0 && met>120. && met<200. && mTmin>110. && pTlll<120.) ||

          //     (j==19 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets>0 && met>120. && met<200. && mTmin>110. && pTlll<120. && pT_j0>70.) ||

          //     (j==20 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets>0 && met>200.) ||

          //     (j==21 && nSignalLeptons>=2 && SFOSpairs.size()>0 && nSignalLeptons==3 && nBaselineLeptons==3 && bjet_veto && pT_l0>25. && pT_l2>20. && mlll>20. && fabs(mSFOS-91.2)<10. && nSignalJets>0 && met>200. && mTmin>110. && mTmin<160.) )

          //     cutFlowVector1[j]++;
          // }

          // cutFlowVector3_str[0] = "All events";
          // cutFlowVector3_str[1] = "2 signal leptons \\& SFOS";
          // cutFlowVector3_str[2] = "B-jet veto";
          // cutFlowVector3_str[3] = "$E_{T}^{miss} > 100 GeV$";
          // cutFlowVector3_str[4] = "2 signal jets";
          // cutFlowVector3_str[5] = "$p_{T}^{j0}, p_{T}^{j1} > 30 GeV$";
          // cutFlowVector3_str[6] = "$81 < m_{Z} < 101 GeV$";
          // cutFlowVector3_str[7] = "$70 < m_{W} < 90 GeV$";
          // cutFlowVector3_str[8] = "$p_{T}^{Z} > 60 GeV$";
          // cutFlowVector3_str[9] = "$\\Delta\\phi(E_{T}^{miss},Z) < 0.8$";
          // cutFlowVector3_str[10] = "$\\Delta\\phi(E_{T}^{miss},W) > 1.5$";
          // cutFlowVector3_str[11] = "$E_{T}^{miss}/p_{T}^{W} < 0.8$";
          // cutFlowVector3_str[12] = "$0.6 < E_{T}^{miss}/p_{T}^{Z} < 1.6$";
          // cutFlowVector3_str[13] = "3-5 signal jets";
          // cutFlowVector3_str[14] = "$p_{T}^{j0}, p_{T}^{j1}, p_{T}^{j2} > 30 GeV$";
          // cutFlowVector3_str[15] = "$81 < m_{Z} < 101 GeV$";
          // cutFlowVector3_str[16] = "$70 < m_{W} < 90 GeV$";
          // cutFlowVector3_str[17] = "$||\\eta (Z)|| < 1.6$";
          // cutFlowVector3_str[18] = "$p_{T}^{Z} > 40 GeV$";
          // cutFlowVector3_str[19] = "$\\Delta\\phi (E_{T}^{miss},ISR) > 2.4$";
          // cutFlowVector3_str[20] = "$\\Delta\\phi (E_{T}^{miss},j1) > 2.6$";
          // cutFlowVector3_str[21] = "$\\Delta\\phi (E_{T}^{miss},W) < 2.2$ ";
          // cutFlowVector3_str[22] = "$0.4 < E_{T}^{miss}/ISR < 0.8$";
          // cutFlowVector3_str[23] = "$\\Delta R(W\\rightarrow 2j) < 2.2$";

     //      cutFlowVector3ATLAS_200_100[0]=20000.;
     //      cutFlowVector3ATLAS_200_100[1]=957.;
     //      cutFlowVector3ATLAS_200_100[2]=880.6;
     //      cutFlowVector3ATLAS_200_100[3]=120.8;
     //      cutFlowVector3ATLAS_200_100[4]=30.2;
     //      cutFlowVector3ATLAS_200_100[5]=20.6;
     //      cutFlowVector3ATLAS_200_100[6]=18.8;
     //      cutFlowVector3ATLAS_200_100[7]=6.2;
     //      cutFlowVector3ATLAS_200_100[8]=5.1;
     //      cutFlowVector3ATLAS_200_100[9]=2.7;
     //      cutFlowVector3ATLAS_200_100[10]=2.7;
     //      cutFlowVector3ATLAS_200_100[11]=2.6;
     //      cutFlowVector3ATLAS_200_100[12]=2.2;
     //      cutFlowVector3ATLAS_200_100[13]=71.7;
     //      cutFlowVector3ATLAS_200_100[14]=47.9;
     //      cutFlowVector3ATLAS_200_100[15]=37.1;
     //      cutFlowVector3ATLAS_200_100[16]=9.3;
     //      cutFlowVector3ATLAS_200_100[17]=7.1;
     //      cutFlowVector3ATLAS_200_100[18]=6.9;
     //      cutFlowVector3ATLAS_200_100[19]=6.3;
     //      cutFlowVector3ATLAS_200_100[20]=5.3;
     //      cutFlowVector3ATLAS_200_100[21]=4.8;
     //      cutFlowVector3ATLAS_200_100[22]=4.0;
     //      cutFlowVector3ATLAS_200_100[23]=3.6;

 //          for (size_t j=0;j<NCUTS3;j++){
 //            if(
 //              (j==0) ||

        //       (j==1 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0) ||

        //       (j==2 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto) ||

        //       (j==3 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100.) ||

        //       (j==4 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets==2) ||

        //       (j==5 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets==2 && pT_j0>30. && pT_j1>30.) ||

        //       (j==6 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets==2 && pT_j0>30. && pT_j1>30. && mll>81. && mll<101.) ||

        //       (j==7 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets==2 && pT_j0>30. && pT_j1>30. && mll>81. && mll<101. && mjj>70. && mjj<90.) ||

        //       (j==8 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets==2 && pT_j0>30. && pT_j1>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.pT()>60.) ||

        //       (j==9 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets==2 && pT_j0>30. && pT_j1>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.pT()>60. && deltaPhi_met_Z<0.8) ||

        //       (j==10 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets==2 && pT_j0>30. && pT_j1>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.pT()>60. && deltaPhi_met_Z<0.8 && deltaPhi_met_W>1.5) ||

        //       (j==11 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets==2 && pT_j0>30. && pT_j1>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.pT()>60. && deltaPhi_met_Z<0.8 && deltaPhi_met_W>1.5 && met/W.pT()<0.8) ||

        //       (j==12 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets==2 && pT_j0>30. && pT_j1>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.pT()>60. && deltaPhi_met_Z<0.8 && deltaPhi_met_W>1.5 && met/W.pT()<0.8 && met/Z.pT()>0.6 && met/Z.pT()<1.6) ||

        //       (j==13 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6) ||

        //       (j==14 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6 && pT_j0>30. && pT_j1>30. && pT_j2>30.) ||

        //       (j==15 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6 && pT_j0>30. && pT_j1>30. && pT_j2>30. && mll>81. && mll<101.) ||

        //       (j==16 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6 && pT_j0>30. && pT_j1>30. && pT_j2>30. && mll>81. && mll<101. && mjj>70. && mjj<90.) ||

        //       (j==17 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6 && pT_j0>30. && pT_j1>30. && pT_j2>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.abseta()<1.6) ||

        //       (j==18 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6 && pT_j0>30. && pT_j1>30. && pT_j2>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.abseta()<1.6 && Z.pT()>40.) ||

        //       (j==19 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6 && pT_j0>30. && pT_j1>30. && pT_j2>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.abseta()<1.6 && Z.pT()>40. && deltaPhi_met_ISR>2.4) ||

        //       (j==20 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6 && pT_j0>30. && pT_j1>30. && pT_j2>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.abseta()<1.6 && Z.pT()>40. && deltaPhi_met_ISR>2.4 && deltaPhi_met_jet0>2.6) ||

        //       (j==21 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6 && pT_j0>30. && pT_j1>30. && pT_j2>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.abseta()<1.6 && Z.pT()>40. && deltaPhi_met_ISR>2.4 && deltaPhi_met_jet0>2.6 && deltaPhi_met_W<2.2) ||

        //       (j==22 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6 && pT_j0>30. && pT_j1>30. && pT_j2>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.abseta()<1.6 && Z.pT()>40. && deltaPhi_met_ISR>2.4 && deltaPhi_met_jet0>2.6 && deltaPhi_met_W<2.2 && met/W_ISR.at(1).pT()>0.4 && met/W_ISR.at(1).pT()<0.8) ||

        //       (j==23 && preselection && pT_l1>25. && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && met>100. && nSignalJets>2 && nSignalJets<6 && pT_j0>30. && pT_j1>30. && pT_j2>30. && mll>81. && mll<101. && mjj>70. && mjj<90. && Z.abseta()<1.6 && Z.pT()>40. && deltaPhi_met_ISR>2.4 && deltaPhi_met_jet0>2.6 && deltaPhi_met_W<2.2 && met/W_ISR.at(1).pT()>0.4 && met/W_ISR.at(1).pT()<0.8 && deltaR_jj<2.2) )

        //       cutFlowVector3[j]++;
        //   }

        // }

        // if ((analysis_name().find("400_200") != string::npos) || (analysis_name().find("500_100") != string::npos)){

 //       cutFlowVector2_str[0] = "All events";
 //       cutFlowVector2_str[1] = "2 signal leptons \\& SFOS";
 //       cutFlowVector2_str[2] = "B-jet veto";
 //       cutFlowVector2_str[3] = "$\\geq$ 2 signal jets";
 //       cutFlowVector2_str[4] = "$p_{T}^{j0}, p_{T}^{j1} > 30 GeV$";
 //       cutFlowVector2_str[5] = "$E_{T}^{miss} > 150 GeV$";
 //       cutFlowVector2_str[6] = "$p_{T}^{Z} > 80 GeV$";
 //       cutFlowVector2_str[7] = "$p_{T}^{W} > 100 GeV$";
 //       cutFlowVector2_str[8] = "$ 81 < m_{Z} < 101 GeV$";
 //       cutFlowVector2_str[9] = "$70 < m_{W} < 100 GeV$";
 //       cutFlowVector2_str[10] = "$m_{T2} > 100 GeV$";
 //       cutFlowVector2_str[11] = "$0.5 < \\Delta\\phi(E_{T}^{miss}, W) < 3.0$";
 //       cutFlowVector2_str[12] = "$\\Delta R(W\\rightarrow jj) <1.5$";
 //       cutFlowVector2_str[13] = "$\\Delta R(Z\\rightarrow ll) <1.8$";

 //          cutFlowVector2ATLAS_400_200[0]=10000.;
 //          cutFlowVector2ATLAS_400_200[1]=83.1;
 //          cutFlowVector2ATLAS_400_200[2]=75.8;
 //          cutFlowVector2ATLAS_400_200[3]=64.7;
 //          cutFlowVector2ATLAS_400_200[4]=53.3;
 //          cutFlowVector2ATLAS_400_200[5]=29.8;
 //          cutFlowVector2ATLAS_400_200[6]=25.0;
 //          cutFlowVector2ATLAS_400_200[7]=20.3;
 //          cutFlowVector2ATLAS_400_200[8]=18.4;
 //          cutFlowVector2ATLAS_400_200[9]=7.7;
 //          cutFlowVector2ATLAS_400_200[10]=5.8;
 //          cutFlowVector2ATLAS_400_200[11]=5.5;
 //          cutFlowVector2ATLAS_400_200[12]=5.4;
 //          cutFlowVector2ATLAS_400_200[13]=5.2;

 //          cutFlowVector2ATLAS_500_100[0]=5000.;
 //          cutFlowVector2ATLAS_500_100[1]=37.9;
 //          cutFlowVector2ATLAS_500_100[2]=33.7;
 //          cutFlowVector2ATLAS_500_100[3]=28.9;
 //          cutFlowVector2ATLAS_500_100[4]=25.3;
 //          cutFlowVector2ATLAS_500_100[5]=20.5;
 //          cutFlowVector2ATLAS_500_100[6]=19.4;
 //          cutFlowVector2ATLAS_500_100[7]=17.5;
 //          cutFlowVector2ATLAS_500_100[8]=15.6;
 //          cutFlowVector2ATLAS_500_100[9]=7.4;
 //          cutFlowVector2ATLAS_500_100[10]=6.7;
 //          cutFlowVector2ATLAS_500_100[11]=5.9;
 //          cutFlowVector2ATLAS_500_100[12]=5.9;
 //          cutFlowVector2ATLAS_500_100[13]=5.9;
        // //maybe add pt_l1>25.
 //          for (size_t j=0;j<NCUTS2;j++){
 //            if(
 //              (j==0) ||

        //       (j==1 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0) ||

        //       (j==2 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto) ||

        //       (j==3 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2) ||

        //       (j==4 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2 && pT_j0>30. && pT_j1>30.) ||

        //       (j==5 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2 && pT_j0>30. && pT_j1>30. && met>150.) ||

        //       (j==6 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2 && pT_j0>30. && pT_j1>30. && met>150. && Z.pT()>80.) ||

        //       (j==7 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2 && pT_j0>30. && pT_j1>30. && met>150. && Z.pT()>80. && W.pT()>100.) ||

        //       (j==8 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2 && pT_j0>30. && pT_j1>30. && met>150. && Z.pT()>80. && W.pT()>100. && mll>81. && mll<101.) ||

        //       (j==9 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2 && pT_j0>30. && pT_j1>30. && met>150. && Z.pT()>80. && W.pT()>100. && mll>81. && mll<101. && mjj>70. && mjj<100.) ||

        //       (j==10 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2 && pT_j0>30. && pT_j1>30. && met>150. && Z.pT()>80. && W.pT()>100. && mll>81. && mll<101. && mjj>70. && mjj<100. && mT2>100.) ||

        //       (j==11 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2 && pT_j0>30. && pT_j1>30. && met>150. && Z.pT()>80. && W.pT()>100. && mll>81. && mll<101. && mjj>70. && mjj<100. && mT2>100. && deltaPhi_met_W>0.5 && deltaPhi_met_W<3.0) ||

        //       (j==12 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2 && pT_j0>30. && pT_j1>30. && met>150. && Z.pT()>80. && W.pT()>100. && mll>81. && mll<101. && mjj>70. && mjj<100. && mT2>100. && deltaPhi_met_W>0.5 && deltaPhi_met_W<3.0 && deltaR_jj<1.5) ||

        //       (j==13 && preselection && nSignalLeptons==2 && SFOSpairs.size()>0 && bjet_veto && nSignalJets>=2 && pT_j0>30. && pT_j1>30. && met>150. && Z.pT()>80. && W.pT()>100. && mll>81. && mll<101. && mjj>70. && mjj<100. && mT2>100. && deltaPhi_met_W>0.5 && deltaPhi_met_W<3.0 && deltaR_jj<1.5 && deltaR_ll<1.8) )

        //       cutFlowVector2[j]++;
        //   }
        // }

        // if (analysis_name().find("800_600") != string::npos){

 //       cutFlowVector4_str[0] = "All events";
 //       cutFlowVector4_str[1] = "3 signal leptons \\& SFOS";
 //       cutFlowVector4_str[2] = "Pass event cleaning";
 //       cutFlowVector4_str[3] = "$m_{T}^{min} > 110 GeV$";
 //       cutFlowVector4_str[4] = "$E_{T}^{miss} > 130 GeV$";
 //       cutFlowVector4_str[5] = "$m^{min}_{SFOS} < 81.2 GeV$";
 //       cutFlowVector4_str[6] = "$20 < p_{T}^{l2} < 30 GeV$";
 //       cutFlowVector4_str[7] = "$p_{T}^{l2} > 30 GeV$";
 //       cutFlowVector4_str[8] = "$m^{min}_{SFOS} > 101.2 GeV$";
 //       cutFlowVector4_str[9] = "$20 < p_{T}^{l2} < 50 GeV$";
 //       cutFlowVector4_str[10] = "$50 < p_{T}^{l2} < 80 GeV$";
 //       cutFlowVector4_str[11] = "$p_{T}^{l2} > 80 GeV$";

 //          // cutFlowVector4ATLAS_800_600[0]=9291.;
 //          // cutFlowVector4ATLAS_800_600[1]=25.13;
 //          // cutFlowVector4ATLAS_800_600[2]=23.54;
 //          // cutFlowVector4ATLAS_800_600[3]=14.43;
 //          // cutFlowVector4ATLAS_800_600[4]=10.22;
 //          // cutFlowVector4ATLAS_800_600[5]=2.10;
 //          // cutFlowVector4ATLAS_800_600[6]=0.11;
 //          // cutFlowVector4ATLAS_800_600[7]=1.99;
 //          // cutFlowVector4ATLAS_800_600[8]=6.8;
 //          // cutFlowVector4ATLAS_800_600[9]=2.53;
 //          // cutFlowVector4ATLAS_800_600[10]=3.01;
 //          // cutFlowVector4ATLAS_800_600[11]=1.25;

 //          for (size_t j=0;j<NCUTS4;j++){
 //            if(
 //              (j==0) ||

        //       (j==1 && nSignalLeptons==3 && SFOSpairs.size()>0) ||

        //       (j==2 && nSignalLeptons==3 && SFOSpairs.size()>0 && preselection && bjet_veto) ||

        //       (j==3 && nSignalLeptons==3 && SFOSpairs.size()>0 && preselection && bjet_veto && mTmin>110.) ||

        //       (j==4 && nSignalLeptons==3 && SFOSpairs.size()>0 && preselection && bjet_veto && mTmin>110. && met>130.) ||

        //       (j==5 && nSignalLeptons==3 && SFOSpairs.size()>0 && preselection && bjet_veto && mTmin>110. && met>130. && mSFOS<81.2) ||

        //       (j==6 && nSignalLeptons==3 && SFOSpairs.size()>0 && preselection && bjet_veto && mTmin>110. && met>130. && mSFOS<81.2 && pT_l2>20. && pT_l2<30.) ||

        //       (j==7 && nSignalLeptons==3 && SFOSpairs.size()>0 && preselection && bjet_veto && mTmin>110. && met>130. && mSFOS<81.2 && pT_l2>30.) ||

        //       (j==8 && nSignalLeptons==3 && SFOSpairs.size()>0 && preselection && bjet_veto && mTmin>110. && met>130. && mSFOS>101.2) ||

        //       (j==9 && nSignalLeptons==3 && SFOSpairs.size()>0 && preselection && bjet_veto && mTmin>110. && met>130. && mSFOS>101.2 && pT_l2>20. && pT_l2<50.) ||

        //       (j==10 && nSignalLeptons==3 && SFOSpairs.size()>0 && preselection && bjet_veto && mTmin>110. && met>130. && mSFOS>101.2 && pT_l2>50. && pT_l2<80.) ||

        //       (j==11 && nSignalLeptons==3 && SFOSpairs.size()>0 && preselection && bjet_veto && mTmin>110. && met>130. && mSFOS>101.2 && pT_l2>80.) )

        //       cutFlowVector4[j]++;
        //   }
        // }

        // if ((analysis_name().find("401_1") != string::npos) || (analysis_name().find("300_150") != string::npos)){

 //       cutFlowVector5_str[0] = "All events";
 //       cutFlowVector5_str[1] = "2 signal leptons \\& OS";
 //       cutFlowVector5_str[2] = "$p_{T}^{l0} > 25 GeV$";
 //       cutFlowVector5_str[3] = "Jet veto";
 //       cutFlowVector5_str[4] = "$m_{ll} > 40 GeV$";
 //       cutFlowVector5_str[5] = "Same flavour";
 //       cutFlowVector5_str[6] = "$m_{ll} > 111 GeV$";
 //       cutFlowVector5_str[7] = "$m_{T2} > 100 GeV$";
 //       cutFlowVector5_str[8] = "Different flavour";
 //       cutFlowVector5_str[9] = "$m_{ll} > 111 GeV$";
 //       cutFlowVector5_str[10] = "$m_{T2} > 100 GeV$";

 //          cutFlowVector5ATLAS_401_1[0]=10000.;
 //          cutFlowVector5ATLAS_401_1[1]=89.7;
 //          cutFlowVector5ATLAS_401_1[2]=89.7;
 //          cutFlowVector5ATLAS_401_1[3]=89.5;
 //          cutFlowVector5ATLAS_401_1[4]=55.7;
 //          cutFlowVector5ATLAS_401_1[5]=55.7;
 //          cutFlowVector5ATLAS_401_1[6]=53.7;
 //          cutFlowVector5ATLAS_401_1[7]=40.4;

 //          cutFlowVector5ATLAS_300_150[0]=25000;
 //          cutFlowVector5ATLAS_300_150[1]=1797.;
 //          cutFlowVector5ATLAS_300_150[2]=1795.3;
 //          cutFlowVector5ATLAS_300_150[3]=1692.1;
 //          cutFlowVector5ATLAS_300_150[4]=1262.;
 //          cutFlowVector5ATLAS_300_150[5]=667.4;
 //          cutFlowVector5ATLAS_300_150[6]=405.;
 //          cutFlowVector5ATLAS_300_150[7]=46.9;
 //          cutFlowVector5ATLAS_300_150[8]=594.5;
 //          cutFlowVector5ATLAS_300_150[9]=363.8;
 //          cutFlowVector5ATLAS_300_150[10]=45.7;

 //          for (size_t j=0;j<NCUTS5;j++){
 //            if(
 //              (j==0) ||

        //       (j==1 && nBaselineLeptons==nSignalLeptons && nSignalLeptons==2 && OSpairs.size()>0) ||

        //       (j==2 && nBaselineLeptons==nSignalLeptons && nSignalLeptons==2 && OSpairs.size()>0 && pT_l0>25.) ||

        //       (j==3 && nBaselineLeptons==nSignalLeptons && nSignalLeptons==2 && OSpairs.size()>0 && pT_l0>25. && central_jet_veto) ||

        //       (j==4 && nBaselineLeptons==nSignalLeptons && nSignalLeptons==2 && OSpairs.size()>0 && pT_l0>25. && central_jet_veto && mll>40.) ||

        //       (j==5 && nBaselineLeptons==nSignalLeptons && nSignalLeptons==2 && OSpairs.size()>0 && pT_l0>25. && central_jet_veto && mll>40. && SFOSpairs.size()>0) ||

        //       (j==6 && nBaselineLeptons==nSignalLeptons && nSignalLeptons==2 && OSpairs.size()>0 && pT_l0>25. && central_jet_veto && mll>40. && SFOSpairs.size()>0 && mll>111.) ||

        //       (j==7 && nBaselineLeptons==nSignalLeptons && nSignalLeptons==2 && OSpairs.size()>0 && pT_l0>25. && central_jet_veto && mll>40. && SFOSpairs.size()>0 && mll>111. && mT2>100.) ||

        //       (j==8 && nBaselineLeptons==nSignalLeptons && nSignalLeptons==2 && OSpairs.size()>0 && pT_l0>25. && central_jet_veto && mll>40. && SFOSpairs.size()==0) ||

        //       (j==9 && nBaselineLeptons==nSignalLeptons && nSignalLeptons==2 && OSpairs.size()>0 && pT_l0>25. && central_jet_veto && mll>40. && SFOSpairs.size()==0 && mll>111.) ||

        //       (j==10 && nBaselineLeptons==nSignalLeptons && nSignalLeptons==2 && OSpairs.size()>0 && pT_l0>25. && central_jet_veto && mll>40. && SFOSpairs.size()==0 && mll>111. && mT2>100.) )

        //       cutFlowVector5[j]++;
        //   }
        // }

      }

      /// Combine the variables of another copy of this analysis (typically on another thread) into this one.
      void combine(const Analysis* other)
      {
        const Analysis_ATLAS_13TeV_MultiLEP_confnote_36invfb* specificOther
                = dynamic_cast<const Analysis_ATLAS_13TeV_MultiLEP_confnote_36invfb*>(other);

        for (auto& pair : _counters) { pair.second += specificOther->_counters.at(pair.first); }

        if (NCUTS1 != specificOther->NCUTS1) NCUTS1 = specificOther->NCUTS1;
        for (size_t j = 0; j < NCUTS1; j++)
        {
          cutFlowVector1[j] += specificOther->cutFlowVector1[j];
          cutFlowVector1_str[j] = specificOther->cutFlowVector1_str[j];
        }

        if (NCUTS2 != specificOther->NCUTS2) NCUTS2 = specificOther->NCUTS2;
        for (size_t j = 0; j < NCUTS2; j++)
        {
          cutFlowVector2[j] += specificOther->cutFlowVector2[j];
          cutFlowVector2_str[j] = specificOther->cutFlowVector2_str[j];
        }

      }


      void collect_results() {

        add_result(SignalRegionData(_counters.at("SR2_SF_loose"), 153., { 133., 22.}));
        add_result(SignalRegionData(_counters.at("SR2_SF_tight"), 9., { 9.8, 2.9}));
        add_result(SignalRegionData(_counters.at("SR2_DF_100"), 78., { 68., 7.}));
        add_result(SignalRegionData(_counters.at("SR2_DF_150"), 11,  { 11.5, 3.1}));
        add_result(SignalRegionData(_counters.at("SR2_DF_200"), 6., { 2.1, 1.9}));
        add_result(SignalRegionData(_counters.at("SR2_DF_300"), 2., { 0.6, 0.6}));
        add_result(SignalRegionData(_counters.at("SR2_int"), 2., { 4.1, 2.6}));
        add_result(SignalRegionData(_counters.at("SR2_high"), 0., { 1.6, 1.6}));
        add_result(SignalRegionData(_counters.at("SR2_low"), 11., { 4.2, 3.8}));
        add_result(SignalRegionData(_counters.at("SR3_slep_a"), 4., { 2.23, 0.79}));
        add_result(SignalRegionData(_counters.at("SR3_slep_b"), 3., { 2.79, 0.43}));
        add_result(SignalRegionData(_counters.at("SR3_slep_c"), 9., { 5.41, 0.93}));
        add_result(SignalRegionData(_counters.at("SR3_slep_d"), 0., { 1.42, 0.38}));
        add_result(SignalRegionData(_counters.at("SR3_slep_e"), 0., { 1.14, 0.23}));
        add_result(SignalRegionData(_counters.at("SR3_WZ_0Ja"), 21., { 21.74, 2.85}));
        add_result(SignalRegionData(_counters.at("SR3_WZ_0Jb"), 1., { 2.68, 0.46}));
        add_result(SignalRegionData(_counters.at("SR3_WZ_0Jc"), 2., { 1.56, 0.33}));
        add_result(SignalRegionData(_counters.at("SR3_WZ_1Ja"), 1., { 2.21, 0.53}));
        add_result(SignalRegionData(_counters.at("SR3_WZ_1Jb"), 3., { 1.82, 0.26}));
        add_result(SignalRegionData(_counters.at("SR3_WZ_1Jc"), 4., { 1.26, 0.34}));

      }


      vector<HEPUtils::P4> get_W_ISR(vector<const HEPUtils::Jet*> jets, HEPUtils::P4 Z, HEPUtils::P4 met) {
        HEPUtils::P4 Z_met_sys=Z+met;
        double deltaR_min=999;
        size_t Wjets_id1;
        size_t Wjets_id2;
        for (size_t i=0;i<jets.size();i++) {
          for (size_t j=0;j<jets.size();j++) {
            if (i!=j) {
              HEPUtils::P4 jj_sys=jets.at(i)->mom()+jets.at(j)->mom();
              double deltaR=fabs(jj_sys.deltaR_eta(Z_met_sys));
              if (deltaR<deltaR_min) {
                deltaR_min=deltaR;
                Wjets_id1=i;
                Wjets_id2=j;
              }
            }
          }
        }
        HEPUtils::P4 W=jets.at(Wjets_id2)->mom()+jets.at(Wjets_id1)->mom();
        HEPUtils::P4 ISR;
        HEPUtils::P4 j0=jets.at(Wjets_id2)->mom();
        HEPUtils::P4 j1=jets.at(Wjets_id1)->mom();
        for (size_t k=0;k<jets.size();k++) {
          if ((k!=Wjets_id1) && (k!=Wjets_id2))ISR+=(jets.at(k)->mom());
        }
        vector<HEPUtils::P4> W_ISR;
        W_ISR.push_back(W);
        W_ISR.push_back(ISR);
        W_ISR.push_back(j0);
        W_ISR.push_back(j1);
        return W_ISR;
      }


    protected:
      void analysis_specific_reset() {

        for (auto& pair : _counters) { pair.second.reset(); }

        std::fill(cutFlowVector1.begin(), cutFlowVector1.end(), 0);
        std::fill(cutFlowVector2.begin(), cutFlowVector2.end(), 0);
        std::fill(cutFlowVector3.begin(), cutFlowVector3.end(), 0);
        std::fill(cutFlowVector4.begin(), cutFlowVector4.end(), 0);
        std::fill(cutFlowVector5.begin(), cutFlowVector5.end(), 0);
      }

    };


    // Factory fn
    DEFINE_ANALYSIS_FACTORY(ATLAS_13TeV_MultiLEP_confnote_36invfb)


  }
}

Updated on 2024-07-18 at 13:53:34 +0000