Skip to content

Commit a55c6dd

Browse files
cnkostercnkoster
authored
Merge pull request #97 from alibuild/alibot-cleanup-16643
[PWGCF] Please consider the following formatting changes to #16643
2 parents 63ab4d1 + 293fded commit a55c6dd

2 files changed

Lines changed: 80 additions & 85 deletions

File tree

‎PWGCF/Flow/TableProducer/zdcQVectors.cxx‎

Lines changed: 19 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -350,7 +350,6 @@ struct ZdcQVectors {
350350
registry.add<TH2>("QA/ZNC_pm2_vs_Centrality", "ZNC_pm2_vs_Centrality", kTH2D, {{100, 0, 100}, {100, 0, 1}});
351351
registry.add<TH2>("QA/ZNC_pm3_vs_Centrality", "ZNC_pm3_vs_Centrality", kTH2D, {{100, 0, 100}, {100, 0, 1}});
352352
registry.add<TH2>("QA/ZNC_pm4_vs_Centrality", "ZNC_pm4_vs_Centrality", kTH2D, {{100, 0, 100}, {100, 0, 1}});
353-
354353
}
355354

356355
// Tower mean energies vs. centrality used for tower gain equalisation
@@ -378,8 +377,8 @@ struct ZdcQVectors {
378377
registry.add<TProfile2D>("CutAnalysis/hvertex_vy", "hvertex_vy", kTProfile2D, {{1, 0., 1.}, {nEventSelections + 5, 0, nEventSelections + 5}});
379378
registry.add<TProfile2D>("CutAnalysis/hvertex_vz", "hvertex_vz", kTProfile2D, {{1, 0., 1.}, {nEventSelections + 5, 0, nEventSelections + 5}});
380379
}
381-
registry.addClone("recentering/before/", "recentering/after/");
382-
registry.addClone("QA/before/", "QA/after/");
380+
registry.addClone("recentering/before/", "recentering/after/");
381+
registry.addClone("QA/before/", "QA/after/");
383382
}
384383
}
385384

@@ -755,7 +754,7 @@ struct ZdcQVectors {
755754
const auto& foundBC = collision.foundBC_as<BCsRun3>();
756755
runnumber = foundBC.runNumber();
757756

758-
if (cfgFillHistRegistry && !cfgFillNothing){
757+
if (cfgFillHistRegistry && !cfgFillNothing) {
759758
registry.fill(HIST("QA/centrality_before"), cent);
760759
}
761760

@@ -841,7 +840,8 @@ struct ZdcQVectors {
841840
}
842841

843842
// load the calibration histos for iteration 0 step 0 (Energy Calibration)
844-
if(!cfgNoGain) loadCalibrations<kEnergyCal>(cfgEnergyCal.value);
843+
if (!cfgNoGain)
844+
loadCalibrations<kEnergyCal>(cfgEnergyCal.value);
845845
// load the calibrations for the mean v
846846
loadCalibrations<kMeanv>(cfgMeanv.value);
847847

@@ -869,25 +869,25 @@ struct ZdcQVectors {
869869

870870
// Now start gain equalisation!
871871
// Fill the list with calibration constants.
872-
if(!cfgNoGain){
872+
if (!cfgNoGain) {
873873
for (int tower = 0; tower < (nTowers + 2); tower++) {
874-
meanEZN[tower] = getCorrection<TProfile2D, kEnergyCal>(namesEcal[tower].Data());
875-
}
876-
}
874+
meanEZN[tower] = getCorrection<TProfile2D, kEnergyCal>(namesEcal[tower].Data());
875+
}
876+
}
877877

878878
// Use the calibration constants but now only loop over towers 1-4
879879
int calibtower = 0;
880880
std::vector<int> towersNocom = {1, 2, 3, 4, 6, 7, 8, 9};
881881

882882
for (const auto& tower : towersNocom) {
883-
if(cfgNoGain) {
884-
e[calibtower] = eZN[calibtower];
885-
} else {
883+
if (cfgNoGain) {
884+
e[calibtower] = eZN[calibtower];
885+
} else {
886886
if (meanEZN[tower] > 0) {
887887
double ecommon = (tower > nTowersPerSide) ? meanEZN[5] : meanEZN[0];
888888
e[calibtower] = eZN[calibtower] * (0.25 * ecommon) / meanEZN[tower];
889889
}
890-
}
890+
}
891891
calibtower++;
892892
}
893893

@@ -918,7 +918,6 @@ struct ZdcQVectors {
918918

919919
registry.fill(HIST("QA/ZNC_pmC_vs_Centrality"), centrality, zdcCol.energyCommonZNC());
920920
registry.fill(HIST("QA/ZNC_pmSUM_vs_Centrality"), centrality, sumZNC);
921-
922921

923922
registry.fill(HIST("QA/ZNA_pm1_vs_Centrality"), centrality, e[0] / sumZNA);
924923
registry.fill(HIST("QA/ZNA_pm2_vs_Centrality"), centrality, e[1] / sumZNA);
@@ -1002,13 +1001,12 @@ struct ZdcQVectors {
10021001
corrQyA.push_back(getCorrection<THnSparse, kRec>(names[0][1].Data(), it, 1));
10031002
corrQxC.push_back(getCorrection<THnSparse, kRec>(names[0][2].Data(), it, 1));
10041003
corrQyC.push_back(getCorrection<THnSparse, kRec>(names[0][3].Data(), it, 1));
1005-
10061004

10071005
if (cfgFillHistRegistry && !cfgFillNothing) {
1008-
registry.get<TH2>(HIST("recentering/QXA_vs_iteration"))->Fill(pb + 1, q[0] - std::accumulate(corrQxA.begin(), corrQxA.end(), 0.0));
1009-
registry.get<TH2>(HIST("recentering/QYA_vs_iteration"))->Fill(pb + 1, q[1] - std::accumulate(corrQyA.begin(), corrQyA.end(), 0.0));
1010-
registry.get<TH2>(HIST("recentering/QXC_vs_iteration"))->Fill(pb + 1, q[2] - std::accumulate(corrQxC.begin(), corrQxC.end(), 0.0));
1011-
registry.get<TH2>(HIST("recentering/QYC_vs_iteration"))->Fill(pb + 1, q[3] - std::accumulate(corrQyC.begin(), corrQyC.end(), 0.0));
1006+
registry.get<TH2>(HIST("recentering/QXA_vs_iteration"))->Fill(pb + 1, q[0] - std::accumulate(corrQxA.begin(), corrQxA.end(), 0.0));
1007+
registry.get<TH2>(HIST("recentering/QYA_vs_iteration"))->Fill(pb + 1, q[1] - std::accumulate(corrQyA.begin(), corrQyA.end(), 0.0));
1008+
registry.get<TH2>(HIST("recentering/QXC_vs_iteration"))->Fill(pb + 1, q[2] - std::accumulate(corrQxC.begin(), corrQxC.end(), 0.0));
1009+
registry.get<TH2>(HIST("recentering/QYC_vs_iteration"))->Fill(pb + 1, q[3] - std::accumulate(corrQyC.begin(), corrQyC.end(), 0.0));
10121010
}
10131011
pb++;
10141012

@@ -1023,7 +1021,7 @@ struct ZdcQVectors {
10231021
registry.get<TH2>(HIST("recentering/QYA_vs_iteration"))->Fill(pb + 1, q[1] - std::accumulate(corrQyA.begin(), corrQyA.end(), 0.0));
10241022
registry.get<TH2>(HIST("recentering/QXC_vs_iteration"))->Fill(pb + 1, q[2] - std::accumulate(corrQxC.begin(), corrQxC.end(), 0.0));
10251023
registry.get<TH2>(HIST("recentering/QYC_vs_iteration"))->Fill(pb + 1, q[3] - std::accumulate(corrQyC.begin(), corrQyC.end(), 0.0));
1026-
}
1024+
}
10271025

10281026
pb++;
10291027
}
@@ -1039,7 +1037,7 @@ struct ZdcQVectors {
10391037
registry.get<TH2>(HIST("recentering/QYA_vs_iteration"))->Fill(pb + 1, q[1] - std::accumulate(corrQyA.begin(), corrQyA.end(), 0.0));
10401038
registry.get<TH2>(HIST("recentering/QXC_vs_iteration"))->Fill(pb + 1, q[2] - std::accumulate(corrQxC.begin(), corrQxC.end(), 0.0));
10411039
registry.get<TH2>(HIST("recentering/QYC_vs_iteration"))->Fill(pb + 1, q[3] - std::accumulate(corrQyC.begin(), corrQyC.end(), 0.0));
1042-
}
1040+
}
10431041
pb++;
10441042
}
10451043
}

‎PWGCF/Flow/Tasks/flowSP.cxx‎

Lines changed: 61 additions & 64 deletions
Original file line numberDiff line numberDiff line change
@@ -258,7 +258,7 @@ struct FlowSP {
258258
double meanPxC = 0;
259259
} spm;
260260

261-
struct ptMaps {
261+
struct ptMaps {
262262
TProfile* meanPTMap = new TProfile("meanPTMap", "meanPTMap", 8, -0.8, 0.8);
263263
TProfile* meanPTMapPos = new TProfile("meanPTMapPos", "meanPTMapPos", 8, -0.8, 0.8);
264264
TProfile* meanPTMapNeg = new TProfile("meanPTMapNeg", "meanPTMapNeg", 8, -0.8, 0.8);
@@ -271,7 +271,7 @@ struct FlowSP {
271271

272272
TProfile* relPxCNeg = new TProfile("relPxCNeg", "relPxCNeg", 8, -0.8, 0.8);
273273
TProfile* relPxCPos = new TProfile("relPxCPos", "relPxCPos", 8, -0.8, 0.8);
274-
} ptmaps;
274+
} ptmaps;
275275

276276
OutputObj<GFWWeights> fWeights{GFWWeights("weights")};
277277
OutputObj<GFWWeights> fWeightsPOS{GFWWeights("weights_positive")};
@@ -700,42 +700,42 @@ struct FlowSP {
700700
fMultCutHigh = new TF1("fMultCutHigh", "[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x", 0, 100);
701701

702702
std::vector<double> paramsMultPVCut;
703-
std::vector<double> paramsMultCut;
704-
int y2023 = 2023;
705-
int y2024 = 2024;
706-
std::array<int> nSigma = {1,2,3};
707-
708-
if(cfg.cUsePredeFinedSigma){
709-
if(cfg.cUsePredeFinedSigmaYear == y2023) {
710-
if(cfg.cUsePredeFinedSigmaNsigma == nSigma[0]){
711-
paramsMultPVCut = {2615.47, -90.5747, 1.25125, -0.00847075, 2.41183e-05, 3399.72, -121.652, 1.84077, -0.0142886, 4.71449e-05};
712-
paramsMultCut = {1716.84, -56.5663, 0.715202, -0.00426007, 1.05075e-05, 2550.82, -87.4873, 1.22205, -0.00852644, 2.54248e-05 };
713-
} else if(cfg.cUsePredeFinedSigmaNsigma == nSigma[1]){
714-
paramsMultPVCut = {2223.49, -75.1444, 0.963572, -0.00570399, 1.34877e-05, 3790.99, -137.064, 2.13044, -0.017122, 5.82834e-05};
715-
paramsMultCut = {1301.56, -41.4615, 0.478224, -0.00239449, 4.46966e-06, 2967.6, -102.927, 1.47488, -0.0106534, 3.28622e-05};
716-
} else if(cfg.cUsePredeFinedSigmaNsigma == nSigma[2]){
717-
paramsMultPVCut = {1837.75, -60.852, 0.724331, -0.00366975, 6.47562e-06, 4182.12, -152.459, 2.41955, -0.0199481, 6.93894e-05};
718-
paramsMultCut = {885.976, -26.3397, 0.240114, -0.000496168, -1.82704e-06, 3384.43, -118.377, 1.72823, -0.0127887, 4.03432e-05};
719-
} else {
720-
LOGF(fatal, "nSigma can only be 1-3 please reset the variable or give the parameters manually and set cfg.cUsePredeFinedSigma to FALSE");
721-
}
722-
} else if(cfg.cUsePredeFinedSigmaYear == y2024){
723-
if(cfg.cUsePredeFinedSigmaNsigma == nSigma[0]){
724-
paramsMultPVCut = {2726.93, -100.128, 1.45046, -0.0099354, 2.71182e-05, 3404.72, -126.569, 1.92500, -0.0142653, 4.31645e-05};
725-
paramsMultCut = {1858.77, -66.6070, 0.929146, -0.00606961, 1.57639e-05 , 2672.43, -96.7708, 1.39109 , -0.00942498, 2.54268e-05};
726-
} else if(cfg.cUsePredeFinedSigmaNsigma == nSigma[1]){
727-
paramsMultPVCut = {2390.04, -87.3154, 1.23176, -0.00806869, 2.06624e-05, 3744.26, -139.927, 2.16863, -0.0165329, 5.17269e-05};
728-
paramsMultCut = {1451.23, -51.4314, 0.694609, -0.00433959, 1.06698e-05, 3080.42, -112.071, 1.63166 , -0.0112533 , 3.10348e-05};
729-
} else if(cfg.cUsePredeFinedSigmaNsigma == nSigma[2]){
730-
paramsMultPVCut = {2053.64, -74.5950, 1.01563, -0.00621473, 1.41276e-05, 4083.79, -153.304, 2.41333, -0.0188198, 6.03974e-05};
731-
paramsMultCut = {1042.50, -35.9374, 0.440681, -0.00222218, 3.20643e-06, 3488.53, -127.396, 1.87339 , -0.0131007, 3.67434e-05};
732-
} else {
733-
LOGF(fatal, "cUsePredeFinedSigmaNsigma can only be 1-3 please reset the variable or give the parameters manually and set cUsePredeFinedSigma to FALSE");
734-
}
735-
} else {
736-
LOGF(fatal, "cUsePredeFinedSigmaYear can only be 2023/2024 please reset the variable or give the parameters manually and set cUsePredeFinedSigma to FALSE");
703+
std::vector<double> paramsMultCut;
704+
int y2023 = 2023;
705+
int y2024 = 2024;
706+
std::array<int> nSigma = {1, 2, 3};
707+
708+
if (cfg.cUsePredeFinedSigma) {
709+
if (cfg.cUsePredeFinedSigmaYear == y2023) {
710+
if (cfg.cUsePredeFinedSigmaNsigma == nSigma[0]) {
711+
paramsMultPVCut = {2615.47, -90.5747, 1.25125, -0.00847075, 2.41183e-05, 3399.72, -121.652, 1.84077, -0.0142886, 4.71449e-05};
712+
paramsMultCut = {1716.84, -56.5663, 0.715202, -0.00426007, 1.05075e-05, 2550.82, -87.4873, 1.22205, -0.00852644, 2.54248e-05};
713+
} else if (cfg.cUsePredeFinedSigmaNsigma == nSigma[1]) {
714+
paramsMultPVCut = {2223.49, -75.1444, 0.963572, -0.00570399, 1.34877e-05, 3790.99, -137.064, 2.13044, -0.017122, 5.82834e-05};
715+
paramsMultCut = {1301.56, -41.4615, 0.478224, -0.00239449, 4.46966e-06, 2967.6, -102.927, 1.47488, -0.0106534, 3.28622e-05};
716+
} else if (cfg.cUsePredeFinedSigmaNsigma == nSigma[2]) {
717+
paramsMultPVCut = {1837.75, -60.852, 0.724331, -0.00366975, 6.47562e-06, 4182.12, -152.459, 2.41955, -0.0199481, 6.93894e-05};
718+
paramsMultCut = {885.976, -26.3397, 0.240114, -0.000496168, -1.82704e-06, 3384.43, -118.377, 1.72823, -0.0127887, 4.03432e-05};
719+
} else {
720+
LOGF(fatal, "nSigma can only be 1-3 please reset the variable or give the parameters manually and set cfg.cUsePredeFinedSigma to FALSE");
721+
}
722+
} else if (cfg.cUsePredeFinedSigmaYear == y2024) {
723+
if (cfg.cUsePredeFinedSigmaNsigma == nSigma[0]) {
724+
paramsMultPVCut = {2726.93, -100.128, 1.45046, -0.0099354, 2.71182e-05, 3404.72, -126.569, 1.92500, -0.0142653, 4.31645e-05};
725+
paramsMultCut = {1858.77, -66.6070, 0.929146, -0.00606961, 1.57639e-05, 2672.43, -96.7708, 1.39109, -0.00942498, 2.54268e-05};
726+
} else if (cfg.cUsePredeFinedSigmaNsigma == nSigma[1]) {
727+
paramsMultPVCut = {2390.04, -87.3154, 1.23176, -0.00806869, 2.06624e-05, 3744.26, -139.927, 2.16863, -0.0165329, 5.17269e-05};
728+
paramsMultCut = {1451.23, -51.4314, 0.694609, -0.00433959, 1.06698e-05, 3080.42, -112.071, 1.63166, -0.0112533, 3.10348e-05};
729+
} else if (cfg.cUsePredeFinedSigmaNsigma == nSigma[2]) {
730+
paramsMultPVCut = {2053.64, -74.5950, 1.01563, -0.00621473, 1.41276e-05, 4083.79, -153.304, 2.41333, -0.0188198, 6.03974e-05};
731+
paramsMultCut = {1042.50, -35.9374, 0.440681, -0.00222218, 3.20643e-06, 3488.53, -127.396, 1.87339, -0.0131007, 3.67434e-05};
732+
} else {
733+
LOGF(fatal, "cUsePredeFinedSigmaNsigma can only be 1-3 please reset the variable or give the parameters manually and set cUsePredeFinedSigma to FALSE");
737734
}
738-
} else {
735+
} else {
736+
LOGF(fatal, "cUsePredeFinedSigmaYear can only be 2023/2024 please reset the variable or give the parameters manually and set cUsePredeFinedSigma to FALSE");
737+
}
738+
} else {
739739
paramsMultPVCut = cfg.cEvSelsMultPv;
740740
paramsMultCut = cfg.cEvSelsMult;
741741
}
@@ -1192,20 +1192,20 @@ struct FlowSP {
11921192
if (track.eta() < 0)
11931193
scale = -1.0;
11941194

1195-
const double invSqrtQQ = 1.0 / std::sqrt(std::fabs(spm.corrQQ));
1195+
const double invSqrtQQ = 1.0 / std::sqrt(std::fabs(spm.corrQQ));
11961196
const double invSqrtQQx = 1.0 / std::sqrt(std::fabs(spm.corrQQx));
11971197
const double invSqrtQQy = 1.0 / std::sqrt(std::fabs(spm.corrQQy));
1198-
const double uqA = spm.uy * spm.qyA + spm.ux * spm.qxA;
1199-
const double uqC = spm.uy * spm.qyC + spm.ux * spm.qxC;
1200-
const double invMeanPtQQ = 1.0 / spm.meanPtWeight;
1198+
const double uqA = spm.uy * spm.qyA + spm.ux * spm.qxA;
1199+
const double uqC = spm.uy * spm.qyC + spm.ux * spm.qxC;
1200+
const double invMeanPtQQ = 1.0 / spm.meanPtWeight;
12011201

12021202
if (cfg.cFillGeneralV1Histos) {
1203-
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnAodd"), track.pt(), track.eta(), spm.centrality, scale * (uqA) * invSqrtQQ, weight);
1204-
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnCodd"), track.pt(), track.eta(), spm.centrality, scale * (uqC) * invSqrtQQ, weight);
1203+
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnAodd"), track.pt(), track.eta(), spm.centrality, scale * (uqA)*invSqrtQQ, weight);
1204+
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnCodd"), track.pt(), track.eta(), spm.centrality, scale * (uqC)*invSqrtQQ, weight);
12051205
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnOdd"), track.pt(), track.eta(), spm.centrality, scale * 0.5 * ((uqA) - (uqC)) * invSqrtQQ, weight);
12061206
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnEven"), track.pt(), track.eta(), spm.centrality, 0.5 * ((uqA) + (uqC)) * invSqrtQQ, weight);
1207-
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnA"), track.pt(), track.eta(), spm.centrality, (uqA) * invSqrtQQ, weight);
1208-
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnC"), track.pt(), track.eta(), spm.centrality, (uqC) * invSqrtQQ, weight);
1207+
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnA"), track.pt(), track.eta(), spm.centrality, (uqA)*invSqrtQQ, weight);
1208+
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnC"), track.pt(), track.eta(), spm.centrality, (uqC)*invSqrtQQ, weight);
12091209
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnCSetPlane"), track.pt(), track.eta(), spm.centrality, (spm.uy + spm.ux) * invSqrtQQ, weight);
12101210
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("vnASetPlane"), track.pt(), track.eta(), spm.centrality, (minusQ * spm.ux - spm.uy) * invSqrtQQ, weight);
12111211
}
@@ -1258,7 +1258,7 @@ struct FlowSP {
12581258
if (!cfg.cFillTrackQA)
12591259
return;
12601260

1261-
double weight = spm.wacc[ct][par] * spm.weff[ct][par] * spm.centWeight;
1261+
double weight = spm.wacc[ct][par] * spm.weff[ct][par] * spm.centWeight;
12621262

12631263
static constexpr std::string_view Time[] = {"before/", "after/"};
12641264
// NOTE: species[kUnidentified] = "" (when nocfg.cTrackSelDo) {
@@ -1492,30 +1492,28 @@ struct FlowSP {
14921492
if (centW <= 0) {
14931493
spm.centWeight = 0.;
14941494
LOGF(fatal, "Centrality weight cannot be negative .. setting to 0. for (%.2f)", spm.centrality);
1495-
} else {
1496-
spm.centWeight = centW;
1495+
} else {
1496+
spm.centWeight = centW;
14971497
}
14981498
}
14991499

15001500
// Always fill centrality histogram after event selections!
15011501
histos.fill(HIST("hCentralityCW"), spm.centrality);
15021502

1503-
15041503
fillEventQA<kAfter>(collision, tracks);
15051504

1506-
ptmaps.meanPTMap->Reset();
1507-
ptmaps.meanPTMapPos->Reset();
1508-
ptmaps.meanPTMapNeg->Reset();
1505+
ptmaps.meanPTMap->Reset();
1506+
ptmaps.meanPTMapPos->Reset();
1507+
ptmaps.meanPTMapNeg->Reset();
15091508

1510-
ptmaps.relPxA->Reset();
1511-
ptmaps.relPxC->Reset();
1509+
ptmaps.relPxA->Reset();
1510+
ptmaps.relPxC->Reset();
15121511

1513-
ptmaps.relPxANeg->Reset();
1514-
ptmaps.relPxAPos->Reset();
1515-
1516-
ptmaps.relPxCNeg->Reset();
1517-
ptmaps.relPxCPos->Reset();
1512+
ptmaps.relPxANeg->Reset();
1513+
ptmaps.relPxAPos->Reset();
15181514

1515+
ptmaps.relPxCNeg->Reset();
1516+
ptmaps.relPxCPos->Reset();
15191517

15201518
double sumPxAEvent = 0;
15211519
double meanPxEventCount = 0;
@@ -1641,11 +1639,10 @@ struct FlowSP {
16411639

16421640
} // end of track loop
16431641

1642+
double meanPxAEvent = 0;
1643+
double meanPxCEvent = 0;
16441644

1645-
double meanPxAEvent = 0;
1646-
double meanPxCEvent = 0;
1647-
1648-
if(meanPxEventCount > 0) {
1645+
if (meanPxEventCount > 0) {
16491646
meanPxAEvent = sumPxAEvent / meanPxEventCount;
16501647
meanPxCEvent = sumPxCEvent / meanPxEventCount;
16511648
}
@@ -2061,7 +2058,7 @@ struct FlowSP {
20612058
auto trackSlice = tracks.sliceBy(trackPerCollision, col.globalIndex());
20622059

20632060
// get filtered tracks that belong to reconstructed collision
2064-
// mcReco uses filtered tracks
2061+
// mcReco uses filtered tracks
20652062
auto filteredTrackSlice = filteredTracks.sliceBy(trackPerCollision, col.globalIndex());
20662063

20672064
spm.centrality = col.centFT0C();

0 commit comments

Comments
 (0)