EvtD0ToKSpi0pi0pi0.hh

Go to the documentation of this file.

#ifndef EVTD0TOKSPI0PI0PI0_HH
#define EVTD0TOKSPI0PI0PI0_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
class EvtParticle;
 
class EvtD0ToKSpi0pi0pi0:public EvtDecayProb  {
 
    public:
 
        EvtD0ToKSpi0pi0pi0() {}
        virtual ~EvtD0ToKSpi0pi0pi0();
 
        void getName(std::string& name);
        EvtDecayBase* clone();
 
        void init();
        void initProbMax();
        void decay(EvtParticle *p);
 
    private:
        double 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];
        int 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