@@ -357,7 +357,7 @@ struct HfTaskCharmPolarisation {
357357 const AxisSpec thnAxisCosThetaStarRandom{configThnAxisCosThetaStarRandom, " cos(#vartheta_{random})" };
358358 const AxisSpec thnAxisCosThetaStarBeam{configThnAxisCosThetaStarBeam, " cos(#vartheta_{beam})" };
359359 const AxisSpec thnAxisCosThetaStarEP{configThnAxisCosThetaStarEP, " cos(#vartheta_{EP})" }; // reaction plane
360- const AxisSpec thnAxisCosPhiStarEP{configThnAxisCosPhiStarEP, " cos(2#varphi_{EP})" }; // reaction plane
360+ const AxisSpec thnAxisCosPhiStarEP{configThnAxisCosPhiStarEP, " cos(2#varphi_{EP})" }; // reaction plane
361361 const AxisSpec thnAxisCosSqThetaStarEP{configThnAxisCosSqThetaStarEP, " cos^{2}(#vartheta_{EP})" }; // reaction plane
362362 const AxisSpec thnAxisMlBkg{configThnAxisMlBkg, " ML bkg" };
363363 const AxisSpec thnAxisMlNonPrompt{configThnAxisMlNonPrompt, " ML non-prompt" };
@@ -901,16 +901,17 @@ struct HfTaskCharmPolarisation {
901901 }; // end init
902902
903903 // / \param ptCharmHad is the pt of the candidate
904- double sampleDeltaPhi (float ptCharmHad) {
904+ double sampleDeltaPhi (float ptCharmHad)
905+ {
905906
906907 double y0 = 6.49e-2 ;
907- double cons = 8.42e-1 ;
908- double coef = -3.78e-1 ;
908+ double cons = 8.42e-1 ;
909+ double coef = -3.78e-1 ;
909910 double x0 = -6.37e-1 ;
910- double norm = -3.64e-1 ;
911- double mu = 4.68e-2 ;
911+ double norm = -3.64e-1 ;
912+ double mu = 4.68e-2 ;
912913 double sigma = 1.513 ;
913-
914+
914915 double sigmoid = y0 + (cons - y0) / (1 + std::exp (-coef * (ptCharmHad - x0)));
915916 double peak = norm * std::exp (-(ptCharmHad - mu) * (ptCharmHad - mu) / (2 * sigma * sigma));
916917
@@ -923,12 +924,13 @@ struct HfTaskCharmPolarisation {
923924
924925 double funcDeltaPhi = 1.0 + 2.0 * v2 * std::cos (2.0 * deltaPhi);
925926
926- if (funcVal < funcDeltaPhi) return deltaPhi;
927+ if (funcVal < funcDeltaPhi)
928+ return deltaPhi;
927929 }
928930
929931 return -1 .;
930932 }
931-
933+
932934 // / \param invMassCharmHad is the invariant-mass of the candidate
933935 // / \param ptCharmHad is the pt of the candidate
934936 // / \param numPvContributors is the number of PV contributors
@@ -1620,7 +1622,7 @@ struct HfTaskCharmPolarisation {
16201622 if (nBkgRotations > 0 ) {
16211623 registry.fill (HIST (" hEP" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], absEtaMin, numItsClsMin, numTpcClsMin, isRotatedCandidate, centrality);
16221624 if (activateTHnSparsePhiAndCosSqStarEP && ptCharmHad > ptMinCosPhiCosSqTheta) {
1623- registry.fill (HIST (" hCosPhiCosSqThetaEP" ), invMassCharmHad, ptCharmHad, std::abs (rapCharmHad), std::cos (2 * phiEuler), cosThetaStar* cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], centrality);
1625+ registry.fill (HIST (" hCosPhiCosSqThetaEP" ), invMassCharmHad, ptCharmHad, std::abs (rapCharmHad), std::cos (2 * phiEuler), cosThetaStar * cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], centrality);
16241626 }
16251627 } else {
16261628 registry.fill (HIST (" hEP" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], absEtaMin, numItsClsMin, numTpcClsMin, centrality);
@@ -1629,7 +1631,7 @@ struct HfTaskCharmPolarisation {
16291631 if (nBkgRotations > 0 ) {
16301632 registry.fill (HIST (" hEP" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], isRotatedCandidate, centrality);
16311633 if (activateTHnSparsePhiAndCosSqStarEP && ptCharmHad > ptMinCosPhiCosSqTheta) {
1632- registry.fill (HIST (" hCosPhiCosSqThetaEP" ), invMassCharmHad, ptCharmHad, std::abs (rapCharmHad), std::cos (2 * phiEuler), cosThetaStar* cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], centrality);
1634+ registry.fill (HIST (" hCosPhiCosSqThetaEP" ), invMassCharmHad, ptCharmHad, std::abs (rapCharmHad), std::cos (2 * phiEuler), cosThetaStar * cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], centrality);
16331635 }
16341636 } else {
16351637 registry.fill (HIST (" hEP" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], centrality);
@@ -1810,7 +1812,7 @@ struct HfTaskCharmPolarisation {
18101812 }
18111813 }
18121814 } else if constexpr (CosThetaStarType == charm_polarisation::CosThetaStarType::EP) { // EP
1813- // Dstar or Lc->pKpi
1815+ // Dstar or Lc->pKpi
18141816 if (origin == RecoDecay::OriginType::Prompt) {
18151817 if (!isPartRecoDstar) {
18161818 registry.fill (HIST (" hGenPromptEP" ), ptCharmHad, numPvContributors, std::abs (rapCharmHad), cosThetaStar, areDausInAcc, resoChannelLc, charge, centrality);
0 commit comments