From 83a60b6f79c5966a57a330018a8686a12b12c5ac Mon Sep 17 00:00:00 2001
From: Raphaelle Bailhache <raphaelle.bailhache@cern.ch>
Date: Fri, 9 May 2025 16:25:20 +0200
Subject: [PATCH] HF cocktail for 13 TeV

---
 .../ini/Pythia8_Beauty_Cocktail_pp13.ini      | 13 ++++++++
 .../PWGEM/ini/Pythia8_Charm_Cocktail_pp13.ini | 13 ++++++++
 .../generator/pythia8_hf_cocktail_pp13.cfg    | 30 +++++++++++++++++++
 3 files changed, 56 insertions(+)
 create mode 100644 MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail_pp13.ini
 create mode 100644 MC/config/PWGEM/ini/Pythia8_Charm_Cocktail_pp13.ini
 create mode 100644 MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp13.cfg

diff --git a/MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail_pp13.ini b/MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail_pp13.ini
new file mode 100644
index 000000000..4a0f54f13
--- /dev/null
+++ b/MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail_pp13.ini
@@ -0,0 +1,13 @@
+[GeneratorExternal]
+fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
+funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332")
+
+[GeneratorPythia8]
+config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp13.cfg
+hooksFileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
+hooksFuncName = pythia8_userhooks_bbbar(-12.,12.)
+includePartonEvent=true
+
+[DecayerPythia8]
+config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
+config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg
diff --git a/MC/config/PWGEM/ini/Pythia8_Charm_Cocktail_pp13.ini b/MC/config/PWGEM/ini/Pythia8_Charm_Cocktail_pp13.ini
new file mode 100644
index 000000000..cad2afe89
--- /dev/null
+++ b/MC/config/PWGEM/ini/Pythia8_Charm_Cocktail_pp13.ini
@@ -0,0 +1,13 @@
+[GeneratorExternal]
+fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
+funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332")
+
+[GeneratorPythia8]
+config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp13.cfg
+hooksFileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
+hooksFuncName = pythia8_userhooks_ccbar(-12.,12.)
+includePartonEvent=true
+
+[DecayerPythia8]
+config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
+config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg
diff --git a/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp13.cfg b/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp13.cfg
new file mode 100644
index 000000000..f7884e068
--- /dev/null
+++ b/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp13.cfg
@@ -0,0 +1,30 @@
+### beams
+Beams:idA 2212			# proton
+Beams:idB 2212                  # proton
+Beams:eCM 13000.                # GeV
+
+### processes
+SoftQCD:inelastic on
+
+### decays
+ParticleDecays:limitTau0 on
+ParticleDecays:tau0Max 10.
+
+### switch on color reconnection in mode 2 (https://arxiv.org/pdf/1505.01681.pdf)
+Tune:pp = 14
+ColourReconnection:mode = 1
+ColourReconnection:allowDoubleJunRem = off
+ColourReconnection:m0 = 0.3
+ColourReconnection:allowJunctions = on
+ColourReconnection:junctionCorrection = 1.20
+ColourReconnection:timeDilationMode = 2
+ColourReconnection:timeDilationPar = 0.18
+StringPT:sigma = 0.335
+StringZ:aLund = 0.36
+StringZ:bLund = 0.56
+StringFlav:probQQtoQ = 0.078
+StringFlav:ProbStoUD = 0.2
+StringFlav:probQQ1toQQ0join = 0.0275,0.0275,0.0275,0.0275
+MultiPartonInteractions:pT0Ref = 2.15
+BeamRemnants:remnantMode = 1
+BeamRemnants:saturation =5
\ No newline at end of file