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