-
Notifications
You must be signed in to change notification settings - Fork 496
Expand file tree
/
Copy pathHitProcessingManager.cxx
More file actions
221 lines (196 loc) · 6.62 KB
/
HitProcessingManager.cxx
File metadata and controls
221 lines (196 loc) · 6.62 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
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
// 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.
#include "Steer/HitProcessingManager.h"
#include <fairlogger/Logger.h>
#include <vector>
#include <map>
#include <iostream>
#include <TFile.h>
#include <TClass.h>
#include <TRandom3.h>
ClassImp(o2::steer::HitProcessingManager);
namespace o2
{
namespace steer
{
bool HitProcessingManager::setupChain()
{
if (mBackgroundFileNames.size() == 0 && mSignalFileNames.size() == 0) {
// there are no files to be analysed/processed;
LOG(warning) << "No files to be analysed";
return false;
}
// count max signal id
int sourcecounter = 0;
int maxsourceid = 0;
for (auto& pair : mSignalFileNames) {
sourcecounter++;
maxsourceid = std::max(pair.first, maxsourceid);
}
if (maxsourceid != sourcecounter) {
LOG(warning) << "max source id " << maxsourceid << " vs " << sourcecounter;
}
LOG(info) << "setting up " << maxsourceid + 1 << " chains";
mSimChains.resize(maxsourceid + 1);
// allocate chains
for (int i = 0; i < mSimChains.size(); ++i) {
// make o2sim a parameter
mSimChains[i] = new TChain("o2sim");
}
// background chain
auto& c = *mSimChains[0];
c.Reset();
for (auto& filename : mBackgroundFileNames) {
c.AddFile(o2::base::NameConf::getMCHeadersFileName(filename.data()).c_str());
}
for (auto& pair : mSignalFileNames) {
const auto& signalid = pair.first;
const auto& filenamevector = pair.second;
auto& chain = *mSimChains[signalid];
for (auto& filename : filenamevector) {
chain.AddFile(o2::base::NameConf::getMCHeadersFileName(filename.data()).c_str());
}
}
return true;
}
void HitProcessingManager::setupRun(int ncollisions)
{
if (!setupChain()) {
return;
}
if (mGeometryFile.size() > 0) {
// load geometry
TGeoManager::Import(mGeometryFile.c_str());
}
//
if (ncollisions != -1) {
mNumberOfCollisions = ncollisions;
} else {
mNumberOfCollisions = mSimChains[0]->GetEntries();
LOG(info) << "Automatic deduction of number of collisions ... will just take number of background entries "
<< mNumberOfCollisions;
}
mDigitizationContext.setNCollisions(mNumberOfCollisions);
sampleCollisionTimes();
// sample collision (background-signal) constituents
sampleCollisionConstituents();
// store prefixes as part of Context
std::vector<std::string> prefixes;
prefixes.emplace_back(mBackgroundFileNames[0]);
for (auto k : mSignalFileNames) {
prefixes.emplace_back(k.second[0]);
}
mDigitizationContext.setSimPrefixes(prefixes);
}
void HitProcessingManager::writeDigitizationContext(const char* filename) const
{
mDigitizationContext.saveToFile(filename);
}
bool HitProcessingManager::setupRunFromExistingContext(const char* filename)
{
auto context = DigitizationContext::loadFromFile(filename);
if (context) {
context->printCollisionSummary();
mDigitizationContext = *context;
return true;
}
LOG(warn) << "NO DIGITIZATIONCONTEXT FOUND";
return false;
}
void HitProcessingManager::sampleCollisionTimes()
{
mDigitizationContext.getEventRecords().resize(mDigitizationContext.getNCollisions());
mInteractionSampler.generateCollisionTimes(mDigitizationContext.getEventRecords());
mDigitizationContext.setBunchFilling(mInteractionSampler.getBunchFilling());
mDigitizationContext.setMuPerBC(mInteractionSampler.getMuPerBC());
}
void HitProcessingManager::sampleCollisionConstituents()
{
TRandom3 rnd(0); // we don't use the global to be in isolation
auto getBackgroundRoundRobin = [this]() {
static int bgcounter = 0;
int numbg = mSimChains[0]->GetEntries();
if (bgcounter == numbg) {
bgcounter = 0;
}
return EventPart(0, bgcounter++);
};
const int nsignalids = mSimChains.size() - 1;
auto getSignalRoundRobin = [this, nsignalids]() {
static int bgcounter = 0;
static int signalid = 0;
static std::vector<int> counter(nsignalids, 0);
if (signalid == nsignalids) {
signalid = 0;
}
const auto realsourceid = signalid + 1;
int numentries = mSimChains[realsourceid]->GetEntries();
if (counter[signalid] == numentries) {
counter[signalid] = 0;
}
EventPart e(realsourceid, counter[signalid]);
counter[signalid]++;
signalid++;
return e;
};
auto getRandomBackground = [this, &rnd]() {
int numbg = mSimChains[0]->GetEntries();
const auto eventID = (int)numbg * rnd.Rndm();
return EventPart(0, eventID);
};
auto getRandomSignal = [this, nsignalids, &rnd]() {
const auto sourceID = 1 + (int)(rnd.Rndm() * nsignalids);
const auto signalID = (int)(rnd.Rndm() * mSimChains[sourceID]->GetEntries());
return EventPart(sourceID, signalID);
};
// we fill mDigitizationContext.mEventParts
auto& eventparts = mDigitizationContext.getEventParts();
eventparts.clear();
eventparts.resize(mDigitizationContext.getEventRecords().size());
for (int i = 0; i < mDigitizationContext.getEventRecords().size(); ++i) {
eventparts[i].clear();
// NOTE: THIS PART WOULD BENEFIT FROM A MAJOR REDESIGN
// WISHFUL ITEMS WOULD BE:
// - ALLOW COLLISION ENGINEERING FROM OUTSIDE (give wanted sequence as file)
// * the outside person can decide what kind of sampling and sequence to use
// - CHECK IF VERTEX IS CONSISTENT
if (mSampleCollisionsRandomly) {
eventparts[i].emplace_back(getRandomBackground());
if (mSimChains.size() > 1) {
eventparts[i].emplace_back(getRandomSignal());
}
} else {
// push any number of constituents?
// for the moment just 2 : one background and one signal
eventparts[i].emplace_back(getBackgroundRoundRobin());
if (mSimChains.size() > 1) {
eventparts[i].emplace_back(getSignalRoundRobin());
}
}
}
// push any number of constituents?
// for the moment just max 2 : one background and one signal
mDigitizationContext.setMaxNumberParts(1);
if (mSimChains.size() > 1) {
mDigitizationContext.setMaxNumberParts(2);
}
mDigitizationContext.printCollisionSummary();
}
void HitProcessingManager::run()
{
setupRun();
// sample other stuff
for (auto& f : mRegisteredRunFunctions) {
f(mDigitizationContext);
}
}
} // end namespace steer
} // end namespace o2