Skip to content

Commit c969012

Browse files
ilikmetaalibuild
ilikmeta
and
alibuild
authored
[PWGCF] Fix FT0 low cut (AliceO2Group#9779)
Co-authored-by: ALICE Action Bot <[email protected]>
1 parent 28d6501 commit c969012

File tree

1 file changed

+67
-73
lines changed

1 file changed

+67
-73
lines changed

‎PWGCF/Flow/Tasks/flowGfwTask.cxx‎

Lines changed: 67 additions & 73 deletions
Original file line numberDiff line numberDiff line change
@@ -61,27 +61,28 @@ struct FlowGfwTask {
6161
O2_DEFINE_CONFIGURABLE(cfgCutChi2prTPCcls, float, 2.5, "Chi2 per TPC clusters")
6262
O2_DEFINE_CONFIGURABLE(cfgCutTPCclu, float, 70.0f, "minimum TPC clusters")
6363
O2_DEFINE_CONFIGURABLE(cfgCutITSclu, float, 5.0f, "minimum ITS clusters")
64-
O2_DEFINE_CONFIGURABLE(cfgUseAdditionalEventCut, bool, false, "Use additional event cut on mult correlations")
65-
O2_DEFINE_CONFIGURABLE(cfgUseAdditionalTrackCut, bool, false, "Use additional track cut on phi")
66-
O2_DEFINE_CONFIGURABLE(cfgUseNch, bool, false, "Use Nch for flow observables")
64+
O2_DEFINE_CONFIGURABLE(cfgMinCentFT0C, float, 0.0f, "Minimum FT0C Centrality")
65+
O2_DEFINE_CONFIGURABLE(cfgMaxCentFT0C, float, 100.0f, "Maximum FT0C Centrality")
66+
O2_DEFINE_CONFIGURABLE(cfgUseAdditionalEventCut, bool, true, "Use additional event cut on mult correlations")
67+
O2_DEFINE_CONFIGURABLE(cfgUseAdditionalTrackCut, bool, true, "Use additional track cut on phi")
68+
O2_DEFINE_CONFIGURABLE(cfgUseNch, bool, true, "Use Nch for flow observables")
6769
O2_DEFINE_CONFIGURABLE(cfgNbootstrap, int, 10, "Number of subsamples")
68-
O2_DEFINE_CONFIGURABLE(cfgOutputNUAWeights, bool, false, "Fill and output NUA weights")
70+
O2_DEFINE_CONFIGURABLE(cfgOutputNUAWeights, bool, true, "Fill and output NUA weights")
6971
O2_DEFINE_CONFIGURABLE(cfgEfficiency, std::string, "", "CCDB path to efficiency object")
7072
O2_DEFINE_CONFIGURABLE(cfgAcceptance, std::string, "", "CCDB path to acceptance object")
7173
O2_DEFINE_CONFIGURABLE(cfgMagnetField, std::string, "GLO/Config/GRPMagField", "CCDB path to Magnet field object")
72-
O2_DEFINE_CONFIGURABLE(cfgCutOccupancyHigh, int, 500, "High cut on TPC occupancy")
74+
O2_DEFINE_CONFIGURABLE(cfgCutOccupancyHigh, int, 3000, "High cut on TPC occupancy")
7375
O2_DEFINE_CONFIGURABLE(cfgCutOccupancyLow, int, 0, "Low cut on TPC occupancy")
7476
O2_DEFINE_CONFIGURABLE(cfgCutDCAz, float, 2, "Custom DCA Z cut")
7577
O2_DEFINE_CONFIGURABLE(cfgCutDCAxy, float, 0.2f, "Custom DCA XY cut")
76-
O2_DEFINE_CONFIGURABLE(cfgTVXinTRD, bool, false, "Use kTVXinTRD (reject TRD triggered events)");
77-
O2_DEFINE_CONFIGURABLE(cfgNoTimeFrameBorder, bool, false, "kNoTimeFrameBorder");
78-
O2_DEFINE_CONFIGURABLE(cfgNoITSROFrameBorder, bool, false, "kNoITSROFrameBorder");
79-
O2_DEFINE_CONFIGURABLE(cfgNoSameBunchPileup, bool, false, "kNoSameBunchPileup");
80-
O2_DEFINE_CONFIGURABLE(cfgIsGoodZvtxFT0vsPV, bool, false, "kIsGoodZvtxFT0vsPV");
81-
O2_DEFINE_CONFIGURABLE(cfgNoCollInTimeRangeStandard, bool, false, "kNoCollInTimeRangeStandard");
82-
O2_DEFINE_CONFIGURABLE(cfgEvSelkIsGoodITSLayersAll, bool, false, "kIsGoodITSLayersAll")
83-
O2_DEFINE_CONFIGURABLE(cfgOccupancy, bool, false, "Bool for event selection on detector occupancy");
84-
O2_DEFINE_CONFIGURABLE(cfgMultCut, bool, false, "Use additional event cut on mult correlations");
78+
O2_DEFINE_CONFIGURABLE(cfgNoTimeFrameBorder, bool, true, "kNoTimeFrameBorder");
79+
O2_DEFINE_CONFIGURABLE(cfgNoITSROFrameBorder, bool, true, "kNoITSROFrameBorder");
80+
O2_DEFINE_CONFIGURABLE(cfgNoSameBunchPileup, bool, true, "kNoSameBunchPileup");
81+
O2_DEFINE_CONFIGURABLE(cfgIsGoodZvtxFT0vsPV, bool, true, "kIsGoodZvtxFT0vsPV");
82+
O2_DEFINE_CONFIGURABLE(cfgNoCollInTimeRangeStandard, bool, true, "kNoCollInTimeRangeStandard");
83+
O2_DEFINE_CONFIGURABLE(cfgEvSelkIsGoodITSLayersAll, bool, true, "kIsGoodITSLayersAll")
84+
O2_DEFINE_CONFIGURABLE(cfgOccupancy, bool, true, "Bool for event selection on detector occupancy");
85+
O2_DEFINE_CONFIGURABLE(cfgMultCut, bool, true, "Use additional event cut on mult correlations");
8586
O2_DEFINE_CONFIGURABLE(FineBinning, bool, false, "Manually change to fine binning")
8687
O2_DEFINE_CONFIGURABLE(cfgTrackSelRun3ITSMatch, bool, false, "System check: Run3ITSMatch")
8788
O2_DEFINE_CONFIGURABLE(cfgTrackSel, bool, false, "System check: track selection")
@@ -98,9 +99,9 @@ struct FlowGfwTask {
9899
ConfigurableAxis axisT0C{"axisT0C", {70, 0, 70000}, "N_{ch} (T0C)"};
99100
ConfigurableAxis axisT0A{"axisT0A", {200, 0, 200000}, "N_{ch} (T0A)"};
100101
ConfigurableAxis axisT0M{"axisT0M", {70, 0, 70000}, "N_{ch} (T0M)"};
101-
ConfigurableAxis axisFT0CAmp{"axisFT0CAmp", {5000, 0, 5000}, "axisFT0CAmp"};
102-
ConfigurableAxis axisFT0AAmp{"axisFT0AAmp", {5000, 0, 5000}, "axisFT0AAmp"};
103-
ConfigurableAxis axisFT0MAmp{"axisFT0MAmp", {10000, 0, 10000}, "axisFT0MAmp"};
102+
ConfigurableAxis axisFT0CAmp{"axisFT0CAmp", {50000, 0, 50000}, "axisFT0CAmp"};
103+
ConfigurableAxis axisFT0AAmp{"axisFT0AAmp", {50000, 0, 50000}, "axisFT0AAmp"};
104+
ConfigurableAxis axisFT0MAmp{"axisFT0MAmp", {50000, 0, 50000}, "axisFT0MAmp"};
104105
ConfigurableAxis axisNchPV{"axisNchPV", {4000, 0, 4000}, "N_{ch} (PV)"};
105106
ConfigurableAxis axisDCAz{"axisDCAz", {200, -2, 2}, "DCA_{z} (cm)"};
106107
ConfigurableAxis axisDCAxy{"axisDCAxy", {200, -1, 1}, "DCA_{xy} (cm)"};
@@ -111,7 +112,7 @@ struct FlowGfwTask {
111112
Configurable<float> maxZP{"maxZP", 3099.5, "Max ZP signal"};
112113
Configurable<float> maxZEM{"maxZEM", 3099.5, "Max ZEM signal"};
113114
Configurable<int> nBinsFit{"nBinsFit", 1000, "nbinsFit"};
114-
Configurable<float> maxMultFT0{"maxMultFT0", 3000, "Max FT0 signal"};
115+
Configurable<float> maxMultFT0{"maxMultFT0", 5000, "Max FT0 signal"};
115116

116117
// Corrections
117118
TH1D* mEfficiency = nullptr;
@@ -171,7 +172,6 @@ struct FlowGfwTask {
171172
kFILTERED,
172173
kSEL8,
173174
kOCCUPANCY,
174-
kTVXINTRD,
175175
kNOTIMEFRAMEBORDER,
176176
kNOITSROFRAMEBORDER,
177177
kNOPSAMEBUNCHPILEUP,
@@ -206,7 +206,6 @@ struct FlowGfwTask {
206206
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(kFILTERED + 1, "Filtered events");
207207
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(kSEL8 + 1, "Sel8");
208208
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(kOCCUPANCY + 1, "Occupancy");
209-
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(kTVXINTRD + 1, "kTVXinTRD");
210209
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(kNOTIMEFRAMEBORDER + 1, "kNoTimeFrameBorder");
211210
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(kNOITSROFRAMEBORDER + 1, "kNoITSROFrameBorder");
212211
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(kNOPSAMEBUNCHPILEUP + 1, "kNoSameBunchPileup");
@@ -257,14 +256,16 @@ struct FlowGfwTask {
257256
registry.add("ZPAcoll", "ZPAcoll; ZPA amplitude; Entries", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZP}}});
258257
registry.add("ZNCcoll", "ZNCcoll; ZNC amplitude; Entries", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZN}}});
259258
registry.add("ZPCcoll", "ZPCcoll; ZPC amplitude; Entries", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZP}}});
260-
registry.add("ZNvsFT0correl", "ZNvsFT0correl", {HistType::kTH2F, {{{nBinsFit, 0., maxMultFT0}, {nBinsAmp, -0.5, 2. * maxZN}}}});
259+
registry.add("ZNvsFT0correl", "ZNvsFT0correl; FT0 amplitude; ZN", {HistType::kTH2F, {{{nBinsFit, 0., maxMultFT0}, {nBinsAmp, -0.5, 2. * maxZN}}}});
261260
registry.add("ZDCAmp", "ZDC Amplitude; ZDC Amplitude; Events", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZP}}});
262-
registry.add("ZNvsZEMcoll", "ZNvsZEMcoll; ZEM; ZNA+ZNC", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
263-
registry.add("ZNvsZEMcoll05", "ZNvsZEMcoll; ZEM; ZNA+ZNC", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
264-
registry.add("ZNvsZEMcoll510", "ZNvsZEMcoll; ZEM; ZNA+ZNC", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
265-
registry.add("ZNvsZEMcoll1020", "ZNvsZEMcoll; ZEM; ZNA+ZNC", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
266-
registry.add("ZNvsZEMcoll2030", "ZNvsZEMcoll; ZEM; ZNA+ZNC", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
267-
registry.add("ZNvsZEMcollrest", "ZNvsZEMcoll; ZEM; ZNA+ZNC", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
261+
registry.add("ZNAmp", "ZNA+ZNC Amplitude; ZN Amplitude; Events", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZN}}});
262+
registry.add("ZPAmp", "ZPA+ZPC Amplitude; ZP Amplitude; Events", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZP}}});
263+
registry.add("ZNvsZEMcoll", "ZNvsZEMcoll; ZEM; ZDC energy (GeV)", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
264+
registry.add("ZNvsZEMcoll05", "ZNvsZEMcoll; ZEM; ZDC energy (GeV)", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
265+
registry.add("ZNvsZEMcoll510", "ZNvsZEMcoll; ZEM; ZDC energy (GeV)", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
266+
registry.add("ZNvsZEMcoll1020", "ZNvsZEMcoll; ZEM; ZDC energy (GeV)", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
267+
registry.add("ZNvsZEMcoll2030", "ZNvsZEMcoll; ZEM; ZDC energy (GeV)", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
268+
registry.add("ZNvsZEMcollrest", "ZNvsZEMcoll; ZEM; ZDC energy (GeV)", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZEM}, {nBinsAmp, -0.5, 2. * maxZN}}}});
268269

269270
// Track plots
270271
registry.add("Events_per_Centrality_Bin", "Events_per_Centrality_Bin;Centrality FT0C;No. of Events", kTH1F, {axisCentrality});
@@ -368,7 +369,7 @@ struct FlowGfwTask {
368369
fWeights->init(true, false);
369370
}
370371

371-
// add in FlowContainer to Get boostrap sample automatically
372+
// add in FlowContainer to Get boostrap sample automatically -- Use post process flow task
372373
TObjArray* oba = new TObjArray();
373374
fFC->SetXAxis(fPtAxis);
374375
fFC->SetName("FlowContainer");
@@ -518,13 +519,6 @@ struct FlowGfwTask {
518519
template <typename TCollision>
519520
bool eventSelected(o2::aod::mult::MultNTracksPV, TCollision collision, const int multTrk, const float centrality)
520521
{
521-
if (cfgTVXinTRD) {
522-
if (collision.alias_bit(kTVXinTRD)) {
523-
// TRD triggered
524-
return false;
525-
}
526-
registry.fill(HIST("hEventCount"), kTVXINTRD);
527-
}
528522
if (cfgNoTimeFrameBorder) {
529523
if (!collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
530524
// reject collisions close to Time Frame borders
@@ -566,7 +560,7 @@ struct FlowGfwTask {
566560
}
567561
if (cfgEvSelkIsGoodITSLayersAll) {
568562
if (cfgEvSelkIsGoodITSLayersAll && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
569-
// no other collisions in this Readout Frame with per-collision multiplicity above threshold
563+
// removes dead staves of ITS
570564
return false;
571565
}
572566
registry.fill(HIST("hEventCount"), kISGOODITSLAYERSALL);
@@ -650,7 +644,7 @@ struct FlowGfwTask {
650644
}
651645

652646
// Apply process filters
653-
Filter collisionFilter = nabs(aod::collision::posZ) < cfgCutVertex;
647+
Filter collisionFilter = nabs(aod::collision::posZ) < cfgCutVertex && (aod::cent::centFT0C > cfgMinCentFT0C) && (aod::cent::centFT0C < cfgMaxCentFT0C);
654648
Filter trackFilter = (nabs(aod::track::eta) < cfgCutEta) && (aod::track::pt > cfgCutPtMin) && (aod::track::pt < cfgCutPtMax) && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t) true)) && (aod::track::tpcChi2NCl < cfgCutChi2prTPCcls) && (nabs(aod::track::dcaZ) < cfgCutDCAz) && (nabs(aod::track::dcaXY) < cfgCutDCAxy);
655649

656650
using Colls = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0Cs, aod::CentFT0As, aod::CentFT0Ms>>; // collisions filter
@@ -668,6 +662,8 @@ struct FlowGfwTask {
668662
if (nTotal < 1)
669663
return;
670664

665+
const auto centrality = collision.centFT0C();
666+
671667
// fill event QA before cuts
672668
registry.fill(HIST("BeforeCut_globalTracks_centT0C"), collision.centFT0C(), tracks.size());
673669
registry.fill(HIST("BeforeCut_PVTracks_centT0C"), collision.centFT0C(), collision.multNTracksPV());
@@ -680,8 +676,6 @@ struct FlowGfwTask {
680676
registry.fill(HIST("BeforeCut_multFT0M_centFT0M"), collision.centFT0M(), collision.multFT0M());
681677
registry.fill(HIST("hEventCount"), kSEL8);
682678

683-
const auto centrality = collision.centFT0C();
684-
685679
if (cfgOccupancy) {
686680
int occupancy = collision.trackOccupancyInTimeRange();
687681
if (occupancy < cfgCutOccupancyLow || occupancy > cfgCutOccupancyHigh)
@@ -693,47 +687,20 @@ struct FlowGfwTask {
693687
return;
694688
}
695689

696-
if (centrality < 0 || centrality >= 70.)
697-
return;
698-
699-
float vtxz = collision.posZ();
700-
float lRandom = fRndm->Rndm();
701-
registry.fill(HIST("hVtxZ"), vtxz);
702-
registry.fill(HIST("hMult"), nTotal);
703-
registry.fill(HIST("hCent"), centrality);
704-
registry.fill(HIST("cent_vs_Nch"), centrality, nTotal);
705-
706-
fGFW->Clear();
707-
708-
auto bc = collision.bc_as<BCsRun3>();
709-
loadCorrections(bc.timestamp());
710-
registry.fill(HIST("hEventCount"), kCENTRALITY);
711-
712-
// fill event QA after cuts
713-
registry.fill(HIST("globalTracks_centT0C_Aft"), collision.centFT0C(), tracks.size());
714-
registry.fill(HIST("PVTracks_centT0C_Aft"), collision.centFT0C(), collision.multNTracksPV());
715-
registry.fill(HIST("globalTracks_PVTracks_Aft"), collision.multNTracksPV(), tracks.size());
716-
registry.fill(HIST("globalTracks_multT0A_Aft"), collision.multFT0A(), tracks.size());
717-
registry.fill(HIST("globalTracks_multV0A_Aft"), collision.multFV0A(), tracks.size());
718-
registry.fill(HIST("multV0A_multT0A_Aft"), collision.multFT0A(), collision.multFV0A());
719-
registry.fill(HIST("multT0C_centT0C_Aft"), collision.centFT0C(), collision.multFT0C());
720-
registry.fill(HIST("multT0A_centT0A_Aft"), collision.centFT0A(), collision.multFT0A());
721-
registry.fill(HIST("multFT0M_centFT0M_Aft"), collision.centFT0M(), collision.multFT0M());
722-
723690
const auto& foundBC = collision.foundBC_as<BCsRun3>();
724691
if (foundBC.has_zdc()) {
725692
registry.fill(HIST("hEventCounterForZDC"), 1);
726693

727694
// FT0 amplitude to use in fine binning
728695
double ft0aAmp = 0;
729696
double ft0cAmp = 0;
697+
double ft0mAmp = 0;
730698

731-
if (collision.has_foundFT0()) {
732-
auto ft0 = collision.foundFT0();
733-
for (const auto& amplitude : ft0.amplitudeA()) {
699+
if (foundBC.has_ft0()) {
700+
for (const auto& amplitude : foundBC.ft0().amplitudeA()) {
734701
ft0aAmp += amplitude;
735702
}
736-
for (const auto& amplitude : ft0.amplitudeC()) {
703+
for (const auto& amplitude : foundBC.ft0().amplitudeC()) {
737704
ft0cAmp += amplitude;
738705
}
739706
} else {
@@ -743,7 +710,7 @@ struct FlowGfwTask {
743710
registry.fill(HIST("FT0AAmp"), ft0aAmp);
744711
registry.fill(HIST("FT0CAmp"), ft0cAmp);
745712

746-
double ft0mAmp = ft0aAmp + ft0cAmp;
713+
ft0mAmp = ft0aAmp + ft0cAmp;
747714
registry.fill(HIST("FT0MAmp"), ft0mAmp);
748715

749716
// ZDC amplitude to use in fine binning
@@ -760,9 +727,12 @@ struct FlowGfwTask {
760727
registry.fill(HIST("ZPAcoll"), aZPA);
761728
registry.fill(HIST("ZPCcoll"), aZPC);
762729

730+
registry.fill(HIST("ZNvsFT0correl"), (ft0aAmp + ft0cAmp) / 100., aZNC + aZNA);
731+
763732
double aZDC = aZNC + aZNA + aZPA + aZPC;
764-
registry.fill(HIST("ZNvsFT0correl"), (ft0aAmp + ft0aAmp) / 100., aZNC + aZNA);
765733
registry.fill(HIST("ZDCAmp"), aZDC);
734+
registry.fill(HIST("ZNAmp"), aZNC + aZNA);
735+
registry.fill(HIST("ZPAmp"), aZPA + aZPC);
766736

767737
registry.fill(HIST("ZNvsZEMcoll"), aZEM1 + aZEM2, aZNA + aZNC);
768738

@@ -779,6 +749,30 @@ struct FlowGfwTask {
779749
}
780750
} // End of ZDC
781751

752+
float vtxz = collision.posZ();
753+
float lRandom = fRndm->Rndm();
754+
registry.fill(HIST("hVtxZ"), vtxz);
755+
registry.fill(HIST("hMult"), nTotal);
756+
registry.fill(HIST("hCent"), centrality);
757+
registry.fill(HIST("cent_vs_Nch"), centrality, nTotal);
758+
759+
fGFW->Clear();
760+
761+
auto bc = collision.bc_as<BCsRun3>();
762+
loadCorrections(bc.timestamp());
763+
registry.fill(HIST("hEventCount"), kCENTRALITY);
764+
765+
// fill event QA after cuts
766+
registry.fill(HIST("globalTracks_centT0C_Aft"), collision.centFT0C(), tracks.size());
767+
registry.fill(HIST("PVTracks_centT0C_Aft"), collision.centFT0C(), collision.multNTracksPV());
768+
registry.fill(HIST("globalTracks_PVTracks_Aft"), collision.multNTracksPV(), tracks.size());
769+
registry.fill(HIST("globalTracks_multT0A_Aft"), collision.multFT0A(), tracks.size());
770+
registry.fill(HIST("globalTracks_multV0A_Aft"), collision.multFV0A(), tracks.size());
771+
registry.fill(HIST("multV0A_multT0A_Aft"), collision.multFT0A(), collision.multFV0A());
772+
registry.fill(HIST("multT0C_centT0C_Aft"), collision.centFT0C(), collision.multFT0C());
773+
registry.fill(HIST("multT0A_centT0A_Aft"), collision.centFT0A(), collision.multFT0A());
774+
registry.fill(HIST("multFT0M_centFT0M_Aft"), collision.centFT0M(), collision.multFT0M());
775+
782776
// track weights
783777
float weff = 1, wacc = 1;
784778
int magnetfield = 0;
@@ -900,7 +894,7 @@ struct FlowGfwTask {
900894
}
901895

902896
} // End of process
903-
PROCESS_SWITCH(FlowGfwTask, processData, "Process analysis for Run 3 data", false);
897+
PROCESS_SWITCH(FlowGfwTask, processData, "Process analysis for Run 3 data", true);
904898

905899
// Filter the Reconstructed tracks
906900
Filter mytrackFilter = (nabs(aod::track::eta) < cfgCutEta) && (aod::track::pt > cfgCutPtMin) && (aod::track::pt < cfgCutPtMax) && (nabs(aod::track::dcaXY) < cfgCutDCAxy);

0 commit comments

Comments
 (0)