Skip to content

Commit 239cd66

Browse files
jesgumjesgum
authored
ALICE 3: Add xicc injector for pp and PbPb (#1688)
* PbPb and pp configurations with xicc gun * Adjust test files * Update generators * Update generators * Update generators * Update generators * Update generators * Edited misprint * Update generators * Update generator
1 parent 26bb981 commit 239cd66

File tree

9 files changed

+1139
-0
lines changed

9 files changed

+1139
-0
lines changed

‎MC/config/ALICE3/ini/tests/xicc_PbPb.C‎

Lines changed: 26 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,26 @@
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 << "\n";
8+
return 1;
9+
}
10+
11+
auto tree = (TTree *)file.Get("o2sim");
12+
if (!tree) {
13+
std::cerr << "Cannot find tree o2sim in file " << path << "\n";
14+
return 1;
15+
}
16+
std::vector<o2::MCTrack> *tracks{};
17+
tree->SetBranchAddress("MCTrack", &tracks);
18+
19+
auto nEvents = tree->GetEntries();
20+
auto nSelected = tree->Scan("MCTrack.GetPdgCode()", "MCTrack.GetPdgCode() == 4422");
21+
if (nSelected == 0) {
22+
std::cerr << "No event of interest\n";
23+
return 1;
24+
}
25+
return 0;
26+
}

‎MC/config/ALICE3/ini/tests/xicc_pp.C‎

Lines changed: 26 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,26 @@
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 << "\n";
8+
return 1;
9+
}
10+
11+
auto tree = (TTree *)file.Get("o2sim");
12+
if (!tree) {
13+
std::cerr << "Cannot find tree o2sim in file " << path << "\n";
14+
return 1;
15+
}
16+
std::vector<o2::MCTrack> *tracks{};
17+
tree->SetBranchAddress("MCTrack", &tracks);
18+
19+
auto nEvents = tree->GetEntries();
20+
auto nSelected = tree->Scan("MCTrack.GetPdgCode()", "MCTrack.GetPdgCode() == 4422");
21+
if (nSelected == 0) {
22+
std::cerr << "No event of interest\n";
23+
return 1;
24+
}
25+
return 0;
26+
}

‎MC/config/ALICE3/ini/xicc_PbPb.ini‎

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
[GeneratorExternal]
2+
fileName=${O2DPG_ROOT}/MC/config/ALICE3/pythia8/generator_pythia8_gun_PbPb.C
3+
funcName=generateNativeXiCC()
4+
5+
[GeneratorPythia8]
6+
config=${O2DPG_ROOT}/MC/config/ALICE3/pythia8/generator/config_custom_xicc.cfg

‎MC/config/ALICE3/ini/xicc_pp.ini‎

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
[GeneratorExternal]
2+
fileName=${O2DPG_ROOT}/MC/config/ALICE3/pythia8/generator_pythia8_gun_pp.C
3+
funcName=generateNativeXiCC()
4+
5+
[GeneratorPythia8]
6+
config=${O2DPG_ROOT}/MC/config/ALICE3/pythia8/generator/config_custom_xicc.cfg

‎MC/config/ALICE3/pythia8/generator/config_custom_xicc.cfg‎

Lines changed: 51 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,51 @@
1+
### Xi_c
2+
#4232:all = Xic+ Xic- 2 3 0 2.46793 0. 0. 0. 10
3+
4422:all = XiCC AntiXiCC 2 6 0 3.6212 0. 0. 0. 10
4+
5+
### Omega_c
6+
#4332:all = Omega_c0 Omega_cbar0 2 0 0 2.69520 0. 0. 0. 10
7+
4432:all = Omega_cc Omega_ccbar 2 3 0 3.746 0. 0. 0. 10
8+
4444:all = Omega_ccc Omega_cccbar 4 6 0 4.797 0. 0. 0. 10
9+
10+
ProcessLevel:all = off
11+
12+
### changing the ctau value in mm/c
13+
### use scaling based on Lambda, Xi, Omega progression
14+
4332:tau0=0.08000000000
15+
4432:tau0=0.04980000000
16+
4444:tau0=0.02496000000
17+
4422:tau0=0.07675
18+
4232:tau0=0.132
19+
20+
### add OmegaC decay absent in PYTHIA8 decay table
21+
4332:addChannel = 1 1 0 3334 211
22+
4432:addChannel = 1 1 0 4332 211
23+
4444:oneChannel = 1 1 0 4432 211
24+
25+
4422:addChannel = 1 1 0 4232 211
26+
4422:onMode = off
27+
4422:onIfMatch = 4232 211
28+
29+
### add OmegaC decay absent in PYTHIA8 decay table
30+
4232:addChannel = 1 0.0001 0 3312 211 211
31+
4232:onMode = off
32+
4232:onIfMatch = 3312 211 211
33+
34+
### force the OmegaC to decay in the Omega_c -> Omega pi channel
35+
4332:onMode = off
36+
4332:onIfMatch = 3334 211
37+
4432:onMode = off
38+
4432:onIfMatch = 4332 211
39+
4444:onMode = off
40+
4444:onIfMatch = 4432 211
41+
42+
#boundless decays
43+
ParticleDecays:limitTau0 off
44+
45+
#do other decays via geant3, etc
46+
310:onMode = off
47+
3122:onMode = off
48+
3112:onMode = off
49+
3222:onMode = off
50+
3312:onMode = off
51+
3334:onMode = off

‎MC/config/ALICE3/pythia8/generator/pythia8_hi.cmnd‎

Lines changed: 21 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,21 @@
1+
### beams
2+
Beams:idA 1000822080 # Pb
3+
Beams:idB 1000822080 # Pb
4+
Beams:eCM 5520. # GeV
5+
6+
### heavy-ion settings (valid for Pb-Pb 5520 only)
7+
HeavyIon:SigFitNGen = 0
8+
HeavyIon:SigFitDefPar = 13.88,1.84,0.22,0.0,0.0,0.0,0.0,0.0
9+
HeavyIon:bWidth = 14.48
10+
11+
### processes (apparently not to be defined)
12+
13+
### decays
14+
ParticleDecays:limitTau0 on
15+
ParticleDecays:tau0Max 0.001
16+
17+
! 2) Seed settings
18+
! Seed is set inside the generator
19+
! If run on the grid, seed is set to job id
20+
! If run locally, seed is set to 0
21+
Random:setSeed = on ! Random seed on

‎MC/config/ALICE3/pythia8/generator/pythia8_pp.cmnd‎

Lines changed: 61 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,61 @@
1+
! WithMPI_MPIBased.cmnd.
2+
! This file contains the More QCD scheme colour reconnection settings scenario
3+
! to be read in for a Pythia8 run.
4+
! Lines beginning with a "!" are comments.
5+
! Author: Renan Acconcia
6+
7+
!pythia default
8+
!Last update: 01/08/2018
9+
!Only Inelastic events: INEL
10+
11+
! 1) Settings related to output in init(), next() and stat().
12+
Init:showChangedSettings = on ! list changed settings
13+
!Init:showChangedParticleData = off ! list changed particle data
14+
!Next:numberCount = 0 ! print message every n events
15+
!Next:numberShowInfo = 0 ! print event information n times
16+
!Next:numberShowProcess = 1 ! print process record n times
17+
!Next:numberShowEvent = 1 ! print event record n times
18+
19+
! 2) Seed settings
20+
! Seed is set inside the generator
21+
! If run on the grid, seed is set to job id
22+
! If run locally, seed is set to 0
23+
Random:setSeed = on ! Random seed on
24+
25+
! 3) Tunes settings
26+
! Tune:preferLHAPDF = 2 ! Using LHAPDF 6
27+
StringPT:sigma = 0.335
28+
StringZ:aLund = 0.36
29+
StringZ:bLund = 0.56
30+
StringFlav:probQQtoQ = 0.078
31+
StringFlav:ProbStoUD = 0.2
32+
StringFlav:probQQ1toQQ0join = 0.0275,0.0275,0.0275,0.0275
33+
34+
MultiPartonInteractions:pT0Ref = 2.12
35+
36+
BeamRemnants:remnantMode = 1
37+
BeamRemnants:saturation = 5
38+
39+
ColourReconnection:mode = 1
40+
ColourReconnection:allowDoubleJunRem = off
41+
ColourReconnection:m0 = 2.9
42+
ColourReconnection:allowJunctions = on
43+
ColourReconnection:junctionCorrection = 1.43
44+
ColourReconnection:timeDilationMode = 0
45+
46+
! 5) Settings for the softQCD-process generation. NSD scenario.
47+
SoftQCD:nonDiffractive = on ! Non Diffractive on
48+
SoftQCD:singleDiffractive = on ! Single Diffractive on
49+
SoftQCD:doubleDiffractive = on ! Double Diffractive on
50+
SoftQCD:centralDiffractive = on ! Central Diffractive on
51+
SoftQCD:elastic = off ! Elastic off
52+
53+
! 6) Decay Supression
54+
ParticleDecays:limitTau0 = on ! Only particles with tau0 < tau0Max are decayed.
55+
ParticleDecays:tau0Max = 0.001 ! Tau0Max = 3. mm/c
56+
57+
! 7) Beam parameter settings. Values below agree with default ones.
58+
!this is default do not need to be specified
59+
Beams:idA = 2212 ! first beam, p = 2212, pbar = -2212
60+
Beams:idB = 2212 ! second beam, p = 2212, pbar = -2212
61+
Beams:eCM = 14000.0 //energy

0 commit comments

Comments
 (0)