Skip to content

Commit 8e2a52d

Browse files
ddobrigkddobrigk
committed
Common: add light ion contamination systems
1 parent ed738ad commit 8e2a52d

File tree

12 files changed

+196
-0
lines changed

12 files changed

+196
-0
lines changed

‎MC/config/common/ini/pythia8_HeO_536.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/common/pythia8/generator/pythia8_HeO_536.cfg

‎MC/config/common/ini/pythia8_NeO_536.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/common/pythia8/generator/pythia8_NeO_536.cfg

‎MC/config/common/ini/pythia8_OHe_536.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/common/pythia8/generator/pythia8_OHe_536.cfg

‎MC/config/common/ini/pythia8_ONe_536.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/common/pythia8/generator/pythia8_ONe_536.cfg

‎MC/config/common/ini/tests/pythia8_HeO_536.C‎

Lines changed: 24 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,24 @@
1+
int External() {
2+
std::string path{"o2sim_Kine.root"};
3+
4+
TFile file(path.c_str(), "READ");
5+
if (file.IsZombie()) {
6+
std::cerr << "Cannot open ROOT file " << path << "\n";
7+
return 1;
8+
}
9+
10+
auto tree = (TTree *)file.Get("o2sim");
11+
if (!tree) {
12+
std::cerr << "Cannot find tree o2sim in file " << path << "\n";
13+
return 1;
14+
}
15+
std::vector<o2::MCTrack> *tracks{};
16+
tree->SetBranchAddress("MCTrack", &tracks);
17+
18+
auto nEvents = tree->GetEntries();
19+
if (nEvents == 0) {
20+
std::cerr << "No event of interest\n";
21+
return 1;
22+
}
23+
return 0;
24+
}

‎MC/config/common/ini/tests/pythia8_NeO_536.C‎

Lines changed: 24 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,24 @@
1+
int External() {
2+
std::string path{"o2sim_Kine.root"};
3+
4+
TFile file(path.c_str(), "READ");
5+
if (file.IsZombie()) {
6+
std::cerr << "Cannot open ROOT file " << path << "\n";
7+
return 1;
8+
}
9+
10+
auto tree = (TTree *)file.Get("o2sim");
11+
if (!tree) {
12+
std::cerr << "Cannot find tree o2sim in file " << path << "\n";
13+
return 1;
14+
}
15+
std::vector<o2::MCTrack> *tracks{};
16+
tree->SetBranchAddress("MCTrack", &tracks);
17+
18+
auto nEvents = tree->GetEntries();
19+
if (nEvents == 0) {
20+
std::cerr << "No event of interest\n";
21+
return 1;
22+
}
23+
return 0;
24+
}

‎MC/config/common/ini/tests/pythia8_OHe_536.C‎

Lines changed: 24 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,24 @@
1+
int External() {
2+
std::string path{"o2sim_Kine.root"};
3+
4+
TFile file(path.c_str(), "READ");
5+
if (file.IsZombie()) {
6+
std::cerr << "Cannot open ROOT file " << path << "\n";
7+
return 1;
8+
}
9+
10+
auto tree = (TTree *)file.Get("o2sim");
11+
if (!tree) {
12+
std::cerr << "Cannot find tree o2sim in file " << path << "\n";
13+
return 1;
14+
}
15+
std::vector<o2::MCTrack> *tracks{};
16+
tree->SetBranchAddress("MCTrack", &tracks);
17+
18+
auto nEvents = tree->GetEntries();
19+
if (nEvents == 0) {
20+
std::cerr << "No event of interest\n";
21+
return 1;
22+
}
23+
return 0;
24+
}

‎MC/config/common/ini/tests/pythia8_ONe_536.C‎

Lines changed: 24 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,24 @@
1+
int External() {
2+
std::string path{"o2sim_Kine.root"};
3+
4+
TFile file(path.c_str(), "READ");
5+
if (file.IsZombie()) {
6+
std::cerr << "Cannot open ROOT file " << path << "\n";
7+
return 1;
8+
}
9+
10+
auto tree = (TTree *)file.Get("o2sim");
11+
if (!tree) {
12+
std::cerr << "Cannot find tree o2sim in file " << path << "\n";
13+
return 1;
14+
}
15+
std::vector<o2::MCTrack> *tracks{};
16+
tree->SetBranchAddress("MCTrack", &tracks);
17+
18+
auto nEvents = tree->GetEntries();
19+
if (nEvents == 0) {
20+
std::cerr << "No event of interest\n";
21+
return 1;
22+
}
23+
return 0;
24+
}

‎MC/config/common/pythia8/generator/pythia8_HeO.cfg‎

Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,16 @@
1+
### OO beams
2+
Beams:idA = 1000020040 # Alpha
3+
Beams:idB = 1000080160 # Oxygen
4+
Beams:frameType 2 # back-to-back beams of different energies and particles
5+
Beams:eA 2680. # 50% of 5360
6+
Beams:eB 2680. # 50% of 5360
7+
8+
ParticleDecays:limitTau0 = on
9+
ParticleDecays:tau0Max = 10. ### match alice: 1cm/c = 10.0mm/c
10+
11+
### Save some CPU at init of jobs
12+
### To avoid refitting, add the following lines to your configuration file:
13+
### HeavyIon:SigFitNGen = 0
14+
### HeavyIon:SigFitDefPar = 2.15,18.42,0.33
15+
16+
Random:setSeed = on

‎MC/config/common/pythia8/generator/pythia8_NeO.cfg‎

Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,16 @@
1+
### OO beams
2+
Beams:idA = 1000100200 # Neon
3+
Beams:idB = 1000080160 # Oxygen
4+
Beams:frameType 2 # back-to-back beams of different energies and particles
5+
Beams:eA 2680. # 50% of 5360
6+
Beams:eB 2680. # 50% of 5360
7+
8+
ParticleDecays:limitTau0 = on
9+
ParticleDecays:tau0Max = 10. ### match alice: 1cm/c = 10.0mm/c
10+
11+
### Save some CPU at init of jobs
12+
### To avoid refitting, add the following lines to your configuration file:
13+
### HeavyIon:SigFitNGen = 0
14+
### HeavyIon:SigFitDefPar = 2.15,18.42,0.33
15+
16+
Random:setSeed = on

0 commit comments

Comments
 (0)