Skip to content

Forced decays + other energies #2004

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
dsekihat merged 1 commit into from
May 20, 2025
Merged

Forced decays + other energies #2004

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
[GeneratorExternal]
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332")
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;511;521;531;5122;5132;5232;5332")

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

[DecayerPythia8]
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg
4 changes: 2 additions & 2 deletions MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail_pp13.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
[GeneratorExternal]
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332")
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;511;521;531;5122;5132;5232;5332")

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

[DecayerPythia8]
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg
13 changes: 13 additions & 0 deletions MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail_pp502.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
[GeneratorExternal]
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;511;521;531;5122;5132;5232;5332")

[GeneratorPythia8]
config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp502.cfg
hooksFileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
hooksFuncName = pythia8_userhooks_bbbar(-12.,12.)
includePartonEvent=true

[DecayerPythia8]
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg
13 changes: 13 additions & 0 deletions MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail_pp536.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
[GeneratorExternal]
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;511;521;531;5122;5132;5232;5332")

[GeneratorPythia8]
config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp536.cfg
hooksFileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
hooksFuncName = pythia8_userhooks_bbbar(-12.,12.)
includePartonEvent=true

[DecayerPythia8]
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg
13 changes: 13 additions & 0 deletions MC/config/PWGEM/ini/Pythia8_Charm_Cocktail_pp502.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
[GeneratorExternal]
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332")

[GeneratorPythia8]
config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp502.cfg
hooksFileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
hooksFuncName = pythia8_userhooks_ccbar(-12.,12.)
includePartonEvent=true

[DecayerPythia8]
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg
13 changes: 13 additions & 0 deletions MC/config/PWGEM/ini/Pythia8_Charm_Cocktail_pp536.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
[GeneratorExternal]
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332")

[GeneratorPythia8]
config = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp536.cfg
hooksFileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/hooks/pythia8_userhooks_qqbar.C
hooksFuncName = pythia8_userhooks_ccbar(-12.,12.)
includePartonEvent=true

[DecayerPythia8]
config[0] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
config[1] = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg
77 changes: 77 additions & 0 deletions MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,77 @@
### semileptonic decays for charm
### D+
411:oneChannel = 1 0.087 0 -311 -11 12
411:addChannel = 1 0.040 0 -321 211 -11 12
411:addChannel = 1 0.037 0 -313 -11 12
### D0
421:oneChannel = 1 0.035 0 -321 -11 12
421:addChannel = 1 0.022 0 -323 -11 12
421:addChannel = 1 0.016 0 -321 111 -11 12
421:addChannel = 1 0.014 0 -311 -211 -11 12
### Ds
431:oneChannel = 1 0.025 0 333 -11 12
431:addChannel = 1 0.027 0 221 -11 12
### Lambdac
4122:oneChannel = 1 0.036 0 3122 -11 12
### chi_{c}^{+}
4232:oneChannel = 1 0.07 0 3322 -11 12
### chi_{c}^{0}
4132:oneChannel = 1 0.014 0 3312 -11 12
### Omega_{c}
4332:oneChannel = 1 0.01224 0 3334 -11 12

### semileptonic decays for beauty
### B0
511:oneChannel = 1 0.0207000 0 12 -11 -411
511:addChannel = 1 0.0570000 0 12 -11 -413
511:addChannel = 1 0.0023000 0 12 -11 -415
511:addChannel = 1 0.0001330 0 12 -11 -211
511:addChannel = 1 0.0002690 0 12 -11 -213
511:addChannel = 1 0.0045000 0 12 -11 -10411
511:addChannel = 1 0.0052000 0 12 -11 -10413
511:addChannel = 1 0.0083000 0 12 -11 -20413

### B+
521:oneChannel = 1 0.0000720 0 12 -11 111
521:addChannel = 1 0.0001450 0 12 -11 113
521:addChannel = 1 0.0000840 0 12 -11 221
521:addChannel = 1 0.0001450 0 12 -11 223
521:addChannel = 1 0.0000840 0 12 -11 331
521:addChannel = 1 0.0224000 0 12 -11 -421
521:addChannel = 1 0.0617000 0 12 -11 -423
521:addChannel = 1 0.0030000 0 12 -11 -425
521:addChannel = 1 0.0049000 0 12 -11 -10421
521:addChannel = 1 0.0056000 0 12 -11 -10423
521:addChannel = 1 0.0090000 0 12 -11 -20423

### Bs
531:oneChannel = 1 0.0002000 0 12 -11 -321
531:addChannel = 1 0.0003000 0 12 -11 -323
531:addChannel = 1 0.0210000 0 12 -11 -431
531:addChannel = 1 0.0490000 0 12 -11 -433
531:addChannel = 1 0.0070000 0 12 -11 -435
531:addChannel = 1 0.0003000 0 12 -11 -10323
531:addChannel = 1 0.0040000 0 12 -11 -10431
531:addChannel = 1 0.0070000 0 12 -11 -10433
531:addChannel = 1 0.0002000 0 12 -11 -20323
531:addChannel = 1 0.0040000 0 12 -11 -20433

### Lambdab
5122:oneChannel = 1 0.0546000 0 -12 11 4122
5122:addChannel = 1 0.0096000 0 -12 11 4124
5122:addChannel = 1 0.0128000 0 -12 11 14122

### Chi_{b}^{-}
5132:oneChannel = 1 0.1080010 0 -12 11 4 3101
5132:addChannel = 1 0.0020000 0 -12 11 2 3101
### Chi_{b}^{0}
5232:oneChannel = 1 0.1080010 0 -12 11 4 3201
5232:addChannel = 1 0.0020000 0 -12 11 2 3201
### Omega_{b}^{-}
5332:oneChannel = 1 0.1080010 1 -12 11 4 3303
5332:oneChannel = 1 0.0020000 1 -12 11 2 3303

# Correct OmegaC decay length (wrong in PYTHIA8 decay table) (mm/c)
4332:tau0 = 0.08000000000
# Correct Lb decay length (wrong in PYTHIA8 decay table)
5122:tau0 = 4.41000e-01
30 changes: 30 additions & 0 deletions MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp502.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,30 @@
### beams
Beams:idA 2212 # proton
Beams:idB 2212 # proton
Beams:eCM 5200. # GeV

### processes
SoftQCD:inelastic on

### decays
ParticleDecays:limitTau0 on
ParticleDecays:tau0Max 10.

### switch on color reconnection in mode 2 (https://arxiv.org/pdf/1505.01681.pdf)
Tune:pp = 14
ColourReconnection:mode = 1
ColourReconnection:allowDoubleJunRem = off
ColourReconnection:m0 = 0.3
ColourReconnection:allowJunctions = on
ColourReconnection:junctionCorrection = 1.20
ColourReconnection:timeDilationMode = 2
ColourReconnection:timeDilationPar = 0.18
StringPT:sigma = 0.335
StringZ:aLund = 0.36
StringZ:bLund = 0.56
StringFlav:probQQtoQ = 0.078
StringFlav:ProbStoUD = 0.2
StringFlav:probQQ1toQQ0join = 0.0275,0.0275,0.0275,0.0275
MultiPartonInteractions:pT0Ref = 2.15
BeamRemnants:remnantMode = 1
BeamRemnants:saturation =5
30 changes: 30 additions & 0 deletions MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail_pp536.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,30 @@
### beams
Beams:idA 2212 # proton
Beams:idB 2212 # proton
Beams:eCM 5360. # GeV

### processes
SoftQCD:inelastic on

### decays
ParticleDecays:limitTau0 on
ParticleDecays:tau0Max 10.

### switch on color reconnection in mode 2 (https://arxiv.org/pdf/1505.01681.pdf)
Tune:pp = 14
ColourReconnection:mode = 1
ColourReconnection:allowDoubleJunRem = off
ColourReconnection:m0 = 0.3
ColourReconnection:allowJunctions = on
ColourReconnection:junctionCorrection = 1.20
ColourReconnection:timeDilationMode = 2
ColourReconnection:timeDilationPar = 0.18
StringPT:sigma = 0.335
StringZ:aLund = 0.36
StringZ:bLund = 0.56
StringFlav:probQQtoQ = 0.078
StringFlav:ProbStoUD = 0.2
StringFlav:probQQ1toQQ0join = 0.0275,0.0275,0.0275,0.0275
MultiPartonInteractions:pT0Ref = 2.15
BeamRemnants:remnantMode = 1
BeamRemnants:saturation =5