Skip to content

Commit 2f69c38

Browse files
klsmith15kKrista Lizbeth Smith
klsmith15k
and
Krista Lizbeth Smith
authored
[PWGHF] Configuration for charm-baryon enriched MC anchored to pp ref run (#1843)
* pp ref eCM * pp ref eCM --------- Co-authored-by: Krista Lizbeth Smith <krista.lizbeth.smith@cern.ch>
1 parent 67a59f0 commit 2f69c38

File tree

2 files changed

+156
-0
lines changed

2 files changed

+156
-0
lines changed

‎MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_CharmBaryons_pp_ref.ini‎

Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,8 @@
1+
### The external generator derives from GeneratorPythia8.
2+
[GeneratorExternal]
3+
fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C
4+
funcName=GeneratorPythia8GapTriggeredCharmAndBeauty(5, -1.5, 1.5, -1.5, 1.5, {4132, 4332, 4232})
5+
6+
[GeneratorPythia8]
7+
config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg
8+
includePartonEvent=true

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

Lines changed: 148 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,148 @@
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

Comments
 (0)