From 4dfdf370d2b5f7b7b3ad3d97cb38528fa763f8c1 Mon Sep 17 00:00:00 2001 From: Laura Gansbartl Date: Thu, 5 Feb 2026 14:44:10 +0100 Subject: [PATCH 1/5] first comit --- .../ini/pythia8_pp_13600_ForcedDalitz.ini | 5 +++++ .../PWGEM/pythia8/decayer/force_dummy.cfg | 0 .../generator/configPythia_ForcedDalitz.cfg | 21 +++++++++++++++++++ 3 files changed, 26 insertions(+) create mode 100644 MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini create mode 100644 MC/config/PWGEM/pythia8/decayer/force_dummy.cfg create mode 100644 MC/config/PWGEM/pythia8/generator/configPythia_ForcedDalitz.cfg diff --git a/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini b/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini new file mode 100644 index 000000000..5d6950977 --- /dev/null +++ b/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini @@ -0,0 +1,5 @@ +[GeneratorPythia8] +config=/data/lgansbartl/alice/O2DPG/MC/config/PWGEM/pythia8/configPythia_ForcedDalitz.cfg + +[DecayerPythia8] +config[0]=/data/lgansbartl/alice/O2DPG/MC/config/PWGEM/pythia8/decayer/force_dummy.cfg diff --git a/MC/config/PWGEM/pythia8/decayer/force_dummy.cfg b/MC/config/PWGEM/pythia8/decayer/force_dummy.cfg new file mode 100644 index 000000000..e69de29bb diff --git a/MC/config/PWGEM/pythia8/generator/configPythia_ForcedDalitz.cfg b/MC/config/PWGEM/pythia8/generator/configPythia_ForcedDalitz.cfg new file mode 100644 index 000000000..79d78de07 --- /dev/null +++ b/MC/config/PWGEM/pythia8/generator/configPythia_ForcedDalitz.cfg @@ -0,0 +1,21 @@ +ProcessLevel:all on + +### beams +Beams:idA 2212 # proton +Beams:idB 2212 # proton +Beams:eCM 13600. # GeV + +### processes +SoftQCD:inelastic on # all inelastic processes + +ParticleDecays:limitTau0 on +ParticleDecays:tau0Max 10 + + +#decay eta -> e+ e- gamma BR = 100% +221:onMode = off +221:onIfMatch = 22 11 -11 +#221:addChannel = 1 1 0 22 11 -11 + +#decay pi0 -> e+ e- gamma BR = 10% +111:addChannel = 1 0.10 0 22 11 -11" \ No newline at end of file From 13a288df708857f61dbd6e0973dabcf84e3c6ab7 Mon Sep 17 00:00:00 2001 From: Laura Gansbartl Date: Thu, 5 Feb 2026 15:23:40 +0100 Subject: [PATCH 2/5] Fix missing new line --- MC/config/PWGEM/pythia8/generator/configPythia_ForcedDalitz.cfg | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/MC/config/PWGEM/pythia8/generator/configPythia_ForcedDalitz.cfg b/MC/config/PWGEM/pythia8/generator/configPythia_ForcedDalitz.cfg index 79d78de07..d51c6f20e 100644 --- a/MC/config/PWGEM/pythia8/generator/configPythia_ForcedDalitz.cfg +++ b/MC/config/PWGEM/pythia8/generator/configPythia_ForcedDalitz.cfg @@ -18,4 +18,4 @@ ParticleDecays:tau0Max 10 #221:addChannel = 1 1 0 22 11 -11 #decay pi0 -> e+ e- gamma BR = 10% -111:addChannel = 1 0.10 0 22 11 -11" \ No newline at end of file +111:addChannel = 1 0.10 0 22 11 -11" From 03c6b4b323b6173f3314ebaacad63af75db3e0ec Mon Sep 17 00:00:00 2001 From: Laura Gansbartl Date: Fri, 6 Feb 2026 14:20:51 +0100 Subject: [PATCH 3/5] Add Pythia8 Forced Dalitz test config for pp 13.6 TeV --- .../ini/tests/pythia8_pp_13600_ForcedDalitz.C | 57 +++++++++++++++++++ 1 file changed, 57 insertions(+) create mode 100644 MC/config/PWGEM/ini/tests/pythia8_pp_13600_ForcedDalitz.C diff --git a/MC/config/PWGEM/ini/tests/pythia8_pp_13600_ForcedDalitz.C b/MC/config/PWGEM/ini/tests/pythia8_pp_13600_ForcedDalitz.C new file mode 100644 index 000000000..be5551f8e --- /dev/null +++ b/MC/config/PWGEM/ini/tests/pythia8_pp_13600_ForcedDalitz.C @@ -0,0 +1,57 @@ +int External() { + std::string path{"o2sim_Kine.root"}; + // Check that file exists, can be opened and has the correct tree + 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); + + // Check if all events are filled + auto nEvents = tree->GetEntries(); + for (Long64_t i = 0; i < nEvents; ++i) + { + tree->GetEntry(i); + if (tracks->empty()) + { + std::cerr << "Empty entry found at event " << i << "\n"; + return 1; + } + } + // // check if each event has at least two oxygen ions + // for (int i = 0; i < nEvents; i++) + // { + // auto check = tree->GetEntry(i); + // int count = 0; + // for (int idxMCTrack = 0; idxMCTrack < tracks->size(); ++idxMCTrack) + // { + // auto track = tracks->at(idxMCTrack); + // if (track.GetPdgCode() == 1000080160) + // { + // count++; + // } + // } + // if (count < 2) + // { + // std::cerr << "Event " << i << " has less than 2 oxygen ions\n"; + // return 1; + // } + // } + + return 0; + } + + int pythia8() + { + return External(); + } + \ No newline at end of file From cbe25fb984103594de965d0228f52f53972c0dbd Mon Sep 17 00:00:00 2001 From: Laura Gansbartl Date: Fri, 6 Feb 2026 15:37:02 +0100 Subject: [PATCH 4/5] update --- .../ini/pythia8_pp_13600_ForcedDalitz.ini | 4 ++-- .../ini/tests/pythia8_pp_13600_ForcedDalitz.C | 20 +------------------ 2 files changed, 3 insertions(+), 21 deletions(-) diff --git a/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini b/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini index 5d6950977..dacad9fdf 100644 --- a/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini +++ b/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini @@ -1,5 +1,5 @@ [GeneratorPythia8] -config=/data/lgansbartl/alice/O2DPG/MC/config/PWGEM/pythia8/configPythia_ForcedDalitz.cfg +config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/configPythia_ForcedDalitz.cfg [DecayerPythia8] -config[0]=/data/lgansbartl/alice/O2DPG/MC/config/PWGEM/pythia8/decayer/force_dummy.cfg +config[0]=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_dummy.cfg diff --git a/MC/config/PWGEM/ini/tests/pythia8_pp_13600_ForcedDalitz.C b/MC/config/PWGEM/ini/tests/pythia8_pp_13600_ForcedDalitz.C index be5551f8e..4829e48d5 100644 --- a/MC/config/PWGEM/ini/tests/pythia8_pp_13600_ForcedDalitz.C +++ b/MC/config/PWGEM/ini/tests/pythia8_pp_13600_ForcedDalitz.C @@ -27,25 +27,7 @@ int External() { return 1; } } - // // check if each event has at least two oxygen ions - // for (int i = 0; i < nEvents; i++) - // { - // auto check = tree->GetEntry(i); - // int count = 0; - // for (int idxMCTrack = 0; idxMCTrack < tracks->size(); ++idxMCTrack) - // { - // auto track = tracks->at(idxMCTrack); - // if (track.GetPdgCode() == 1000080160) - // { - // count++; - // } - // } - // if (count < 2) - // { - // std::cerr << "Event " << i << " has less than 2 oxygen ions\n"; - // return 1; - // } - // } + return 0; } From 7b8122c2e083112f715a719eaf8bfed2d0a91567 Mon Sep 17 00:00:00 2001 From: Laura Gansbartl Date: Fri, 6 Feb 2026 15:48:57 +0100 Subject: [PATCH 5/5] update --- MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini b/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini index dacad9fdf..17caba124 100644 --- a/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini +++ b/MC/config/PWGEM/ini/pythia8_pp_13600_ForcedDalitz.ini @@ -1,5 +1,5 @@ [GeneratorPythia8] -config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/configPythia_ForcedDalitz.cfg +config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/configPythia_ForcedDalitz.cfg [DecayerPythia8] config[0]=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_dummy.cfg