@@ -100,15 +100,16 @@ struct HfCandidateSelectorXicToPKPi {
100100 selectorKaon = selectorPion;
101101 selectorProton = selectorPion;
102102
103- if (activateQA) {
103+ if (activateQA) {
104104 constexpr int kNBinsSelections = 1 + aod::SelectionStep::NSelectionSteps;
105105 std::string labels[kNBinsSelections ];
106106 labels[0 ] = " No selection" ;
107107 labels[1 + aod::SelectionStep::RecoSkims] = " Skims selection" ;
108108 labels[1 + aod::SelectionStep::RecoTopol] = " Skims & Topological selections" ;
109109 labels[1 + aod::SelectionStep::RecoPID] = " Skims & Topological & PID selections" ;
110110 labels[1 + aod::SelectionStep::RecoMl] = " ML selection" ;
111- static const AxisSpec axisSelections = {kNBinsSelections , 0.5 , kNBinsSelections + 0.5 , " " }; registry.add (" hSelections" , " Selections;;#it{p}_{T} (GeV/#it{c})" , {HistType::kTH2F , {axisSelections, {(std::vector<double >)binsPt, " #it{p}_{T} (GeV/#it{c})" }}});
111+ static const AxisSpec axisSelections = {kNBinsSelections , 0.5 , kNBinsSelections + 0.5 , " " };
112+ registry.add (" hSelections" , " Selections;;#it{p}_{T} (GeV/#it{c})" , {HistType::kTH2F , {axisSelections, {(std::vector<double >)binsPt, " #it{p}_{T} (GeV/#it{c})" }}});
112113 for (int iBin = 0 ; iBin < kNBinsSelections ; ++iBin) {
113114 registry.get <TH2>(HIST (" hSelections" ))->GetXaxis ()->SetBinLabel (iBin + 1 , labels[iBin].data ());
114115 }
@@ -242,12 +243,12 @@ struct HfCandidateSelectorXicToPKPi {
242243 hfMlXicToPKPiCandidate (outputMlXicToPKPi, outputMlXicToPiKP);
243244 }
244245 if (activateQA) {
245- registry.fill (HIST (" hSelections" ), 1 , ptCand);
246+ registry.fill (HIST (" hSelections" ), 1 , ptCand);
246247 }
247248 continue ;
248249 }
249- if (activateQA) {
250- registry.fill (HIST (" hSelections" ), 2 + aod::SelectionStep::RecoSkims, ptCand);
250+ if (activateQA) {
251+ registry.fill (HIST (" hSelections" ), 2 + aod::SelectionStep::RecoSkims, ptCand);
251252 }
252253 SETBIT (statusXicToPKPi, aod::SelectionStep::RecoSkims);
253254 SETBIT (statusXicToPiKP, aod::SelectionStep::RecoSkims);
@@ -285,7 +286,7 @@ struct HfCandidateSelectorXicToPKPi {
285286 if (topolXicToPiKP) {
286287 SETBIT (statusXicToPiKP, aod::SelectionStep::RecoTopol);
287288 }
288-
289+
289290 if (activateQA) {
290291 registry.fill (HIST (" hSelections" ), 2 + aod::SelectionStep::RecoTopol, candidate.pt ());
291292 }
@@ -299,25 +300,25 @@ struct HfCandidateSelectorXicToPKPi {
299300 pidXicToPiKP = 1 ;
300301 } else {
301302 // track-level PID selection
302- TrackSelectorPID::Status pidTrackPos1Proton = TrackSelectorPID::Accepted;
303- TrackSelectorPID::Status pidTrackPos2Proton = TrackSelectorPID::Accepted;
304- TrackSelectorPID::Status pidTrackPos1Pion = TrackSelectorPID::Accepted;
305- TrackSelectorPID::Status pidTrackPos2Pion = TrackSelectorPID::Accepted;
306- TrackSelectorPID::Status pidTrackNegKaon = TrackSelectorPID::Accepted;
307- if (usePidTpcAndTof){
308-
309- pidTrackPos1Proton = selectorProton.statusTpcAndTof (trackPos1);
310- pidTrackPos2Proton = selectorProton.statusTpcAndTof (trackPos2);
311- pidTrackPos1Pion = selectorPion.statusTpcAndTof (trackPos1);
312- pidTrackPos2Pion = selectorPion.statusTpcAndTof (trackPos2);
313- pidTrackNegKaon = selectorKaon.statusTpcAndTof (trackNeg);
314- } else {
315- pidTrackPos1Proton = selectorProton.statusTpcOrTof (trackPos1);
316- pidTrackPos2Proton = selectorProton.statusTpcOrTof (trackPos2);
317- pidTrackPos1Pion = selectorPion.statusTpcOrTof (trackPos1);
318- pidTrackPos2Pion = selectorPion.statusTpcOrTof (trackPos2);
319- pidTrackNegKaon = selectorKaon.statusTpcOrTof (trackNeg);
320- }
303+ TrackSelectorPID::Status pidTrackPos1Proton = TrackSelectorPID::Accepted;
304+ TrackSelectorPID::Status pidTrackPos2Proton = TrackSelectorPID::Accepted;
305+ TrackSelectorPID::Status pidTrackPos1Pion = TrackSelectorPID::Accepted;
306+ TrackSelectorPID::Status pidTrackPos2Pion = TrackSelectorPID::Accepted;
307+ TrackSelectorPID::Status pidTrackNegKaon = TrackSelectorPID::Accepted;
308+ if (usePidTpcAndTof) {
309+
310+ pidTrackPos1Proton = selectorProton.statusTpcAndTof (trackPos1);
311+ pidTrackPos2Proton = selectorProton.statusTpcAndTof (trackPos2);
312+ pidTrackPos1Pion = selectorPion.statusTpcAndTof (trackPos1);
313+ pidTrackPos2Pion = selectorPion.statusTpcAndTof (trackPos2);
314+ pidTrackNegKaon = selectorKaon.statusTpcAndTof (trackNeg);
315+ } else {
316+ pidTrackPos1Proton = selectorProton.statusTpcOrTof (trackPos1);
317+ pidTrackPos2Proton = selectorProton.statusTpcOrTof (trackPos2);
318+ pidTrackPos1Pion = selectorPion.statusTpcOrTof (trackPos1);
319+ pidTrackPos2Pion = selectorPion.statusTpcOrTof (trackPos2);
320+ pidTrackNegKaon = selectorKaon.statusTpcOrTof (trackNeg);
321+ }
321322
322323 if (pidTrackPos1Proton == TrackSelectorPID::Accepted &&
323324 pidTrackNegKaon == TrackSelectorPID::Accepted &&
@@ -384,7 +385,7 @@ struct HfCandidateSelectorXicToPKPi {
384385 if (isSelectedMlXicToPiKP) {
385386 SETBIT (statusXicToPiKP, aod::SelectionStep::RecoMl);
386387 }
387- if (activateQA) {
388+ if (activateQA) {
388389 registry.fill (HIST (" hSelections" ), 2 + aod::SelectionStep::RecoMl, candidate.pt ());
389390 }
390391 }
0 commit comments