@@ -554,8 +554,8 @@ struct EbyeMaker {
554554
555555 // tracking variables QA
556556 histos.add <TH2>(" QA/tpcCRvsCls" , " ;#it{N}_{TPCCR};#it{N}_{TPCcls}" , HistType::kTH2F , {nTpcAxis, nTpcAxis});
557- histos.add <TH2>(" QA/dcaxyVsPt" , " ;#it{p}_{T} (GeV/#it{c});DCA_{#it{xy}} (cm)" , HistType::kTH2F , {momAxis , dcaAxis});
558- histos.add <TH2>(" QA/dcazVsPt" , " ;#it{p}_{T} (GeV/#it{c});DCA_{#it{z}} (cm)" , HistType::kTH2F , {momAxis , dcaAxis});
557+ histos.add <TH2>(" QA/dcaxyVsPt" , " ;#it{p}_{T} (GeV/#it{c});DCA_{#it{xy}} (cm)" , HistType::kTH2F , {signMomAxis , dcaAxis});
558+ histos.add <TH2>(" QA/dcazVsPt" , " ;#it{p}_{T} (GeV/#it{c});DCA_{#it{z}} (cm)" , HistType::kTH2F , {signMomAxis , dcaAxis});
559559 histos.add <TH3>(" QA/phiVsPtVsCls" , " ;#it{N}^{TPC}_{CR};#it{p}_{T} (GeV/#it{c});#phi (rad)" , HistType::kTH3F , {nTpcAxis, signMomAxis, phiAxis});
560560
561561 ptMin = std::array<float , kNpart >{antipPtMin, antidPtMin};
@@ -604,12 +604,12 @@ struct EbyeMaker {
604604 auto trackPt = trackParCov.getPt ();
605605 auto trackPhi = trackParCov.getPhi ();
606606 auto trackEta = trackParCov.getEta ();
607- histos.fill (HIST (" QA/dcaxyVsPt" ), track.pt (), dcaInfo[0 ]); // TODO: this should rather be trackPt (likely small effect for 7 ITS clusters?)
608- histos.fill (HIST (" QA/dcazVsPt" ), track.pt (), dcaInfo[1 ]);
609- if (std::abs (dcaInfo[0 ]) > dcaSigma (track.pt (), " dcaxy" ) || std::abs (dcaInfo[1 ]) > dcaSigma (track.pt (), " dcaz" ) || dca > dcaSigma (track.pt (), " dca" ))
610- continue ;
611607 if (track.tpcNClsFound () < trackNclusTpcCut || track.tpcNClsCrossedRows () < trackNcrossedRows)
612608 continue ;
609+ histos.fill (HIST (" QA/dcaxyVsPt" ), track.sign () > 0 . ? track.pt () : -track.pt (), dcaInfo[0 ]); // TODO: this should rather be trackPt (likely small effect for 7 ITS clusters?)
610+ histos.fill (HIST (" QA/dcazVsPt" ), track.sign () > 0 . ? track.pt () : -track.pt (), dcaInfo[1 ]);
611+ if (std::abs (dcaInfo[0 ]) > dcaSigma (track.pt (), " dcaxy" ) || std::abs (dcaInfo[1 ]) > dcaSigma (track.pt (), " dcaz" ) || dca > dcaSigma (track.pt (), " dca" ))
612+ continue ;
613613 histos.fill (HIST (" QA/phiVsPtVsCls" ), track.tpcNClsFound (), track.sign () > 0 . ? trackPt : -trackPt, trackPhi);
614614 histos.fill (HIST (" QA/tpcSignal" ), track.tpcInnerParam (), track.tpcSignal ());
615615 if (trackPt > ptMin[0 ] && trackPt < ptMax[0 ] && ((track.sign () < 0 && countOnlyLSTrk == TracksCharge::kNegative ) || (track.sign () > 0 && countOnlyLSTrk == TracksCharge::kPositive ) || (countOnlyLSTrk == TracksCharge::kAll ))) {
0 commit comments