@@ -82,6 +82,7 @@ struct HfTaskSingleElectron {
8282 Configurable<float > tpcNClsFoundOverFindableMin{" tpcNClsFoundOverFindableMin" , 0.8 , " min # of TPC found/findable clusters" };
8383 Configurable<float > tpcChi2perNClMax{" tpcChi2perNClMax" , 4 ., " min # of tpc chi2 per clusters" };
8484 Configurable<int > itsIBClsMin{" itsIBClsMin" , 3 , " min # of its clusters in IB" };
85+ Configurable<float > itsChi2perNClMax{" itsChi2perNClMax" , 6 ., " min # of tpc chi2 per clusters" };
8586 Configurable<float > dcaxyMax{" dcaxyMax" , 1 ., " max of track dca in xy" };
8687 Configurable<float > dcazMax{" dcazMax" , 2 ., " max of track dca in z" };
8788 Configurable<float > tofNSigmaMax{" tofNSigmaMax" , 3 ., " max of tof nsigma" };
@@ -95,7 +96,7 @@ struct HfTaskSingleElectron {
9596
9697 // using declarations
9798 using MyCollisions = soa::Join<aod::Collisions, aod::EvSels>;
98- using TracksEl = soa::Join<aod::Tracks, aod::TrackExtra , aod::TracksDCA, aod::pidTOFFullEl, aod::pidTPCFullEl>;
99+ using TracksEl = soa::Join<aod::Tracks, aod::TrackSelection, aod::TrackSelectionExtension, aod::TracksExtra , aod::TracksDCA, aod::pidTOFFullEl, aod::pidTPCFullEl>;
99100 using McTracksEl = soa::Join<aod::Tracks, aod::TrackExtra, aod::TracksDCA, aod::pidTOFFullEl, aod::pidTPCFullEl, aod::McTrackLabels>;
100101
101102 // Filter
@@ -131,6 +132,7 @@ struct HfTaskSingleElectron {
131132 histos.add (" tpcFoundFindableTrack" , " " , kTH1D , {{10 , 0 , 1 }});
132133 histos.add (" tpcChi2Track" , " " , kTH1D , {{100 , 0 , 10 }});
133134 histos.add (" itsIBClsTrack" , " " , kTH1D , {{10 , 0 , 10 }});
135+ histos.add (" itsChi2Track" , " " , kTH1D , {{50 , 0 , 50 }});
134136 histos.add (" dcaXYTrack" , " " , kTH1D , {{600 , -3 , 3 }});
135137 histos.add (" dcaZTrack" , " " , kTH1D , {{600 , -3 , 3 }});
136138
@@ -168,9 +170,11 @@ struct HfTaskSingleElectron {
168170 if (track.tpcNClsCrossedRows () < tpcNCrossedRowMin) {
169171 return false ;
170172 }
173+
171174 if (track.tpcCrossedRowsOverFindableCls () < tpcNClsFoundOverFindableMin) {
172175 return false ;
173176 }
177+
174178 if (track.tpcChi2NCl () > tpcChi2perNClMax) {
175179 return false ;
176180 }
@@ -179,9 +183,14 @@ struct HfTaskSingleElectron {
179183 return false ;
180184 }
181185
186+ if (track.itsChi2NCl () > itsChi2perNClMax) {
187+ return false ;
188+ }
189+
182190 if (std::abs (track.dcaXY ()) > dcaxyMax) {
183191 return false ;
184192 }
193+
185194 if (std::abs (track.dcaZ ()) > dcazMax) {
186195 return false ;
187196 }
@@ -405,13 +414,18 @@ struct HfTaskSingleElectron {
405414 continue ;
406415 }
407416
417+ if (!(track.passedITSRefit () && track.passedTPCRefit ())) {
418+ continue ;
419+ }
420+
408421 histos.fill (HIST (" etaTrack" ), track.eta ());
409422 histos.fill (HIST (" ptTrack" ), track.pt ());
410423
411424 histos.fill (HIST (" tpcNClsTrack" ), track.tpcNClsCrossedRows ());
412425 histos.fill (HIST (" tpcFoundFindableTrack" ), track.tpcCrossedRowsOverFindableCls ());
413426 histos.fill (HIST (" tpcChi2Track" ), track.tpcChi2NCl ());
414427 histos.fill (HIST (" itsIBClsTrack" ), track.itsNClsInnerBarrel ());
428+ histos.fill (HIST (" itsChi2Track" ), track.itsChi2NCl ());
415429 histos.fill (HIST (" dcaXYTrack" ), track.dcaXY ());
416430 histos.fill (HIST (" dcaZTrack" ), track.dcaZ ());
417431
0 commit comments