Skip to content

Commit ee2a012

Browse files
MartijnLaarhovenMartijnLaarhoven
authored
Merge branch 'AliceO2Group:master' into master
2 parents d6fd502 + 7898348 commit ee2a012

23 files changed

+2962
-1832
lines changed

‎ALICE3/Core/FastTracker.cxx‎

Lines changed: 55 additions & 46 deletions
Original file line numberDiff line numberDiff line change
@@ -236,6 +236,42 @@ void FastTracker::AddTPC(float phiResMean, float zResMean)
236236
}
237237
}
238238

239+
std::map<std::string, std::map<std::string, std::string>> FastTracker::parseTEnvConfiguration(std::string filename)
240+
{
241+
std::map<std::string, std::map<std::string, std::string>> configMap;
242+
243+
TEnv env(filename.c_str());
244+
THashList* table = env.GetTable();
245+
std::vector<std::string> layers;
246+
for (int i = 0; i < table->GetEntries(); ++i) {
247+
const std::string key = table->At(i)->GetName();
248+
// key should contain exactly one dot
249+
if (key.find('.') == std::string::npos || key.find('.') != key.rfind('.')) {
250+
LOG(fatal) << "Key " << key << " does not contain exactly one dot";
251+
continue;
252+
}
253+
const std::string firstPart = key.substr(0, key.find('.'));
254+
if (std::find(layers.begin(), layers.end(), firstPart) == layers.end()) {
255+
layers.push_back(firstPart);
256+
}
257+
}
258+
env.Print();
259+
260+
// Layers
261+
for (const auto& layer : layers) {
262+
LOG(info) << " Reading layer " << layer;
263+
for (int i = 0; i < table->GetEntries(); ++i) {
264+
const std::string key = table->At(i)->GetName();
265+
if (key.find(layer + ".") == 0) {
266+
const std::string paramName = key.substr(key.find('.') + 1);
267+
const std::string value = env.GetValue(key.c_str(), "");
268+
configMap[layer][paramName] = value;
269+
}
270+
}
271+
}
272+
return configMap;
273+
}
274+
239275
void FastTracker::AddGenericDetector(std::string filename, o2::ccdb::BasicCCDBManager* ccdbManager)
240276
{
241277
LOG(info) << " Adding generic detector from file " << filename;
@@ -263,54 +299,27 @@ void FastTracker::AddGenericDetector(std::string filename, o2::ccdb::BasicCCDBMa
263299
return;
264300
}
265301

266-
TEnv env(filename.c_str());
267-
THashList* table = env.GetTable();
268-
std::vector<std::string> layers;
269-
for (int i = 0; i < table->GetEntries(); ++i) {
270-
const std::string key = table->At(i)->GetName();
271-
// key should contain exactly one dot
272-
if (key.find('.') == std::string::npos || key.find('.') != key.rfind('.')) {
273-
LOG(fatal) << "Key " << key << " does not contain exactly one dot";
274-
continue;
275-
}
276-
const std::string firstPart = key.substr(0, key.find('.'));
277-
if (std::find(layers.begin(), layers.end(), firstPart) == layers.end()) {
278-
layers.push_back(firstPart);
279-
}
280-
}
281-
// env.Print();
302+
std::map<std::string, std::map<std::string, std::string>> configMap = parseTEnvConfiguration(filename);
282303
// Layers
283-
for (const auto& layer : layers) {
284-
LOG(info) << " Reading layer " << layer;
285-
286-
auto getKey = [&layer, &env](const std::string& name, const bool required = true) {
287-
std::string key = layer + "." + name;
288-
if (!env.Defined(key.c_str())) {
289-
if (required) {
290-
LOG(fatal) << "Key " << key << " not defined in configuration file";
291-
}
292-
LOG(debug) << "Key " << key << " not defined in configuration file, getting the default value";
293-
}
294-
LOG(debug) << " Getting key " << key << " from configuration file";
295-
return key;
296-
};
297-
const float r = env.GetValue(getKey("r").c_str(), -1.0f);
298-
LOG(info) << " Layer " << layer << " has radius " << r;
299-
const float z = env.GetValue(getKey("z").c_str(), -1.0f);
300-
const float x0 = env.GetValue(getKey("x0").c_str(), 0.0f);
301-
const float xrho = env.GetValue(getKey("xrho").c_str(), 0.0f);
302-
const float resRPhi = env.GetValue(getKey("resRPhi").c_str(), 0.0f);
303-
const float resZ = env.GetValue(getKey("resZ").c_str(), 0.0f);
304-
const float eff = env.GetValue(getKey("eff").c_str(), 0.0f);
305-
const int type = env.GetValue(getKey("type").c_str(), 0);
306-
const char* deadPhiRegions = env.GetValue(getKey("deadPhiRegions", false).c_str(), "");
304+
for (const auto& layer : configMap) {
305+
LOG(info) << " Reading layer " << layer.first;
306+
const float r = std::stof(layer.second.at("r"));
307+
LOG(info) << " Layer " << layer.first << " has radius " << r;
308+
const float z = std::stof(layer.second.at("z"));
309+
const float x0 = std::stof(layer.second.at("x0"));
310+
const float xrho = std::stof(layer.second.at("xrho"));
311+
const float resRPhi = std::stof(layer.second.at("resRPhi"));
312+
const float resZ = std::stof(layer.second.at("resZ"));
313+
const float eff = std::stof(layer.second.at("eff"));
314+
const int type = std::stoi(layer.second.at("type"));
315+
const std::string deadPhiRegions = layer.second.at("deadPhiRegions");
307316

308317
// void AddLayer(TString name, float r, float z, float x0, float xrho, float resRPhi = 0.0f, float resZ = 0.0f, float eff = 0.0f, int type = 0);
309-
LOG(info) << " Adding layer " << layer << " r=" << r << " z=" << z << " x0=" << x0 << " xrho=" << xrho << " resRPhi=" << resRPhi << " resZ=" << resZ << " eff=" << eff << " type=" << type << " deadPhiRegions=" << deadPhiRegions;
318+
LOG(info) << " Adding layer " << layer.first << " r=" << r << " z=" << z << " x0=" << x0 << " xrho=" << xrho << " resRPhi=" << resRPhi << " resZ=" << resZ << " eff=" << eff << " type=" << type << " deadPhiRegions=" << deadPhiRegions;
310319

311-
DetLayer* addedLayer = AddLayer(layer.c_str(), r, z, x0, xrho, resRPhi, resZ, eff, type);
312-
if (strlen(deadPhiRegions) > 0) { // Taking it as ccdb path or local file
313-
// Check if it begins with ccdb:
320+
DetLayer* addedLayer = AddLayer(layer.first.c_str(), r, z, x0, xrho, resRPhi, resZ, eff, type);
321+
if (!deadPhiRegions.empty()) { // Taking it as ccdb path or local file
322+
// Check if it begins with ccdb:
314323
if (std::string(deadPhiRegions).rfind("ccdb:", 0) == 0) {
315324
std::string ccdbPath = std::string(deadPhiRegions).substr(5); // remove "ccdb:" prefix
316325
if (ccdbManager == nullptr) {
@@ -321,7 +330,7 @@ void FastTracker::AddGenericDetector(std::string filename, o2::ccdb::BasicCCDBMa
321330
addedLayer->setDeadPhiRegions(g);
322331
} else {
323332
// Taking it as local file
324-
TFile infile(deadPhiRegions, "READ");
333+
TFile infile(deadPhiRegions.c_str(), "READ");
325334
if (!infile.IsOpen()) {
326335
LOG(fatal) << "Cannot open dead phi regions file " << deadPhiRegions;
327336
return;
@@ -331,7 +340,7 @@ void FastTracker::AddGenericDetector(std::string filename, o2::ccdb::BasicCCDBMa
331340
addedLayer->setDeadPhiRegions(g);
332341
}
333342
} else {
334-
LOG(debug) << " No dead phi regions for layer " << layer;
343+
LOG(debug) << " No dead phi regions for layer " << layer.first;
335344
}
336345
}
337346
}

‎ALICE3/Core/FastTracker.h‎

Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -68,6 +68,14 @@ class FastTracker
6868
void AddSiliconALICE3v2(std::vector<float> pixelResolution);
6969
void AddSiliconALICE3(float scaleX0VD, std::vector<float> pixelResolution);
7070
void AddTPC(float phiResMean, float zResMean);
71+
72+
/**
73+
* @brief Parses a TEnv configuration file and returns the key-value pairs split per entry
74+
* @param filename Path to the TEnv configuration file
75+
* @return A map where each key is a layer name and the value is another map of key-value pairs for that layer
76+
*/
77+
std::map<std::string, std::map<std::string, std::string>> parseTEnvConfiguration(std::string filename);
78+
7179
/**
7280
* @brief Adds a generic detector configuration from the specified file.
7381
*

‎ALICE3/TableProducer/OTF/onTheFlyTrackerPid.cxx‎

Lines changed: 25 additions & 26 deletions
Original file line numberDiff line numberDiff line change
@@ -136,38 +136,37 @@ class ToTLUT
136136

137137
bool load(int pdg, const std::string& filename)
138138
{
139-
if (!filename.empty() && strncmp(filename.c_str(), "ccdb:", 5) == 0) {
140-
std::string basePath = std::string(filename).substr(5);
141-
std::string path = basePath + "/PDG_" + std::to_string(pdg);
142-
const std::string outPath = "/tmp/ToTLUTs/";
143-
144-
std::string localFilename = Form("%s/lut_tot_%d.root", outPath.c_str(), pdg);
139+
if (filename.empty()) {
140+
LOG(warning) << "Provided filename is empty for PDG " << pdg;
141+
return false;
142+
}
143+
if (strncmp(filename.c_str(), "ccdb:", 5) == 0) { // Check if filename starts with "ccdb:"
144+
const std::string basePath = std::string(filename).substr(5);
145+
const std::string outPath = "/tmp/ToTLUTs/" + basePath;
146+
const std::string localFilename = outPath + "/snapshot.root";
145147
std::ifstream checkFile(localFilename);
146-
if (!checkFile.is_open()) {
148+
if (!checkFile.is_open()) { // File is not found, need to download it from CCDB
147149
if (!mCcdbManager) {
148150
LOG(fatal) << "CCDB manager not set. Please set it before loading LUT from CCDB.";
149151
}
150152
std::map<std::string, std::string> metadata;
151-
mCcdbManager->getCCDBAccessor().retrieveBlob(path, outPath, metadata, 1);
152-
153-
std::string foundFile = Form("%s/%s/snapshot.root", outPath.c_str(), path.c_str());
154-
std::ifstream testFile(foundFile);
153+
mCcdbManager->getCCDBAccessor().retrieveBlob(basePath, outPath, metadata, 1);
154+
std::ifstream testFile(localFilename);
155155
if (!testFile.is_open()) {
156-
LOG(error) << "Could not find downloaded CCDB file for PDG " << pdg;
156+
LOG(fatal) << "Could not find downloaded CCDB file for PDG " << pdg;
157157
return false;
158158
}
159159
testFile.close();
160-
161-
return load(pdg, foundFile);
162-
} else {
160+
return load(pdg, localFilename);
161+
} else { // File is found, proceed to load it
163162
checkFile.close();
164163
return load(pdg, localFilename);
165164
}
166165
}
167-
166+
// In case the file is already available locally
168167
TFile* f = TFile::Open(filename.c_str());
169168
if (!f || f->IsZombie()) {
170-
LOG(error) << "Failed to open LUT file: " << filename;
169+
LOG(fatal) << "Failed to open LUT file: " << filename;
171170
return false;
172171
}
173172

@@ -397,15 +396,15 @@ struct OnTheFlyTrackerPid {
397396

398397
HistogramRegistry histos{"histos", {}, OutputObjHandlingPolicy::AnalysisObject};
399398

400-
Configurable<std::string> lutTotEl{"lutTotEl", "ccdb:Users/h/hfribert/ToT_LUTs", "ToT LUT for electrons"};
401-
Configurable<std::string> lutTotMu{"lutTotMu", "ccdb:Users/h/hfribert/ToT_LUTs", "ToT LUT for muons"};
402-
Configurable<std::string> lutTotPi{"lutTotPi", "ccdb:Users/h/hfribert/ToT_LUTs", "ToT LUT for pions"};
403-
Configurable<std::string> lutTotKa{"lutTotKa", "ccdb:Users/h/hfribert/ToT_LUTs", "ToT LUT for kaons"};
404-
Configurable<std::string> lutTotPr{"lutTotPr", "ccdb:Users/h/hfribert/ToT_LUTs", "ToT LUT for protons"};
405-
Configurable<std::string> lutTotDe{"lutTotDe", "ccdb:Users/h/hfribert/ToT_LUTs", "ToT LUT for deuteron"};
406-
Configurable<std::string> lutTotTr{"lutTotTr", "ccdb:Users/h/hfribert/ToT_LUTs", "ToT LUT for triton"};
407-
Configurable<std::string> lutTotHe{"lutTotHe", "ccdb:Users/h/hfribert/ToT_LUTs", "ToT LUT for helium-3"};
408-
Configurable<std::string> lutTotAl{"lutTotAl", "ccdb:Users/h/hfribert/ToT_LUTs", "ToT LUT for alphas"};
399+
Configurable<std::string> lutTotEl{"lutTotEl", "ccdb:Users/h/hfribert/ToT_LUTs/PDG_11/", "ToT LUT for electrons"};
400+
Configurable<std::string> lutTotMu{"lutTotMu", "ccdb:Users/h/hfribert/ToT_LUTs/PDG_13/", "ToT LUT for muons"};
401+
Configurable<std::string> lutTotPi{"lutTotPi", "ccdb:Users/h/hfribert/ToT_LUTs/PDG_211/", "ToT LUT for pions"};
402+
Configurable<std::string> lutTotKa{"lutTotKa", "ccdb:Users/h/hfribert/ToT_LUTs/PDG_321/", "ToT LUT for kaons"};
403+
Configurable<std::string> lutTotPr{"lutTotPr", "ccdb:Users/h/hfribert/ToT_LUTs/PDG_2212/", "ToT LUT for protons"};
404+
Configurable<std::string> lutTotDe{"lutTotDe", "ccdb:Users/h/hfribert/ToT_LUTs/PDG_1000010020/", "ToT LUT for deuteron"};
405+
Configurable<std::string> lutTotTr{"lutTotTr", "ccdb:Users/h/hfribert/ToT_LUTs/PDG_1000010030/", "ToT LUT for triton"};
406+
Configurable<std::string> lutTotHe{"lutTotHe", "ccdb:Users/h/hfribert/ToT_LUTs/PDG_1000020030/", "ToT LUT for helium-3"};
407+
Configurable<std::string> lutTotAl{"lutTotAl", "ccdb:Users/h/hfribert/ToT_LUTs/PDG_1000020040/", "ToT LUT for alphas"};
409408

410409
Configurable<float> dBz{"dBz", 20, "magnetic field (kilogauss) for track propagation"};
411410
Configurable<int> maxBarrelLayers{"maxBarrelLayers", 11, "Maximum number of barrel layers"};

‎ALICE3/Tasks/alice3-qa-singleparticle.cxx‎

Lines changed: 30 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -15,13 +15,13 @@
1515
/// \brief Task to monitor the single particle QA, at the particle and track level, showing the tracked and the origin of particles
1616
///
1717

18-
// O2 includes
19-
#include "Framework/AnalysisTask.h"
20-
#include "Framework/runDataProcessing.h"
21-
#include "Framework/HistogramRegistry.h"
22-
#include "Framework/O2DatabasePDGPlugin.h"
23-
#include "TDatabasePDG.h"
24-
#include "TMCProcess.h"
18+
#include <Framework/AnalysisTask.h>
19+
#include <Framework/HistogramRegistry.h>
20+
#include <Framework/O2DatabasePDGPlugin.h>
21+
#include <Framework/runDataProcessing.h>
22+
23+
#include <TDatabasePDG.h>
24+
#include <TMCProcess.h>
2525

2626
using namespace o2;
2727
using namespace o2::framework;
@@ -97,11 +97,9 @@ struct Alice3SingleParticle {
9797
const AxisSpec axisProdz{prodBinsZ, prodMinZ, prodMaxZ, "Prod. Vertex Z (cm)"};
9898
const AxisSpec axisProdRadius{prodBins, 0., 2. * prodMax, "Prod. Vertex Radius (cm)"};
9999

100-
if (!doprocessParticleOnly) {
101-
histos.add("event/VtxX", "Vertex X", kTH1D, {axisVx});
102-
histos.add("event/VtxY", "Vertex Y", kTH1D, {axisVy});
103-
histos.add("event/VtxZ", "Vertex Z", kTH1D, {axisVz});
104-
}
100+
histos.add("event/VtxX", "Vertex X", kTH1D, {axisVx});
101+
histos.add("event/VtxY", "Vertex Y", kTH1D, {axisVy});
102+
histos.add("event/VtxZ", "Vertex Z", kTH1D, {axisVz});
105103

106104
histos.add("particle/PDGs", "Particle PDGs", kTH2D, {axisPDGs, axisCharge});
107105
histos.add("particle/PDGsPrimaries", "Particle PDGs of Primaries", kTH2D, {axisPDGs, axisCharge});
@@ -145,6 +143,7 @@ struct Alice3SingleParticle {
145143
histos.add("particle/Py", "Particle Py " + tit, kTH1D, {axisPy});
146144
histos.add("particle/Pz", "Particle Pz " + tit, kTH1D, {axisPz});
147145

146+
histos.add("particle/daughters/Number", "Number of Daughters " + tit, kTH1D, {{20, -0.5, 19.5}});
148147
histos.add("particle/daughters/PDGs", "Daughters PDGs " + tit, kTH2D, {axisPDGs, axisCharge});
149148
histos.add("particle/daughters/PDGsPrimaries", "Daughters PDGs Primaries of " + tit, kTH2D, {axisPDGs, axisCharge});
150149
histos.add("particle/daughters/PDGsSecondaries", "Daughters PDGs Secondaries of " + tit, kTH2D, {axisPDGs, axisCharge});
@@ -158,6 +157,7 @@ struct Alice3SingleParticle {
158157
histos.add("particle/daughters/prodRadiusVsPt", "Daughters Prod. Vertex Radius " + tit, kTH2D, {axisPt, axisProdRadius});
159158
histos.add("particle/daughters/prodRadius3DVsPt", "Daughters Prod. Vertex Radius XYZ " + tit, kTH2D, {axisPt, axisProdRadius});
160159

160+
histos.add("particle/mothers/Number", "Number of Mothers " + tit, kTH1D, {{20, -0.5, 19.5}});
161161
histos.add("particle/mothers/PDGs", "Mothers PDGs " + tit, kTH2D, {axisPDGs, axisCharge});
162162
histos.add("particle/mothers/PDGsPrimaries", "Mothers PDGs Primaries of " + tit, kTH2D, {axisPDGs, axisCharge});
163163
histos.add("particle/mothers/PDGsSecondaries", "Mothers PDGs Secondaries of " + tit, kTH2D, {axisPDGs, axisCharge});
@@ -167,6 +167,8 @@ struct Alice3SingleParticle {
167167
histos.add("particle/mothers/prodRadiusVsPt", "Mothers Prod. Vertex Radius " + tit, kTH2D, {axisPt, axisProdRadius});
168168
histos.add("particle/mothers/prodRadius3DVsPt", "Mothers Prod. Vertex Radius XYZ " + tit, kTH2D, {axisPt, axisProdRadius});
169169

170+
// Go up one generation
171+
histos.add("particle/mothers/mothers/Number", "Number of Mothers mothers " + tit, kTH1D, {{20, -0.5, 19.5}});
170172
histos.add("particle/mothers/mothers/PDGs", "Mothers mothers PDGs " + tit, kTH2D, {axisPDGs, axisCharge});
171173
histos.add("particle/mothers/mothers/PDGsPrimaries", "Mothers mothers PDGs Primaries of " + tit, kTH2D, {axisPDGs, axisCharge});
172174
histos.add("particle/mothers/mothers/PDGsSecondaries", "Mothers mothers PDGs Secondaries of " + tit, kTH2D, {axisPDGs, axisCharge});
@@ -272,6 +274,7 @@ struct Alice3SingleParticle {
272274
histos.fill(HIST("particle/prodVz"), mcParticle.vz());
273275
if (mcParticle.has_daughters()) {
274276
auto daughters = mcParticle.daughters_as<aod::McParticles>();
277+
histos.fill(HIST("particle/daughters/Number"), daughters.size());
275278
for (const auto& daughter : daughters) {
276279
const auto& pdgStringDau = getPdgCodeString(daughter);
277280
const auto& pdgChargeDau = getCharge(daughter);
@@ -293,6 +296,8 @@ struct Alice3SingleParticle {
293296
histos.fill(HIST("particle/daughters/prodRadiusVsPt"), mcParticle.pt(), std::sqrt(daughter.vx() * daughter.vx() + daughter.vy() * daughter.vy()));
294297
histos.fill(HIST("particle/daughters/prodRadius3DVsPt"), mcParticle.pt(), std::sqrt(daughter.vx() * daughter.vx() + daughter.vy() * daughter.vy() + daughter.vz() * daughter.vz()));
295298
}
299+
} else {
300+
histos.fill(HIST("particle/daughters/Number"), 0.f);
296301
}
297302
if (mcParticle.has_mothers()) {
298303
const auto& mothers = mcParticle.mothers_as<aod::McParticles>();
@@ -415,8 +420,14 @@ struct Alice3SingleParticle {
415420
}
416421
PROCESS_SWITCH(Alice3SingleParticle, processStandard, "Process IU tracks", true);
417422

418-
void processParticleOnly(const aod::McParticles& mcParticles)
423+
void processParticleOnly(const o2::aod::McCollisions& colls,
424+
const aod::McParticles& mcParticles)
419425
{
426+
for (const auto& col : colls) {
427+
histos.fill(HIST("event/VtxX"), col.posX());
428+
histos.fill(HIST("event/VtxY"), col.posY());
429+
histos.fill(HIST("event/VtxZ"), col.posZ());
430+
}
420431
for (const auto& mcParticle : mcParticles) {
421432
const auto& pdgString = getPdgCodeString(mcParticle);
422433
const auto& pdgCharge = getCharge(mcParticle);
@@ -467,6 +478,7 @@ struct Alice3SingleParticle {
467478
histos.fill(HIST("particle/prodVz"), mcParticle.vz());
468479
if (mcParticle.has_daughters()) {
469480
auto daughters = mcParticle.daughters_as<aod::McParticles>();
481+
histos.fill(HIST("particle/daughters/Number"), daughters.size());
470482
for (const auto& daughter : daughters) {
471483
const auto& pdgStringDau = getPdgCodeString(daughter);
472484
const auto& pdgChargeDau = getCharge(daughter);
@@ -484,6 +496,8 @@ struct Alice3SingleParticle {
484496
histos.fill(HIST("particle/daughters/prodRadiusVsPt"), mcParticle.pt(), std::sqrt(mcParticle.vx() * mcParticle.vx() + mcParticle.vy() * mcParticle.vy()));
485497
histos.fill(HIST("particle/daughters/prodRadius3DVsPt"), mcParticle.pt(), std::sqrt(mcParticle.vx() * mcParticle.vx() + mcParticle.vy() * mcParticle.vy() + mcParticle.vz() * mcParticle.vz()));
486498
}
499+
} else {
500+
histos.fill(HIST("particle/daughters/Number"), 0.f);
487501
}
488502
if (mcParticle.has_mothers()) {
489503
auto mothers = mcParticle.mothers_as<aod::McParticles>();
@@ -586,6 +600,7 @@ struct Alice3SingleParticle {
586600
histos.fill(HIST("particle/prodVz"), mcParticle.vz());
587601
if (mcParticle.has_daughters()) {
588602
auto daughters = mcParticle.daughters_as<aod::McParticles>();
603+
histos.fill(HIST("particle/daughters/Number"), daughters.size());
589604
for (const auto& daughter : daughters) {
590605
const auto& pdgStringDau = getPdgCodeString(daughter);
591606
const auto& pdgChargeDau = getCharge(daughter);
@@ -603,6 +618,8 @@ struct Alice3SingleParticle {
603618
histos.fill(HIST("particle/daughters/prodRadiusVsPt"), mcParticle.pt(), std::sqrt(mcParticle.vx() * mcParticle.vx() + mcParticle.vy() * mcParticle.vy()));
604619
histos.fill(HIST("particle/daughters/prodRadius3DVsPt"), mcParticle.pt(), std::sqrt(mcParticle.vx() * mcParticle.vx() + mcParticle.vy() * mcParticle.vy() + mcParticle.vz() * mcParticle.vz()));
605620
}
621+
} else {
622+
histos.fill(HIST("particle/daughters/Number"), 0.f);
606623
}
607624
if (mcParticle.has_mothers()) {
608625
auto mothers = mcParticle.mothers_as<aod::McParticles>();

0 commit comments

Comments
 (0)