message: PWGCF: AliFemtoModelCorrFctnDEtaDPhiStar - Cleaned implementation of LCMS q calculation
diff: @@ -120,28 +120,21 @@ Qcms(const AliFemtoLorentzVector &p1, const AliFemtoLorentzVector &p2)
  , d = p1 - p2
  ;
  - Double_t dx = d.x()
  - , dy = d.y()
  - , dz = d.z();
  -
  - Double_t xt = p.x()
  - , yt = p.y()
  - , k1 = p.Perp()
  - , k2 = dx*xt + dy*yt;
  + Double_t k1 = p.Perp(),
  + k2 = d.x()*p.x() + d.y()*p.y();
  // relative momentum out component in lab frame
  Double_t qout = (k1 == 0) ? 0.0 : k2/k1;
  -
  // relative momentum side component in lab frame
  Double_t qside = (k1 == 0) ? 0.0 : 2.0 * (p2.x()*p1.y() - p1.x()*p2.y())/k1;
  // relative momentum component in lab frame
  - double beta = p.z()/p.t();
  - double gamma = 1.0 / TMath::Sqrt((1.0-beta)*(1.0+beta));
  + Double_t beta = p.z()/p.t(),
  + gamma = 1.0 / TMath::Sqrt((1.0-beta)*(1.0+beta));
  - double qlong = gamma * (dz - beta*d.t());
  + Double_t qlong = gamma * (d.z() - beta*d.t());
  // double qlong = (p.t()*d.z() - p.z()*d.t()) / TMath::Sqrt(p.t()*p.t() - p.z()*p.z());
  