From 609476977f45ef8c0bdda0a122c090dba5bfe801 Mon Sep 17 00:00:00 2001 From: stefanopolitano Date: Tue, 15 Apr 2025 17:26:32 +0200 Subject: [PATCH 1/3] [PWGHF] Add config and ini file for correlated bkg production in pp --- ...D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.ini | 8 + ...hia8_charmhadronic_with_decays_CorrBkg.cfg | 159 ++++++++++++++++++ 2 files changed, 167 insertions(+) create mode 100644 MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.ini create mode 100644 MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg diff --git a/MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.ini b/MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.ini new file mode 100644 index 000000000..978cf5499 --- /dev/null +++ b/MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.ini @@ -0,0 +1,8 @@ +### The external generator derives from GeneratorPythia8. +[GeneratorExternal] +fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C +funcName=GeneratorPythia8GapTriggeredCharmAndBeauty(1, -1.5, 1.5) + +[GeneratorPythia8] +config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg +includePartonEvent=true diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg new file mode 100644 index 000000000..54756cc71 --- /dev/null +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg @@ -0,0 +1,159 @@ +### authors: Stefano Politano (stefano.politano@cern.ch) +### last update: April 2025 + +### beams +Beams:idA 2212 # proton +Beams:idB 2212 # proton +Beams:eCM 13600. # GeV + +### processes +SoftQCD:inelastic on # all inelastic processes + +### decays +ParticleDecays:limitTau0 on +ParticleDecays:tau0Max 10. + +### switching on Pythia Mode2 +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 + +# Correct decay lengths (wrong in PYTHIA8 decay table) +# Lb +5122:tau0 = 0.4390 +# Xic0 +4132:tau0 = 0.0455 +# OmegaC +4332:tau0 = 0.0803 + +### HF decays +### BR are set to yield the same amount of counts for each final state +### +### D0 decays +421:oneChannel = 1 0.170 0 -321 211 ### D0 -> K- π+ +421:addChannel = 1 0.166 0 -321 211 111 ### D0 -> K- π+ π0 +421:addChannel = 1 0.166 0 -211 211 ### D0 -> π- π+ +421:addChannel = 1 0.166 0 -211 211 111 ### D0 -> π- π+ π0 +421:addChannel = 1 0.166 0 -321 321 ### D0 -> K- K+ +421:addChannel = 1 0.166 0 -321 321 111 ### D0 -> K- K+ π0 + +### D+ decays +411:oneChannel = 1 0.17 0 -321 211 211 ### D+ -> K- π+ π+ +411:addChannel = 1 0.33 0 -211 211 211 ### D+ -> π- π+ π+ +411:addChannel = 1 0.17 0 313 211 ### D+ -> K* π+ -> K- π+ π+ +411:addChannel = 1 0.33 0 333 211 ### D+ -> φ π+ -> K- K+ π+ + +### Ds+ decays +431:oneChannel = 1 0.125 0 333 211 ### Ds+ -> φ π+ -> K- K+ π+ +431:addChannel = 1 0.125 0 313 321 ### Ds+ -> K* K+ -> K- K+ π+ +431:addChannel = 1 0.250 0 -211 211 211 ### Ds+ -> π- π+ π+ +431:addChannel = 1 0.250 0 -211 321 211 ### Ds+ -> π- K+ π+ +431:addChannel = 1 0.250 0 -321 321 321 ### Ds+ -> K- K+ K+ + +## Lc decays +4122:oneChannel = 1 0.0825 0 2212 -321 211 ### Λc+ -> p K- π+ +4122:addChannel = 1 0.0825 100 2212 -313 ### Λc+ -> p K* -> p K- π+ +4122:addChannel = 1 0.0825 100 2224 -321 ### Λc+ -> Delta++ K- -> p K- π+ +4122:addChannel = 1 0.0825 100 102134 211 ### Λc+ -> Lambda(1520) K- -> p K- π+ +4122:addChannel = 1 0.1650 0 2212 211 211 ### Λc+ -> p π+ π+ +4122:addChannel = 1 0.1650 0 2212 311 ### Λc+ -> p K0s -> p π+ π- +4122:addChannel = 1 0.3300 0 2212 -321 321 ### Λc+ -> p K- K+ + +### Ξc+ decays +4232:oneChannel = 1 0.34 0 2212 -321 211 ### Ξc+ -> p K- π+ +4232:addChannel = 1 0.33 0 -3312 211 211 ### Ξc+ -> Ξ- π+ π+ +4232:addChannel = 1 0.33 0 2212 333 ### Ξc+ -> p φ -> p K- K+ + +# Allow the decay of resonances in the decay chain +### K* -> K π +313:onMode = off +313:onIfAll = 321 211 +### for Ds -> φ π+ +333:onMode = off +333:onIfAll = 321 321 +### for D0 -> rho0 π+ k- +113:onMode = off +113:onIfAll = 211 211 +### for Λc -> Delta++ K- +2224:onMode = off +2224:onIfAll = 2212 211 +### for Λc -> Lambda(1520) K- +102134:onMode = off +102134:onIfAll = 2212 321 + +### Switch off all decay channels +411:onMode = off +421:onMode = off +431:onMode = off +4122:onMode = off +4232:onMode = off +4332:onMode = off + +# Allow the decay of HF +### D0 -> K π +421:onIfMatch = 321 211 +### D0 -> K π π0 +421:onIfMatch = 321 211 111 +### D0 -> π π +421:onIfMatch = 211 211 +### D0 -> π π π0 +421:onIfMatch = 211 211 111 +### D0 -> K K +421:onIfMatch = 321 321 +### D0 -> K K π0 +421:onIfMatch = 321 321 111 + +### D+/- -> K π π +411:onIfMatch = 321 211 211 +### D+/- -> K* π +411:onIfMatch = 313 211 +### D+/- -> φ π +411:onIfMatch = 333 211 +### D+/- -> π π π +411:onIfMatch = 211 211 211 + +### Ds -> φ π +431:onIfMatch = 333 211 +### Ds -> K* K +431:onIfMatch = 321 313 +### Ds -> π π π +431:onIfMatch = 211 211 211 +### Ds -> K π π +431:onIfMatch = 321 211 211 +### Ds -> K K K +431:onIfMatch = 321 321 321 + +### Λc -> p K π +4122:onIfMatch = 2212 321 211 +### Λc -> p K* +4122:onIfMatch = 2212 313 +### Λc -> Delta++ K +4122:onIfMatch = 2224 321 +### Λc -> Lambda(1520) π +4122:onIfMatch = 102134 211 +### Λc -> p π π +4122:onIfMatch = 2212 211 211 +### Λc -> pK0s +4122:onIfMatch = 2212 311 +### Λc -> p K K +4122:onIfMatch = 2212 321 321 + +### Ξc+ -> p K- π+ +4232:onIfMatch = 2212 321 211 +### Ξc+ -> Ξ- π+ π+ +4232:onIfMatch = 3312 211 211 +### Ξc+ -> p φ +4232:onIfMatch = 2212 333 From 507b91dc7f58056dd224ae0b154e25277040d150 Mon Sep 17 00:00:00 2001 From: stefanopolitano Date: Tue, 15 Apr 2025 17:29:13 +0200 Subject: [PATCH 2/3] [PWGHF] minor fix in ini file --- .../ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.ini | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.ini b/MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.ini index 978cf5499..49102a6b2 100644 --- a/MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.ini +++ b/MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.ini @@ -1,7 +1,7 @@ ### The external generator derives from GeneratorPythia8. [GeneratorExternal] fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C -funcName=GeneratorPythia8GapTriggeredCharmAndBeauty(1, -1.5, 1.5) +funcName=GeneratorPythia8GapTriggeredCharmAndBeauty(5, -1.5, 1.5) [GeneratorPythia8] config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg From e4cbba817c4160379d01c82df9c90a3fef1c977b Mon Sep 17 00:00:00 2001 From: stefanopolitano Date: Tue, 15 Apr 2025 17:38:55 +0200 Subject: [PATCH 3/3] [PWGHF] minor fix in config file --- .../generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg | 1 - 1 file changed, 1 deletion(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg index 54756cc71..4e39ac46f 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg @@ -100,7 +100,6 @@ BeamRemnants:saturation 5 431:onMode = off 4122:onMode = off 4232:onMode = off -4332:onMode = off # Allow the decay of HF ### D0 -> K π