Skip to content

Commit 460bdf8

Browse files
axti98axti98
authored
[ALICE3,PWGDQ] Implement Multiplicity and Centrality in ALICE 3 DQ Framework (#17606)
1 parent daf20ea commit 460bdf8

6 files changed

Lines changed: 106 additions & 31 deletions

File tree

‎ALICE3/DataModel/ReducedTablesAlice3.h‎

Lines changed: 9 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -23,6 +23,8 @@
2323
#include "ALICE3/DataModel/OTFPIDTrk.h"
2424
#include "ALICE3/DataModel/OTFRICH.h"
2525
#include "ALICE3/DataModel/OTFTOF.h"
26+
#include "ALICE3/DataModel/collisionAlice3.h"
27+
#include "Common/DataModel/Centrality.h"
2628

2729
#include <Framework/ASoA.h>
2830
#include <Framework/AnalysisDataModel.h>
@@ -41,20 +43,16 @@ DECLARE_SOA_COLUMN(MCPosY, mcPosY, float); //! MC event position Y
4143
DECLARE_SOA_COLUMN(MCPosZ, mcPosZ, float); //! MC event position Z
4244
} // namespace reducedeventalice3
4345

44-
namespace reducedeventmcalice3
45-
{
46-
DECLARE_SOA_COLUMN(MultMCNParticlesEta05, multMCNParticlesEta05, float);
47-
DECLARE_SOA_COLUMN(MultMCNParticlesEta08, multMCNParticlesEta08, float);
48-
DECLARE_SOA_COLUMN(MultMCNParticlesEta10, multMCNParticlesEta10, float);
49-
DECLARE_SOA_COLUMN(MultMCNParticlesEta20, multMCNParticlesEta20, float);
50-
DECLARE_SOA_COLUMN(MultMCNParticlesEta40, multMCNParticlesEta40, float);
51-
} // namespace reducedeventmcalice3
52-
5346
DECLARE_SOA_TABLE_STAGED(ReA3Events, "REA3EVENT", //! Main event information table
5447
o2::soa::Index<>,
5548
collision::PosX, collision::PosY, collision::PosZ, collision::NumContrib,
5649
collision::CollisionTime, collision::CollisionTimeRes, reducedeventalice3::MultDensity);
5750

51+
DECLARE_SOA_TABLE(ReA3EventsExtended, "AOD", "REA3EVENTEXT", //! Event vertex covariance matrix
52+
mult::MultNTracksPV, mult::MultNTracksPVeta1, mult::MultNTracksPVetaHalf,
53+
mult::MultNTracksGlobal, mult::MultNGlobalTracksPV,
54+
cent::CentRun2V0M);
55+
5856
DECLARE_SOA_TABLE(ReducedA3EventsVtxCov, "AOD", "REA3VTXCOV", //! Event vertex covariance matrix
5957
collision::CovXX, collision::CovXY, collision::CovXZ,
6058
collision::CovYY, collision::CovYZ, collision::CovZZ, collision::Chi2);
@@ -65,7 +63,8 @@ DECLARE_SOA_TABLE(ReducedA3EventsInfo, "AOD", "REA3EVENTINFO", //! Main event
6563
DECLARE_SOA_TABLE(ReA3MCEvents, "AOD", "REA3MCEVENT", //! Event level MC truth information
6664
o2::soa::Index<>,
6765
mccollision::GeneratorsID, reducedeventalice3::MCPosX, reducedeventalice3::MCPosY, reducedeventalice3::MCPosZ,
68-
mccollision::T, mccollision::Weight, mccollision::ImpactParameter);
66+
mccollision::T, mccollision::Weight, mccollision::ImpactParameter,
67+
mcmult_alice3::MultMC, mcmult_alice3::MultMC25, mcmult_alice3::MultMC125, mcmult_alice3::MultMC09);
6968

7069
using ReducedA3MCEvent = ReA3MCEvents::iterator;
7170
using ReA3Event = ReA3Events::iterator;

‎ALICE3/TableProducer/alice3DqTableMaker.cxx‎

Lines changed: 9 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -63,8 +63,8 @@ using MyBarrelTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA,
6363
aod::UpgradeTofs, aod::UpgradeRichs, aod::UpgradeRichSignals,
6464
aod::McTrackLabels>;
6565

66-
using MyEvents = soa::Join<aod::Collisions, aod::CollisionsAlice3, aod::McCollisionLabels>;
67-
using MyEventsMC = aod::McCollisions;
66+
using MyEvents = soa::Join<aod::Collisions, aod::CollisionsAlice3, aod::McCollisionLabels, aod::PVMults, aod::MultsGlobal, aod::CentRun2V0Ms>;
67+
using MyEventsMC = soa::Join<aod::McCollisions, aod::MultsMCAlice3>;
6868

6969
constexpr static uint32_t GkEventFillMap = VarManager::ObjTypes::Collision;
7070
constexpr static uint32_t GkEventMcFillMap = VarManager::ObjTypes::CollisionMC;
@@ -84,6 +84,7 @@ struct Alice3DqTableMaker {
8484
Produces<ReA3MCTracks> trackMC;
8585

8686
Produces<ReA3Events> event;
87+
Produces<ReA3EventsExtended> eventExtended;
8788
Produces<ReducedA3EventsVtxCov> eventVtxCov;
8889
Produces<ReducedA3MCEventLabels> eventMClabels;
8990

@@ -345,7 +346,8 @@ struct Alice3DqTableMaker {
345346
fHistMan->FillHistClass("Event_MCTruth", dqefficiency_helpers::varValues());
346347

347348
eventMC(mcCollision.generatorsID(), mcCollision.posX(), mcCollision.posY(), mcCollision.posZ(),
348-
mcCollision.t(), mcCollision.weight(), mcCollision.impactParameter()); // TODO: Determine and fill multiplicity values
349+
mcCollision.t(), mcCollision.weight(), mcCollision.impactParameter(),
350+
mcCollision.multMC(), mcCollision.multMC25(), mcCollision.multMC125(), mcCollision.multMC09());
349351
}
350352
}
351353

@@ -448,6 +450,10 @@ struct Alice3DqTableMaker {
448450
event(collision.posX(), collision.posY(), collision.posZ(), collision.numContrib(),
449451
collision.collisionTime(), collision.collisionTimeRes(), collision.multDensity());
450452

453+
eventExtended(collision.multNTracksPV(), collision.multNTracksPVeta1(), collision.multNTracksPVetaHalf(),
454+
collision.multNTracksGlobal(), collision.multNGlobalTracksPV(),
455+
collision.centRun2V0M());
456+
451457
eventVtxCov(collision.covXX(), collision.covXY(), collision.covXZ(), collision.covYY(), collision.covYZ(), collision.covZZ(), collision.chi2());
452458
eventMClabels(collision.mcCollisionId(), collision.mcMask());
453459

‎ALICE3/Tasks/alice3DqEfficiency.cxx‎

Lines changed: 9 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -111,10 +111,9 @@ DECLARE_SOA_TABLE(OniaMCTruth, "AOD", "MCTRUTHONIA", dqanalysisflags::OniaPt, dq
111111

112112
// TODO: USE PROPER TABLES
113113

114-
using MyEvents = soa::Join<aod::ReA3Events, aod::ReducedA3MCEventLabels>;
115-
using MyEventsSelected = soa::Join<aod::ReA3Events, aod::EventCuts, aod::ReducedA3MCEventLabels>;
116-
using MyEventsVtxCov = soa::Join<aod::ReA3Events, aod::ReducedA3EventsVtxCov, aod::ReducedA3MCEventLabels>;
117-
using MyEventsVtxCovSelected = soa::Join<aod::ReA3Events, aod::ReducedA3EventsVtxCov, aod::EventCuts, aod::ReducedA3MCEventLabels>;
114+
using MyEvents = soa::Join<aod::ReA3Events, aod::ReA3EventsExtended, aod::ReducedA3MCEventLabels>;
115+
using MyEventsVtxCov = soa::Join<aod::ReA3Events, aod::ReA3EventsExtended, aod::ReducedA3EventsVtxCov, aod::ReducedA3MCEventLabels>;
116+
using MyEventsVtxCovSelected = soa::Join<aod::ReA3Events, aod::ReA3EventsExtended, aod::ReducedA3EventsVtxCov, aod::EventCuts, aod::ReducedA3MCEventLabels>;
118117

119118
using MyBarrelAssocs = soa::Join<aod::ReducedA3TracksAssoc, aod::BarrelTrackCuts>;
120119
using MyBarrelAssocsPrefilter = soa::Join<aod::ReducedA3TracksAssoc, aod::BarrelTrackCuts, aod::Prefilter>;
@@ -598,10 +597,10 @@ struct Alice3DqEfficiencyAnalysisPrefilterSelection {
598597
// get the list of cuts that were computed in the barrel track-selection task and create a bit mask
599598
// to mark just the ones we want to apply a prefilter on
600599
string trackCuts;
601-
getTaskOptionValue<string>(context, "analysis-track-selection", "cfgTrackCuts", trackCuts, false);
600+
getTaskOptionValue<string>(context, "alice3-dq-efficiency-analysis-track-selection", "cfgTrackCuts", trackCuts, false);
602601
TString allTrackCutsStr = trackCuts;
603602
// check also the cuts added via JSON and add them to the string of cuts
604-
getTaskOptionValue<string>(context, "analysis-track-selection", "cfgBarrelTrackCutsJSON", trackCuts, false);
603+
getTaskOptionValue<string>(context, "alice3-dq-efficiency-analysis-track-selection", "cfgBarrelTrackCutsJSON", trackCuts, false);
605604
TString addTrackCutsStr = trackCuts;
606605
if (addTrackCutsStr != "") {
607606
std::vector<AnalysisCut*> addTrackCuts = dqcuts::GetCutsFromJSON(addTrackCutsStr.Data());
@@ -852,10 +851,10 @@ struct Alice3DqEfficiencyAnalysisSameEventPairing {
852851

853852
// get the barrel track selection cuts
854853
string tempCuts;
855-
getTaskOptionValue<string>(context, "analysis-track-selection", "cfgTrackCuts", tempCuts, false);
854+
getTaskOptionValue<string>(context, "alice3-dq-efficiency-analysis-track-selection", "cfgTrackCuts", tempCuts, false);
856855
TString tempCutsStr = tempCuts;
857856
// check also the cuts added via JSON and add them to the string of cuts
858-
getTaskOptionValue<string>(context, "analysis-track-selection", "cfgBarrelTrackCutsJSON", tempCuts, false);
857+
getTaskOptionValue<string>(context, "alice3-dq-efficiency-analysis-track-selection", "cfgBarrelTrackCutsJSON", tempCuts, false);
859858
TString addTrackCutsStr = tempCuts;
860859
if (addTrackCutsStr != "") {
861860
std::vector<AnalysisCut*> addTrackCuts = dqcuts::GetCutsFromJSON(addTrackCutsStr.Data());
@@ -1649,10 +1648,10 @@ struct Alice3DqEfficiencyAnalysisAsymmetricPairing {
16491648

16501649
// Get the barrel track selection cuts
16511650
string tempCuts;
1652-
getTaskOptionValue<string>(context, "analysis-track-selection", "cfgTrackCuts", tempCuts, false);
1651+
getTaskOptionValue<string>(context, "alice3-dq-efficiency-analysis-track-selection", "cfgTrackCuts", tempCuts, false);
16531652
TString tempCutsStr = tempCuts;
16541653
// check also the cuts added via JSON and add them to the string of cuts
1655-
getTaskOptionValue<string>(context, "analysis-track-selection", "cfgBarrelTrackCutsJSON", tempCuts, false);
1654+
getTaskOptionValue<string>(context, "alice3-dq-efficiency-analysis-track-selection", "cfgBarrelTrackCutsJSON", tempCuts, false);
16561655
TString addTrackCutsStr = tempCuts;
16571656
if (addTrackCutsStr != "") {
16581657
std::vector<AnalysisCut*> addTrackCuts = dqcuts::GetCutsFromJSON(addTrackCutsStr.Data());

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 19 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -110,6 +110,14 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
110110
hm->AddHistogram(histClass, "CentFT0C_MultTPC", "CentFT0C vs MultTPC", false, 100, 0., 100., VarManager::kCentFT0C, 100, 0., 50000., VarManager::kMultTPC);
111111
hm->AddHistogram(histClass, "CentFT0C_Run", "Cent FT0C", true, 1, -0.5, 0.5, VarManager::kRunNo, 100, 0., 100., VarManager::kCentFT0C, 1, 0, 1, VarManager::kNothing, "", "", "", VarManager::kNothing, VarManager::kNothing, false, true);
112112
}
113+
if (subGroupStr.Contains("alice3multcent")) {
114+
hm->AddHistogram(histClass, "Cent", "Cent", false, 100, 0.0, 100.0, VarManager::kCent);
115+
hm->AddHistogram(histClass, "MultPV", "MultPV", false, 100, 0.0, 10000.0, VarManager::kMultPV);
116+
hm->AddHistogram(histClass, "MultPVeta1", "MultPVeta1", false, 100, 0.0, 10000.0, VarManager::kMultPVeta1);
117+
hm->AddHistogram(histClass, "MultPVetaHalf", "MultPVetaHalf", false, 100, 0.0, 10000.0, VarManager::kMultPVetaHalf);
118+
hm->AddHistogram(histClass, "MultGlobal", "MultGlobal", false, 100, 0.0, 10000.0, VarManager::kMultGlobalTracks);
119+
hm->AddHistogram(histClass, "MultGlobalPV", "MultGlobalPV", false, 100, 0.0, 10000.0, VarManager::kMultGlobalTracksPV);
120+
}
113121
if (subGroupStr.Contains("randomplane")) {
114122
hm->AddHistogram(histClass, "random plane angle", "random Psi", false, 100, -TMath::Pi() / 2, TMath::Pi() / 2, VarManager::kRandomPsi2);
115123
}
@@ -1395,6 +1403,17 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
13951403
}
13961404
hm->AddHistogram(histClass, "Mass_PtLong", "", false, 500, alice3MassBins.data(), VarManager::kMass, 20, pTBins.data(), VarManager::kPt);
13971405
}
1406+
if (subGroupStr.Contains("alice3pbpb")) {
1407+
hm->AddHistogram(histClass, "Mass_Cent", "", false, 125, 0.0, 5.0, VarManager::kMass, 20, 0.0, 100.0, VarManager::kCent);
1408+
hm->AddHistogram(histClass, "Pt_CentFT0C", "", false, 100, 0.0, 10.0, VarManager::kPt, 20, 0.0, 100.0, VarManager::kCent);
1409+
hm->AddHistogram(histClass, "Mass_Pt_CentFT0C", "", false, 75, 1.5, 4.5, VarManager::kMass, 20, 0.0, 10.0, VarManager::kPt, 10, 0.0, 100.0, VarManager::kCent);
1410+
}
1411+
if (subGroupStr.Contains("alice3mult")) {
1412+
hm->AddHistogram(histClass, "Mass_Pt_MultGlobal", "Mass vs pT vs multGlobal", false, 200, 0.0, 5.0, VarManager::kMass, 40, 0.0, 40.0, VarManager::kPt, 100, 0.0, 10000.0, VarManager::kMultGlobalTracks);
1413+
hm->AddHistogram(histClass, "Mass_MultGlobal", "Mass vs Global Tracks", false, 200, 0.0, 5.0, VarManager::kMass, 100, 0, 10000.0, VarManager::kMultGlobalTracks);
1414+
hm->AddHistogram(histClass, "Mass_MultGlobal_VtxZ", "Mass vs Global Tracks vs VtxZ", false, 200, 2.0, 5.0, VarManager::kMass, 100, 0, 10000.0, VarManager::kMultGlobalTracks, 20, -10.0, 10.0, VarManager::kVtxZ);
1415+
hm->AddHistogram(histClass, "VtxZ_MultGlobal", "VtxZ vs VtxNcontribReal", false, 240, -12.0, 12.0, VarManager::kVtxZ, 100, 0, 10000.0, VarManager::kMultGlobalTracks);
1416+
}
13981417
if (subGroupStr.Contains("dielectron-polarization-he-pbpb")) {
13991418
std::array<int, 5> varsHEpbpb = {VarManager::kMass, VarManager::kPt, VarManager::kCentFT0C, VarManager::kCosThetaHE, VarManager::kPhiHE};
14001419
std::array<int, 5> binspT = {100, 30, 10, 10, 10};

‎PWGDQ/Core/VarManager.cxx‎

Lines changed: 30 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -1900,10 +1900,26 @@ void VarManager::SetDefaultVarNames()
19001900
fgVariableUnits[kBGFDDCpf] = "";
19011901
fgVariableNames[kMultDensity] = "dNdeta ALICE3";
19021902
fgVariableUnits[kMultDensity] = "";
1903-
fgVariableNames[kMultMCNParticlesEta40] = "Multiplicity_eta40";
1904-
fgVariableUnits[kMultMCNParticlesEta40] = "";
1905-
fgVariableNames[kMultMCNParticlesEta20] = "Multiplicity_eta20";
1906-
fgVariableUnits[kMultMCNParticlesEta20] = "";
1903+
fgVariableNames[kCent] = "Centrality ALICE3";
1904+
fgVariableUnits[kCent] = "";
1905+
fgVariableNames[kMultPV] = "Mult PV ALICE3";
1906+
fgVariableUnits[kMultPV] = "";
1907+
fgVariableNames[kMultPVeta1] = "Mult PV ALICE3 in |eta| < 1";
1908+
fgVariableUnits[kMultPVeta1] = "";
1909+
fgVariableNames[kMultPVetaHalf] = "Mult PV ALICE3 in |eta| < 0.5";
1910+
fgVariableUnits[kMultPVetaHalf] = "";
1911+
fgVariableNames[kMultGlobalTracks] = "Mult global tracks ALICE3";
1912+
fgVariableUnits[kMultGlobalTracks] = "";
1913+
fgVariableNames[kMultGlobalTracksPV] = "Mult global tracks that are PV ALICE3";
1914+
fgVariableUnits[kMultGlobalTracksPV] = "";
1915+
fgVariableNames[kMultMCNParticlesAll] = "Multiplicity MC all";
1916+
fgVariableUnits[kMultMCNParticlesAll] = "";
1917+
fgVariableNames[kMultMCNParticlesEta25] = "Multiplicity MC |eta| < 2.5";
1918+
fgVariableUnits[kMultMCNParticlesEta25] = "";
1919+
fgVariableNames[kMultMCNParticlesEta125] = "Multiplicity MC |eta| < 1.25";
1920+
fgVariableUnits[kMultMCNParticlesEta125] = "";
1921+
fgVariableNames[kMultMCNParticlesEta09] = "Multiplicity MC |eta| < 0.9";
1922+
fgVariableUnits[kMultMCNParticlesEta09] = "";
19071923
fgVariableNames[kIsReconstructed] = "is track reconstructed";
19081924
fgVariableUnits[kIsReconstructed] = "";
19091925
fgVariableNames[kNSiliconHits] = "Number of hits in silicon layers";
@@ -2834,8 +2850,16 @@ void VarManager::SetDefaultVarNames()
28342850
fgVarNamesMap["kBBFDDCpf"] = kBBFDDCpf;
28352851
fgVarNamesMap["kBGFDDCpf"] = kBGFDDCpf;
28362852
fgVarNamesMap["kMultDensity"] = kMultDensity;
2837-
fgVarNamesMap["kMultMCNParticlesEta40"] = kMultMCNParticlesEta40;
2838-
fgVarNamesMap["kMultMCNParticlesEta20"] = kMultMCNParticlesEta20;
2853+
fgVarNamesMap["kCent"] = kCent;
2854+
fgVarNamesMap["kMultPV"] = kMultPV;
2855+
fgVarNamesMap["kMultPVeta1"] = kMultPVeta1;
2856+
fgVarNamesMap["kMultPVetaHalf"] = kMultPVetaHalf;
2857+
fgVarNamesMap["kMultGlobalTracks"] = kMultGlobalTracks;
2858+
fgVarNamesMap["kMultGlobalTracksPV"] = kMultGlobalTracksPV;
2859+
fgVarNamesMap["kMultMCNParticlesAll"] = kMultMCNParticlesAll;
2860+
fgVarNamesMap["kMultMCNParticlesEta25"] = kMultMCNParticlesEta25;
2861+
fgVarNamesMap["kMultMCNParticlesEta125"] = kMultMCNParticlesEta125;
2862+
fgVarNamesMap["kMultMCNParticlesEta09"] = kMultMCNParticlesEta09;
28392863
fgVarNamesMap["kIsReconstructed"] = kIsReconstructed;
28402864
fgVarNamesMap["kNSiliconHits"] = kNSiliconHits;
28412865
fgVarNamesMap["kNTPCHits"] = kNTPCHits;

‎PWGDQ/Core/VarManager.h‎

Lines changed: 30 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1117,8 +1117,16 @@ class VarManager : public TObject
11171117

11181118
// ALICE 3 Variables
11191119
kMultDensity,
1120-
kMultMCNParticlesEta40,
1121-
kMultMCNParticlesEta20,
1120+
kCent,
1121+
kMultPV,
1122+
kMultPVeta1,
1123+
kMultPVetaHalf,
1124+
kMultGlobalTracks,
1125+
kMultGlobalTracksPV,
1126+
kMultMCNParticlesAll,
1127+
kMultMCNParticlesEta25,
1128+
kMultMCNParticlesEta125,
1129+
kMultMCNParticlesEta09,
11221130
kIsReconstructed,
11231131
kNSiliconHits,
11241132
kNTPCHits,
@@ -7197,6 +7205,12 @@ void VarManager::FillEventAlice3(T const& event, float* values)
71977205
values[kVtxChi2] = event.chi2();
71987206
values[kCollisionTime] = event.collisionTime();
71997207
values[kCollisionTimeRes] = event.collisionTimeRes();
7208+
values[kCent] = event.centRun2V0M();
7209+
values[kMultPV] = event.multNTracksPV();
7210+
values[kMultPVeta1] = event.multNTracksPVeta1();
7211+
values[kMultPVetaHalf] = event.multNTracksPVetaHalf();
7212+
values[kMultGlobalTracks] = event.multNTracksGlobal();
7213+
values[kMultGlobalTracksPV] = event.multNGlobalTracksPV();
72007214
}
72017215

72027216
if constexpr ((fillMap & ReducedEvent) > 0) {
@@ -7207,6 +7221,12 @@ void VarManager::FillEventAlice3(T const& event, float* values)
72077221
values[kVtxNcontrib] = event.numContrib();
72087222
values[kCollisionTime] = event.collisionTime();
72097223
values[kCollisionTimeRes] = event.collisionTimeRes();
7224+
values[kCent] = event.centRun2V0M();
7225+
values[kMultPV] = event.multNTracksPV();
7226+
values[kMultPVeta1] = event.multNTracksPVeta1();
7227+
values[kMultPVetaHalf] = event.multNTracksPVetaHalf();
7228+
values[kMultGlobalTracks] = event.multNTracksGlobal();
7229+
values[kMultGlobalTracksPV] = event.multNGlobalTracksPV();
72107230
}
72117231
if constexpr ((fillMap & ReducedEventVtxCov) > 0) {
72127232
values[kVtxCovXX] = event.covXX();
@@ -7227,6 +7247,10 @@ void VarManager::FillEventAlice3(T const& event, float* values)
72277247
values[kMCEventTime] = event.t();
72287248
values[kMCEventWeight] = event.weight();
72297249
values[kMCEventImpParam] = event.impactParameter();
7250+
values[kMultMCNParticlesAll] = event.multMC();
7251+
values[kMultMCNParticlesEta25] = event.multMC25();
7252+
values[kMultMCNParticlesEta125] = event.multMC125();
7253+
values[kMultMCNParticlesEta09] = event.multMC09();
72307254
}
72317255

72327256
if constexpr ((fillMap & ReducedEventMC) > 0) {
@@ -7238,6 +7262,10 @@ void VarManager::FillEventAlice3(T const& event, float* values)
72387262
values[kMCEventTime] = event.t();
72397263
values[kMCEventWeight] = event.weight();
72407264
values[kMCEventImpParam] = event.impactParameter();
7265+
values[kMultMCNParticlesAll] = event.multMC();
7266+
values[kMultMCNParticlesEta25] = event.multMC25();
7267+
values[kMultMCNParticlesEta125] = event.multMC125();
7268+
values[kMultMCNParticlesEta09] = event.multMC09();
72417269
}
72427270
}
72437271

0 commit comments

Comments
 (0)