EvtDToKSpieta.hh

Go to the documentation of this file.

#ifndef EVTDTOKSPIETA_HH
#define EVTDTOKSPIETA_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
 
class EvtParticle;
 
class EvtDToKSpieta:public EvtDecayProb  {
 
public:
 
   EvtDToKSpieta() {}
   virtual ~EvtDToKSpieta();
 
   void getName(std::string& name);
   EvtDecayBase* clone();
 
   void init();
   void initProbMax();
   void decay(EvtParticle *p);
 
private:
    double calPDF(double Ks[], double Pic[], 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(int l, double sa, double sb, double sc, double r, double mass);
    double Barrier(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 propagatorFlatte(double mass, double width, double sa, double sb, double sc, int r, double prop[2]);
    void propagatorGS(double mass2, 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]);
    double DDalitz(double P1[4], double P2[4], double P3[4], int Ang, double mass);
   
   double G[4][4],E[4][4][4][4];
   double rho[24];
   double phi[24];
   int spin[10];
   double modetype[20];
   double mKstr0, GKstr0;
   double ma0, Ga0, mK1270, GK1270, mK1400, GK1400, mKn_1430, GKn_1430;
   double rho_omega, phi_omega;
 
   double mD; 
   double mass_Ks; 
   double mass_Eta; 
   double rD; 
   double metap;
   double mkstr;
   double mk0;
   double mass_Kaon;
   double mass_Pion;
   double mass_Pi0;
   double math_pi;
 
   double pi;
   double mpi;
   double g1;
   double g2;
 
};
 
#endif