EvtBSemiTauonicHelicityAmplitudeCalculator.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#pragma once
 
#include "EvtGenBase/EvtComplex.hh"
 
namespace Belle2 {
  class EvtBSemiTauonicHelicityAmplitudeCalculator {
 
  public:
 
    EvtBSemiTauonicHelicityAmplitudeCalculator();
 
    EvtBSemiTauonicHelicityAmplitudeCalculator(const double rho12, const double rhoA12, const double ffR11, const double ffR21,
                                               const double AS1, const double AR3,
                                               const double bottomMass, const double charmMass,
                                               const EvtComplex& CV1, const EvtComplex& CV2, const EvtComplex& CS1, const EvtComplex& CS2, const EvtComplex& CT,
                                               const double parentMass, const double DMass, const double DstarMass);
 
    virtual ~EvtBSemiTauonicHelicityAmplitudeCalculator() {};
 
    EvtComplex helAmp(double mtau, int tauhel, int Dhel, double w, double costau) const;
 
  public:
 
    EvtComplex helAmp(const EvtComplex& CV1, const EvtComplex& CV2, const EvtComplex& CS1, const EvtComplex& CS2, const EvtComplex& CT,
                      double mtau, int tauhel, int Dhel, double w, double costau) const;
 
    double Lep(const double mtau, int tauhel, int whel, double q2, double costau) const; // Vector
 
    double Lep(const double mtau, int tauhel, double q2, double costau) const; // Scalar
 
    double Lep(const double mtau, int tauhel, int whel1, int whel2, double q2, double costau) const; // Tensor
 
    double HadV1(int Dhel, int whel, double w) const; // V-A
 
    double HadV2(int Dhel, int whel, double w) const; // V+A
 
    double HadS1(int Dhel, double w) const; // S+P
 
    double HadS2(int Dhel, double w) const; // S-P
 
    double HadT(int Dhel, int whel1, int whel2, double w) const; // Tensor
 
    double helampSM(double mtau, int tauhel, int Dhel, double w, double costau) const; // SM
 
    double helampV1(double mtau, int tauhel, int Dhel, double w, double costau) const; // V-A
 
    double helampV2(double mtau, int tauhel, int Dhel, double w, double costau) const; // V+A
 
    double helampS1(double mtau, int tauhel, int Dhel, double w, double costau) const; // S+P
 
    double helampS2(double mtau, int tauhel, int Dhel, double w, double costau) const; // S-P
 
    double helampT(double mtau, int tauhel, int Dhel, double w, double costau) const;  // Tensor
 
    double eta(int whel) const {return (whel == 2) ? -1 : 1;}
 
 
 
    // Vector and axial-vector form factors in terms of CLN form factors
 
    double hp(double w) const;
 
    double hm(double w) const;
 
    double hA1(double w) const;
 
    double hV(double w) const;
 
    double hA2(double w) const;
 
    double hA3(double w) const;
 
    // Scalar form factors in terms of V and A form factors
 
    double hS(double w) const;
 
    double hP(double w) const;
 
    // Tensor form factors in terms of V and A form factors
 
    double hT(double w) const;
 
    double hT1(double w) const;
 
    double hT2(double w) const;
 
    double hT3(double w) const;
 
    double z(double w) const;
 
    double ffV1(double w) const;
 
    double ffS1(double w) const;
 
    double ffA1(double w) const;
 
    double ffR1(double w) const;
 
    double ffR2(double w) const;
 
    double ffR3(double w) const;
 
    double ffV11() const {return 1.;} // cancels in R(D)
 
    double ffA11() const {return 1.;} // cancels in R(Ds)
 
    double dS1(double w) const;
 
    double dR3(double w) const;
 
    double aS1() const {return getAS1();} // {return 1.;}
 
    double aR3() const {return getAR3();} // {return 1.;}
 
    double mD(int Dhel) const;
 
    double r(int Dhel) const {return mD(Dhel) / m_mB;} // meson mass ratio
 
 
    double rq() const {return m_mCharm / m_mBottom;} // quark mass ratio
 
    double v(double mtau, double q2) const;
 
    double q2(int Dhel, double w) const;
 
    double qh2(int Dhel, double w) const;
 
    // range of q2 and w
 
    double q2min(double mtau) const {return mtau * mtau;}
 
    double q2max(int Dhel) const {return (m_mB - mD(Dhel)) * (m_mB - mD(Dhel));}
 
    double wfunc(int Dhel, double q2) const {return (1. + r(Dhel) * r(Dhel) - q2 / m_mB / m_mB) / 2. / r(Dhel);}
 
    double wmin() const {return 1.;};
 
    double wmax(double mtau, int Dhel) const {return wfunc(Dhel, q2min(mtau));}
 
    double getRho12() const {return m_rho12;}
 
    double getRhoA12() const {return m_rhoA12;}
 
    double getR11() const {return m_ffR11;}
 
    double getR21() const {return m_ffR21;}
 
    double getAS1() const {return m_aS1;}
 
    double getAR3() const {return m_aR3;}
 
    double getMB() const {return m_mB;}
 
    double getMD() const {return m_mD;}
 
    double getMDst() const {return m_mDst;}
 
    double getMBottom() const {return m_mBottom;}
 
    double getMCharm() const {return m_mCharm;}
 
    EvtComplex getCV1() const {return m_CV1;}
 
    EvtComplex getCV2() const {return m_CV2;}
 
    EvtComplex getCS1() const {return m_CS1;}
 
    EvtComplex getCS2() const {return m_CS2;}
 
    EvtComplex getCT() const {return m_CT;}
 
    void setRho12(double v) {m_rho12 = v;}
 
    void setRhoA12(double v) {m_rhoA12 = v;}
 
    void setR11(double v) {m_ffR11 = v;}
 
    void setR21(double v) {m_ffR21 = v;}
 
    void setAS1(double v) {m_aS1 = v;}
 
    void setAR3(double v) {m_aR3 = v;}
 
    void setMB(double m)  {m_mB = m;}
 
    void setMD(double m)  {m_mD = m;}
 
    void setMDst(double m) {m_mDst = m;}
 
    void setMBottom(double m) {m_mBottom = m;}
 
    void setMCharm(double m) {m_mCharm = m;}
 
    void setCV1(const EvtComplex& v) {m_CV1 = v;}
 
    void setCV2(const EvtComplex& v) {m_CV2 = v;}
 
    void setCS1(const EvtComplex& v) {m_CS1 = v;}
 
    void setCS2(const EvtComplex& v) {m_CS2 = v;}
 
    void setCT(const EvtComplex& v) {m_CT = v;}
 
  private:
    // Parameters
 
    // physics constant <-- not affect distributions
    //double hbar;
    //double GF;
    //double Vcb;
 
    double m_rho12;
 
    double m_rhoA12;
 
    double m_ffR11;
 
    double m_ffR21;
 
    double m_aS1;
 
    double m_aR3;
 
    double m_mB;
 
    double m_mD;
 
    double m_mDst;
 
    double m_mBottom;
 
    double m_mCharm;
 
    EvtComplex m_CV1;
 
    EvtComplex m_CV2;
 
    EvtComplex m_CS1;
 
    EvtComplex m_CS2;
 
    EvtComplex m_CT;
 
    bool chkDhel(int Dhel) const;
 
    bool chkwhel(int whel) const;
 
    bool chktauhel(int tauhel) const;
 
    //bool chkcostau(double costau) const;
 
  };
 
} // Belle 2 Namespace