Bsmumu.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 SusyFit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef BSMUMU_H
#define BSMUMU_H
 
class StandardModel;
class EvolBsmm;
#include <memory>
#include "ThObservable.h"
#include "QCD.h"
#include "OrderScheme.h"
#include "DGamma.h"
 
class Bsmumu : public ThObservable {
public:   
    Bsmumu(const StandardModel& SM_i, int obsFlag, QCD::lepton lep_i=QCD::MU);
    
    double computeThValue();
    double computeAmumu(orders order);
    double computeSmumu(orders order);
    
    
protected:
    
    void computeAmpSq(orders order, orders_qed order_qed, double mu);
    void computeObs(orders order, orders_qed order_qed);
    
private:
    
    QCD::lepton lep; 
    double beta;
    double mBs;
    double mW;
    double mlep;
    double mb;
    double ms;
    double chiral;
    double absP;
    double argP;
    double absS;
    double argS;
    double ampSq;
    double Amumu;
    double Smumu;
    double phiNP;
    double timeInt;
//    double convertFromSingletoDoubleGF;
    double coupling;
    int obs;
    int leptonindex;
    gslpp::complex C_10;
    gslpp::complex C_10p;
    gslpp::complex C_S;
    gslpp::complex C_Sp;
    gslpp::complex C_P;
    gslpp::complex C_Pp;
    
    std::unique_ptr<EvolBsmm> evolbsmm;
    gslpp::vector<gslpp::complex> ** allcoeff;
 
 
};
 
#endif  /* BSMUMU_H */