add cluster cut

f3sch · Nov 20, 2023 · 7b9907b · 7b9907b
1 parent 8de0701
commit 7b9907b
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 0 deletions.
diff --git a/Detectors/Vertexing/include/DetectorsVertexing/SVertexerParams.h b/Detectors/Vertexing/include/DetectorsVertexing/SVertexerParams.h
@@ -99,6 +99,7 @@ struct SVertexerParams : public o2::conf::ConfigurableParamHelper<SVertexerParam
   float mTPCTrackMaxX = constants::geom::XTPCInnerRef; // don't use TPC standalone tracks with X exceeding this;
   float mTPCTrack2Beam = 21.f;                         // straight line for TPC track back to beamline
   bool mTPCTrackPhotonTune = true;                     // use TPC-only photon tuning
+  int mTPCTrackMinNClusters = 25;                      // minimum number of clusters
   float mTPCTrackXY2Radius = 21.f;                     // check for conversions close to material
   float mTPCTrackD2R = 30.f;                           // check for maximal distance of pairs and their radii
   float minTPCdEdx = 250;                              // starting from this dEdx value, tracks with p > minMomTPCdEdx are always accepted

diff --git a/Detectors/Vertexing/src/SVertexer.cxx b/Detectors/Vertexing/src/SVertexer.cxx
@@ -1247,6 +1247,8 @@ bool SVertexer::processTPCTrack(const o2::tpc::TrackTPC& trTPC, GIndex gid, int
 
   if (mSVParams->mTPCTrackPhotonTune) {
     bool discard = false;
+    // require minimum of tpc clusters
+    discard = trTPC.getNClusters() < mSVParams->mTPCTrackMinNClusters;
     // check track z cuts
     discard = std::abs(trLoc.getX() * trLoc.getTgl() - trLoc.getZ() - vtx.getZ()) > mSVParams->mTPCTrack2Beam;
     // check track transveres cuts