EvtDToKpienu.hh

Go to the documentation of this file.

 
/***********************************************************************
* Copyright 1998-2020 CERN for the benefit of the EvtGen authors       *
*                                                                      *
* This file is part of EvtGen.                                         *
*                                                                      *
* EvtGen is free software: you can redistribute it and/or modify       *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation, either version 3 of the License, or    *
* (at your option) any later version.                                  *
*                                                                      *
* EvtGen is distributed in the hope that it will be useful,            *
* but WITHOUT ANY WARRANTY; without even the implied warranty of       *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the        *
* GNU General Public License for more details.                         *
*                                                                      *
* You should have received a copy of the GNU General Public License    *
* along with EvtGen.  If not, see <https://www.gnu.org/licenses/>.     *
***********************************************************************/
 
#ifndef EVTDTOKPIENU_HH
#define EVTDTOKPIENU_HH
 
#include "EvtGenBase/EvtComplex.hh"
#include "EvtGenBase/EvtDecayProb.hh"
#include "EvtGenBase/EvtVector4R.hh"
 
#include <array>
 
class EvtParticle;
 
class EvtDToKpienu : public EvtDecayProb {
  public:
    std::string getName() const override;
    EvtDecayBase* clone() const override;
 
    void init() override;
    void initProbMax() override;
    void decay( EvtParticle* p ) override;
 
  private:
    void KinVGen( const EvtVector4R& vp4_K, const EvtVector4R& vp4_Pi,
                  const EvtVector4R& vp4_Lep, const EvtVector4R& vp4_Nu,
                  const int charm, double& m2, double& q2, double& cosV,
                  double& cosL, double& chi ) const;
    double calPDF( const double m2, const double q2, const double cosV,
                   const double cosL, const double chi ) const;
    void ResonanceP( const double m, const double q, const double mV,
                     const double mA, const double V_0, const double A1_0,
                     const double A2_0, const double m0, const double width0,
                     const double rBW, double& amplitude, double& delta,
                     EvtComplex& F11, EvtComplex& F21, EvtComplex& F31 ) const;
    void NRS( const double m, const double q, const double rS, const double rS1,
              const double a_delta, const double b_delta, const double mA,
              const double m0, const double width0, double& amplitude,
              double& delta, EvtComplex& F10 ) const;
    void ResonanceD( const double m, const double q, const double mV,
                     const double mA, const double TV_0, const double T1_0,
                     const double T2_0, const double m0, const double width0,
                     const double rBW, double& amplitude, double& delta,
                     EvtComplex& F12, EvtComplex& F22, EvtComplex& F32 ) const;
    double getPStar( const double m, const double m1, const double m2 ) const;
    double getF1( const double m, const double m0, const double m_c1,
                  const double m_c2, const double rBW ) const;
    double getF2( const double m, const double m0, const double m_c1,
                  const double m_c2, const double rBW ) const;
    double getWidth0( const double m, const double m0, const double m_c1,
                      const double m_c2, const double width0 ) const;
    double getWidth1( const double m, const double m0, const double m_c1,
                      const double m_c2, const double width0,
                      const double rBW ) const;
    double getWidth2( const double m, const double m0, const double m_c1,
                      const double m_c2, const double width0,
                      const double rBW ) const;
    EvtComplex getCoef( const double rho, const double phi ) const;
 
    int m_nAmps;
    std::array<int, 5> m_type;
 
    double m_rS;
    double m_rS1;
    double m_a_delta;
    double m_b_delta;
    double m_m0_1430_S;
    double m_width0_1430_S;
 
    double m_mV;
    double m_mA;
    double m_V_0;
    double m_A1_0;
    double m_A2_0;
    double m_m0;
    double m_width0;
    double m_rBW;
    double m_rho;
    double m_phi;
    double m_m0_1410;
    double m_width0_1410;
    double m_rho_1410;
    double m_phi_1410;
    double m_TV_0;
    double m_T1_0;
    double m_T2_0;
    double m_m0_1430;
    double m_width0_1430;
    double m_rho_1430;
    double m_phi_1430;
 
    double m_mD;
    double m_mPi;
    double m_mK;
    double m_Pi;
    double m_root2;
    double m_root2d3;
    double m_root1d2;
    double m_root3d2;
};
 
#endif