D0ToKSpipipi0.h

Go to the documentation of this file.

#ifndef D0TOKSPIPIPI0_H
#define D0TOKSPIPIPI0_H
 
#include <vector>
#include <complex>
#include <map>
 
using namespace std;
 
class D0ToKSpipipi0 {
 
public:
  
  D0ToKSpipipi0() {}
  virtual ~D0ToKSpipipi0();
 
  void init();
 
    complex<double> Amp(vector<double> ks0, vector<double> pip, vector<double> pim, vector<double> pi0); 
private:
 
  //double AmplitudeSquare();
 
  void readInputFile();
  void readInputCoeff();
  void initPar();
 
  // Functions
  void addPartialWave(complex<double> amp, double mag, double pha);
  void addPartialWave(complex<double> amp1, complex<double> amp2, double mag, double pha);
  void addPartialWave(double t1, complex<double> resX, complex<double> resY, double mag, double pha);
  void addPartialWave(double t1, complex<double> resX1, complex<double> resY1, double t2, complex<double> resX2, complex<double> resY2, double fact, double mag, double pha); // For iso-spin conjugate
 
  void createPropagator(vector<double> ks0, vector<double> pip, vector<double> pim, vector<double> pi0);
  void createSpinfactor(vector<double> ks0, vector<double> pip, vector<double> pim, vector<double> pi0);
 
  // Spin Factor
  vector<double> Proj(vector<double> pa, int rank);
  double contract_11_0(vector<double> pa, vector<double> pb);
  double contract_22_0(vector<double> pa, vector<double> pb);
  vector<double> contract_21_1(vector<double> pa, vector<double> pb);
  vector<double> contract_22_2(vector<double> pa, vector<double> pb);
  vector<double> contract_31_2(vector<double> pa, vector<double> pb);
  vector<double> contract_41_3(vector<double> pa, vector<double> pb);
  vector<double> contract_42_2(vector<double> pa, vector<double> pb);
 
  vector<double> OrbitalTensors(vector<double> pa, vector<double> pb, vector<double> pc, double r, int rank);
 
  //------------------------------------------------------------//
  double D2PP_P2SP();
  double D2PP_P2VP(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2VP_V2VP(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2AP_A2SP(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2AP_A2VP(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2AP_A2TP(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2TP_T2VP(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2TP_T2TP(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2PTP_PT2SP(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2PTP_PT2VP(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2PTP_PT2TP(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2SS();
  double D2VS(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2VV(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2TS(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2TV(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
  double D2TT(vector<double> p1, vector<double> p2, vector<double> p3, vector<double> p4, int l);
 
 
  // Propagator
  double fundecaymomentum(double mr2, double m1_2, double m2_2);
  double fundecaymomentum2(double mr2, double m1_2, double m2_2);
  double wid(double mass, double sa, double sb, double sc, double r, int l);
  complex<double> BW(double mx2, double mr, double wr);
  complex<double> RBW(double mx2, double mr, double wr, double m1_2, double m2_2, double r, int l);
  complex<double> LASS(double mx2, double m1_2, double m2_2);
  
  double h(double m, double q);
  double dh(double m0, double q0);
  double f(double m0, double sx, double q0, double q);
  double d(double m0, double q0);
  complex<double> GS(double mx2, double mr, double wr, double m1_2, double m2_2, double r, int l);
 
  double rho4pi(double s);
  double rho2pi(double s);
  complex<double> sigma(double mx2, double mr, double gf);
 
  complex<double> irho(double mr2, double m1_2, double m2_2);
  complex<double> Flatte2(double mx2, double mr2, double g1, double m1a, double m1b, double g2, double m2a, double m2b);
 
  complex<double> create_RBW_propagator(string name, double mx2, double m1_2, double m2_2, int l);
  complex<double> create_BW_propagator(string name, double mx2);
  complex<double> create_GS_propagator(string name, double mx2); //for rho0
  complex<double> create_KPiSLASS_propagator(string name, double mx2, double m1_2, double m2_2);
  complex<double> create_sigma_propagator(string name, double mx2);
  complex<double> create_Flatte2_propagator(string name, double mx2, double, double, double, double);
 
 
  int _nd;
 
  double rRes;
  double rD;
 
  double math_pi;
  double mpip;
  double mpim;
  double mpion;
  double mkaon;
 
  vector<double> m_ks;
  vector<double> m_pip;
  vector<double> m_pim;
  vector<double> m_pi0;
 
  vector<double> g;
  vector<double> epsilon;
 
  complex<double> totalAmp;
 
  std::map<string,double> resonance_par;
  std::map<string,double> coefficient;
  std::map<string,double> spinfactor;
  std::map< string,complex<double> > propagator;
 
  std::vector<string> VwaveName;
  std::vector<string> VcoeffName;
  std::vector<string> Vres1;
  std::vector<string> Vres2;
  std::vector<string> Vmodetype;
  std::vector<string> Vcom1;
  std::vector<int>    Vl;
  std::vector<int>    Vconj;
};
 
#endif