From 12509ad1cac1314423b8dcd9625e0ebed43f5f0d Mon Sep 17 00:00:00 2001 From: romainschotter Date: Fri, 22 May 2026 11:43:15 +0200 Subject: [PATCH 1/3] Add configuration for Pb-Pb with hadronic rescattering --- .../ini/pythia8_PbPb_rescattering_536.ini | 9 +++++++ .../ini/tests/pythia8_PbPb_rescattering_536.C | 24 +++++++++++++++++++ .../pythia8_PbPb_rescattering_536.cfg | 22 +++++++++++++++++ 3 files changed, 55 insertions(+) create mode 100644 MC/config/common/ini/pythia8_PbPb_rescattering_536.ini create mode 100644 MC/config/common/ini/tests/pythia8_PbPb_rescattering_536.C create mode 100644 MC/config/common/pythia8/generator/pythia8_PbPb_rescattering_536.cfg diff --git a/MC/config/common/ini/pythia8_PbPb_rescattering_536.ini b/MC/config/common/ini/pythia8_PbPb_rescattering_536.ini new file mode 100644 index 000000000..c1cd7fb4f --- /dev/null +++ b/MC/config/common/ini/pythia8_PbPb_rescattering_536.ini @@ -0,0 +1,9 @@ +[Diamond] +width[2]=6.0 + +[GeneratorExternal] +fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/ALICE3/pythia8/generator_pythia8_ALICE3.C +funcName=generator_pythia8_ALICE3() + +[GeneratorPythia8] +config=${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/generator/pythia8_PbPb_rescattering_536.cfg diff --git a/MC/config/common/ini/tests/pythia8_PbPb_rescattering_536.C b/MC/config/common/ini/tests/pythia8_PbPb_rescattering_536.C new file mode 100644 index 000000000..1c28040e2 --- /dev/null +++ b/MC/config/common/ini/tests/pythia8_PbPb_rescattering_536.C @@ -0,0 +1,24 @@ +int External() { + std::string path{"o2sim_Kine.root"}; + + TFile file(path.c_str(), "READ"); + if (file.IsZombie()) { + std::cerr << "Cannot open ROOT file " << path << "\n"; + return 1; + } + + auto tree = (TTree *)file.Get("o2sim"); + if (!tree) { + std::cerr << "Cannot find tree o2sim in file " << path << "\n"; + return 1; + } + std::vector *tracks{}; + tree->SetBranchAddress("MCTrack", &tracks); + + auto nEvents = tree->GetEntries(); + if (nEvents == 0) { + std::cerr << "No event of interest\n"; + return 1; + } + return 0; +} diff --git a/MC/config/common/pythia8/generator/pythia8_PbPb_rescattering_536.cfg b/MC/config/common/pythia8/generator/pythia8_PbPb_rescattering_536.cfg new file mode 100644 index 000000000..42469b62d --- /dev/null +++ b/MC/config/common/pythia8/generator/pythia8_PbPb_rescattering_536.cfg @@ -0,0 +1,22 @@ +### OO beams +Beams:idA 1000822080 # Pb +Beams:idB 1000822080 # Pb +Beams:eCM = 5360.0 ### energy + +Beams:frameType = 1 +ParticleDecays:limitTau0 = on +ParticleDecays:tau0Max = 10. ### match alice: 1cm/c = 10.0mm/c + +### Save some CPU at init of jobs +### To avoid refitting, add the following lines to your configuration file: +HeavyIon:SigFitNGen = 0 +HeavyIon:SigFitDefPar = 2.15,18.42,0.33 + +### enable hadronic rescattering +HadronLevel:Rescatter = on # default = off +Fragmentation:setVertices = on # default = off +PartonVertex:setVertex = on # default = off +Rescattering:nearestNeighbours = off # default = on (but "require a larger retuning effort") +Rescattering:inelastic = on # default = on + +Random:setSeed = on From 2438f7c40c4b45cf6e4c1d527cc9002472ccfb06 Mon Sep 17 00:00:00 2001 From: romainschotter Date: Fri, 22 May 2026 12:01:47 +0200 Subject: [PATCH 2/3] Update PbPb Angantyr settings --- .../generator/pythia8_PbPb_rescattering_536.cfg | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/MC/config/common/pythia8/generator/pythia8_PbPb_rescattering_536.cfg b/MC/config/common/pythia8/generator/pythia8_PbPb_rescattering_536.cfg index 42469b62d..d4158bdcf 100644 --- a/MC/config/common/pythia8/generator/pythia8_PbPb_rescattering_536.cfg +++ b/MC/config/common/pythia8/generator/pythia8_PbPb_rescattering_536.cfg @@ -7,10 +7,12 @@ Beams:frameType = 1 ParticleDecays:limitTau0 = on ParticleDecays:tau0Max = 10. ### match alice: 1cm/c = 10.0mm/c -### Save some CPU at init of jobs -### To avoid refitting, add the following lines to your configuration file: -HeavyIon:SigFitNGen = 0 -HeavyIon:SigFitDefPar = 2.15,18.42,0.33 +### Initialize the Angantyr model to fit the total and semi-includive +### cross sections in Pythia within some tolerance. +HeavyIon:SigFitErr = 0.02,0.02,0.1,0.05,0.05,0.0,0.1,0.0 + +### These parameters are typicall suitable for sqrt(S_NN)=5TeV +HeavyIon:SigFitDefPar = 17.24,2.15,0.33,0.0,0.0,0.0,0.0,0.0 ### enable hadronic rescattering HadronLevel:Rescatter = on # default = off From f32b1979bdc3943baf6c7b89855ed1b924ed53e3 Mon Sep 17 00:00:00 2001 From: romainschotter Date: Sun, 24 May 2026 11:49:59 +0200 Subject: [PATCH 3/3] Add #NEV_TEST> 1 --- MC/config/common/ini/pythia8_PbPb_rescattering_536.ini | 1 + 1 file changed, 1 insertion(+) diff --git a/MC/config/common/ini/pythia8_PbPb_rescattering_536.ini b/MC/config/common/ini/pythia8_PbPb_rescattering_536.ini index c1cd7fb4f..1ed1278ef 100644 --- a/MC/config/common/ini/pythia8_PbPb_rescattering_536.ini +++ b/MC/config/common/ini/pythia8_PbPb_rescattering_536.ini @@ -1,3 +1,4 @@ +#NEV_TEST> 1 [Diamond] width[2]=6.0