EvtDsToEtapi2pi0.hh

Go to the documentation of this file.

#ifndef EVTDSTOETAPI2PI0_HH
#define EVTDSTOETAPI2PI0_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
 
class EvtParticle;
 
class EvtDsToEtapi2pi0:public  EvtDecayProb  {
 
public:
 
  EvtDsToEtapi2pi0() {}
  virtual ~EvtDsToEtapi2pi0();
 
  void getName(std::string& name);
  EvtDecayBase* clone();
 
  void init();
  void initProbMax();
  void decay(EvtParticle *p);
 
private:
  void calEvaMy(double* Pic, double* Pi01, double* Pi02, 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);
  double barrier(int l, double sa, double sb, double sc, double r, double mass);
  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[3];
  double mass2[3];
  double width1[3];
  double width2[3];
  int  modetype[3];
  double rho[3];
  double phi[3];
 
  double mDs;
  double rRes;
  double rRes2;
  double rD;
  double rD2;
  double metap;
  double mkstr;
  double mk0;
  double meta;
  double mass_Ks;
  double mass_Kaon;
  double mass_Pion;
  double mass_Pi0;
  double math_pi;
  double mrho770_a1_S;
  double ma0980;
  double ma0980_S1; 
  double ma11260_S;
  double mrho770;
  double mzero_a01;
  double Grho770_a1_S;
  double Ga0980;
  double Ga0980_S1;
  double Ga11260_S;
  double Grho770;
  double Gzero_a01;
  double GS1;
  double GS2;
  double GS3;
  double GS4;
  int G[4][4];
  int E[4][4][4][4];
};
 
#endif