Skip to content

Commit 64670ba

Browse files
wuctlbywuctlby
committed
Add pp reference with correlated bkg
1 parent 5f56534 commit 64670ba

File tree

1 file changed

+244
-0
lines changed

1 file changed

+244
-0
lines changed

‎MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref_CorrBkg.cfg‎

Lines changed: 244 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,244 @@
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+
### HF decays
45+
### BR are set to yield 50% of signal in golden channel and uniform abundance of corr. bkg channels (weighted by BR from PDG)
46+
###
47+
### D0 decays
48+
### D0 -> K- π+ (50%)
49+
421:oneChannel = 1 0.50000 0 -321 211 ### D0 -> K- π+ 3.94%
50+
### D0 -> K- π+ π0 (12.50%)
51+
421:addChannel = 1 0.00865 0 -321 211 111 ### D0 -> K- π+ π0 (non-resonant) 1.15% (e.g. 115/(115+231+195+1120)*0.2)
52+
421:addChannel = 1 0.08429 0 213 -321 ### D0 -> rho+ K- 11.2%
53+
421:addChannel = 1 0.01467 0 -313 111 ### D0 -> antiK*0(892) π0 1.95%
54+
421:addChannel = 1 0.01739 0 -323 211 ### D0 -> K*-(892) π+ 2.31%
55+
### D0 -> π- π+ (12.50%)
56+
421:addChannel = 1 0.12500 0 -211 211 ### D0 -> π- π+ (non-resonant) 1.0e-4
57+
### D0 -> π- π+ π0 (12.50%)
58+
421:addChannel = 1 0.12341 0 213 -211 ### D0 -> rho+ π- 1.01%
59+
421:addChannel = 1 0.00159 0 -211 211 111 ### D0 -> π- π+ π0 (non-resonant) 1.3e-4
60+
### D0 -> K- K+ (12.50%)
61+
421:addChannel = 1 0.12500 0 -321 321 ### D0 -> K- K+ (non-resonant) 4.08e-3
62+
63+
### D+ decays
64+
### D+ -> K- π+ π+ (50%)
65+
411:oneChannel = 1 0.50000 0 -321 211 211 ### D+ -> K- π+ π+ 9.38%
66+
### D+ -> K- π+ π+ π0 (small, 3%)
67+
411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25%
68+
### D+ -> K- K+ π+ (23.5%)
69+
411:addChannel = 1 0.06045 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3
70+
411:addChannel = 1 0.04418 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3
71+
411:addChannel = 1 0.06530 0 333 211 ### D+ -> φ π+ 2.69e-3
72+
411:addChannel = 1 0.06506 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3
73+
### D+ -> π- π+ π+ (23.5%)
74+
411:addChannel = 1 0.14100 0 113 211 ### D+ -> rho0 π+ 8.4e-4
75+
411:addChannel = 1 0.07721 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4
76+
411:addChannel = 1 0.01679 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4
77+
78+
### Ds+ decays
79+
### Ds+ -> K- K+ π+ (50%)
80+
431:oneChannel = 1 0.24776 0 333 211 ### Ds+ -> φ(1020) π+ 2.21%
81+
431:addChannel = 1 0.25224 0 -313 321 ### Ds+ -> antiK*(892) K+ 2.58%
82+
### Ds+ -> K- K+ π+ π0 (small, 2%)
83+
431:addChannel = 1 0.02000 0 333 213 ### Ds+ -> φ(1020) rho 5.50%
84+
### Ds+ -> π- π+ π+ (16.00%)
85+
431:addChannel = 1 0.00188 0 113 211 ### Ds+ -> rho π+ 1.1e-4
86+
431:addChannel = 1 0.00239 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3
87+
431:addChannel = 1 0.15573 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave)
88+
### Ds+ -> π- K+ π+ (16.00%)
89+
431:addChannel = 1 0.04310 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3
90+
431:addChannel = 1 0.03097 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3
91+
431:addChannel = 1 0.05600 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3
92+
431:addChannel = 1 0.02993 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3
93+
### Ds+ -> π+ π- π+ π0 (16.00%)
94+
431:addChannel = 1 0.16000 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel)
95+
96+
## Λc decays
97+
### Λc+ -> p K- π+ (50%, BR set up to 50%)
98+
4122:oneChannel = 1 0.06472 0 2212 -321 211 ### Λc+ -> p K- π+ (non-resonant) 3.5%
99+
4122:addChannel = 1 0.03624 100 2212 -313 ### Λc+ -> p K*0(892) 1.96%
100+
4122:addChannel = 1 0.01997 100 2224 -321 ### Λc+ -> Delta++ K- 1.08%
101+
4122:addChannel = 1 0.00407 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3
102+
### Λc+ -> p K- π+ π0 (small, 5%)
103+
4122:addChannel = 1 0.05000 0 2212 -321 211 111 ### Λc+ -> p K- π+ π0 (non-resonant) 4.6%
104+
### Λc+ -> p π- π+ (23.50%)
105+
4122:addChannel = 1 0.23500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59%
106+
### Λc+ -> p K- K+ (23.50%)
107+
4122:addChannel = 1 0.23500 0 2212 333 ### Λc+ -> p phi 1.06%
108+
109+
### Ξc+ decays (uniformly distributed)
110+
### Ξc+ -> p K- π+ (34%)
111+
4232:oneChannel = 1 0.17 0 2212 -321 211 ### Ξc+ -> p K- π+ 6.2e-3
112+
4232:addChannel = 1 0.17 0 2212 -313 ### Ξc+ -> p antiK*0(892) 3.3e-3
113+
### Ξc+ -> p K- K+ (33%)
114+
4232:addChannel = 1 0.33 0 2212 333 ### Ξc+ -> p φ
115+
### Ξc+ -> sigma+ π+ π- (33%)
116+
4232:addChannel = 1 0.33 0 3222 -211 211 ### Ξc+ -> sigma+ π- π+
117+
118+
# Allow the decay of resonances in the decay chain
119+
### K*0(892) -> K- π+
120+
313:onMode = off
121+
313:onIfAll = 321 211
122+
### K*(892)+ -> K- π0
123+
323:onMode = off
124+
323:onIfAll = 321 111
125+
### K*(1430)0 -> K- π+
126+
10311:onMode = off
127+
10311:onIfAll = 321 211
128+
### rho+ -> π+ π0
129+
213:onMode = off
130+
213:onIfAll = 211 111
131+
### φ -> K+ K-
132+
333:onMode = off
133+
333:onIfAll = 321 321
134+
### rho0 -> π+ π-
135+
113:onMode = off
136+
113:onIfAll = 211 211
137+
### f2(1270) -> π+ π-
138+
225:onMode = off
139+
225:onIfAll = 211 211
140+
### f0(1370) -> π+ π-
141+
10221:onMode = off
142+
10221:onIfAll = 211 211
143+
### eta -> π+ π-
144+
221:onMode = off
145+
221:onIfAll = 111 211 211
146+
### for Λc -> Delta++ K-
147+
2224:onMode = off
148+
2224:onIfAll = 2212 211
149+
### for Λc -> Lambda(1520) K-
150+
102134:onMode = off
151+
102134:onIfAll = 2212 321
152+
153+
### Switch off all decay channels
154+
411:onMode = off
155+
421:onMode = off
156+
431:onMode = off
157+
4122:onMode = off
158+
4232:onMode = off
159+
160+
# Allow the decay of HF
161+
### D0 -> K π
162+
421:onIfMatch = 321 211
163+
### D0 -> K π π0
164+
421:onIfMatch = 321 211 111
165+
### D0 -> rho K
166+
421:onIfMatch = 213 321
167+
### D0 -> antiK*0 π0
168+
421:onIfMatch = 313 111
169+
### D0 -> K*- π+
170+
421:onIfMatch = 323 211
171+
### D0 -> π π
172+
421:onIfMatch = 211 211
173+
### D0 -> rho+ π-
174+
421:onIfMatch = 213 -211
175+
### D0 -> π π π0
176+
421:onIfMatch = 211 211 111
177+
### D0 -> K K
178+
421:onIfMatch = 321 321
179+
180+
### D+/- -> K π π
181+
411:onIfMatch = 321 211 211
182+
### D+/- -> K π π
183+
411:onIfMatch = 321 211 211 111
184+
### D+/- -> K* K
185+
411:onIfMatch = 313 321
186+
### D+/- -> antiK* π
187+
411:onIfMatch = 10311 321
188+
### D+/- -> φ π
189+
411:onIfMatch = 333 211
190+
### D+/- -> K K π
191+
411:onIfMatch = 321 321 211
192+
### D+/- -> f2(1270) π
193+
411:onIfMatch = 225 211
194+
### D+/- -> rho π
195+
411:onIfMatch = 113 211
196+
### D+/- -> π π π
197+
411:onIfMatch = 211 211 211
198+
199+
### Ds -> φ π
200+
431:onIfMatch = 333 211
201+
### Ds -> K* K
202+
431:onIfMatch = 313 321
203+
### Ds -> φ rho+
204+
431:onIfMatch = 333 213
205+
### Ds -> rho π
206+
431:onIfMatch = 113 211
207+
### Ds -> f2 π
208+
431:onIfMatch = 225 211
209+
### Ds -> π π π
210+
431:onIfMatch = 211 211 211
211+
### Ds -> K*(892)0 π+
212+
431:onIfMatch = 313 211
213+
### Ds -> f0(1370) K+
214+
431:onIfMatch = 10221 321
215+
### Ds -> rho0 K+
216+
431:onIfMatch = 113 321
217+
### Ds -> K π π
218+
431:onIfMatch = 321 211 211
219+
### Ds -> eta π+
220+
431:onIfMatch = 221 211
221+
222+
### Λc -> p K π
223+
4122:onIfMatch = 2212 321 211
224+
### Λc -> p K*
225+
4122:onIfMatch = 2212 313
226+
### Λc -> Delta++ K
227+
4122:onIfMatch = 2224 321
228+
### Λc -> Lambda(1520) π
229+
4122:onIfMatch = 102134 211
230+
### Λc -> p K- π+ π0
231+
4122:onIfMatch = 2212 321 211 111
232+
### Λc -> p π π
233+
4122:onIfMatch = 2212 211 211
234+
### Λc -> p K K
235+
4122:onIfMatch = 2212 333
236+
237+
### Ξc+ -> p K- π+
238+
4232:onIfMatch = 2212 321 211
239+
### Ξc+ -> p antiK*0(892)
240+
4232:onIfMatch = 2212 313
241+
### Ξc+ -> p φ
242+
4232:onIfMatch = 2212 333
243+
### Ξc+ -> sigma- π+ π+
244+
4232:onIfMatch = 3222 211 211

0 commit comments

Comments
 (0)