THDMcache.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

#ifndef THDMCACHE_H
#define THDMCACHE_H

#include <cmath>
#include "PVfunctions.h"
//#include "THDM.h"

#include <stdexcept>
#include "gslpp.h"

class THDMcache {
    
public:

    THDMcache();
    
    static const int CacheSize = 5;

    int CacheCheck(const gslpp::complex cache[][CacheSize],
                   const int NumPar, const double params[]) const;

    void CacheShift(gslpp::complex cache[][CacheSize], const int NumPar,
                    const double params[], const gslpp::complex newResult) const; 

    const PVfunctions getPV() const {
        return PV;
    }

    gslpp::matrix<double> readTable(std::string filename, int rowN, int colN);

    double interpolate (gslpp::matrix<double> arrayTab, double x);

    double interpolate2D (gslpp::matrix<double> arrayTab, double x, double y);

    void read();

    gslpp::matrix<double> br_tt, br_bb, br_tautau, br_cc, br_mumu, br_ZZ, br_WW;

    gslpp::matrix<double> pc_ggF, pc_VBF, pc_WH, pc_ZH, pc_ttH;

    gslpp::matrix<double> GammaHtotSM;

    gslpp::matrix<double> cs_ggH, cs_ggH_tt, cs_ggH_bb, cs_ggA, cs_ggA_tt, cs_ggA_bb;

    gslpp::matrix<double> cs_bbFtoHP;

    gslpp::matrix<double> ATLAS_ggF_phi_gaga, ATLAS_ggF_phi_tautau, ATLAS_bbF_phi_tautau, ATLAS_ggF_A_hZ_tautauZ, ATLAS_ggF_A_hZ_bbZ, ATLAS_pp_phi_tt, ATLAS_ggF_H_WW, ATLAS_VBF_H_WW, ATLAS_ggF_H_hh;

    gslpp::matrix<double> CMS_pp_H_ZZ;

    gslpp::matrix<double> CMS_ggF_A_hZ_bbll, CMS_pp_H_hh_gagabb, CMS_pp_H_hh_bbbb, CMS_bbF_phi_bb, CMS_ggF_phi_tautau, CMS_bbF_phi_tautau, CMS_ggF_phi_gaga, CMS_ggF_H_hh_bbtautau, CMS_ggF_A_hZ_tautaull;

    gslpp::matrix<double> arraybsgamma;

    double ip_Br_HPtott(double mass);

    double ip_Br_HPtobb(double mass);

    double ip_Br_HPtotautau(double mass);

    double ip_Br_HPtocc(double mass);

    double ip_Br_HPtomumu(double mass);

    double ip_Br_HPtoZZ(double mass);

    double ip_Br_HPtoWW(double mass);

    double ip_pc_ggFtoHP(double mass);

    double ip_pc_VBFtoHP(double mass);

    double ip_pc_WHP_HP(double mass);

    double ip_pc_ZHP_HP(double mass);

    double ip_pc_ttFtoHP(double mass);

    double ip_GammaHPtotSM(double mass);

    double ip_cs_ggFtoHP(double mass);

    double ip_cs_ggHP_tt(double mass);

    double ip_cs_ggHP_bb(double mass);

    double ip_cs_ggA(double mass);

    double ip_cs_ggA_tt(double mass);

    double ip_cs_ggA_bb(double mass);

    double ip_cs_bbFtoHP(double mass);

    double ip_ex_ggF_phi_gaga_ATLAS(double mass);

    double ip_ex_ggF_phi_tautau_ATLAS(double mass);

    double ip_ex_bbF_phi_tautau_ATLAS(double mass);

    double ip_ex_ggF_A_hZ_tautauZ_ATLAS(double mass);

    double ip_ex_ggF_A_hZ_bbZ_ATLAS(double mass);

    double ip_ex_pp_phi_tt_ATLAS(double mass);

    double ip_ex_ggF_H_WW_ATLAS(double mass);

    double ip_ex_VBF_H_WW_ATLAS(double mass);

    double ip_ex_ggF_H_hh_ATLAS(double mass);

    double ip_ex_pp_H_ZZ_CMS(double mass);

    double ip_ex_ggF_A_hZ_bbll_CMS(double mass);

    double ip_ex_pp_phi_hh_gagabb_CMS(double mass);

    double ip_ex_pp_phi_hh_bbbb_CMS(double mass);

    double ip_ex_bbF_phi_bb_CMS(double mass);

    double ip_ex_ggF_phi_tautau_CMS(double mass);

    double ip_ex_bbF_phi_tautau_CMS(double mass);

    double ip_ex_ggF_phi_gaga_CMS(double mass);

    double ip_ex_ggF_H_hh_bbtautau_CMS(double mass);

    double ip_ex_ggF_A_hZ_tautaull_CMS(double mass);

    double ip_ex_bsgamma(double logtb, double logmHp);

    /*One-loop functions*/

    gslpp::complex B0_MZ2_0_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const;
    gslpp::complex B0_MZ2_0_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const;
    gslpp::complex B0_MZ2_0_MZ2_mHh2(const double MZ2, const double mHh2) const;
    gslpp::complex B0_MZ2_0_MZ2_mHl2(const double MZ2, const double mHl2) const;
    gslpp::complex B0_MZ2_MW2_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const;
    gslpp::complex B0_MZ2_MW2_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const;
    gslpp::complex B0_MZ2_MZ2_MZ2_mHh2(const double MZ2, const double mHh2) const;
    gslpp::complex B0_MZ2_MZ2_MZ2_mHl2(const double MZ2, const double mHl2) const;

    gslpp::complex B00_MZ2_0_mA2_mHp2(const double MZ2, const double mA2, const double mHp2) const;
    gslpp::complex B00_MZ2_0_mHh2_mA2(const double MZ2, const double mHh2, const double mA2) const;
    gslpp::complex B00_MZ2_0_mHh2_mHp2(const double MZ2, const double mHh2, const double mHp2) const;
    gslpp::complex B00_MZ2_0_mHl2_mA2(const double MZ2, const double mHl2, const double mA2) const;
    gslpp::complex B00_MZ2_0_mHl2_mHp2(const double MZ2, const double mHl2, const double mHp2) const;
    gslpp::complex B00_MZ2_0_mHp2_mHp2(const double MZ2, const double mHp2) const;
    gslpp::complex B00_MZ2_0_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const;
    gslpp::complex B00_MZ2_0_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const;
    gslpp::complex B00_MZ2_0_MZ2_mHh2(const double MZ2, const double mHh2) const;
    gslpp::complex B00_MZ2_0_MZ2_mHl2(const double MZ2, const double mHl2) const;
    gslpp::complex B00_MZ2_MW2_mA2_mHp2(const double MZ2, const double MW2, const double mA2, const double mHp2) const;
    gslpp::complex B00_MZ2_MW2_mHh2_mHp2(const double MZ2, const double MW2, const double mHh2, const double mHp2) const;
    gslpp::complex B00_MZ2_MW2_mHl2_mHp2(const double MZ2, const double MW2, const double mHl2, const double mHp2) const;
    gslpp::complex B00_MZ2_MW2_mHp2_mHp2(const double MZ2, const double MW2, const double mHp2) const;
    gslpp::complex B00_MZ2_MW2_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const;
    gslpp::complex B00_MZ2_MW2_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const;
    gslpp::complex B00_MZ2_MZ2_mHh2_mA2(const double MZ2, const double mHh2, const double mA2) const;
    gslpp::complex B00_MZ2_MZ2_mHl2_mA2(const double MZ2, const double mHl2, const double mA2) const;
    gslpp::complex B00_MZ2_MZ2_mHp2_mHp2(const double MZ2, const double mHp2) const;
    gslpp::complex B00_MZ2_MZ2_MZ2_mHh2(const double MZ2, const double mHh2) const;
    gslpp::complex B00_MZ2_MZ2_MZ2_mHl2(const double MZ2, const double mHl2) const;


private:

    const PVfunctions PV;
    //const THDM myTHDM;

    //Caches
    
    /*One-loop functions*/

    mutable gslpp::complex B0_MZ2_0_MW2_mHh2_cache[4][CacheSize];
    mutable gslpp::complex B0_MZ2_0_MW2_mHl2_cache[4][CacheSize];
    mutable gslpp::complex B0_MZ2_0_MZ2_mHh2_cache[3][CacheSize];
    mutable gslpp::complex B0_MZ2_0_MZ2_mHl2_cache[3][CacheSize];
    mutable gslpp::complex B0_MZ2_MW2_MW2_mHh2_cache[4][CacheSize];
    mutable gslpp::complex B0_MZ2_MW2_MW2_mHl2_cache[4][CacheSize];
    mutable gslpp::complex B0_MZ2_MZ2_MZ2_mHh2_cache[3][CacheSize];
    mutable gslpp::complex B0_MZ2_MZ2_MZ2_mHl2_cache[3][CacheSize];

    mutable gslpp::complex B00_MZ2_0_mA2_mHp2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_0_mHh2_mA2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_0_mHh2_mHp2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_0_mHl2_mA2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_0_mHl2_mHp2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_0_mHp2_mHp2_cache[3][CacheSize];
    mutable gslpp::complex B00_MZ2_0_MW2_mHh2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_0_MW2_mHl2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_0_MZ2_mHh2_cache[3][CacheSize];
    mutable gslpp::complex B00_MZ2_0_MZ2_mHl2_cache[3][CacheSize];
    mutable gslpp::complex B00_MZ2_MW2_mA2_mHp2_cache[5][CacheSize];
    mutable gslpp::complex B00_MZ2_MW2_mHh2_mHp2_cache[5][CacheSize];
    mutable gslpp::complex B00_MZ2_MW2_mHl2_mHp2_cache[5][CacheSize];
    mutable gslpp::complex B00_MZ2_MW2_mHp2_mHp2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_MW2_MW2_mHh2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_MW2_MW2_mHl2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_MZ2_mHh2_mA2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_MZ2_mHl2_mA2_cache[4][CacheSize];
    mutable gslpp::complex B00_MZ2_MZ2_mHp2_mHp2_cache[3][CacheSize];
    mutable gslpp::complex B00_MZ2_MZ2_MZ2_mHh2_cache[3][CacheSize];
    mutable gslpp::complex B00_MZ2_MZ2_MZ2_mHl2_cache[3][CacheSize];

};

#endif  /* THDMCACHE_H */