EvtBSemiTauonicHelicityAmplitudeCalculator.h

Go to the documentation of this file.

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2013 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Koji Hara                                                *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#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