@@ -456,7 +456,7 @@ struct FemtoUniversePairTaskTrackPhi {
456456 }
457457
458458 template <bool isMC, typename PartitionType, typename PartType, typename MCParticles = std::nullptr_t >
459- void doSameEvent (const PartitionType& groupPartsTrack, const PartitionType& groupPartsPhi, const PartType& parts, float magFieldTesla, int multCol, [[maybe_unused]] MCParticles mcParts = nullptr )
459+ void doSameEvent (const PartitionType& groupPartsTrack, const PartitionType& groupPartsPhi, const PartType& parts, float magFieldTesla, int multCol, [[maybe_unused]] const MCParticles& mcParts = nullptr )
460460 {
461461 for (auto const & phicandidate : groupPartsPhi) {
462462 // TODO: add phi meson minv cut here
@@ -584,7 +584,7 @@ struct FemtoUniversePairTaskTrackPhi {
584584 }
585585
586586 template <bool isMC, typename PartitionType, typename PartType, typename MCParticles = std::nullptr_t >
587- void doMixedEvent (const PartitionType& groupPartsTrack, const PartitionType& groupPartsPhi, const PartType& parts, float magFieldTesla, int multCol, [[maybe_unused]] MCParticles mcParts = nullptr )
587+ void doMixedEvent (const PartitionType& groupPartsTrack, const PartitionType& groupPartsPhi, const PartType& parts, float magFieldTesla, int multCol, [[maybe_unused]] const MCParticles& mcParts = nullptr )
588588 {
589589 for (auto const & [track, phicandidate] : combinations (CombinationsFullIndexPolicy (groupPartsTrack, groupPartsPhi))) {
590590 if (confTrackIsIdentified) {
@@ -675,8 +675,9 @@ struct FemtoUniversePairTaskTrackPhi {
675675 void processMCTruth (aod::FDParticles const & parts)
676676 {
677677 for (auto const & part : parts) {
678- if (part.partType () != uint8_t (aod::femtouniverseparticle::ParticleType::kMCTruthTrack ))
678+ if (part.partType () != uint8_t (aod::femtouniverseparticle::ParticleType::kMCTruthTrack )) {
679679 continue ;
680+ }
680681
681682 int pdgCode = static_cast <int >(part.pidCut ());
682683 const auto & pdgParticle = pdgMC->GetParticle (pdgCode);
@@ -733,8 +734,9 @@ struct FemtoUniversePairTaskTrackPhi {
733734 {
734735 for (auto const & part : parts) {
735736 auto mcPartId = part.fdMCParticleId ();
736- if (mcPartId == kInvalidMCPartId )
737+ if (mcPartId == kInvalidMCPartId ) {
737738 continue ; // no MC particle
739+ }
738740 const auto & mcpart = mcparts.iteratorAt (mcPartId);
739741
740742 if (mcpart.pdgMCTruth () == confTrackPDGCode && (part.pt () > confTrackPtLow) && (part.pt () < confTrackPtHigh) && isParticleNSigmaAccepted (part.p (), trackCuts.getNsigmaTPC (part, o2::track::PID ::Proton), trackCuts.getNsigmaTOF (part, o2::track::PID ::Proton), trackCuts.getNsigmaTPC (part, o2::track::PID ::Pion), trackCuts.getNsigmaTOF (part, o2::track::PID ::Pion), trackCuts.getNsigmaTPC (part, o2::track::PID ::Kaon), trackCuts.getNsigmaTOF (part, o2::track::PID ::Kaon))) {
@@ -748,8 +750,9 @@ struct FemtoUniversePairTaskTrackPhi {
748750 registryMCpT.fill (HIST (" MCReco/C_phi_pT" ), part.pt (), weightPhi);
749751 }
750752
751- if (isParticleNSigmaAccepted (part.p (), trackCuts.getNsigmaTPC (part, o2::track::PID ::Proton), trackCuts.getNsigmaTOF (part, o2::track::PID ::Proton), trackCuts.getNsigmaTPC (part, o2::track::PID ::Pion), trackCuts.getNsigmaTOF (part, o2::track::PID ::Pion), trackCuts.getNsigmaTPC (part, o2::track::PID ::Kaon), trackCuts.getNsigmaTOF (part, o2::track::PID ::Kaon)))
753+ if (isParticleNSigmaAccepted (part.p (), trackCuts.getNsigmaTPC (part, o2::track::PID ::Proton), trackCuts.getNsigmaTOF (part, o2::track::PID ::Proton), trackCuts.getNsigmaTPC (part, o2::track::PID ::Pion), trackCuts.getNsigmaTOF (part, o2::track::PID ::Pion), trackCuts.getNsigmaTPC (part, o2::track::PID ::Kaon), trackCuts.getNsigmaTOF (part, o2::track::PID ::Kaon))) {
752754 hTrackDCA.fillQA <true , true >(part);
755+ }
753756 if ((part.partType () == aod::femtouniverseparticle::ParticleType::kPhi ) && (mcpart.pdgMCTruth () == o2::constants::physics::Pdg::kPhi ) && (mcpart.partOriginMCTruth () == aod::femtouniverse_mc_particle::ParticleOriginMCTruth::kPrimary )) {
754757 registryMCreco.fill (HIST (" MCrecoPhi" ), mcpart.pt (), mcpart.eta ()); // phi
755758 registryMCreco.fill (HIST (" MCrecoPhiPt" ), mcpart.pt ());
0 commit comments