Skip to content

Commit 22da986

Browse files
authored
Add LF generator for MC with coalescence production for De Tr He3 (#2207)
* add generators for nuclei MC with coalescence De Tr He3 * Add LF generator for MC with coalescence production for De Tr He3 * rm old cfg
1 parent 35954dd commit 22da986

File tree

3 files changed

+53
-0
lines changed

3 files changed

+53
-0
lines changed
Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,7 @@
1+
[GeneratorExternal]
2+
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGLF/pythia8/generator_pythia8_coalescence.C
3+
funcName = generateCoalescence({1000010020, 1000010030, 1000020030}, 1, 0.239)
4+
5+
[GeneratorPythia8]
6+
config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGLF/pythia8/generator/pythia8_inel_pp536tev.cfg
7+
Lines changed: 28 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,28 @@
1+
int External()
2+
{
3+
std::string path{"o2sim_Kine.root"};
4+
5+
TFile file(path.c_str(), "READ");
6+
if (file.IsZombie())
7+
{
8+
std::cerr << "Cannot open ROOT file " << path << "\n";
9+
return 1;
10+
}
11+
12+
auto tree = (TTree *)file.Get("o2sim");
13+
if (!tree)
14+
{
15+
std::cerr << "Cannot find tree o2sim in file " << path << "\n";
16+
return 1;
17+
}
18+
std::vector<o2::MCTrack> *tracks{};
19+
tree->SetBranchAddress("MCTrack", &tracks);
20+
21+
auto nEvents = tree->GetEntries();
22+
if (nEvents < 1)
23+
{
24+
std::cerr << "No events actually generated: not OK!";
25+
return 1;
26+
}
27+
return 0;
28+
}
Lines changed: 18 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,18 @@
1+
### beams
2+
Beams:idA = 2212 # proton
3+
Beams:idB = 2212 # proton
4+
Beams:eCM = 5360. # GeV
5+
6+
### processes
7+
SoftQCD:inelastic = on # all inelastic processes
8+
9+
### decays
10+
ParticleDecays:limitTau0 = on
11+
ParticleDecays:tau0Max = 10.
12+
13+
### phase space cuts
14+
PhaseSpace:pTHatMin = 0.000000
15+
PhaseSpace:pTHatMax = -1.000000
16+
17+
Random:setSeed = on
18+
Random:seed = 0

0 commit comments

Comments
 (0)