#include "BEAM/Main/Beam_Base.H" #include "ATOOLS/Org/Message.H" #include "ATOOLS/Org/My_MPI.H" using namespace BEAM; Beam_Base::Beam_Base(std::string _type,const ATOOLS::Flavour _beam, const double _energy,const double _polarisation, const int _dir) : m_type(_type), m_beam(_beam), m_bunch(m_beam), m_dir(_dir), m_energy(_energy), m_polarisation(_polarisation), m_x(1.), m_Q2(0.), m_weight(1.) { Init(); } bool Beam_Base::Init(int mode) { double disc = mode?1.0:1.0-ATOOLS::sqr(m_beam.Mass()/m_energy); if (disc<0) { msg_Error()<<"Error in Beam_Base :"<