EvtD0ToKSKK.hh

Go to the documentation of this file.

#ifndef EVTD0TOKSKK_HH
#define EVTD0TOKSKK_HH
 
#include "../EvtGenBase/EvtDecayProb.hh"
 
class EvtParticle;
 
class EvtD0ToKSKK : public EvtDecayProb {
 
public:
  EvtD0ToKSKK() {}
  virtual ~EvtD0ToKSKK();
 
  void          getName( std::string& name );
  EvtDecayBase* clone();
 
  void init();
  void initProbMax();
  void decay( EvtParticle* p );
 
private:
  double EvaluateAmp( EvtVector4R& p4_d1, EvtVector4R& p4_d2,
                      EvtVector4R& p4_d3 ); // double EvaluateAmp(EvtVector4R& p4_d1,
                                            // EvtVector4R& p4_d2, EvtVector4R& p4_d3, double
                                            // res[]);
  double helicityAngle( double M, double m, double m2, double mSpec, double invMassSqA,
                        double invMassSqB );
  std::complex<double> AmpRelBreitWignerRes( double mSq, double mR, double ma, double mb,
                                             double width, unsigned int J,
                                             double mesonRadius );
  std::complex<double> widthToCoupling( double mSq, double mR, double width, double ma,
                                        double mb, double spin, double mesonRadius );
  std::complex<double> AmpFlatteRes( double mSq, double mR, double massA1, double massA2,
                                     double gA, double massB1, double massB2, double couplingB,
                                     unsigned int J, double mesonRadius );
  std::complex<double> AmpFlatteRes( double mSq, double mR, double massA1, double massA2,
                                     double gA, double massB1, double massB2, double couplingB,
                                     double massC1, double massC2, double couplingC,
                                     unsigned int J, double mesonRadius );
  double FormFactor( double sqrtS, double ma, double mb, double spin, double mesonRadius );
  std::complex<double> phspFactor( double sqrtS, double ma, double mb );
  std::complex<double> qValue( double sqrtS, double ma, double mb );
  double               qSqValue( double sqrtS, double ma, double mb );
  std::complex<double> GetCoefficient( double mag, double phase ) const;
  double Wigner_d( unsigned int __j, unsigned int __m, unsigned int __n, double __beta );
  int    factorial( int i );
 
  double MD;
  double m1;
  double m2;
  double m3;
 
  double       mag_a00;
  double       phase_a00;
  unsigned int spin_a00;
  double       mass_a00;
  double       gA_a00;
  double       gB_a00;
  double       gC_a00;
  double       massB1_a00;
  double       massB2_a00;
  double       massC1_a00;
  double       massC2_a00;
 
  double       mag_a0p;
  double       phase_a0p;
  unsigned int spin_a0p;
  double       mass_a0p;
  double       gA_a0p;
  double       gB_a0p;
  double       massB1_a0p;
  double       massB2_a0p;
 
  double       mag_phi;
  double       phase_phi;
  unsigned int spin_phi;
  double       mass_phi;
  double       width_phi;
 
  double       mag_a2p;
  double       phase_a2p;
  unsigned int spin_a2p;
  double       mass_a2p;
  double       width_a2p;
 
  double       mag_a2m;
  double       phase_a2m;
  unsigned int spin_a2m;
  double       mass_a2m;
  double       width_a2m;
 
  double       mag_a0_1450m;
  double       phase_a0_1450m;
  unsigned int spin_a0_1450m;
  double       mass_a0_1450m;
  double       width_a0_1450m;
 
  double mesonRadius;
  double motherRadius;
};
 
#endif