EvtD0ToKSLKK.hh

Go to the documentation of this file.

#ifndef EVTD0TOKSLKK_HH 
#define EVTD0TOKSLKK_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
#include "TComplex.h"
 
class EvtParticle;
 
class EvtD0ToKSLKK:public  EvtDecayProb  {
 
public:
   EvtD0ToKSLKK() {}
   virtual ~EvtD0ToKSLKK();
 
   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]);
   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]);
   void propagatorFlatte(double mass, double width, double sa, double prop[2]);
   void propagatorGS(double mass, double width, double sa, double sb, double sc, double r2, double prop[2]);
   void calEva(double* K, double* Pi1, double* Pi2 , double *mass1, double *width1, double *amp, double *phase,int* g0, int* spin,int* modetype, double *r0, double *r1, double & Result, int first, int last, int charge, bool SorL);
   double SCADot(double a1[4], double a2[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 r2, int l);
   double widl1(double mass2, double mass, double sa, double sb, double sc, double r2);
   double DDalitz(double P1[4], double P2[4], double P3[4], int Ang, double mass);
   void propagatorRBW(double mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2]);
     void propagatorKstr1430(double mass, double sx, double *sb, double *sc, double prop[2]);
     void Flatte_rhoab(double sa, double sb, double sc, double rho[2]);
     void propagatorsigma500(double sa, 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 KPiSLASS(double sa, double sb, double sc, double prop[2]);
     void propagatora0980p(double mass, double sx, double *sb, double *sc, double prop[2]); //a0980p Flatte
     void propagatora0980pfloated(double mass, double sx, double *sb, double *sc, double gKK, double prop[2]); //a0980p Flatte
   void propagatora0980wm(double mass, double width, double sx, double sb, double sc, double prop[2]); //a0980wm Flatte
     void propagatora09800(double mass, double sx, double *sb, double *sc, double prop[2]); //a09800 Flatte
     void propagatora09800floated(double mass, double sx, double *sb, double *sc, double gKK, double prop[2]); //a09800 Flatte
     void propagatora098002channel(double mass, double sx, double *sb, double *sc, double prop[2]); //a09800 Flatte
     void getprop(double sa, double sb, double sc, double mass, double width, double prop[2]);//rho propagator
     void calPDF(double *Ks0, double *K1, double *K2, double* mass1,   double* width1, double* amp, double* phase, int* g0, int* spin, int* modetype,   double* r0, double* r1, int first, int last, double PDF[2]);
 
 
     int Narg;
     int Uspin;
     int Daug0Id;
     int charge;
     bool SorL;
   int modetype[8];
   double mass[8];
   double width[8];
   double rho[8];
   double phi[8];
   double spin[8];
     double mDM;
     double mK0;
     double mKa;
     double mPi;
     double mEta;
     double mKa2;//0.49368^2;
     double mPi2;//0.13957^2;
     double mEta2;//0.547862^2;
     double mass_EtaP;
     double mass_Kaon;
     double mass_KS;
    
     double math_pi;
     double mass_Pion2;
     double mass_2Pion;
     double math_2pi;
     double rD2;    // 5*5
     double rRes2;   // 3*3
     double g2;             //K*0(1430)
   double rho_omega;
     double phi_omega; 
 
 
   double GS1;
   double GS2;
   double GS3;
   double GS4;
   int G[4][4];
};
 
#endif