YFS/Main/Coulomb.H

0001 #ifndef YFS_Main_Coulomb_H
0002 #define YFS_Main_Coulomb_H
0003
0004 #include "YFS/Main/YFS_Base.H"
0005
0006 namespace YFS{
0007
0008   class Coulomb : public YFS_Base
0009   {
0010   public:
0011     Coulomb();
0012     ~Coulomb();
0013
0014     ATOOLS::Vec4D p1,p2; // final state momentum eg ttbar or W+W-
0015     double m_weight, m_p, m_s1, m_s2, m_MW, m_GW, m_cms, m_threshold;
0016     int m_subtract;
0017     void Calculate(const ATOOLS::Vec4D p1, const ATOOLS::Vec4D p2);
0018     inline double GetWeight() { return m_weight; }
0019     inline void SetAlphaQED(double alp) { m_alpha = alp; }
0020     void Subtract();
0021   };
0022 }
0023
0024 #endif