tb leptonic: redefine label for UE; add option for doDY; fix the mm c…

…hannel
MiT-HEP · Oct 25, 2018 · d5c9373 · d5c9373
1 parent aeab466
commit d5c9373
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 3 deletions.
diff --git a/interface/AnalysisChargedHiggsTopBottom.hpp b/interface/AnalysisChargedHiggsTopBottom.hpp
@@ -22,10 +22,11 @@ class ChargedHiggsTopBottom:  virtual public AnalysisBase
     ChargedHiggsTopBottom() : AnalysisBase () {}
     virtual ~ChargedHiggsTopBottom () {}
 
-    bool doJECSources = true;
+    bool doJECSources = false;
 
     bool doSynch = false;
     bool doQCD = false;
+    bool doDY = false;
     bool doBDTSyst = true;
     bool doFinal = true; // true for final ; no need to compute fancy stuff
     bool writeTree = false; // false for final

diff --git a/src/AnalysisChargedHiggsTopBottom.cpp b/src/AnalysisChargedHiggsTopBottom.cpp
@@ -3310,6 +3310,10 @@ int ChargedHiggsTopBottom::analyze(Event*e,string systname)
     if ( label == "WZZ") label = "EWK";
     if ( label == "ZZZ") label = "EWK";
 
+    // redefine plots
+    if ( label == "TT_TuneCUETP8M2T4up_13TeV-powheg-pythia8" ) label = "TT_TuneCUETP8M2T4_13TeV-powheg-pythia8";
+    if ( label == "TT_TuneCUETP8M2T4down_13TeV-powheg-pythia8" ) label = "TT_TuneCUETP8M2T4_13TeV-powheg-pythia8";
+
     /*
     std::cout << e->GetName() << std::endl;
     if(e->GetName().find("SingleMuon")!=string::npos) std::cout << "this should be SingleMuon" << std::endl;
@@ -3386,9 +3390,10 @@ int ChargedHiggsTopBottom::analyze(Event*e,string systname)
             if(doQCD and do2lAnalysis) pairSign = (it->Charge()*leadLep->Charge()>0);
 
             if( pairSign and (it->GetP4() + leadLep->GetP4()).M() > 12 ) {
-                if(pairSign || (it->IsElectron() && leadLep->IsElectron()) ) {
+                if(((it->IsMuon() && leadLep->IsMuon())  || (it->IsElectron() && leadLep->IsElectron()))) {
                     // same flavour
-                    if( (std::abs((it->GetP4() + leadLep->GetP4()).M()-91) > 15) && (it->GetP4() + leadLep->GetP4()).M()>50 ) { trailLep=it; nOSLepPair++; }
+                    if( not doDY && (std::abs((it->GetP4() + leadLep->GetP4()).M()-91) > 15) && (it->GetP4() + leadLep->GetP4()).M()>50 ) { trailLep=it; nOSLepPair++; }
+                    if( doDY && (std::abs((it->GetP4() + leadLep->GetP4()).M()-91) < 15)) { trailLep=it; nOSLepPair++; }
                 } else { trailLep=it; nOSLepPair++; }
             } // and OS
         }
@@ -3661,6 +3666,7 @@ int ChargedHiggsTopBottom::analyze(Event*e,string systname)
     ////$$$$$$$ Apply SF
     ////$$$$$$$
 
+    /* // do not apply topPT in the central value
     if (not e->IsRealData() and
         ((label.find("TTTo2L2Nu_TuneCUETP8M2")!=string::npos) || (label.find("TTToSemilepton_TuneCUETP8M2")!=string::npos)
          || (label.find("TTTo2L2Nu_ttbbFilter")!=string::npos) || (label.find("TTToSemilepton_ttbbFilter")!=string::npos)
@@ -3671,6 +3677,7 @@ int ChargedHiggsTopBottom::analyze(Event*e,string systname)
          || (label.find("TT_TuneCUETP8M2T4_13TeV-powheg-isrup")!=string::npos)
          || (label.find("TTJets")!=string::npos)
          )) { e->ApplyTopReweight(); }
+    */
 
     ////$$$$$$$
     ////$$$$$$$ Fill Cut Flow