@@ -123,13 +123,13 @@ std::tuple<std::vector<int>, std::vector<int>> MatchJetsGeometricallyImpl(
123123 const std::vector<T>& jetsBaseEta,
124124 std::vector<T> jetsBasePhiForMatching,
125125 std::vector<T> jetsBaseEtaForMatching,
126- const std::vector<std::size_t > jetMapBaseToJetIndex,
126+ const std::vector<std::size_t >& jetMapBaseToJetIndex,
127127 const std::vector<T>& jetsTagPhi,
128128 const std::vector<T>& jetsTagEta,
129129 std::vector<T> jetsTagPhiForMatching,
130130 std::vector<T> jetsTagEtaForMatching,
131- const std::vector<std::size_t > jetMapTagToJetIndex,
132- double maxMatchingDistance)
131+ const std::vector<std::size_t >& jetMapTagToJetIndex,
132+ const double maxMatchingDistance)
133133{
134134 // Validation
135135 // If no jets in either collection, then return immediately.
@@ -544,8 +544,10 @@ void MatchPt(T const& jetsBasePerCollision, U const& jetsTagPerCollision, std::v
544544 auto jetTagClusters = getConstituents (jetTag, clustersTag);
545545 auto jetTagCandidates = getConstituents (jetTag, candidatesTag);
546546
547- ptSumBase = getPtSum < jetfindingutilities::isEMCALClusterTable<N>() || jetfindingutilities::isEMCALClusterTable<Q>(), (jetcandidateutilities::isCandidateTable<M>() || jetcandidateutilities::isCandidateMcTable<M>()) && (jetcandidateutilities::isCandidateTable<P>() || jetcandidateutilities::isCandidateMcTable<P>()), jetsBaseIsMc, jetsTagIsMc > (jetBaseTracks, jetBaseCandidates, jetBaseClusters, jetTagTracks, jetTagCandidates, jetTagClusters, tracksBase, tracksTag);
548- ptSumTag = getPtSum < jetfindingutilities::isEMCALClusterTable<N>() || jetfindingutilities::isEMCALClusterTable<Q>(), (jetcandidateutilities::isCandidateTable<M>() || jetcandidateutilities::isCandidateMcTable<M>()) && (jetcandidateutilities::isCandidateTable<P>() || jetcandidateutilities::isCandidateMcTable<P>()), jetsTagIsMc, jetsBaseIsMc > (jetTagTracks, jetTagCandidates, jetTagClusters, jetBaseTracks, jetBaseCandidates, jetBaseClusters, tracksTag, tracksBase);
547+ constexpr bool IsEMCAL{jetfindingutilities::isEMCALClusterTable<N>() || jetfindingutilities::isEMCALClusterTable<Q>()};
548+ constexpr bool IsCandidate{(jetcandidateutilities::isCandidateTable<M>() || jetcandidateutilities::isCandidateMcTable<M>()) && (jetcandidateutilities::isCandidateTable<P>() || jetcandidateutilities::isCandidateMcTable<P>())};
549+ ptSumBase = getPtSum<IsEMCAL, IsCandidate, jetsBaseIsMc, jetsTagIsMc>(jetBaseTracks, jetBaseCandidates, jetBaseClusters, jetTagTracks, jetTagCandidates, jetTagClusters, tracksBase, tracksTag);
550+ ptSumTag = getPtSum<IsEMCAL, IsCandidate, jetsTagIsMc, jetsBaseIsMc>(jetTagTracks, jetTagCandidates, jetTagClusters, jetBaseTracks, jetBaseCandidates, jetBaseClusters, tracksTag, tracksBase);
549551 if (ptSumBase > jetBase.pt () * minPtFraction) {
550552 baseToTagMatchingPt[jetBase.globalIndex ()].push_back (jetTag.globalIndex ());
551553 }
0 commit comments