diff --git a/PWGEM/Dilepton/TableProducer/skimmerPrimaryMFTTrack.cxx b/PWGEM/Dilepton/TableProducer/skimmerPrimaryMFTTrack.cxx index 28055a15057..977061cee75 100644 --- a/PWGEM/Dilepton/TableProducer/skimmerPrimaryMFTTrack.cxx +++ b/PWGEM/Dilepton/TableProducer/skimmerPrimaryMFTTrack.cxx @@ -198,13 +198,16 @@ struct skimmerPrimaryMFTTrack { // As minimal cuts, following cuts are applied. The cut values are hardcoded on the purpose for consistent bit operation. // Ncls MFT >= 5 - // chi2/ndf MFT < 5 + // chi2/ndf MFT < 4 // |dcaXY| < 0.06 cm - if (mfttrack.nClusters() < 6 || mfttrack.chi2() / ndf > 5.f || std::fabs(dcaXY) > 0.05) { + if (mfttrack.nClusters() < 5 || mfttrack.chi2() / ndf > 4.f || std::fabs(dcaXY) > 0.06) { return; } + if (mfttrack.nClusters() >= 6) { + trackBit |= static_cast(RefMFTTrackBit::kNclsMFT6); + } if (mfttrack.nClusters() >= 7) { trackBit |= static_cast(RefMFTTrackBit::kNclsMFT7); } @@ -212,13 +215,16 @@ struct skimmerPrimaryMFTTrack { trackBit |= static_cast(RefMFTTrackBit::kNclsMFT8); } - if (mfttrack.chi2() / ndf < 4.f) { - trackBit |= static_cast(RefMFTTrackBit::kChi2MFT4); - } if (mfttrack.chi2() / ndf < 3.f) { trackBit |= static_cast(RefMFTTrackBit::kChi2MFT3); } + if (mfttrack.chi2() / ndf < 2.f) { + trackBit |= static_cast(RefMFTTrackBit::kChi2MFT2); + } + if (std::fabs(dcaXY) < 0.05) { + trackBit |= static_cast(RefMFTTrackBit::kDCAxy005cm); + } if (std::fabs(dcaXY) < 0.04) { trackBit |= static_cast(RefMFTTrackBit::kDCAxy004cm); } diff --git a/PWGEM/Dilepton/Utils/EMTrackUtilities.h b/PWGEM/Dilepton/Utils/EMTrackUtilities.h index 54bc28eb012..88daecccb62 100644 --- a/PWGEM/Dilepton/Utils/EMTrackUtilities.h +++ b/PWGEM/Dilepton/Utils/EMTrackUtilities.h @@ -47,14 +47,16 @@ enum class RefTrackBit : uint16_t { // This is not for leptons, but charged trac }; enum class RefMFTTrackBit : uint16_t { // This is not for leptons, but charged tracks for reference flow. - kNclsMFT7 = 1, // default is 6 - kNclsMFT8 = 2, - kChi2MFT4 = 4, // default is 5 - kChi2MFT3 = 8, - kDCAxy004cm = 16, // default is 0.05 cm - kDCAxy003cm = 32, - kDCAxy002cm = 64, - kDCAxy001cm = 128, + kNclsMFT6 = 1, // default is 5 + kNclsMFT7 = 2, + kNclsMFT8 = 4, + kChi2MFT3 = 8, // default is 4 + kChi2MFT2 = 16, + kDCAxy005cm = 32, // default is 0.06 cm + kDCAxy004cm = 64, + kDCAxy003cm = 128, + kDCAxy002cm = 256, + kDCAxy001cm = 512, }; //_______________________________________________________________________