EvtDToKmPipPipPi0.hh

Go to the documentation of this file.

#ifndef EVTDTOKMPIPPIPPI0_HH
#define EVTDTOKMPIPPIPPI0_HH
 
#include "../EvtGenBase/EvtDecayProb.hh"
#include "TComplex.h"
 
class EvtParticle;
 
class EvtDToKmPipPipPi0 : public EvtDecayProb {
 
public:
  EvtDToKmPipPipPi0() {}
  virtual ~EvtDToKmPipPipPi0();
 
  void          getName( std::string& name );
  EvtDecayBase* clone();
 
  void init();
  void initProbMax();
  void decay( EvtParticle* p );
 
private:
  double calDalEva( double Km[], double Pip1[], double Pip2[], double Pi0[] );
  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] );
  void   calt1( double daug1[4], double daug2[4], double t1[4] );
  void   calt2( double daug1[4], double daug2[4], double t2[4][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 r, int l );
  double widl1( double mass2, double mass, double sa, double sb, double sc, double r );
  void propagatorRBW( double mass2, double mass, double width, double sa, double sb, double sc,
                      double r, int l, double prop[2] );
  void propagatorRBW_K1( double mass2, double mass, double width, double sa, double sb,
                         double sc, double r, int l, double prop[2] );
  void propagatorNBW( double mass2, double mass, double width, double sa, double sb, double sc,
                      double r, int l, double prop[2] );
  void propagatorRBWl1( double mass2, double mass, double width, double sa, double sb,
                        double sc, double r, double prop[2] );
  void propagatorGS( double mass2, double mass, double width, double sa, double sb, double sc,
                     double r, double prop[2] );
  void KPiSLASS( double sa, double sb, double sc, double prop[2] );
 
  void calEvaMy( double* Km, double* Pip1, double* Pip2, double* Pi0, double* mass1,
                 double* mass2, double* width1, double* width2, double* amp, double* phase,
                 int* g0, int* g1, int* g2, int* modetype, int nstates, double& Result );
 
  TComplex ResonanceSkm( double& m2 );
  double   CalRho4pi( double s );
 
  double rho[15];
  double phi[15];
  double mD;
  double mKst0;
  double mrho;
  double mrho1450;
  double mK1400;
  double mK1460;
  double mK1270;
  double mK1650;
  double mK1680;
  // double a1_width[2000];
  // double K1270_width[2000];
  // double K1400_width[2000];
 
  double GKst0;
  double Grho;
  double Grho1450;
  double GK1400;
  double GK1460;
  double GK1270;
  double GK1650;
  double GK1680;
 
  double mass_Pion;
  double mass_Pion_N;
  double mass_Eta;
  double mass_Kaon;
  double math_pi;
  double rD2;   // 5*5
  double rRes1; // 3*3
  double rRes2; // 3*3
 
  double GS1;
  double GS2;
  double GS3; //  1/(2*math_2pi)
  double GS4; //  mass_Pion2/math_pi
  int    G[4][4];
  int    E[4][4][4][4];
};
 
#endif