EvtDToKppipipi0.hh

Go to the documentation of this file.

#ifndef EVTDTOKPPIPIPI0_HH
#define EVTDTOKPPIPIPI0_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
class EvtParticle;
 
class EvtDToKppipipi0:public EvtDecayProb  {
 
public:
 
   EvtDToKppipipi0() {}
   virtual ~EvtDToKppipipi0();
 
   void getName(std::string& name);
   EvtDecayBase* clone();
 
   void init();
   void initProbMax();
   void decay(EvtParticle *p);
 
private:
   double calPDF(double Kp[], double Pip[], 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 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 KPiSLASS(double sa, double sb, double sc, double prop[2]);
 
   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;
   double ma1, Ga1, mK1270, GK1270, mK1400, GK1400, mK1460, GK1460, mK1650, GK1650;
   double mD; 
   double rD; 
   double math_pi;
 
   double pi;
};
 
#endif