EvtDTopipi0pi0.hh

Go to the documentation of this file.

#ifndef EVTDTOPIPI0PI0_HH 
#define EVTDTOPIPI0PI0_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
 
class EvtParticle;
 
class EvtDTopipi0pi0:public  EvtDecayProb  {
 
public:
   EvtDTopipi0pi0() {}
   virtual ~EvtDTopipi0pi0();
 
   void getName(std::string& name);
   EvtDecayBase* clone();
 
   void init();
   void initProbMax();
   void decay(EvtParticle *p);
 
private:
 
   void propagator(double mass2, double mass, double width, double sx, double prop[2]);
   void Com_Multi(double a1[2], double a2[2], double res[2]);
   void Com_Divide(double a1[2], double a2[2], double res[2]);
   void calt1(double daug1[4], double daug2[4], double t1[4]);
   void calt2(double daug1[4], double daug2[4], double t2[4][4]);
   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]);
   void PiPiSWASS(double sa, double sb, double sc, double prop[2]);
   void propagatorFlatte(double mass, double width, double sa, double sb, double sc, int r, double prop[2]);
   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 Flatte_rhoab(double sa, double sb, double rho[2]);
   void propagator980(double mass, double sx, double *sb, double prop[2]);
   void calEva(double* Pic, double* Pi01, double* Pi02, double *mass1, double *width1, double *amp, double *phase,int* g0,int* spin, int* modetype, int nstates, double & Result ,double* r0 , double* r1);
   double CalRho4pi(double_t s);
   double SCADot(double a1[4], double a2[4]);
   double barrier(int l, double sa, double sb, double sc, double r, double mass);
   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);
   double DDalitz(double P1[4], double P2[4], double P3[4], int Ang, double mass, double rRES0, double rRES1);
   void rhoMTX(int i, int j, double s, double Rho[2]);
   void KMTX(int i, int j, double s, double KM[2]);
   void IMTX(int i, int j, double IMTX[2]);
   void FMTX(double Kijx, double Kijy, double rhojjx, double rhojjy, int i, int j, double FM[2]);
   void FINVMTX(double s, double *FINVx, double *FINVy);
   void PVTR(int ID, double s, double PV[2]);
   void Fvector(double sa, double s0, double Fv[2]);
   TComplex ResonanceSkm(Double_t& m2);
 
   int modetype[4];
   double mass[4];
   double width[4];
   double rho[4];
   double phi[4];
   double spin[4];
   double Gf21270;
   double Gf0980;
   double Gf0500;
   double Gf01370;
   double Grho;
   double GKstm;
   double Gpipi0VP;
   double mf21270;
   double mf0980;
   double mf0500;
   double mf01370;
   double mrho;
   double mKstm;
   double mpipi0VP;
 
   double mD;
   double mDs;
   double rRes;
   double rD;
   double metap;
   double mkstr;
   double mk0;
   double mass_Kaon;
   double mass_Pion;
 
double mass_Eta ;
double mass_Etap;
 
   double mass_Pi0;
   double math_pi;
   double meta;
   double ma0;
   double Ga0;
   double GS1;
   double GS2;
   double GS3;
   double GS4;
   int G[4][4];
   double pi180inv;
   double ma[5], g[5][5];
 
};
 
#endif