SimplePhaseSpace.h

// -*- C++ -*-
//
// SimplePhaseSpace.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef ThePEG_SimplePhaseSpace_H
#define ThePEG_SimplePhaseSpace_H
 
#include "ThePEG/Config/ThePEG.h"
 
#include "ThePEG/Vectors/LorentzRotation.h"
#include "ThePEG/Vectors/LorentzRotation.h"
#include "ThePEG/EventRecord/Particle.h"
#include "ThePEG/EventRecord/ParticleTraits.h"
#include "ThePEG/Repository/UseRandom.h"
#include "SimplePhaseSpace.xh"
#include <numeric>
 
namespace ThePEG {
 
namespace SimplePhaseSpace {
 
  template <typename PType>
  void CMS(Energy2 s, PType & p1, PType & p2);
 
  template <typename PType>
  void CMS(PType & p1, PType & p2, Energy2 s,
                  double cosTheta, double phi);
 
  template <typename PType>
  void CMS(PType & p1, PType & p2, Energy2 s, Energy2 t, double phi,
                  const PType & p0);
 
  template <typename PType>
  void CMS(PType & p1, PType & p2, Energy2 s);
 
  template <typename PPairType>
  void CMS(const PPairType & p, Energy2 s)
  {
    CMS(*p.first, *p.second, s);
  }
 
  template <typename PType>
  void CMS(PType & p1, PType & p2, PType & p3, Energy2 s,
                  double x1, double x3);
 
  template <typename PType>
  void CMS(PType & p1, PType & p2, PType & p3, Energy2 s,
                  double x1, double x3, double phii = 0.0,
                  double theta = 0.0, double phi = 0.0);
 
  Energy getMagnitude(Energy2 s, Energy m1, Energy m2);
 
//  static Energy checkMagnitude(Energy2 s, Energy m1, Energy m2);
 
  inline Momentum3 polar3Vector(Energy p, double costheta, double phi)
  {
    return Momentum3(p*sqrt(1.0 - sqr(costheta))*sin(phi),
                     p*sqrt(1.0 - sqr(costheta))*cos(phi),
                     p*costheta);
  }
 
  vector<LorentzMomentum>
  CMSn(Energy m0, const vector<Energy> & m);
 
  template <typename Container>
  void CMSn(Container & particles, Energy m0);
 
}
 
}
 
#ifndef ThePEG_TEMPLATES_IN_CC_FILE
#include "SimplePhaseSpace.tcc"
#endif
 
#endif /* ThePEG_SimplePhaseSpace_H */