// Return a modified copy of the four-vector FourVector FourVector::boost_z(const double velocity) { // natural units: c = 1 double gamma = 1.0 / sqrt(1.0 - velocity * velocity); double t2 = gamma * (t - velocity * z); double z2 = gamma * (z - velocity * t); return FourVector(t2, x, y, z2); // return a FourVector }
void Decay_gammagamma::do_decay() { do_decay(FourVector(Model.higgs.m(), 0., 0., 0.)); }