[PWGHF] taskCharmPolarisation: Fill the THsparse when not apply rotation background (#15935)

Mingze129 · web-flow · commit 0174bdc46bfa · 2026-04-23T16:23:52.000+02:00
diff --git a/PWGHF/D2H/Tasks/taskCharmPolarisation.cxx b/PWGHF/D2H/Tasks/taskCharmPolarisation.cxx
@@ -1626,6 +1626,9 @@ struct HfTaskCharmPolarisation {
                 }
               } else {
                 registry.fill(HIST("hEP"), invMassCharmHad, ptCharmHad, numPvContributors, std::abs(rapCharmHad), invMassD0, cosThetaStar, outputMl[0], /*outputMl[1],*/ outputMl[2], absEtaMin, numItsClsMin, numTpcClsMin, centrality);
+                if (activateTHnSparsePhiAndCosSqStarEP && ptCharmHad > ptMinCosPhiCosSqTheta) {
+                  registry.fill(HIST("hCosPhiCosSqThetaEP"), invMassCharmHad, ptCharmHad, std::abs(rapCharmHad), std::cos(2 * phiEuler), cosThetaStar * cosThetaStar, outputMl[0], /*outputMl[1],*/ outputMl[2], centrality);
+                }
               }
             } else {
               if (nBkgRotations > 0) {
@@ -1635,6 +1638,9 @@ struct HfTaskCharmPolarisation {
                 }
               } else {
                 registry.fill(HIST("hEP"), invMassCharmHad, ptCharmHad, numPvContributors, std::abs(rapCharmHad), invMassD0, cosThetaStar, outputMl[0], /*outputMl[1],*/ outputMl[2], centrality);
+                if (activateTHnSparsePhiAndCosSqStarEP && ptCharmHad > ptMinCosPhiCosSqTheta) {
+                  registry.fill(HIST("hCosPhiCosSqThetaEP"), invMassCharmHad, ptCharmHad, std::abs(rapCharmHad), std::cos(2 * phiEuler), cosThetaStar * cosThetaStar, outputMl[0], /*outputMl[1],*/ outputMl[2], centrality);
+                }
               }
             }
           }

Original file line number	Diff line number	Diff line change
`@@ -1626,6 +1626,9 @@ struct HfTaskCharmPolarisation {`
`1626`	`1626`	`}`
`1627`	`1627`	`} else {`
`1628`	`1628`	`registry.fill(HIST("hEP"), invMassCharmHad, ptCharmHad, numPvContributors, std::abs(rapCharmHad), invMassD0, cosThetaStar, outputMl[0], /outputMl[1],/ outputMl[2], absEtaMin, numItsClsMin, numTpcClsMin, centrality);`
	`1629`	`+ if (activateTHnSparsePhiAndCosSqStarEP && ptCharmHad > ptMinCosPhiCosSqTheta) {`
	`1630`	`+ registry.fill(HIST("hCosPhiCosSqThetaEP"), invMassCharmHad, ptCharmHad, std::abs(rapCharmHad), std::cos(2 * phiEuler), cosThetaStar * cosThetaStar, outputMl[0], /outputMl[1],/ outputMl[2], centrality);`
	`1631`	`+ }`
`1629`	`1632`	`}`
`1630`	`1633`	`} else {`
`1631`	`1634`	`if (nBkgRotations > 0) {`
`@@ -1635,6 +1638,9 @@ struct HfTaskCharmPolarisation {`
`1635`	`1638`	`}`
`1636`	`1639`	`} else {`
`1637`	`1640`	`registry.fill(HIST("hEP"), invMassCharmHad, ptCharmHad, numPvContributors, std::abs(rapCharmHad), invMassD0, cosThetaStar, outputMl[0], /outputMl[1],/ outputMl[2], centrality);`
	`1641`	`+ if (activateTHnSparsePhiAndCosSqStarEP && ptCharmHad > ptMinCosPhiCosSqTheta) {`
	`1642`	`+ registry.fill(HIST("hCosPhiCosSqThetaEP"), invMassCharmHad, ptCharmHad, std::abs(rapCharmHad), std::cos(2 * phiEuler), cosThetaStar * cosThetaStar, outputMl[0], /outputMl[1],/ outputMl[2], centrality);`
	`1643`	`+ }`
`1638`	`1644`	`}`
`1639`	`1645`	`}`
`1640`	`1646`	`}`