@@ -43,6 +43,9 @@ struct flowTest {
4343 Configurable<float > maxB{" maxB" , 20 .0f , " max impact parameter" };
4444 Configurable<int > pdgSelection{" pdgSelection" , 0 , " pdg code selection for tracking study (0: no selection)" };
4545
46+ Configurable<int > analysisMinimumITSClusters{" analysisMinimumITSClusters" , 5 , " minimum ITS clusters for analysis track category" };
47+ Configurable<int > analysisMinimumTPCClusters{" analysisMinimumTPCClusters" , 70 , " minimum TPC clusters for analysis track category" };
48+
4649 ConfigurableAxis axisB{" axisB" , {100 , 0 .0f , 20 .0f }, " " };
4750 ConfigurableAxis axisPhi{" axisPhi" , {100 , 0 .0f , 2 .0f * TMath::Pi ()}, " " };
4851 ConfigurableAxis axisNch{" axisNch" , {300 , 0 .0f , 3000 .0f }, " Nch in |eta|<0.8" };
@@ -60,6 +63,8 @@ struct flowTest {
6063 histos.add <TH3>(" hBVsPtVsPhiGenerated" , " hBVsPtVsPhiGenerated" , HistType::kTH3D , {axisB, axisPhi, axisPt});
6164 histos.add <TH3>(" hBVsPtVsPhiGlobal" , " hBVsPtVsPhiGlobal" , HistType::kTH3D , {axisB, axisPhi, axisPt});
6265 histos.add <TH3>(" hBVsPtVsPhiGlobalFake" , " hBVsPtVsPhiGlobalFake" , HistType::kTH3D , {axisB, axisPhi, axisPt});
66+ histos.add <TH3>(" hBVsPtVsPhiAnalysis" , " hBVsPtVsPhiAnalysis" , HistType::kTH3D , {axisB, axisPhi, axisPt});
67+ histos.add <TH3>(" hBVsPtVsPhiAnalysisFake" , " hBVsPtVsPhiAnalysisFake" , HistType::kTH3D , {axisB, axisPhi, axisPt});
6368 histos.add <TH3>(" hBVsPtVsPhiAny" , " hBVsPtVsPhiAny" , HistType::kTH3D , {axisB, axisPhi, axisPt});
6469 histos.add <TH3>(" hBVsPtVsPhiTPCTrack" , " hBVsPtVsPhiTPCTrack" , HistType::kTH3D , {axisB, axisPhi, axisPt});
6570 histos.add <TH3>(" hBVsPtVsPhiITSTrack" , " hBVsPtVsPhiITSTrack" , HistType::kTH3D , {axisB, axisPhi, axisPt});
@@ -126,6 +131,8 @@ struct flowTest {
126131 bool validITSTrackFake = false ;
127132 bool validITSABTrack = false ;
128133 bool validITSABTrackFake = false ;
134+ bool validAnalysisTrack = false ;
135+ bool validAnalysisTrackFake = false ;
129136 if (mcParticle.has_tracks ()) {
130137 auto const & tracks = mcParticle.tracks_as <recoTracks>();
131138 for (auto const & track : tracks) {
@@ -137,6 +144,12 @@ struct flowTest {
137144 }
138145 }
139146
147+ if (track.tpcNClsFound () >= analysisMinimumTPCClusters && track.itsNCls () >= analysisMinimumITSClusters) {
148+ validAnalysisTrack = true ;
149+ if (isITSFake) {
150+ validAnalysisTrackFake = true ;
151+ }
152+ }
140153 if (track.hasTPC () && track.hasITS ()) {
141154 validGlobal = true ;
142155 if (isITSFake) {
@@ -172,6 +185,12 @@ struct flowTest {
172185 if (validGlobalFake) {
173186 histos.fill (HIST (" hBVsPtVsPhiGlobalFake" ), imp, deltaPhi, mcParticle.pt ());
174187 }
188+ if (validAnalysisTrack) {
189+ histos.fill (HIST (" hBVsPtVsPhiAnalysis" ), imp, deltaPhi, mcParticle.pt ());
190+ }
191+ if (validAnalysisTrackFake) {
192+ histos.fill (HIST (" hBVsPtVsPhiAnalysisFake" ), imp, deltaPhi, mcParticle.pt ());
193+ }
175194 // if any track present, fill
176195 if (validTrack)
177196 histos.fill (HIST (" hBVsPtVsPhiAny" ), imp, deltaPhi, mcParticle.pt ());
0 commit comments