Skip to content

Commit 4d53cc8

Browse files
dsekihatdsekihat
committed
MC/PWGEM: add cfg and init for vm2ll in pp/PbPb
1 parent 3cd6b38 commit 4d53cc8

File tree

9 files changed

+267
-0
lines changed

9 files changed

+267
-0
lines changed

‎MC/config/PWGEM/ini/pythia8_PbPb_5360_VM2ll.ini‎

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

‎MC/config/PWGEM/ini/pythia8_pp_13600_VM2ll.ini‎

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

‎MC/config/PWGEM/ini/pythia8_pp_5360_VM2ll.ini‎

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

‎MC/config/PWGEM/ini/tests/pythia8_PbPb_5360_VM2ll.C‎

Lines changed: 56 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,56 @@
1+
int External() {
2+
std::string path{"o2sim_Kine.root"};
3+
// Check that file exists, can be opened and has the correct tree
4+
TFile file(path.c_str(), "READ");
5+
if (file.IsZombie())
6+
{
7+
std::cerr << "Cannot open ROOT file " << path << "\n";
8+
return 1;
9+
}
10+
auto tree = (TTree *)file.Get("o2sim");
11+
if (!tree)
12+
{
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+
// Check if all events are filled
20+
auto nEvents = tree->GetEntries();
21+
for (Long64_t i = 0; i < nEvents; ++i)
22+
{
23+
tree->GetEntry(i);
24+
if (tracks->empty())
25+
{
26+
std::cerr << "Empty entry found at event " << i << "\n";
27+
return 1;
28+
}
29+
}
30+
// check if each event has at least two oxygen ions
31+
for (int i = 0; i < nEvents; i++)
32+
{
33+
auto check = tree->GetEntry(i);
34+
int count = 0;
35+
for (int idxMCTrack = 0; idxMCTrack < tracks->size(); ++idxMCTrack)
36+
{
37+
auto track = tracks->at(idxMCTrack);
38+
if (track.GetPdgCode() == 1000080160)
39+
{
40+
count++;
41+
}
42+
}
43+
if (count < 2)
44+
{
45+
std::cerr << "Event " << i << " has less than 2 oxygen ions\n";
46+
return 1;
47+
}
48+
}
49+
50+
return 0;
51+
}
52+
53+
int pythia8()
54+
{
55+
return External();
56+
}

‎MC/config/PWGEM/ini/tests/pythia8_pp_13600_VM2ll.C‎

Lines changed: 56 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,56 @@
1+
int External() {
2+
std::string path{"o2sim_Kine.root"};
3+
// Check that file exists, can be opened and has the correct tree
4+
TFile file(path.c_str(), "READ");
5+
if (file.IsZombie())
6+
{
7+
std::cerr << "Cannot open ROOT file " << path << "\n";
8+
return 1;
9+
}
10+
auto tree = (TTree *)file.Get("o2sim");
11+
if (!tree)
12+
{
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+
// Check if all events are filled
20+
auto nEvents = tree->GetEntries();
21+
for (Long64_t i = 0; i < nEvents; ++i)
22+
{
23+
tree->GetEntry(i);
24+
if (tracks->empty())
25+
{
26+
std::cerr << "Empty entry found at event " << i << "\n";
27+
return 1;
28+
}
29+
}
30+
// check if each event has at least two oxygen ions
31+
for (int i = 0; i < nEvents; i++)
32+
{
33+
auto check = tree->GetEntry(i);
34+
int count = 0;
35+
for (int idxMCTrack = 0; idxMCTrack < tracks->size(); ++idxMCTrack)
36+
{
37+
auto track = tracks->at(idxMCTrack);
38+
if (track.GetPdgCode() == 1000080160)
39+
{
40+
count++;
41+
}
42+
}
43+
if (count < 2)
44+
{
45+
std::cerr << "Event " << i << " has less than 2 oxygen ions\n";
46+
return 1;
47+
}
48+
}
49+
50+
return 0;
51+
}
52+
53+
int pythia8()
54+
{
55+
return External();
56+
}

‎MC/config/PWGEM/ini/tests/pythia8_pp_5360_VM2ll.C‎

Lines changed: 56 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,56 @@
1+
int External() {
2+
std::string path{"o2sim_Kine.root"};
3+
// Check that file exists, can be opened and has the correct tree
4+
TFile file(path.c_str(), "READ");
5+
if (file.IsZombie())
6+
{
7+
std::cerr << "Cannot open ROOT file " << path << "\n";
8+
return 1;
9+
}
10+
auto tree = (TTree *)file.Get("o2sim");
11+
if (!tree)
12+
{
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+
// Check if all events are filled
20+
auto nEvents = tree->GetEntries();
21+
for (Long64_t i = 0; i < nEvents; ++i)
22+
{
23+
tree->GetEntry(i);
24+
if (tracks->empty())
25+
{
26+
std::cerr << "Empty entry found at event " << i << "\n";
27+
return 1;
28+
}
29+
}
30+
// check if each event has at least two oxygen ions
31+
for (int i = 0; i < nEvents; i++)
32+
{
33+
auto check = tree->GetEntry(i);
34+
int count = 0;
35+
for (int idxMCTrack = 0; idxMCTrack < tracks->size(); ++idxMCTrack)
36+
{
37+
auto track = tracks->at(idxMCTrack);
38+
if (track.GetPdgCode() == 1000080160)
39+
{
40+
count++;
41+
}
42+
}
43+
if (count < 2)
44+
{
45+
std::cerr << "Event " << i << " has less than 2 oxygen ions\n";
46+
return 1;
47+
}
48+
}
49+
50+
return 0;
51+
}
52+
53+
int pythia8()
54+
{
55+
return External();
56+
}

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

Lines changed: 23 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,23 @@
1+
### Specify beams
2+
Beams:idA = 1000822080
3+
Beams:idB = 1000822080
4+
Beams:eCM = 5360.0 ### energy
5+
6+
Beams:frameType = 1
7+
ParticleDecays:limitTau0 = on
8+
ParticleDecays:tau0Max = 10. ### match alice: 1cm/c = 10.0mm/c
9+
10+
### Initialize the Angantyr model to fit the total and semi-includive
11+
### cross sections in Pythia within some tolerance.
12+
HeavyIon:SigFitErr = 0.02,0.02,0.1,0.05,0.05,0.0,0.1,0.0
13+
14+
### These parameters are typicall suitable for sqrt(S_NN)=5TeV
15+
HeavyIon:SigFitDefPar = 17.24,2.15,0.33,0.0,0.0,0.0,0.0,0.0
16+
17+
Random:setSeed = on
18+
19+
# change omega, phi meson's BR below
20+
223:oneChannel = 1 0.5 0 -11 11
21+
223:addChannel = 1 0.5 0 -13 13
22+
333:oneChannel = 1 0.5 0 -11 11
23+
333:addChannel = 1 0.5 0 -13 13

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

Lines changed: 23 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,23 @@
1+
### Specify beams
2+
Beams:idA = 2212
3+
Beams:idB = 2212
4+
Beams:eCM = 13600. ### energy
5+
6+
Beams:frameType = 1
7+
ParticleDecays:limitTau0 = on
8+
ParticleDecays:tau0Max = 10. ### match alice: 1cm/c = 10.0mm/c
9+
10+
### processes
11+
SoftQCD:inelastic = on # all inelastic processes
12+
13+
# default: do nothing, Monash 2013 will do its thing
14+
Tune:pp = 14
15+
16+
Random:setSeed = on
17+
18+
# change eta, omega, phi meson's BR below
19+
221:oneChannel = 1 0.5 0 -13 13
20+
223:oneChannel = 1 0.5 0 -11 11
21+
223:addChannel = 1 0.5 0 -13 13
22+
333:oneChannel = 1 0.5 0 -11 11
23+
333:addChannel = 1 0.5 0 -13 13

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

Lines changed: 23 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,23 @@
1+
### Specify beams
2+
Beams:idA = 2212
3+
Beams:idB = 2212
4+
Beams:eCM = 5360. ### energy
5+
6+
Beams:frameType = 1
7+
ParticleDecays:limitTau0 = on
8+
ParticleDecays:tau0Max = 10. ### match alice: 1cm/c = 10.0mm/c
9+
10+
### processes
11+
SoftQCD:inelastic = on # all inelastic processes
12+
13+
# default: do nothing, Monash 2013 will do its thing
14+
Tune:pp = 14
15+
16+
Random:setSeed = on
17+
18+
# change eta, omega, phi meson's BR below
19+
221:oneChannel = 1 0.5 0 -13 13
20+
223:oneChannel = 1 0.5 0 -11 11
21+
223:addChannel = 1 0.5 0 -13 13
22+
333:oneChannel = 1 0.5 0 -11 11
23+
333:addChannel = 1 0.5 0 -13 13

0 commit comments

Comments
 (0)