[PWGDQ] Update VarManager.h to correct the cos theta definition in PP frame

KangkanGoswami · web-flow · commit 1ea95779d278 · 2025-06-03T11:00:31.000+05:30
Update VarManager.h to correct the cos theta definition in the production plane
diff --git a/PWGDQ/Core/VarManager.h b/PWGDQ/Core/VarManager.h
@@ -2910,7 +2910,7 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
       ROOT::Math::XYZVector yaxis_PP{(v12.Vect()).Unit()};
       ROOT::Math::XYZVector xaxis_PP{(yaxis_PP.Cross(zaxis_PP)).Unit()};
       if (fgUsedVars[kCosThetaPP]) {
-        values[kCosThetaPP] = zaxis_PP.Dot(v_CM);
+        values[kCosThetaPP] = zaxis_PP.Dot(v_CM)/std::sqrt(zaxis_PP.Mag2());
       }
       if (fgUsedVars[kPhiPP]) {
         values[kPhiPP] = TMath::ATan2(yaxis_PP.Dot(v_CM), xaxis_PP.Dot(v_CM));

Original file line number	Diff line number	Diff line change
`@@ -2910,7 +2910,7 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)`
`2910`	`2910`	`ROOT::Math::XYZVector yaxis_PP{(v12.Vect()).Unit()};`
`2911`	`2911`	`ROOT::Math::XYZVector xaxis_PP{(yaxis_PP.Cross(zaxis_PP)).Unit()};`
`2912`	`2912`	`if (fgUsedVars[kCosThetaPP]) {`
`2913`		`- values[kCosThetaPP] = zaxis_PP.Dot(v_CM);`
	`2913`	`+ values[kCosThetaPP] = zaxis_PP.Dot(v_CM)/std::sqrt(zaxis_PP.Mag2());`
`2914`	`2914`	`}`
`2915`	`2915`	`if (fgUsedVars[kPhiPP]) {`
`2916`	`2916`	`values[kPhiPP] = TMath::ATan2(yaxis_PP.Dot(v_CM), xaxis_PP.Dot(v_CM));`