New updates for making run_plot faster, usage of two count histograms for NLO V+jets, update 2018 configs

README.md

@@ -8,3 +8,19 @@ Setup & instructions
 Full instructions on how to setup the code and run the analysis can be found on the Twiki page:
 
 https://twiki.cern.ch/twiki/bin/viewauth/CMS/HiggsBBcodeETH
+
+
+**Branch: HIG-20-001-NLO2017_2018** <br />
+**Used for NLO 2017 and 2018 STXS analysis** <br />
+
+**Final ntuples:**<br />
+**2017**: <br />
+**Zll** : root://t3dcachedb03.psi.ch:1094//pnfs/psi.ch/cms/trivcat/store/user/creissel/VHbb/Zll/VHbbPostNano2017_V11/2022_02_06/eval/STXSunc/ <br />
+**Wlv** : root://t3dcachedb03.psi.ch:1094//pnfs/psi.ch/cms/trivcat/store/user/creissel/VHbb/Wlv/VHbbPostNano2017_V11/2022_02_06/eval/STXSunc/ <br />
+**Zvv** : root://t3dcachedb03.psi.ch:1094//pnfs/psi.ch/cms/trivcat/store/user/creissel/VHbb/Zvv/VHbbPostNano2017_V11/2022_02_06/eval/STXSunc/ <br />
+
+**2018**: <br />
+**Zll** : root://t3dcachedb03.psi.ch:1094//pnfs/psi.ch/cms/trivcat/store/user/krgedia/VHbb/Zll/VHbbPostNano2018/mva/13may22_vjets_vptnlowgtcorrect_  <br />
+**Wlv** : root://t3dcachedb03.psi.ch:1094//pnfs/psi.ch/cms/trivcat/store/user/krgedia/VHbb/VHbbPostNano2018/V12/Wlv/mva/13may22_vjets_vptnlowgtcorrect_ <br />
+**Zvv** : root://t3dcachedb03.psi.ch:1094//pnfs/psi.ch/cms/trivcat/store/user/krgedia/VHbb/VHbbPostNano2018/V12/Zvv/eval/13may22_vjets_vptnlowgtcorrect_ <br />
+

interface/VHbbNameSpace.h

@@ -258,6 +258,12 @@ double SoverSBWeight(double BDT, int channel) {
     return TMath::Sqrt(deta*deta + dphi*dphi);
   }
 
+
+  double deltaR_v2(double deta, double dphi) {
+    return TMath::Sqrt(deta*deta + dphi*dphi);
+  }
+
+
   double HJetVarByCSV(double CSVj1, double Varj1, double CSVj2, double Varj2, bool pickhighestBtag) {
     if ((CSVj1 > CSVj2) && pickhighestBtag) {
       return Varj1;

python/VHbb2017commonconfig/general.ini

@@ -14,6 +14,7 @@ Run_D = (run>=302030&&run<=303434)
 Run_E = (run>=303435&&run<=304826)
 Run_F = (run>=304911&&run<=306462)
 
+VZ_optimize=False 
 Debug = True
 merge = True
 weightexpression = 1
@@ -89,6 +90,30 @@ diJetEta = <!General|Hbtag!>_eta
 Jet_tight = (Jet_puId>6||Jet_Pt>50)&&Jet_jetId>4&&Jet_lepFilter 
 Jet_additional = Iteration$!=<!General|btagidx0!>&&Iteration$!=<!General|btagidx1!>
 
+resolvedBtagWeightApplies  = (<!Cuts|isResolved!>)
+DeepCSVweightApplies = <!.|resolvedBtagWeightApplies!>
+
+weightFxLHEScaleWeight6    = (<!Weights|weightF!>*LHEScaleWeight[6])
+weightFxLHEScaleWeight4    = (<!Weights|weightF!>*LHEScaleWeight[4])
+weightFxLHEScaleWeight3    = (<!Weights|weightF!>*LHEScaleWeight[3])
+weightFxLHEScaleWeight1    = (<!Weights|weightF!>*LHEScaleWeight[1])
+weightnom                  = (<!Weights|weightF!>)
+
+LHEScalefactor2AppliesmuF    = ((( <!General|weightFxLHEScaleWeight3!> < <!General|weightnom!> )&&( <!General|weightFxLHEScaleWeight4!> < <!General|weightnom!> )) || ((<!General|weightFxLHEScaleWeight3!> > <!General|weightnom!>)&&(<!General|weightFxLHEScaleWeight4!> > <!General|weightnom!>)))
+LHEScalefactor2AppliesmuR    = ((( <!General|weightFxLHEScaleWeight1!> < <!General|weightnom!> )&&( <!General|weightFxLHEScaleWeight6!> < <!General|weightnom!> )) || ((<!General|weightFxLHEScaleWeight1!> > <!General|weightnom!>)&&(<!General|weightFxLHEScaleWeight6!> > <!General|weightnom!>)))
+
+LHEScalefactor2AppliesEye    =( ((sampleIndex==13000)||(sampleIndex==13200)||(sampleIndex==13400)||(sampleIndex==13600)||(sampleIndex==13300)||(sampleIndex==13500)||(sampleIndex==13700)||(sampleIndex==18100)||(sampleIndex==18200)||(sampleIndex==6000)||(sampleIndex==6200)||(sampleIndex==6400)||(sampleIndex==6600)||(sampleIndex==6800)||(sampleIndex==7000)||(sampleIndex==6500)||(sampleIndex==6700)||(sampleIndex==6900)||(sampleIndex==7100)||(sampleIndex==19100)||(sampleIndex==19200)||(sampleIndex==17000)||(sampleIndex==17200)||(sampleIndex==17400)||(sampleIndex==17600)||(sampleIndex==17100)||(sampleIndex==17300)||(sampleIndex==17500)||(sampleIndex==17700)) || ((sampleIndex==13001)||(sampleIndex==13201)||(sampleIndex==13401)||(sampleIndex==13601)||(sampleIndex==13301)||(sampleIndex==13501)||(sampleIndex==13701)||(sampleIndex==18101)||(sampleIndex==18201)||(sampleIndex==6001)||(sampleIndex==6201)||(sampleIndex==6401)||(sampleIndex==6601)||(sampleIndex==6801)||(sampleIndex==7001)||(sampleIndex==6501)||(sampleIndex==6701)||(sampleIndex==6901)||(sampleIndex==7101)||(sampleIndex==19101)||(sampleIndex==19201)||(sampleIndex==17001)||(sampleIndex==17201)||(sampleIndex==17401)||(sampleIndex==17601)||(sampleIndex==17101)||(sampleIndex==17301)||(sampleIndex==17501)||(sampleIndex==17701)) ||  ((sampleIndex==13002)||(sampleIndex==13202)||(sampleIndex==13402)||(sampleIndex==13602)||(sampleIndex==13302)||(sampleIndex==13502)||(sampleIndex==13702)||(sampleIndex==18102)||(sampleIndex==18202)||(sampleIndex==6002)||(sampleIndex==6202)||(sampleIndex==6402)||(sampleIndex==6602)||(sampleIndex==6802)||(sampleIndex==7002)||(sampleIndex==6502)||(sampleIndex==6702)||(sampleIndex==6902)||(sampleIndex==7102)||(sampleIndex==19102)||(sampleIndex==19202)||(sampleIndex==17002)||(sampleIndex==17202)||(sampleIndex==17402)||(sampleIndex==17602)||(sampleIndex==17102)||(sampleIndex==17302)||(sampleIndex==17502)||(sampleIndex==17702)) ||  ((sampleIndex==13004)||(sampleIndex==13204)||(sampleIndex==13404)||(sampleIndex==13604)||(sampleIndex==13304)||(sampleIndex==13504)||(sampleIndex==13704)||(sampleIndex==18104)||(sampleIndex==18204)||(sampleIndex==6004)||(sampleIndex==6204)||(sampleIndex==6404)||(sampleIndex==6604)||(sampleIndex==6804)||(sampleIndex==7004)||(sampleIndex==6504)||(sampleIndex==6704)||(sampleIndex==6904)||(sampleIndex==7104)||(sampleIndex==19104)||(sampleIndex==19204)||(sampleIndex==17004)||(sampleIndex==17204)||(sampleIndex==17404)||(sampleIndex==17604)||(sampleIndex==17104)||(sampleIndex==17304)||(sampleIndex==17504)||(sampleIndex==17704)) ) 
+
+
+LHE1ScaleWeightforW2Jets = ( (V_pt>=100.0&&V_pt<150.0)*(1.32) + (V_pt>=150.0&&V_pt<250.0)*(1.3) + (V_pt>=250.0&&V_pt<400.0)*(1.28) + (V_pt>=400.0)*(1.26) )
+
+LHE3ScaleWeightforW2Jets = ( (V_pt>=100.0&&V_pt<150.0)*(1.02) + (V_pt>=150.0&&V_pt<250.0)*(1.04) + (V_pt>=250.0&&V_pt<400.0)*(1.06) + (V_pt>=400.0)*(1.08) )
+
+LHE5ScaleWeightforW2Jets = ( (V_pt>=100.0&&V_pt<150.0)*(0.98) + (V_pt>=150.0&&V_pt<250.0)*(0.96) + (V_pt>=250.0&&V_pt<400.0)*(0.96) + (V_pt>=400.0)*(0.94) )
+
+LHE7ScaleWeightforW2Jets = ( (V_pt>=100.0&&V_pt<150.0)*(0.78) + (V_pt>=150.0&&V_pt<250.0)*(0.8) + (V_pt>=250.0&&V_pt<400.0)*(0.8) + (V_pt>=400.0)*(0.8) )
+
+
 #-------------------------------------------------
 # Scale Factors
 
@@ -461,6 +486,13 @@ sys_cut_include = []
 
 [Cuts]
 
+dRbb_gt0p5 = ((-99+Alt$((hJidx[0]>-1)*(hJidx[1]>-1)*(99+(VHbb::deltaR_v2(hJets_FSRrecovered_dEta,hJets_FSRrecovered_dPhi))),0)) > 0.5)
+dRbb_gt0p6 = ((-99+Alt$((hJidx[0]>-1)*(hJidx[1]>-1)*(99+(VHbb::deltaR_v2(hJets_FSRrecovered_dEta,hJets_FSRrecovered_dPhi))),0)) > 0.6)
+dRbb_gt0p7 = ((-99+Alt$((hJidx[0]>-1)*(hJidx[1]>-1)*(99+(VHbb::deltaR_v2(hJets_FSRrecovered_dEta,hJets_FSRrecovered_dPhi))),0)) > 0.7)
+dRbb_gt0p8 = ((-99+Alt$((hJidx[0]>-1)*(hJidx[1]>-1)*(99+(VHbb::deltaR_v2(hJets_FSRrecovered_dEta,hJets_FSRrecovered_dPhi))),0)) > 0.8)
+dRbb_gt0p9 = ((-99+Alt$((hJidx[0]>-1)*(hJidx[1]>-1)*(99+(VHbb::deltaR_v2(hJets_FSRrecovered_dEta,hJets_FSRrecovered_dPhi))),0)) > 0.9)
+dRbb_gt1p0 = ((-99+Alt$((hJidx[0]>-1)*(hJidx[1]>-1)*(99+(VHbb::deltaR_v2(hJets_FSRrecovered_dEta,hJets_FSRrecovered_dPhi))),0)) > 1.0)
+
 # --------------------------------------------------------------------------------------------------
 # MVA training 
 # --------------------------------------------------------------------------------------------------
@@ -483,24 +515,46 @@ boostedJets  = (Hbb_fjidx>-1)
 # isResolved should be required for all resolved regions and isBoosted for all boosted regions
 #isResolved_BOOSTED  = (<!.|resolvedJets!> && (!<!Cuts|all_BOOST!>))
 
+## flag isBoosted = <!Cuts|all_BOOST!>
+#isResolved_BOOSTED  = (<!.|resolvedJets!> && !(boostedCR||boostedSR))
+#isResolved_RESOLVED = (<!.|resolvedJets!>)
+#isResolved_SR_A     = (<!.|resolvedJets!> && (hJidx[0]>-1&&hJidx[1]>-1) && (!(boostedCR||boostedSR)||(resolvedSR&&boostedCR)))
+#isResolved_SR_B     = (<!.|resolvedJets!> && !(boostedSR&&!resolvedSR))
+#isResolved_SR_B_VZ     = (<!.|resolvedJets!> && !(boostedSR_VZ&&!resolvedSR_VZ))
+#isResolved          = <!.|isResolved_<!.|OverlapResolvedBoosted!>!>
+#isResolved_VZ          = (<!.|resolvedJets!> && !(boostedSR_VZ&&!resolvedSR_VZ))
+#
+#isBoosted_BOOSTED   = (<!.|boostedJets!>)
+#isBoosted_RESOLVED  = (<!.|boostedJets!> && !(resolvedCR||resolvedSR))
+#isBoosted_SR_A      = (<!.|boostedJets!> && !(!(boostedCR||boostedSR)||(resolvedSR&&boostedCR)))
+#isBoosted_SR_B      = (<!.|boostedJets!> && (!(resolvedSR||resolvedCR)||(boostedSR&&resolvedCR)))
+#isBoosted           = <!.|isBoosted_<!.|OverlapResolvedBoosted!>!>
+#isBoosted_VZ           = (<!.|boostedJets!> && (!(resolvedSR_VZ||resolvedCR_VZ)||(boostedSR_VZ&&resolvedCR_VZ)))
+#
+#
+#isBoosted_VZ          = (<!.|boostedJets!> && (!(resolvedSR_VZ||resolvedCR_VZ)||(boostedSR_VZ&&resolvedCR_VZ)))
+
+
 # flag isBoosted = <!Cuts|all_BOOST!>
 isResolved_BOOSTED  = (<!.|resolvedJets!> && !(boostedCR||boostedSR))
 isResolved_RESOLVED = (<!.|resolvedJets!>)
 isResolved_SR_A     = (<!.|resolvedJets!> && (hJidx[0]>-1&&hJidx[1]>-1) && (!(boostedCR||boostedSR)||(resolvedSR&&boostedCR)))
 isResolved_SR_B     = (<!.|resolvedJets!> && !(boostedSR&&!resolvedSR))
-isResolved_SR_B_VZ     = (<!.|resolvedJets!> && !(boostedSR_VZ&&!resolvedSR_VZ))
-isResolved          = <!.|isResolved_<!.|OverlapResolvedBoosted!>!>
-isResolved_VZ          = (<!.|resolvedJets!> && !(boostedSR_VZ&&!resolvedSR_VZ))
+isResolved_VZ_SR_B  = (<!.|resolvedJets!> && !(boostedSR_VZ&&!resolvedSR_VZ))  
+isResolved_False    = <!.|isResolved_<!.|OverlapResolvedBoosted!>!>
+isResolved_True     = <!.|isResolved_VZ_<!.|OverlapResolvedBoosted!>!>
+isResolved          = (<!.|isResolved_<!General|VZ_optimize!>!>)
 
 isBoosted_BOOSTED   = (<!.|boostedJets!>)
 isBoosted_RESOLVED  = (<!.|boostedJets!> && !(resolvedCR||resolvedSR))
 isBoosted_SR_A      = (<!.|boostedJets!> && !(!(boostedCR||boostedSR)||(resolvedSR&&boostedCR)))
 isBoosted_SR_B      = (<!.|boostedJets!> && (!(resolvedSR||resolvedCR)||(boostedSR&&resolvedCR)))
-isBoosted           = <!.|isBoosted_<!.|OverlapResolvedBoosted!>!>
-isBoosted_VZ           = (<!.|boostedJets!> && (!(resolvedSR_VZ||resolvedCR_VZ)||(boostedSR_VZ&&resolvedCR_VZ)))
+isBoosted_VZ_SR_B   = (<!.|boostedJets!> && (!(resolvedSR_VZ||resolvedCR_VZ)||(boostedSR_VZ&&resolvedCR_VZ)))
+isBoosted_False     = <!.|isBoosted_<!.|OverlapResolvedBoosted!>!>
+isBoosted_True      = <!.|isBoosted_VZ_<!.|OverlapResolvedBoosted!>!>
+isBoosted           = (<!.|isBoosted_<!General|VZ_optimize!>!>)
 
 
-isBoosted_VZ          = (<!.|boostedJets!> && (!(resolvedSR_VZ||resolvedCR_VZ)||(boostedSR_VZ&&resolvedCR_VZ)))
 
 # essential branches which are always kept when skimming, even though they are not strictly needed for given cuts/variables
 [Branches]