-
Notifications
You must be signed in to change notification settings - Fork 496
Expand file tree
/
Copy pathGenCosmics.C
More file actions
81 lines (72 loc) · 2.54 KB
/
GenCosmics.C
File metadata and controls
81 lines (72 loc) · 2.54 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
//< Macro to run QED background generator, us it as e.g.
//< o2-sim -n10000 -m PIPE ITS TPC -g extgen --configKeyValues "GeneratorExternal.fileName=$O2_ROOT/share/Generators/external/GenCosmicsLoader.C"
R__LOAD_LIBRARY(libGeneratorCosmics)
using namespace o2::eventgen;
o2::eventgen::GeneratorTGenerator* GenCosmics()
{
auto genCosm = new GeneratorCosmics();
auto& cosmParam = GenCosmicsParam::Instance();
if (cosmParam.param == GenCosmicsParam::ParamMI) {
genCosm->setParamMI();
} else if (cosmParam.param == GenCosmicsParam::ParamACORDE) {
genCosm->setParamACORDE();
} else if (cosmParam.param == GenCosmicsParam::ParamTPC) {
genCosm->setParamTPC();
} else {
LOG(fatal) << "Unknown cosmics param type " << cosmParam.param;
}
genCosm->setNPart(cosmParam.nPart);
genCosm->setPRange(cosmParam.pmin, cosmParam.pmax);
switch (cosmParam.accept) {
case GenCosmicsParam::ITS0:
genCosm->requireITS0();
break;
case GenCosmicsParam::ITS1:
genCosm->requireITS1();
break;
case GenCosmicsParam::ITS2:
genCosm->requireITS2();
break;
case GenCosmicsParam::ITS3:
genCosm->requireITS3();
break;
case GenCosmicsParam::ITS4:
genCosm->requireITS4();
break;
case GenCosmicsParam::ITS5:
genCosm->requireITS5();
break;
case GenCosmicsParam::ITS6:
genCosm->requireITS6();
break;
case GenCosmicsParam::TPC:
genCosm->requireTPC();
break;
case GenCosmicsParam::Custom:
genCosm->requireXZAccepted(cosmParam.customAccX, cosmParam.customAccZ);
break;
default:
LOG(fatal) << "Unknown cosmics acceptance type " << cosmParam.accept;
break;
}
genCosm->Init();
// instance and configure TGenerator interface
auto tgen = new o2::eventgen::GeneratorTGenerator();
tgen->setMomentumUnit(1.); // [GeV/c]
tgen->setEnergyUnit(1.); // [GeV/c]
tgen->setPositionUnit(1.); // [cm]
tgen->setTimeUnit(1.); // [s]
tgen->setTGenerator(genCosm);
fg = tgen;
return tgen;
}