BDlnuDecay.h

Go to the documentation of this file.

#ifndef BDLNUDECAY_H
#define BDLNUDECAY_H
 
#include "DecayParent.h"
#include "ObsWilsonProxy.h"
#include "ObsQCDProxy.h"
#include "DefaultConfig.h"
 
struct BDlnuConfig : public DecayConfig {
    enum class B_Charge {B_0, B_PLUS};
    B_Charge charge {B_Charge::B_PLUS};
};
 
struct BDlnuDecayCache {
    double G_F;
    double m_e, m_tau;
    double m_B, m_D;
    double tau_B;
    double V11, rho_D2, Delta;
 
    complex_t C_V, C_S, C_T;
    bool C_V_flag, C_S_flag, C_T_flag;
 
    double r_D, r_e, r_tau, r_qp, r_qm;
    double w_e, w_tau;
    double BR_pref;
    double Gamma_p, Gamma_m;
};
 
class BDlnuDecay : public DecayParentConfigurable<BDlnuConfig> {
private:
    BDlnuDecayCache cache;
    BDlnuConfig cfg;
 
protected:
    // Utility
    void set_cfg_flags(BDlnuConfig::B_Charge charge);
    void fill_wilson_cache();
    void load_cfg_dep_params();
 
    double t(double w);
    double lambda_D(double w);
    double x_l(double rl, double w);
    double phi(double rl, double w);
    double w_max(double rD, double r_l);
 
    double V_1(double w);
    double S_1(double w);
 
    double H_V0 (double w);
    double H_Vt (double w);
    double H_S  (double w);
    double H_T  (double w);
 
    double F_V0_1  (double r_l, double w_m);
    double F_V0_2  (double r_l, double w_m);
    double F_Vt    (double r_l, double w_m);
    double F_S     (double r_l, double w_m);
    double F_T_1   (double r_l, double w_m);
    double F_T_2   (double r_l, double w_m);
    double G_V0_Vt (double r_l, double w_m);
    double G_V0_S  (double r_l, double w_m);
    double G_V0_T  (double r_l, double w_m);
    double G_Vt_S  (double r_l, double w_m);
    double G_Vt_T  (double r_l, double w_m);
    double G_S_T   (double r_l, double w_m);
 
    double gamma_p(double r_l, double w_m);
    double gamma_m(double r_l, double w_m);
 
    double BR();
    double A_FB();
    double R_D();
    double P_tau();
 
public:
    BDlnuDecay(QCDOrder order, double matching_scale, double hadronic_scale, ObservablePortsConfig& ports) : DecayParentConfigurable(DecayMapper::to_id(Decays::B__D_l_nu), matching_scale, hadronic_scale, order, ports) {
        this->w_config.groups = {GroupMapper::to_id(WGroup::CC_bc)};
        this->max_order = QCDOrder::LO;
    }
 
    void load_params() override;
    std::vector<ObservableValue> compute_observable(Observables obs) override;
    std::vector<ObservableValue> compute_observable(ObservableId obs) override;
    void set_config_spe(BDlnuConfig config) override {this->cfg = config;}
};
 
#endif // __BDLNUDECAY_H__