From b705e8a404c181f6221b99842178d0901d14b2d9 Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Wed, 29 Jul 2026 17:33:29 +0200 Subject: [PATCH 1/4] Add configs to force decays of Xic for MC productions for triggered data --- .../decayer/force_hadronic_Xic_trigger.cfg | 31 +++++++++ .../geant4_externaldecayer_Xic_trigger.in | 68 +++++++++++++++++++ 2 files changed, 99 insertions(+) create mode 100644 MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg create mode 100644 MC/config/PWGHF/pythia8/decayer/geant4_externaldecayer_Xic_trigger.in diff --git a/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg b/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg new file mode 100644 index 000000000..34a192f5b --- /dev/null +++ b/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg @@ -0,0 +1,31 @@ +# Decay configuration taken from MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +Init:showChangedParticleData = on +4132:tau0 = 0.0455 # Xic0 + +# switch off all decay channels +4232:onMode = off +4132:onMode = off + +## Xic decays +### Ξc+ -> Ξ- π+ π+ (70%) +4232:addChannel = 1 0.35000 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% +### Ξc+ -> Ξ*0 π+ (30%) +4232:addChannel = 1 0.12500 0 3324 211 + +### Ξc+ -> Ξ*0 π+, Ξ*0 -> Ξ- π+ +4232:onIfMatch = 3324 211 +### Ξc+ -> Ξ- π+ π+ +4232:onIfMatch = 3312 211 211 + +## Xic0 decays +### add Xic0 decays absent in PYTHIA8 decay table +4132:oneChannel = 1 0.0143 0 3312 211 ### Xi_c()0 --> Xi- pi+ 0.01432524810113947 + +### Xic0 -> Xi- pi+ +4132:onIfMatch = 3312 211 +# Matching Exclusive Decay Channels: +# ------------------------------------------------------- +# Xi_c()0 --> Xi- pi+ 0.01432524810113947 +# -> PDG Codes: 3312, 211 +# Xi_c()0 --> Xi- pi+ pi+ pi- 0.04775082700379833 +# -> PDG Codes: 3312, 211, 211, -211 diff --git a/MC/config/PWGHF/pythia8/decayer/geant4_externaldecayer_Xic_trigger.in b/MC/config/PWGHF/pythia8/decayer/geant4_externaldecayer_Xic_trigger.in new file mode 100644 index 000000000..0ecc07cee --- /dev/null +++ b/MC/config/PWGHF/pythia8/decayer/geant4_externaldecayer_Xic_trigger.in @@ -0,0 +1,68 @@ + +/control/verbose 2 +/mcVerbose/all 1 +/mcVerbose/geometryManager 1 +/mcVerbose/opGeometryManager 1 +/mcTracking/loopVerbose 1 +/mcVerbose/composedPhysicsList 2 +/mcVerbose/runAction 2 # For looping thresholds control +#/tracking/verbose 1 +#//control/cout/ignoreThreadsExcept 0 + +/mcPhysics/rangeCuts 0.001 mm +/mcRegions/setRangePrecision 5 +/mcTracking/skipNeutrino true +/mcDet/setIsMaxStepInLowDensityMaterials true +/mcDet/setMaxStepInLowDensityMaterials 10 m +/mcMagField/setConstDistance 1 mm +/mcDet/setIsZeroMagField true +/mcControl/useRootRandom true # couple G4 random seed to gRandom + +# optical + +/process/optical/verbose 0 +/process/optical/processActivation Scintillation 0 +/process/optical/processActivation OpWLS 0 +/process/optical/processActivation OpMieHG 0 +/process/optical/cerenkov/setTrackSecondariesFirst false +/mcMagField/stepperType NystromRK4 + +# PAI for TRD +# Geant4 VMC >= v3.2 +/mcPhysics/emModel/setEmModel PAI +/mcPhysics/emModel/setRegions TRD_Gas-mix +/mcPhysics/emModel/setParticles all +/mcPrimaryGenerator/skipUnknownParticles true # don't crash when seeing unknown ion etc. (issue warning) + +# +# Precise Msc for EMCAL +# +# Geant4 VMC >= v3.2 +/mcPhysics/emModel/setEmModel SpecialUrbanMsc +/mcPhysics/emModel/setRegions EMC_Lead$ EMC_Scintillator$ +/mcPhysics/emModel/setParticles e- e+ + +# combined transportation + Msc mode is currently broken for ALICE (Geant 10.2.0) +/process/em/transportationWithMsc Disabled + +# +# Adding extra lines for fixing tracking bias +# +/mcMagField/setDeltaIntersection 1.0e-05 mm +/mcMagField/setMinimumEpsilonStep 0.5e-05 +/mcMagField/setMaximumEpsilonStep 1.0e-05 +/mcMagField/printParameters + +# Change default parameters for killing looping particles +# +/mcPhysics/useHighLooperThresholds +/mcRun/setLooperThresholdImportantEnergy 100. MeV + +# Define media with the INCLXX physics list; here basically in all ITS media +#/mcVerbose/biasingConfigurationManager 3 +/mcPhysics/biasing/setModel inclxx +/mcPhysics/biasing/setRegions ITS_AIR$ ITS_WATER$ ITS_COPPER$ ITS_KAPTON(POLYCH2)$ ITS_GLUE_IBFPC$ ITS_CERAMIC$ ITS_K13D2U2k$ ITS_K13D2U120$ ITS_F6151B05M$ ITS_M60J3K$ ITS_M55J6K$ ITS_FGS003$ ITS_CarbonFleece$ ITS_PEEKCF30$ ITS_GLUE$ ITS_ALUMINUM$ ITS_INOX304$ ALPIDE_METALSTACK$ ALPIDE_SI$ +/mcPhysics/biasing/setParticles proton neutron pi+ pi- + +# external decayer +/mcPhysics/setExtDecayerSelection xi_c0 anti_xi_c0 xi_c+ xi_c- From dbcfcfb43e99bd5c4c3a63c09b1a920232f2e143 Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Thu, 30 Jul 2026 10:37:46 +0200 Subject: [PATCH 2/4] Modify decay tables for 2025 and 2026 charm baryon MCs --- .../decayer/force_hadronic_charmbaryon.cfg | 22 +++--- .../force_hadronic_charmbaryon_trigger.cfg | 70 +++++++++++++++++++ 2 files changed, 83 insertions(+), 9 deletions(-) create mode 100644 MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon_trigger.cfg diff --git a/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon.cfg b/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon.cfg index 9e5e79e45..03534d1c2 100644 --- a/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon.cfg +++ b/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon.cfg @@ -9,17 +9,17 @@ Init:showChangedParticleData = on 4132:onMode = off ## Xic decays -### Ξc+ -> p K- π+ (35%) -4232:oneChannel = 1 0.17500 0 2212 -321 211 ### Ξc+ -> p K- π+ 6.18e-3 -4232:addChannel = 1 0.17500 0 2212 -313 ### Ξc+ -> p antiK*0(892) -### Ξc+ -> Ξ- π+ π+ (35%) (set the same as Ξc+ -> p K- π+) -4232:addChannel = 1 0.35000 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% +### Ξc+ -> p K- π+ (30%) +4232:oneChannel = 1 0.1500 0 2212 -321 211 ### Ξc+ -> p K- π+ 6.18e-3 +4232:addChannel = 1 0.1500 0 2212 -313 ### Ξc+ -> p antiK*0(892) +### Ξc+ -> Ξ- π+ π+ (40%) (set the same as Ξc+ -> p K- π+) +4232:addChannel = 1 0.4000 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% ### Ξc+ -> p φ (10%) 4232:addChannel = 1 0.10000 0 2212 333 ### Ξc+ -> p φ ### Ξc+ -> sigma+ π+ π- (10%) -4232:addChannel = 1 0.12500 0 3222 -211 211 ### Ξc+ -> sigma+ π- π+ 1.37% +4232:addChannel = 1 0.10000 0 3222 -211 211 ### Ξc+ -> sigma+ π- π+ 1.37% ### Ξc+ -> Ξ*0 π+ (10%) -4232:addChannel = 1 0.12500 0 3324 211 +4232:addChannel = 1 0.10000 0 3324 211 ### Ξc+ -> p K- π+ 4232:onIfMatch = 2212 321 211 @@ -49,8 +49,9 @@ Init:showChangedParticleData = on ## OmegaC decays ### add custom OmegaC decays absent in PYTHIA8 decay table -4332:oneChannel = 1 0.5 0 3334 211 -4332:addChannel = 1 0.5 0 3312 211 +4332:oneChannel = 1 0.7 0 3334 211 +4332:addChannel = 1 0.1 0 3312 211 +4332:addChannel = 1 0.2 0 3334 321 ### Omega_c -> Omega pi 4332:onIfMatch = 3334 211 @@ -72,6 +73,9 @@ Init:showChangedParticleData = on # Omega_c()0 --> Xi- K- 2 pi+ 0.625932203389831 # -> PDG Codes: 3312, -321, [2], 211 +### Omega_c -> Omega Kaon +4332:onIfMatch = 3334 321 + ## Allow the decay of resonances in the decay chain ### K*0(892) -> K- π+ 313:onMode = off diff --git a/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon_trigger.cfg b/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon_trigger.cfg new file mode 100644 index 000000000..fa5f60330 --- /dev/null +++ b/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon_trigger.cfg @@ -0,0 +1,70 @@ +# Decay configuration taken from MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +Init:showChangedParticleData = on +4332:tau0 = 0.0803 # OmegaC +4132:tau0 = 0.0455 # Xic0 + +# switch off all decay channels +4232:onMode = off +4332:onMode = off +4132:onMode = off + +## Xic decays +### Ξc+ -> Ξ- π+ π+ (40%) (set the same as Ξc+ -> p K- π+) +4232:oneChannel = 1 0.8000 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% +### Ξc+ -> Ξ*0 π+ (10%) +4232:addChannel = 1 0.20000 0 3324 211 + +### Ξc+ -> Ξ*0 π+, Ξ*0 -> Ξ- π+ +4232:onIfMatch = 3324 211 +### Ξc+ -> Ξ- π+ π+ +4232:onIfMatch = 3312 211 211 + +## Xic0 decays +### add Xic0 decays absent in PYTHIA8 decay table +4132:oneChannel = 1 0.0143 0 3312 211 ### Xi_c()0 --> Xi- pi+ 0.01432524810113947 + +### Xic0 -> Xi- pi+ +4132:onIfMatch = 3312 211 +# Matching Exclusive Decay Channels: +# ------------------------------------------------------- +# Xi_c()0 --> Xi- pi+ 0.01432524810113947 +# -> PDG Codes: 3312, 211 +# Xi_c()0 --> Xi- pi+ pi+ pi- 0.04775082700379833 +# -> PDG Codes: 3312, 211, 211, -211 + +## OmegaC decays +### add custom OmegaC decays absent in PYTHIA8 decay table +4332:oneChannel = 1 0.7 0 3334 211 +4332:addChannel = 1 0.1 0 3312 211 +4332:addChannel = 1 0.2 0 3334 321 + +### Omega_c -> Omega pi +4332:onIfMatch = 3334 211 +# Matching Exclusive Decay Channels: +# ------------------------------------------------------- +# Omega_c()0 --> Omega- pi+ None +# -> PDG Codes: 3334, 211 +# Omega_c()0 --> Omega- pi+ pi0 1.79939678284182 +# -> PDG Codes: 3334, 211, 111 +# Omega_c()0 --> Omega- pi- 2 pi+ 0.30954954954955 +# -> PDG Codes: 3334, -211, [2], 211 + +### Omega_c -> Xi pi +4332:onIfMatch = 3312 211 +# Matching Exclusive Decay Channels: +# ------------------------------------------------------- +# Omega_c()0 --> Xi- Kbar0 pi+ 2.12 +# -> PDG Codes: 3312, -311, 211 +# Omega_c()0 --> Xi- K- 2 pi+ 0.625932203389831 +# -> PDG Codes: 3312, -321, [2], 211 + +### Omega_c -> Omega Kaon +4332:onIfMatch = 3334 321 + +## Allow the decay of resonances in the decay chain +### and Omega_c -> pi Xi -> pi pi Lambda -> pi pi pi p +3312:onMode = off +3312:onIfAll = 3122 -211 +### for Omega_c -> pi Omega -> pi K Lambda -> pi K pi p +3334:onMode = off +3334:onIfAll = 3122 -321 From 34cdabb5e57b7609b63f99b13442529175619eb7 Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Thu, 30 Jul 2026 10:45:49 +0200 Subject: [PATCH 3/4] Fix typo --- MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg b/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg index 34a192f5b..de59765a4 100644 --- a/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg +++ b/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg @@ -8,7 +8,7 @@ Init:showChangedParticleData = on ## Xic decays ### Ξc+ -> Ξ- π+ π+ (70%) -4232:addChannel = 1 0.35000 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% +4232:oneChannel = 1 0.35000 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% ### Ξc+ -> Ξ*0 π+ (30%) 4232:addChannel = 1 0.12500 0 3324 211 From d3a4ba8388046e4553e0a61f62014d55e7c0338f Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Thu, 30 Jul 2026 11:26:38 +0200 Subject: [PATCH 4/4] Add Xic->OmegaK --- .../PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg | 4 +++- .../PWGHF/pythia8/decayer/force_hadronic_charmbaryon.cfg | 8 +++++--- .../decayer/force_hadronic_charmbaryon_trigger.cfg | 8 +++++--- 3 files changed, 13 insertions(+), 7 deletions(-) diff --git a/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg b/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg index de59765a4..404f675bc 100644 --- a/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg +++ b/MC/config/PWGHF/pythia8/decayer/force_hadronic_Xic_trigger.cfg @@ -19,7 +19,8 @@ Init:showChangedParticleData = on ## Xic0 decays ### add Xic0 decays absent in PYTHIA8 decay table -4132:oneChannel = 1 0.0143 0 3312 211 ### Xi_c()0 --> Xi- pi+ 0.01432524810113947 +4132:oneChannel = 1 0.7 0 3312 211 ### Xi_c()0 --> Xi- pi+ +4132:addChannel = 1 0.3 0 3334 321 ### Xi_c()0 --> Omega- K+ ### Xic0 -> Xi- pi+ 4132:onIfMatch = 3312 211 @@ -29,3 +30,4 @@ Init:showChangedParticleData = on # -> PDG Codes: 3312, 211 # Xi_c()0 --> Xi- pi+ pi+ pi- 0.04775082700379833 # -> PDG Codes: 3312, 211, 211, -211 +4132:onIfMatch = 3334 321 diff --git a/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon.cfg b/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon.cfg index 03534d1c2..a6128cb23 100644 --- a/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon.cfg +++ b/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon.cfg @@ -36,7 +36,8 @@ Init:showChangedParticleData = on ## Xic0 decays ### add Xic0 decays absent in PYTHIA8 decay table -4132:oneChannel = 1 0.0143 0 3312 211 ### Xi_c()0 --> Xi- pi+ 0.01432524810113947 +4132:oneChannel = 1 0.7 0 3312 211 ### Xi_c()0 --> Xi- pi+ +4132:addChannel = 1 0.3 0 3334 321 ### Xi_c()0 --> Omega- K+ ### Xic0 -> Xi- pi+ 4132:onIfMatch = 3312 211 @@ -46,12 +47,13 @@ Init:showChangedParticleData = on # -> PDG Codes: 3312, 211 # Xi_c()0 --> Xi- pi+ pi+ pi- 0.04775082700379833 # -> PDG Codes: 3312, 211, 211, -211 +4132:onIfMatch = 3334 321 ## OmegaC decays ### add custom OmegaC decays absent in PYTHIA8 decay table -4332:oneChannel = 1 0.7 0 3334 211 +4332:oneChannel = 1 0.8 0 3334 211 4332:addChannel = 1 0.1 0 3312 211 -4332:addChannel = 1 0.2 0 3334 321 +4332:addChannel = 1 0.1 0 3334 321 ### Omega_c -> Omega pi 4332:onIfMatch = 3334 211 diff --git a/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon_trigger.cfg b/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon_trigger.cfg index fa5f60330..5983a65f1 100644 --- a/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon_trigger.cfg +++ b/MC/config/PWGHF/pythia8/decayer/force_hadronic_charmbaryon_trigger.cfg @@ -21,7 +21,8 @@ Init:showChangedParticleData = on ## Xic0 decays ### add Xic0 decays absent in PYTHIA8 decay table -4132:oneChannel = 1 0.0143 0 3312 211 ### Xi_c()0 --> Xi- pi+ 0.01432524810113947 +4132:oneChannel = 1 0.7 0 3312 211 ### Xi_c()0 --> Xi- pi+ +4132:addChannel = 1 0.3 0 3334 321 ### Xi_c()0 --> Omega- K+ ### Xic0 -> Xi- pi+ 4132:onIfMatch = 3312 211 @@ -31,12 +32,13 @@ Init:showChangedParticleData = on # -> PDG Codes: 3312, 211 # Xi_c()0 --> Xi- pi+ pi+ pi- 0.04775082700379833 # -> PDG Codes: 3312, 211, 211, -211 +4132:onIfMatch = 3334 321 ## OmegaC decays ### add custom OmegaC decays absent in PYTHIA8 decay table -4332:oneChannel = 1 0.7 0 3334 211 +4332:oneChannel = 1 0.8 0 3334 211 4332:addChannel = 1 0.1 0 3312 211 -4332:addChannel = 1 0.2 0 3334 321 +4332:addChannel = 1 0.1 0 3334 321 ### Omega_c -> Omega pi 4332:onIfMatch = 3334 211