Skip to content

Commit 6094769

Browse files
stefanopolitanostefanopolitano
committed
[PWGHF] Add config and ini file for correlated bkg production in pp
1 parent e09374d commit 6094769

File tree

2 files changed

+167
-0
lines changed

2 files changed

+167
-0
lines changed

‎MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkg.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(1, -1.5, 1.5)
5+
6+
[GeneratorPythia8]
7+
config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_CorrBkg.cfg
8+
includePartonEvent=true

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

Lines changed: 159 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,159 @@
1+
### authors: Stefano Politano (stefano.politano@cern.ch)
2+
### last update: April 2025
3+
4+
### beams
5+
Beams:idA 2212 # proton
6+
Beams:idB 2212 # proton
7+
Beams:eCM 13600. # GeV
8+
9+
### processes
10+
SoftQCD:inelastic on # all inelastic processes
11+
12+
### decays
13+
ParticleDecays:limitTau0 on
14+
ParticleDecays:tau0Max 10.
15+
16+
### switching on Pythia Mode2
17+
ColourReconnection:mode 1
18+
ColourReconnection:allowDoubleJunRem off
19+
ColourReconnection:m0 0.3
20+
ColourReconnection:allowJunctions on
21+
ColourReconnection:junctionCorrection 1.20
22+
ColourReconnection:timeDilationMode 2
23+
ColourReconnection:timeDilationPar 0.18
24+
StringPT:sigma 0.335
25+
StringZ:aLund 0.36
26+
StringZ:bLund 0.56
27+
StringFlav:probQQtoQ 0.078
28+
StringFlav:ProbStoUD 0.2
29+
StringFlav:probQQ1toQQ0join 0.0275,0.0275,0.0275,0.0275
30+
MultiPartonInteractions:pT0Ref 2.15
31+
BeamRemnants:remnantMode 1
32+
BeamRemnants:saturation 5
33+
34+
# Correct decay lengths (wrong in PYTHIA8 decay table)
35+
# Lb
36+
5122:tau0 = 0.4390
37+
# Xic0
38+
4132:tau0 = 0.0455
39+
# OmegaC
40+
4332:tau0 = 0.0803
41+
42+
### HF decays
43+
### BR are set to yield the same amount of counts for each final state
44+
###
45+
### D0 decays
46+
421:oneChannel = 1 0.170 0 -321 211 ### D0 -> K- π+
47+
421:addChannel = 1 0.166 0 -321 211 111 ### D0 -> K- π+ π0
48+
421:addChannel = 1 0.166 0 -211 211 ### D0 -> π- π+
49+
421:addChannel = 1 0.166 0 -211 211 111 ### D0 -> π- π+ π0
50+
421:addChannel = 1 0.166 0 -321 321 ### D0 -> K- K+
51+
421:addChannel = 1 0.166 0 -321 321 111 ### D0 -> K- K+ π0
52+
53+
### D+ decays
54+
411:oneChannel = 1 0.17 0 -321 211 211 ### D+ -> K- π+ π+
55+
411:addChannel = 1 0.33 0 -211 211 211 ### D+ -> π- π+ π+
56+
411:addChannel = 1 0.17 0 313 211 ### D+ -> K* π+ -> K- π+ π+
57+
411:addChannel = 1 0.33 0 333 211 ### D+ -> φ π+ -> K- K+ π+
58+
59+
### Ds+ decays
60+
431:oneChannel = 1 0.125 0 333 211 ### Ds+ -> φ π+ -> K- K+ π+
61+
431:addChannel = 1 0.125 0 313 321 ### Ds+ -> K* K+ -> K- K+ π+
62+
431:addChannel = 1 0.250 0 -211 211 211 ### Ds+ -> π- π+ π+
63+
431:addChannel = 1 0.250 0 -211 321 211 ### Ds+ -> π- K+ π+
64+
431:addChannel = 1 0.250 0 -321 321 321 ### Ds+ -> K- K+ K+
65+
66+
## Lc decays
67+
4122:oneChannel = 1 0.0825 0 2212 -321 211 ### Λc+ -> p K- π+
68+
4122:addChannel = 1 0.0825 100 2212 -313 ### Λc+ -> p K* -> p K- π+
69+
4122:addChannel = 1 0.0825 100 2224 -321 ### Λc+ -> Delta++ K- -> p K- π+
70+
4122:addChannel = 1 0.0825 100 102134 211 ### Λc+ -> Lambda(1520) K- -> p K- π+
71+
4122:addChannel = 1 0.1650 0 2212 211 211 ### Λc+ -> p π+ π+
72+
4122:addChannel = 1 0.1650 0 2212 311 ### Λc+ -> p K0s -> p π+ π-
73+
4122:addChannel = 1 0.3300 0 2212 -321 321 ### Λc+ -> p K- K+
74+
75+
### Ξc+ decays
76+
4232:oneChannel = 1 0.34 0 2212 -321 211 ### Ξc+ -> p K- π+
77+
4232:addChannel = 1 0.33 0 -3312 211 211 ### Ξc+ -> Ξ- π+ π+
78+
4232:addChannel = 1 0.33 0 2212 333 ### Ξc+ -> p φ -> p K- K+
79+
80+
# Allow the decay of resonances in the decay chain
81+
### K* -> K π
82+
313:onMode = off
83+
313:onIfAll = 321 211
84+
### for Ds -> φ π+
85+
333:onMode = off
86+
333:onIfAll = 321 321
87+
### for D0 -> rho0 π+ k-
88+
113:onMode = off
89+
113:onIfAll = 211 211
90+
### for Λc -> Delta++ K-
91+
2224:onMode = off
92+
2224:onIfAll = 2212 211
93+
### for Λc -> Lambda(1520) K-
94+
102134:onMode = off
95+
102134:onIfAll = 2212 321
96+
97+
### Switch off all decay channels
98+
411:onMode = off
99+
421:onMode = off
100+
431:onMode = off
101+
4122:onMode = off
102+
4232:onMode = off
103+
4332:onMode = off
104+
105+
# Allow the decay of HF
106+
### D0 -> K π
107+
421:onIfMatch = 321 211
108+
### D0 -> K π π0
109+
421:onIfMatch = 321 211 111
110+
### D0 -> π π
111+
421:onIfMatch = 211 211
112+
### D0 -> π π π0
113+
421:onIfMatch = 211 211 111
114+
### D0 -> K K
115+
421:onIfMatch = 321 321
116+
### D0 -> K K π0
117+
421:onIfMatch = 321 321 111
118+
119+
### D+/- -> K π π
120+
411:onIfMatch = 321 211 211
121+
### D+/- -> K* π
122+
411:onIfMatch = 313 211
123+
### D+/- -> φ π
124+
411:onIfMatch = 333 211
125+
### D+/- -> π π π
126+
411:onIfMatch = 211 211 211
127+
128+
### Ds -> φ π
129+
431:onIfMatch = 333 211
130+
### Ds -> K* K
131+
431:onIfMatch = 321 313
132+
### Ds -> π π π
133+
431:onIfMatch = 211 211 211
134+
### Ds -> K π π
135+
431:onIfMatch = 321 211 211
136+
### Ds -> K K K
137+
431:onIfMatch = 321 321 321
138+
139+
### Λc -> p K π
140+
4122:onIfMatch = 2212 321 211
141+
### Λc -> p K*
142+
4122:onIfMatch = 2212 313
143+
### Λc -> Delta++ K
144+
4122:onIfMatch = 2224 321
145+
### Λc -> Lambda(1520) π
146+
4122:onIfMatch = 102134 211
147+
### Λc -> p π π
148+
4122:onIfMatch = 2212 211 211
149+
### Λc -> pK0s
150+
4122:onIfMatch = 2212 311
151+
### Λc -> p K K
152+
4122:onIfMatch = 2212 321 321
153+
154+
### Ξc+ -> p K- π+
155+
4232:onIfMatch = 2212 321 211
156+
### Ξc+ -> Ξ- π+ π+
157+
4232:onIfMatch = 3312 211 211
158+
### Ξc+ -> p φ
159+
4232:onIfMatch = 2212 333

0 commit comments

Comments
 (0)