Update charm-hadron decay lenghts + force Xic and Omegac decay channels

Author: fcatalan92

MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays.cfg

@@ -1,5 +1,7 @@
-### author: Fabrizio Grosa (fabrizio.grosa@cern.ch)
-### since:  January 2022
+### authors: Fabrizio Grosa (fabrizio.grosa@cern.ch)
+###          Cristina Terrevoli (cristina.terrevoli@cern.ch)
+###          Fabio Catalano (fabio.catalano@cern.ch)
+### last update:  November 2023
 
 ### beams
 Beams:idA 2212			# proton
@@ -13,8 +15,15 @@ SoftQCD:inelastic on		# all inelastic processes
 ParticleDecays:limitTau0 on	
 ParticleDecays:tau0Max 10.
 
-### Force golden charm hadrons decay modes 
+# Correct decay lengths (wrong in PYTHIA8 decay table)
+# Lb
+5122:tau0 = 0.4390
+# Xic0
+4132:tau0 = 0.0455
+# OmegaC
+4332:tau0 = 0.0803
 
+### Force golden charm hadrons decay modes for D2H studies
 ### add D+ decays absent in PYTHIA8 decay table and set BRs from PDG for other
 411:oneChannel = 1 0.0752 0 -321 211 211
 411:addChannel = 1 0.0104 0 -313 211
@@ -32,7 +41,9 @@ ParticleDecays:tau0Max 10.
 4232:addChannel = 1 0.2 0 3324 211
 4232:addChannel = 1 0.2 0 3312 211 211
 ### add Xic0 decays absent in PYTHIA8 decay table
-4132:addChannel = 1 0.2 0 3312 211
+4132:addChannel = 1 0.0143 0 3312 211
+### add OmegaC decay absent in PYTHIA8 decay table
+4332:addChannel = 1 1. 0 3334 211
 
 ### K* -> K pi
 313:onMode = off
@@ -49,7 +60,14 @@ ParticleDecays:tau0Max 10.
 ### for Lambda_c -> Lambda(1520) K-
 102134:onMode = off
 102134:onIfAll = 2212 321
-
+### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p
+3312:onMode = off
+3312:onIfAll = 3122 -211
+3122:onMode = off
+3122:onIfAll = 2212 -211
+### for Omega_c -> pi Omega -> pi K Lambda -> pi K pi p
+3334:onMode = off
+3334:onIfAll = 3122 -321
 
 ### switch off all decay channels
 411:onMode = off
@@ -85,7 +103,6 @@ ParticleDecays:tau0Max 10.
 4122:onIfMatch = 3124 211
 ### Lambda_c -> p K pi
 4122:onIfMatch = 2212 321 211
-
 ### Lambda_c -> pK0s
 4122:onIfMatch = 2212 311 
 
@@ -103,9 +120,3 @@ ParticleDecays:tau0Max 10.
 
 ### Omega_c -> Omega pi
 4332:onIfMatch = 3334 211
-
-# Correct Lb decay length (wrong in PYTHIA8 decay table)
-5122:tau0 = 4.41000e-01
-
-# Correct OmegaC decay length (wrong in PYTHIA8 decay table)
-4332:tau0 = 0.08000000000

MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg

@@ -1,5 +1,7 @@
-### author: Fabrizio Grosa (fabrizio.grosa@cern.ch)
-### since:  January 2022
+### authors: Fabrizio Grosa (fabrizio.grosa@cern.ch)
+###          Cristina Terrevoli (cristina.terrevoli@cern.ch)
+###          Fabio Catalano (fabio.catalano@cern.ch)
+### last update:  November 2023
 
 ### beams
 Beams:idA 2212			# proton
@@ -30,10 +32,16 @@ StringFlav:probQQ1toQQ0join 0.0275,0.0275,0.0275,0.0275
 MultiPartonInteractions:pT0Ref 2.15
 BeamRemnants:remnantMode 1
 BeamRemnants:saturation 5
- 
 
-### Force golden charm hadrons decay modes for trigger studies
+# Correct decay lengths (wrong in PYTHIA8 decay table)
+# Lb
+5122:tau0 = 0.4390
+# Xic0
+4132:tau0 = 0.0455
+# OmegaC
+4332:tau0 = 0.0803
 
+### Force golden charm hadrons decay modes for D2H studies
 ### add D+ decays absent in PYTHIA8 decay table and set BRs from PDG for other
 411:oneChannel = 1 0.0752 0 -321 211 211
 411:addChannel = 1 0.0104 0 -313 211
@@ -51,7 +59,9 @@ BeamRemnants:saturation 5
 4232:addChannel = 1 0.2 0 3324 211
 4232:addChannel = 1 0.2 0 3312 211 211
 ### add Xic0 decays absent in PYTHIA8 decay table
-4132:addChannel = 1 0.2 0 3312 211
+4132:addChannel = 1 0.0143 0 3312 211
+### add OmegaC decay absent in PYTHIA8 decay table
+4332:addChannel = 1 1. 0 3334 211
 
 ### K* -> K pi
 313:onMode = off
@@ -68,6 +78,14 @@ BeamRemnants:saturation 5
 ### for Lambda_c -> Lambda(1520) K-
 102134:onMode = off
 102134:onIfAll = 2212 321
+### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p
+3312:onMode = off
+3312:onIfAll = 3122 -211
+3122:onMode = off
+3122:onIfAll = 2212 -211
+### for Omega_c -> pi Omega -> pi K Lambda -> pi K pi p
+3334:onMode = off
+3334:onIfAll = 3122 -321
 
 ### switch off all decay channels
 411:onMode = off
@@ -120,9 +138,3 @@ BeamRemnants:saturation 5
 
 ### Omega_c -> Omega pi
 4332:onIfMatch = 3334 211
-
-# Correct Lb decay length (wrong in PYTHIA8 decay table)
-5122:tau0 = 4.41000e-01
-
-# Correct OmegaC decay length (wrong in PYTHIA8 decay table)
-4332:tau0 = 0.08000000000