Keep only charm baryon decay modes

Updated decay modes for charm baryons and added new channels for Omega_c decays.

Author: FDUEnrich

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

@@ -45,62 +45,10 @@ BeamRemnants:saturation 5
 # OmegaC
 4332:tau0 = 0.0803
 
-### Force golden charm hadrons decay modes for D2H studies
-# HF decays
-### BR are set to yield 50% of signal in golden channel and uniform abundance of corr. bkg channels (weighted by BR from PDG) 
-###
-### D0 decays
-### D0 -> K- π+ (50%)
-421:oneChannel = 1 0.50000 0 -321 211              ### D0 -> K- π+ 3.94% 
-### D0 -> K- π+ π0 (12.50%)
-421:addChannel = 1 0.00625 0 -321 211 111          ### D0 -> K- π+ π0 (non-resonant) 1.15% (e.g. 115/(115+231+195+1120)*0.2)
-421:addChannel = 1 0.08750 0  213 -321             ### D0 -> rho+ K- 11.2%
-421:addChannel = 1 0.01250 0 -313 111              ### D0 -> antiK*0(892) π0 1.95%
-421:addChannel = 1 0.01875 0 -323 211              ### D0 -> K*-(892) π+ 2.31%
-### D0 -> π- π+ (12.50%)
-421:addChannel = 1 0.12500 0 -211 211              ### D0 -> π- π+ (non-resonant) 1.0e-4
-### D0 -> π- π+ π0 (12.50%)
-421:addChannel = 1 0.08750 0  213 -211             ### D0 -> rho+ π- 1.01%
-421:addChannel = 1 0.03750 0 -211 211 111          ### D0 -> π- π+ π0 (non-resonant) 1.3e-4
-### D0 -> K- K+ (12.50%)
-421:addChannel = 1 0.12500 0 -321 321              ### D0 -> K- K+ (non-resonant) 4.08e-3
+### Force charm baryon decay modes (Ξc and Ωc)
+### D mesons use default PYTHIA8 branching ratios
 
-### D+ decays
-### D+ -> K- π+ π+ (50%)
-411:oneChannel = 1 0.40189 0 -321 211 211          ### D+ -> K- π+ π+ 9.38%
-411:addChannel = 1 0.05356 0 -10311 211            ### D+ -> antiK*0(1430) π+ 1.25%
-411:addChannel = 1 0.04455 0 -313 211              ### D+ -> K*0(892) pi+ 1.04%
-### D+ -> K- π+ π+ π0 (small, 3%)
-411:addChannel = 1 0.03000 0 -321 211 211 111      ### D+ -> K- π+ π+ π0 6.25%
-### D+ -> K- K+ π+ (36.00%, set 25% for D+ -> φ π+, 11% for the rest)
-411:addChannel = 1 0.25000 0  333 211              ### D+ -> φ π+ 2.69e-3 !needed for signal
-411:addChannel = 1 0.03929 0 -313 321              ### D+ -> K*0(892) K+ 2.49e-3
-411:addChannel = 1 0.02865 0 -10311 321            ### D+ -> antiK*0(1430) K+ 1.82e-3
-411:addChannel = 1 0.04206 0 -321 321 211          ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 
-### D+ -> π- π+ π+ (11.00%)
-411:addChannel = 1 0.02911 0  113 211               ### D+ -> rho0 π+ 8.4e-4
-411:addChannel = 1 0.01618 0  225 211               ### D+ -> f2(1270) π+ 4.6e-4
-411:addChannel = 1 0.06471 0 -211 211 211           ### D+ -> π- π+ π+ (non-resonant) 1.0e-4
-
-### Ds+ decays
-### Ds+ -> K- K+ π+ (50%)
-431:oneChannel = 1 0.50000 0  333 211             ### Ds+ -> φ(1020) π+ 2.21%
-431:addChannel = 1 0.15000 0 -313 321             ### Ds+ -> antiK*(892) K+ 2.58%
-### Ds+ -> K- K+ π+ π0 (small, 2%)
-431:addChannel = 1 0.02000 0  333 213              ### Ds+ -> φ(1020) rho 5.50%
-### Ds+ -> π- π+ π+ (11.00%)
-431:addChannel = 1 0.00220 0  113 211             ### Ds+ -> rho π+ 1.1e-4
-431:addChannel = 1 0.00220 0  225 211             ### Ds+ -> f2(1270) π+ 1.4e-3
-431:addChannel = 1 0.10560 0 -211 211 211         ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave)
-### Ds+ -> π- K+ π+ (11.00%)
-431:addChannel = 1 0.03080 0  313 211             ### Ds+ -> K*(892)0 π+ 1.67e-3
-431:addChannel = 1 0.02200 0  10221 321           ### Ds+ -> f0(1370) K+ 1.2e-3
-431:addChannel = 1 0.03960 0  113 321             ### Ds+ -> rho0 K+ 2.17e-3
-431:addChannel = 1 0.01760 0 -211 321 211         ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3
-### Ds+ -> π+ π- π+ π0 (11.00%)
-431:addChannel = 1 0.11000 0  221 211             ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel)
-
-## Λc decays
+## Λc decays (optional, can be kept for reference)
 ### Λc+ -> p K- π+ (36%)
 4122:oneChannel = 1 0.14400 0 2212 -321 211       ### Λc+ -> p K- π+ (non-resonant) 3.5%
 4122:addChannel = 1 0.08100 100 2212 -313         ### Λc+ -> p K*0(892) 1.96%
@@ -135,40 +83,19 @@ BeamRemnants:saturation 5
 4332:oneChannel = 1 0.5 0 3334 211
 4332:addChannel = 1 0.5 0 3312 211
 
-# Allow the decay of resonances in the decay chain
-### K*0(892) -> K- π+
-313:onMode = off
-313:onIfAll = 321 211
-### K*(892)+ -> K- π0
-323:onMode = off
-323:onIfAll = 321 111
-### K*(1430)0 -> K- π+
-10311:onMode = off
-10311:onIfAll = 321 211
-### rho+ -> π+ π0
-213:onMode = off
-213:onIfAll = 211 111
-### φ -> K+ K-
-333:onMode = off
-333:onIfAll = 321 321
-### rho0 -> π+ π-
-113:onMode = off
-113:onIfAll = 211 211
-### f2(1270) -> π+ π-
-225:onMode = off
-225:onIfAll = 211 211
-### f0(1370) -> π+ π-
-10221:onMode = off
-10221:onIfAll = 211 211
-### eta -> π+ π-
-221:onMode = off
-221:onIfAll =  111 211 211
+# Allow the decay of resonances in the decay chain (only for charm baryons)
 ### for Λc -> Delta++ K-
 2224:onMode = off
 2224:onIfAll = 2212 211
 ### for Λc -> Lambda(1520) K-
 102134:onMode = off
 102134:onIfAll = 2212 321
+### for Λc -> p K*0(892)
+313:onMode = off
+313:onIfAll = 321 211
+### for Λc -> p phi
+333:onMode = off
+333:onIfAll = 321 321
 ### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p 
 ### and Omega_c -> pi Xi -> pi pi Lambda -> pi pi pi p
 3312:onMode = off
@@ -179,82 +106,13 @@ BeamRemnants:saturation 5
 3334:onMode = off
 3334:onIfAll = 3122 -321
 
-# Switch off all decay channels
-411:onMode = off
-421:onMode = off
-431:onMode = off
+# Switch off charm baryon decay channels (D mesons use PYTHIA8 defaults)
 4122:onMode = off
 4232:onMode = off
 4132:onMode = off
-443: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 -> rho K
-421:onIfMatch = 213 321
-### D0 -> antiK*0 π0
-421:onIfMatch = 313 111
-### D0 -> K*- π+
-421:onIfMatch = 323 211 
-### D0 -> π π
-421:onIfMatch = 211 211
-### D0 -> rho+ π-
-421:onIfMatch = 213 -211
-### D0 -> π π π0
-421:onIfMatch = 211 211 111
-### D0 -> K K
-421:onIfMatch = 321 321
-
-### D+/- -> K π π 
-411:onIfMatch = 321 211 211
-### D+/- -> K π π 
-411:onIfMatch = 321 211 211 111
-### D+/- -> K* π
-411:onIfMatch = 313 211
-### D+/- -> K* K
-411:onIfMatch = 313 321
-### D+/- -> antiK* π
-411:onIfMatch = 10311 211
-### D+/- -> antiK* K
-411:onIfMatch = 10311 321
-### D+/- -> φ π
-411:onIfMatch = 333 211
-### D+/- -> K K π 
-411:onIfMatch = 321 321 211
-### D+/- -> f2(1270) π
-411:onIfMatch = 225 211
-### D+/- -> rho π
-411:onIfMatch = 113 211
-### D+/- -> π π π
-411:onIfMatch = 211 211 211
-
-### Ds -> φ π
-431:onIfMatch = 333 211
-### Ds -> K* K
-431:onIfMatch = 313 321
-### Ds -> φ rho+
-431:onIfMatch = 333 213
-### Ds -> rho π
-431:onIfMatch = 113 211
-### Ds -> f2 π
-431:onIfMatch = 225 211
-### Ds -> π π π
-431:onIfMatch = 211 211 211
-### Ds -> K*(892)0 π+
-431:onIfMatch = 313 211
-### Ds -> f0(1370) K+
-431:onIfMatch = 10221 321
-### Ds -> rho0 K+
-431:onIfMatch = 113 321
-### Ds -> K π π
-431:onIfMatch =  321 211 211
-### Ds -> eta π+
-431:onIfMatch = 221 211
-
+# Allow the decay of charm baryons
 ### Λc -> pK0s
 4122:onIfMatch = 2212 311 
 ### Λc -> p K- π+ π0
@@ -294,3 +152,5 @@ BeamRemnants:saturation 5
 4332:onIfMatch = 3334 211
 ### Omega_c -> Xi pi
 4332:onIfMatch = 3312 211
+### Omega_c -> Omega- Ka+
+4332:onIfMatch = 3334 321