D0Topipipi0.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of models developed at BES collaboration
//      based on the EvtGen framework.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/BesCopyright
//      Copyright (A) 2006      Ping Rong-Gang @IHEP
//
// Module:  EvtD0Topipipi0.cc
//
// Description: Model provided by user, see the mannual
//
// Modification history:
//
//    Liaoyuan Dong       Aug. 8, 2022       Module created
//
//------------------------------------------------------------------------
 
#ifndef D0TOPIPIPI0_H
#define D0TOPIPIPI0_H
 
#include <vector>
#include <complex>
 
using namespace std;
 
class D0Topipipi0{
 
public:
  
  D0Topipipi0() {}
  virtual ~D0Topipipi0();
 
  void init();
 
    complex<double> Amp(vector<double> Pip, vector<double> Pim, vector<double> Pi0);
 
private:
 
  int _nd;
  int charm;
  int tagmode;
 
    void setInput(double* pip, double* pim, double* pi0);
  int CalAmp();
 
  complex<double> m_AmpD0;
  complex<double> m_AmpDb;
  complex<double> m_RAmp;
  complex<double> Get_AmpD0(){return m_AmpD0;}
  complex<double> Get_AmpDb(){return m_AmpDb;}
 
  double Get_strongPhase();
  double arg(complex<double> x);
  double mag2(complex<double> x);
    double mag2_itf(complex<double> x, complex<double> y);
 
  complex<double> CalD0Amp();
  complex<double> CalDbAmp();
 
  vector<double> sum_tensor(vector<double> pa, vector<double> pb);
  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);
  vector<double> ProjectionTensors(vector<double> pa, int rank);
 
  double fundecaymomentum(double mr2, double m1_2, double m2_2);
  double fundecaymomentum2(double mr2, double m1_2, double m2_2);
  complex<double>  breitwigner(double mx2, double mr, double wr);
  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);
  double wid(double mass, double sa, double sb, double sc, double r, int l);
 
  complex<double> GS(double mx2, double mr, double wr, double m1_2, double m2_2, double r, int l);
    complex<double> irho(double mr2, double m1_2, double m2_2);
    complex<double> Flatte(double mx2, double mr, double g1, double g2, double m1a, double m1b, double m2a, double m2b);
    complex<double> RBW(double mx2, double mr, double wr, double m1_2, double m2_2, double r, int l);
    double rho22(double sc);
    complex<double> rhoMTX(int i, int j, double s);
    complex<double> KMTX(int i, int j, double s);
    complex<double> IMTX(int i, int j);
    complex<double> FMTX(double Kijx, double Kijy, double rhojjx, double rhojjy, int i, int j);
    double FINVMTX(double s, double *FINVx, double *FINVy);
    complex<double> PVTR(int ID, double s);
    complex<double> Fvector(double sa, double s0, int l);
 
    vector<double> p4_Pip;
    vector<double> p4_Pim;
    vector<double> p4_Pi0;
  
    vector<double> g_uv;
  vector<double> epsilon_uvmn;
  vector< complex<double> > fitpara;
 
  double math_pi;
  double mass_Pion;
 
  double rRes;
  double rD;
    double m_Pi;
    double m2_Pi;
    double m_Pi0;
    double m2_Pi0;
 
    double m0_rho7700;
    double w0_rho7700;
 
    double m0_rho770p;
    double w0_rho770p;
 
    double m0_f21270;
    double w0_f21270;
 
    double s0_prod;
 
 
  double delta;
 
};
 
#endif