EvtDToKSKSK.hh

Go to the documentation of this file.

#ifndef EVTDTOKSKSK_HH 
#define EVTDTOKSKSK_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
#include "TComplex.h"
 
class EvtParticle;
 
class EvtDToKSKSK:public  EvtDecayProb  {
 
    public:
        EvtDToKSKSK() {}
        virtual ~EvtDToKSKSK();
 
        void getName(std::string& name);
        EvtDecayBase* clone();
 
        void init();
        void initProbMax();
        void decay(EvtParticle *p);
 
    private:
        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(int l, double sa, double sb, double sc, double r, double mass);
        void calt1(double daug1[4], double daug2[4], double t1[4]);
        void calt2(double daug1[4], double daug2[4], double t2[4][4]);
        void propagatorCBW(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 mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2]);
        void propagatorFlatte(double mass, double width, double sa, double sb, double sc, double prop[2]);
        void propagatorGS(double mass, double width, double sa, double sb, double sc, double r2, double prop[2]);
        void PiPiSWASS(double sa, double sb, double sc, double prop[2]);
        void KPiSLASS(double sa, double sb, double sc, double prop[2]);
        void Flatte_rhoab(double sa, double sb, double sc, double rho[2]);
        void propagatorKstr1430(double mass, double sx, double *sb, double *sc, 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 getprop(double sa, double sb, double sc, double mass, double width, double prop[2]);
        double DDalitz(double P1[4], double P2[4], double P3[4], int Ang, double mass);
        void calEva(double* Ks01, double* Ks02, double*  Kp, double *mass1, double *width1, double *amp, double *phase,int* g0,int* spin, int* modetype, int nstates, double & Result);   
 
 
 
        int modetype[3];
        double mass[3];
        double width[3];
        double rho[3];
        double phi[3];
//      double spin[3];
        double Ga0_980;
        double Grho_1450;
        double ma0_980;
        double mrho_1450;
        double mD0  = 1.86966;
        double mK0  = 0.497611;
        double mKa  = 0.493677;
        double mPi  = 0.13957;
        double mK02 = 0.237616707;
        double mPi2 = 0.01947978;
        double mass_EtaP = 0.95778;
        double mass_Kaon = 0.49368;
 
        double math_pi = 3.1415926;
        double mass_Pion = 0.13957;
        double mass_Pion2 = 0.0194797849;
        double mass_2Pion = 0.27914;
        double math_2pi = 6.2831852;
        double rD2 = 25.0;    
        double rRes2 = 9.0;   
        double g2  = 0.23;            
        double GS1 = 0.636619783;
        double GS2 = 0.23459086; 
        double GS3 = 0.1591549458;   
        double GS4 = 0.078196955;  
 
        double rho_omega = 0.00294;
        double phi_omega = -0.02;
        double math_K = 0.493677;
 
        int G[4][4];
};
 
#endif