From 90e9aa25c2e2d01094876ab3b7d8b62e532b5005 Mon Sep 17 00:00:00 2001 From: shahoian Date: Fri, 23 May 2025 17:39:22 +0200 Subject: [PATCH 1/2] Fix condition to disable V0/Cascade hypothesis in the SVertexer --- .../Vertexing/include/DetectorsVertexing/SVertexHypothesis.h | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/Detectors/Vertexing/include/DetectorsVertexing/SVertexHypothesis.h b/Detectors/Vertexing/include/DetectorsVertexing/SVertexHypothesis.h index c3fd74aa7eeff..bc6673288f752 100644 --- a/Detectors/Vertexing/include/DetectorsVertexing/SVertexHypothesis.h +++ b/Detectors/Vertexing/include/DetectorsVertexing/SVertexHypothesis.h @@ -60,7 +60,7 @@ class SVertexHypothesis bool check(float p2Pos, float p2Neg, float p2V0, float ptV0) const { // check if given mass and pt is matching to hypothesis - return mPars[SigmaM] > 0 && check(calcMass(p2Pos, p2Neg, p2V0), ptV0); + return mPars[SigmaM] >= 0.f && check(calcMass(p2Pos, p2Neg, p2V0), ptV0); } bool check(float mass, float pt) const { // check if given mass and pt is matching to hypothesis @@ -151,7 +151,7 @@ class SVertex3Hypothesis bool check(float p2Pos, float p2Neg, float p2Bach, float p2Tot, float ptV0) const { // check if given mass and pt is matching to hypothesis - return mPars[SigmaM] > 0 && check(calcMass(p2Pos, p2Neg, p2Bach, p2Tot), ptV0); + return mPars[SigmaM] >= 0.f && check(calcMass(p2Pos, p2Neg, p2Bach, p2Tot), ptV0); } bool check(float mass, float pt) const From 8cd9960a6f4142bdea76f0846233bc9e0df06c32 Mon Sep 17 00:00:00 2001 From: shahoian Date: Fri, 23 May 2025 17:40:46 +0200 Subject: [PATCH 2/2] Fix some rejection status assignments --- Common/DCAFitter/include/DCAFitter/DCAFitterN.h | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/Common/DCAFitter/include/DCAFitter/DCAFitterN.h b/Common/DCAFitter/include/DCAFitter/DCAFitterN.h index aac451f1f8978..df732bd4bde63 100644 --- a/Common/DCAFitter/include/DCAFitter/DCAFitterN.h +++ b/Common/DCAFitter/include/DCAFitter/DCAFitterN.h @@ -142,6 +142,8 @@ class DCAFitterN FailInv2ndDeriv, // inversion of 2nd derivatives failed FailCorrTracks, // correction of tracks to updated x failed FailCloserAlt, // alternative PCA is closer + // + NStatusesDefined }; static constexpr int getNProngs() { return N; } @@ -989,7 +991,7 @@ GPUd() bool DCAFitterN::minimizeChi2() } if (mMaxDZIni > 0 && !roughDZCut()) { // apply rough cut on tracks Z difference - mFitStatus[mCurHyp] = FitStatus::RejTrackX; + mFitStatus[mCurHyp] = FitStatus::RejTrackRoughZ; return false; } @@ -1063,7 +1065,7 @@ GPUd() bool DCAFitterN::minimizeChi2NoErr() setTrackPos(mTrPos[mCurHyp][i], mCandTr[mCurHyp][i]); // prepare positions } if (mMaxDZIni > 0 && !roughDZCut()) { // apply rough cut on tracks Z difference - mFitStatus[mCurHyp] = FitStatus::RejTrackX; + mFitStatus[mCurHyp] = FitStatus::RejTrackRoughZ; return false; }