use tgl instead of eta

f3sch · Nov 25, 2023 · 708e852 · 708e852
1 parent 6ee3adf
commit 708e852
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 6 deletions.
diff --git a/Detectors/Vertexing/include/DetectorsVertexing/SVertexerParams.h b/Detectors/Vertexing/include/DetectorsVertexing/SVertexerParams.h
@@ -59,8 +59,8 @@ struct SVertexerParams : public o2::conf::ConfigurableParamHelper<SVertexerParam
   float maxDCAXYToMeanVertex = 0.2;       ///< max DCA of V0 from beam line (mean vertex) for prompt V0 candidates
   float maxDCAXYToMeanVertexV0Casc = 2.;  ///< max DCA of V0 from beam line (mean vertex) for cascade V0 candidates
   float maxDCAXYToMeanVertex3bodyV0 = 2.; ///< max DCA of V0 from beam line (mean vertex) for 3body V0 candidates
-  float maxV0EtaAbsDiff = 1.;             ///< max absolute difference in eta for V0 before rotate into same Frame and x
-  float maxV0PhiAbsDiff = 2.;             ///< max absolute difference in phi for V0 before rotate into same Frame and x
+  float maxV0TglAbsDiff = 0.4;            ///< max absolute difference in eta for V0 before rotate into same Frame and x
+  float maxV0PhiAbsDiff = 1.4;            ///< max absolute difference in phi for V0 before rotate into same Frame and x
   float minPtV0 = 0.01;                   ///< v0 minimum pT
   float minPtV0FromCascade = 0.3;         ///< v0 minimum pT for v0 to be used in cascading (lowest pT Run 2 lambda: 0.4)
   float maxTglV0 = 2.;                    ///< maximum tgLambda of V0

diff --git a/Detectors/Vertexing/src/SVertexer.cxx b/Detectors/Vertexing/src/SVertexer.cxx
@@ -576,9 +576,9 @@ bool SVertexer::checkV0(const TrackCand& seedP, const TrackCand& seedN, int iP,
   // Fast rough cuts on pairs before feeding to DCAFitter, tracks are not in the same Frame or at same X
   bool usesTPCOnly = (seedP.hasTPC && seedP.nITSclu < 1) || (seedN.hasTPC && seedN.nITSclu < 1);
   if (usesTPCOnly && mSVParams->mTPCTrackPhotonTune) {
-    auto deta = std::abs(seedP.getEta() - seedN.getEta()), dphi = std::abs(seedP.getPhi() - seedN.getPhi());
-    if (deta > mSVParams->maxV0EtaAbsDiff || dphi > mSVParams->maxV0PhiAbsDiff) {
-      LOG(debug) << "RejEtaPhi " << deta << " " << dphi;
+    auto dTgl = std::abs(seedP.getTgl() - seedN.getTgl()), dPhi = std::abs(seedP.getPhi() - seedN.getPhi());
+    if (dTgl > mSVParams->maxV0TglAbsDiff || dPhi > mSVParams->maxV0PhiAbsDiff) {
+      LOG(debug) << "RejTglPhi " << dTgl << " " << dPhi;
       return false;
     }
     float sna, csa;
@@ -645,7 +645,7 @@ bool SVertexer::checkV0(const TrackCand& seedP, const TrackCand& seedN, int iP,
 
   bool goodHyp = false;
   std::array<bool, NHypV0> hypCheckStatus{};
-  for (int ipid = 0; ipid < NHypV0 && mSVParams->checkV0Hypothesis; ipid++) {
+  for (int ipid = 0; (ipid < ((usesTPCOnly) ? NHypV0 : (Photon + 1))) && mSVParams->checkV0Hypothesis; ipid++) {
     if (mV0Hyps[ipid].check(p2Pos, p2Neg, p2V0, ptV0)) {
       goodHyp = hypCheckStatus[ipid] = true;
     }