Skip to content

Commit 88ee20d

Browse files
rbailhacrbailhac
authored
Forced decays + other energies (#2004)
1 parent 3247645 commit 88ee20d

9 files changed

+193
-4
lines changed

‎MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail.ini‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,6 @@
11
[GeneratorExternal]
22
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
3-
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332")
3+
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;511;521;531;5122;5132;5232;5332")
44

55
[GeneratorPythia8]
66
config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail.cfg
@@ -10,4 +10,4 @@ includePartonEvent=true
1010

1111
[DecayerPythia8]
1212
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
13-
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg
13+
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg

‎MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail_pp13.ini‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,6 @@
11
[GeneratorExternal]
22
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
3-
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332")
3+
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;511;521;531;5122;5132;5232;5332")
44

55
[GeneratorPythia8]
66
config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp13.cfg
@@ -10,4 +10,4 @@ includePartonEvent=true
1010

1111
[DecayerPythia8]
1212
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
13-
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg
13+
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg

‎MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail_pp502.ini‎

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
[GeneratorExternal]
2+
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
3+
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;511;521;531;5122;5132;5232;5332")
4+
5+
[GeneratorPythia8]
6+
config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp502.cfg
7+
hooksFileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
8+
hooksFuncName = pythia8_userhooks_bbbar(-12.,12.)
9+
includePartonEvent=true
10+
11+
[DecayerPythia8]
12+
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
13+
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg

‎MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail_pp536.ini‎

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
[GeneratorExternal]
2+
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
3+
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;511;521;531;5122;5132;5232;5332")
4+
5+
[GeneratorPythia8]
6+
config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp536.cfg
7+
hooksFileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
8+
hooksFuncName = pythia8_userhooks_bbbar(-12.,12.)
9+
includePartonEvent=true
10+
11+
[DecayerPythia8]
12+
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
13+
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg

‎MC/config/PWGEM/ini/Pythia8_Charm_Cocktail_pp502.ini‎

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
[GeneratorExternal]
2+
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
3+
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332")
4+
5+
[GeneratorPythia8]
6+
config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp502.cfg
7+
hooksFileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
8+
hooksFuncName = pythia8_userhooks_ccbar(-12.,12.)
9+
includePartonEvent=true
10+
11+
[DecayerPythia8]
12+
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
13+
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg

‎MC/config/PWGEM/ini/Pythia8_Charm_Cocktail_pp536.ini‎

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
[GeneratorExternal]
2+
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
3+
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332")
4+
5+
[GeneratorPythia8]
6+
config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp536.cfg
7+
hooksFileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
8+
hooksFuncName = pythia8_userhooks_ccbar(-12.,12.)
9+
includePartonEvent=true
10+
11+
[DecayerPythia8]
12+
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
13+
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg

‎MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg‎

Lines changed: 77 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,77 @@
1+
### semileptonic decays for charm
2+
### D+
3+
411:oneChannel = 1 0.087 0 -311 -11 12
4+
411:addChannel = 1 0.040 0 -321 211 -11 12
5+
411:addChannel = 1 0.037 0 -313 -11 12
6+
### D0
7+
421:oneChannel = 1 0.035 0 -321 -11 12
8+
421:addChannel = 1 0.022 0 -323 -11 12
9+
421:addChannel = 1 0.016 0 -321 111 -11 12
10+
421:addChannel = 1 0.014 0 -311 -211 -11 12
11+
### Ds
12+
431:oneChannel = 1 0.025 0 333 -11 12
13+
431:addChannel = 1 0.027 0 221 -11 12
14+
### Lambdac
15+
4122:oneChannel = 1 0.036 0 3122 -11 12
16+
### chi_{c}^{+}
17+
4232:oneChannel = 1 0.07 0 3322 -11 12
18+
### chi_{c}^{0}
19+
4132:oneChannel = 1 0.014 0 3312 -11 12
20+
### Omega_{c}
21+
4332:oneChannel = 1 0.01224 0 3334 -11 12
22+
23+
### semileptonic decays for beauty
24+
### B0
25+
511:oneChannel = 1 0.0207000 0 12 -11 -411
26+
511:addChannel = 1 0.0570000 0 12 -11 -413
27+
511:addChannel = 1 0.0023000 0 12 -11 -415
28+
511:addChannel = 1 0.0001330 0 12 -11 -211
29+
511:addChannel = 1 0.0002690 0 12 -11 -213
30+
511:addChannel = 1 0.0045000 0 12 -11 -10411
31+
511:addChannel = 1 0.0052000 0 12 -11 -10413
32+
511:addChannel = 1 0.0083000 0 12 -11 -20413
33+
34+
### B+
35+
521:oneChannel = 1 0.0000720 0 12 -11 111
36+
521:addChannel = 1 0.0001450 0 12 -11 113
37+
521:addChannel = 1 0.0000840 0 12 -11 221
38+
521:addChannel = 1 0.0001450 0 12 -11 223
39+
521:addChannel = 1 0.0000840 0 12 -11 331
40+
521:addChannel = 1 0.0224000 0 12 -11 -421
41+
521:addChannel = 1 0.0617000 0 12 -11 -423
42+
521:addChannel = 1 0.0030000 0 12 -11 -425
43+
521:addChannel = 1 0.0049000 0 12 -11 -10421
44+
521:addChannel = 1 0.0056000 0 12 -11 -10423
45+
521:addChannel = 1 0.0090000 0 12 -11 -20423
46+
47+
### Bs
48+
531:oneChannel = 1 0.0002000 0 12 -11 -321
49+
531:addChannel = 1 0.0003000 0 12 -11 -323
50+
531:addChannel = 1 0.0210000 0 12 -11 -431
51+
531:addChannel = 1 0.0490000 0 12 -11 -433
52+
531:addChannel = 1 0.0070000 0 12 -11 -435
53+
531:addChannel = 1 0.0003000 0 12 -11 -10323
54+
531:addChannel = 1 0.0040000 0 12 -11 -10431
55+
531:addChannel = 1 0.0070000 0 12 -11 -10433
56+
531:addChannel = 1 0.0002000 0 12 -11 -20323
57+
531:addChannel = 1 0.0040000 0 12 -11 -20433
58+
59+
### Lambdab
60+
5122:oneChannel = 1 0.0546000 0 -12 11 4122
61+
5122:addChannel = 1 0.0096000 0 -12 11 4124
62+
5122:addChannel = 1 0.0128000 0 -12 11 14122
63+
64+
### Chi_{b}^{-}
65+
5132:oneChannel = 1 0.1080010 0 -12 11 4 3101
66+
5132:addChannel = 1 0.0020000 0 -12 11 2 3101
67+
### Chi_{b}^{0}
68+
5232:oneChannel = 1 0.1080010 0 -12 11 4 3201
69+
5232:addChannel = 1 0.0020000 0 -12 11 2 3201
70+
### Omega_{b}^{-}
71+
5332:oneChannel = 1 0.1080010 1 -12 11 4 3303
72+
5332:oneChannel = 1 0.0020000 1 -12 11 2 3303
73+
74+
# Correct OmegaC decay length (wrong in PYTHIA8 decay table) (mm/c)
75+
4332:tau0 = 0.08000000000
76+
# Correct Lb decay length (wrong in PYTHIA8 decay table)
77+
5122:tau0 = 4.41000e-01

‎MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp502.cfg‎

Lines changed: 30 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,30 @@
1+
### beams
2+
Beams:idA 2212 # proton
3+
Beams:idB 2212 # proton
4+
Beams:eCM 5200. # GeV
5+
6+
### processes
7+
SoftQCD:inelastic on
8+
9+
### decays
10+
ParticleDecays:limitTau0 on
11+
ParticleDecays:tau0Max 10.
12+
13+
### switch on color reconnection in mode 2 (https://arxiv.org/pdf/1505.01681.pdf)
14+
Tune:pp = 14
15+
ColourReconnection:mode = 1
16+
ColourReconnection:allowDoubleJunRem = off
17+
ColourReconnection:m0 = 0.3
18+
ColourReconnection:allowJunctions = on
19+
ColourReconnection:junctionCorrection = 1.20
20+
ColourReconnection:timeDilationMode = 2
21+
ColourReconnection:timeDilationPar = 0.18
22+
StringPT:sigma = 0.335
23+
StringZ:aLund = 0.36
24+
StringZ:bLund = 0.56
25+
StringFlav:probQQtoQ = 0.078
26+
StringFlav:ProbStoUD = 0.2
27+
StringFlav:probQQ1toQQ0join = 0.0275,0.0275,0.0275,0.0275
28+
MultiPartonInteractions:pT0Ref = 2.15
29+
BeamRemnants:remnantMode = 1
30+
BeamRemnants:saturation =5

‎MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp536.cfg‎

Lines changed: 30 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,30 @@
1+
### beams
2+
Beams:idA 2212 # proton
3+
Beams:idB 2212 # proton
4+
Beams:eCM 5360. # GeV
5+
6+
### processes
7+
SoftQCD:inelastic on
8+
9+
### decays
10+
ParticleDecays:limitTau0 on
11+
ParticleDecays:tau0Max 10.
12+
13+
### switch on color reconnection in mode 2 (https://arxiv.org/pdf/1505.01681.pdf)
14+
Tune:pp = 14
15+
ColourReconnection:mode = 1
16+
ColourReconnection:allowDoubleJunRem = off
17+
ColourReconnection:m0 = 0.3
18+
ColourReconnection:allowJunctions = on
19+
ColourReconnection:junctionCorrection = 1.20
20+
ColourReconnection:timeDilationMode = 2
21+
ColourReconnection:timeDilationPar = 0.18
22+
StringPT:sigma = 0.335
23+
StringZ:aLund = 0.36
24+
StringZ:bLund = 0.56
25+
StringFlav:probQQtoQ = 0.078
26+
StringFlav:ProbStoUD = 0.2
27+
StringFlav:probQQ1toQQ0join = 0.0275,0.0275,0.0275,0.0275
28+
MultiPartonInteractions:pT0Ref = 2.15
29+
BeamRemnants:remnantMode = 1
30+
BeamRemnants:saturation =5

0 commit comments

Comments
 (0)