Skip to content

Commit c33782e

Browse files
[PWGHF] Fix charm-hadron EsE filling (#17323)
1 parent 55ada5f commit c33782e

1 file changed

Lines changed: 6 additions & 8 deletions

File tree

PWGHF/D2H/Tasks/taskFlowCharmHadrons.cxx

Lines changed: 6 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -608,8 +608,8 @@ struct HfTaskFlowCharmHadrons {
608608
if constexpr (StoreInfo == RunMode::kEsE || StoreInfo == RunMode::kSPEsE) {
609609
qVecRedComps = getEseQvec(collision, qVecRedDetector.value);
610610
}
611-
float xRedQVec = qVecRedComps[0];
612-
float yRedQVec = qVecRedComps[1];
611+
float const xRedQVec = qVecRedComps[0];
612+
float const yRedQVec = qVecRedComps[1];
613613
float const amplRedQVec = qVecRedComps[2];
614614

615615
for (const auto& candidate : candidates) {
@@ -811,12 +811,10 @@ struct HfTaskFlowCharmHadrons {
811811
// subtract daughters' contribution from the (normalized) Q-vector
812812
const float redQVecXDaugSubtr = xRedQVec - std::accumulate(tracksRedQx.begin(), tracksRedQx.end(), 0.0);
813813
const float redQVecYDaugSubtr = yRedQVec - std::accumulate(tracksRedQy.begin(), tracksRedQy.end(), 0.0);
814-
if (qVecRedDetector.value == QvecEstimator::TPCTot || qVecRedDetector.value == QvecEstimator::TPCPos || qVecRedDetector.value == QvecEstimator::TPCNeg) {
815-
// Correct for track multiplicity
816-
redQVec = std::hypot(redQVecXDaugSubtr, redQVecYDaugSubtr) * std::sqrt(amplRedQVec) / std::sqrt(amplRedQVec - tracksRedQx.size());
817-
} else {
818-
redQVec = std::hypot(xRedQVec, yRedQVec);
819-
}
814+
// Correct for track multiplicity
815+
redQVec = std::hypot(redQVecXDaugSubtr, redQVecYDaugSubtr) * std::sqrt(amplRedQVec) / std::sqrt(amplRedQVec - tracksRedQx.size());
816+
} else {
817+
redQVec = std::hypot(xRedQVec, yRedQVec);
820818
}
821819
if (storeRedQVec) {
822820
rowCandFlowEsE(massCand, ptCand, outputMl[0], outputMl[1], scalprodCand, cent, redQVec);

0 commit comments

Comments
 (0)