From 6ce4d2969df35e837e7d41d3b8081a3628c12f4d Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Sat, 20 Apr 2024 20:08:21 +0200 Subject: [PATCH 1/2] PWGHF: fix typo in soft pion rotation --- PWGHF/D2H/Tasks/taskCharmPolarisation.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PWGHF/D2H/Tasks/taskCharmPolarisation.cxx b/PWGHF/D2H/Tasks/taskCharmPolarisation.cxx index 3ff9529b43f..beeb137c8f3 100644 --- a/PWGHF/D2H/Tasks/taskCharmPolarisation.cxx +++ b/PWGHF/D2H/Tasks/taskCharmPolarisation.cxx @@ -562,7 +562,7 @@ struct TaskPolarisationCharmHadrons { massDau = massPi; // (*) const float bkgRotAngle = bkgRotationAngleStep * bkgRotationId; - std::array threeVecSoftPi{candidate.pxSoftPi() * std::cos(bkgRotAngle) - candidate.pxSoftPi() * std::sin(bkgRotAngle), candidate.pxSoftPi() * std::cos(bkgRotAngle) + candidate.pxSoftPi() * std::sin(bkgRotAngle), candidate.pzSoftPi()}; // we rotate the soft pion + std::array threeVecSoftPi{candidate.pxSoftPi() * std::cos(bkgRotAngle) - candidate.pySoftPi() * std::sin(bkgRotAngle), candidate.pxSoftPi() * std::cos(bkgRotAngle) + candidate.pySoftPi() * std::sin(bkgRotAngle), candidate.pzSoftPi()}; // we rotate the soft pion std::array threeVecD0Prong0{candidate.pxProng0(), candidate.pxProng0(), candidate.pzProng0()}; std::array threeVecD0Prong1{candidate.pxProng1(), candidate.pyProng1(), candidate.pzProng1()}; if (bkgRotationId > 0) { From fae46e23021cd6a931801ed64ef8d425b3f4f4f6 Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Sat, 20 Apr 2024 20:10:43 +0200 Subject: [PATCH 2/2] Fix second bug --- PWGHF/D2H/Tasks/taskCharmPolarisation.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PWGHF/D2H/Tasks/taskCharmPolarisation.cxx b/PWGHF/D2H/Tasks/taskCharmPolarisation.cxx index beeb137c8f3..cbbe2904cd4 100644 --- a/PWGHF/D2H/Tasks/taskCharmPolarisation.cxx +++ b/PWGHF/D2H/Tasks/taskCharmPolarisation.cxx @@ -562,7 +562,7 @@ struct TaskPolarisationCharmHadrons { massDau = massPi; // (*) const float bkgRotAngle = bkgRotationAngleStep * bkgRotationId; - std::array threeVecSoftPi{candidate.pxSoftPi() * std::cos(bkgRotAngle) - candidate.pySoftPi() * std::sin(bkgRotAngle), candidate.pxSoftPi() * std::cos(bkgRotAngle) + candidate.pySoftPi() * std::sin(bkgRotAngle), candidate.pzSoftPi()}; // we rotate the soft pion + std::array threeVecSoftPi{candidate.pxSoftPi() * std::cos(bkgRotAngle) - candidate.pySoftPi() * std::sin(bkgRotAngle), candidate.pxSoftPi() * std::sin(bkgRotAngle) + candidate.pySoftPi() * std::cos(bkgRotAngle), candidate.pzSoftPi()}; // we rotate the soft pion std::array threeVecD0Prong0{candidate.pxProng0(), candidate.pxProng0(), candidate.pzProng0()}; std::array threeVecD0Prong1{candidate.pxProng1(), candidate.pyProng1(), candidate.pzProng1()}; if (bkgRotationId > 0) {