EvtDToKSKmpippip.hh

Go to the documentation of this file.

#ifndef EVTDTOKSKMPIPPIP_HH
#define EVTDTOKSKMPIPPIP_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
class EvtParticle;
 
class EvtDToKSKmpippip:public EvtDecayProb  {
 
    public:
 
        EvtDToKSKmpippip() {}
        virtual ~EvtDToKSKmpippip();
 
        void getName(std::string& name);
        EvtDecayBase* clone();
 
        void init();
        void initProbMax();
        void decay(EvtParticle *p);
 
    private:
                double calEva(double* KS,double* Km,double* Pip1,double* Pip2,double* mass1,double* mass2,double* width1,double* width2,double* rho,double* phi,int* g0,int* g1,int* g2,int* modetype,double value);
 
        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);
        double barrier(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 Flatte_rhoab(double sa, double sb, double sc, double rho[2]);
        void propagator980(double mass2, double mass, double sx, double *sb, double *sc, double prop[2]);
        void propagatora0980(double mass2, double mass, double sx, double *sb, double *sc, double prop[2]);
        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 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];
        int modetype[20];
   double mK1400, m1510, meta1405, meta1475, mf1285, mKstp, mKst0, ma0_980;
   double GK1400, G1510, Geta1405, Geta1475, Gf1285, GKstp, GKst0, Ga0_980;
 
   double mD, mPi, mKa;
   double math_pi;
   double GS1;
   double GS2;
   double GS3;
   double GS4;
   double mass_Pion;
   double mass_Pion_N;
   double mass_Eta;
   double rRes2;
   double rD2;
};
 
#endif