EvtD0ToKKpi0.hh

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:  EvtD0ToKKpi0.cc
//
// Description: Model provided by user, see the mannual
//
// Modification history:
//
//    Liaoyuan Dong       Aug. 8, 2022       Module created
//
//------------------------------------------------------------------------
 
#ifndef EVTD0TOKKPI0_HH
#define EVTD0TOKKPI0_HH
 
#include "EvtGenBase/EvtDecayProb.hh"
#include <vector>
#include <complex>
 
class EvtParticle;
 
class EvtD0ToKKpi0:public EvtDecayProb {
 
public:
  
  EvtD0ToKKpi0() {}
  virtual ~EvtD0ToKKpi0();
 
  void getName(std::string& name);
  EvtDecayBase* clone();
 
  void init();
  void initProbMax();
  void decay(EvtParticle *p); 
 
private:
 
  double AmplitudeSquare(int charm, int tagmode);
 
  EvtVector4R GetDaugMomLab(int i){return _p4Lab[i];}
  EvtVector4R GetDaugMomCM(int i) {return _p4CM[i];}
  EvtVector4R GetDaugMomHel(int i){return _p4Hel[i];}
 
  EvtVector4R _p4Lab[10],_p4CM[10],_p4Hel[10];
  int _nd;
  int charm;
  int tagmode;
 
    void setInput(double* pip, double* pim, double* pi0);
  int CalAmp();
 
  complex<double> Get_AmpD0(){return m_AmpD0;}
  complex<double> Get_AmpDb(){return m_AmpDb;}
 
  double Get_strongPhase();
    complex<double> Amp(vector<double> Pip, vector<double> Pim, vector<double> Pi0);
  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> KPiSLass(double mx2, double mr, double wr, double a, double r, double R, double phi_R, double phi_F);
    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_Kap;
    vector<double> p4_Kam;
    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 m_Ka;
    double m2_Ka;
 
    double m0_Kst892;
    double w0_Kst892;
 
    double m0_Kst1410;
    double w0_Kst1410;
 
    double m0_phi1020;
    double w0_phi1020;
 
    double m0_f0980;
    double g1_f0980;
    double g2_f0980;
 
    double m0_K01430;
    double w0_K01430;
    double a_lass;
    double r_lass;
 
    double s0_prod;
 
  complex<double> m_AmpD0;
  complex<double> m_AmpDb;
  complex<double> m_RAmp;
 
  double delta;
 
};
class EvtD0ToKKpi0:public EvtDecayProb {…};
 
#endif