HADRONS%2B%2B/ME_Library/B_K_Semileptonic.H

0001 #ifndef HADRONS_ME_Library_B_K_Semileptonic_H
0002 #define HADRONS_ME_Library_B_K_Semileptonic_H
0003
0004 #include "HADRONS++/ME_Library/HD_ME_Base.H"
0005 #include "ATOOLS/Phys/Flavour.H"
0006 #include "ATOOLS/Math/Vector.H"
0007 #include "ATOOLS/Math/MathTools.H"
0008 #include "ATOOLS/Math/Histogram.H"
0009 #include "ATOOLS/Org/Message.H"
0010 #include "HADRONS++/Main/Tools.H"
0011
0012 namespace HADRONS {
0013
0014   class B_K_Semileptonic : public HD_ME_Base {
0015     // indices: 0=B, 1=K, 2=lepton-, 3=lepton+
0016     double m_global;
0017     bool m_LD;
0018     Complex m_cR_T1, m_cR_T2, m_cL_T1, m_cL_T2;
0019     double m_C1, m_C2, m_C3, m_C4, m_C5, m_C6, m_C7eff, m_C9, m_C10;
0020     double m_mc, m_mb, m_ms, m_Vtb, m_Vub;
0021     public:
0022     B_K_Semileptonic(const ATOOLS::Flavour_Vector& flavs,
0023                      const std::vector<int>& decayindices,
0024                      const std::string& name):
0025       HD_ME_Base(flavs,decayindices,name) {};
0026     void Calculate(const ATOOLS::Vec4D_Vector& momenta, bool anti=false);
0027     void SetModelParameters( GeneralModel _md );
0028     Complex C9sehgal(double sHat);
0029     Complex sehgalld(double sHat);
0030     Complex g(double shat);
0031     Complex g0(double shat);
0032     Complex gc(double shat);
0033     double Theta(double x);
0034     Complex gSD(double mhat, double shat);
0035   };
0036 }
0037
0038 #endif