Skip to content

Commit 506cec9

Browse files
stefanopolitanostefanopolitano
authored
[PWGHF] Add config for correlated bkg production in pp (#1965)
* [PWGHF] Add config and ini file for correlated bkg production in pp * [PWGHF] minor fix in ini file * [PWGHF] minor fix in config file
1 parent ab53525 commit 506cec9

File tree

2 files changed

+166
-0
lines changed

2 files changed

+166
-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(5, -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: 158 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,158 @@
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+
104+
# Allow the decay of HF
105+
### D0 -> K π
106+
421:onIfMatch = 321 211
107+
### D0 -> K π π0
108+
421:onIfMatch = 321 211 111
109+
### D0 -> π π
110+
421:onIfMatch = 211 211
111+
### D0 -> π π π0
112+
421:onIfMatch = 211 211 111
113+
### D0 -> K K
114+
421:onIfMatch = 321 321
115+
### D0 -> K K π0
116+
421:onIfMatch = 321 321 111
117+
118+
### D+/- -> K π π
119+
411:onIfMatch = 321 211 211
120+
### D+/- -> K* π
121+
411:onIfMatch = 313 211
122+
### D+/- -> φ π
123+
411:onIfMatch = 333 211
124+
### D+/- -> π π π
125+
411:onIfMatch = 211 211 211
126+
127+
### Ds -> φ π
128+
431:onIfMatch = 333 211
129+
### Ds -> K* K
130+
431:onIfMatch = 321 313
131+
### Ds -> π π π
132+
431:onIfMatch = 211 211 211
133+
### Ds -> K π π
134+
431:onIfMatch = 321 211 211
135+
### Ds -> K K K
136+
431:onIfMatch = 321 321 321
137+
138+
### Λc -> p K π
139+
4122:onIfMatch = 2212 321 211
140+
### Λc -> p K*
141+
4122:onIfMatch = 2212 313
142+
### Λc -> Delta++ K
143+
4122:onIfMatch = 2224 321
144+
### Λc -> Lambda(1520) π
145+
4122:onIfMatch = 102134 211
146+
### Λc -> p π π
147+
4122:onIfMatch = 2212 211 211
148+
### Λc -> pK0s
149+
4122:onIfMatch = 2212 311
150+
### Λc -> p K K
151+
4122:onIfMatch = 2212 321 321
152+
153+
### Ξc+ -> p K- π+
154+
4232:onIfMatch = 2212 321 211
155+
### Ξc+ -> Ξ- π+ π+
156+
4232:onIfMatch = 3312 211 211
157+
### Ξc+ -> p φ
158+
4232:onIfMatch = 2212 333

0 commit comments

Comments
 (0)