Skip to content

Commit c22d11f

Browse files
cnkosteralibuild
cnkoster
and
alibuild
authored
[PWGCF] flowSP: add EBE mean px/pt and overall mean px/pt (#14957)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent 0d08253 commit c22d11f

File tree

2 files changed

+78
-50
lines changed

2 files changed

+78
-50
lines changed

‎PWGCF/Flow/TableProducer/zdcQVectors.cxx‎

Lines changed: 23 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -99,6 +99,7 @@ std::vector<float> v(3, 0); // vx, vy, vz
9999
bool isSelected = true;
100100
std::vector<float> cents; // centrality estimaters
101101
uint64_t timestamp = 0;
102+
double rsTimestamp = 0;
102103

103104
} // namespace o2::analysis::qvectortask
104105

@@ -529,7 +530,7 @@ struct ZdcQVectors {
529530
}
530531

531532
template <FillType ft>
532-
inline void fillCommonRegistry(double qxa, double qya, double qxc, double qyc, std::vector<float> v, double centrality)
533+
inline void fillCommonRegistry(double qxa, double qya, double qxc, double qyc, std::vector<float> v, double centrality, double rsTimestamp)
533534
{
534535
// loop for filling multiple histograms with different naming patterns
535536
// Always fill the uncentered "raw" Q-vector histos!
@@ -580,15 +581,15 @@ struct ZdcQVectors {
580581
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQYA_QXC_vs_vz"), v[2], qya * qxc);
581582
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQXA_QYC_vs_vz"), v[2], qxa * qyc);
582583

583-
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQXA_vs_timestamp"), rescaleTimestamp(timestamp, runnumber), qxa);
584-
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQYA_vs_timestamp"), rescaleTimestamp(timestamp, runnumber), qya);
585-
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQXC_vs_timestamp"), rescaleTimestamp(timestamp, runnumber), qxc);
586-
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQYC_vs_timestamp"), rescaleTimestamp(timestamp, runnumber), qyc);
584+
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQXA_vs_timestamp"), rsTimestamp, qxa);
585+
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQYA_vs_timestamp"), rsTimestamp, qya);
586+
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQXC_vs_timestamp"), rsTimestamp, qxc);
587+
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQYC_vs_timestamp"), rsTimestamp, qyc);
587588

588-
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQXA_QXC_vs_timestamp"), rescaleTimestamp(timestamp, runnumber), qxa * qxc);
589-
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQYA_QYC_vs_timestamp"), rescaleTimestamp(timestamp, runnumber), qya * qyc);
590-
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQYA_QXC_vs_timestamp"), rescaleTimestamp(timestamp, runnumber), qya * qxc);
591-
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQXA_QYC_vs_timestamp"), rescaleTimestamp(timestamp, runnumber), qxa * qyc);
589+
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQXA_QXC_vs_timestamp"), rsTimestamp, qxa * qxc);
590+
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQYA_QYC_vs_timestamp"), rsTimestamp, qya * qyc);
591+
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQYA_QXC_vs_timestamp"), rsTimestamp, qya * qxc);
592+
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/hQXA_QYC_vs_timestamp"), rsTimestamp, qxa * qyc);
592593

593594
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/ZNA_Qx_vs_Centrality"), centrality, qxa);
594595
registry.fill(HIST("QA/") + HIST(Time[ft]) + HIST("/ZNA_Qy_vs_Centrality"), centrality, qya);
@@ -680,7 +681,7 @@ struct ZdcQVectors {
680681
bin = h->GetXaxis()->FindBin(TString::Format("%i", runnumber));
681682
}
682683
if (name.Contains("timestamp")) {
683-
bin = h->GetXaxis()->FindBin(rescaleTimestamp(timestamp, runnumber));
684+
bin = h->GetXaxis()->FindBin(rsTimestamp);
684685
}
685686
calibConstant = h->GetBinContent(bin);
686687
} else if (hist->InheritsFrom("THnSparse")) {
@@ -762,6 +763,7 @@ struct ZdcQVectors {
762763
registry.fill(HIST("hEventCount"), evSel_FilteredEvent);
763764

764765
timestamp = foundBC.timestamp();
766+
rsTimestamp = rescaleTimestamp(timestamp, runnumber);
765767

766768
if (!foundBC.has_zdc()) {
767769
isSelected = false;
@@ -1027,15 +1029,15 @@ struct ZdcQVectors {
10271029

10281030
if (cal.atIteration == 0) {
10291031
if (isSelected && cfgFillHistRegistry)
1030-
fillCommonRegistry<kBefore>(q[0], q[1], q[2], q[3], v, centrality);
1032+
fillCommonRegistry<kBefore>(q[0], q[1], q[2], q[3], v, centrality, rsTimestamp);
10311033

10321034
spTableZDC(runnumber, cents, v, foundBC.timestamp(), q[0], q[1], q[2], q[3], isSelected, eventSelectionFlags);
10331035
counter++;
10341036
lastRunNumber = runnumber;
10351037
return;
10361038
} else {
10371039
if (cfgFillHistRegistry)
1038-
fillCommonRegistry<kBefore>(q[0], q[1], q[2], q[3], v, centrality);
1040+
fillCommonRegistry<kBefore>(q[0], q[1], q[2], q[3], v, centrality, rsTimestamp);
10391041

10401042
// vector of 4
10411043
std::vector<double> corrQxA;
@@ -1091,15 +1093,6 @@ struct ZdcQVectors {
10911093
}
10921094
}
10931095

1094-
if (isSelected && cfgFillHistRegistry && !cfgFillNothing) {
1095-
fillCommonRegistry<kAfter>(qRec[0], qRec[1], qRec[2], qRec[3], v, centrality);
1096-
registry.fill(HIST("QA/centrality_after"), centrality);
1097-
registry.get<TProfile>(HIST("QA/after/ZNA_Qx"))->Fill(Form("%d", runnumber), qRec[0]);
1098-
registry.get<TProfile>(HIST("QA/after/ZNA_Qy"))->Fill(Form("%d", runnumber), qRec[1]);
1099-
registry.get<TProfile>(HIST("QA/after/ZNC_Qx"))->Fill(Form("%d", runnumber), qRec[2]);
1100-
registry.get<TProfile>(HIST("QA/after/ZNC_Qy"))->Fill(Form("%d", runnumber), qRec[3]);
1101-
}
1102-
11031096
// do shift for psi.
11041097
double psiZDCA = 1.0 * std::atan2(qRec[1], qRec[0]);
11051098
double psiZDCC = 1.0 * std::atan2(qRec[3], qRec[2]);
@@ -1182,6 +1175,15 @@ struct ZdcQVectors {
11821175
double qXcShift = std::hypot(qRec[2], qRec[3]) * std::cos(psiZDCCshift);
11831176
double qYcShift = std::hypot(qRec[2], qRec[3]) * std::sin(psiZDCCshift);
11841177

1178+
if (isSelected && cfgFillHistRegistry && !cfgFillNothing) {
1179+
fillCommonRegistry<kAfter>(qRec[0], qRec[1], qRec[2], qRec[3], v, centrality, rsTimestamp);
1180+
registry.fill(HIST("QA/centrality_after"), centrality);
1181+
registry.get<TProfile>(HIST("QA/after/ZNA_Qx"))->Fill(Form("%d", runnumber), qRec[0]);
1182+
registry.get<TProfile>(HIST("QA/after/ZNA_Qy"))->Fill(Form("%d", runnumber), qRec[1]);
1183+
registry.get<TProfile>(HIST("QA/after/ZNC_Qx"))->Fill(Form("%d", runnumber), qRec[2]);
1184+
registry.get<TProfile>(HIST("QA/after/ZNC_Qy"))->Fill(Form("%d", runnumber), qRec[3]);
1185+
}
1186+
11851187
spTableZDC(runnumber, cents, v, foundBC.timestamp(), qXaShift, qYaShift, qXcShift, qYcShift, isSelected, eventSelectionFlags);
11861188
qRec.clear();
11871189

‎PWGCF/Flow/Tasks/flowSP.cxx‎

Lines changed: 55 additions & 29 deletions
Original file line numberDiff line numberDiff line change
@@ -72,6 +72,7 @@ struct FlowSP {
7272
// event selection configurable group
7373
O2_DEFINE_CONFIGURABLE(cfgEvSelsUseAdditionalEventCut, bool, true, "Bool to enable Additional Event Cut");
7474
O2_DEFINE_CONFIGURABLE(cfgEvSelsMaxOccupancy, int, 10000, "Maximum occupancy of selected events");
75+
O2_DEFINE_CONFIGURABLE(cfgEvSelsMinOccupancy, int, 0, "Minimum occupancy of selected events");
7576
O2_DEFINE_CONFIGURABLE(cfgEvSelsNoSameBunchPileupCut, bool, true, "kNoSameBunchPileupCut");
7677
O2_DEFINE_CONFIGURABLE(cfgEvSelsIsGoodZvtxFT0vsPV, bool, true, "kIsGoodZvtxFT0vsPV");
7778
O2_DEFINE_CONFIGURABLE(cfgEvSelsNoCollInTimeRangeStandard, bool, true, "kNoCollInTimeRangeStandard");
@@ -88,6 +89,7 @@ struct FlowSP {
8889
O2_DEFINE_CONFIGURABLE(cfgFillEventPlaneQA, bool, false, "Fill histograms for Event Plane QA");
8990
O2_DEFINE_CONFIGURABLE(cfgFillQABefore, bool, false, "Fill QA histograms before cuts, only for processData");
9091
O2_DEFINE_CONFIGURABLE(cfgFillMeanPT, bool, false, "Fill histograms for mean PX/PT");
92+
O2_DEFINE_CONFIGURABLE(cfgUseCentAveragePt, bool, false, "Use <pt> in 1% centrality intervals and not ecent average");
9193
// Flags to make and fill histograms
9294
O2_DEFINE_CONFIGURABLE(cfgFillGeneralV1Histos, bool, true, "Fill histograms for vn analysis");
9395
O2_DEFINE_CONFIGURABLE(cfgFillMixedHarmonics, bool, true, "Flag to make and fill histos for mixed harmonics");
@@ -139,6 +141,7 @@ struct FlowSP {
139141
O2_DEFINE_CONFIGURABLE(cfgCCDB_NUE, std::string, "Users/c/ckoster/flowSP/LHC23_PbPb_pass5/NUE/Default", "ccdb dir for NUE corrections (pt)");
140142
O2_DEFINE_CONFIGURABLE(cfgCCDB_NUE2D, std::string, "Users/c/ckoster/flowSP/LHC23_PbPb_pass5/NUE/2D", "ccdb dir for NUE 2D corrections (eta, pt)");
141143
O2_DEFINE_CONFIGURABLE(cfgCCDBdir_centrality, std::string, "", "ccdb dir for Centrality corrections");
144+
O2_DEFINE_CONFIGURABLE(cfgCCDBdir_meanPt, std::string, "", "ccdb dir for Mean Pt corrections");
142145
// Confogirable axis
143146
ConfigurableAxis axisCentrality{"axisCentrality", {20, 0, 100}, "Centrality bins for vn "};
144147
ConfigurableAxis axisNch = {"axisNch", {400, 0, 4000}, "Global N_{ch}"};
@@ -187,19 +190,20 @@ struct FlowSP {
187190
TProfile* hcorrQQy = nullptr;
188191
TProfile* hEvPlaneRes = nullptr;
189192
TH1D* hCentrality = nullptr;
190-
TProfile2D* hRelPt = nullptr;
193+
TProfile2D* hMeanPt = nullptr;
191194

192195
bool clQQ = false;
193196
bool clEvPlaneRes = false;
194197
bool clCentrality = false;
195-
bool clRelPt = false;
198+
bool clMeanPt = false;
196199

197200
} cfg;
198201

199202
struct SPMvars {
200203
std::vector<std::map<int, float>> wacc = {{{0, 1.0}, {1, 1.0}, {2, 1.0}, {3, 1.0}}, {{0, 1.0}, {1, 1.0}, {2, 1.0}, {3, 1.0}}, {{0, 1.0}, {1, 1.0}, {2, 1.0}, {3, 1.0}}}; // int for part species, float for weight vector for kIncl, kPos, kNeg
201204
std::vector<std::map<int, float>> weff = {{{0, 1.0}, {1, 1.0}, {2, 1.0}, {3, 1.0}}, {{0, 1.0}, {1, 1.0}, {2, 1.0}, {3, 1.0}}, {{0, 1.0}, {1, 1.0}, {2, 1.0}, {3, 1.0}}}; // int for part species, float for weight vector for kIncl, kPos, kNeg
202205
double centWeight = 1.0;
206+
double meanPtWeight = 1.0;
203207
double ux = 0;
204208
double uy = 0;
205209
double uxMH = 0;
@@ -527,10 +531,10 @@ struct FlowSP {
527531
registry.add<TProfile2D>("incl/meanPT/meanRelPtA", "", kTProfile2D, {axisEtaVn, axisCentrality});
528532
registry.add<TProfile2D>("incl/meanPT/meanRelPtC", "", kTProfile2D, {axisEtaVn, axisCentrality});
529533

530-
registry.add<TProfile2D>("incl/meanPT/hRelEtaPt", "", kTProfile2D, {axisEtaVn, axisCentrality});
531-
registry.add<TProfile2D>("incl/meanPT/ptV1A", "", kTProfile2D, {axisEtaVn, axisCentrality});
532-
registry.add<TProfile2D>("incl/meanPT/ptV1C", "", kTProfile2D, {axisEtaVn, axisCentrality});
533-
registry.add<TProfile>("incl/meanPT/meanPT", "", kTProfile, {axisCentrality});
534+
registry.add<TProfile2D>("incl/meanPT/hMeanPtEtaCent", "", kTProfile2D, {axisEtaVn, axisCent});
535+
registry.add<TProfile2D>("incl/meanPT/ptV1A", "", kTProfile2D, {axisEtaVn, axisCent});
536+
registry.add<TProfile2D>("incl/meanPT/ptV1C", "", kTProfile2D, {axisEtaVn, axisCent});
537+
registry.add<TProfile>("incl/meanPT/hMeanPtCent", "", kTProfile, {axisCent});
534538
}
535539
if (cfgFillPID) {
536540
registry.add<TProfile3D>("incl/pion/vnC", "", kTProfile3D, {axisPt, axisEtaVn, axisCentrality});
@@ -848,7 +852,7 @@ struct FlowSP {
848852
// Occupancy
849853
if (cfgEvSelsDoOccupancySel) {
850854
auto occupancy = collision.trackOccupancyInTimeRange();
851-
if (occupancy > cfgEvSelsMaxOccupancy) {
855+
if (occupancy > cfgEvSelsMaxOccupancy || occupancy < cfgEvSelsMinOccupancy) {
852856
return 0;
853857
}
854858
histos.fill(HIST("hEventCount"), evSel_occupancy);
@@ -1071,10 +1075,10 @@ struct FlowSP {
10711075
}
10721076

10731077
if (cfgFillMeanPT) {
1074-
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("meanPT/hRelEtaPt"), track.eta(), spm.centrality, track.pt(), weight);
1075-
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("meanPT/meanPT"), spm.centrality, track.pt(), weight);
1076-
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("meanPT/ptV1A"), track.eta(), spm.centrality, track.pt() * ((spm.uy * spm.qyA + spm.ux * spm.qxA) / std::sqrt(std::fabs(spm.corrQQ))), weight);
1077-
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("meanPT/ptV1C"), track.eta(), spm.centrality, track.pt() * ((spm.uy * spm.qyC + spm.ux * spm.qxC) / std::sqrt(std::fabs(spm.corrQQ))), weight);
1078+
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("meanPT/hMeanPtEtaCent"), track.eta(), spm.centrality, track.pt(), weight);
1079+
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("meanPT/hMeanPtCent"), spm.centrality, track.pt(), weight);
1080+
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("meanPT/ptV1A"), track.eta(), spm.centrality, track.pt() * ((spm.uy * spm.qyA + spm.ux * spm.qxA) / (std::sqrt(std::fabs(spm.corrQQ)) * spm.meanPtWeight)), weight);
1081+
registry.fill(HIST(Charge[ct]) + HIST(Species[pt]) + HIST("meanPT/ptV1C"), track.eta(), spm.centrality, track.pt() * ((spm.uy * spm.qyC + spm.ux * spm.qxC) / (std::sqrt(std::fabs(spm.corrQQ)) * spm.meanPtWeight)), weight);
10781082
}
10791083
}
10801084

@@ -1305,14 +1309,19 @@ struct FlowSP {
13051309
evPlaneRes = std::sqrt(evPlaneRes);
13061310
}
13071311

1312+
spm.centWeight = 1.;
13081313
if (cfgCCDBdir_centrality.value.empty() == false) {
13091314
if (!cfg.clCentrality) {
13101315
cfg.hCentrality = ccdb->getForTimeStamp<TH1D>(cfgCCDBdir_centrality.value, bc.timestamp());
13111316
cfg.clCentrality = true;
13121317
}
1313-
spm.centWeight = cfg.hCentrality->GetBinContent(cfg.hCentrality->FindBin(spm.centrality));
1314-
if (spm.centWeight < 0)
1315-
LOGF(fatal, "Centrality weight cannot be negative.. abort for (%.2f)", spm.centrality);
1318+
double centW = cfg.hCentrality->GetBinContent(cfg.hCentrality->FindBin(spm.centrality));
1319+
if (centW < 0) {
1320+
spm.centWeight = 1. / centW;
1321+
} else {
1322+
LOGF(fatal, "Centrality weight cannot be negative .. setting to 0. for (%.2f)", spm.centrality);
1323+
spm.centWeight = 0.;
1324+
}
13161325
}
13171326

13181327
fillEventQA<kAfter>(collision, tracks);
@@ -1366,6 +1375,23 @@ struct FlowSP {
13661375
if (!trackSelected(track, field))
13671376
continue;
13681377

1378+
spm.meanPtWeight = 1.0;
1379+
if (cfgCCDBdir_meanPt.value.empty() == false) {
1380+
if (!cfg.clMeanPt) {
1381+
cfg.hMeanPt = ccdb->getForTimeStamp<TProfile2D>(cfgCCDBdir_meanPt.value, bc.timestamp());
1382+
cfg.clMeanPt = true;
1383+
}
1384+
int etaBin = cfg.hMeanPt->GetXaxis()->FindBin(track.eta());
1385+
int centBin = cfg.hMeanPt->GetYaxis()->FindBin(spm.centrality);
1386+
double weight = cfg.hMeanPt->GetBinContent(etaBin, centBin);
1387+
if (weight > 0) {
1388+
spm.meanPtWeight = 1.0 / weight;
1389+
} else {
1390+
LOGF(info, "MeanPt cannot be negative or 0.. (weight = 1/meanPt) -> setting to 0. for (%.2f, %.2f)", track.eta(), spm.centrality);
1391+
spm.meanPtWeight = 0.0;
1392+
}
1393+
}
1394+
13691395
// constrain angle to 0 -> [0,0+2pi]
13701396
auto phi = RecoDecay::constrainAngle(track.phi(), 0);
13711397

@@ -1491,20 +1517,20 @@ struct FlowSP {
14911517
double drelPxA = track.pt() * ((spm.uy * spm.qyA + spm.ux * spm.qxA) / std::sqrt(std::fabs(spm.corrQQ)));
14921518
double drelPxC = track.pt() * ((spm.uy * spm.qyC + spm.ux * spm.qxC) / std::sqrt(std::fabs(spm.corrQQ)));
14931519

1494-
meanPTMap->Fill(track.eta(), track.pt(), spm.wacc[kInclusive][kUnidentified] * spm.weff[kInclusive][kUnidentified]);
1495-
relPxA->Fill(track.eta(), drelPxA, spm.wacc[kInclusive][kUnidentified] * spm.weff[kInclusive][kUnidentified]);
1496-
relPxC->Fill(track.eta(), drelPxC, spm.wacc[kInclusive][kUnidentified] * spm.weff[kInclusive][kUnidentified]);
1520+
meanPTMap->Fill(track.eta(), track.pt(), spm.wacc[kInclusive][kUnidentified] * spm.weff[kInclusive][kUnidentified] * spm.centWeight);
1521+
relPxA->Fill(track.eta(), drelPxA, spm.wacc[kInclusive][kUnidentified] * spm.weff[kInclusive][kUnidentified] * spm.centWeight);
1522+
relPxC->Fill(track.eta(), drelPxC, spm.wacc[kInclusive][kUnidentified] * spm.weff[kInclusive][kUnidentified] * spm.centWeight);
14971523

14981524
if (spm.charge == kPositive) {
1499-
meanPTMapPos->Fill(track.eta(), track.pt(), spm.wacc[kPositive][kUnidentified] * spm.weff[kPositive][kUnidentified]);
1500-
relPxAPos->Fill(track.eta(), drelPxA, spm.wacc[kPositive][kUnidentified] * spm.weff[kPositive][kUnidentified]);
1501-
relPxCPos->Fill(track.eta(), drelPxC, spm.wacc[kPositive][kUnidentified] * spm.weff[kPositive][kUnidentified]);
1525+
meanPTMapPos->Fill(track.eta(), track.pt(), spm.wacc[kPositive][kUnidentified] * spm.weff[kPositive][kUnidentified] * spm.centWeight);
1526+
relPxAPos->Fill(track.eta(), drelPxA, spm.wacc[kPositive][kUnidentified] * spm.weff[kPositive][kUnidentified] * spm.centWeight);
1527+
relPxCPos->Fill(track.eta(), drelPxC, spm.wacc[kPositive][kUnidentified] * spm.weff[kPositive][kUnidentified] * spm.centWeight);
15021528
}
15031529

15041530
if (spm.charge == kNegative) {
1505-
meanPTMapNeg->Fill(track.eta(), track.pt(), spm.wacc[kNegative][kUnidentified] * spm.weff[kNegative][kUnidentified]);
1506-
relPxANeg->Fill(track.eta(), drelPxA, spm.wacc[kNegative][kUnidentified] * spm.weff[kNegative][kUnidentified]);
1507-
relPxCNeg->Fill(track.eta(), drelPxC, spm.wacc[kNegative][kUnidentified] * spm.weff[kNegative][kUnidentified]);
1531+
meanPTMapNeg->Fill(track.eta(), track.pt(), spm.wacc[kNegative][kUnidentified] * spm.weff[kNegative][kUnidentified] * spm.centWeight);
1532+
relPxANeg->Fill(track.eta(), drelPxA, spm.wacc[kNegative][kUnidentified] * spm.weff[kNegative][kUnidentified] * spm.centWeight);
1533+
relPxCNeg->Fill(track.eta(), drelPxC, spm.wacc[kNegative][kUnidentified] * spm.weff[kNegative][kUnidentified] * spm.centWeight);
15081534
}
15091535

15101536
} // end of track loop
@@ -1529,12 +1555,12 @@ struct FlowSP {
15291555
double drelPxCPos = relPxCPos->GetBinContent(etabin);
15301556

15311557
if (meanPt != 0) {
1532-
registry.fill(HIST("incl/meanPT/meanRelPtA"), eta, spm.centrality, drelPxA / meanPt, 1);
1533-
registry.fill(HIST("neg/meanPT/meanRelPtA"), eta, spm.centrality, drelPxANeg / meanPtNeg, 1);
1534-
registry.fill(HIST("pos/meanPT/meanRelPtA"), eta, spm.centrality, drelPxAPos / meanPtPos, 1);
1535-
registry.fill(HIST("incl/meanPT/meanRelPtC"), eta, spm.centrality, drelPxC / meanPt, 1);
1536-
registry.fill(HIST("neg/meanPT/meanRelPtC"), eta, spm.centrality, drelPxCNeg / meanPtNeg, 1);
1537-
registry.fill(HIST("pos/meanPT/meanRelPtC"), eta, spm.centrality, drelPxCPos / meanPtPos, 1);
1558+
registry.fill(HIST("incl/meanPT/meanRelPtA"), eta, spm.centrality, drelPxA / meanPt, spm.centWeight);
1559+
registry.fill(HIST("neg/meanPT/meanRelPtA"), eta, spm.centrality, drelPxANeg / meanPtNeg, spm.centWeight);
1560+
registry.fill(HIST("pos/meanPT/meanRelPtA"), eta, spm.centrality, drelPxAPos / meanPtPos, spm.centWeight);
1561+
registry.fill(HIST("incl/meanPT/meanRelPtC"), eta, spm.centrality, drelPxC / meanPt, spm.centWeight);
1562+
registry.fill(HIST("neg/meanPT/meanRelPtC"), eta, spm.centrality, drelPxCNeg / meanPtNeg, spm.centWeight);
1563+
registry.fill(HIST("pos/meanPT/meanRelPtC"), eta, spm.centrality, drelPxCPos / meanPtPos, spm.centWeight);
15381564
}
15391565
}
15401566
}

0 commit comments

Comments
 (0)