Skip to content

Commit 7436b66

Browse files
stefanopolitanostefanopolitano
committed
[PWGHF] Updating config for pp 13.6, 5.36 TeV and PbPb 5.36 TeV
1 parent 05fde30 commit 7436b66

File tree

3 files changed

+322
-76
lines changed

3 files changed

+322
-76
lines changed

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

Lines changed: 107 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -41,45 +41,64 @@ BeamRemnants:saturation 5
4141
# OmegaC
4242
4332:tau0 = 0.0803
4343

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
44+
### HF decays
45+
### BR are set to yield 50% the golden channel and 50% the correlated bkg ones
46+
###
47+
### D0 decays
48+
421:oneChannel = 1 0.500 0 -321 211 ### D0 -> K- π+
49+
421:addChannel = 1 0.100 0 -321 211 111 ### D0 -> K- π+ π0
50+
421:addChannel = 1 0.100 0 -211 211 ### D0 -> π- π+
51+
421:addChannel = 1 0.100 0 -211 211 111 ### D0 -> π- π+ π0
52+
421:addChannel = 1 0.100 0 -321 321 ### D0 -> K- K+
53+
421:addChannel = 1 0.100 0 -321 321 111 ### D0 -> K- K+ π0
54+
55+
### D+ decays
56+
411:oneChannel = 1 0.50 0 -321 211 211 ### D+ -> K- π+ π+
57+
411:addChannel = 1 0.17 0 -211 211 211 ### D+ -> π- π+ π+
58+
411:addChannel = 1 0.16 0 313 211 ### D+ -> K* π+ -> K- π+ π+
59+
411:addChannel = 1 0.17 0 333 211 ### D+ -> φ π+ -> K- K+ π+
60+
61+
### Ds+ decays
62+
431:oneChannel = 1 0.500 0 333 211 ### Ds+ -> φ π+ -> K- K+ π+
63+
431:addChannel = 1 0.125 0 313 321 ### Ds+ -> K* K+ -> K- K+ π+
64+
431:addChannel = 1 0.125 0 -211 211 211 ### Ds+ -> π- π+ π+
65+
431:addChannel = 1 0.125 0 -211 321 211 ### Ds+ -> π- K+ π+
66+
431:addChannel = 1 0.125 0 -321 321 321 ### Ds+ -> K- K+ K+
67+
68+
## Lc decays
69+
4122:oneChannel = 1 0.2500 2212 -321 211 ### Λc+ -> p K- π+
70+
4122:addChannel = 1 0.0560 100 2212 -313 ### Λc+ -> p K* -> p K- π+
71+
4122:addChannel = 1 0.0550 100 2224 -321 ### Λc+ -> Delta++ K- -> p K- π+
72+
4122:addChannel = 1 0.0550 100 102134 211 ### Λc+ -> Lambda(1520) K- -> p K- π+
73+
4122:addChannel = 1 0.1670 0 2212 211 211 ### Λc+ -> p π+ π+
74+
4122:addChannel = 1 0.2500 0 2212 311 ### Λc+ -> p K0s -> p π+ π-
75+
4122:addChannel = 1 0.1670 0 2212 -321 321 ### Λc+ -> p K- K+
76+
77+
### Ξc+ decays
78+
4232:oneChannel = 1 0.50 0 2212 -321 211 ### Ξc+ -> p K- π+
79+
4232:addChannel = 1 0.25 0 -3312 211 211 ### Ξc+ -> Ξ- π+ π+
80+
4232:addChannel = 1 0.25 0 2212 333 ### Ξc+ -> p φ -> p K- K+
81+
6482
### add Xic0 decays absent in PYTHIA8 decay table
6583
4132:addChannel = 1 0.0143 0 3312 211
6684
### add OmegaC decays absent in PYTHIA8 decay table
6785
4332:addChannel = 1 0.5 0 3334 211
6886
4332:addChannel = 1 0.5 0 3312 211
6987

70-
### K* -> K pi
88+
# Allow the decay of resonances in the decay chain
89+
### K* -> K π
7190
313:onMode = off
7291
313:onIfAll = 321 211
73-
### for Ds -> Phi pi+
92+
### for Ds -> φ π+
7493
333:onMode = off
7594
333:onIfAll = 321 321
76-
### for D0 -> rho0 pi+ k-
95+
### for D0 -> rho0 π+ k-
7796
113:onMode = off
7897
113:onIfAll = 211 211
79-
### for Lambda_c -> Delta++ K-
98+
### for Λc -> Delta++ K-
8099
2224:onMode = off
81100
2224:onIfAll = 2212 211
82-
### for Lambda_c -> Lambda(1520) K-
101+
### for Λc -> Lambda(1520) K-
83102
102134:onMode = off
84103
102134:onIfAll = 2212 321
85104
### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p
@@ -92,6 +111,69 @@ BeamRemnants:saturation 5
92111
3334:onMode = off
93112
3334:onIfAll = 3122 -321
94113

114+
### Switch off all decay channels
115+
411:onMode = off
116+
421:onMode = off
117+
431:onMode = off
118+
4122:onMode = off
119+
4232:onMode = off
120+
121+
# Allow the decay of HF
122+
### D0 -> K π
123+
421:onIfMatch = 321 211
124+
### D0 -> K π π0
125+
421:onIfMatch = 321 211 111
126+
### D0 -> π π
127+
421:onIfMatch = 211 211
128+
### D0 -> π π π0
129+
421:onIfMatch = 211 211 111
130+
### D0 -> K K
131+
421:onIfMatch = 321 321
132+
### D0 -> K K π0
133+
421:onIfMatch = 321 321 111
134+
135+
### D+/- -> K π π
136+
411:onIfMatch = 321 211 211
137+
### D+/- -> K* π
138+
411:onIfMatch = 313 211
139+
### D+/- -> φ π
140+
411:onIfMatch = 333 211
141+
### D+/- -> π π π
142+
411:onIfMatch = 211 211 211
143+
144+
### Ds -> φ π
145+
431:onIfMatch = 333 211
146+
### Ds -> K* K
147+
431:onIfMatch = 321 313
148+
### Ds -> π π π
149+
431:onIfMatch = 211 211 211
150+
### Ds -> K π π
151+
431:onIfMatch = 321 211 211
152+
### Ds -> K K K
153+
431:onIfMatch = 321 321 321
154+
155+
### Λc -> p K π
156+
4122:onIfMatch = 2212 321 211
157+
### Λc -> p K*
158+
4122:onIfMatch = 2212 313
159+
### Λc -> Delta++ K
160+
4122:onIfMatch = 2224 321
161+
### Λc -> Lambda(1520) π
162+
4122:onIfMatch = 102134 211
163+
### Λc -> p π π
164+
4122:onIfMatch = 2212 211 211
165+
### Λc -> pK0s
166+
4122:onIfMatch = 2212 311
167+
### Λc -> p K K
168+
4122:onIfMatch = 2212 321 321
169+
170+
### Ξc+ -> p K- π+
171+
4232:onIfMatch = 2212 321 211
172+
### Ξc+ -> Ξ- π+ π+
173+
4232:onIfMatch = 3312 211 211
174+
### Ξc+ -> p φ
175+
4232:onIfMatch = 2212 333
176+
95177
### switch off all decay channels
96178
411:onMode = off
97179
421:onMode = off

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

Lines changed: 107 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -45,45 +45,64 @@ BeamRemnants:saturation 5
4545
# OmegaC
4646
4332:tau0 = 0.0803
4747

48-
### Force golden charm hadrons decay modes for D2H studies
49-
### add D0 decays absent in PYTHIA8 decay table and set BRs from PDG for other
50-
421:oneChannel = 1 0.0389 0 -321 211
51-
421:addChannel = 1 0.00389 0 -321 211 111
52-
### add D+ decays absent in PYTHIA8 decay table and set BRs from PDG for other
53-
411:oneChannel = 1 0.0752 0 -321 211 211
54-
411:addChannel = 1 0.0104 0 -313 211
55-
411:addChannel = 1 0.0156 0 311 211
56-
411:addChannel = 1 0.0752 0 333 211 # to have the same amount of D+->KKpi and D+->Kpipi
57-
## add Lc decays absent in PYTHIA8 decay table and set BRs from PDG for other
58-
4122:oneChannel = 1 0.0196 100 2212 -313
59-
4122:addChannel = 1 0.0108 100 2224 -321
60-
4122:addChannel = 1 0.022 100 102134 211
61-
4122:addChannel = 1 0.035 0 2212 -321 211
62-
4122:addChannel = 1 0.0159 0 2212 311
63-
### add Xic+ decays absent in PYTHIA8 decay table
64-
4232:addChannel = 1 0.2 0 2212 -313
65-
4232:addChannel = 1 0.2 0 2212 -321 211
66-
4232:addChannel = 1 0.2 0 3324 211
67-
4232:addChannel = 1 0.2 0 3312 211 211
48+
### HF decays
49+
### BR are set to yield 50% the golden channel and 50% the correlated bkg ones
50+
###
51+
### D0 decays
52+
421:oneChannel = 1 0.500 0 -321 211 ### D0 -> K- π+
53+
421:addChannel = 1 0.100 0 -321 211 111 ### D0 -> K- π+ π0
54+
421:addChannel = 1 0.100 0 -211 211 ### D0 -> π- π+
55+
421:addChannel = 1 0.100 0 -211 211 111 ### D0 -> π- π+ π0
56+
421:addChannel = 1 0.100 0 -321 321 ### D0 -> K- K+
57+
421:addChannel = 1 0.100 0 -321 321 111 ### D0 -> K- K+ π0
58+
59+
### D+ decays
60+
411:oneChannel = 1 0.50 0 -321 211 211 ### D+ -> K- π+ π+
61+
411:addChannel = 1 0.17 0 -211 211 211 ### D+ -> π- π+ π+
62+
411:addChannel = 1 0.16 0 313 211 ### D+ -> K* π+ -> K- π+ π+
63+
411:addChannel = 1 0.17 0 333 211 ### D+ -> φ π+ -> K- K+ π+
64+
65+
### Ds+ decays
66+
431:oneChannel = 1 0.500 0 333 211 ### Ds+ -> φ π+ -> K- K+ π+
67+
431:addChannel = 1 0.125 0 313 321 ### Ds+ -> K* K+ -> K- K+ π+
68+
431:addChannel = 1 0.125 0 -211 211 211 ### Ds+ -> π- π+ π+
69+
431:addChannel = 1 0.125 0 -211 321 211 ### Ds+ -> π- K+ π+
70+
431:addChannel = 1 0.125 0 -321 321 321 ### Ds+ -> K- K+ K+
71+
72+
## Lc decays
73+
4122:oneChannel = 1 0.2500 2212 -321 211 ### Λc+ -> p K- π+
74+
4122:addChannel = 1 0.0560 100 2212 -313 ### Λc+ -> p K* -> p K- π+
75+
4122:addChannel = 1 0.0550 100 2224 -321 ### Λc+ -> Delta++ K- -> p K- π+
76+
4122:addChannel = 1 0.0550 100 102134 211 ### Λc+ -> Lambda(1520) K- -> p K- π+
77+
4122:addChannel = 1 0.1670 0 2212 211 211 ### Λc+ -> p π+ π+
78+
4122:addChannel = 1 0.2500 0 2212 311 ### Λc+ -> p K0s -> p π+ π-
79+
4122:addChannel = 1 0.1670 0 2212 -321 321 ### Λc+ -> p K- K+
80+
81+
### Ξc+ decays
82+
4232:oneChannel = 1 0.50 0 2212 -321 211 ### Ξc+ -> p K- π+
83+
4232:addChannel = 1 0.25 0 -3312 211 211 ### Ξc+ -> Ξ- π+ π+
84+
4232:addChannel = 1 0.25 0 2212 333 ### Ξc+ -> p φ -> p K- K+
85+
6886
### add Xic0 decays absent in PYTHIA8 decay table
6987
4132:addChannel = 1 0.0143 0 3312 211
7088
### add OmegaC decays absent in PYTHIA8 decay table
7189
4332:addChannel = 1 0.5 0 3334 211
7290
4332:addChannel = 1 0.5 0 3312 211
7391

74-
### K* -> K pi
92+
# Allow the decay of resonances in the decay chain
93+
### K* -> K π
7594
313:onMode = off
7695
313:onIfAll = 321 211
77-
### for Ds -> Phi pi+
96+
### for Ds -> φ π+
7897
333:onMode = off
7998
333:onIfAll = 321 321
80-
### for D0 -> rho0 pi+ k-
99+
### for D0 -> rho0 π+ k-
81100
113:onMode = off
82101
113:onIfAll = 211 211
83-
### for Lambda_c -> Delta++ K-
102+
### for Λc -> Delta++ K-
84103
2224:onMode = off
85104
2224:onIfAll = 2212 211
86-
### for Lambda_c -> Lambda(1520) K-
105+
### for Λc -> Lambda(1520) K-
87106
102134:onMode = off
88107
102134:onIfAll = 2212 321
89108
### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p
@@ -96,6 +115,69 @@ BeamRemnants:saturation 5
96115
3334:onMode = off
97116
3334:onIfAll = 3122 -321
98117

118+
### Switch off all decay channels
119+
411:onMode = off
120+
421:onMode = off
121+
431:onMode = off
122+
4122:onMode = off
123+
4232:onMode = off
124+
125+
# Allow the decay of HF
126+
### D0 -> K π
127+
421:onIfMatch = 321 211
128+
### D0 -> K π π0
129+
421:onIfMatch = 321 211 111
130+
### D0 -> π π
131+
421:onIfMatch = 211 211
132+
### D0 -> π π π0
133+
421:onIfMatch = 211 211 111
134+
### D0 -> K K
135+
421:onIfMatch = 321 321
136+
### D0 -> K K π0
137+
421:onIfMatch = 321 321 111
138+
139+
### D+/- -> K π π
140+
411:onIfMatch = 321 211 211
141+
### D+/- -> K* π
142+
411:onIfMatch = 313 211
143+
### D+/- -> φ π
144+
411:onIfMatch = 333 211
145+
### D+/- -> π π π
146+
411:onIfMatch = 211 211 211
147+
148+
### Ds -> φ π
149+
431:onIfMatch = 333 211
150+
### Ds -> K* K
151+
431:onIfMatch = 321 313
152+
### Ds -> π π π
153+
431:onIfMatch = 211 211 211
154+
### Ds -> K π π
155+
431:onIfMatch = 321 211 211
156+
### Ds -> K K K
157+
431:onIfMatch = 321 321 321
158+
159+
### Λc -> p K π
160+
4122:onIfMatch = 2212 321 211
161+
### Λc -> p K*
162+
4122:onIfMatch = 2212 313
163+
### Λc -> Delta++ K
164+
4122:onIfMatch = 2224 321
165+
### Λc -> Lambda(1520) π
166+
4122:onIfMatch = 102134 211
167+
### Λc -> p π π
168+
4122:onIfMatch = 2212 211 211
169+
### Λc -> pK0s
170+
4122:onIfMatch = 2212 311
171+
### Λc -> p K K
172+
4122:onIfMatch = 2212 321 321
173+
174+
### Ξc+ -> p K- π+
175+
4232:onIfMatch = 2212 321 211
176+
### Ξc+ -> Ξ- π+ π+
177+
4232:onIfMatch = 3312 211 211
178+
### Ξc+ -> p φ
179+
4232:onIfMatch = 2212 333
180+
99181
### switch off all decay channels
100182
411:onMode = off
101183
421:onMode = off

0 commit comments

Comments
 (0)