EvtDToKKpi.hh

Go to the documentation of this file.

#ifndef EVTDTOKKPI_HH
#define EVTDTOKKPI_HH
 
#include "../EvtGenBase/EvtDecayProb.hh"
 
class EvtParticle;
 
class EvtDToKKpi : public EvtDecayProb {
 
public:
  EvtDToKKpi() {}
  virtual ~EvtDToKKpi();
 
  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 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 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 KPiSLASS( double sa, double sb, double sc, double prop[2] );
  void calEva( double* K01, double* K02, double* Pi, double* mass1, double* width1,
               double* amp, double* phase, int* g0, int* spin, int* modetype, int nstates,
               double& Result, double* r0, double* r1 );
  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 );
  int    modetype[8];
  double mass[8];
  double width[8];
  double rho[8];
  double phi[8];
  int    spin[8];
  double GKsta0;
  double GK1680;
  double Gf0500;
  double Gf0980;
  double Gf01370;
  double Gf21270;
  double mKsta0;
  double mK1680;
  double mf0500;
  double mf0980;
  double mf01370;
  double mf21270;
 
  double mDp;
  double mK0;
  double mKa;
  double mPi;
  double mK02;
  double mPi2;
  double mass_EtaP;
  double mass_Kaon;
  double math_pi;
  double mass_Pion;
  double mass_2Pion;
  double mass_Pion2;
  double math_2pi;
  double rD2;
  double rRes2;
  // double g1;
  double g2;
  double GS1;
  double GS2;
  double GS3;
  double GS4;
  int    G[4][4];
};
 
#endif