EvtDToKmPipPipPi0.hh

Go to the documentation of this file.

#ifndef EVTDTOKMPIPPIPPI0_HH
#define EVTDTOKMPIPPIPPI0_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
#include "TComplex.h"
 
class EvtParticle;
 
class EvtDToKmPipPipPi0:public  EvtDecayProb  {
 
    public:
 
        EvtDToKmPipPipPi0() {}
        virtual ~EvtDToKmPipPipPi0();
 
        void getName(std::string& name);
        EvtDecayBase* clone();
 
        void init();
        void initProbMax();
        void decay(EvtParticle *p);
 
    private:
        double calDalEva(double Km[], double Pip1[], double Pip2[], double Pi0[]);
        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]);
        void calt1(double daug1[4], double daug2[4], double t1[4]);
        void calt2(double daug1[4], double daug2[4], double t2[4][4]);
 
        double Barrier(int l, double sa, double sb, double sc, double r, double mass);
        double wid(double mass2,double mass, double sa, double sb, double sc, double r, int l);
        double widl1(double mass2, double mass, double sa, double sb, double sc, double r); 
        void propagatorRBW(double mass2, double mass, double width, double sa, double sb, double sc, double r, int l, double prop[2]);
        void propagatorRBW_K1(double mass2, double mass, double width, double sa, double sb, double sc, double r, int l, double prop[2]);
        void propagatorNBW(double mass2, double mass, double width, double sa, double sb, double sc, double r, int l, double prop[2]);
        void propagatorRBWl1(double mass2, double mass, double width, double sa, double sb, double sc, double r, double prop[2]);
        void propagatorGS(double mass2, double mass, double width, double sa, double sb, double sc, double r, double prop[2]);
        void KPiSLASS(double sa, double sb, double sc, double prop[2]);
 
        void calEvaMy(double* Km,  double* Pip1, double* Pip2,   double* Pi0,  double *mass1,  double *mass2, double *width1, double *width2, double *amp,  double *phase,int* g0,int* g1,int* g2,int* modetype,  int nstates,  double  & Result);
 
        TComplex ResonanceSkm(double & m2);
        double CalRho4pi(double s);
 
        double rho[15];
        double phi[15];
        double  mD;  
        double  mKst0; 
        double  mrho;   
        double  mrho1450;   
        double  mK1400; 
        double  mK1460; 
        double  mK1270;
        double  mK1650;
        double  mK1680;
        //double a1_width[2000];
        //double K1270_width[2000];
        //double K1400_width[2000];
 
        double  GKst0;  
        double  Grho;  
        double  Grho1450;  
        double  GK1400; 
        double  GK1460; 
        double  GK1270; 
        double  GK1650; 
        double  GK1680; 
 
        double mass_Pion;
        double mass_Pion_N;
        double mass_Eta;
        double mass_Kaon;
        double math_pi;
        double rD2;    // 5*5
        double rRes1;   // 3*3
        double rRes2;   // 3*3
 
        double GS1;
        double GS2;
        double GS3;   //  1/(2*math_2pi)
        double GS4;  //  mass_Pion2/math_pi
        int G[4][4];
        int E[4][4][4][4];
};
 
#endif