Boost the space part of phase to PRF

Author: alibuild

PWGCF/FemtoUniverse/Core/FemtoUniverseMath.h

@@ -269,8 +269,9 @@ class FemtoUniverseMath
       const double mRS = (y1_lcms - y2_lcms) / 0.197327;
       const double mRL = (z1_lcms - z2_lcms) / 0.197327;
       const double mDT = (t1_lcms - t2_lcms) / 0.197327;
-
-      const double phase_inv = fDKOutPRF * mROutSmeared + fDKSidePRF * mRSideSmeared + fDKLongPRF * mRLongSmeared;
+      const double rOutPRF  = gammaOut * (mRO - betaOut * mDT);
+      
+      const double phase_inv = fDKOutPRF * mROut + fDKSidePRF * mRS + fDKLongPRF * mRL;
       double quantumweight = 1.0 + TMath::Cos(2.0*phase_inv);
 
       vect.push_back(quantumweight);