diff --git a/PWGEM/Dilepton/Utils/ElectronModule.h b/PWGEM/Dilepton/Utils/ElectronModule.h index 9ccb58e6299..a1a259163a9 100644 --- a/PWGEM/Dilepton/Utils/ElectronModule.h +++ b/PWGEM/Dilepton/Utils/ElectronModule.h @@ -697,6 +697,10 @@ class ElectronModule void fillMapMLPID(TCollision const& collision, TTrack const& track) { if (fElectronCut.usePIDML) { + if (!isElectron_TOFif(track, collision, fElectronCut)) { // minimal n sigma cut is taken from the main electron cut. + return; + } + o2::dataformats::DCA mDcaInfoCov; mDcaInfoCov.set(999, 999, 999, 999, 999); auto trackParCov = getTrackParCov(track); @@ -706,10 +710,6 @@ class ElectronModule return; } - if (!isElectron_TOFif(track, collision, fElectronCut)) { // minimal n sigma cut is taken from the main electron cut. - return; - } - o2::analysis::pwgem::dilepton::mlpid::candidate candidate; candidate.tpcInnerParam = track.tpcInnerParam(); candidate.tpcNClsFound = track.tpcNClsFound(); @@ -726,7 +726,7 @@ class ElectronModule } candidate.meanClusterSizeITSobCosTgl = static_cast(total_cluster_size_ob) / static_cast(nl_ob) * std::cos(std::atan(trackParCov.getTgl())); - std::vector inputFeatures = mlResponsePID.getInputFeatures(candidate); + auto inputFeatures = mlResponsePID.getInputFeatures(candidate); float binningFeature = mlResponsePID.getBinningFeature(candidate); int pbin = lower_bound(fElectronCut.binsMl.value.begin(), fElectronCut.binsMl.value.end(), binningFeature) - fElectronCut.binsMl.value.begin() - 1; @@ -1671,7 +1671,7 @@ class ElectronModule candidate.signedMassLH = electron.sign() * hadron.sign() * eTpair.mass; candidate.tpcNSigmaKa = hadron.tpcNSigmaKa(); - std::vector inputFeatures = mlResponseSCTeT.getInputFeatures(candidate); + auto inputFeatures = mlResponseSCTeT.getInputFeatures(candidate); float binningFeature = mlResponseSCTeT.getBinningFeature(candidate); int pbin = lower_bound(fConfigDFeT.binsMl.value.begin(), fConfigDFeT.binsMl.value.end(), binningFeature) - fConfigDFeT.binsMl.value.begin() - 1; @@ -1721,7 +1721,7 @@ class ElectronModule candidate.ptL = trackParCov.getPt(); candidate.massH = o2::constants::physics::MassK0Short; - std::vector inputFeatures = mlResponseSCTeV0.getInputFeatures(candidate); + auto inputFeatures = mlResponseSCTeV0.getInputFeatures(candidate); float binningFeature = mlResponseSCTeV0.getBinningFeature(candidate); int pbin = lower_bound(fConfigDFeV0.binsMl.value.begin(), fConfigDFeV0.binsMl.value.end(), binningFeature) - fConfigDFeV0.binsMl.value.begin() - 1; @@ -1772,7 +1772,7 @@ class ElectronModule candidate.ptL = trackParCov.getPt(); candidate.massH = o2::constants::physics::MassLambda; - std::vector inputFeatures = mlResponseSCTeV0.getInputFeatures(candidate); + auto inputFeatures = mlResponseSCTeV0.getInputFeatures(candidate); float binningFeature = mlResponseSCTeV0.getBinningFeature(candidate); int pbin = lower_bound(fConfigDFeV0.binsMl.value.begin(), fConfigDFeV0.binsMl.value.end(), binningFeature) - fConfigDFeV0.binsMl.value.begin() - 1; @@ -1823,7 +1823,7 @@ class ElectronModule candidate.ptL = trackParCov.getPt(); candidate.massH = o2::constants::physics::MassXiMinus; - std::vector inputFeatures = mlResponseSCTeC.getInputFeatures(candidate); + auto inputFeatures = mlResponseSCTeC.getInputFeatures(candidate); float binningFeature = mlResponseSCTeC.getBinningFeature(candidate); int pbin = lower_bound(fConfigDFeC.binsMl.value.begin(), fConfigDFeC.binsMl.value.end(), binningFeature) - fConfigDFeC.binsMl.value.begin() - 1; @@ -1874,7 +1874,7 @@ class ElectronModule candidate.ptL = trackParCov.getPt(); candidate.massH = o2::constants::physics::MassOmegaMinus; - std::vector inputFeatures = mlResponseSCTeC.getInputFeatures(candidate); + auto inputFeatures = mlResponseSCTeC.getInputFeatures(candidate); float binningFeature = mlResponseSCTeC.getBinningFeature(candidate); int pbin = lower_bound(fConfigDFeC.binsMl.value.begin(), fConfigDFeC.binsMl.value.end(), binningFeature) - fConfigDFeC.binsMl.value.begin() - 1; @@ -1925,7 +1925,7 @@ class ElectronModule candidate.ptL = trackParCov.getPt(); candidate.massH = o2::constants::physics::MassLambda; - std::vector inputFeatures = mlResponseSCTeV0.getInputFeatures(candidate); + auto inputFeatures = mlResponseSCTeV0.getInputFeatures(candidate); float binningFeature = mlResponseSCTeV0.getBinningFeature(candidate); int pbin = lower_bound(fConfigDFeV0.binsMl.value.begin(), fConfigDFeV0.binsMl.value.end(), binningFeature) - fConfigDFeV0.binsMl.value.begin() - 1; @@ -1976,7 +1976,7 @@ class ElectronModule candidate.ptL = trackParCov.getPt(); candidate.massH = o2::constants::physics::MassXiMinus; - std::vector inputFeatures = mlResponseSCTeC.getInputFeatures(candidate); + auto inputFeatures = mlResponseSCTeC.getInputFeatures(candidate); float binningFeature = mlResponseSCTeC.getBinningFeature(candidate); int pbin = lower_bound(fConfigDFeC.binsMl.value.begin(), fConfigDFeC.binsMl.value.end(), binningFeature) - fConfigDFeC.binsMl.value.begin() - 1; @@ -2027,7 +2027,7 @@ class ElectronModule candidate.ptL = trackParCov.getPt(); candidate.massH = o2::constants::physics::MassOmegaMinus; - std::vector inputFeatures = mlResponseSCTeC.getInputFeatures(candidate); + auto inputFeatures = mlResponseSCTeC.getInputFeatures(candidate); float binningFeature = mlResponseSCTeC.getBinningFeature(candidate); int pbin = lower_bound(fConfigDFeC.binsMl.value.begin(), fConfigDFeC.binsMl.value.end(), binningFeature) - fConfigDFeC.binsMl.value.begin() - 1;