EvtDToKSKpi0.hh

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of models developed at BES collaboration
//      based on the EvtGen framework.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/BesCopyright
//      Copyright (A) 2006      Ping Rong-Gang @IHEP
//
// Module:  EvtDToKSKpi0.hh
//
// Description: Model provided by user, see the mannual
//
// Modification history:
//
//    Liaoyuan Dong       Aug  8, 2022       Module created
//
//------------------------------------------------------------------------
 
#ifndef EVTDTOKSKPI0_HH
#define EVTDTOKSKPI0_HH
 
#include "../EvtGenBase/EvtDecayProb.hh"
 
class EvtParticle;
 
class EvtDToKSKpi0 : public EvtDecayProb {
 
public:
  EvtDToKSKpi0() {}
  virtual ~EvtDToKSKpi0();
 
  void          getName( std::string& name );
  EvtDecayBase* clone();
 
  void init();
  void initProbMax();
  void decay( EvtParticle* p );
 
private:
  double AmplitudeSquare();
 
  EvtVector4R GetDaugMomLab( int i ) { return _p4Lab[i]; }
  EvtVector4R GetDaugMomCM( int i ) { return _p4CM[i]; }
  EvtVector4R GetDaugMomHel( int i ) { return _p4Hel[i]; }
 
  EvtVector4R _p4Lab[10], _p4CM[10], _p4Hel[10];
  int         _nd;
  EvtVector4R _pd[3];
 
  double twoBodyCMmom( double rMassSq, double d1m, double d2m );
  double dampingFactorSquare( const double& cmmom, const int& spin, const double& mRadius );
  double spinFactor( int spin, double motherMass, double daug1Mass, double daug2Mass,
                     double daug3Mass, double m12, double m13, double m23 );
  double getM23( double massPZ, double massPM ) {
    return ( _mDp2 + KsMass * KsMass + KpMass * KpMass + pi0Mass * pi0Mass - massPZ - massPM );
  }
  EvtComplex RBW( int id, double resmass, double reswidth, int spin );
  EvtComplex LASS( int id, double resmass, double reswidth );
  EvtComplex Flatte( int id, double resmass, double g1, double rg2og1 );
 
  double _mDp;
  double c_motherMass;
  double _mDp2;
  double _mDp2inv;
  double KsMass;
  double KpMass;
  double pi0Mass;
  double etamass;
  double pipMass;
  double c_meson_radius_inter;
  double c_meson_radius_Dp;
};
 
#endif