EvtDsToEtapipi0.hh

Go to the documentation of this file.

#ifndef EVTDSTOETAPIPI0_HH
#define EVTDSTOETAPIPI0_HH
 
#include "../EvtGenBase/EvtDecayProb.hh"
#include "TComplex.h"
 
class EvtParticle;
 
class EvtDsToEtapipi0 : public EvtDecayProb {
 
public:
  EvtDsToEtapipi0() {}
  virtual ~EvtDsToEtapipi0();
 
  void          getName( std::string& name );
  EvtDecayBase* clone();
 
  void init();
  void initProbMax();
  void decay( EvtParticle* p );
 
private:
  Double_t calDalEva( Double_t P1[], Double_t P2[], Double_t P3[] );
  TComplex Spin_factor( Double_t P1[], Double_t P2[], Double_t P3[], Int_t spin, Int_t flag );
  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   barrierNeo( int l, double sa, double sb, double sc, double r, double mR );
  void     calt1( double daug1[4], double daug2[4], double t1[4] );
  void     calt2( double daug1[4], double daug2[4], double t2[4][4] );
  double   wid( double mass, double sa, double sb, double sc, double r, int l );
  void     Flatte_rhoab( double sa, double sb, double sc, double rho[2] );
  TComplex Flatte_rhoab( Double_t sa, Double_t sb, Double_t sc );
  TComplex propagatorRBW( Double_t mass, Double_t width, Double_t sa, Double_t sb, Double_t sc,
                          Double_t r, Int_t l );
 
  double mass[4];
  double width[4];
  double rho[24];
  double phi[24];
  double mD;
  double mDs;
  double rRes;
  double rD;
  double metap;
  double mkstr;
  double mk0;
  double mass_Kaon;
  double mass_Pion;
  double mass_Pi0;
  double math_pi;
  double mrho;
  double Grho;
  double ma0;
  double Ga0;
  double meta;
  int    G[4][4];
};
 
#endif