diff --git a/Detectors/Vertexing/include/DetectorsVertexing/SVertexerParams.h b/Detectors/Vertexing/include/DetectorsVertexing/SVertexerParams.h index 4e4f66e1747ac..8221bbd1749db 100644 --- a/Detectors/Vertexing/include/DetectorsVertexing/SVertexerParams.h +++ b/Detectors/Vertexing/include/DetectorsVertexing/SVertexerParams.h @@ -95,11 +95,11 @@ struct SVertexerParams : public o2::conf::ConfigurableParamHelper minMomTPCdEdx are always accepted float minMomTPCdEdx = 0.8; // minimum p for tracks with dEdx > mMinTPCdEdx to be accepted - uint8_t mITSSAminNcluCascades = 6; // require at least this many ITS clusters if no TPC info present + uint8_t mITSSAminNcluCascades = 6; // require at least this many ITS clusters if no TPC info present bool mRequireTPCforCascBaryons = true; // require that baryon daughter of cascade has TPC - + // percent deviation from expected proton dEdx - to be replaced - estimated sigma from TPC for now 6%; a 6*sigma cut is therefore 36% = 0.36f. Any value above 1.0f will be ignored manually when checking. - float mFractiondEdxforCascBaryons = 0.36f; + float mFractiondEdxforCascBaryons = 0.36f; // cuts on different V0 PID params bool checkV0Hypothesis = true; diff --git a/Detectors/Vertexing/src/SVertexer.cxx b/Detectors/Vertexing/src/SVertexer.cxx index 216fd44a7e691..024d55f93c3ac 100644 --- a/Detectors/Vertexing/src/SVertexer.cxx +++ b/Detectors/Vertexing/src/SVertexer.cxx @@ -302,7 +302,7 @@ void SVertexer::updateTimeDependentParams() ft.setBz(bz); } - // TPC dE/dx usage information + // TPC dE/dx usage information // default: average 2023 from C. Sonnabend, Nov 2023: ([0.217553 4.02762 0.00850178 2.33324 0.880904 ]) // to-do: grab from CCDB when available -> discussion with TPC experts, not available yet const float bbPars[] = {0.217553, 4.02762, 0.00850178, 2.33324, 0.880904}; @@ -499,7 +499,7 @@ void SVertexer::buildT2V(const o2::globaltracking::RecoContainer& recoData) // a auto protonId = o2::track::PID::Proton; float dEdxExpected = mPIDresponse.getExpectedSignal(tpcTrack, protonId); float fracDevProton = std::abs((dEdxTPC - dEdxExpected) / dEdxExpected); - if(fracDevProtonmFractiondEdxforCascBaryons){ + if (fracDevProton < mSVParams->mFractiondEdxforCascBaryons) { compatibleWithProton = true; } } @@ -605,7 +605,7 @@ bool SVertexer::checkV0(const TrackCand& seedP, const TrackCand& seedN, int iP, } // check tight lambda mass only bool goodLamForCascade = false, goodALamForCascade = false; - if( ptV0 > mSVParams->minPtV0FromCascade && seedP.nITSclu >= mSVParams->mITSSAminNcluCascades && seedN.nITSclu >= mSVParams->mITSSAminNcluCascades ){ + if (ptV0 > mSVParams->minPtV0FromCascade && seedP.nITSclu >= mSVParams->mITSSAminNcluCascades && seedN.nITSclu >= mSVParams->mITSSAminNcluCascades) { if (mV0Hyps[Lambda].checkTight(p2Pos, p2Neg, p2V0, ptV0) && (!mSVParams->mRequireTPCforCascBaryons || seedP.hasTPC) && seedP.compatibleProton) { goodLamForCascade = true; } @@ -614,7 +614,6 @@ bool SVertexer::checkV0(const TrackCand& seedP, const TrackCand& seedN, int iP, } } - // apply mass selections for 3-body decay bool good3bodyV0Hyp = false; for (int ipid = 2; ipid < 4; ipid++) {