marley_kinematics.hh

//
// This file is part of MARLEY (Model of Argon Reaction Low Energy Yields)
//
// MARLEY is free software: you can redistribute it and/or modify it under the
// terms of version 3 of the GNU General Public License as published by the
// Free Software Foundation.
//
// For the full text of the license please see COPYING or
// visit http://opensource.org/licenses/GPL-3.0
//
// Please respect the MCnet academic usage guidelines. See GUIDELINES
// or visit https://www.montecarlonet.org/GUIDELINES for details.
 
#pragma once
#include <cmath>
 
#include "marley/Error.hh"
 
namespace marley {
  class Particle;
}
 
namespace marley_kinematics {
 
  // Lorentz boost a particle, replacing its energy and momentum with the
  // boosted versions
  void lorentz_boost(double beta_x, double beta_y, double beta_z,
    marley::Particle& particle_to_boost);
 
  void two_body_decay(const marley::Particle& initial_particle,
    marley::Particle& first_product, marley::Particle& second_product,
    double cos_theta_first, double phi_first);
 
  // Rotates a particle's 3-momentum so that it points in the (x, y, z)
  // direction
  void rotate_momentum_vector(double x, double y, double z,
    marley::Particle& particle_to_rotate);
 
  // Gets the square of the total energy in the center of momentum frame for
  // two particles
  double get_mandelstam_s(const marley::Particle& p1,
    const marley::Particle& p2);
 
  // Boost two particles into their mutual center-of-momentum frame
  void boost_to_cm_frame(marley::Particle& p1, marley::Particle& p2);
}