@@ -233,15 +233,11 @@ struct RadialFlowDecorr {
233233
234234 TProfile3D* pmeanTruNchEtabinPtbinStep2 = nullptr ;
235235 TProfile3D* pmeanRecoNchEtabinPtbinStep2 = nullptr ;
236- // TProfile3D* pmeanRecoMatchedNchEtabinPtbinStep2 = nullptr;
237236 TProfile3D* pmeanRecoEffcorrNchEtabinPtbinStep2 = nullptr ;
238- // TProfile3D* pmeanRecoMatchedEffcorrNchEtabinPtbinStep2 = nullptr;
239237
240238 TProfile3D* pmeanEtTruNchEtabinPtbinStep2 = nullptr ;
241239 TProfile3D* pmeanEtRecoNchEtabinPtbinStep2 = nullptr ;
242- // TProfile3D* pmeanEtRecoMatchedNchEtabinPtbinStep2 = nullptr;
243240 TProfile3D* pmeanEtRecoEffcorrNchEtabinPtbinStep2 = nullptr ;
244- // TProfile3D* pmeanEtRecoMatchedEffcorrNchEtabinPtbinStep2 = nullptr;
245241
246242 TProfile3D* pmeanMultTruNchEtabinPtbinStep2 = nullptr ;
247243 TProfile3D* pmeanMultRecoNchEtabinPtbinStep2 = nullptr ;
@@ -310,7 +306,6 @@ struct RadialFlowDecorr {
310306 auto * pd = pdg->GetParticle (particle.pdgCode ());
311307 if (!pd)
312308 return false ;
313- // if (dpt::isStrangeBaryonPDG(particle.pdgCode())) return false;
314309 if (std::abs (pd->Charge ()) == 0 )
315310 return false ;
316311 if (particle.pt () < cfgCutPtLower || particle.pt () > cfgCutPtUpper || std::abs (particle.eta ()) > cfgCutEta)
@@ -655,18 +650,14 @@ struct RadialFlowDecorr {
655650
656651 histos.add <TProfile3D>(" pmeanTruNchEtabinPtbin" , " ;N_{PV};#eta; p_{T}" , kTProfile3D , {{nChAxis}, {KNEta + 1 , -KBinOffset, KNEta + KBinOffset}, {KNpT + 1 , -KBinOffset, KNpT + KBinOffset}});
657652 histos.add <TProfile3D>(" pmeanRecoNchEtabinPtbin" , " ;N_{PV};#eta; p_{T}" , kTProfile3D , {{nChAxis}, {KNEta + 1 , -KBinOffset, KNEta + KBinOffset}, {KNpT + 1 , -KBinOffset, KNpT + KBinOffset}});
658- // histos.add<TProfile3D>("pmeanRecoMatchedNchEtabinPtbin", ";N_{PV};#eta; p_{T}", kTProfile3D, {{nChAxis}, {KNEta + 1, -KBinOffset, KNEta + KBinOffset}, {KNpT + 1, -KBinOffset, KNpT + KBinOffset}});
659653 histos.add <TProfile3D>(" pmeanRecoEffcorrNchEtabinPtbin" , " ;N_{PV};#eta; p_{T}" , kTProfile3D , {{nChAxis}, {KNEta + 1 , -KBinOffset, KNEta + KBinOffset}, {KNpT + 1 , -KBinOffset, KNpT + KBinOffset}});
660- // histos.add<TProfile3D>("pmeanRecoMatchedEffcorrNchEtabinPtbin", ";N_{PV};#eta; p_{T}", kTProfile3D, {{nChAxis}, {KNEta + 1, -KBinOffset, KNEta + KBinOffset}, {KNpT + 1, -KBinOffset, KNpT + KBinOffset}});
661654
662655 histos.add (" MCGen/Prof_Cent_MeanEt" , " ;cent;#LT E_{T}#GT" , kTProfile , {centAxis1Per});
663656 histos.add (" MCGen/Prof_Mult_MeanEt" , " ;N_{PV};#LT E_{T}#GT" , kTProfile , {nChAxis});
664657
665658 histos.add <TProfile3D>(" pmeanEtTruNchEtabinPtbin" , " ;N_{PV};#eta; p_{T}" , kTProfile3D , {{nChAxis}, {KNEta + 1 , -KBinOffset, KNEta + KBinOffset}, {KNpT + 1 , -KBinOffset, KNpT + KBinOffset}});
666659 histos.add <TProfile3D>(" pmeanEtRecoNchEtabinPtbin" , " ;N_{PV};#eta; p_{T}" , kTProfile3D , {{nChAxis}, {KNEta + 1 , -KBinOffset, KNEta + KBinOffset}, {KNpT + 1 , -KBinOffset, KNpT + KBinOffset}});
667- // histos.add<TProfile3D>("pmeanEtRecoMatchedNchEtabinPtbin", ";N_{PV};#eta; p_{T}", kTProfile3D, {{nChAxis}, {KNEta + 1, -KBinOffset, KNEta + KBinOffset}, {KNpT + 1, -KBinOffset, KNpT + KBinOffset}});
668660 histos.add <TProfile3D>(" pmeanEtRecoEffcorrNchEtabinPtbin" , " ;N_{PV};#eta; p_{T}" , kTProfile3D , {{nChAxis}, {KNEta + 1 , -KBinOffset, KNEta + KBinOffset}, {KNpT + 1 , -KBinOffset, KNpT + KBinOffset}});
669- // histos.add<TProfile3D>("pmeanEtRecoMatchedEffcorrNchEtabinPtbin", ";N_{PV};#eta; p_{T}", kTProfile3D, {{nChAxis}, {KNEta + 1, -KBinOffset, KNEta + KBinOffset}, {KNpT + 1, -KBinOffset, KNpT + KBinOffset}});
670661
671662 histos.add <TProfile3D>(" pmeanMultTruNchEtabinPtbin" , " ;N_{PV};#eta; p_{T}" , kTProfile3D , {{nChAxis}, {KNEta + 1 , -KBinOffset, KNEta + KBinOffset}, {KNpT + 1 , -KBinOffset, KNpT + KBinOffset}});
672663 histos.add <TProfile3D>(" pmeanMultRecoNchEtabinPtbin" , " ;N_{PV};#eta; p_{T}" , kTProfile3D , {{nChAxis}, {KNEta + 1 , -KBinOffset, KNEta + KBinOffset}, {KNpT + 1 , -KBinOffset, KNpT + KBinOffset}});
@@ -801,7 +792,6 @@ struct RadialFlowDecorr {
801792 }
802793 auto hWMap = reinterpret_cast <THnSparseF*>(hRaw->Clone (mapName));
803794 hWMap->SetTitle (Form (" Flattening Weight Map %s (w_{#phi} = <N_{#phi}> / N_{#phi})" , mapName));
804- // hWMap->SetDirectory(nullptr);
805795 hWMap->Reset ();
806796 auto axV = hRaw->GetAxis (0 ); // Vz
807797 auto axChg = hRaw->GetAxis (1 ); // Charge
@@ -922,16 +912,12 @@ struct RadialFlowDecorr {
922912 if (cfgRunMCMean) {
923913 declareMCMeanHists ();
924914 histos.addClone (" MCGen/" , " MCReco/" );
925- // histos.addClone("MCGen/", "MCRecoMatched/");
926915 histos.addClone (" MCGen/" , " MCRecoEffCorr/" );
927- // histos.addClone("MCGen/", "MCRecoMatchedEffCorr/");
928916 }
929917 if (cfgRunMCFluc) {
930918 declareMCFlucHists ();
931919 histos.addClone (" MCGen/" , " MCReco/" );
932- // histos.addClone("MCGen/", "MCRecoMatched/");
933920 histos.addClone (" MCGen/" , " MCRecoEffCorr/" );
934- // histos.addClone("MCGen/", "MCRecoMatchedEffCorr/");
935921 }
936922 if (cfgRunGetDataFlat) {
937923 declareDataGetFlatHists ();
@@ -1066,15 +1052,11 @@ struct RadialFlowDecorr {
10661052 if (lstMCMean) {
10671053 loadTProfile3DFromList (lstMCMean, " pmeanTruNchEtabinPtbin" , pmeanTruNchEtabinPtbinStep2);
10681054 loadTProfile3DFromList (lstMCMean, " pmeanRecoNchEtabinPtbin" , pmeanRecoNchEtabinPtbinStep2);
1069- // loadTProfile3DFromList(lstMCMean, "pmeanRecoMatchedNchEtabinPtbin", pmeanRecoMatchedNchEtabinPtbinStep2);
10701055 loadTProfile3DFromList (lstMCMean, " pmeanRecoEffcorrNchEtabinPtbin" , pmeanRecoEffcorrNchEtabinPtbinStep2);
1071- // loadTProfile3DFromList(lstMCMean, "pmeanRecoMatchedEffcorrNchEtabinPtbin", pmeanRecoMatchedEffcorrNchEtabinPtbinStep2);
10721056
10731057 loadTProfile3DFromList (lstMCMean, " pmeanEtTruNchEtabinPtbin" , pmeanEtTruNchEtabinPtbinStep2);
10741058 loadTProfile3DFromList (lstMCMean, " pmeanEtRecoNchEtabinPtbin" , pmeanEtRecoNchEtabinPtbinStep2);
1075- // loadTProfile3DFromList(lstMCMean, "pmeanEtRecoMatchedNchEtabinPtbin", pmeanEtRecoMatchedNchEtabinPtbinStep2);
10761059 loadTProfile3DFromList (lstMCMean, " pmeanEtRecoEffcorrNchEtabinPtbin" , pmeanEtRecoEffcorrNchEtabinPtbinStep2);
1077- // loadTProfile3DFromList(lstMCMean, "pmeanEtRecoMatchedEffcorrNchEtabinPtbin", pmeanEtRecoMatchedEffcorrNchEtabinPtbinStep2);
10781060
10791061 loadTProfile3DFromList (lstMCMean, " pmeanMultTruNchEtabinPtbin" , pmeanMultTruNchEtabinPtbinStep2);
10801062 loadTProfile3DFromList (lstMCMean, " pmeanMultRecoNchEtabinPtbin" , pmeanMultRecoNchEtabinPtbinStep2);
@@ -1328,7 +1310,6 @@ struct RadialFlowDecorr {
13281310 histos.fill (HIST (" TruthTracKVz" ), mcCollision.posZ (), col.posZ ());
13291311 histos.fill (HIST (" vzResolution" ), mcCollision.posZ (), mcCollision.posZ () - col.posZ ());
13301312
1331- // LOGF(info, "Event Check: cent = %.1f, nTracks = %d", cent, (int)trackSlice.size());
13321313 memset (sumWiTruth, 0 , sizeof (sumWiTruth));
13331314 memset (sumWiptiTruth, 0 , sizeof (sumWiptiTruth));
13341315 memset (sumWiReco, 0 , sizeof (sumWiReco));
@@ -1612,7 +1593,6 @@ struct RadialFlowDecorr {
16121593 auto colSlice = collisions.sliceBy (colPerMcCollision, mcCollision.globalIndex ());
16131594 if (colSlice.size () != 1 )
16141595 continue ;
1615- // histos.fill(HIST("MCGen/hVtxZ"), mcCollision.posZ());
16161596 for (const auto & col : colSlice) {
16171597
16181598 auto trackSlice = tracks.sliceBy (trackPerCollision, col.globalIndex ());
0 commit comments