EvtD0ToKSKppimpi0.hh

Go to the documentation of this file.

#ifndef EVTD0TOKSKPPIMPI0_HH
#define EVTD0TOKSKPPIMPI0_HH
 
#include "../EvtGenBase/EvtDecayProb.hh"
class EvtParticle;
 
class EvtD0ToKSKppimpi0 : public EvtDecayProb {
 
public:
  EvtD0ToKSKppimpi0() {}
  virtual ~EvtD0ToKSKppimpi0();
 
  void          getName( std::string& name );
  EvtDecayBase* clone();
 
  void init();
  void initProbMax();
  void decay( EvtParticle* p );
 
private:
  void calPDF( double KS[], double Kp[], double Pim[], double Pi0[], double& Result );
 
  void   Com_Multi( double a1[2], double a2[2], double res[2] );
  void   Com_Divide( double a1[2], double a2[2], double res[2] );
  double SCADot( double a1[4], double a2[4] );
  double Barrier( double mass2, int l, double sa, double sb, double sc, double r2 );
  void   calt1( double daug1[4], double daug2[4], double t1[4] );
  void   calt2( double daug1[4], double daug2[4], double t2[4][4] );
  void   propagator( double mass2, double mass, double width, double sx, double prop[2] );
  double wid( double mass2, double mass, double sa, double sb, double sc, double r2, int l );
  double widl1( double mass2, double mass, double sa, double sb, double sc, double r2 );
  void propagatorRBW( double mass2, double mass, double width, double sa, double sb, double sc,
                      double r2, int l, double prop[2] );
  void propagatorRBWl1( double mass2, double mass, double width, double sa, double sb,
                        double sc, double r2, double prop[2] );
  void propagatorRBW_a1( double mass2, double mass, double width, double sa, double sb,
                         double sc, double r2, int l, double prop[2] );
  void propagatorRBW_K1( double mass2, double mass, double width, double sa, double sb,
                         double sc, double r2, int l, double prop[2] );
  void propagatorGS( double mass2, double mass, double width, double sa, double sb, double sc,
                     double r2, double prop[2] );
  void rhoab( double sa, double sb, double sc, double res[2] );
  void rho4Pi( double sa, double res[2] );
  void propagatorsigma500( double sa, double sb, double sc, double prop[2] );
  void Flatte_rhoab( double sa, double sb, double sc, double rho[2] );
  void propagator980( double mass, double sx, double* sb, double* sc, double prop[2] );
  void propagatorFlatte( double mass, double width, double sa, double prop[2] );
  void KPiSLASS( double sa, double sb, double sc, double prop[2] );
  // double a1_width[2000], K1270_width[2000], K1400_width[2000];
 
  int    mother_c;
  double G[4][4], E[4][4][4][4];
  double rho[24];
  double phi[24];
  double modetype[20];
  double mrhop, Grhop, mrho0, Grho0, mKstr0, GKstr0, mKstrp, GKstrp, msigma, Gsigma, mOmega,
      GOmega, ma980, Ga980;
  double ma1_1260, Ga1_1260, mK1_1270, GK1_1270, mK1_1400, GK1_1400, mK1460, GK1460, mK1650,
      GK1650, meta1475, Geta1475, mD, mPi, mKa;
  double math_pi;
};
 
#endif