@@ -195,11 +195,11 @@ struct UccZdc {
195195
196196 registry.add (" zPos" , " ;;Entries;" , kTH1F , {axisZpos});
197197 registry.add (" T0Ccent" , " ;;Entries" , kTH1F , {axisCent});
198- registry.add (" NchUncorrected" , " ;#it{N}_{ch} (|#eta| < 0.8);Entries;" , kTH1F , {{300 , 0 ., 3000 . }});
198+ registry.add (" NchUncorrected" , " ;#it{N}_{ch} (|#eta| < 0.8);Entries;" , kTH1F , {{nBinsNch, minNch, maxNch }});
199199 registry.add (" hEventCounter" , " ;;Events" , kTH1F , {axisEvent});
200200 registry.add (" ZNamp" , " ;ZNA+ZNC;Entries;" , kTH1F , {{nBinsZN, -0.5 , maxZN}});
201201 registry.add (" ExcludedEvtVsFT0M" , " ;T0A+T0C (#times 1/100, -3.3 < #eta < -2.1 and 3.5 < #eta < 4.9);Entries;" , kTH1F , {{nBinsAmpFT0, 0 ., maxAmpFT0}});
202- registry.add (" ExcludedEvtVsNch" , " ;#it{N}_{ch} (|#eta|<0.8);Entries;" , kTH1F , {{300 , 0 , 3000 }});
202+ registry.add (" ExcludedEvtVsNch" , " ;#it{N}_{ch} (|#eta|<0.8);Entries;" , kTH1F , {{nBinsNch, minNch, maxNch }});
203203 registry.add (" Nch" , " ;#it{N}_{ch} (|#eta| < 0.8, Corrected);" , kTH1F , {{nBinsNch, minNch, maxNch}});
204204 registry.add (" NchVsOneParCorr" , " ;#it{N}_{ch} (|#eta| < 0.8, Corrected);#LT[#it{p}_{T}^{(1)}]#GT (GeV/#it{c})" , kTProfile , {{nBinsNch, minNch, maxNch}});
205205 registry.add (" EtaVsPhi" , " ;#eta;#varphi" , kTH2F , {{{axisEta}, {100 , -0.1 * PI, +2.1 * PI}}});
@@ -299,24 +299,24 @@ struct UccZdc {
299299 if (doprocessQA) {
300300 registry.add (" Debunch" , " ;t_{ZDC}-t_{ZDA};t_{ZDC}+t_{ZDA}" , kTH2F , {{{nBinsTDC, minTdc, maxTdc}, {nBinsTDC, minTdc, maxTdc}}});
301301 registry.add (" NchVsFT0M" , " ;T0A+T0C (#times 1/100, -3.3 < #eta < -2.1 and 3.5 < #eta < 4.9);#it{N}_{ch} (|#eta|<0.8);" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsNch, minNch, maxNch}}});
302- registry.add (" NchVsFT0A" , " ;T0A (#times 1/100, 3.5 < #eta < 4.9);#it{N}_{ch} (|#eta|<0.8);" , kTH2F , {{{80 , 0 ., 1800 . }, {nBinsNch, minNch, maxNch}}});
303- registry.add (" NchVsFT0C" , " ;T0C (#times 1/100, -3.3 < #eta < -2.1);#it{N}_{ch} (|#eta|<0.8);" , kTH2F , {{{80 , 0 ., 600 . }, {nBinsNch, minNch, maxNch}}});
302+ registry.add (" NchVsFT0A" , " ;T0A (#times 1/100, 3.5 < #eta < 4.9);#it{N}_{ch} (|#eta|<0.8);" , kTH2F , {{{nBinsAmpFT0 , 0 ., maxAmpFT0 }, {nBinsNch, minNch, maxNch}}});
303+ registry.add (" NchVsFT0C" , " ;T0C (#times 1/100, -3.3 < #eta < -2.1);#it{N}_{ch} (|#eta|<0.8);" , kTH2F , {{{nBinsAmpFT0 , 0 ., maxAmpFT0 }, {nBinsNch, minNch, maxNch}}});
304304 registry.add (" NchVsFV0A" , " ;V0A (#times 1/100, 2.2 < #eta < 5);#it{N}_{ch} (|#eta|<0.8);" , kTH2F , {{{80 , 0 ., maxAmpFV0}, {nBinsNch, minNch, maxNch}}});
305305 registry.add (" NchVsEt" , " ;#it{E}_{T} (|#eta|<0.8);#LTITS+TPC tracks#GT (|#eta|<0.8);" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsNch, minNch, maxNch}}});
306- registry.add (" NchVsNPV" , " ;#it{N}_{PV} (|#eta|<1);ITS+TPC tracks (|#eta|<0.8);" , kTH2F , {{{300 , - 0.5 , 5999.5 }, {nBinsNch, minNch, maxNch}}});
307- registry.add (" NchVsITStracks" , " ;ITS tracks nCls >= 5;TITS+TPC tracks (|#eta|<0.8);" , kTH2F , {{{300 , - 0.5 , 5999.5 }, {nBinsNch, minNch, maxNch}}});
308- registry.add (" ZNVsFT0A" , " ;T0A (#times 1/100);ZNA+ZNC;" , kTH2F , {{{80 , 0 ., 1800 . }, {nBinsZN, -0.5 , maxZN}}});
309- registry.add (" ZNVsFT0C" , " ;T0C (#times 1/100);ZNA+ZNC;" , kTH2F , {{{80 , 0 ., 600 . }, {nBinsZN, -0.5 , maxZN}}});
306+ registry.add (" NchVsNPV" , " ;#it{N}_{PV} (|#eta|<1);ITS+TPC tracks (|#eta|<0.8);" , kTH2F , {{{nBinsITSTrack, minITSTrack, maxITSTrack }, {nBinsNch, minNch, maxNch}}});
307+ registry.add (" NchVsITStracks" , " ;ITS tracks nCls >= 5;TITS+TPC tracks (|#eta|<0.8);" , kTH2F , {{{nBinsITSTrack, minITSTrack, maxITSTrack }, {nBinsNch, minNch, maxNch}}});
308+ registry.add (" ZNVsFT0A" , " ;T0A (#times 1/100);ZNA+ZNC;" , kTH2F , {{{nBinsAmpFT0 , 0 ., maxAmpFT0 }, {nBinsZN, -0.5 , maxZN}}});
309+ registry.add (" ZNVsFT0C" , " ;T0C (#times 1/100);ZNA+ZNC;" , kTH2F , {{{nBinsAmpFT0 , 0 ., maxAmpFT0 }, {nBinsZN, -0.5 , maxZN}}});
310310 registry.add (" ZNVsFT0M" , " ;T0A+T0C (#times 1/100);ZNA+ZNC;" , kTH2F , {{{nBinsAmpFT0, 0 ., maxAmpFT0}, {nBinsZN, -0.5 , maxZN}}});
311311 registry.add (" ZNAamp" , " ;ZNA;Entries;" , kTH1F , {{nBinsZN, -0.5 , maxZN}});
312312 registry.add (" ZPAamp" , " ;ZPA;Entries;" , kTH1F , {{nBinsZP, -0.5 , maxZP}});
313313 registry.add (" ZNCamp" , " ;ZNC;Entries;" , kTH1F , {{nBinsZN, -0.5 , maxZN}});
314314 registry.add (" ZPCamp" , " ;ZPC;Entries;" , kTH1F , {{nBinsZP, -0.5 , maxZP}});
315- registry.add (" ZNAVsZNC" , " ;ZNC;ZNA" , kTH2F , {{{30 , -0.5 , maxZN}, {30 , -0.5 , maxZN}}});
316- registry.add (" ZPAVsZPC" , " ;ZPC;ZPA;" , kTH2F , {{{100 , -0.5 , maxZP}, {100 , -0.5 , maxZP}}});
317- registry.add (" ZNAVsZPA" , " ;ZPA;ZNA;" , kTH2F , {{{20 , -0.5 , maxZP}, {30 , -0.5 , maxZN}}});
318- registry.add (" ZNCVsZPC" , " ;ZPC;ZNC;" , kTH2F , {{{20 , -0.5 , maxZP}, {30 , -0.5 , maxZN}}});
319- registry.add (" ZNVsZEM" , " ;ZEM;ZNA+ZNC;" , kTH2F , {{{60 , -0.5 , maxZEM}, {60 , -0.5 , maxZN}}});
315+ registry.add (" ZNAVsZNC" , " ;ZNC;ZNA" , kTH2F , {{{nBinsZN , -0.5 , maxZN}, {nBinsZN , -0.5 , maxZN}}});
316+ registry.add (" ZPAVsZPC" , " ;ZPC;ZPA;" , kTH2F , {{{nBinsZP , -0.5 , maxZP}, {nBinsZP , -0.5 , maxZP}}});
317+ registry.add (" ZNAVsZPA" , " ;ZPA;ZNA;" , kTH2F , {{{nBinsZP , -0.5 , maxZP}, {nBinsZN , -0.5 , maxZN}}});
318+ registry.add (" ZNCVsZPC" , " ;ZPC;ZNC;" , kTH2F , {{{nBinsZP , -0.5 , maxZP}, {nBinsZN , -0.5 , maxZN}}});
319+ registry.add (" ZNVsZEM" , " ;ZEM;ZNA+ZNC;" , kTH2F , {{{nBinsZN , -0.5 , maxZEM}, {nBinsZN , -0.5 , maxZN}}});
320320 registry.add (" ZNCVsNch" , " ;#it{N}_{ch} (|#eta|<0.8);ZNC;" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZN, minZN, maxZN}}});
321321 registry.add (" ZNAVsNch" , " ;#it{N}_{ch} (|#eta|<0.8);ZNA;" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZN, minZN, maxZN}}});
322322 registry.add (" ZNVsNch" , " ;#it{N}_{ch} (|#eta|<0.8);ZNA+ZNC;" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZN, minZN, maxZN}}});
@@ -464,7 +464,6 @@ struct UccZdc {
464464 return ;
465465 }
466466 registry.fill (HIST (" hEventCounter" ), EvCutLabel::TZero);
467-
468467 if (foundBC.has_fv0a ()) {
469468 for (const auto & amplitude : foundBC.fv0a ().amplitude ()) {
470469 aV0A += amplitude;
@@ -492,7 +491,6 @@ struct UccZdc {
492491 zpC /= kCollEnergy ;
493492 float sumZNs{znA + znC};
494493 float sumZEMs{aZEM1 + aZEM2};
495-
496494 // TDC cut
497495 if (isTDCcut) {
498496 if (std::sqrt (std::pow (tZDCdif, 2 .) + std::pow (tZDCsum, 2 .)) > tdcCut) {
@@ -543,15 +541,14 @@ struct UccZdc {
543541 const double nSigmaSelection{nSigmaNchCut * sigmaNch};
544542 const double diffMeanNch{meanNch - glbTracks};
545543
546- if (!( std::abs (diffMeanNch) < nSigmaSelection) ) {
544+ if (std::abs (diffMeanNch) > nSigmaSelection) {
547545 registry.fill (HIST (" ExcludedEvtVsFT0M" ), normT0M);
548546 registry.fill (HIST (" ExcludedEvtVsNch" ), glbTracks);
549- } else {
550547 skipEvent = true ;
551548 }
552549 }
553550
554- if (! skipEvent) {
551+ if (useMidRapNchSel && skipEvent) {
555552 return ;
556553 }
557554
@@ -602,9 +599,8 @@ struct UccZdc {
602599 registry.fill (HIST (" ZNCVsNch" ), glbTracks, znC);
603600 registry.fill (HIST (" ZNVsNch" ), glbTracks, sumZNs);
604601 registry.fill (HIST (" ZNDifVsNch" ), glbTracks, znA - znC);
605- if (glbTracks >= minNchSel) {
602+ if (glbTracks >= minNchSel)
606603 registry.fill (HIST (" NchVsOneParCorr" ), glbTracks, meanpt / glbTracks);
607- }
608604 }
609605 PROCESS_SWITCH (UccZdc, processQA, " Process QA" , true );
610606
0 commit comments