add bbyy normalisation and adapt to splitting b and tau-tags in k4sim…

…delphes

examples/FCChh/ggHH_bbyy/analysis_final.py

@@ -5,14 +5,23 @@
 outputDir  = "outputs/FCChh/ggHH_bbyy/final/"
 
 processList = {
-    'pwp8_pp_hh_5f_hhbbyy':{},#Run over the full statistics from stage2 input file <inputDir>/p8_ee_ZZ_ecm240.root. Keep the same output name as input
+    # 'pwp8_pp_hh_5f_hhbbyy':{},#Run over the full statistics from stage2 input file <inputDir>/p8_ee_ZZ_ecm240.root. Keep the same output name as input
+    'pwp8_pp_hh_5f_hhbbyy_split_HF_tau_tags':{},#Run over the full statistics from stage2 input file <inputDir>/p8_ee_ZZ_ecm240.root. Keep the same output name as input
 }
 
 #Link to the dictonary that contains all the cross section informations etc...
 procDict = "FCCee_procDict_spring2021_IDEA.json" # will need an FCC-hh one!
 
 #Add MySample_p8_ee_ZH_ecm240 as it is not an offical process TO UPDATE 
-procDictAdd={"pwp8_pp_hh_5f_hhbbyy":{"numberOfEvents": 10000000, "sumOfWeights": 10000000, "crossSection": 1.0, "kfactor": 1.0, "matchingEfficiency": 1.0}}
+# procDictAdd={"pwp8_pp_hh_5f_hhbbyy":{"numberOfEvents": 10000000, "sumOfWeights": 10000000, "crossSection": 1.0, "kfactor": 1.0, "matchingEfficiency": 1.0}}
+# procDictAdd={"pwp8_pp_hh_5f_hhbbyy_split_HF_tau_tags":{"numberOfEvents": 10000000, "sumOfWeights": 10000000, "crossSection": 1.0, "kfactor": 1.0, "matchingEfficiency": 1.0}}
+procDictAdd={"pwp8_pp_hh_5f_hhbbyy_split_HF_tau_tags": {"numberOfEvents": 4980000, "sumOfWeights": 4980000.0, "crossSection": 0.0029844128399999998, "kfactor": 1.075363, "matchingEfficiency": 1.0}}
+
+# Expected integrated luminosity
+intLumi = 30e+06  # pb-1
+
+# Whether to scale to expected integrated luminosity
+doScale = True
 
 #Number of CPUs to use
 nCPUS = 2

examples/FCChh/ggHH_bbyy/analysis_stage1.py

@@ -28,7 +28,8 @@ def __init__(self, cmdline_args):
             # # - <name of process> needs to correspond either the name of the input .root file, or the name of a directory containing root files 
             # # If you want to process only part of the events, split the output into chunks or give a different name to the output use the optional arguments
             # # or leave blank to use defaults = run the full statistics in one output file named the same as the process:
-            'pwp8_pp_hh_5f_hhbbyy': {},
+            'pwp8_pp_hh_5f_hhbbyy_split_HF_tau_tags': {},
+            # 'pwp8_pp_hh_5f_hhbbyy': {},
         }
 
         # Mandatory: Input directory where to find the samples, or a production tag when running over the centrally produced
@@ -112,7 +113,7 @@ def analyzers(self, dframe):
 
             #b-tagged jets:
             #b tagged jets
-            .Define("bjets", "AnalysisFCChh::get_tagged_jets(Jet, Jet_tags, _Jet_tags_particle, _Jet_tags_parameters, 1)") #bit 1 = medium WP, see: https://github.com/delphes/delphes/blob/master/cards/FCC/scenarios/FCChh_I.tcl
+            .Define("bjets", "AnalysisFCChh::get_tagged_jets(Jet, Jet_HF_tags, _Jet_HF_tags_particle, _Jet_HF_tags_parameters, 1)") #bit 1 = medium WP, see: https://github.com/delphes/delphes/blob/master/cards/FCC/scenarios/FCChh_I.tcl
             .Define("selpt_bjets", "FCCAnalyses::ReconstructedParticle::sel_pt(30.)(bjets)")
             .Define("sel_bjets_unsort", "FCCAnalyses::ReconstructedParticle::sel_eta(4)(selpt_bjets)")
             .Define("sel_bjets", "AnalysisFCChh::SortParticleCollection(sel_bjets_unsort)")
@@ -134,12 +135,12 @@ def analyzers(self, dframe):
             ########################################### APPLY PRE-SELECTION ########################################### 
 
             #require at least two b-jets and two photons, both with invariant masses compatible with the Higgs mass
-            # .Filter("sel_bjets.size()>1")
-            # .Filter("sel_gamma.size()>1") 
-            # .Filter("m_bb[0] < 200.") 
-            # .Filter("m_bb[0] > 80.") 
-            # .Filter("m_yy[0] < 180.") 
-            # .Filter("m_yy[0] > 100.")   
+            .Filter("sel_bjets.size()>1")
+            .Filter("sel_gamma.size()>1") 
+            .Filter("m_bb[0] < 200.") 
+            .Filter("m_bb[0] > 80.") 
+            .Filter("m_yy[0] < 180.") 
+            .Filter("m_yy[0] > 100.")   
 
         )
         return dframe2