Skip to content

Commit 401ea7a

Browse files
authored
[PWGHF,Trigger] Revert to symmetric NsigmaDe cut for charm femto trigger (#11353)
1 parent 4ec9c7a commit 401ea7a

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

EventFiltering/PWGHF/HFFilterHelpers.h

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -278,7 +278,7 @@ constexpr float massJPsi = o2::constants::physics::MassJPsi;
278278

279279
static const o2::framework::AxisSpec ptAxis{50, 0.f, 50.f};
280280
static const o2::framework::AxisSpec pAxis{50, 0.f, 10.f};
281-
static const o2::framework::AxisSpec kstarAxis{100, 0.f, 1.f};
281+
static const o2::framework::AxisSpec kstarAxis{200, 0.f, 2.f};
282282
static const o2::framework::AxisSpec etaAxis{30, -1.5f, 1.5f};
283283
static const o2::framework::AxisSpec nSigmaAxis{100, -10.f, 10.f};
284284
static const o2::framework::AxisSpec alphaAxis{100, -1.f, 1.f};
@@ -958,19 +958,19 @@ inline bool HfFilterHelper::isSelectedTrack4Femto(const T1& track, const T2& tra
958958
// Apply different PID strategy in different pt range
959959
// one side selection only
960960
if (pt <= ptThresholdPidStrategy) {
961-
if (NSigmaTPC < -nSigmaCuts[0] || NSigmaITS < -nSigmaCuts[3]) { // Use TPC and ITS below the threshold, NSigmaITS for deuteron with a lower limit
961+
if (std::fabs(NSigmaTPC) > nSigmaCuts[0] || NSigmaITS < -nSigmaCuts[3]) { // Use TPC and ITS below the threshold, NSigmaITS for deuteron with a lower limit
962962
return false;
963963
}
964964
} else {
965-
if (NSigmaTOF < -nSigmaCuts[1] || NSigmaTPC < -nSigmaCuts[0]) { // Use combined TPC and TOF above the threshold
965+
if (std::fabs(NSigmaTOF) > nSigmaCuts[1] || std::fabs(NSigmaTPC) > nSigmaCuts[0]) { // Use combined TPC and TOF above the threshold
966966
return false;
967967
}
968968
}
969969
}
970970

971971
if (activateQA > 1) {
972972
hTPCPID->Fill(track.p(), NSigmaTPC);
973-
if ((forceTof || track.hasTOF())) {
973+
if ((!forceTof || track.hasTOF())) {
974974
if (trackSpecies == kProtonForFemto)
975975
hTOFPID->Fill(track.p(), NSigmaTOF);
976976
else if (trackSpecies == kDeuteronForFemto && pt > ptThresholdPidStrategy)

0 commit comments

Comments
 (0)