//
Double_t beta = TMath::TanH(fDyBoost);
- Double_t gamma = 1./TMath::Sqrt(1.-beta*beta);
+ Double_t gamma = 1./TMath::Sqrt((1.-beta)*(1.+beta));
Double_t gb = gamma * beta;
// printf("\n Boosting particles to lab frame %f %f %f", fDyBoost, beta, gamma);