EvtDsToKSKmPipPip.hh

Go to the documentation of this file.

#ifndef EVTDSTOKSKMPIPPIP_HH
#define EVTDSTOKSKMPIPPIP_HH
 
#include "../EvtGenBase/EvtDecayProb.hh"
 
class EvtParticle;
 
class EvtDsToKSKmPipPip : public EvtDecayProb {
 
public:
  EvtDsToKSKmPipPip() {}
  virtual ~EvtDsToKSKmPipPip();
 
  void          getName( std::string& name );
  EvtDecayBase* clone();
 
  void init();
  void initProbMax();
  void decay( EvtParticle* p );
 
private:
  double calDalEva( double Ks[], double Km[], double Pip[], double Pip2[] );
  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 r2 );
  void   calt1( double daug1[4], double daug2[4], double t1[4] );
  void   calt2( double daug1[4], double daug2[4], double t2[4][4] );
  double wid( double mass, double sa, double sb, double sc, double r, int l );
  void   Flatte_rhoab( double sa, double sb, double sc, double rho[2] );
  void   propagatora0980( double mass, double sx, double* sb, double* sc, double prop[2] );
  void   propagatorRBW( double mass, double width, double sa, double sb, double sc, double r2,
                        int l, double prop[2] );
  void   KPiSLASS( double sa, double sb, double sc, double prop[2] );
  void   calEvaMy( double* Ks, double* Km, double* Pip, double* Pip2, double* mass1,
                   double* width1, double* amp, double* phase, int* g0, int* modetype,
                   int nstates, double& Result );
 
  double mass[20];
  double width[20];
  int    modetype[20];
  double rho[20];
  double phi[20];
  double amp[20];
  double phase[20];
 
  double mD;
  double metap;
  double mk0;
  double mass_Kaon;
  double mass_Eta;
  double mass_Pion;
  double math_pi;
  double mass_Pion2;
  double mass_2Pion;
  double math_2pi;
 
  double G1510;
  double GKst0;
  double GKstp;
  double Ga0_980;
  double Geta1475;
  double Gf1285;
  double m1510;
  double mKst0;
  double mKstp;
  double ma0_980;
  double meta1475;
  double mf1285;
  double rD2;
  double rRes2;
  double g1;
  double g2;
 
  double GS1;
  double GS2;
  double GS3;
  double GS4;
 
  int G[4][4];
  int E[4][4][4][4];
};
 
#endif