D0ToKSLKK.h

Go to the documentation of this file.

#ifndef D0TOKSLKK_HH
#define D0TOKSLKK_HH
 
#include <complex>
#include <vector>
 
using namespace std;
 
class D0ToKSLKK {
 
public:
  D0ToKSLKK() {}
  virtual ~D0ToKSLKK();
 
  void init( int Daug0Id, int Uspin );
 
  complex<double> Amp( vector<double> k0l, vector<double> kp, vector<double> km, int Daug0Id,
                       int Uspin );
 
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] );
  void calt1( double daug1[4], double daug2[4], double t1[4] );
  void calt2( double daug1[4], double daug2[4], double t2[4][4] );
  void propagatorCBW( double mass, double width, double sx, double prop[2] );
  void propagatorFlatte( double mass, double width, double sa, double prop[2] );
  void propagatorGS( double mass, double width, double sa, double sb, double sc, double r2,
                     double prop[2] );
  void calEva( double* K, double* Pi1, double* Pi2, double* mass1, double* width1, double* amp,
               double* phase, int* g0, int* spin, int* modetype, double* r0, double* r1,
               double& Result, int first, int last, int charge, bool SorL );
  void calEva_QC( double* K, double* Pi1, double* Pi2, double* mass1, double* width1,
                  double* amp, double* phase, int* g0, int* spin, int* modetype, double* r0,
                  double* r1, double& Result, int first, int last, int charge, bool SorL );
  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 );
  void   propagatorRBW( double mass, double width, double sa, double sb, double sc, double r2,
                        int l, double prop[2] );
  void   propagatorKstr1430( double mass, double sx, double* sb, double* sc, double prop[2] );
  void   Flatte_rhoab( double sa, double sb, double sc, double rho[2] );
  void   propagatorsigma500( double sa, double sb, double sc, double prop[2] );
  void   rhoab( double sa, double sb, double sc, double res[2] );
  void   rho4Pi( double sa, double res[2] );
  void   KPiSLASS( double sa, double sb, double sc, double prop[2] );
  void   propagatora0980p( double mass, double sx, double* sb, double* sc,
                           double prop[2] ); // a0980p Flatte
  void   propagatora0980pfloated( double mass, double sx, double* sb, double* sc, double gKK,
                                  double prop[2] ); // a0980p Flatte
  void   propagatora0980wm( double mass, double width, double sx, double sb, double sc,
                            double prop[2] ); // a0980wm Flatte
  void   propagatora09800( double mass, double sx, double* sb, double* sc,
                           double prop[2] ); // a09800 Flatte
  void   propagatora09800floated( double mass, double sx, double* sb, double* sc, double gKK,
                                  double prop[2] ); // a09800 Flatte
  void   propagatora098002channel( double mass, double sx, double* sb, double* sc,
                                   double prop[2] ); // a09800 Flatte
  void   getprop( double sa, double sb, double sc, double mass, double width,
                  double prop[2] ); // rho propagator
  void calPDF( double* Ks0, double* K1, double* K2, double* mass1, double* width1, double* amp,
               double* phase, int* g0, int* spin, int* modetype, double* r0, double* r1,
               int first, int last, double PDF[2] );
 
  int Narg;
  // int Uspin;
  // int Daug0Id;
  int    charge;
  bool   SorL;
  int    modetype[8];
  double mass[8];
  double width[8];
  double rho[8];
  double phi[8];
  double spin[8];
  double mDM;
  double mK0;
  double mKa;
  double mPi;
  double mEta;
  double mKa2;  // 0.49368^2;
  double mPi2;  // 0.13957^2;
  double mEta2; // 0.547862^2;
  double mass_EtaP;
  double mass_Kaon;
  double mass_KS;
 
  double math_pi;
  double mass_Pion2;
  double mass_2Pion;
  double math_2pi;
  double rD2;   // 5*5
  double rRes2; // 3*3
  double g2;    // K*0(1430)
  double rho_omega;
  double phi_omega;
 
  double GS1;
  double GS2;
  double GS3;
  double GS4;
  int    G[4][4];
};
 
#endif