ZFitterWrapper.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */
 
/*
 *  Documents of ZFITTER:
 *    hep-ph/9412201
 *    hep-ph/9908433 for v6.21
 *    hep-ph/0507146 for v6.42
 *
 *  Constants defined in the ZFITTER Fortran codes:
 *   - G_F^mu and alpha(0) are set in CONST1() in zfbib6_40.f and EWINIT()
 *     in zfbib6_40.f (and other places?). 
 *   - All lepton and quark masses (u,d,s,c,b) are set in CONST1().
 *   - Since it is hard to access some of these constants from the current
 *     class, the constants G_F^mu and alpha(0), which are necessary for
 *     the computations of the oblique parameters, are also set here in
 *     GF() and alpha().
 *
 *  Usage:
 *   1. call a constructor
 *        ZFitter()
 *   2. set flags (if necessary)
 *        setAllFlags() or setFlag()
 *   3. set cuts (if necessary)
 *   4. compute observables
 *
 *  Important flags in ZFITTER:
 *       AMT4 = 6 : with complete two-loop corrections to Mw and fermionic
 *                  two-loop corrections to sin^2(theta_eff^lept)
 *       ALEM = 2 : DAL5H is supplied by the user as input. 
 *       FOT2 = 3 : controls Initial-state radiation
 *       FINR = 0 : controls Final-state radiation 
 *       INTF = 0 : ISR/FSR interference is ignored. 
 *                  (INTF=1: with O(alpha), INTF=2: with the exponentiation)
 *       ISPP = 2 : with ISR pairs (ISPP=0: without the pairs)
 *       ICOL = -1: controls cuts
 *
 *  Fitted parameters in ZFITTER:
 *    ZMASS: the pole mass of Z in GeV
 *    TMASS: the pole mass of the top quark in GeV [10-400]
 *    HMASS: the pole mass of the Higgs in GeV [10-1000]
 *    ALFAS: the strong coupling constant at Mz
 *    DAL5H: the 5-quark flavor hadronic vacuum polarization at Mz
 *    UMASS: the constituent mass of the up quark (fixed)
 *    DMASS: the constituent mass of the down quark (fixed)
 *    V_TB: the CKM matrix elements V_{tb} (fixed)
 *
 *  INDF = Fermion index:
 *    INDF==0 neutrinos
 *    INDF==1 electron
 *    INDF==2 mu
 *    INDF==3 tau
 *    INDF==4 u
 *    INDF==5 d
 *    INDF==6 c
 *    INDF==7 s
 *    INDF==8 t
 *    INDF==9 b
 *    INDF==10 hadronic
 *    INDF==11 total
 * 
 */
 
#ifndef ZFITTERWRAPPER_H
#define ZFITTERWRAPPER_H
 
#include <iostream>
#include <cstring>
#include <gslpp.h>
 
 
#include <StandardModel.h>
 
using namespace gslpp;
 
 
class ZFitterWrapper {
public:
 
    ZFitterWrapper(const StandardModel& mySM);
 
 
    // Flags
    
    void setFlag(const std::string CHFLAG, const int IVALUE) const;
 
    void setAllFlags(const int flags[46], const int flagPrint) const;    
    
    void FlagInfo() const;
 
    
    // Cuts
    
    void setCuts(const int INDF, const int ICUT, const double ACOL,
                 const double EMIN, const double S_PR, const double ANG0,
                 const double ANG1, const double SIPP) const;
 
    void setAllCuts(const int ICUT[12], const double ACOL[12],
                    const double EMIN[12], const double S_PR[12],
                    const double ANG0[12], const double ANG1[12],
                    const double SPP[12], const int flagPrint) const;
 
    void setSprimeCut(const int INDF, const double s) const;
    
    void CutInfo() const;
    
    
    // Effective couplings
    
    complex rhoZ_f(const int INDF) const;
    
    complex kappaZ_f(const int INDF) const;
 
    complex gZ_f(const int INDF) const;
 
    
    // LR-asymmetry parameter
    
    double Af(const int INDF) const;
    
    
    // Z-pole observables
    
    double alphaMZ() const;
    
    double Mw() const;
    
    double Gamma_W() const;
 
    double sw2() const;
 
    double s2teff_f(const int INDF) const;
    
    double Gamma_f(const int INDF) const;
 
    double Gamma_inv() const;
 
    double Gamma_had() const;
 
    double Gamma_Z() const;
    
    
    
    void calcXS_AFB(const int INDF, const double SQRS, double *XS, double *AFB) const;
 
    void calcDXS(const int INDF, const double SQRS, const double CSA, double *DXS) const;
 
    void calcTauPol(const double SQRS, double *TAUPOL, double *TAUAFB) const;
 
    void calcALR(const int INDF, const double SQRS, const double POL,
                 double *XSPL, double *XSMI) const;
 
    void calcAPV(double *C1U, double *C1D, double *C2U, double *C2D) const;
 
    void calcXS(const int INDF, const double SQRS, const double GAMZ0,
                const double GAMEE, const double GAMFF, double *XS) const;
 
    void calcXS_AFB_2(const int INDF, const double SQRS, const double GAMZ0,
                      const int MODE, const double GVE, const double XE,
                      const double GVF, const double XF, double *XS, double *AFB) const;
 
    void calcXS_AFB_3(const int INDF, const double SQRS, const double GAMZ0,
                      const int MODE, const double GV2, const double X2,
                      double *XS, double *AFB) const;
 
    void calcXS_AFB_4(const int INDF, const double SQRS, const double GAMZ0,
                      const double PFOUR, const double PVAE2, const double PVAF2,
                      double *XS, double *AFB) const;
 
    void calcTauPol_2(const double SQRS, const double GAMZ0, const int MODE,
                      const double GVE, const double XE, const double GVF,
                      const double XF, double *TAUPOL, double *TAUAFB) const;
 
    
 
    void printConstants() const;
 
    void printInputs() const;
 
    void printIntermediateResults() const;
 
    std::string convertINDF(const int INDF) const;
    
    void test(const int IMISC) const;
    
    const StandardModel& getModel() const {
        return SM;
    }
 
    
    
private:
    const StandardModel& SM;
    
    mutable bool IsFlagChanged;
 
    // Caches
    mutable double MzCache, MtCache, MhCache, AlsMzCache, DAle5MzCache;
 
    // Constants (set in the constructor)
    double V_TB;
    double UMASS; // the constituent mass of the up quark
    double DMASS; // the constituent mass of the down quark 
 
    
    
    void init(const int IPRINT) const;
 
    void calcCommonBlocks() const;
 
        
    // get variables from Common blocks of ZFITTER
 
    double getCommonALPHST() const;
 
    double getCommonSIN2TW() const;
 
    double getCommonS2TEFF(const int INDF) const;
 
    double getCommonALLCH(const int INDF) const;
 
    double getCommonALLMS(const int INDF) const;
    
    double getCommonWIDTHS(const int INDF) const;
    
    double getCommonPARTZ(const int INDF) const;
   
    double getCommonPARTW(const int i) const;   
    
 
    // get variables from Common blocks of DIZET
 
    double getCommonARROFZ(const int INDF) const;
 
    double getCommonARKAFZ(const int INDF) const;
 
    double getCommonARVEFZ(const int INDF) const;
 
    double getCommonARSEFZ(const int INDF) const;
 
    double getCommonAROTFZ(const int INDF) const;
 
    double getCommonAIROFZ(const int INDF) const;
 
    double getCommonAIKAFZ(const int INDF) const;
 
    double getCommonAIVEFZ(const int INDF) const;    
    
    
};
    
#endif  /* ZFITTERWRAPPER_H */