EvtBSemiTauonicDecayRateCalculator.h

/**************************************************************************
 * 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 <cstddef>
 
namespace Belle2 {
  class EvtBSemiTauonicHelicityAmplitudeCalculator;
 
  class EvtBSemiTauonicDecayRateCalculator {
 
  public:
    EvtBSemiTauonicDecayRateCalculator(): m_BSTD(NULL) {}
 
    ~EvtBSemiTauonicDecayRateCalculator() {}
 
    double dGammadwdcostau(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, double mtau, int tauhel, int Dhel, double w,
                           double costau);
 
    double dGammadw(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, double mtau, int tauhel, int Dhel, double w);
 
    double dGammadcostau(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, double mtau, int tauhel, int Dhel, double costau);
 
    double Gamma(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, double mtau, int tauhel, int Dhel);
 
    double GammaD(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, double mtau);
 
    double GammaDstar(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, double mtau);
 
    double GammaSMD(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, const double mlep = 0.0005110);
 
    double GammaSMDstar(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, const double mlep = 0.0005110);
 
    double RGammaD(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, const double mtau, const double mlep = 0.0005110);
 
    double RGammaDstar(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, const double mtau, const double mlep = 0.0005110);
 
    double PtauD(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, const double mtau);
 
    double PtauDstar(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, const double mtau);
 
    double PDstar(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, const double mtau);
 
    double pf(const EvtBSemiTauonicHelicityAmplitudeCalculator& BSTD, double mtau, int Dhel, double w);
 
  private:
    double EvaluateByW(double* x, double* param);   // x=w
 
    double EvaluateByCostau(double* x, double* param); // x=costau
 
    double EvaluateBy2D(double* x, double* param);  // x={w,costau}
 
    const EvtBSemiTauonicHelicityAmplitudeCalculator* m_BSTD;
  };
 
} // namespace Belle2