forW-MEc.h

#ifndef _forW_MEc_h_included_
#define _forW_MEc_h_included_
#include <complex>
using std::complex;
 
namespace Photospp {
 
  class PhotosMEforW {
  public:
    static void   PHOBWnlo(double* WT);
 
  private:
    static double WDecayEikonalSqrKS_1ph(double p3[4], double p1[4], double p2[4], double k[4]);
    static double WDecayBornAmpSqrKS_1ph(double p3[4], double p1[4], double p2[4]);
    static double WDecayAmplitudeSqrKS_1ph(double p3[4], double p1[4], double p2[4], double k[4]);
    static double SANC_WT(double PW[4], double PNE[4], double PMU[4], double PPHOT[4], double B_PW[4], double B_PNE[4],
                          double B_PMU[4]);
    static void   SANC_INIT1(double QB0, double QF20, double MF10, double MF20, double MB0);
    static void   SANC_INIT(double ALPHA, int PHLUN);
 
  private:
    static complex<double> InProd_zero(double p1[4], int l1, double p2[4], int l2);
    static double          InSqrt(double p[4], double q[4]);
    static complex<double> InProd_mass(double p1[4], double m1, int l1, double p2[4], double m2, int l2);
    static complex<double> BsFactor(int s, double k[4], double p[4], double m);
    static complex<double> WDecayEikonalKS_1ph(double p3[4], double p1[4], double p2[4], double k[4], int s);
    static complex<double> SoftFactor(int s, double k[4], double p1[4], double m1, double p2[4], double m2, double Gmass2);
    static complex<double> TrMatrix_zero(double p1[4], double m1, int l1, double k[4], int s, double p2[4], double m2, int l2);
    static complex<double> TrMatrix_mass(double p1[4], double m1, int l1, double k[4], double m, int s, double p2[4], double m2,
                                         int l2);
    static complex<double> WDecayBornAmpKS_1ph(double p3[4], int l3, double p1[4], int l1, double p2[4], int l2);
    static complex<double> WDecayAmplitudeKS_1ph(double p3[4], int l3, double p1[4], int l1, double p2[4], int l2, double k[4], int s);
 
  private:
    // COMMON /Kleiss_Stirling/spV,bet
    static double spV[4], bet[4];
 
    // COMMON /mc_parameters/pi,sw,cw,alphaI,qb,mb,mf1,mf2,qf1,qf2,vf,af,mcLUN
    static double pi, sw, cw, alphaI, qb, mb, mf1, mf2, qf1, qf2, vf, af, mcLUN;
  };
 
} // namespace Photospp
#endif