Skip to content

Commit ed09c35

Browse files
somabhattasomabhatta
committed
Update radial flow decorr task toinclude PID (2)
1 parent 53c7673 commit ed09c35

File tree

1 file changed

+23
-62
lines changed

1 file changed

+23
-62
lines changed

‎PWGCF/EbyEFluctuations/Tasks/radialFlowDecorr.cxx‎

Lines changed: 23 additions & 62 deletions
Original file line numberDiff line numberDiff line change
@@ -550,10 +550,10 @@ struct RadialFlowDecorr {
550550
aod::TracksDCA,
551551
aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr,
552552
aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr>;
553-
// Filter trackFilter = nabs(aod::track::eta) < cfgEtaCut &&
554-
// aod::track::pt > cfgPtMin&&
555-
// aod::track::pt < cfgPtMax&&
556-
// nabs(aod::track::dcaXY) < cfgDCAXY&& nabs(aod::track::dcaZ) < cfgDCAZ;
553+
Filter trackFilter = nabs(aod::track::eta) < cfgEtaCut &&
554+
aod::track::pt > cfgPtMin&&
555+
aod::track::pt < cfgPtMax&&
556+
nabs(aod::track::dcaXY) < cfgDCAXY&& nabs(aod::track::dcaZ) < cfgDCAZ;
557557
using AodTracksSel = soa::Filtered<UnfilteredTracks>;
558558
using TCs = soa::Join<UnfilteredTracks, aod::McTrackLabels>;
559559
using FilteredTCs = soa::Filtered<TCs>;
@@ -564,14 +564,16 @@ struct RadialFlowDecorr {
564564
aod::CentFT0As, aod::CentFT0Cs, aod::CentFT0Ms, aod::CentFV0As,
565565
aod::CentNGlobals, aod::McCollisionLabels>;
566566

567-
// using MyMCTracks = soa::Join<
568-
// aod::Tracks, aod::TrackSelection, aod::TracksExtra, aod::TracksDCA,
569-
// aod::McTrackLabels,
570-
// aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr,
571-
// aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr>;
567+
using MyMCTracks = soa::Join<
568+
aod::Tracks, aod::TrackSelection, aod::TracksExtra, aod::TracksDCA,
569+
aod::McTrackLabels,
570+
aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr,
571+
aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr>;
572572

573573
PresliceUnsorted<aod::McParticles> partPerMcCollision = aod::mcparticle::mcCollisionId;
574574
PresliceUnsorted<MyRun3MCCollisions> colPerMcCollision = aod::mccollisionlabel::mcCollisionId;
575+
PresliceUnsorted<TCs> trackPerMcParticle = aod::mctracklabel::mcParticleId;
576+
Preslice<MyMCTracks> perCollision = aod::track::collisionId;
575577
Preslice<FilteredTCs> trackPerCollision = aod::track::collisionId;
576578

577579
void declareCommonQA()
@@ -606,7 +608,7 @@ struct RadialFlowDecorr {
606608

607609
histos.add("h_RecoUnMatchedToPrimary", ";p_{T}", kTH1F, {{KNbinsPtRes, cfgPtMin, cfgPtMax}});
608610
}
609-
/*
611+
610612
void declareMCGetFlatHists()
611613
{
612614
for (const auto& suf : pidSuffix) {
@@ -640,12 +642,10 @@ struct RadialFlowDecorr {
640642
histos.add<TProfile3D>("pmean_cent_id_eta_FT0", ";cent;id;#eta", kTProfile3D, {{centAxis1Per}, {100, -0.5, 99.5}, {100, -5.0, 5.0}});
641643
histos.add("h3_cent_id_eta_FT0", ";cent;id;#eta", kTH3F, {{centAxis1Per}, {100, -0.5, 99.5}, {100, -5.0, 5.0}});
642644

643-
// for (const char* lvl : {"MCGen/", "MCReco/", "MCRecoEffCorr/"}) {
644645
histos.add<TProfile2D>("MCGen/Prof_Cent_Nsp_Nchrec", ";cent;isp", kTProfile2D, {{centAxis1Per}, {KNsp, -0.5, KNsp - 0.5}});
645646
histos.add<TProfile2D>("MCGen/Prof_Mult_Nsp_Nchrec", ";mult;isp", kTProfile2D, {{nChAxis}, {KNsp, -0.5, KNsp - 0.5}});
646647
histos.add<TProfile2D>("MCGen/Prof_Cent_Nsp_MeanpT", ";cent;isp", kTProfile2D, {{centAxis1Per}, {KNsp, -0.5, KNsp - 0.5}});
647648
histos.add<TProfile2D>("MCGen/Prof_Mult_Nsp_MeanpT", ";mult;isp", kTProfile2D, {{nChAxis}, {KNsp, -0.5, KNsp - 0.5}});
648-
// }
649649

650650
histos.add<TProfile3D>("pmeanTru_nch_etabin_spbin", ";mult;eta;isp", kTProfile3D, {{nChAxis}, {KNEta, 0.5, KNEta + 0.5}, {KNsp, -0.5, KNsp - 0.5}});
651651
histos.add<TProfile3D>("pmeanReco_nch_etabin_spbin", ";mult;eta;isp", kTProfile3D, {{nChAxis}, {KNEta, 0.5, KNEta + 0.5}, {KNsp, -0.5, KNsp - 0.5}});
@@ -751,7 +751,7 @@ struct RadialFlowDecorr {
751751
histos.add<TProfile3D>("Prof_CovFT0C2D_Cent_etaA_etaC" + suf, ";cent;#eta_{A};#eta_{B}", kTProfile3D, {{centAxis1Per}, {etaAxis}, {etaAxis}});
752752
}
753753
}
754-
*/
754+
755755
THnSparseF* buildWeightMapFromRaw(THnSparseF* hRaw, const char* mapName)
756756
{
757757
if (!hRaw) {
@@ -870,7 +870,6 @@ struct RadialFlowDecorr {
870870
if (cfgRunMCMean || cfgRunMCFluc || cfgRunGetEff) {
871871
declareMCCommonHists();
872872
}
873-
/*
874873
if (cfgRunMCMean) {
875874
declareMCMeanHists();
876875
histos.addClone("MCGen/", "MCReco/");
@@ -895,7 +894,7 @@ struct RadialFlowDecorr {
895894
if (cfgRunDataFluc) {
896895
declareDataFlucHists();
897896
}
898-
*/
897+
899898
if (!cfgRunGetEff && (cfgEff)) {
900899
TList* lst = ccdb->getForTimeStamp<TList>(pathEff, now);
901900

@@ -1070,7 +1069,6 @@ struct RadialFlowDecorr {
10701069

10711070
void processGetEffHists(aod::McCollisions const& mcColl, MyRun3MCCollisions const& collisions, TCs const& tracks, FilteredTCs const& /*filteredTracks*/, aod::McParticles const& mcParticles)
10721071
{
1073-
/*
10741072
for (const auto& mcCollision : mcColl) {
10751073
auto colSlice = collisions.sliceBy(colPerMcCollision, mcCollision.globalIndex());
10761074
if (colSlice.size() != 1)
@@ -1299,12 +1297,11 @@ struct RadialFlowDecorr {
12991297
} // end track loop
13001298
}
13011299
}
1302-
*/
13031300
LOGF(info, "FINISHED RUNNING processGetEffHists");
13041301
}
13051302
PROCESS_SWITCH(RadialFlowDecorr, processGetEffHists, "process MC to calculate EffWeights", cfgRunGetEff);
1306-
/*
1307-
void processMCFlat(aod::McCollisions const& mcColl, MyRun3MCCollisions const& collisions, TCs const& tracks, FilteredTCs const& , aod::McParticles const& mcParticles)
1303+
1304+
void processMCFlat(aod::McCollisions const& mcColl, MyRun3MCCollisions const& collisions, /*soa::SmallGroups<MyRun3MCCollisions> const& collisions,*/ TCs const& tracks, FilteredTCs const& /*filteredTracks*/, aod::McParticles const& mcParticles)
13081305
{
13091306
for (const auto& mcCollision : mcColl) {
13101307
auto colSlice = collisions.sliceBy(colPerMcCollision, mcCollision.globalIndex());
@@ -1398,9 +1395,8 @@ struct RadialFlowDecorr {
13981395
}
13991396
PROCESS_SWITCH(RadialFlowDecorr, processMCFlat, "process MC to calculate FlatWeights", cfgRunGetMCFlat);
14001397

1401-
void processMCMean(aod::McCollisions const& mcColl, MyRun3MCCollisions const& collisions, TCs const& tracks, FilteredTCs const& , aod::FT0s const&, aod::McParticles const& mcParticles)
1398+
void processMCMean(aod::McCollisions const& mcColl, MyRun3MCCollisions const& collisions, TCs const& tracks, FilteredTCs const& /*filteredTracks*/, aod::FT0s const&, aod::McParticles const& mcParticles)
14021399
{
1403-
// Track-sum arrays using KNsp index (isp=0: Incl, 1: Pi, 2: Ka, 3: Pr)
14041400
double sumWiTruth[KNsp][KNEta]{}, sumWiptiTruth[KNsp][KNEta]{};
14051401
double sumWiReco[KNsp][KNEta]{}, sumWiptiReco[KNsp][KNEta]{};
14061402
double sumWiRecoEffCorr[KNsp][KNEta]{}, sumWiptiRecoEffCorr[KNsp][KNEta]{};
@@ -1425,15 +1421,13 @@ struct RadialFlowDecorr {
14251421
float multPV = col.multNTracksPV();
14261422
float vz = col.posZ();
14271423

1428-
// Reset local event sum
14291424
memset(sumWiTruth, 0, sizeof(sumWiTruth));
14301425
memset(sumWiptiTruth, 0, sizeof(sumWiptiTruth));
14311426
memset(sumWiReco, 0, sizeof(sumWiReco));
14321427
memset(sumWiptiReco, 0, sizeof(sumWiptiReco));
14331428
memset(sumWiRecoEffCorr, 0, sizeof(sumWiRecoEffCorr));
14341429
memset(sumWiptiRecoEffCorr, 0, sizeof(sumWiptiRecoEffCorr));
14351430

1436-
// --- 1. Truth Loop ---
14371431
for (const auto& particle : partSlice) {
14381432
if (!isParticleSelected(particle) || !particle.isPhysicalPrimary())
14391433
continue;
@@ -1469,7 +1463,6 @@ struct RadialFlowDecorr {
14691463
}
14701464
}
14711465

1472-
// Full Event MCGen
14731466
for (int isp = 0; isp < KNsp; ++isp) {
14741467
histos.fill(HIST("MCGen/Prof_Cent_Nsp_Nchrec"), cent, isp, sumWiTruth[isp][0]);
14751468
histos.fill(HIST("MCGen/Prof_Mult_Nsp_Nchrec"), multPV, isp, sumWiTruth[isp][0]);
@@ -1479,7 +1472,6 @@ struct RadialFlowDecorr {
14791472
}
14801473
}
14811474

1482-
// --- 2. Reco Loop ---
14831475
for (const auto& track : trackSlice) {
14841476
if (!isTrackSelected(track))
14851477
continue;
@@ -1605,8 +1597,6 @@ struct RadialFlowDecorr {
16051597
for (int ietaA = 0; ietaA < KNEta; ++ietaA) {
16061598
for (int ietaC = 0; ietaC < KNEta; ++ietaC) {
16071599
for (int isp = 0; isp < KNsp; ++isp) {
1608-
1609-
// 1. Truth Sub-event Mean
16101600
float nTruAB = sumWiTruth[isp][ietaA] + sumWiTruth[isp][ietaC];
16111601
float nRecoAB = sumWiReco[isp][ietaA] + sumWiReco[isp][ietaC];
16121602
float nCorrAB = sumWiRecoEffCorr[isp][ietaA] + sumWiRecoEffCorr[isp][ietaC];
@@ -1684,7 +1674,7 @@ struct RadialFlowDecorr {
16841674
histos.fill(HIST("Prof2D_MeanpTSub_RecoEffCorr_AllPr"), cent, ietaA, ietaC, mptsubRecoEffCorr);
16851675
}
16861676

1687-
if (ietaA == ietaC) { // only fill once per eta bin
1677+
if (ietaA == ietaC) {
16881678
if (sumWiTruth[isp][ietaA] > 0) {
16891679
float val = sumWiptiTruth[isp][ietaA] / sumWiTruth[isp][ietaA];
16901680
histos.fill(HIST("pmeanTru_nch_etabin_spbin"), multPV, ietaA, isp, val);
@@ -1734,15 +1724,13 @@ struct RadialFlowDecorr {
17341724
}
17351725
PROCESS_SWITCH(RadialFlowDecorr, processMCMean, "process MC to calculate mean pt and Eff Hists", cfgRunMCMean);
17361726

1737-
void processMCFluc(aod::McCollisions const& mcColl, MyRun3MCCollisions const& collisions, TCs const& tracks, FilteredTCs const& , aod::FT0s const&, aod::McParticles const& mcParticles)
1727+
void processMCFluc(aod::McCollisions const& mcColl, MyRun3MCCollisions const& collisions, TCs const& tracks, FilteredTCs const& /*filteredTracks*/, aod::FT0s const&, aod::McParticles const& mcParticles)
17381728
{
17391729
if (!pmeanTruNchEtabinSpbinStep2 || !pmeanRecoNchEtabinSpbinStep2 || !pmeanRecoEffcorrNchEtabinSpbinStep2 ||
17401730
!pmeanMultTruNchEtabinSpbinStep2 || !pmeanMultRecoNchEtabinSpbinStep2 || !pmeanMultRecoEffcorrNchEtabinSpbinStep2) {
17411731
LOGF(warning, "MC fluc: Unified Mean pT or Mult map missing");
17421732
return;
17431733
}
1744-
1745-
// Expanded with KNsp index (isp=0: Incl, 1: Pi, 2: Ka, 3: Pr)
17461734
double sumPmwkTru[KNsp][KNEta][KIntM][KIntK]{};
17471735
double sumWkTru[KNsp][KNEta][KIntK]{};
17481736
double sumPmwkReco[KNsp][KNEta][KIntM][KIntK]{};
@@ -1779,8 +1767,6 @@ struct RadialFlowDecorr {
17791767
if (cent > KCentMax)
17801768
continue;
17811769
float multPV = col.multNTracksPV();
1782-
1783-
// Reset local arrays
17841770
memset(sumPmwkTru, 0, sizeof(sumPmwkTru));
17851771
memset(sumWkTru, 0, sizeof(sumWkTru));
17861772
memset(sumPmwkReco, 0, sizeof(sumPmwkReco));
@@ -1809,7 +1795,6 @@ struct RadialFlowDecorr {
18091795

18101796
double p1kBarFt0A = 0.0, p1kBarFt0C = 0.0;
18111797

1812-
// --- 1. Truth Loop ---
18131798
for (const auto& particle : partSlice) {
18141799
if (!isParticleSelected(particle) || !particle.isPhysicalPrimary())
18151800
continue;
@@ -1849,8 +1834,6 @@ struct RadialFlowDecorr {
18491834
}
18501835
}
18511836
} // end truth loop
1852-
1853-
// --- 2. Reco Loop ---
18541837
float vz = col.posZ();
18551838
for (const auto& track : trackSlice) {
18561839
if (!isTrackSelected(track))
@@ -2029,10 +2012,8 @@ struct RadialFlowDecorr {
20292012
}
20302013
}
20312014

2032-
// --- 3. Fill 1D Profiles: Gen, Reco, and EffCorr Levels ---
20332015
for (int ieta = 0; ieta < KNEta; ++ieta) {
20342016
for (int isp = 0; isp < KNsp; ++isp) {
2035-
// --- MCGen (Truth) ---
20362017
if (std::isfinite(meanTru[isp][ieta])) {
20372018
histos.fill(HIST("MCGen/Prof_MeanpT_Cent_etabin_spbin"), cent, ieta, isp, meanTru[isp][ieta]);
20382019
histos.fill(HIST("MCGen/Prof_MeanpT_Mult_etabin_spbin"), col.multNTracksPV(), ieta, isp, meanTru[isp][ieta]);
@@ -2041,7 +2022,6 @@ struct RadialFlowDecorr {
20412022
histos.fill(HIST("MCGen/Prof_C2_Cent_etabin_spbin"), cent, ieta, isp, c2Tru[isp][ieta]);
20422023
histos.fill(HIST("MCGen/Prof_C2_Mult_etabin_spbin"), col.multNTracksPV(), ieta, isp, c2Tru[isp][ieta]);
20432024
}
2044-
// --- MCReco ---
20452025
if (std::isfinite(meanReco[isp][ieta])) {
20462026
histos.fill(HIST("MCReco/Prof_MeanpT_Cent_etabin_spbin"), cent, ieta, isp, meanReco[isp][ieta]);
20472027
histos.fill(HIST("MCReco/Prof_MeanpT_Mult_etabin_spbin"), col.multNTracksPV(), ieta, isp, meanReco[isp][ieta]);
@@ -2050,7 +2030,6 @@ struct RadialFlowDecorr {
20502030
histos.fill(HIST("MCReco/Prof_C2_Cent_etabin_spbin"), cent, ieta, isp, c2Reco[isp][ieta]);
20512031
histos.fill(HIST("MCReco/Prof_C2_Mult_etabin_spbin"), col.multNTracksPV(), ieta, isp, c2Reco[isp][ieta]);
20522032
}
2053-
// --- MCRecoEffCorr ---
20542033
if (std::isfinite(meanRecoEffCor[isp][ieta])) {
20552034
histos.fill(HIST("MCRecoEffCorr/Prof_MeanpT_Cent_etabin_spbin"), cent, ieta, isp, meanRecoEffCor[isp][ieta]);
20562035
histos.fill(HIST("MCRecoEffCorr/Prof_MeanpT_Mult_etabin_spbin"), col.multNTracksPV(), ieta, isp, meanRecoEffCor[isp][ieta]);
@@ -2065,7 +2044,6 @@ struct RadialFlowDecorr {
20652044
p1kBarFt0A = amplFT0A - pmeanFT0AmultpvStep2->GetBinContent(pmeanFT0AmultpvStep2->GetXaxis()->FindBin(col.multNTracksPV()));
20662045
p1kBarFt0C = amplFT0C - pmeanFT0CmultpvStep2->GetBinContent(pmeanFT0CmultpvStep2->GetXaxis()->FindBin(col.multNTracksPV()));
20672046

2068-
// --- 4. Symmetric Sub-Event (1D) Covariances ---
20692047
for (int ietaA = 1; ietaA <= (KNEta - 1) / 2; ++ietaA) {
20702048
int ietaC = KNEta - ietaA;
20712049
for (int isp = 0; isp < KNsp; ++isp) {
@@ -2138,11 +2116,9 @@ struct RadialFlowDecorr {
21382116
}
21392117
}
21402118

2141-
// --- 5. Full 2D Covariances & GapSum2D Profiles ---
21422119
for (int ietaA = 1; ietaA < KNEta; ++ietaA) {
21432120
for (int ietaC = 1; ietaC < KNEta; ++ietaC) {
21442121

2145-
// Gap and Sum calculations
21462122
float etaValA = (etaLw[ietaA] + etaUp[ietaA]) / 2.0f;
21472123
float etaValB = (etaLw[ietaC] + etaUp[ietaC]) / 2.0f;
21482124
float gap = etaValA - etaValB;
@@ -2524,7 +2500,7 @@ struct RadialFlowDecorr {
25242500
}
25252501
PROCESS_SWITCH(RadialFlowDecorr, processMCFluc, "process MC to calculate pt fluc", cfgRunMCFluc);
25262502

2527-
void processGetDataFlat(AodCollisionsSel::iterator const& coll, BCsRun3 const& , aod::Zdcs const& , AodTracksSel const& tracks)
2503+
void processGetDataFlat(AodCollisionsSel::iterator const& coll, BCsRun3 const& /*bcs*/, aod::Zdcs const& /*zdcsData*/, AodTracksSel const& tracks)
25282504
{
25292505
histos.fill(HIST("hVtxZ"), coll.posZ());
25302506
if (!isEventSelected(coll))
@@ -2573,8 +2549,6 @@ struct RadialFlowDecorr {
25732549
for (int isp = 0; isp < KNsp; ++isp) {
25742550
if (!isSpecies[isp])
25752551
continue;
2576-
2577-
// Fetch efficiency specifically for this particle species
25782552
float eff = getEfficiency(coll.multNTracksPV(), pt, eta, static_cast<PID>(isp), 0, cfgEff);
25792553
float fake = getEfficiency(coll.multNTracksPV(), pt, eta, static_cast<PID>(isp), 1, cfgEff);
25802554
float w = (1.0 - fake) / eff;
@@ -2641,7 +2615,7 @@ struct RadialFlowDecorr {
26412615
}
26422616
PROCESS_SWITCH(RadialFlowDecorr, processGetDataFlat, "process data to calculate Flattening maps", cfgRunGetDataFlat);
26432617

2644-
void processDataMean(AodCollisionsSel::iterator const& coll, BCsRun3 const&, aod::Zdcs const& , aod::FT0s const&, AodTracksSel const& tracks)
2618+
void processDataMean(AodCollisionsSel::iterator const& coll, BCsRun3 const& /*bcs*/, aod::Zdcs const& /*zdcsData*/, aod::FT0s const&, AodTracksSel const& tracks)
26452619
{
26462620
double sumWi[KNsp][KNEta]{}, sumWipti[KNsp][KNEta]{};
26472621

@@ -2752,7 +2726,6 @@ struct RadialFlowDecorr {
27522726
histos.fill(HIST("hEtaPhiRecoWtd_AllPr"), vz, sign, pt, eta, phi, w);
27532727
}
27542728

2755-
// Accumulate sum
27562729
for (int ieta = 0; ieta < KNEta; ++ieta) {
27572730
if (eta <= etaLw[ieta] || eta > etaUp[ieta])
27582731
continue;
@@ -2762,15 +2735,12 @@ struct RadialFlowDecorr {
27622735
}
27632736
}
27642737

2765-
// Full Event Means
27662738
for (int isp = 0; isp < KNsp; ++isp) {
27672739
histos.fill(HIST("Prof_Cent_Nsp_Nchrec"), cent, isp, sumWi[isp][0]);
27682740
histos.fill(HIST("Prof_Mult_Nsp_Nchrec"), coll.multNTracksPV(), isp, sumWi[isp][0]);
27692741
if (sumWi[isp][0] > 1.0f)
27702742
histos.fill(HIST("Prof_Cent_Nsp_MeanpT"), cent, isp, sumWipti[isp][0] / sumWi[isp][0]);
27712743
}
2772-
2773-
// Kinematic Bin Means (1D and 2D Sub-event)
27742744
for (int ietaA = 0; ietaA < KNEta; ++ietaA) {
27752745
for (int ietaC = 0; ietaC < KNEta; ++ietaC) {
27762746
for (int isp = 0; isp < KNsp; ++isp) {
@@ -2840,7 +2810,7 @@ struct RadialFlowDecorr {
28402810
}
28412811
PROCESS_SWITCH(RadialFlowDecorr, processDataMean, "process data to calculate mean pT", cfgRunDataMean);
28422812

2843-
void processDataFluc(AodCollisionsSel::iterator const& coll, BCsRun3 const& , aod::Zdcs const& , aod::FT0s const&, AodTracksSel const& tracks)
2813+
void processDataFluc(AodCollisionsSel::iterator const& coll, BCsRun3 const& /*bcs*/, aod::Zdcs const& /*zdcsData*/, aod::FT0s const&, AodTracksSel const& tracks)
28442814
{
28452815
if (!isEventSelected(coll))
28462816
return;
@@ -2853,7 +2823,6 @@ struct RadialFlowDecorr {
28532823
return;
28542824
}
28552825

2856-
// 2. Safety Check: Correction Maps (Looping over Inclusive, Pi, Ka, Pr)
28572826
for (int isp = 0; isp < KNsp; ++isp) {
28582827
auto pid = static_cast<PID>(isp);
28592828
if (!hEff[pid] || !hFake[pid] || !hFlatWeight[pid]) {
@@ -2862,7 +2831,6 @@ struct RadialFlowDecorr {
28622831
}
28632832
}
28642833

2865-
// Expanded arrays to handle KNsp species (0: Incl, 1: Pi, 2: Ka, 3: Pr)
28662834
double sumpmwk[KNsp][KNEta][KIntM][KIntK]{};
28672835
double sumwk[KNsp][KNEta][KIntK]{};
28682836

@@ -2916,8 +2884,6 @@ struct RadialFlowDecorr {
29162884

29172885
if (!std::isfinite(w) || w <= KFloatEpsilon || eff <= KFloatEpsilon)
29182886
continue;
2919-
2920-
// Kinematic Bin sum
29212887
for (int ieta = 0; ieta < KNEta; ++ieta) {
29222888
if (eta <= etaLw[ieta] || eta > etaUp[ieta])
29232889
continue;
@@ -2946,7 +2912,6 @@ struct RadialFlowDecorr {
29462912
double p1kBarFt0A = amplFT0A - pmeanFT0AmultpvStep2->GetBinContent(pmeanFT0AmultpvStep2->GetXaxis()->FindBin(coll.multNTracksPV()));
29472913
double p1kBarFt0C = amplFT0C - pmeanFT0CmultpvStep2->GetBinContent(pmeanFT0CmultpvStep2->GetXaxis()->FindBin(coll.multNTracksPV()));
29482914

2949-
// --- 2. Step 2 Means and 1D Fluc Variables ---
29502915
for (int ieta = 0; ieta < KNEta; ++ieta) {
29512916
const int ibx = pmeanNchEtabinSpbinStep2->GetXaxis()->FindBin(coll.multNTracksPV());
29522917
const int iby = ieta + 1;
@@ -2967,8 +2932,6 @@ struct RadialFlowDecorr {
29672932
p1kBarMult[isp][ieta] = meanMult[isp][ieta] - mmMult;
29682933
}
29692934
}
2970-
2971-
// --- 3. Fill 1D Profiles ---
29722935
for (int ieta = 0; ieta < KNEta; ++ieta) {
29732936
for (int isp = 0; isp < KNsp; ++isp) {
29742937
if (std::isfinite(mean[isp][ieta])) {
@@ -3014,8 +2977,6 @@ struct RadialFlowDecorr {
30142977
}
30152978
}
30162979
}
3017-
3018-
// --- 5. Full 2D Covariances & GapSum Profiles ---
30192980
for (int ietaA = 1; ietaA < KNEta; ++ietaA) {
30202981
for (int ietaC = 1; ietaC < KNEta; ++ietaC) {
30212982

@@ -3146,7 +3107,7 @@ struct RadialFlowDecorr {
31463107
}
31473108
}
31483109
PROCESS_SWITCH(RadialFlowDecorr, processDataFluc, "process data to calculate fluc pT", cfgRunDataFluc);
3149-
*/
3110+
31503111
};
31513112

31523113
WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)

0 commit comments

Comments
 (0)