PhotosMEforZ Class Reference

Static Public Member Functions
static double	phwtnlo ()
static void	set_VakPol (double(*fun)(double[4], double[4], double[4], double[4], double[4], int, int))

Static Private Member Functions
static void	GIVIZO (int idferm, int ihelic, double *sizo3, double *charge, int *kolor)
static double	PHBORNM (double svar, double costhe, double T3e, double qe, double T3f, double qf, int Ncf)
static double	AFBCALC (double SVAR, int IDEE, int IDFF)
static int	GETIDEE (int IDE)
static double	PHASYZ (double SVAR, int IDE, int IDF)
static double	Zphwtnlo (double svar, double xk, int IDHEP3, int IREP, double qp[4], double qm[4], double ph[4], double pp[4], double pm[4], double COSTHG, double BETA, double th1, int IDE, int IDF)
static double	VakPol (double qp[4], double qm[4], double ph[4], double pp[4], double pm[4], int IDE, int IDF)
static double	default_VakPol (double qp[4], double qm[4], double ph[4], double pp[4], double pm[4], int IDE, int IDF)

Static Private Attributes
static double(*	currentVakPol )(double[4], double[4], double[4], double[4], double[4], int, int)

Detailed Description

Definition at line 6 of file forZ-MEc.h.

Member Function Documentation

PHBORNM()

double PHBORNM ( double  svar, double  costhe, double  T3e, double  qe, double  T3f, double  qf, int  NCf  )

staticprivate

                                                                  //

This routine provides unsophisticated Born differential cross section // at the crude x-section level, with Z and gamma s-chanel exchange. //

Definition at line 62 of file forZ-MEc.cc.

  {
 
    double   s, Sw2, MZ, MZ2, GammZ, AlfInv, GFermi; // t,MW,MW2,
    double   Ve, Ae, thresh;        //  sum,deno,
    double   xe, yf, xf, ye, ff0, ff1, amx2, amfin, Vf, Af;
    double   ReChiZ, SqChiZ, RaZ;   //,RaW,ReChiW,SqChiW;
    double   Born;                  //, BornS;
    //  int  KeyZet,HadMin,KFbeam;
    //  int  i,ke,KFfin,kf,IsGenerated,iKF;
    int  KeyWidFix;
 
    AlfInv = 137.0359895;
    GFermi = 1.16639e-5;
 
    //--------------------------------------------------------------------
    s = svar;
    //------------------------------
    //     EW paratemetrs taken from BornV
    MZ = 91.187;
    GammZ = 2.50072032;
    Sw2 = .22276773;
    //------------------------------
    // Z and gamma couplings to beams (electrons)
    // Z and gamma couplings to final fermions
    // Loop over all flavours defined in m_xpar(400+i)
 
 
    //------ incoming fermion
    Ve =  2 * T3e - 4 * qe * Sw2;
    Ae =  2 * T3e;
    //------ final fermion couplings
    amfin = 0.000511; //  m_xpar(kf+6)
    Vf =  2 * T3f - 4 * qf * Sw2;
    Af =  2 * T3f;
    if (fabs(costhe) > 1.0) {
      cout << "+++++STOP in PHBORN: costhe>0 =" << costhe << endl;
      exit(-1);
    }
    MZ2  = MZ * MZ;
    RaZ  = (GFermi * MZ2 * AlfInv) / (sqrt(2.0) * 8.0 * PI); //
    RaZ  = 1 / (16.0 * Sw2 * (1.0 - Sw2));
    KeyWidFix = 1;       // fixed width
    KeyWidFix = 0;       // variable width
    if (KeyWidFix == 0) {
      ReChiZ = (s - MZ2) * s / ((s - MZ2) * (s - MZ2) + (GammZ * s / MZ) * (GammZ * s / MZ)) * RaZ; // variable width
      SqChiZ =      s * s / ((s - MZ2) * (s - MZ2) + (GammZ * s / MZ) * (GammZ * s / MZ)) * RaZ * RaZ; // variable width
    } else {
      ReChiZ = (s - MZ2) * s / ((s - MZ2) * (s - MZ2) + (GammZ * MZ) * (GammZ * MZ)) * RaZ; // fixed width
      SqChiZ =      s * s / ((s - MZ2) * (s - MZ2) + (GammZ * MZ) * (GammZ * MZ)) * RaZ * RaZ; // fixed width
    }
    xe = Ve * Ve + Ae * Ae;
    xf = Vf * Vf + Af * Af;
    ye = 2 * Ve * Ae;
    yf = 2 * Vf * Af;
    ff0 = qe * qe * qf * qf + 2 * ReChiZ * qe * qf * Ve * Vf + SqChiZ * xe * xf;
    ff1 =             +2 * ReChiZ * qe * qf * Ae * Af + SqChiZ * ye * yf;
    Born    = (1.0 + costhe * costhe) * ff0 + 2.0 * costhe * ff1;
    // Colour factor
    Born = NCf * Born;
    // Crude method of correcting threshold, cos(theta) depencence incorrect!!!
    if (svar <  4.0 * amfin * amfin) {
      thresh = 0.0;
    } else if (svar < 16.0 * amfin * amfin) {
      amx2 = 4.0 * amfin * amfin / svar;
      thresh = sqrt(1.0 - amx2) * (1.0 + amx2 / 2.0);
    } else {
      thresh = 1.0;
    }
 
    Born = Born * thresh;
    return Born;
  }

Zphwtnlo()

double Zphwtnlo ( double  svar, double  xk, int  IDHEP3, int  IREP, double  qp[4], double  qm[4], double  ph[4], double  pp[4], double  pm[4], double  COSTHG, double  BETA, double  th1, int  IDE, int  IDF  )

staticprivate

waga=waga*wagan2 waga=waga*(1-delta) /wt2 ! sprawdzone ze to jest =2/(1.D0+COSTHG*BETA)

% * svar/4./xkap*(1.D0-COSTHG*BETA)*sqrt(1.0-xk)

fprintf(PHLUN,"%f %f %f %f %f") 'pneutr= ',phomom.pneutr[0],phomom.pneutr[1],phomom.pneutr[2],phomom.pneutr[3],phomom.pneutr[4];

Definition at line 243 of file forZ-MEc.cc.

---

The documentation for this class was generated from the following files:

• generators/kkmc/photospp/forZ-MEc.h
• generators/kkmc/photospp/forZ-MEc.cc