TopQuarkObservables.h

Go to the documentation of this file.

/* 
 * File:   TopQuarkObservables.h
 * Author: silvest
 *
 * Created on 21 settembre 2023, 16.45
 */
 
#ifndef TOPQUARKOBSERVABLES_H
#define TOPQUARKOBSERVABLES_H
 
#include "NPSMEFTd6General.h"
#include "ThObservable.h"
 
#include "PVfunctions.h"
 
 
class TopQuarkObservables {
public:
 
    
    
    TopQuarkObservables(const NPSMEFTd6General& NP_i);
//    TopQuarkObservables(const TopQuarkObservables& orig);
    virtual ~TopQuarkObservables(){};
    
 
    const NPSMEFTd6General& GetmyNPSMEFTd6General() const {
        return myNPSMEFTd6General;
    }
 
    
    /*
    inline double getSMEFTCoeffEW(const std::string name) const
    {
        return GetmyNPSMEFTd6General().getSMEFTCoeffEW(name);
    }
    
    inline double getSMEFTCoeffEW(const std::string name, int i, int j) const
    {
        return GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j);
    }
    
    inline double getSMEFTCoeffEW(const std::string name, int i, int j, int k, int l) const
    {
        return GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l);
    }
    */
    
protected:
        
 
 
private:
    
    const NPSMEFTd6General& myNPSMEFTd6General;
  
};
 
    
 
 
 
    class FB_asymmetry_Tevatron_tt_diff_mtt_LO: public ThObservable {
    public:
 
    FB_asymmetry_Tevatron_tt_diff_mtt_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
        
//        const PVfunctions PV;
    };
 
    
    
    
    class sigma_tt_diff_mtt_LO_CMS_181106625: public ThObservable {
    public:
 
    sigma_tt_diff_mtt_LO_CMS_181106625(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
 
 
    
    class sigma_tt_diff_mtt_CMS_LO: public ThObservable {
    public:
 
    sigma_tt_diff_mtt_CMS_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    class sigma_norm_tt_diff_mtt_ATLAS_LO: public ThObservable {
    public:
 
    sigma_norm_tt_diff_mtt_ATLAS_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    class sigma_tt_13_LO: public ThObservable {
    public:
 
    sigma_tt_13_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    class R_tt_8_o_7_LO: public ThObservable {
    public:
 
    R_tt_8_o_7_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    class R_tt_13_o_8_LO: public ThObservable {
    public:
 
    R_tt_13_o_8_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    class charge_asymmetry_tt_diff_mtt_LO: public ThObservable {
    public:
 
    charge_asymmetry_tt_diff_mtt_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    class sigma_ttz_diff_LO_CMS_190711270: public ThObservable {
    public:
 
    sigma_ttz_diff_LO_CMS_190711270(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    
    
    
    
    class sigma_ttz_diff_LO_ATLAS_210312603: public ThObservable {
    public:
 
    sigma_ttz_diff_LO_ATLAS_210312603(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    class sigma_ttz_diff_LO_ATLAS_231204450: public ThObservable {
    public:
 
    sigma_ttz_diff_LO_ATLAS_231204450(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    
    class sigma_tta_diff_LO_ATLAS_emu_200706946: public ThObservable {
    public:
 
    sigma_tta_diff_LO_ATLAS_emu_200706946(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
 
    
    
    
    
    
    class sigma_tta_diff_LO_CMS_semileptonic_210701508: public ThObservable {
    public:
 
    sigma_tta_diff_LO_CMS_semileptonic_210701508(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    class sigma_tta_diff_LO_CMS_dilepton_220107301: public ThObservable {
    public:
 
    sigma_tta_diff_LO_CMS_dilepton_220107301(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
     
        double b_min = 0.;
        double b_max = 0.;
        
        
        const TopQuarkObservables mytopobs;
 
    };
    
 
    
    
    
    
    
    
        
    class sigma_tb_8_LO: public ThObservable {
    public:
 
    sigma_tb_8_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    class sigma_tb_13_LO: public ThObservable {
    public:
 
    sigma_tb_13_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    class sigma_tq_7_LO: public ThObservable {
    public:
 
    sigma_tq_7_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    class sigma_tq_8_LO: public ThObservable {
    public:
 
    sigma_tq_8_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    class sigma_tq_13_LO: public ThObservable {
    public:
 
    sigma_tq_13_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    class sigma_tq_top_13_LO: public ThObservable {
    public:
 
    sigma_tq_top_13_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    
    
    
    
    
    class sigma_tq_antitop_13_LO: public ThObservable {
    public:
 
    sigma_tq_antitop_13_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
 
    
    
    
    
    
    class sigma_taq_LO_CMS: public ThObservable {
    public:
 
    sigma_taq_LO_CMS(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    class sigma_taq_LO_ATLAS: public ThObservable {
    public:
 
    sigma_taq_LO_ATLAS(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    class sigma_tzq_LO: public ThObservable {
    public:
 
    sigma_tzq_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    class sigma_tw_7_LO: public ThObservable {
    public:
 
    sigma_tw_7_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    
    
    
    class sigma_tw_8_LO: public ThObservable {
    public:
 
    sigma_tw_8_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    
    class sigma_tw_13_LO: public ThObservable {
    public:
 
    sigma_tw_13_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
 
    
    
    
    
    class sigma_ttw_LO: public ThObservable {
    public:
 
    sigma_ttw_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
 
    class R_ttw_LO: public ThObservable {
    public:
 
    R_ttw_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    
    
    
    
    class F0_LO: public ThObservable {
    public:
 
    F0_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
 
    const TopQuarkObservables mytopobs;
 
    };
 
 
    class FL_LO: public ThObservable {
    public:
 
    FL_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
    
    
    private:
 
    const TopQuarkObservables mytopobs;
    
    
    
 
    };
 
 
    
    
    
    
    
    class FR_LO: public ThObservable {
    public:
 
    FR_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
    
    
    private:
 
    const TopQuarkObservables mytopobs;
    
 
    };
    
    
    
    
    
    
    class sigma_tttt_13_LO: public ThObservable {
    public:
 
    sigma_tttt_13_LO(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    
    class sigma_ttbb_13_LO_dilepton: public ThObservable {
    public:
 
    sigma_ttbb_13_LO_dilepton(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
    class sigma_ttbb_13_LO_lepjet: public ThObservable {
    public:
 
    sigma_ttbb_13_LO_lepjet(const StandardModel& SM_i);
 
    double computeThValue();
    
    
    
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double ewgc(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double ewgc(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
 
    private:
        
        const TopQuarkObservables mytopobs;
 
    };
    
    
    
    
    
 
 
    class Test_direct: public ThObservable {
    public:
 
    Test_direct(const StandardModel& SM_i);
 
    double computeThValue();
 
    //We need to multiply by the square of the scale used to obtain the parametrisation (1000 GeV)
    inline double getSMEFTCoeffEW(const std::string name) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
    }
    
    inline double getSMEFTCoeffEW(const std::string name, int i, int j) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
    }
    
    inline double getSMEFTCoeffEW(const std::string name, int i, int j, int k, int l) const
    {
        return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
    }
    
    
    private:
 
    const TopQuarkObservables mytopobs;
 
    
 
    };
    
    
 
 
 
#endif /* TOPQUARKOBSERVABLES_H */