EvtDsTopipipipi0pi0.hh

Go to the documentation of this file.

#ifndef EVTDSTOPIPIPIPI0PI0_HH
#define EVTDSTOPIPIPIPI0PI0_HH
 
#include "../EvtGenBase/EvtDecayProb.hh"
#include "TComplex.h"
 
class EvtParticle;
 
class EvtDsTopipipipi0pi0 : public EvtDecayProb {
 
public:
  EvtDsTopipipipi0pi0() {}
  virtual ~EvtDsTopipipipi0pi0();
 
  void          getName( std::string& name );
  EvtDecayBase* clone();
 
  void init();
  void initProbMax();
  void decay( EvtParticle* p );
 
private:
  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 r2 );
  double BarrierN( 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   calG2( double daug1[4], double daug2[4], double g2[4][4] );
  void   propagator( double mass2, double mass, double width, double sx, double prop[2] );
  double widttt( double mass2, double mass, double sa, double sb, double sc, double r2, int l,
                 double ttt[3] );
  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 propagatorRBW( double mass2, double mass, double width, double sa, double sb, double sc,
                      double r2, int l, double prop[2] );
  void propagatorRBWttt( double mass2, double mass, double width, double sa, double sb,
                         double sc, double r2, int l, double ttt[6], double prop[2] );
  void propagatorRBWl1( double mass2, double mass, double width, double sa, double sb,
                        double sc, double r2, double prop[2] );
  void propagatorGS( double mass2, double mass, double width, double sa, double sb, double sc,
                     double r2, double prop[2] );
  void propagatorGS_Omg( double mass2, double mass, double width, double sa, double sb,
                         double sc, double r2, double amp_Omg, double phs_Omg,
                         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 KPiSLASS( double sa, double sb, double sc, double prop[2] );
  void calEva( double* Pip1, double* Pip2, double* Pi01, double* Pi02, double* Pim,
               double* mass1, double* mass2, double* width1, double* width2, double* amp,
               double* phase, int* g0, int* g1, int* g2, int* modetype, double nstates,
               double& Result );
  // void calEva(double* Pip1, double* Pip2, double* Pi01, double* Pi02, double* Pim, double*
  // mass1,  double* mass2, double* width1, double* width2, double* amp, double* phase, int*
  // g0, int* g1, int* g2, int* modetype, double nstates, double &Result, double* PPP1, double*
  // PPP2);
  void propagatorRBW_a1( double mass2, double mass, double width, double sa, double sb,
                         double sc, double r2, int l, double prop[2] );
 
  int charge;
 
  int    modetype[50];
  double mass1[50];
  double mass2[50];
  double width1[50];
  double width2[50];
  double rho[50];
  double phi[50];
  double spin[50];
  double g0[50];
  double g1[50];
  double Gb1;
  double Gomega;
  double Gphi;
  double Grho;
  double Grho1450;
  double mb1;
  double momega;
  double mphi;
  double mrho;
  double mrho1450;
 
  double mDsM;
  double mD;
  double metap;
  double mkstr;
  double mk0;
  double mass_Kaon;
  double mass_Pion;
  double math_pi;
  double mass_Pion2;
  double mass_2Pion;
  double math_2pi;
  double rD2;   // 5*5
  double rDs2;  // 5*5
  double rRes2; // 3*3
  double gg1;
  double gg2; // K*0(1430)
  double GS1;
  double GS2;
  double GS3; //  1/(2*math_2pi)
  double GS4;
 
  double G[4][4];
 
  double E[4][4][4][4];
};
 
#endif