@@ -64,8 +64,6 @@ struct FlowZdcTask {
6464 Configurable<int > nBinsAmp{" nBinsAmp" , 1025 , " nbinsAmp" };
6565 Configurable<int > nBinsCent{" nBinsCent" , 90 , " nBinsCent" };
6666 Configurable<float > maxZn{" maxZn" , 125.5 , " Max ZN signal" };
67- Configurable<float > vtxRange{" vtxRange" , 10 .0f , " Vertex Z range to consider" };
68- Configurable<float > etaRange{" etaRange" , 1 .0f , " Eta range to consider" };
6967 // configs for process QA
7068 Configurable<int > nBinsNch{" nBinsNch" , 2501 , " N bins Nch (|eta|<0.8)" };
7169 Configurable<int > nBinsAmpFT0{" nBinsAmpFT0" , 100 , " N bins FT0 amp" };
@@ -74,8 +72,6 @@ struct FlowZdcTask {
7472 Configurable<int > nBinsAmpFV0{" nBinsAmpFV0" , 100 , " N bins FV0 amp" };
7573 Configurable<float > maxAmpFV0{" maxAmpFV0" , 2000 , " Max FV0 amp" };
7674 Configurable<int > nBinsZDC{" nBinsZDC" , 400 , " nBinsZDC" };
77- Configurable<int > nBinsZN{" nBinsZN" , 400 , " N bins ZN" };
78- Configurable<int > nBinsZP{" nBinsZP" , 160 , " N bins ZP" };
7975 Configurable<float > minNch{" minNch" , 0 , " Min Nch (|eta|<0.8)" };
8076 Configurable<float > maxNch{" maxNch" , 2500 , " Max Nch (|eta|<0.8)" };
8177 Configurable<int > nBinsTDC{" nBinsTDC" , 150 , " nbinsTDC" };
@@ -94,7 +90,6 @@ struct FlowZdcTask {
9490 Configurable<bool > isTDCcut{" isTDCcut" , false , " Use TDC cut?" };
9591 Configurable<bool > useMidRapNchSel{" useMidRapNchSel" , false , " Use mid-rapidity Nch selection" };
9692 Configurable<bool > applyEff{" applyEff" , true , " Apply track-by-track efficiency correction" };
97- Configurable<bool > applyFD{" applyFD" , false , " Apply track-by-track feed down correction" };
9893 Configurable<bool > correctNch{" correctNch" , true , " Correct also Nch" };
9994
10095 Configurable<float > nSigmaNchCut{" nSigmaNchCut" , 1 ., " nSigma Nch selection" };
@@ -107,16 +102,10 @@ struct FlowZdcTask {
107102 Configurable<float > maxPtSpectra{" maxPtSpectra" , 50 ., " maximum pt of the tracks" };
108103 Configurable<float > zemCut{" zemCut" , 100 ., " ZEM cut" };
109104 // axis configs
110- ConfigurableAxis axisPhi{" axisPhi" , {60 , 0.0 , constants::math::TwoPI}, " phi axis for histograms" };
111105 ConfigurableAxis axisMultiplicity{" axisMultiplicity" , {3500 , 0 , 3500 }, " centrality axis for histograms" };
112106 ConfigurableAxis axisZN{" axisZN" , {5000 , 0 , 500 }, " axisZN" };
113107 ConfigurableAxis axisZP{" axisZP" , {5000 , 0 , 500 }, " axisZP" };
114- ConfigurableAxis axisFT0CAmp{" axisFT0CAmp" , {5000 , 0 , 5000 }, " axisFT0CAmp" };
115- ConfigurableAxis axisFT0AAmp{" axisFT0AAmp" , {5000 , 0 , 5000 }, " axisFT0AAmp" };
116- ConfigurableAxis axisFT0MAmp{" axisFT0MAmp" , {10000 , 0 , 10000 }, " axisFT0MAmp" };
117- ConfigurableAxis multHistBin{" multHistBin" , {501 , -0.5 , 500.5 }, " " };
118108 ConfigurableAxis axisCent{" axisCent" , {10 , 0 , 100 }, " axisCent" };
119- ConfigurableAxis ft0cMultHistBin{" ft0cMultHistBin" , {501 , -0.5 , 500.5 }, " " };
120109 ConfigurableAxis binsPt{" binsPt" , {VARIABLE_WIDTH, 0.0 , 0.1 , 0.12 }, " pT binning" };
121110 Configurable<float > posZcut{" posZcut" , +10.0 , " z-vertex position cut" };
122111 Configurable<float > minEta{" minEta" , -0.8 , " minimum eta" };
@@ -176,21 +165,16 @@ struct FlowZdcTask {
176165 void init (InitContext const &)
177166 {
178167 // define axes
179- const AxisSpec axisCounter{1 , 0 , +1 , " " };
180168 const AxisSpec axisEvent{18 , 0.5 , 18.5 , " " };
181169 const AxisSpec axisZpos{48 , -12 ., 12 ., " Vtx_{z} (cm)" };
182170 const AxisSpec axisEta{40 , -1 ., +1 ., " #eta" };
183171 const AxisSpec axisPt{binsPt, " #it{p}_{T} (GeV/#it{c})" };
184172
185173 AxisSpec axisVtxZ{40 , -20 , 20 , " Vertex Z" , " VzAxis" };
186- AxisSpec axisMult = {multHistBin, " Mult" , " MultAxis" };
187- AxisSpec axisFT0CMult = {ft0cMultHistBin, " ft0c" , " FT0CMultAxis" };
188174
189175 // create histograms
190176 histos.add (" hEventCounter" , " Event counter" , kTH1F , {axisEvent});
191177 histos.add (" zPos" , " ;;Entries;" , kTH1F , {axisZpos});
192-
193- histos.add (" eventCounter" , " eventCounter" , kTH1F , {axisCounter});
194178 histos.add (" hZNvsFT0Ccent" ,
195179 " ZN Energy vs FT0C Centrality" ,
196180 kTH2F ,
@@ -204,17 +188,13 @@ struct FlowZdcTask {
204188 {{{nBinsNch, -0.5 , maxNch}, {nBinsNch, -0.5 , maxNch}}});
205189 histos.add (" T0Ccent" , " ;;Entries" , kTH1F , {axisCent});
206190 histos.add (" NchUncorrected" , " ;#it{N}_{ch} (|#eta| < 0.8);Entries;" , kTH1F , {{300 , 0 ., 3000 .}});
207- histos.add (" ZNamp" , " ;ZNA+ZNC;Entries;" , kTH1F , {{nBinsZN , -0.5 , maxZn}});
191+ histos.add (" ZNamp" , " ;ZNA+ZNC;Entries;" , kTH1F , {{nBinsZDC , -0.5 , maxZn}});
208192 histos.add (" ExcludedEvtVsFT0M" , " ;T0A+T0C (#times 1/100, -3.3 < #eta < -2.1 and 3.5 < #eta < 4.9);Entries;" , kTH1F , {{nBinsAmpFT0, 0 ., maxAmpFT0}});
209193 histos.add (" ExcludedEvtVsNch" , " ;#it{N}_{ch} (|#eta|<0.8);Entries;" , kTH1F , {{300 , 0 , 3000 }});
210194 histos.add (" Nch" , " ;#it{N}_{ch} (|#eta| < 0.8, Corrected);" , kTH1F , {{nBinsNch, minNch, maxNch}});
211- histos.add (" NchVsOneParCorr" , " ;#it{N}_{ch} (|#eta| < 0.8, Corrected);#LT[#it{p}_{T}^{(1)}]#GT (GeV/#it{c})" , kTProfile , {{nBinsNch, minNch, maxNch}});
212195 histos.add (" EtaVsPhi" , " ;#eta;#varphi" , kTH2F , {{{axisEta}, {100 , -0.1 * PI, +2.1 * PI}}});
213196 histos.add (" ZposVsEta" , " " , kTProfile , {axisZpos});
214- histos.add (" sigma1Pt" , " ;;#sigma(p_{T})/p_{T};" , kTProfile , {axisPt});
215197 histos.add (" dcaXYvspT" , " ;DCA_{xy} (cm);;" , kTH2F , {{{50 , -1 ., 1 .}, {axisPt}}});
216- histos.add (" GlobalMult_vs_FT0C" , " GlobalMult_vs_FT0C" , kTH2F , {axisMult, axisFT0CMult});
217- histos.add (" VtxZHist" , " VtxZHist" , kTH1D , {axisVtxZ});
218198
219199 // event selection steps
220200 histos.add (" eventSelectionSteps" , " eventSelectionSteps" , kTH1D , {axisEvent});
@@ -254,16 +234,31 @@ struct FlowZdcTask {
254234 histos.add (" ZPAenergy" , " common sum ZP energy side a" , kTH1F , {axisZP});
255235 histos.add (" ZNenergy" , " common sum zn (a + c sides) energy" , kTH1F , {axisZN});
256236 histos.add (" ZPenergy" , " common sum zp energy (a + c sides)" , kTH1F , {axisZP});
257- histos.add (" hZNvsFT0CAmp" , " ZN Energy vs FT0C Amplitude" , kTH2F , {axisFT0CAmp , axisZN});
258- histos.add (" hZPvsFT0CAmp" , " ZP Energy vs FT0C Amplitude" , kTH2F , {axisFT0CAmp , axisZP});
237+ histos.add (" hZNvsFT0CAmp" , " ZN Energy vs FT0C Amplitude" , kTH2F , {{nBinsAmpFT0, 0 ., maxAmpFT0} , axisZN});
238+ histos.add (" hZPvsFT0CAmp" , " ZP Energy vs FT0C Amplitude" , kTH2F , {{nBinsAmpFT0, 0 ., maxAmpFT0} , axisZP});
259239 histos.add (" hZNvsMult" , " ZN Energy vs Multiplicity" , kTH2F , {axisMultiplicity, axisZN});
260240 histos.add (" hZPvsMult" , " ZP Energy vs Multiplicity" , kTH2F , {axisMultiplicity, axisZP});
261241 }
262242
263243 if (doprocessQA) {
264- histos.add (" ZNVsFT0A" , " ;T0A (#times 1/100);ZNA+ZNC;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZn}}});
265- histos.add (" ZNVsFT0C" , " ;T0C (#times 1/100);ZNA+ZNC;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZn}}});
266- histos.add (" ZNVsFT0M" , " ;T0A+T0C (#times 1/100);ZNA+ZNC;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0M}, {nBinsZDC, -0.5 , maxZn}}});
244+ histos.add (" ZNVsFT0A" , " ;T0A (#times 1/100);ZNA+ZNC Sector Energy;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZn}}});
245+ histos.add (" ZNVsFT0C" , " ;T0C (#times 1/100);ZNA+ZNC Sector Energy;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZn}}});
246+ histos.add (" ZNVsFT0M" , " ;T0A+T0C (#times 1/100);ZNA+ZNC Sector Energy;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0M}, {nBinsZDC, -0.5 , maxZn}}});
247+ histos.add (" ZPVsFT0A" , " ;T0A (#times 1/100);ZPA+ZPC Sector Energy;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZp}}});
248+ histos.add (" ZPVsFT0C" , " ;T0C (#times 1/100);ZPA+ZPC Sector Energy;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZp}}});
249+ histos.add (" ZPVsFT0M" , " ;T0A+T0C (#times 1/100);ZPA+ZPC Sector Energy;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0M}, {nBinsZDC, -0.5 , maxZp}}});
250+ histos.add (" ZNAVsFT0A" , " ;T0A (#times 1/100);ZNA Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZn}}});
251+ histos.add (" ZNAVsFT0C" , " ;T0C (#times 1/100);ZNA Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZn}}});
252+ histos.add (" ZNAVsFT0M" , " ;T0A+T0C (#times 1/100);ZNA Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0M}, {nBinsZDC, -0.5 , maxZn}}});
253+ histos.add (" ZNCVsFT0A" , " ;T0A (#times 1/100);ZNC Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZn}}});
254+ histos.add (" ZNCVsFT0C" , " ;T0C (#times 1/100);ZNC Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZn}}});
255+ histos.add (" ZNCVsFT0M" , " ;T0A+T0C (#times 1/100);ZNC Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0M}, {nBinsZDC, -0.5 , maxZn}}});
256+ histos.add (" ZPAVsFT0A" , " ;T0A (#times 1/100);ZPA Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZp}}});
257+ histos.add (" ZPAVsFT0C" , " ;T0C (#times 1/100);ZPA Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZp}}});
258+ histos.add (" ZPAVsFT0M" , " ;T0A+T0C (#times 1/100);ZPA Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0M}, {nBinsZDC, -0.5 , maxZp}}});
259+ histos.add (" ZPCVsFT0A" , " ;T0A (#times 1/100);ZPC Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZp}}});
260+ histos.add (" ZPCVsFT0C" , " ;T0C (#times 1/100);ZPC Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZDC, -0.5 , maxZp}}});
261+ histos.add (" ZPCVsFT0M" , " ;T0A+T0C (#times 1/100);ZPC Sector;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0M}, {nBinsZDC, -0.5 , maxZp}}});
267262 histos.add (" ZN" , " ;ZNA+ZNC;Entries;" , kTH1F , {{nBinsZDC, -0.5 , maxZn}});
268263 histos.add (" ZNA" , " ;ZNA Amplitude;Entries;" , kTH1F , {{nBinsZDC, -0.5 , maxZn}});
269264 histos.add (" ZPA" , " ;ZPA Amplitude;Entries;" , kTH1F , {{nBinsZDC, -0.5 , maxZp}});
@@ -293,7 +288,6 @@ struct FlowZdcTask {
293288 histos.add (" ZEM1Vstdc" , " ;t_{ZEM1};ZEM1;" , kTH2F , {{{30 , -15 ., 15 .}, {30 , -0.5 , 2000.5 }}});
294289 histos.add (" ZEM2Vstdc" , " ;t_{ZEM2};ZEM2;" , kTH2F , {{{30 , -15 ., 15 .}, {30 , -0.5 , 2000.5 }}});
295290 histos.add (" debunch" , " ;t_{ZDC}-t_{ZDA};t_{ZDC}+t_{ZDA}" , kTH2F , {{{nBinsTDC, minTdc, maxTdc}, {nBinsTDC, minTdc, maxTdc}}});
296-
297291 histos.add (" GlbTracks" , " Nch" , kTH1F , {{nBinsNch, minNch, maxNch}});
298292 histos.add (" ampFT0C" , " ;T0C (#times 1/100);" , kTH1F , {{nBinsAmpFT0, 0 ., maxAmpFT0}});
299293 histos.add (" ampFT0A" , " ;T0A (#times 1/100);" , kTH1F , {{nBinsAmpFT0, 0 ., maxAmpFT0}});
@@ -317,6 +311,12 @@ struct FlowZdcTask {
317311 histos.add (" pZPCvsFT0Ccent" , " ;FT0C centrality;ZPC Amplitude" , kTProfile , {{nBinsCent, minT0CcentCut, maxT0CcentCut}});
318312 histos.add (" pZPAvsGlbTrack" , " ;Global Tracks (ITS + TPC);ZPA Amplitude" , kTProfile , {{nBinsNch, minNch, maxNch}});
319313 histos.add (" pZPCvsGlbTrack" , " ;Global Tracks (ITS + TPC);ZPC Amplitude" , kTProfile , {{nBinsNch, minNch, maxNch}});
314+ histos.add (" hZPASectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);ZPA Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZp}}});
315+ histos.add (" hZPCSectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);ZPC Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZp}}});
316+ histos.add (" hZNASectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);ZNA Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZn}}});
317+ histos.add (" hZNCSectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);ZNC Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZn}}});
318+ histos.add (" hZPSectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);(ZPA + ZPC) Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZp}}});
319+ histos.add (" hZNSectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);(ZNA + ZNC) Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZn}}});
320320 }
321321 if (doprocessZdc) {
322322 histos.add (" commonZna" , " ;ZNA Common Energy;Entries;" , kTH1F , {{nBinsZDC, -0.5 , maxZn}});
@@ -481,6 +481,8 @@ struct FlowZdcTask {
481481 float sumZNA = ((zdc.energySectorZNA ())[0 ] + (zdc.energySectorZNA ())[1 ] + (zdc.energySectorZNA ())[2 ] + (zdc.energySectorZNA ())[3 ]) / cfgCollisionEnergy;
482482 float sumZPC = ((zdc.energySectorZPC ())[0 ] + (zdc.energySectorZPC ())[1 ] + (zdc.energySectorZPC ())[2 ] + (zdc.energySectorZPC ())[3 ]) / cfgCollisionEnergy;
483483 float sumZPA = ((zdc.energySectorZPA ())[0 ] + (zdc.energySectorZPA ())[1 ] + (zdc.energySectorZPA ())[2 ] + (zdc.energySectorZPA ())[3 ]) / cfgCollisionEnergy;
484+ float sumSectZN = (sumZNC + sumZNA);
485+ float sumSectZP = (sumZPC + sumZPA);
484486
485487 // TDC cut
486488 if (isTDCcut) {
@@ -572,9 +574,24 @@ struct FlowZdcTask {
572574 histos.fill (HIST (" ZNCVstdc" ), tZNC, znC);
573575 histos.fill (HIST (" ZNAVstdc" ), tZNA, znA);
574576 histos.fill (HIST (" ZPCVstdc" ), tZPC, zpC);
575- histos.fill (HIST (" ZNVsFT0A" ), aT0A / 100 ., sumZNs);
576- histos.fill (HIST (" ZNVsFT0C" ), aT0C / 100 ., sumZNs);
577- histos.fill (HIST (" ZNVsFT0M" ), (aT0A + aT0C) / 100 ., sumZNs);
577+ histos.fill (HIST (" ZNVsFT0A" ), aT0A / 100 ., sumSectZN);
578+ histos.fill (HIST (" ZNVsFT0C" ), aT0C / 100 ., sumSectZN);
579+ histos.fill (HIST (" ZNVsFT0M" ), (aT0A + aT0C) / 100 ., sumSectZN);
580+ histos.fill (HIST (" ZPVsFT0A" ), aT0A / 100 ., sumSectZP);
581+ histos.fill (HIST (" ZPVsFT0C" ), aT0C / 100 ., sumSectZP);
582+ histos.fill (HIST (" ZPVsFT0M" ), (aT0A + aT0C) / 100 ., sumSectZP);
583+ histos.fill (HIST (" ZPAVsFT0A" ), aT0A / 100 ., sumZPA);
584+ histos.fill (HIST (" ZPAVsFT0C" ), aT0C / 100 ., sumZPA);
585+ histos.fill (HIST (" ZPAVsFT0M" ), (aT0A + aT0C) / 100 ., sumZPA);
586+ histos.fill (HIST (" ZPCVsFT0A" ), aT0A / 100 ., sumZPC);
587+ histos.fill (HIST (" ZPCVsFT0C" ), aT0C / 100 ., sumZPC);
588+ histos.fill (HIST (" ZPCVsFT0M" ), (aT0A + aT0C) / 100 ., sumZPC);
589+ histos.fill (HIST (" ZNCVsFT0A" ), aT0A / 100 ., sumZNC);
590+ histos.fill (HIST (" ZNCVsFT0C" ), aT0C / 100 ., sumZNC);
591+ histos.fill (HIST (" ZNCVsFT0M" ), (aT0A + aT0C) / 100 ., sumZNC);
592+ histos.fill (HIST (" ZNAVsFT0A" ), aT0A / 100 ., sumZNA);
593+ histos.fill (HIST (" ZNAVsFT0C" ), aT0C / 100 ., sumZNA);
594+ histos.fill (HIST (" ZNAVsFT0M" ), (aT0A + aT0C) / 100 ., sumZNA);
578595 histos.fill (HIST (" ZPA" ), zpA);
579596 histos.fill (HIST (" ZPACommon" ), commonSumZpa);
580597 histos.fill (HIST (" ZPASector" ), sumZPA);
@@ -588,11 +605,17 @@ struct FlowZdcTask {
588605 if (std::isfinite (zpA) && !std::isnan (zpA) && cent >= minT0CcentCut && cent < maxT0CcentCut && glbTracks >= minNch && glbTracks < maxNch) {
589606 histos.fill (HIST (" pZPAvsFT0Ccent" ), cent, zpA);
590607 histos.fill (HIST (" pZPAvsGlbTrack" ), glbTracks, zpA);
608+ histos.fill (HIST (" hZPASectorvsGlbTrack" ), glbTracks, sumZPA);
591609 }
592610 if (std::isfinite (zpC) && !std::isnan (zpC) && cent >= minT0CcentCut && cent < maxT0CcentCut && glbTracks >= minNch && glbTracks < maxNch) {
593611 histos.fill (HIST (" pZPCvsFT0Ccent" ), cent, zpC);
594612 histos.fill (HIST (" pZPCvsGlbTrack" ), glbTracks, zpC);
613+ histos.fill (HIST (" hZPCSectorvsGlbTrack" ), glbTracks, sumZPC);
595614 }
615+ histos.fill (HIST (" hZNASectorvsGlbTrack" ), glbTracks, sumZNA);
616+ histos.fill (HIST (" hZNCSectorvsGlbTrack" ), glbTracks, sumZNC);
617+ histos.fill (HIST (" hZPSectorvsGlbTrack" ), glbTracks, sumSectZP);
618+ histos.fill (HIST (" hZNSectorvsGlbTrack" ), glbTracks, sumSectZN);
596619 // ZDC Correlations
597620 histos.fill (HIST (" ZNAVsZNC" ), znC, znA);
598621 histos.fill (HIST (" ZNAVsZPA" ), zpA, znA);
@@ -619,6 +642,7 @@ struct FlowZdcTask {
619642 histos.fill (HIST (" NchVsFT0A" ), aT0A / 100 ., glbTracks);
620643 histos.fill (HIST (" NchVsFT0C" ), aT0C / 100 ., glbTracks);
621644 histos.fill (HIST (" NchVsFT0M" ), (aT0A + aT0C) / 100 ., glbTracks);
645+ histos.fill (HIST (" hNchvsNPV" ), collision.multNTracksPVeta1 (), tracks.size ());
622646 histos.fill (HIST (" NchVsEt" ), et, glbTracks);
623647 histos.fill (HIST (" NchVsITStracks" ), itsTracks, glbTracks);
624648 if (glbTracks >= minNchSel) {
@@ -786,7 +810,6 @@ struct FlowZdcTask {
786810 histos.fill (HIST (" hZPvsFT0CAmp" ), ft0cAmp, sumZP);
787811 histos.fill (HIST (" hZNvsMult" ), nTot, sumZN);
788812 histos.fill (HIST (" hZPvsMult" ), nTot, sumZP);
789- histos.fill (HIST (" hNchvsNPV" ), collision.multNTracksPVeta1 (), nTot);
790813 histos.fill (HIST (" Nch" ), nchMult);
791814 histos.fill (HIST (" ZNamp" ), sumZNs);
792815 histos.fill (HIST (" NchVsZN" ), nchMult, sumZNs);
0 commit comments