Skip to content

Commit 868db81

Browse files
njacazionjacazio
authored
Add pythia CLR-BLC Modes 0,2,3 (#2201)
- Based on https://link.springer.com/article/10.1007/JHEP08(2015)003
1 parent 57566cc commit 868db81

9 files changed

+182
-0
lines changed

‎MC/config/ALICE3/ini/pythia8_pp_clr-blc0_13tev.ini‎

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,9 @@
1+
[Diamond]
2+
width[2]=6.0
3+
4+
[GeneratorExternal]
5+
fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/ALICE3/pythia8/generator_pythia8_ALICE3.C
6+
funcName=generator_pythia8_ALICE3()
7+
8+
[GeneratorPythia8]
9+
config=${O2DPG_MC_CONFIG_ROOT}/MC/config/ALICE3/pythia8/generator/pythia8_pp_clr-blc0_13tev.cfg

‎MC/config/ALICE3/ini/pythia8_pp_clr-blc2_13tev.ini‎

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,9 @@
1+
[Diamond]
2+
width[2]=6.0
3+
4+
[GeneratorExternal]
5+
fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/ALICE3/pythia8/generator_pythia8_ALICE3.C
6+
funcName=generator_pythia8_ALICE3()
7+
8+
[GeneratorPythia8]
9+
config=${O2DPG_MC_CONFIG_ROOT}/MC/config/ALICE3/pythia8/generator/pythia8_pp_clr-blc2_13tev.cfg

‎MC/config/ALICE3/ini/pythia8_pp_clr-blc3_13tev.ini‎

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,9 @@
1+
[Diamond]
2+
width[2]=6.0
3+
4+
[GeneratorExternal]
5+
fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/ALICE3/pythia8/generator_pythia8_ALICE3.C
6+
funcName=generator_pythia8_ALICE3()
7+
8+
[GeneratorPythia8]
9+
config=${O2DPG_MC_CONFIG_ROOT}/MC/config/ALICE3/pythia8/generator/pythia8_pp_clr-blc3_13tev.cfg

‎MC/config/ALICE3/ini/tests/pythia8_pp_clr-blc0_13tev.C‎

Lines changed: 19 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,19 @@
1+
int External()
2+
{
3+
std::string path{"o2sim_Kine.root"};
4+
5+
TFile file(path.c_str(), "read");
6+
if (file.IsZombie()) {
7+
std::cerr << "Cannot open ROOT file " << path << std::endl;
8+
return 1;
9+
}
10+
11+
TTree* tree = (TTree*)file.Get("o2sim");
12+
13+
if (!tree) {
14+
std::cerr << "Cannot find tree o2sim in file " << path << "\n";
15+
return 1;
16+
}
17+
18+
return 0;
19+
}

‎MC/config/ALICE3/ini/tests/pythia8_pp_clr-blc2_13tev.C‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1 @@
1+
pythia8_pp_clr-blc0_13tev.C

‎MC/config/ALICE3/ini/tests/pythia8_pp_clr-blc3_13tev.C‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1 @@
1+
pythia8_pp_clr-blc0_13tev.C

‎MC/config/ALICE3/pythia8/generator/pythia8_pp_clr-blc0_13tev.cfg‎

Lines changed: 44 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,44 @@
1+
### ============================================
2+
### PYTHIA 8 — CLR-BLC Mode 0 configuration
3+
### Based on table: Monash vs CLR-BLC (Mode 0) https://link.springer.com/article/10.1007/JHEP08(2015)003 in Appendix A
4+
### ============================================
5+
6+
### Beams
7+
Beams:idA = 2212
8+
Beams:idB = 2212
9+
Beams:eCM = 13000.
10+
Beams:frameType = 1
11+
12+
### Decays (match ALICE)
13+
ParticleDecays:limitTau0 = on
14+
ParticleDecays:tau0Max = 10.0
15+
16+
### Processes and base tune
17+
SoftQCD:inelastic = on # all inelastic processes
18+
Tune:pp = 14 # Monash 2013 as baseline
19+
20+
### String fragmentation
21+
StringPT:sigma = 0.335
22+
StringZ:aLund = 0.36
23+
StringZ:bLund = 0.56
24+
StringFlav:probQQtoQ = 0.078 # probQQtoQ → Q (Mode 0,2,3)
25+
StringFlav:probStoUD = 0.20 # s/u,d suppression (Mode 0,2,3)
26+
27+
# Diquark–diquark join probabilities (vector of 4 numbers)
28+
StringFlav:probQQ1toQQ0join = 0.0275,0.0275,0.0275,0.0275 # (Mode 0,2,3)
29+
### Multi-parton interactions
30+
MultiPartonInteractions:pT0Ref = 2.12 # Mode 0 pT0Ref
31+
32+
### Beam remnants
33+
BeamRemnants:remnantMode = 1 # required for CR mode 1
34+
BeamRemnants:saturation = 5
35+
36+
### Colour reconnection: CLR-BLC Mode 0
37+
ColourReconnection:mode = 1 # new CR scheme
38+
ColourReconnection:allowDoubleJunRem = off # as in table
39+
ColourReconnection:m0 = 2.9 # Mode 0 mass scale
40+
ColourReconnection:allowJunctions = on
41+
ColourReconnection:junctionCorrection = 1.43 # Mode 0
42+
ColourReconnection:timeDilationMode = 0 # Mode 0
43+
44+
Random:setSeed = on

‎MC/config/ALICE3/pythia8/generator/pythia8_pp_clr-blc2_13tev.cfg‎

Lines changed: 45 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,45 @@
1+
### ============================================
2+
### PYTHIA 8 — CLR-BLC Mode 2 configuration
3+
### Based on table: Monash vs CLR-BLC (Mode 2) https://link.springer.com/article/10.1007/JHEP08(2015)003 in Appendix A
4+
### ============================================
5+
6+
### Beams
7+
Beams:idA = 2212
8+
Beams:idB = 2212
9+
Beams:eCM = 13000.
10+
Beams:frameType = 1
11+
12+
### Decays (match ALICE)
13+
ParticleDecays:limitTau0 = on
14+
ParticleDecays:tau0Max = 10.0
15+
16+
### Processes and base tune
17+
SoftQCD:inelastic = on # all inelastic processes
18+
Tune:pp = 14 # Monash 2013 as baseline
19+
20+
### String fragmentation
21+
StringPT:sigma = 0.335
22+
StringZ:aLund = 0.36
23+
StringZ:bLund = 0.56
24+
StringFlav:probQQtoQ = 0.078 # probQQtoQ → Q (Mode 0,2,3)
25+
StringFlav:probStoUD = 0.20 # s/u,d suppression (Mode 0,2,3)
26+
27+
# Diquark–diquark join probabilities (vector of 4 numbers)
28+
StringFlav:probQQ1toQQ0join = 0.0275,0.0275,0.0275,0.0275 # (Mode 0,2,3)
29+
### Multi-parton interactions
30+
MultiPartonInteractions:pT0Ref = 2.15 # Mode 2 pT0Ref
31+
32+
### Beam remnants
33+
BeamRemnants:remnantMode = 1 # required for CR mode 1
34+
BeamRemnants:saturation = 5
35+
36+
### Colour reconnection: CLR-BLC Mode 2
37+
ColourReconnection:mode = 1 # new CR scheme
38+
ColourReconnection:allowDoubleJunRem = off # as in table
39+
ColourReconnection:m0 = 0.3 # Mode 2 mass scale
40+
ColourReconnection:allowJunctions = on
41+
ColourReconnection:junctionCorrection = 1.2 # Mode 2
42+
ColourReconnection:timeDilationMode = 2 # Mode 2
43+
ColourReconnection:timeDilationPar = 0.18 # Mode 2
44+
45+
Random:setSeed = on

‎MC/config/ALICE3/pythia8/generator/pythia8_pp_clr-blc3_13tev.cfg‎

Lines changed: 45 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,45 @@
1+
### ============================================
2+
### PYTHIA 8 — CLR-BLC Mode 3 configuration
3+
### Based on table: Monash vs CLR-BLC (Mode 3) https://link.springer.com/article/10.1007/JHEP08(2015)003 in Appendix A
4+
### ============================================
5+
6+
### Beams
7+
Beams:idA = 2212
8+
Beams:idB = 2212
9+
Beams:eCM = 13000.
10+
Beams:frameType = 1
11+
12+
### Decays (match ALICE)
13+
ParticleDecays:limitTau0 = on
14+
ParticleDecays:tau0Max = 10.0
15+
16+
### Processes and base tune
17+
SoftQCD:inelastic = on # all inelastic processes
18+
Tune:pp = 14 # Monash 2013 as baseline
19+
20+
### String fragmentation
21+
StringPT:sigma = 0.335
22+
StringZ:aLund = 0.36
23+
StringZ:bLund = 0.56
24+
StringFlav:probQQtoQ = 0.078 # probQQtoQ → Q (Mode 0,2,3)
25+
StringFlav:probStoUD = 0.20 # s/u,d suppression (Mode 0,2,3)
26+
27+
# Diquark–diquark join probabilities (vector of 4 numbers)
28+
StringFlav:probQQ1toQQ0join = 0.0275,0.0275,0.0275,0.0275 # (Mode 0,2,3)
29+
### Multi-parton interactions
30+
MultiPartonInteractions:pT0Ref = 2.05 # Mode 3 pT0Ref
31+
32+
### Beam remnants
33+
BeamRemnants:remnantMode = 1 # required for CR mode 1
34+
BeamRemnants:saturation = 5
35+
36+
### Colour reconnection: CLR-BLC Mode 3
37+
ColourReconnection:mode = 1 # new CR scheme
38+
ColourReconnection:allowDoubleJunRem = off # as in table
39+
ColourReconnection:m0 = 0.3 # Mode 3 mass scale
40+
ColourReconnection:allowJunctions = on
41+
ColourReconnection:junctionCorrection = 1.15 # Mode 3
42+
ColourReconnection:timeDilationMode = 3 # Mode 3
43+
ColourReconnection:timeDilationPar = 0.073 # Mode 3
44+
45+
Random:setSeed = on

0 commit comments

Comments
 (0)