D0ToKSLKK.h

Go to the documentation of this file.

#ifndef D0TOKSLKK_HH 
#define D0TOKSLKK_HH
 
#include <vector>
#include <complex>
 
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