EvtDsToEta3pi.hh

Go to the documentation of this file.

#ifndef EVTDSTOETA3PI_HH
#define EVTDSTOETA3PI_HH
 
#include "../EvtGenBase/EvtDecayProb.hh"
 
class EvtParticle;
 
class EvtDsToEta3pi : public EvtDecayProb {
 
public:
  EvtDsToEta3pi() {}
  virtual ~EvtDsToEta3pi();
 
  void          getName( std::string& name );
  EvtDecayBase* clone();
 
  void init();
  void initProbMax();
  void decay( EvtParticle* p );
 
private:
  void   calEvaMy( double* Pip1, double* Pip2, double* Pim, double* Eta, double* mass1,
                   double* mass2, double* width1, double* width2, double* amp, double* phase,
                   int* g0, int* g1, int* g2, int* modetype, int nstates, 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 );
  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   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   propagatora0_980( double mass2, double mass, double width, double sa, double sb,
                           double sc, int charge, double prop[2] );
  void   Flatte_rhoab( double sa, double sb, double sc, double rho[2] );
  void   propagator980( double mass, double sx, double* sb, double* sc, double prop[2] );
  void   propagatora0980( double mass, double sx, double* sb, double* sc, double prop[2] );
  void propagatorRBW( double mass2, double mass, double width, double sa, double sb, double sc,
                      double r2, int l, double prop[2] );
  void propagatorGS( double mass2, double mass, double width, double sa, double sb, double sc,
                     double r2, double prop[2] );
 
  double mass1[11];
  double mass2[11];
  double width1[11];
  double width2[11];
  int    modetype[11];
  double rho[11];
  double phi[11];
 
  double mDs;
  double rRes;
  double rRes2;
  double rD;
  double rD2;
  double metap;
  double mkstr;
  double mk0;
  double mass_Ks;
  double mass_Kaon;
  double mass_Pion;
  double mass_Pi0;
  double math_pi;
  double mrho;
  double mf0980;
  double msigma;
  double ma0980;
  double ma11260;
  double meta;
  double mf11285;
  double meta1475;
  double mzero;
  double mf11420;
  double Gf11420;
  double Grho;
  double Gsigma;
  double Ga0980;
  double Ga11260;
  double Gf0980;
  double Gf11510;
  double Geta1475;
  double Gzero;
  double GS1;
  double GS2;
  double GS3;
  double GS4;
  int    G[4][4];
  int    E[4][4][4][4];
};
 
#endif