Skip to content

Commit 6013121

Browse files
wuctlbywuctlby
authored
[PWGHF] Add the correlated bkg channel for PbPb@5.36TeV, pp@13.6TeV, and pp reference @5.36TeV (#2030)
* Add the correlated bkg channel for PbPb@5.36TeV, pp@13.6TeV, and ppref@5.36TeV * correct the abundance ratio of Lc+, and add the Lambdac+ -> p k0s * Enlarge some ratio of the mian analysed channel * correct the ratio of corr_bkgs Ds * Add the decay channel D+ -> K*0(892) pi+ * Turn on the match of Λc -> pK0s and add the D+ -> K*0(1430) π+ 1.25% * Set 50% for Ds+ -> φ(1020) π+ * correct the decay channel D+ -> antiK*0(1430) π+ 1.25% * correct the comment of 321 for K and 211 for pion * add onIfMatch for `D+/- -> K* π` * Set the `D+ -> φ π+` to 25% and average the rest, move the `D+ -> K*0(892) pi+` to where the right final state is
1 parent dbe3a8d commit 6013121

File tree

3 files changed

+563
-137
lines changed

3 files changed

+563
-137
lines changed

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

Lines changed: 188 additions & 46 deletions
Original file line numberDiff line numberDiff line change
@@ -42,44 +42,127 @@ BeamRemnants:saturation 5
4242
4332:tau0 = 0.0803
4343

4444
### 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
45+
# HF decays
46+
### BR are set to yield 50% of signal in golden channel and uniform abundance of corr. bkg channels (weighted by BR from PDG)
47+
###
48+
### D0 decays
49+
### D0 -> K- π+ (50%)
50+
421:oneChannel = 1 0.50000 0 -321 211 ### D0 -> K- π+ 3.94%
51+
### D0 -> K- π+ π0 (12.50%)
52+
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)
53+
421:addChannel = 1 0.08750 0 213 -321 ### D0 -> rho+ K- 11.2%
54+
421:addChannel = 1 0.01250 0 -313 111 ### D0 -> antiK*0(892) π0 1.95%
55+
421:addChannel = 1 0.01875 0 -323 211 ### D0 -> K*-(892) π+ 2.31%
56+
### D0 -> π- π+ (12.50%)
57+
421:addChannel = 1 0.12500 0 -211 211 ### D0 -> π- π+ (non-resonant) 1.0e-4
58+
### D0 -> π- π+ π0 (12.50%)
59+
421:addChannel = 1 0.08750 0 213 -211 ### D0 -> rho+ π- 1.01%
60+
421:addChannel = 1 0.03750 0 -211 211 111 ### D0 -> π- π+ π0 (non-resonant) 1.3e-4
61+
### D0 -> K- K+ (12.50%)
62+
421:addChannel = 1 0.12500 0 -321 321 ### D0 -> K- K+ (non-resonant) 4.08e-3
63+
64+
### D+ decays
65+
### D+ -> K- π+ π+ (50%)
66+
411:oneChannel = 1 0.40189 0 -321 211 211 ### D+ -> K- π+ π+ 9.38%
67+
411:addChannel = 1 0.05356 0 -10311 211 ### D+ -> antiK*0(1430) π+ 1.25%
68+
411:addChannel = 1 0.04455 0 -313 211 ### D+ -> K*0(892) pi+ 1.04%
69+
### D+ -> K- π+ π+ π0 (small, 3%)
70+
411:addChannel = 1 0.03000 0 -321 211 211 111 ### D+ -> K- π+ π+ π0 6.25%
71+
### D+ -> K- K+ π+ (36.00%, set 25% for D+ -> φ π+, 11% for the rest)
72+
411:addChannel = 1 0.25000 0 333 211 ### D+ -> φ π+ 2.69e-3 !needed for signal
73+
411:addChannel = 1 0.03929 0 -313 321 ### D+ -> K*0(892) K+ 2.49e-3
74+
411:addChannel = 1 0.02865 0 -10311 321 ### D+ -> antiK*0(1430) K+ 1.82e-3
75+
411:addChannel = 1 0.04206 0 -321 321 211 ### D+ -> K- K+ π+ (non-resonant) 2.68e-3
76+
### D+ -> π- π+ π+ (11.00%)
77+
411:addChannel = 1 0.02911 0 113 211 ### D+ -> rho0 π+ 8.4e-4
78+
411:addChannel = 1 0.01618 0 225 211 ### D+ -> f2(1270) π+ 4.6e-4
79+
411:addChannel = 1 0.06471 0 -211 211 211 ### D+ -> π- π+ π+ (non-resonant) 1.0e-4
80+
81+
### Ds+ decays
82+
### Ds+ -> K- K+ π+ (50%)
83+
431:oneChannel = 1 0.50000 0 333 211 ### Ds+ -> φ(1020) π+ 2.21%
84+
431:addChannel = 1 0.15000 0 -313 321 ### Ds+ -> antiK*(892) K+ 2.58%
85+
### Ds+ -> K- K+ π+ π0 (small, 2%)
86+
431:addChannel = 1 0.02000 0 333 213 ### Ds+ -> φ(1020) rho 5.50%
87+
### Ds+ -> π- π+ π+ (11.00%)
88+
431:addChannel = 1 0.00220 0 113 211 ### Ds+ -> rho π+ 1.1e-4
89+
431:addChannel = 1 0.00220 0 225 211 ### Ds+ -> f2(1270) π+ 1.4e-3
90+
431:addChannel = 1 0.10560 0 -211 211 211 ### Ds+ -> π- π+ π+ 9.12e-3 (s-wave)
91+
### Ds+ -> π- K+ π+ (11.00%)
92+
431:addChannel = 1 0.03080 0 313 211 ### Ds+ -> K*(892)0 π+ 1.67e-3
93+
431:addChannel = 1 0.02200 0 10221 321 ### Ds+ -> f0(1370) K+ 1.2e-3
94+
431:addChannel = 1 0.03960 0 113 321 ### Ds+ -> rho0 K+ 2.17e-3
95+
431:addChannel = 1 0.01760 0 -211 321 211 ### Ds+ -> π- K+ π+ (non-resonant) 1.16-3
96+
### Ds+ -> π+ π- π+ π0 (11.00%)
97+
431:addChannel = 1 0.11000 0 221 211 ### Ds+ -> eta π+ -> π0 π+ π+ π- (affects D+ golden channel)
98+
99+
## Λc decays
100+
### Λc+ -> p K- π+ (36%)
101+
4122:oneChannel = 1 0.14400 0 2212 -321 211 ### Λc+ -> p K- π+ (non-resonant) 3.5%
102+
4122:addChannel = 1 0.08100 100 2212 -313 ### Λc+ -> p K*0(892) 1.96%
103+
4122:addChannel = 1 0.04500 100 2224 -321 ### Λc+ -> Delta++ K- 1.08%
104+
4122:addChannel = 1 0.09000 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3
105+
### Λc+ -> p K0S (36%)
106+
4122:addChannel = 1 0.36000 0 2212 311 ### Λc+ -> p K0S 1.59%
107+
### Λc+ -> p K- π+ π0 (small, 3%)
108+
4122:addChannel = 1 0.03000 0 2212 -321 211 111 ### Λc+ -> p K- π+ π0 (non-resonant) 4.6%
109+
### Λc+ -> p π- π+ (12.50%)
110+
4122:addChannel = 1 0.12500 0 2212 -211 211 ### Λc+ -> p π+ π+ 4.59%
111+
### Λc+ -> p K- K+ (12.50%)
112+
4122:addChannel = 1 0.12500 0 2212 333 ### Λc+ -> p phi 1.06%
113+
114+
## Xic decays
115+
### Ξc+ -> p K- π+ (35%)
116+
4232:oneChannel = 1 0.17500 0 2212 -321 211 ### Ξc+ -> p K- π+ 6.18e-3
117+
4232:addChannel = 1 0.17500 0 2212 -313 ### Ξc+ -> p antiK*0(892)
118+
### Ξc+ -> Ξ- π+ π+ (35%) (set the same as Ξc+ -> p K- π+)
119+
4232:addChannel = 1 0.35000 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86%
120+
### Ξc+ -> p φ (10%)
121+
4232:addChannel = 1 0.10000 0 2212 333 ### Ξc+ -> p φ
122+
### Ξc+ -> sigma+ π+ π- (10%)
123+
4232:addChannel = 1 0.12500 0 3222 -211 211 ### Ξc+ -> sigma+ π- π+ 1.37%
124+
### Ξc+ -> Ξ*0 π+ (10%)
125+
4232:addChannel = 1 0.12500 0 3324 211
126+
64127
### add Xic0 decays absent in PYTHIA8 decay table
65128
4132:addChannel = 1 0.0143 0 3312 211
129+
66130
### add OmegaC decays absent in PYTHIA8 decay table
67131
4332:addChannel = 1 0.5 0 3334 211
68132
4332:addChannel = 1 0.5 0 3312 211
69133

70-
### K* -> K pi
134+
# Allow the decay of resonances in the decay chain
135+
### K*0(892) -> K- π+
71136
313:onMode = off
72137
313:onIfAll = 321 211
73-
### for Ds -> Phi pi+
138+
### K*(892)+ -> K- π0
139+
323:onMode = off
140+
323:onIfAll = 321 111
141+
### K*(1430)0 -> K- π+
142+
10311:onMode = off
143+
10311:onIfAll = 321 211
144+
### rho+ -> π+ π0
145+
213:onMode = off
146+
213:onIfAll = 211 111
147+
### φ -> K+ K-
74148
333:onMode = off
75149
333:onIfAll = 321 321
76-
### for D0 -> rho0 pi+ k-
150+
### rho0 -> π+ π-
77151
113:onMode = off
78152
113:onIfAll = 211 211
79-
### for Lambda_c -> Delta++ K-
153+
### f2(1270) -> π+ π-
154+
225:onMode = off
155+
225:onIfAll = 211 211
156+
### f0(1370) -> π+ π-
157+
10221:onMode = off
158+
10221:onIfAll = 211 211
159+
### eta -> π+ π-
160+
221:onMode = off
161+
221:onIfAll = 111 211 211
162+
### for Λc -> Delta++ K-
80163
2224:onMode = off
81164
2224:onIfAll = 2212 211
82-
### for Lambda_c -> Lambda(1520) K-
165+
### for Λc -> Lambda(1520) K-
83166
102134:onMode = off
84167
102134:onIfAll = 2212 321
85168
### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p
@@ -92,7 +175,7 @@ BeamRemnants:saturation 5
92175
3334:onMode = off
93176
3334:onIfAll = 3122 -321
94177

95-
### switch off all decay channels
178+
# switch off all decay channels
96179
411:onMode = off
97180
421:onMode = off
98181
431:onMode = off
@@ -102,41 +185,100 @@ BeamRemnants:saturation 5
102185
443:onMode = off
103186
4332:onMode = off
104187

105-
### D0 -> K pi
188+
# Allow the decay of HF
189+
### D0 -> K π
106190
421:onIfMatch = 321 211
107-
### D0 -> K pi pi0
191+
### D0 -> K π π0
108192
421:onIfMatch = 321 211 111
193+
### D0 -> rho K
194+
421:onIfMatch = 213 321
195+
### D0 -> antiK*0 π0
196+
421:onIfMatch = 313 111
197+
### D0 -> K*- π+
198+
421:onIfMatch = 323 211
199+
### D0 -> π π
200+
421:onIfMatch = 211 211
201+
### D0 -> rho+ π-
202+
421:onIfMatch = 213 -211
203+
### D0 -> π π π0
204+
421:onIfMatch = 211 211 111
205+
### D0 -> K K
206+
421:onIfMatch = 321 321
109207

110-
### D+/- -> K pi pi
208+
### D+/- -> K π π
111209
411:onIfMatch = 321 211 211
112-
### D+/- -> K* pi
210+
### D+/- -> K π π
211+
411:onIfMatch = 321 211 211 111
212+
### D+/- -> K* π
113213
411:onIfMatch = 313 211
114-
### D+/- -> phi pi
214+
### D+/- -> K* K
215+
411:onIfMatch = 313 321
216+
### D+/- -> antiK* π
217+
411:onIfMatch = 10311 211
218+
### D+/- -> antiK* K
219+
411:onIfMatch = 10311 321
220+
### D+/- -> φ π
115221
411:onIfMatch = 333 211
222+
### D+/- -> K K π
223+
411:onIfMatch = 321 321 211
224+
### D+/- -> f2(1270) π
225+
411:onIfMatch = 225 211
226+
### D+/- -> rho π
227+
411:onIfMatch = 113 211
228+
### D+/- -> π π π
229+
411:onIfMatch = 211 211 211
116230

117-
### D_s -> K K*
118-
431:onIfMatch = 321 313
119-
### D_s -> Phi pi
231+
### Ds -> φ π
120232
431:onIfMatch = 333 211
233+
### Ds -> K* K
234+
431:onIfMatch = 313 321
235+
### Ds -> φ rho+
236+
431:onIfMatch = 333 213
237+
### Ds -> rho π
238+
431:onIfMatch = 113 211
239+
### Ds -> f2 π
240+
431:onIfMatch = 225 211
241+
### Ds -> π π π
242+
431:onIfMatch = 211 211 211
243+
### Ds -> K*(892)0 π+
244+
431:onIfMatch = 313 211
245+
### Ds -> f0(1370) K+
246+
431:onIfMatch = 10221 321
247+
### Ds -> rho0 K+
248+
431:onIfMatch = 113 321
249+
### Ds -> K π π
250+
431:onIfMatch = 321 211 211
251+
### Ds -> eta π+
252+
431:onIfMatch = 221 211
121253

122-
### Lambda_c -> p K*
254+
### Λc -> pK0s
255+
4122:onIfMatch = 2212 311
256+
### Λc -> p K- π+ π0
257+
4122:onIfMatch = 2212 321 211
258+
### Λc -> p K*
123259
4122:onIfMatch = 2212 313
124-
### Lambda_c -> Delta K
260+
### Λc -> Delta++ K
125261
4122:onIfMatch = 2224 321
126-
### Lambda_c -> Lambda(1520) pi
262+
### Λc -> Lambda(1520) π
127263
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
264+
### Λc -> p K- π+ π0
265+
4122:onIfMatch = 2212 321 211 111
266+
### Λc -> p π π
267+
4122:onIfMatch = 2212 211 211
268+
### Λc -> p K K
269+
4122:onIfMatch = 2212 333
132270

133-
### Xic+ -> pK*0
134-
4232:onIfMatch = 2212 313
135-
### Xic+ -> p K- pi+
271+
### Ξc+ -> p K- π+
136272
4232:onIfMatch = 2212 321 211
137-
### Xic+ -> Xi*0 pi+, Xi*->Xi- pi+
273+
### Ξc+ -> p antiK*0(892)
274+
4232:onIfMatch = 2212 313
275+
### Ξc+ -> p φ
276+
4232:onIfMatch = 2212 333
277+
### Ξc+ -> sigma- π+ π+
278+
4232:onIfMatch = 3222 211 211
279+
### Ξc+ -> Ξ*0 π+, Ξ*0 -> Ξ- π+
138280
4232:onIfMatch = 3324 211
139-
### Xic+ -> Xi- pi+ pi+
281+
### Ξc+ -> Ξ- π+ π+
140282
4232:onIfMatch = 3312 211 211
141283

142284
### Xic0 -> Xi- pi+
@@ -145,4 +287,4 @@ BeamRemnants:saturation 5
145287
### Omega_c -> Omega pi
146288
4332:onIfMatch = 3334 211
147289
### Omega_c -> Xi pi
148-
4332:onIfMatch = 3312 211
290+
4332:onIfMatch = 3312 211

0 commit comments

Comments
 (0)