From 096182adf3e56ab72c3b897e66637a96f624d5a3 Mon Sep 17 00:00:00 2001 From: wuctlby Date: Thu, 12 Jun 2025 17:35:33 +0200 Subject: [PATCH 01/11] Add the correlated bkg channel for PbPb@5.36TeV, pp@13.6TeV, and ppref@5.36TeV --- ...ythia8_charmhadronic_with_decays_Mode2.cfg | 224 ++++++++++++++---- ...adronic_with_decays_Mode2_hardQCD_5TeV.cfg | 224 ++++++++++++++---- ...charmhadronic_with_decays_Mode2_pp_ref.cfg | 221 +++++++++++++---- 3 files changed, 532 insertions(+), 137 deletions(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg index aebd1cc63..aaf85b86e 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg @@ -42,44 +42,123 @@ BeamRemnants:saturation 5 4332:tau0 = 0.0803 ### Force golden charm hadrons decay modes for D2H studies -### add D0 decays absent in PYTHIA8 decay table and set BRs from PDG for other -421:oneChannel = 1 0.0389 0 -321 211 -421:addChannel = 1 0.00389 0 -321 211 111 -### 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 -411:addChannel = 1 0.0156 0 311 211 -411:addChannel = 1 0.0752 0 333 211 # to have the same amount of D+->KKpi and D+->Kpipi -## add Lc decays absent in PYTHIA8 decay table and set BRs from PDG for other -4122:oneChannel = 1 0.0196 100 2212 -313 -4122:addChannel = 1 0.0108 100 2224 -321 -4122:addChannel = 1 0.022 100 102134 211 -4122:addChannel = 1 0.035 0 2212 -321 211 -4122:addChannel = 1 0.0159 0 2212 311 -### add Xic+ decays absent in PYTHIA8 decay table -4232:addChannel = 1 0.2 0 2212 -313 -4232:addChannel = 1 0.2 0 2212 -321 211 -4232:addChannel = 1 0.2 0 3324 211 -4232:addChannel = 1 0.2 0 3312 211 211 +### 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 + +### D+ decays +### D+ -> K- π+ π+ (50%) +411:oneChannel = 1 0.50000 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% +### D+ -> K- π+ π+ π0 (small, 3%) +411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% +### D+ -> K- K+ π+ (23.50%) +411:addChannel = 1 0.06053 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 +411:addChannel = 1 0.04415 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 +411:addChannel = 1 0.06552 0 333 211 ### D+ -> φ π+ 2.69e-3 +411:addChannel = 1 0.06480 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 +### D+ -> π- π+ π+ (23.50%) +411:addChannel = 1 0.06221 0 113 211 ### D+ -> rho0 π+ 8.4e-4 +411:addChannel = 1 0.03456 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 +411:addChannel = 1 0.13824 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 + +### Ds+ decays +### Ds+ -> K- K+ π+ (50%) +431:oneChannel = 1 0.25000 0 333 211 ### Ds+ -> φ(1020) π+ 2.21% +431:addChannel = 1 0.25000 0 -313 321 ### Ds+ -> antiK*(892) K+ 2.58% +### Ds+ -> K- K+ π+ π0 (small, 3%) +431:addChannel = 1 0.03000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% +### Ds+ -> π- π+ π+ (15.50%) +431:addChannel = 1 0.00310 0 113 211 ### Ds+ -> rho π+ 1.1e-4 +431:addChannel = 1 0.00310 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 +431:addChannel = 1 0.14880 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) +### Ds+ -> π- K+ π+ (15.50%) +431:addChannel = 1 0.04340 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 +431:addChannel = 1 0.03100 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 +431:addChannel = 1 0.05580 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 +431:addChannel = 1 0.02480 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 +### Ds+ -> π+ π- π+ π0 (15.50%) +431:addChannel = 1 0.15500 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) + +## Λc decays +### Λc+ -> p K- π+ (50%, BR set up to 50%) +4122:oneChannel = 1 0.06472 0 2212 -321 211 ### Λc+ -> p K- π+ (non-resonant) 3.5% +4122:addChannel = 1 0.03624 100 2212 -313 ### Λc+ -> p K*0(892) 1.96% +4122:addChannel = 1 0.01997 100 2224 -321 ### Λc+ -> Delta++ K- 1.08% +4122:addChannel = 1 0.00407 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3 +### Λc+ -> p K- π+ π0 (small, 3%) +4122:addChannel = 1 0.03000 0 2212 -321 211 111 ### Λc+ -> p K- π+ π0 (non-resonant) 4.6% +### Λc+ -> p π- π+ (23.50%) +4122:addChannel = 1 0.23500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59% +### Λc+ -> p K- K+ (23.50%) +4122:addChannel = 1 0.23500 0 2212 333 ### Λc+ -> p phi 1.06% + +## Xic decays +### Ξc+ -> p K- π+ (50%) +4232:oneChannel = 1 0.25 0 2212 -321 211 ### Ξc+ -> p K- π+ 6.18e-3 +4232:addChannel = 1 0.25 0 2212 -313 ### Ξc+ -> p antiK*0(892) +### Ξc+ -> p φ (12.5%) +4232:addChannel = 1 0.125 0 2212 333 ### Ξc+ -> p φ +### Ξc+ -> sigma+ π+ π- (12.5%) +4232:addChannel = 1 0.125 0 3222 -211 211 ### Ξc+ -> sigma+ π- π+ 1.37% +### Ξc+ -> Ξ*0 π+ (12.5%) +4232:addChannel = 1 0.125 0 3324 211 +### Ξc+ -> Ξ- π+ π+ (12.5%) +4232:addChannel = 1 0.125 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% + ### add Xic0 decays absent in PYTHIA8 decay table 4132:addChannel = 1 0.0143 0 3312 211 + ### add OmegaC decays absent in PYTHIA8 decay table 4332:addChannel = 1 0.5 0 3334 211 4332:addChannel = 1 0.5 0 3312 211 -### K* -> K pi +# Allow the decay of resonances in the decay chain +### K*0(892) -> K- π+ 313:onMode = off 313:onIfAll = 321 211 -### for Ds -> Phi pi+ +### 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 -### for D0 -> rho0 pi+ k- +### rho0 -> π+ π- 113:onMode = off 113:onIfAll = 211 211 -### for Lambda_c -> Delta++ K- +### 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 +### for Λc -> Delta++ K- 2224:onMode = off 2224:onIfAll = 2212 211 -### for Lambda_c -> Lambda(1520) K- +### for Λc -> Lambda(1520) K- 102134:onMode = off 102134:onIfAll = 2212 321 ### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p @@ -102,41 +181,94 @@ BeamRemnants:saturation 5 443:onMode = off 4332:onMode = off -### D0 -> K pi +# Allow the decay of HF +### D0 -> K π 421:onIfMatch = 321 211 -### D0 -> K pi pi0 +### 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 pi pi +### D+/- -> K π π 411:onIfMatch = 321 211 211 -### D+/- -> K* pi -411:onIfMatch = 313 211 -### D+/- -> phi pi +### D+/- -> K π π +411:onIfMatch = 321 211 211 111 +### D+/- -> K* K +411:onIfMatch = 313 321 +### D+/- -> antiK* π +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 -### D_s -> K K* -431:onIfMatch = 321 313 -### D_s -> Phi pi +### 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 -### Lambda_c -> p K* +### Λc -> p K π +4122:onIfMatch = 2212 321 211 +### Λc -> p K* 4122:onIfMatch = 2212 313 -### Lambda_c -> Delta K +### Λc -> Delta++ K 4122:onIfMatch = 2224 321 -### Lambda_c -> Lambda(1520) pi +### Λc -> Lambda(1520) π 4122:onIfMatch = 102134 211 -### Lambda_c -> p K pi -4122:onIfMatch = 2212 321 211 -### Lambda_c -> pK0s -4122:onIfMatch = 2212 311 +### Λc -> p K- π+ π0 +4122:onIfMatch = 2212 321 211 111 +### Λc -> p π π +4122:onIfMatch = 2212 211 211 +### Λc -> p K K +4122:onIfMatch = 2212 333 -### Xic+ -> pK*0 -4232:onIfMatch = 2212 313 -### Xic+ -> p K- pi+ +### Ξc+ -> p K- π+ 4232:onIfMatch = 2212 321 211 -### Xic+ -> Xi*0 pi+, Xi*->Xi- pi+ +### Ξc+ -> p antiK*0(892) +4232:onIfMatch = 2212 313 +### Ξc+ -> p φ +4232:onIfMatch = 2212 333 +### Ξc+ -> sigma- π+ π+ +4232:onIfMatch = 3222 211 211 +### Ξc+ -> Ξ*0 π+, Ξ*0 -> Ξ- π+ 4232:onIfMatch = 3324 211 -### Xic+ -> Xi- pi+ pi+ +### Ξc+ -> Ξ- π+ π+ 4232:onIfMatch = 3312 211 211 ### Xic0 -> Xi- pi+ @@ -145,4 +277,4 @@ BeamRemnants:saturation 5 ### Omega_c -> Omega pi 4332:onIfMatch = 3334 211 ### Omega_c -> Xi pi -4332:onIfMatch = 3312 211 +4332:onIfMatch = 3312 211 \ No newline at end of file diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg index cbce46503..24058f40d 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg @@ -46,44 +46,123 @@ BeamRemnants:saturation 5 4332:tau0 = 0.0803 ### Force golden charm hadrons decay modes for D2H studies -### add D0 decays absent in PYTHIA8 decay table and set BRs from PDG for other -421:oneChannel = 1 0.0389 0 -321 211 -421:addChannel = 1 0.00389 0 -321 211 111 -### 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 -411:addChannel = 1 0.0156 0 311 211 -411:addChannel = 1 0.0752 0 333 211 # to have the same amount of D+->KKpi and D+->Kpipi -## add Lc decays absent in PYTHIA8 decay table and set BRs from PDG for other -4122:oneChannel = 1 0.0196 100 2212 -313 -4122:addChannel = 1 0.0108 100 2224 -321 -4122:addChannel = 1 0.022 100 102134 211 -4122:addChannel = 1 0.035 0 2212 -321 211 -4122:addChannel = 1 0.0159 0 2212 311 -### add Xic+ decays absent in PYTHIA8 decay table -4232:addChannel = 1 0.2 0 2212 -313 -4232:addChannel = 1 0.2 0 2212 -321 211 -4232:addChannel = 1 0.2 0 3324 211 -4232:addChannel = 1 0.2 0 3312 211 211 +### 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 + +### D+ decays +### D+ -> K- π+ π+ (50%) +411:oneChannel = 1 0.50000 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% +### D+ -> K- π+ π+ π0 (small, 3%) +411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% +### D+ -> K- K+ π+ (23.50%) +411:addChannel = 1 0.06053 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 +411:addChannel = 1 0.04415 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 +411:addChannel = 1 0.06552 0 333 211 ### D+ -> φ π+ 2.69e-3 +411:addChannel = 1 0.06480 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 +### D+ -> π- π+ π+ (23.50%) +411:addChannel = 1 0.06221 0 113 211 ### D+ -> rho0 π+ 8.4e-4 +411:addChannel = 1 0.03456 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 +411:addChannel = 1 0.13824 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 + +### Ds+ decays +### Ds+ -> K- K+ π+ (50%) +431:oneChannel = 1 0.25000 0 333 211 ### Ds+ -> φ(1020) π+ 2.21% +431:addChannel = 1 0.25000 0 -313 321 ### Ds+ -> antiK*(892) K+ 2.58% +### Ds+ -> K- K+ π+ π0 (small, 3%) +431:addChannel = 1 0.03000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% +### Ds+ -> π- π+ π+ (15.50%) +431:addChannel = 1 0.00310 0 113 211 ### Ds+ -> rho π+ 1.1e-4 +431:addChannel = 1 0.00310 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 +431:addChannel = 1 0.14880 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) +### Ds+ -> π- K+ π+ (15.50%) +431:addChannel = 1 0.04340 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 +431:addChannel = 1 0.03100 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 +431:addChannel = 1 0.05580 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 +431:addChannel = 1 0.02480 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 +### Ds+ -> π+ π- π+ π0 (15.50%) +431:addChannel = 1 0.15500 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) + +## Λc decays +### Λc+ -> p K- π+ (50%, BR set up to 50%) +4122:oneChannel = 1 0.06472 0 2212 -321 211 ### Λc+ -> p K- π+ (non-resonant) 3.5% +4122:addChannel = 1 0.03624 100 2212 -313 ### Λc+ -> p K*0(892) 1.96% +4122:addChannel = 1 0.01997 100 2224 -321 ### Λc+ -> Delta++ K- 1.08% +4122:addChannel = 1 0.00407 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3 +### Λc+ -> p K- π+ π0 (small, 3%) +4122:addChannel = 1 0.03000 0 2212 -321 211 111 ### Λc+ -> p K- π+ π0 (non-resonant) 4.6% +### Λc+ -> p π- π+ (23.50%) +4122:addChannel = 1 0.23500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59% +### Λc+ -> p K- K+ (23.50%) +4122:addChannel = 1 0.23500 0 2212 333 ### Λc+ -> p phi 1.06% + +## Xic decays +### Ξc+ -> p K- π+ (50%) +4232:oneChannel = 1 0.25 0 2212 -321 211 ### Ξc+ -> p K- π+ 6.18e-3 +4232:addChannel = 1 0.25 0 2212 -313 ### Ξc+ -> p antiK*0(892) +### Ξc+ -> p φ (12.5%) +4232:addChannel = 1 0.125 0 2212 333 ### Ξc+ -> p φ +### Ξc+ -> sigma+ π+ π- (12.5%) +4232:addChannel = 1 0.125 0 3222 -211 211 ### Ξc+ -> sigma+ π- π+ 1.37% +### Ξc+ -> Ξ*0 π+ (12.5%) +4232:addChannel = 1 0.125 0 3324 211 +### Ξc+ -> Ξ- π+ π+ (12.5%) +4232:addChannel = 1 0.125 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% + ### add Xic0 decays absent in PYTHIA8 decay table 4132:addChannel = 1 0.0143 0 3312 211 + ### add OmegaC decays absent in PYTHIA8 decay table 4332:addChannel = 1 0.5 0 3334 211 4332:addChannel = 1 0.5 0 3312 211 -### K* -> K pi +# Allow the decay of resonances in the decay chain +### K*0(892) -> K- π+ 313:onMode = off 313:onIfAll = 321 211 -### for Ds -> Phi pi+ +### 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 -### for D0 -> rho0 pi+ k- +### rho0 -> π+ π- 113:onMode = off 113:onIfAll = 211 211 -### for Lambda_c -> Delta++ K- +### 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 +### for Λc -> Delta++ K- 2224:onMode = off 2224:onIfAll = 2212 211 -### for Lambda_c -> Lambda(1520) K- +### for Λc -> Lambda(1520) K- 102134:onMode = off 102134:onIfAll = 2212 321 ### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p @@ -106,41 +185,94 @@ BeamRemnants:saturation 5 443:onMode = off 4332:onMode = off -### D0 -> K pi +# Allow the decay of HF +### D0 -> K π 421:onIfMatch = 321 211 -### D0 -> K pi pi0 +### 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 pi pi +### D+/- -> K π π 411:onIfMatch = 321 211 211 -### D+/- -> K* pi -411:onIfMatch = 313 211 -### D+/- -> phi pi +### D+/- -> K π π +411:onIfMatch = 321 211 211 111 +### D+/- -> K* K +411:onIfMatch = 313 321 +### D+/- -> antiK* π +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 -### D_s -> K K* -431:onIfMatch = 321 313 -### D_s -> Phi pi +### 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 -### Lambda_c -> p K* +### Λc -> p K π +4122:onIfMatch = 2212 321 211 +### Λc -> p K* 4122:onIfMatch = 2212 313 -### Lambda_c -> Delta K +### Λc -> Delta++ K 4122:onIfMatch = 2224 321 -### Lambda_c -> Lambda(1520) pi +### Λc -> Lambda(1520) π 4122:onIfMatch = 102134 211 -### Lambda_c -> p K pi -4122:onIfMatch = 2212 321 211 -### Lambda_c -> pK0s -4122:onIfMatch = 2212 311 +### Λc -> p K- π+ π0 +4122:onIfMatch = 2212 321 211 111 +### Λc -> p π π +4122:onIfMatch = 2212 211 211 +### Λc -> p K K +4122:onIfMatch = 2212 333 -### Xic+ -> pK*0 -4232:onIfMatch = 2212 313 -### Xic+ -> p K- pi+ +### Ξc+ -> p K- π+ 4232:onIfMatch = 2212 321 211 -### Xic+ -> Xi*0 pi+, Xi*->Xi- pi+ +### Ξc+ -> p antiK*0(892) +4232:onIfMatch = 2212 313 +### Ξc+ -> p φ +4232:onIfMatch = 2212 333 +### Ξc+ -> sigma- π+ π+ +4232:onIfMatch = 3222 211 211 +### Ξc+ -> Ξ*0 π+, Ξ*0 -> Ξ- π+ 4232:onIfMatch = 3324 211 -### Xic+ -> Xi- pi+ pi+ +### Ξc+ -> Ξ- π+ π+ 4232:onIfMatch = 3312 211 211 ### Xic0 -> Xi- pi+ @@ -149,4 +281,4 @@ BeamRemnants:saturation 5 ### Omega_c -> Omega pi 4332:onIfMatch = 3334 211 ### Omega_c -> Xi pi -4332:onIfMatch = 3312 211 +4332:onIfMatch = 3312 211 \ No newline at end of file diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index 63ce139a5..330fc840d 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -42,44 +42,122 @@ BeamRemnants:saturation 5 4332:tau0 = 0.0803 ### Force golden charm hadrons decay modes for D2H studies -### add D0 decays absent in PYTHIA8 decay table and set BRs from PDG for other -421:oneChannel = 1 0.0389 0 -321 211 -421:addChannel = 1 0.00389 0 -321 211 111 -### 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 -411:addChannel = 1 0.0156 0 311 211 -411:addChannel = 1 0.0752 0 333 211 # to have the same amount of D+->KKpi and D+->Kpipi -## add Lc decays absent in PYTHIA8 decay table and set BRs from PDG for other -4122:oneChannel = 1 0.0196 100 2212 -313 -4122:addChannel = 1 0.0108 100 2224 -321 -4122:addChannel = 1 0.022 100 102134 211 -4122:addChannel = 1 0.035 0 2212 -321 211 -4122:addChannel = 1 0.0159 0 2212 311 -### add Xic+ decays absent in PYTHIA8 decay table -4232:addChannel = 1 0.2 0 2212 -313 -4232:addChannel = 1 0.2 0 2212 -321 211 -4232:addChannel = 1 0.2 0 3324 211 -4232:addChannel = 1 0.2 0 3312 211 211 +### 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 + +### D+ decays +### D+ -> K- π+ π+ (50%) +411:oneChannel = 1 0.50000 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% +### D+ -> K- π+ π+ π0 (small, 3%) +411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% +### D+ -> K- K+ π+ (23.50%) +411:addChannel = 1 0.06053 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 +411:addChannel = 1 0.04415 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 +411:addChannel = 1 0.06552 0 333 211 ### D+ -> φ π+ 2.69e-3 +411:addChannel = 1 0.06480 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 +### D+ -> π- π+ π+ (23.50%) +411:addChannel = 1 0.06221 0 113 211 ### D+ -> rho0 π+ 8.4e-4 +411:addChannel = 1 0.03456 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 +411:addChannel = 1 0.13824 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 + +### Ds+ decays +### Ds+ -> K- K+ π+ (50%) +431:oneChannel = 1 0.25000 0 333 211 ### Ds+ -> φ(1020) π+ 2.21% +431:addChannel = 1 0.25000 0 -313 321 ### Ds+ -> antiK*(892) K+ 2.58% +### Ds+ -> K- K+ π+ π0 (small, 3%) +431:addChannel = 1 0.03000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% +### Ds+ -> π- π+ π+ (15.50%) +431:addChannel = 1 0.00310 0 113 211 ### Ds+ -> rho π+ 1.1e-4 +431:addChannel = 1 0.00310 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 +431:addChannel = 1 0.14880 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) +### Ds+ -> π- K+ π+ (15.50%) +431:addChannel = 1 0.04340 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 +431:addChannel = 1 0.03100 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 +431:addChannel = 1 0.05580 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 +431:addChannel = 1 0.02480 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 +### Ds+ -> π+ π- π+ π0 (15.50%) +431:addChannel = 1 0.15500 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) + +## Λc decays +### Λc+ -> p K- π+ (50%, BR set up to 50%) +4122:oneChannel = 1 0.06472 0 2212 -321 211 ### Λc+ -> p K- π+ (non-resonant) 3.5% +4122:addChannel = 1 0.03624 100 2212 -313 ### Λc+ -> p K*0(892) 1.96% +4122:addChannel = 1 0.01997 100 2224 -321 ### Λc+ -> Delta++ K- 1.08% +4122:addChannel = 1 0.00407 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3 +### Λc+ -> p K- π+ π0 (small, 3%) +4122:addChannel = 1 0.03000 0 2212 -321 211 111 ### Λc+ -> p K- π+ π0 (non-resonant) 4.6% +### Λc+ -> p π- π+ (23.50%) +4122:addChannel = 1 0.23500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59% +### Λc+ -> p K- K+ (23.50%) +4122:addChannel = 1 0.23500 0 2212 333 ### Λc+ -> p phi 1.06% + +## Xic decays +### Ξc+ -> p K- π+ (50%) +4232:oneChannel = 1 0.25 0 2212 -321 211 ### Ξc+ -> p K- π+ 6.18e-3 +4232:addChannel = 1 0.25 0 2212 -313 ### Ξc+ -> p antiK*0(892) +### Ξc+ -> p φ (12.5%) +4232:addChannel = 1 0.125 0 2212 333 ### Ξc+ -> p φ +### Ξc+ -> sigma+ π+ π- (12.5%) +4232:addChannel = 1 0.125 0 3222 -211 211 ### Ξc+ -> sigma+ π- π+ 1.37% +### Ξc+ -> Ξ*0 π+ (12.5%) +4232:addChannel = 1 0.125 0 3324 211 +### Ξc+ -> Ξ- π+ π+ (12.5%) +4232:addChannel = 1 0.125 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% + ### add Xic0 decays absent in PYTHIA8 decay table 4132:addChannel = 1 0.0143 0 3312 211 ### add OmegaC decays absent in PYTHIA8 decay table 4332:addChannel = 1 0.5 0 3334 211 4332:addChannel = 1 0.5 0 3312 211 -### K* -> K pi +# Allow the decay of resonances in the decay chain +### K*0(892) -> K- π+ 313:onMode = off 313:onIfAll = 321 211 -### for Ds -> Phi pi+ +### 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 -### for D0 -> rho0 pi+ k- +### rho0 -> π+ π- 113:onMode = off 113:onIfAll = 211 211 -### for Lambda_c -> Delta++ K- +### 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 +### for Λc -> Delta++ K- 2224:onMode = off 2224:onIfAll = 2212 211 -### for Lambda_c -> Lambda(1520) K- +### for Λc -> Lambda(1520) K- 102134:onMode = off 102134:onIfAll = 2212 321 ### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p @@ -102,41 +180,94 @@ BeamRemnants:saturation 5 443:onMode = off 4332:onMode = off -### D0 -> K pi +# Allow the decay of HF +### D0 -> K π 421:onIfMatch = 321 211 -### D0 -> K pi pi0 +### 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 pi pi +### D+/- -> K π π 411:onIfMatch = 321 211 211 -### D+/- -> K* pi -411:onIfMatch = 313 211 -### D+/- -> phi pi +### D+/- -> K π π +411:onIfMatch = 321 211 211 111 +### D+/- -> K* K +411:onIfMatch = 313 321 +### D+/- -> antiK* π +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 -### D_s -> K K* -431:onIfMatch = 321 313 -### D_s -> Phi pi +### 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 -### Lambda_c -> p K* +### Λc -> p K π +4122:onIfMatch = 2212 321 211 +### Λc -> p K* 4122:onIfMatch = 2212 313 -### Lambda_c -> Delta K +### Λc -> Delta++ K 4122:onIfMatch = 2224 321 -### Lambda_c -> Lambda(1520) pi +### Λc -> Lambda(1520) π 4122:onIfMatch = 102134 211 -### Lambda_c -> p K pi -4122:onIfMatch = 2212 321 211 -### Lambda_c -> pK0s -4122:onIfMatch = 2212 311 +### Λc -> p K- π+ π0 +4122:onIfMatch = 2212 321 211 111 +### Λc -> p π π +4122:onIfMatch = 2212 211 211 +### Λc -> p K K +4122:onIfMatch = 2212 333 -### Xic+ -> pK*0 -4232:onIfMatch = 2212 313 -### Xic+ -> p K- pi+ +### Ξc+ -> p K- π+ 4232:onIfMatch = 2212 321 211 -### Xic+ -> Xi*0 pi+, Xi*->Xi- pi+ +### Ξc+ -> p antiK*0(892) +4232:onIfMatch = 2212 313 +### Ξc+ -> p φ +4232:onIfMatch = 2212 333 +### Ξc+ -> sigma- π+ π+ +4232:onIfMatch = 3222 211 211 +### Ξc+ -> Ξ*0 π+, Ξ*0 -> Ξ- π+ 4232:onIfMatch = 3324 211 -### Xic+ -> Xi- pi+ pi+ +### Ξc+ -> Ξ- π+ π+ 4232:onIfMatch = 3312 211 211 ### Xic0 -> Xi- pi+ From c3f8714351a79c43bd38643284837ec9521ba3ad Mon Sep 17 00:00:00 2001 From: wuctlby Date: Thu, 12 Jun 2025 19:10:12 +0200 Subject: [PATCH 02/11] correct the abundance ratio of Lc+, and add the Lambdac+ -> p k0s --- ...ythia8_charmhadronic_with_decays_Mode2.cfg | 20 ++++++++++--------- ...adronic_with_decays_Mode2_hardQCD_5TeV.cfg | 20 ++++++++++--------- ...charmhadronic_with_decays_Mode2_pp_ref.cfg | 20 ++++++++++--------- 3 files changed, 33 insertions(+), 27 deletions(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg index aaf85b86e..ba7933667 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg @@ -95,17 +95,19 @@ BeamRemnants:saturation 5 431:addChannel = 1 0.15500 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) ## Λc decays -### Λc+ -> p K- π+ (50%, BR set up to 50%) -4122:oneChannel = 1 0.06472 0 2212 -321 211 ### Λc+ -> p K- π+ (non-resonant) 3.5% -4122:addChannel = 1 0.03624 100 2212 -313 ### Λc+ -> p K*0(892) 1.96% -4122:addChannel = 1 0.01997 100 2224 -321 ### Λc+ -> Delta++ K- 1.08% -4122:addChannel = 1 0.00407 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3 +### Λ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% +4122:addChannel = 1 0.04500 100 2224 -321 ### Λc+ -> Delta++ K- 1.08% +4122:addChannel = 1 0.09000 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3 +### Λc+ -> p K0S (36%) +4122:addChannel = 1 0.36000 0 2212 311 ### Λc+ -> p K0S 1.59% ### Λc+ -> p K- π+ π0 (small, 3%) 4122:addChannel = 1 0.03000 0 2212 -321 211 111 ### Λc+ -> p K- π+ π0 (non-resonant) 4.6% -### Λc+ -> p π- π+ (23.50%) -4122:addChannel = 1 0.23500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59% -### Λc+ -> p K- K+ (23.50%) -4122:addChannel = 1 0.23500 0 2212 333 ### Λc+ -> p phi 1.06% +### Λc+ -> p π- π+ (12.50%) +4122:addChannel = 1 0.12500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59% +### Λc+ -> p K- K+ (12.50%) +4122:addChannel = 1 0.12500 0 2212 333 ### Λc+ -> p phi 1.06% ## Xic decays ### Ξc+ -> p K- π+ (50%) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg index 24058f40d..8740e891f 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg @@ -99,17 +99,19 @@ BeamRemnants:saturation 5 431:addChannel = 1 0.15500 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) ## Λc decays -### Λc+ -> p K- π+ (50%, BR set up to 50%) -4122:oneChannel = 1 0.06472 0 2212 -321 211 ### Λc+ -> p K- π+ (non-resonant) 3.5% -4122:addChannel = 1 0.03624 100 2212 -313 ### Λc+ -> p K*0(892) 1.96% -4122:addChannel = 1 0.01997 100 2224 -321 ### Λc+ -> Delta++ K- 1.08% -4122:addChannel = 1 0.00407 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3 +### Λ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% +4122:addChannel = 1 0.04500 100 2224 -321 ### Λc+ -> Delta++ K- 1.08% +4122:addChannel = 1 0.09000 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3 +### Λc+ -> p K0S (36%) +4122:addChannel = 1 0.36000 0 2212 311 ### Λc+ -> p K0S 1.59% ### Λc+ -> p K- π+ π0 (small, 3%) 4122:addChannel = 1 0.03000 0 2212 -321 211 111 ### Λc+ -> p K- π+ π0 (non-resonant) 4.6% -### Λc+ -> p π- π+ (23.50%) -4122:addChannel = 1 0.23500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59% -### Λc+ -> p K- K+ (23.50%) -4122:addChannel = 1 0.23500 0 2212 333 ### Λc+ -> p phi 1.06% +### Λc+ -> p π- π+ (12.50%) +4122:addChannel = 1 0.12500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59% +### Λc+ -> p K- K+ (12.50%) +4122:addChannel = 1 0.12500 0 2212 333 ### Λc+ -> p phi 1.06% ## Xic decays ### Ξc+ -> p K- π+ (50%) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index 330fc840d..321e77102 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -95,17 +95,19 @@ BeamRemnants:saturation 5 431:addChannel = 1 0.15500 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) ## Λc decays -### Λc+ -> p K- π+ (50%, BR set up to 50%) -4122:oneChannel = 1 0.06472 0 2212 -321 211 ### Λc+ -> p K- π+ (non-resonant) 3.5% -4122:addChannel = 1 0.03624 100 2212 -313 ### Λc+ -> p K*0(892) 1.96% -4122:addChannel = 1 0.01997 100 2224 -321 ### Λc+ -> Delta++ K- 1.08% -4122:addChannel = 1 0.00407 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3 +### Λ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% +4122:addChannel = 1 0.04500 100 2224 -321 ### Λc+ -> Delta++ K- 1.08% +4122:addChannel = 1 0.09000 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3 +### Λc+ -> p K0S (36%) +4122:addChannel = 1 0.36000 0 2212 311 ### Λc+ -> p K0S 1.59% ### Λc+ -> p K- π+ π0 (small, 3%) 4122:addChannel = 1 0.03000 0 2212 -321 211 111 ### Λc+ -> p K- π+ π0 (non-resonant) 4.6% -### Λc+ -> p π- π+ (23.50%) -4122:addChannel = 1 0.23500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59% -### Λc+ -> p K- K+ (23.50%) -4122:addChannel = 1 0.23500 0 2212 333 ### Λc+ -> p phi 1.06% +### Λc+ -> p π- π+ (12.50%) +4122:addChannel = 1 0.12500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59% +### Λc+ -> p K- K+ (12.50%) +4122:addChannel = 1 0.12500 0 2212 333 ### Λc+ -> p phi 1.06% ## Xic decays ### Ξc+ -> p K- π+ (50%) From 3befe66ffdbca05b1d7c5fdbb232854273b4ae73 Mon Sep 17 00:00:00 2001 From: wuctlby Date: Fri, 13 Jun 2025 20:28:57 +0200 Subject: [PATCH 03/11] Enlarge some ratio of the mian analysed channel --- ...ythia8_charmhadronic_with_decays_Mode2.cfg | 44 +++++++++---------- ...adronic_with_decays_Mode2_hardQCD_5TeV.cfg | 44 +++++++++---------- ...charmhadronic_with_decays_Mode2_pp_ref.cfg | 44 +++++++++---------- 3 files changed, 66 insertions(+), 66 deletions(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg index ba7933667..bf5c676d4 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg @@ -66,20 +66,20 @@ BeamRemnants:saturation 5 411:oneChannel = 1 0.50000 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% -### D+ -> K- K+ π+ (23.50%) -411:addChannel = 1 0.06053 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 -411:addChannel = 1 0.04415 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 -411:addChannel = 1 0.06552 0 333 211 ### D+ -> φ π+ 2.69e-3 -411:addChannel = 1 0.06480 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 -### D+ -> π- π+ π+ (23.50%) -411:addChannel = 1 0.06221 0 113 211 ### D+ -> rho0 π+ 8.4e-4 -411:addChannel = 1 0.03456 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 -411:addChannel = 1 0.13824 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 +### D+ -> K- K+ π+ (25.00%) +411:addChannel = 1 0.06439 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 +411:addChannel = 1 0.04697 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 +411:addChannel = 1 0.06970 0 333 211 ### D+ -> φ π+ 2.69e-3 +411:addChannel = 1 0.06894 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 +### D+ -> π- π+ π+ (22.00%) +411:addChannel = 1 0.05824 0 113 211 ### D+ -> rho0 π+ 8.4e-4 +411:addChannel = 1 0.03235 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 +411:addChannel = 1 0.12941 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 ### Ds+ decays ### Ds+ -> K- K+ π+ (50%) -431:oneChannel = 1 0.25000 0 333 211 ### Ds+ -> φ(1020) π+ 2.21% -431:addChannel = 1 0.25000 0 -313 321 ### Ds+ -> antiK*(892) K+ 2.58% +431:oneChannel = 1 0.35000 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, 3%) 431:addChannel = 1 0.03000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% ### Ds+ -> π- π+ π+ (15.50%) @@ -110,17 +110,17 @@ BeamRemnants:saturation 5 4122:addChannel = 1 0.12500 0 2212 333 ### Λc+ -> p phi 1.06% ## Xic decays -### Ξc+ -> p K- π+ (50%) -4232:oneChannel = 1 0.25 0 2212 -321 211 ### Ξc+ -> p K- π+ 6.18e-3 -4232:addChannel = 1 0.25 0 2212 -313 ### Ξc+ -> p antiK*0(892) -### Ξc+ -> p φ (12.5%) -4232:addChannel = 1 0.125 0 2212 333 ### Ξc+ -> p φ -### Ξc+ -> sigma+ π+ π- (12.5%) -4232:addChannel = 1 0.125 0 3222 -211 211 ### Ξc+ -> sigma+ π- π+ 1.37% -### Ξc+ -> Ξ*0 π+ (12.5%) -4232:addChannel = 1 0.125 0 3324 211 -### Ξc+ -> Ξ- π+ π+ (12.5%) -4232:addChannel = 1 0.125 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% +### Ξ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 φ (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% +### Ξc+ -> Ξ*0 π+ (10%) +4232:addChannel = 1 0.12500 0 3324 211 ### add Xic0 decays absent in PYTHIA8 decay table 4132:addChannel = 1 0.0143 0 3312 211 diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg index 8740e891f..b6fd3b837 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg @@ -70,20 +70,20 @@ BeamRemnants:saturation 5 411:oneChannel = 1 0.50000 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% -### D+ -> K- K+ π+ (23.50%) -411:addChannel = 1 0.06053 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 -411:addChannel = 1 0.04415 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 -411:addChannel = 1 0.06552 0 333 211 ### D+ -> φ π+ 2.69e-3 -411:addChannel = 1 0.06480 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 -### D+ -> π- π+ π+ (23.50%) -411:addChannel = 1 0.06221 0 113 211 ### D+ -> rho0 π+ 8.4e-4 -411:addChannel = 1 0.03456 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 -411:addChannel = 1 0.13824 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 +### D+ -> K- K+ π+ (25.00%) +411:addChannel = 1 0.06439 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 +411:addChannel = 1 0.04697 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 +411:addChannel = 1 0.06970 0 333 211 ### D+ -> φ π+ 2.69e-3 +411:addChannel = 1 0.06894 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 +### D+ -> π- π+ π+ (22.00%) +411:addChannel = 1 0.05824 0 113 211 ### D+ -> rho0 π+ 8.4e-4 +411:addChannel = 1 0.03235 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 +411:addChannel = 1 0.12941 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 ### Ds+ decays ### Ds+ -> K- K+ π+ (50%) -431:oneChannel = 1 0.25000 0 333 211 ### Ds+ -> φ(1020) π+ 2.21% -431:addChannel = 1 0.25000 0 -313 321 ### Ds+ -> antiK*(892) K+ 2.58% +431:oneChannel = 1 0.35000 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, 3%) 431:addChannel = 1 0.03000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% ### Ds+ -> π- π+ π+ (15.50%) @@ -114,17 +114,17 @@ BeamRemnants:saturation 5 4122:addChannel = 1 0.12500 0 2212 333 ### Λc+ -> p phi 1.06% ## Xic decays -### Ξc+ -> p K- π+ (50%) -4232:oneChannel = 1 0.25 0 2212 -321 211 ### Ξc+ -> p K- π+ 6.18e-3 -4232:addChannel = 1 0.25 0 2212 -313 ### Ξc+ -> p antiK*0(892) -### Ξc+ -> p φ (12.5%) -4232:addChannel = 1 0.125 0 2212 333 ### Ξc+ -> p φ -### Ξc+ -> sigma+ π+ π- (12.5%) -4232:addChannel = 1 0.125 0 3222 -211 211 ### Ξc+ -> sigma+ π- π+ 1.37% -### Ξc+ -> Ξ*0 π+ (12.5%) -4232:addChannel = 1 0.125 0 3324 211 -### Ξc+ -> Ξ- π+ π+ (12.5%) -4232:addChannel = 1 0.125 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% +### Ξ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 φ (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% +### Ξc+ -> Ξ*0 π+ (10%) +4232:addChannel = 1 0.12500 0 3324 211 ### add Xic0 decays absent in PYTHIA8 decay table 4132:addChannel = 1 0.0143 0 3312 211 diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index 321e77102..e0755ff07 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -66,20 +66,20 @@ BeamRemnants:saturation 5 411:oneChannel = 1 0.50000 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% -### D+ -> K- K+ π+ (23.50%) -411:addChannel = 1 0.06053 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 -411:addChannel = 1 0.04415 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 -411:addChannel = 1 0.06552 0 333 211 ### D+ -> φ π+ 2.69e-3 -411:addChannel = 1 0.06480 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 -### D+ -> π- π+ π+ (23.50%) -411:addChannel = 1 0.06221 0 113 211 ### D+ -> rho0 π+ 8.4e-4 -411:addChannel = 1 0.03456 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 -411:addChannel = 1 0.13824 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 +### D+ -> K- K+ π+ (25.00%) +411:addChannel = 1 0.06439 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 +411:addChannel = 1 0.04697 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 +411:addChannel = 1 0.06970 0 333 211 ### D+ -> φ π+ 2.69e-3 +411:addChannel = 1 0.06894 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 +### D+ -> π- π+ π+ (22.00%) +411:addChannel = 1 0.05824 0 113 211 ### D+ -> rho0 π+ 8.4e-4 +411:addChannel = 1 0.03235 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 +411:addChannel = 1 0.12941 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 ### Ds+ decays ### Ds+ -> K- K+ π+ (50%) -431:oneChannel = 1 0.25000 0 333 211 ### Ds+ -> φ(1020) π+ 2.21% -431:addChannel = 1 0.25000 0 -313 321 ### Ds+ -> antiK*(892) K+ 2.58% +431:oneChannel = 1 0.35000 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, 3%) 431:addChannel = 1 0.03000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% ### Ds+ -> π- π+ π+ (15.50%) @@ -110,17 +110,17 @@ BeamRemnants:saturation 5 4122:addChannel = 1 0.12500 0 2212 333 ### Λc+ -> p phi 1.06% ## Xic decays -### Ξc+ -> p K- π+ (50%) -4232:oneChannel = 1 0.25 0 2212 -321 211 ### Ξc+ -> p K- π+ 6.18e-3 -4232:addChannel = 1 0.25 0 2212 -313 ### Ξc+ -> p antiK*0(892) -### Ξc+ -> p φ (12.5%) -4232:addChannel = 1 0.125 0 2212 333 ### Ξc+ -> p φ -### Ξc+ -> sigma+ π+ π- (12.5%) -4232:addChannel = 1 0.125 0 3222 -211 211 ### Ξc+ -> sigma+ π- π+ 1.37% -### Ξc+ -> Ξ*0 π+ (12.5%) -4232:addChannel = 1 0.125 0 3324 211 -### Ξc+ -> Ξ- π+ π+ (12.5%) -4232:addChannel = 1 0.125 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86% +### Ξ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 φ (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% +### Ξc+ -> Ξ*0 π+ (10%) +4232:addChannel = 1 0.12500 0 3324 211 ### add Xic0 decays absent in PYTHIA8 decay table 4132:addChannel = 1 0.0143 0 3312 211 From ffc4e7a3b25a75bbaa646a3ca029c39e1f6f308b Mon Sep 17 00:00:00 2001 From: wuctlby Date: Fri, 13 Jun 2025 20:35:19 +0200 Subject: [PATCH 04/11] correct the ratio of corr_bkgs Ds --- ...ythia8_charmhadronic_with_decays_Mode2.cfg | 26 +++++++++---------- ...adronic_with_decays_Mode2_hardQCD_5TeV.cfg | 26 +++++++++---------- ...charmhadronic_with_decays_Mode2_pp_ref.cfg | 26 +++++++++---------- 3 files changed, 39 insertions(+), 39 deletions(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg index bf5c676d4..3996ff514 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg @@ -80,19 +80,19 @@ BeamRemnants:saturation 5 ### Ds+ -> K- K+ π+ (50%) 431:oneChannel = 1 0.35000 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, 3%) -431:addChannel = 1 0.03000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% -### Ds+ -> π- π+ π+ (15.50%) -431:addChannel = 1 0.00310 0 113 211 ### Ds+ -> rho π+ 1.1e-4 -431:addChannel = 1 0.00310 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 -431:addChannel = 1 0.14880 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) -### Ds+ -> π- K+ π+ (15.50%) -431:addChannel = 1 0.04340 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 -431:addChannel = 1 0.03100 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 -431:addChannel = 1 0.05580 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 -431:addChannel = 1 0.02480 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 -### Ds+ -> π+ π- π+ π0 (15.50%) -431:addChannel = 1 0.15500 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) +### Ds+ -> K- K+ π+ π0 (small, 2%) +431:addChannel = 1 0.02000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% +### Ds+ -> π- π+ π+ (16.00%) +431:addChannel = 1 0.00320 0 113 211 ### Ds+ -> rho π+ 1.1e-4 +431:addChannel = 1 0.00320 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 +431:addChannel = 1 0.15360 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) +### Ds+ -> π- K+ π+ (16.00%) +431:addChannel = 1 0.04480 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 +431:addChannel = 1 0.03200 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 +431:addChannel = 1 0.05760 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 +431:addChannel = 1 0.02560 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 +### Ds+ -> π+ π- π+ π0 (16.00%) +431:addChannel = 1 0.16000 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) ## Λc decays ### Λc+ -> p K- π+ (36%) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg index b6fd3b837..c4c7c85fc 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg @@ -84,19 +84,19 @@ BeamRemnants:saturation 5 ### Ds+ -> K- K+ π+ (50%) 431:oneChannel = 1 0.35000 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, 3%) -431:addChannel = 1 0.03000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% -### Ds+ -> π- π+ π+ (15.50%) -431:addChannel = 1 0.00310 0 113 211 ### Ds+ -> rho π+ 1.1e-4 -431:addChannel = 1 0.00310 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 -431:addChannel = 1 0.14880 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) -### Ds+ -> π- K+ π+ (15.50%) -431:addChannel = 1 0.04340 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 -431:addChannel = 1 0.03100 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 -431:addChannel = 1 0.05580 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 -431:addChannel = 1 0.02480 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 -### Ds+ -> π+ π- π+ π0 (15.50%) -431:addChannel = 1 0.15500 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) +### Ds+ -> K- K+ π+ π0 (small, 2%) +431:addChannel = 1 0.02000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% +### Ds+ -> π- π+ π+ (16.00%) +431:addChannel = 1 0.00320 0 113 211 ### Ds+ -> rho π+ 1.1e-4 +431:addChannel = 1 0.00320 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 +431:addChannel = 1 0.15360 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) +### Ds+ -> π- K+ π+ (16.00%) +431:addChannel = 1 0.04480 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 +431:addChannel = 1 0.03200 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 +431:addChannel = 1 0.05760 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 +431:addChannel = 1 0.02560 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 +### Ds+ -> π+ π- π+ π0 (16.00%) +431:addChannel = 1 0.16000 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) ## Λc decays ### Λc+ -> p K- π+ (36%) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index e0755ff07..b7d6e0663 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -80,19 +80,19 @@ BeamRemnants:saturation 5 ### Ds+ -> K- K+ π+ (50%) 431:oneChannel = 1 0.35000 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, 3%) -431:addChannel = 1 0.03000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% -### Ds+ -> π- π+ π+ (15.50%) -431:addChannel = 1 0.00310 0 113 211 ### Ds+ -> rho π+ 1.1e-4 -431:addChannel = 1 0.00310 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 -431:addChannel = 1 0.14880 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) -### Ds+ -> π- K+ π+ (15.50%) -431:addChannel = 1 0.04340 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 -431:addChannel = 1 0.03100 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 -431:addChannel = 1 0.05580 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 -431:addChannel = 1 0.02480 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 -### Ds+ -> π+ π- π+ π0 (15.50%) -431:addChannel = 1 0.15500 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) +### Ds+ -> K- K+ π+ π0 (small, 2%) +431:addChannel = 1 0.02000 0 333 213 ### Ds+ -> φ(1020) rho 5.50% +### Ds+ -> π- π+ π+ (16.00%) +431:addChannel = 1 0.00320 0 113 211 ### Ds+ -> rho π+ 1.1e-4 +431:addChannel = 1 0.00320 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 +431:addChannel = 1 0.15360 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) +### Ds+ -> π- K+ π+ (16.00%) +431:addChannel = 1 0.04480 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 +431:addChannel = 1 0.03200 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 +431:addChannel = 1 0.05760 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 +431:addChannel = 1 0.02560 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 +### Ds+ -> π+ π- π+ π0 (16.00%) +431:addChannel = 1 0.16000 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) ## Λc decays ### Λc+ -> p K- π+ (36%) From ae42de08b8c062de6b15a76e67723f9a18a881d5 Mon Sep 17 00:00:00 2001 From: wuctlby Date: Mon, 16 Jun 2025 11:49:20 +0200 Subject: [PATCH 05/11] Add the decay channel D+ -> K*0(892) pi+ --- .../pythia8_charmhadronic_with_decays_Mode2.cfg | 9 +++++---- ...hia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg | 9 +++++---- .../pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg | 9 +++++---- 3 files changed, 15 insertions(+), 12 deletions(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg index 3996ff514..e3e63268b 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg @@ -67,10 +67,11 @@ BeamRemnants:saturation 5 ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% ### D+ -> K- K+ π+ (25.00%) -411:addChannel = 1 0.06439 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 -411:addChannel = 1 0.04697 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 -411:addChannel = 1 0.06970 0 333 211 ### D+ -> φ π+ 2.69e-3 -411:addChannel = 1 0.06894 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 +411:addChannel = 1 0.06243 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 +411:addChannel = 1 0.00764 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 +411:addChannel = 1 0.04553 0 -313 211 ### D+ -> K*0(892) pi+ +411:addChannel = 1 0.06757 0 333 211 ### D+ -> φ π+ 2.69e-3 +411:addChannel = 1 0.06683 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 ### D+ -> π- π+ π+ (22.00%) 411:addChannel = 1 0.05824 0 113 211 ### D+ -> rho0 π+ 8.4e-4 411:addChannel = 1 0.03235 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg index c4c7c85fc..3ce89b7a8 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg @@ -71,10 +71,11 @@ BeamRemnants:saturation 5 ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% ### D+ -> K- K+ π+ (25.00%) -411:addChannel = 1 0.06439 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 -411:addChannel = 1 0.04697 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 -411:addChannel = 1 0.06970 0 333 211 ### D+ -> φ π+ 2.69e-3 -411:addChannel = 1 0.06894 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 +411:addChannel = 1 0.06243 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 +411:addChannel = 1 0.00764 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 +411:addChannel = 1 0.04553 0 -313 211 ### D+ -> K*0(892) pi+ +411:addChannel = 1 0.06757 0 333 211 ### D+ -> φ π+ 2.69e-3 +411:addChannel = 1 0.06683 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 ### D+ -> π- π+ π+ (22.00%) 411:addChannel = 1 0.05824 0 113 211 ### D+ -> rho0 π+ 8.4e-4 411:addChannel = 1 0.03235 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index b7d6e0663..e120f3693 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -67,10 +67,11 @@ BeamRemnants:saturation 5 ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% ### D+ -> K- K+ π+ (25.00%) -411:addChannel = 1 0.06439 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 -411:addChannel = 1 0.04697 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 -411:addChannel = 1 0.06970 0 333 211 ### D+ -> φ π+ 2.69e-3 -411:addChannel = 1 0.06894 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 +411:addChannel = 1 0.06243 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 +411:addChannel = 1 0.00764 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 +411:addChannel = 1 0.04553 0 -313 211 ### D+ -> K*0(892) pi+ +411:addChannel = 1 0.06757 0 333 211 ### D+ -> φ π+ 2.69e-3 +411:addChannel = 1 0.06683 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 ### D+ -> π- π+ π+ (22.00%) 411:addChannel = 1 0.05824 0 113 211 ### D+ -> rho0 π+ 8.4e-4 411:addChannel = 1 0.03235 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 From 119852c8d5f605af63d2a1513a696a2bcadfbc6f Mon Sep 17 00:00:00 2001 From: wuctlby Date: Mon, 16 Jun 2025 12:18:07 +0200 Subject: [PATCH 06/11] =?UTF-8?q?Turn=20on=20the=20match=20of=20=CE=9Bc=20?= =?UTF-8?q?->=20pK0s=20and=20add=20the=20D+=20->=20K*0(1430)=20=CF=80+=201?= =?UTF-8?q?.25%?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .../generator/pythia8_charmhadronic_with_decays_Mode2.cfg | 7 +++++-- ...thia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg | 7 +++++-- .../pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg | 8 ++++++-- 3 files changed, 16 insertions(+), 6 deletions(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg index e3e63268b..f51a4968d 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg @@ -63,7 +63,8 @@ BeamRemnants:saturation 5 ### D+ decays ### D+ -> K- π+ π+ (50%) -411:oneChannel = 1 0.50000 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% +411:oneChannel = 1 0.44120 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% +411:oneChannel = 1 0.05880 0 10311 211 ### D+ -> K*0(1430) π+ 1.25% ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% ### D+ -> K- K+ π+ (25.00%) @@ -246,7 +247,9 @@ BeamRemnants:saturation 5 ### Ds -> eta π+ 431:onIfMatch = 221 211 -### Λc -> p K π +### Λc -> pK0s +4122:onIfMatch = 2212 311 +### Λc -> p K- π+ π0 4122:onIfMatch = 2212 321 211 ### Λc -> p K* 4122:onIfMatch = 2212 313 diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg index 3ce89b7a8..32f040e31 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg @@ -67,7 +67,8 @@ BeamRemnants:saturation 5 ### D+ decays ### D+ -> K- π+ π+ (50%) -411:oneChannel = 1 0.50000 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% +411:oneChannel = 1 0.44120 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% +411:oneChannel = 1 0.05880 0 10311 211 ### D+ -> K*0(1430) π+ 1.25% ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% ### D+ -> K- K+ π+ (25.00%) @@ -250,7 +251,9 @@ BeamRemnants:saturation 5 ### Ds -> eta π+ 431:onIfMatch = 221 211 -### Λc -> p K π +### Λc -> pK0s +4122:onIfMatch = 2212 311 +### Λc -> p K- π+ π0 4122:onIfMatch = 2212 321 211 ### Λc -> p K* 4122:onIfMatch = 2212 313 diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index e120f3693..d01dd4c8a 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -63,7 +63,8 @@ BeamRemnants:saturation 5 ### D+ decays ### D+ -> K- π+ π+ (50%) -411:oneChannel = 1 0.50000 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% +411:oneChannel = 1 0.44120 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% +411:oneChannel = 1 0.05880 0 10311 211 ### D+ -> K*0(1430) π+ 1.25% ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% ### D+ -> K- K+ π+ (25.00%) @@ -125,6 +126,7 @@ BeamRemnants:saturation 5 ### add Xic0 decays absent in PYTHIA8 decay table 4132:addChannel = 1 0.0143 0 3312 211 + ### add OmegaC decays absent in PYTHIA8 decay table 4332:addChannel = 1 0.5 0 3334 211 4332:addChannel = 1 0.5 0 3312 211 @@ -245,7 +247,9 @@ BeamRemnants:saturation 5 ### Ds -> eta π+ 431:onIfMatch = 221 211 -### Λc -> p K π +### Λc -> pK0s +4122:onIfMatch = 2212 311 +### Λc -> p K- π+ π0 4122:onIfMatch = 2212 321 211 ### Λc -> p K* 4122:onIfMatch = 2212 313 From 9708fd8fe81e521d8be168d265ea31b31bd3b020 Mon Sep 17 00:00:00 2001 From: wuctlby Date: Mon, 16 Jun 2025 12:23:30 +0200 Subject: [PATCH 07/11] =?UTF-8?q?Set=2050%=20for=20Ds+=20->=20=CF=86(1020)?= =?UTF-8?q?=20=CF=80+?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- ...ythia8_charmhadronic_with_decays_Mode2.cfg | 26 +++++++++---------- ...adronic_with_decays_Mode2_hardQCD_5TeV.cfg | 26 +++++++++---------- ...charmhadronic_with_decays_Mode2_pp_ref.cfg | 26 +++++++++---------- 3 files changed, 39 insertions(+), 39 deletions(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg index f51a4968d..3f0725175 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg @@ -80,21 +80,21 @@ BeamRemnants:saturation 5 ### Ds+ decays ### Ds+ -> K- K+ π+ (50%) -431:oneChannel = 1 0.35000 0 333 211 ### Ds+ -> φ(1020) π+ 2.21% +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+ -> π- π+ π+ (16.00%) -431:addChannel = 1 0.00320 0 113 211 ### Ds+ -> rho π+ 1.1e-4 -431:addChannel = 1 0.00320 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 -431:addChannel = 1 0.15360 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) -### Ds+ -> π- K+ π+ (16.00%) -431:addChannel = 1 0.04480 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 -431:addChannel = 1 0.03200 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 -431:addChannel = 1 0.05760 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 -431:addChannel = 1 0.02560 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 -### Ds+ -> π+ π- π+ π0 (16.00%) -431:addChannel = 1 0.16000 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) +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+ -> p K- π+ (36%) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg index 32f040e31..965d39e1a 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg @@ -84,21 +84,21 @@ BeamRemnants:saturation 5 ### Ds+ decays ### Ds+ -> K- K+ π+ (50%) -431:oneChannel = 1 0.35000 0 333 211 ### Ds+ -> φ(1020) π+ 2.21% +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+ -> π- π+ π+ (16.00%) -431:addChannel = 1 0.00320 0 113 211 ### Ds+ -> rho π+ 1.1e-4 -431:addChannel = 1 0.00320 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 -431:addChannel = 1 0.15360 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) -### Ds+ -> π- K+ π+ (16.00%) -431:addChannel = 1 0.04480 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 -431:addChannel = 1 0.03200 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 -431:addChannel = 1 0.05760 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 -431:addChannel = 1 0.02560 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 -### Ds+ -> π+ π- π+ π0 (16.00%) -431:addChannel = 1 0.16000 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) +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+ -> p K- π+ (36%) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index d01dd4c8a..9fcd921dc 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -80,21 +80,21 @@ BeamRemnants:saturation 5 ### Ds+ decays ### Ds+ -> K- K+ π+ (50%) -431:oneChannel = 1 0.35000 0 333 211 ### Ds+ -> φ(1020) π+ 2.21% +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+ -> π- π+ π+ (16.00%) -431:addChannel = 1 0.00320 0 113 211 ### Ds+ -> rho π+ 1.1e-4 -431:addChannel = 1 0.00320 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3 -431:addChannel = 1 0.15360 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave) -### Ds+ -> π- K+ π+ (16.00%) -431:addChannel = 1 0.04480 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3 -431:addChannel = 1 0.03200 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3 -431:addChannel = 1 0.05760 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3 -431:addChannel = 1 0.02560 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3 -### Ds+ -> π+ π- π+ π0 (16.00%) -431:addChannel = 1 0.16000 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel) +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+ -> p K- π+ (36%) From 735f013256d230457da38b40079cea9b1f87f081 Mon Sep 17 00:00:00 2001 From: wuctlby Date: Mon, 16 Jun 2025 15:06:40 +0200 Subject: [PATCH 08/11] =?UTF-8?q?correct=20the=20decay=20channel=20D+=20->?= =?UTF-8?q?=20antiK*0(1430)=20=CF=80+=201.25%?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .../generator/pythia8_charmhadronic_with_decays_Mode2.cfg | 4 +++- .../pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg | 4 +++- .../pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg | 2 +- 3 files changed, 7 insertions(+), 3 deletions(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg index 3f0725175..11739ad2d 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg @@ -64,7 +64,7 @@ BeamRemnants:saturation 5 ### D+ decays ### D+ -> K- π+ π+ (50%) 411:oneChannel = 1 0.44120 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% -411:oneChannel = 1 0.05880 0 10311 211 ### D+ -> K*0(1430) π+ 1.25% +411:addChannel = 1 0.05880 0 -10311 211 ### D+ -> antiK*0(1430) π+ 1.25% ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% ### D+ -> K- K+ π+ (25.00%) @@ -212,6 +212,8 @@ BeamRemnants:saturation 5 ### 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 diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg index 965d39e1a..221b62fd8 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg @@ -68,7 +68,7 @@ BeamRemnants:saturation 5 ### D+ decays ### D+ -> K- π+ π+ (50%) 411:oneChannel = 1 0.44120 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% -411:oneChannel = 1 0.05880 0 10311 211 ### D+ -> K*0(1430) π+ 1.25% +411:addChannel = 1 0.05880 0 -10311 211 ### D+ -> antiK*0(1430) π+ 1.25% ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% ### D+ -> K- K+ π+ (25.00%) @@ -216,6 +216,8 @@ BeamRemnants:saturation 5 ### 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 diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index 9fcd921dc..41905463e 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -64,7 +64,7 @@ BeamRemnants:saturation 5 ### D+ decays ### D+ -> K- π+ π+ (50%) 411:oneChannel = 1 0.44120 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% -411:oneChannel = 1 0.05880 0 10311 211 ### D+ -> K*0(1430) π+ 1.25% +411:addChannel = 1 0.05880 0 -10311 211 ### D+ -> antiK*0(1430) π+ 1.25% ### D+ -> K- π+ π+ π0 (small, 3%) 411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25% ### D+ -> K- K+ π+ (25.00%) From 8c2e824e0ba5876df87633c30ad16b07ebeebe02 Mon Sep 17 00:00:00 2001 From: wuctlby Date: Mon, 16 Jun 2025 15:08:13 +0200 Subject: [PATCH 09/11] correct the comment of 321 for K and 211 for pion --- .../pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg | 2 ++ 1 file changed, 2 insertions(+) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index 41905463e..6a2375eff 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -212,6 +212,8 @@ BeamRemnants:saturation 5 ### 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 From d8ae175697b730e95be01869f1108e8a42879909 Mon Sep 17 00:00:00 2001 From: wuctlby Date: Mon, 16 Jun 2025 15:19:14 +0200 Subject: [PATCH 10/11] =?UTF-8?q?add=20onIfMatch=20for=20`D+/-=20->=20K*?= =?UTF-8?q?=20=CF=80`?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .../generator/pythia8_charmhadronic_with_decays_Mode2.cfg | 2 ++ .../pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg | 2 ++ .../pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg | 2 ++ 3 files changed, 6 insertions(+) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg index 11739ad2d..2d0df6d61 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg @@ -209,6 +209,8 @@ BeamRemnants:saturation 5 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* π diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg index 221b62fd8..51279c032 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg @@ -213,6 +213,8 @@ BeamRemnants:saturation 5 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* π diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index 6a2375eff..c3fe0c0c2 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -209,6 +209,8 @@ BeamRemnants:saturation 5 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* π From db13521b4aec616d83fcb84c7694400b9297e1b7 Mon Sep 17 00:00:00 2001 From: wuctlby Date: Mon, 16 Jun 2025 18:19:43 +0200 Subject: [PATCH 11/11] =?UTF-8?q?Set=20the=20`D+=20->=20=CF=86=20=CF=80+`?= =?UTF-8?q?=20to=2025%=20and=20average=20the=20rest,=20move=20the=20`D+=20?= =?UTF-8?q?->=20K*0(892)=20pi+`=20to=20where=20the=20right=20final=20state?= =?UTF-8?q?=20is?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- ...ythia8_charmhadronic_with_decays_Mode2.cfg | 28 +++++++++---------- ...adronic_with_decays_Mode2_hardQCD_5TeV.cfg | 28 +++++++++---------- ...charmhadronic_with_decays_Mode2_pp_ref.cfg | 28 +++++++++---------- 3 files changed, 42 insertions(+), 42 deletions(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg index 2d0df6d61..95402f146 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg @@ -42,7 +42,7 @@ BeamRemnants:saturation 5 4332:tau0 = 0.0803 ### Force golden charm hadrons decay modes for D2H studies -### HF decays +# 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 @@ -63,20 +63,20 @@ BeamRemnants:saturation 5 ### D+ decays ### D+ -> K- π+ π+ (50%) -411:oneChannel = 1 0.44120 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% -411:addChannel = 1 0.05880 0 -10311 211 ### D+ -> antiK*0(1430) π+ 1.25% +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+ π+ (25.00%) -411:addChannel = 1 0.06243 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 -411:addChannel = 1 0.00764 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 -411:addChannel = 1 0.04553 0 -313 211 ### D+ -> K*0(892) pi+ -411:addChannel = 1 0.06757 0 333 211 ### D+ -> φ π+ 2.69e-3 -411:addChannel = 1 0.06683 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 -### D+ -> π- π+ π+ (22.00%) -411:addChannel = 1 0.05824 0 113 211 ### D+ -> rho0 π+ 8.4e-4 -411:addChannel = 1 0.03235 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 -411:addChannel = 1 0.12941 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 +### 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%) @@ -175,7 +175,7 @@ BeamRemnants:saturation 5 3334:onMode = off 3334:onIfAll = 3122 -321 -### switch off all decay channels +# switch off all decay channels 411:onMode = off 421:onMode = off 431:onMode = off diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg index 51279c032..d75c500ac 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg @@ -46,7 +46,7 @@ BeamRemnants:saturation 5 4332:tau0 = 0.0803 ### Force golden charm hadrons decay modes for D2H studies -### HF decays +# 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 @@ -67,20 +67,20 @@ BeamRemnants:saturation 5 ### D+ decays ### D+ -> K- π+ π+ (50%) -411:oneChannel = 1 0.44120 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% -411:addChannel = 1 0.05880 0 -10311 211 ### D+ -> antiK*0(1430) π+ 1.25% +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+ π+ (25.00%) -411:addChannel = 1 0.06243 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 -411:addChannel = 1 0.00764 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 -411:addChannel = 1 0.04553 0 -313 211 ### D+ -> K*0(892) pi+ -411:addChannel = 1 0.06757 0 333 211 ### D+ -> φ π+ 2.69e-3 -411:addChannel = 1 0.06683 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 -### D+ -> π- π+ π+ (22.00%) -411:addChannel = 1 0.05824 0 113 211 ### D+ -> rho0 π+ 8.4e-4 -411:addChannel = 1 0.03235 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 -411:addChannel = 1 0.12941 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 +### 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%) @@ -179,7 +179,7 @@ BeamRemnants:saturation 5 3334:onMode = off 3334:onIfAll = 3122 -321 -### switch off all decay channels +# Switch off all decay channels 411:onMode = off 421:onMode = off 431:onMode = off diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg index c3fe0c0c2..00743a912 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg @@ -42,7 +42,7 @@ BeamRemnants:saturation 5 4332:tau0 = 0.0803 ### Force golden charm hadrons decay modes for D2H studies -### HF decays +# 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 @@ -63,20 +63,20 @@ BeamRemnants:saturation 5 ### D+ decays ### D+ -> K- π+ π+ (50%) -411:oneChannel = 1 0.44120 0 -321 211 211 ### D+ -> K- π+ π+ 9.38% -411:addChannel = 1 0.05880 0 -10311 211 ### D+ -> antiK*0(1430) π+ 1.25% +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+ π+ (25.00%) -411:addChannel = 1 0.06243 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3 -411:addChannel = 1 0.00764 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3 -411:addChannel = 1 0.04553 0 -313 211 ### D+ -> K*0(892) pi+ -411:addChannel = 1 0.06757 0 333 211 ### D+ -> φ π+ 2.69e-3 -411:addChannel = 1 0.06683 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 -### D+ -> π- π+ π+ (22.00%) -411:addChannel = 1 0.05824 0 113 211 ### D+ -> rho0 π+ 8.4e-4 -411:addChannel = 1 0.03235 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4 -411:addChannel = 1 0.12941 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4 +### 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%) @@ -175,7 +175,7 @@ BeamRemnants:saturation 5 3334:onMode = off 3334:onIfAll = 3122 -321 -### switch off all decay channels +# switch off all decay channels 411:onMode = off 421:onMode = off 431:onMode = off