diff --git a/PWGLF/Tasks/Strangeness/phiStrangeCorrelation.cxx b/PWGLF/Tasks/Strangeness/phiStrangeCorrelation.cxx index da93877f57f..7226974f530 100644 --- a/PWGLF/Tasks/Strangeness/phiStrangeCorrelation.cxx +++ b/PWGLF/Tasks/Strangeness/phiStrangeCorrelation.cxx @@ -78,11 +78,13 @@ static constexpr std::array phiMassRegionLabels{"Signal", " enum ParticleOfInterest { Phi = 0, K0S, - PionTPC, - PionTPCTOF + Pion, + /*PionTPC, + PionTPCTOF*/ + ParticleOfInterestSize }; -static constexpr std::array particleOfInterestLabels{"Phi", "K0S", "PionTPC", "PionTPCTOF"}; +static constexpr std::array particleOfInterestLabels{"Phi", "K0S", "Pion" /*"PionTPC", "PionTPCTOF"*/}; /* #define LIST_OF_PARTICLES_OF_INTEREST \ @@ -180,8 +182,8 @@ struct PhiStrangenessCorrelation { Configurable forceTOF{"forceTOF", false, "force the TOF signal for the PID"}; Configurable tofPIDThreshold{"tofPIDThreshold", 0.5, "minimum pT after which TOF PID is applicable"}; Configurable> trkPIDspecies{"trkPIDspecies", std::vector{o2::track::PID::Pion, o2::track::PID::Kaon, o2::track::PID::Proton}, "Trk sel: Particles species for PID, proton, pion, kaon"}; - Configurable> pidTPCMax{"pidTPCMax", std::vector{2., 0., 0.}, "maximum nSigma TPC"}; - Configurable> pidTOFMax{"pidTOFMax", std::vector{2., 0., 0.}, "maximum nSigma TOF"}; + Configurable> pidTPCMax{"pidTPCMax", std::vector{2.0f, 2.0f, 2.0f}, "maximum nSigma TPC"}; + Configurable> pidTOFMax{"pidTOFMax", std::vector{2.0f, 2.0f, 2.0f}, "maximum nSigma TOF"}; } trackConfigs; // Configurables on phi selection @@ -261,7 +263,7 @@ struct PhiStrangenessCorrelation { using V0DauMCTracks = soa::Join; // Defining the type of the tracks for data and MC - using FullTracks = soa::Join; + using FullTracks = soa::Join; using FullMCTracks = soa::Join; // using FilteredTracks = soa::Filtered; @@ -286,7 +288,7 @@ struct PhiStrangenessCorrelation { std::shared_ptr effMapPionTPC = nullptr; std::shared_ptr effMapPionTPCTOF = nullptr;*/ - std::array, 4> effMaps{}; + std::array, ParticleOfInterestSize> effMaps{}; void init(InitContext&) { @@ -344,7 +346,7 @@ struct PhiStrangenessCorrelation { ccdb->setLocalObjectValidityChecking(); ccdb->setFatalWhenNull(false); - for (int i = 0; i < 4; ++i) { + for (int i = 0; i < ParticleOfInterestSize; ++i) { loadEfficiencyMapFromCCDB(static_cast(i)); } } @@ -443,7 +445,7 @@ struct PhiStrangenessCorrelation { if (std::abs(nSigmaTPC) >= trackConfigs.pidTPCMax->at(speciesIndex)) { return false; // TPC check failed } - if (trackConfigs.forceTOF || (track.pt() > trackConfigs.tofPIDThreshold && track.hasTOF())) { + if (trackConfigs.forceTOF || (track.pt() >= trackConfigs.tofPIDThreshold && track.hasTOF())) { auto nSigmaTOF = aod::pidutils::tofNSigma(pid, track); if (std::abs(nSigmaTOF) >= trackConfigs.pidTOFMax->at(speciesIndex)) { return false; // TOF check failed @@ -553,7 +555,7 @@ struct PhiStrangenessCorrelation { if (!selectionPion(track)) continue; - auto Pion = track.pt() < trackConfigs.tofPIDThreshold ? PionTPC : PionTPCTOF; + // auto Pion = track.pt() < trackConfigs.tofPIDThreshold ? PionTPC : PionTPCTOF; float weightPhiPion = computeWeight(BoundEfficiencyMap(effMaps[Phi], multiplicity, phiCand.pt(), phiCand.y()), BoundEfficiencyMap(effMaps[Pion], multiplicity, track.pt(), track.rapidity(massPi)));