1+ # ## authors: Fabrizio Grosa (fabrizio.grosa@cern.ch)
2+ # ## Cristina Terrevoli (cristina.terrevoli@cern.ch)
3+ # ## Fabio Catalano (fabio.catalano@cern.ch)
4+ # ## last update: November 2023
5+
6+ # ## beams
7+ Beams:idA 2212 # proton
8+ Beams:idB 2212 # proton
9+ Beams:eCM 5360. # GeV
10+
11+ # ## processes
12+ SoftQCD:inelastic on # all inelastic processes
13+
14+ # ## decays
15+ ParticleDecays:limitTau0 on
16+ ParticleDecays:tau0Max 10.
17+
18+ # ## switching on Pythia Mode2
19+ ColourReconnection:mode 1
20+ ColourReconnection:allowDoubleJunRem off
21+ ColourReconnection:m0 0.3
22+ ColourReconnection:allowJunctions on
23+ ColourReconnection:junctionCorrection 1.20
24+ ColourReconnection:timeDilationMode 2
25+ ColourReconnection:timeDilationPar 0.18
26+ StringPT:sigma 0.335
27+ StringZ:aLund 0.36
28+ StringZ:bLund 0.56
29+ StringFlav:probQQtoQ 0.078
30+ StringFlav:ProbStoUD 0.2
31+ StringFlav:probQQ1toQQ0join 0.0275,0.0275,0.0275,0.0275
32+ MultiPartonInteractions:pT0Ref 2.15
33+ BeamRemnants:remnantMode 1
34+ BeamRemnants:saturation 5
35+
36+ # Correct decay lengths (wrong in PYTHIA8 decay table)
37+ # Lb
38+ 5122:tau0 = 0.4390
39+ # Xic0
40+ 4132:tau0 = 0.0455
41+ # OmegaC
42+ 4332:tau0 = 0.0803
43+
44+ # ## Force golden charm hadrons decay modes for D2H studies
45+ # ## add D0 decays absent in PYTHIA8 decay table and set BRs from PDG for other
46+ 421:oneChannel = 1 0.0389 0 -321 211
47+ 421:addChannel = 1 0.00389 0 -321 211 111
48+ # ## add D+ decays absent in PYTHIA8 decay table and set BRs from PDG for other
49+ 411:oneChannel = 1 0.0752 0 -321 211 211
50+ 411:addChannel = 1 0.0104 0 -313 211
51+ 411:addChannel = 1 0.0156 0 311 211
52+ 411:addChannel = 1 0.0752 0 333 211 # to have the same amount of D+->KKpi and D+->Kpipi
53+ # # add Lc decays absent in PYTHIA8 decay table and set BRs from PDG for other
54+ 4122:oneChannel = 1 0.0196 100 2212 -313
55+ 4122:addChannel = 1 0.0108 100 2224 -321
56+ 4122:addChannel = 1 0.022 100 102134 211
57+ 4122:addChannel = 1 0.035 0 2212 -321 211
58+ 4122:addChannel = 1 0.0159 0 2212 311
59+ # ## add Xic+ decays absent in PYTHIA8 decay table
60+ 4232:addChannel = 1 0.2 0 2212 -313
61+ 4232:addChannel = 1 0.2 0 2212 -321 211
62+ 4232:addChannel = 1 0.2 0 3324 211
63+ 4232:addChannel = 1 0.2 0 3312 211 211
64+ # ## add Xic0 decays absent in PYTHIA8 decay table
65+ 4132:addChannel = 1 0.0143 0 3312 211
66+ # ## add OmegaC decays absent in PYTHIA8 decay table
67+ 4332:addChannel = 1 0.5 0 3334 211
68+ 4332:addChannel = 1 0.5 0 3312 211
69+
70+ # ## K* -> K pi
71+ 313:onMode = off
72+ 313:onIfAll = 321 211
73+ # ## for Ds -> Phi pi+
74+ 333:onMode = off
75+ 333:onIfAll = 321 321
76+ # ## for D0 -> rho0 pi+ k-
77+ 113:onMode = off
78+ 113:onIfAll = 211 211
79+ # ## for Lambda_c -> Delta++ K-
80+ 2224:onMode = off
81+ 2224:onIfAll = 2212 211
82+ # ## for Lambda_c -> Lambda(1520) K-
83+ 102134:onMode = off
84+ 102134:onIfAll = 2212 321
85+ # ## for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p
86+ # ## and Omega_c -> pi Xi -> pi pi Lambda -> pi pi pi p
87+ 3312:onMode = off
88+ 3312:onIfAll = 3122 -211
89+ 3122:onMode = off
90+ 3122:onIfAll = 2212 -211
91+ # ## for Omega_c -> pi Omega -> pi K Lambda -> pi K pi p
92+ 3334:onMode = off
93+ 3334:onIfAll = 3122 -321
94+
95+ # ## switch off all decay channels
96+ 411:onMode = off
97+ 421:onMode = off
98+ 431:onMode = off
99+ 4122:onMode = off
100+ 4232:onMode = off
101+ 4132:onMode = off
102+ 443:onMode = off
103+ 4332:onMode = off
104+
105+ # ## D0 -> K pi
106+ 421:onIfMatch = 321 211
107+ # ## D0 -> K pi pi0
108+ 421:onIfMatch = 321 211 111
109+
110+ # ## D+/- -> K pi pi
111+ 411:onIfMatch = 321 211 211
112+ # ## D+/- -> K* pi
113+ 411:onIfMatch = 313 211
114+ # ## D+/- -> phi pi
115+ 411:onIfMatch = 333 211
116+
117+ # ## D_s -> K K*
118+ 431:onIfMatch = 321 313
119+ # ## D_s -> Phi pi
120+ 431:onIfMatch = 333 211
121+
122+ # ## Lambda_c -> p K*
123+ 4122:onIfMatch = 2212 313
124+ # ## Lambda_c -> Delta K
125+ 4122:onIfMatch = 2224 321
126+ # ## Lambda_c -> Lambda(1520) pi
127+ 4122:onIfMatch = 102134 211
128+ # ## Lambda_c -> p K pi
129+ 4122:onIfMatch = 2212 321 211
130+ # ## Lambda_c -> pK0s
131+ 4122:onIfMatch = 2212 311
132+
133+ # ## Xic+ -> pK*0
134+ 4232:onIfMatch = 2212 313
135+ # ## Xic+ -> p K- pi+
136+ 4232:onIfMatch = 2212 321 211
137+ # ## Xic+ -> Xi*0 pi+, Xi*->Xi- pi+
138+ 4232:onIfMatch = 3324 211
139+ # ## Xic+ -> Xi- pi+ pi+
140+ 4232:onIfMatch = 3312 211 211
141+
142+ # ## Xic0 -> Xi- pi+
143+ 4132:onIfMatch = 3312 211
144+
145+ # ## Omega_c -> Omega pi
146+ 4332:onIfMatch = 3334 211
147+ # ## Omega_c -> Xi pi
148+ 4332:onIfMatch = 3312 211
0 commit comments