Rely on GeneratorPythia8 for HF particle decays

Author: fgrosa

MC/config/PWGHF/ini/GeneratorHFTrigger_bbbar.ini

@@ -1,24 +1,17 @@
-### The setup uses an external event generator
-### This part sets the path of the file and the function call to retrieve it
-
-[GeneratorExternal]
-fileName = ${O2DPG_ROOT}/MC/config/PWGHF/external/generator/GeneratorHF.C
-funcName = GeneratorHF_bbbar()
-
 ### The external generator derives from GeneratorPythia8.
 ### This part configures the bits of the interface: configuration and user hooks
 
 [GeneratorPythia8]
-config = ${O2DPG_ROOT}/MC/config/common/pythia8/generator/pythia8_hf.cfg
+config = ${O2DPG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmtriggers_with_decays.cfg
 hooksFileName = ${O2DPG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
 hooksFuncName = pythia8_userhooks_bbbar(-1.5,1.5)
 
 ### The setup uses the base configuration of the decayer which is loaded from the file specified by config[0].
 ### On top of the base configuration, two more sets of settings are loaded sequentially from config[1] and [2].
 
-[DecayerPythia8]
-config[0] = ${O2DPG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
-config[1] = ${O2DPG_ROOT}/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmHadTrigger.cfg
+# [DecayerPythia8]
+# config[0] = ${O2DPG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
+# config[1] = ${O2DPG_ROOT}/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmHadTrigger.cfg
 
 ### The setup forces some particles to be decayed by the external decayer instead of Geant.
 ### The PDG list of the particles is specified below.

MC/config/PWGHF/ini/GeneratorHFTrigger_ccbar.ini

@@ -1,24 +1,17 @@
-### The setup uses an external event generator
-### This part sets the path of the file and the function call to retrieve it
-
-[GeneratorExternal]
-fileName = ${O2DPG_ROOT}/MC/config/PWGHF/external/generator/GeneratorHF.C
-funcName = GeneratorHF_ccbar()
-
 ### The external generator derives from GeneratorPythia8.
 ### This part configures the bits of the interface: configuration and user hooks
 
 [GeneratorPythia8]
-config = ${O2DPG_ROOT}/MC/config/common/pythia8/generator/pythia8_hf.cfg
+config = ${O2DPG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmtriggers_with_decays.cfg
 hooksFileName = ${O2DPG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
 hooksFuncName = pythia8_userhooks_ccbar(-1.5,1.5)
 
 ### The setup uses the base configuration of the decayer which is loaded from the file specified by config[0].
 ### On top of the base configuration, two more sets of settings are loaded sequentially from config[1] and [2].
 
-[DecayerPythia8]
-config[0] = ${O2DPG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
-config[1] = ${O2DPG_ROOT}/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmHadTrigger.cfg
+# [DecayerPythia8]
+# config[0] = ${O2DPG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
+# config[1] = ${O2DPG_ROOT}/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmHadTrigger.cfg
 
 ### The setup forces some particles to be decayed by the external decayer instead of Geant.
 ### The PDG list of the particles is specified below.

…cayer/force_hadronic_charmHadTrigger.cfg …or/pythia8_charmtriggers_with_decays.cfgMC/config/PWGHF/pythia8/decayer/force_hadronic_charmHadTrigger.cfg renamed to MC/config/PWGHF/pythia8/generator/pythia8_charmtriggers_with_decays.cfg MC/config/PWGHF/pythia8/decayer/force_hadronic_charmHadTrigger.cfg renamed to MC/config/PWGHF/pythia8/generator/pythia8_charmtriggers_with_decays.cfg

@@ -1,6 +1,19 @@
 ### author: Fabrizio Grosa (fabrizio.grosa@cern.ch)
 ### since:  January 2022
 
+### beams
+Beams:idA 2212			# proton
+Beams:idB 2212 			# proton
+Beams:eCM 13500. 		# GeV
+
+### processes
+HardQCD:hardccbar on		# scatterings g-g / q-qbar -> c-cbar
+HardQCD:hardbbbar on		# scatterings g-g / q-qbar -> b-bbar
+
+### decays
+ParticleDecays:limitTau0 on	
+ParticleDecays:tau0Max 10. # 1 micron it's enough for all the open charm and beauty hadrons
+
 ### Force golden chamr hadrons decay modes for trigger studies
 
 ### add D+ decays absent in PYTHIA8 decay table and set BRs from PDG for other
@@ -45,6 +58,11 @@
 4122:onMode = off
 4232:onMode = off
 4132:onMode = off
+443:onMode = off
+4332:onMode = off
+
+### D0 -> K pi
+421:onIfMatch = 321 211
 
 ### D+/- -> K pi pi 
 411:onIfMatch = 321 211 211
@@ -53,11 +71,6 @@
 ### D+/- -> phi pi
 411:onIfMatch = 333 211
 
-### D0 -> K pi
-421:onIfMatch = 321 211
-
-### D_s -> K K*
-431:onIfMatch = 321 313
 ### D_s -> Phi pi
 431:onIfMatch = 333 211
 
@@ -69,8 +82,6 @@
 4122:onIfMatch = 3124 211
 ### Lambda_c -> p K pi
 4122:onIfMatch = 2212 321 211
-### Lambda_c -> Lambda pi 
-4122:onIfMatch = 3122 211
 
 ### Xic+ -> pK*0
 4232:onIfMatch = 2212 313
@@ -85,9 +96,7 @@
 4132:onIfMatch = 3312 211
 
 ### Omega_c -> Omega pi
-4332:onMode = off
 4332:onIfMatch = 3334 211
 
 ### Jpsi -> ee
-443:onMode = off
-4132:onIfMatch = 11 11
+443:onIfMatch = 11 11

MC/run/PWGHF/run_pp_HFtriggers_bbbar_userhook.sh

@@ -23,7 +23,7 @@ NTIMEFRAMES=${NTIMEFRAMES:-1}
 # create workflow
 
 #ccbar filter
-${O2DPG_ROOT}/MC/bin/o2dpg_sim_workflow.py -eCM 13500 -col pp -gen pythia8 -proc "inel" -j ${NWORKERS} -ns ${NSIGEVENTS} -tf ${NTIMEFRAMES} -interactionRate 500000 -confKey "Diamond.width[2]=6." -e TGeant4 -mod "--skipModules ZDC" \
+${O2DPG_ROOT}/MC/bin/o2dpg_sim_workflow.py -eCM 13500 -col pp -gen pythia8 -proc "inel" -j ${NWORKERS} -ns ${NSIGEVENTS} -tf ${NTIMEFRAMES} -interactionRate 500000 -confKey "Diamond.width[2]=6.;GeneratorPythia8.config=${O2DPG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmtriggers_with_decays.cfg" -e TGeant4 -mod "--skipModules ZDC" \
         -ini $O2DPG_ROOT/MC/config/PWGHF/ini/GeneratorHFTrigger_bbbar.ini \
 
 # run workflow

MC/run/PWGHF/run_pp_HFtriggers_ccbar_userhook.sh

@@ -23,7 +23,7 @@ NTIMEFRAMES=${NTIMEFRAMES:-1}
 # create workflow
 
 #ccbar filter
-${O2DPG_ROOT}/MC/bin/o2dpg_sim_workflow.py -eCM 13500 -col pp -gen pythia8 -proc "inel" -j ${NWORKERS} -ns ${NSIGEVENTS} -tf ${NTIMEFRAMES} -interactionRate 500000 -confKey "Diamond.width[2]=6." -e TGeant4 -mod "--skipModules ZDC" \
+${O2DPG_ROOT}/MC/bin/o2dpg_sim_workflow.py -eCM 13500 -col pp -gen pythia8 -proc "inel" -j ${NWORKERS} -ns ${NSIGEVENTS} -tf ${NTIMEFRAMES} -interactionRate 500000 -confKey "Diamond.width[2]=6.;GeneratorPythia8.config=${O2DPG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmtriggers_with_decays.cfg" -e TGeant4 -mod "--skipModules ZDC" \
         -ini $O2DPG_ROOT/MC/config/PWGHF/ini/GeneratorHFTrigger_ccbar.ini \
 
 # run workflow