EvtDsToKpPipPimPi0.hh

Go to the documentation of this file.

#ifndef EVTDSTOKPPIPPIMPI0_HH
#define EVTDSTOKPPIPPIMPI0_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
#include "TComplex.h"
 
class EvtParticle;
 
class EvtDsToKpPipPimPi0:public  EvtDecayProb  {
 
    public:
 
        EvtDsToKpPipPimPi0() {}
        virtual ~EvtDsToKpPipPimPi0();
 
        void getName(std::string& name);
        EvtDecayBase* clone();
 
        void init();
        void initProbMax();
        void decay(EvtParticle *p);
 
    private:
        double calDalEva(double Kp[], double Pip[], double Pim[], 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 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* Kp,  double* Pip, double* Pim,   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, double & q1270);
 
        TComplex ResonanceSkm(double & m2);
        double CalRho4pi(double s);
 
        double rho[11];
        double phi[11];
        double  mDsM;  
        double  mKst0; 
        double  mKstp;
        double  mKstp2;
        double  mrho;   
        double  mrho0;   
        double  mK1400; 
        double  mK1410; 
        double  mK1650; 
        double  mA1;    
        double  mK1270; 
        double  mOmega;
        double  GOmega;
        double  GKst0;  
        double  GKstp;  
        double  Grho;  
        double  Grho0;  
        double  GK1400; 
        double  GK1410; 
        double  GK1650; 
        double  GA1;    
        double  GK1270; 
 
        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