EvtDalitzReso Class Reference

#include <EvtDalitzReso.hh>

Public Types
enum	NumType {   NBW = 0 , RBW_ZEMACH = 1 , RBW_KUEHN = 2 , RBW_CLEO = 3 ,   RBW_ZEMACH2 = 4 , GS_CLEO = 5 , K_MATRIX = 6 , RBW_CLEO_ZEMACH = 7 ,   GS_CLEO_ZEMACH = 8 , LASS = 9 , K_MATRIX_I = 10 , K_MATRIX_II = 11 ,   GAUSS_CLEO = 12 , GAUSS_CLEO_ZEMACH = 13 }
enum	CouplingType {   Undefined =0 , PicPic =1 , PizPiz , PiPi ,   KcKc , KzKz , KK , EtaPic ,   EtaPiz , PicPicKK , WA76 }

Public Member Functions
	EvtDalitzReso ()
	EvtDalitzReso (const EvtDalitzPlot &dp, EvtCyclic3::Pair pairAng, EvtCyclic3::Pair pairRes, EvtSpinType::spintype spin, double m0, double g0, NumType typeN)
	EvtDalitzReso (const EvtDalitzPlot &dp, EvtCyclic3::Pair pairAng, EvtCyclic3::Pair pairRes, EvtSpinType::spintype spin, double m0, double g0, NumType typeN, double m0_mix, double g0_mix, double delta_mix, EvtComplex amp_mix)
	EvtDalitzReso (const EvtDalitzPlot &dp, EvtCyclic3::Pair pairAng, EvtCyclic3::Pair pairRes, EvtSpinType::spintype spin, double m0, NumType typeN, double g1, double g2, CouplingType coupling2)
	EvtDalitzReso (const EvtDalitzPlot &dp, EvtCyclic3::Pair pairRes, std::string nameIndex, NumType typeN, EvtComplex fr12prod, EvtComplex fr13prod, EvtComplex fr14prod, EvtComplex fr15prod, double s0prod)
	EvtDalitzReso (const EvtDalitzPlot &dp, EvtCyclic3::Pair pairRes, double m0, double g0, double a, double r, double B, double phiB, double R, double phiR)
	EvtDalitzReso (const EvtDalitzReso &other)
	~EvtDalitzReso ()
EvtDalitzReso *	clone () const
EvtComplex	evaluate (const EvtDalitzPoint &p)
void	set_fd (double R)
void	set_fb (double R)

Detailed Description

Definition at line 34 of file EvtDalitzReso.hh.

Member Enumeration Documentation

CouplingType

enum EvtDalitzReso::CouplingType

Enumerator
Undefined
PicPic
PizPiz
PiPi
KcKc
KzKz
KK
EtaPic
EtaPiz
PicPicKK
WA76

Definition at line 51 of file EvtDalitzReso.hh.

{Undefined=0,PicPic=1,PizPiz,PiPi,KcKc,KzKz,KK,EtaPic,EtaPiz,PicPicKK,WA76};

NumType

enum EvtDalitzReso::NumType

Enumerator
NBW
RBW_ZEMACH
RBW_KUEHN
RBW_CLEO
RBW_ZEMACH2
GS_CLEO
K_MATRIX
RBW_CLEO_ZEMACH
GS_CLEO_ZEMACH
LASS
K_MATRIX_I
K_MATRIX_II
GAUSS_CLEO
GAUSS_CLEO_ZEMACH

Definition at line 38 of file EvtDalitzReso.hh.

               { NBW            = 0 , RBW_ZEMACH        = 1 , RBW_KUEHN  = 2 , RBW_CLEO        = 3 ,
         RBW_ZEMACH2    = 4 , GS_CLEO           = 5 , K_MATRIX   = 6 , RBW_CLEO_ZEMACH = 7 ,
         GS_CLEO_ZEMACH = 8 , LASS              = 9 , K_MATRIX_I = 10, K_MATRIX_II     = 11,
         GAUSS_CLEO     = 12, GAUSS_CLEO_ZEMACH = 13 };

Constructor & Destructor Documentation

EvtDalitzReso() [1/7]

EvtDalitzReso::EvtDalitzReso ( )

inline

Definition at line 53 of file EvtDalitzReso.hh.

{};

Referenced by clone().

EvtDalitzReso() [2/7]

EvtDalitzReso::EvtDalitzReso	(	const EvtDalitzPlot &	dp,
		EvtCyclic3::Pair	pairAng,
		EvtCyclic3::Pair	pairRes,
		EvtSpinType::spintype	spin,
		double	m0,
		double	g0,
		NumType	typeN )

Definition at line 38 of file EvtDalitzReso.cc.

  : _dp(dp),
    _pairAng(pairAng),
    _pairRes(pairRes),
    _spin(spin),
    _typeN(typeN),
    _m0(m0),_g0(g0),
    _massFirst(dp.m(first(pairRes))),_massSecond(dp.m(second(pairRes))),
    _m0_mix(-1.),_g0_mix(0.),_delta_mix(0.),_amp_mix(0.,0.),
    _g1(-1.),_g2(-1.),_coupling2(Undefined),
    _kmatrix_index(-1),_fr12prod(0.,0.),_fr13prod(0.,0.),_fr14prod(0.,0.),_fr15prod(0.,0.),_s0prod(0.),
    _a(0.),_r(0.),_Blass(0.),_phiB(0.),_R(0.),_phiR(0.)
{
  _vb = EvtTwoBodyVertex(_m0,_dp.m(EvtCyclic3::other(_pairRes)),_dp.bigM(),_spin); 
  _vd = EvtTwoBodyVertex(_massFirst,_massSecond,_m0,_spin);
  _vb.set_f( 0.0 ); // Default values for Blatt-Weisskopf factors.
  _vd.set_f( 1.5 );
  assert(_typeN != K_MATRIX && _typeN != K_MATRIX_I && _typeN != K_MATRIX_II);  // single BW cannot be K-matrix
}

EvtDalitzReso() [3/7]

EvtDalitzReso::EvtDalitzReso	(	const EvtDalitzPlot &	dp,
		EvtCyclic3::Pair	pairAng,
		EvtCyclic3::Pair	pairRes,
		EvtSpinType::spintype	spin,
		double	m0,
		double	g0,
		NumType	typeN,
		double	m0_mix,
		double	g0_mix,
		double	delta_mix,
		EvtComplex	amp_mix )

Definition at line 61 of file EvtDalitzReso.cc.

  : _dp(dp),
    _pairAng(pairAng),
    _pairRes(pairRes),
    _spin(spin),
    _typeN(typeN),
    _m0(m0),_g0(g0),
    _massFirst(dp.m(first(pairRes))),_massSecond(dp.m(second(pairRes))),
    _m0_mix(m0_mix),_g0_mix(g0_mix),_delta_mix(delta_mix),_amp_mix(amp_mix),
    _g1(-1.),_g2(-1.),_coupling2(Undefined),
    _kmatrix_index(-1),_fr12prod(0.,0.),_fr13prod(0.,0.),_fr14prod(0.,0.),_fr15prod(0.,0.),_s0prod(0.),
    _a(0.),_r(0.),_Blass(0.),_phiB(0.),_R(0.),_phiR(0.)
{
  _vb = EvtTwoBodyVertex(_m0,_dp.m(EvtCyclic3::other(_pairRes)),_dp.bigM(),_spin); 
  _vd = EvtTwoBodyVertex(_massFirst,_massSecond,_m0,_spin);
  _vb.set_f( 0.0 ); // Default values for Blatt-Weisskopf factors.
  _vd.set_f( 1.5 );
  // single BW (with electromagnetic mixing) cannot be K-matrix
  assert(_typeN != K_MATRIX && _typeN != K_MATRIX_I && _typeN != K_MATRIX_II);
}

EvtDalitzReso() [4/7]

EvtDalitzReso::EvtDalitzReso	(	const EvtDalitzPlot &	dp,
		EvtCyclic3::Pair	pairAng,
		EvtCyclic3::Pair	pairRes,
		EvtSpinType::spintype	spin,
		double	m0,
		NumType	typeN,
		double	g1,
		double	g2,
		CouplingType	coupling2 )

Definition at line 85 of file EvtDalitzReso.cc.

  : _dp(dp),
    _pairAng(pairAng),
    _pairRes(pairRes),
    _spin(spin),
    _typeN(typeN),
    _m0(m0),_g0(-1.),
    _massFirst(dp.m(first(pairRes))),_massSecond(dp.m(second(pairRes))),
    _m0_mix(-1.),_g0_mix(0.),_delta_mix(0.),_amp_mix(0.,0.),
    _g1(g1),_g2(g2),_coupling2(coupling2),
    _kmatrix_index(-1),_fr12prod(0.,0.),_fr13prod(0.,0.),_fr14prod(0.,0.),_fr15prod(0.,0.),_s0prod(0.),
    _a(0.),_r(0.),_Blass(0.),_phiB(0.),_R(0.),_phiR(0.)
{
  _vb = EvtTwoBodyVertex(_m0,_dp.m(EvtCyclic3::other(_pairRes)),_dp.bigM(),_spin);   
  _vd = EvtTwoBodyVertex(_massFirst,_massSecond,_m0,_spin);
  _vb.set_f( 0.0 ); // Default values for Blatt-Weisskopf factors.
  _vd.set_f( 1.5 );
  assert(_coupling2 != Undefined);
  assert(_typeN != K_MATRIX && _typeN != K_MATRIX_I && _typeN != K_MATRIX_II); // coupled BW cannot be K-matrix
  assert(_typeN != LASS);     // coupled BW cannot be LASS
  assert(_typeN != NBW);      // for coupled BW, only relativistic BW 
}

EvtDalitzReso() [5/7]

EvtDalitzReso::EvtDalitzReso	(	const EvtDalitzPlot &	dp,
		EvtCyclic3::Pair	pairRes,
		std::string	nameIndex,
		NumType	typeN,
		EvtComplex	fr12prod,
		EvtComplex	fr13prod,
		EvtComplex	fr14prod,
		EvtComplex	fr15prod,
		double	s0prod )

Definition at line 111 of file EvtDalitzReso.cc.

  : _dp(dp),
    _pairRes(pairRes),
    _typeN(typeN),
    _m0(0.),_g0(0.),
    _massFirst(dp.m(first(pairRes))),_massSecond(dp.m(second(pairRes))),
    _m0_mix(-1.),_g0_mix(0.),_delta_mix(0.),_amp_mix(0.,0.),
    _g1(-1.),_g2(-1.),_coupling2(Undefined),
    _kmatrix_index(-1),_fr12prod(fr12prod),_fr13prod(fr13prod),_fr14prod(fr14prod),_fr15prod(fr15prod),_s0prod(s0prod),
    _a(0.),_r(0.),_Blass(0.),_phiB(0.),_R(0.),_phiR(0.)
{
  assert(_typeN==K_MATRIX || _typeN==K_MATRIX_I || _typeN==K_MATRIX_II);
  _spin=EvtSpinType::SCALAR;
  if (nameIndex=="Pole1") _kmatrix_index=1;
  else if (nameIndex=="Pole2") _kmatrix_index=2;
  else if (nameIndex=="Pole3") _kmatrix_index=3;
  else if (nameIndex=="Pole4") _kmatrix_index=4;
  else if (nameIndex=="Pole5") _kmatrix_index=5;
  else if (nameIndex=="f11prod") _kmatrix_index=6;
  else assert(0);
}

EvtDalitzReso() [6/7]

EvtDalitzReso::EvtDalitzReso	(	const EvtDalitzPlot &	dp,
		EvtCyclic3::Pair	pairRes,
		double	m0,
		double	g0,
		double	a,
		double	r,
		double	B,
		double	phiB,
		double	R,
		double	phiR )

Definition at line 136 of file EvtDalitzReso.cc.

  : _dp(dp),
    _pairRes(pairRes),
    _typeN(LASS),
    _m0(m0),_g0(g0),
    _massFirst(dp.m(first(pairRes))),_massSecond(dp.m(second(pairRes))),
    _m0_mix(-1.),_g0_mix(0.),_delta_mix(0.),_amp_mix(0.,0.),
    _g1(-1.),_g2(-1.),_coupling2(Undefined),
    _kmatrix_index(-1),_fr12prod(0.,0.),_fr13prod(0.,0.),_fr14prod(0.,0.),_fr15prod(0.,0.),_s0prod(0.),
    _a(a),_r(r),_Blass(B),_phiB(phiB),_R(R),_phiR(phiR)
{
  _spin=EvtSpinType::SCALAR;
  _vd = EvtTwoBodyVertex(_massFirst,_massSecond,_m0,_spin);
  _vd.set_f( 1.5 ); // Default values for Blatt-Weisskopf factors.
}

EvtDalitzReso() [7/7]

EvtDalitzReso::EvtDalitzReso

(

const EvtDalitzReso &

other

)

Definition at line 155 of file EvtDalitzReso.cc.

  : _dp(other._dp),
    _pairAng(other._pairAng),
    _pairRes(other._pairRes),
    _spin(other._spin),
    _typeN(other._typeN),
    _m0(other._m0),_g0(other._g0),
    _vb(other._vb),_vd(other._vd),
    _massFirst(other._massFirst),_massSecond(other._massSecond),
    _m0_mix(other._m0_mix),_g0_mix(other._g0_mix),_delta_mix(other._delta_mix),_amp_mix(other._amp_mix),
    _g1(other._g1),_g2(other._g2),_coupling2(other._coupling2),
    _kmatrix_index(other._kmatrix_index),
    _fr12prod(other._fr12prod),_fr13prod(other._fr13prod),_fr14prod(other._fr14prod),_fr15prod(other._fr15prod),
    _s0prod(other._s0prod),
    _a(other._a),_r(other._r),_Blass(other._Blass),_phiB(other._phiB),_R(other._R),_phiR(other._phiR)
{}

~EvtDalitzReso()

EvtDalitzReso::~EvtDalitzReso

(

)

Definition at line 173 of file EvtDalitzReso.cc.

{}

Member Function Documentation

clone()

EvtDalitzReso * EvtDalitzReso::clone ( ) const

inline

Definition at line 77 of file EvtDalitzReso.hh.

{ return new EvtDalitzReso(*this); }

evaluate()

EvtComplex EvtDalitzReso::evaluate

(

const EvtDalitzPoint &

p

)

Definition at line 177 of file EvtDalitzReso.cc.

{
  double m = sqrt(x.q(_pairRes));
 
  // do use always hash table (speed up fitting)
  if (_typeN==K_MATRIX || _typeN==K_MATRIX_I || _typeN==K_MATRIX_II)
    return Fvector( m*m, _kmatrix_index );
 
  if (_typeN==LASS)
    return lass(m*m);
 
  EvtComplex amp(1.0,0.0);
 
  if (_dp.bigM() != x.bigM()) _vb = EvtTwoBodyVertex(_m0,_dp.m(EvtCyclic3::other(_pairRes)),x.bigM(),_spin); 
  EvtTwoBodyKine vb(m,x.m(EvtCyclic3::other(_pairRes)),x.bigM());
  EvtTwoBodyKine vd(_massFirst,_massSecond,m);   
 
  EvtComplex prop(0,0);
  if (_typeN==NBW) {
    prop = propBreitWigner(_m0,_g0,m);
  } else if (_typeN==GAUSS_CLEO || _typeN==GAUSS_CLEO_ZEMACH) {
    prop = propGauss(_m0,_g0,m);
  } else {
    if (_coupling2==Undefined) {  
      // single BW
      double g = (_g0<=0. || _vd.pD()<=0.)? -_g0 : _g0*_vd.widthFactor(vd);  // running width
      if (_typeN==GS_CLEO || _typeN==GS_CLEO_ZEMACH) {
    // Gounaris-Sakurai (GS)
    assert(_massFirst==_massSecond);
    prop = propGounarisSakurai(_m0,fabs(_g0),_vd.pD(),m,g,vd.p());
      } else {
    // standard relativistic BW
    prop = propBreitWignerRel(_m0,g,m);
      }
    } else {    
      // coupled width BW
      EvtComplex G1,G2;
      switch (_coupling2) { 
      case PicPic: {
    G1 = _g1*_g1*psFactor(_massFirst,_massSecond,m);
    static double mPic = EvtPDL::getMass( EvtPDL::getId( "pi+" ) );
    G2 = _g2*_g2*psFactor(mPic,mPic,m);
    break;
      }
      case PizPiz: {
    G1 = _g1*_g1*psFactor(_massFirst,_massSecond,m);
    static double mPiz = EvtPDL::getMass( EvtPDL::getId( "pi0" ) );
    G2 = _g2*_g2*psFactor(mPiz,mPiz,m);
    break;
      }
      case PiPi: {
    G1 = _g1*_g1*psFactor(_massFirst,_massSecond,m);
    static double mPic = EvtPDL::getMass( EvtPDL::getId( "pi+" ) );
    static double mPiz = EvtPDL::getMass( EvtPDL::getId( "pi0" ) );
    G2 = _g2*_g2*psFactor(mPic,mPic,mPiz,mPiz,m);
    break;
      }
      case KcKc: {
    G1 = _g1*_g1*psFactor(_massFirst,_massSecond,m);
    static double mKc = EvtPDL::getMass( EvtPDL::getId( "K+" ) );
    G2 = _g2*_g2*psFactor(mKc,mKc,m);
    break;
      }
      case KzKz: {
    G1 = _g1*_g1*psFactor(_massFirst,_massSecond,m);
    static double mKz = EvtPDL::getMass( EvtPDL::getId( "K0" ) );
    G2 = _g2*_g2*psFactor(mKz,mKz,m);
    break;
      }
      case KK: {
    G1 = _g1*_g1*psFactor(_massFirst,_massSecond,m);
    static double mKc = EvtPDL::getMass( EvtPDL::getId( "K+" ) );
    static double mKz = EvtPDL::getMass( EvtPDL::getId( "K0" ) );
    G2 = _g2*_g2*psFactor(mKc,mKc,mKz,mKz,m);
    break;
      }
      case EtaPic: {
    G1 = _g1*_g1*psFactor(_massFirst,_massSecond,m);
    static double mEta = EvtPDL::getMass( EvtPDL::getId( "eta" ) );
    static double mPic = EvtPDL::getMass( EvtPDL::getId( "pi+" ) );
    G2 = _g2*_g2*psFactor(mEta,mPic,m);
    break;
      }
      case EtaPiz: {
    G1 = _g1*_g1*psFactor(_massFirst,_massSecond,m);
    static double mEta = EvtPDL::getMass( EvtPDL::getId( "eta" ) );
    static double mPiz = EvtPDL::getMass( EvtPDL::getId( "pi0" ) );
    G2 = _g2*_g2*psFactor(mEta,mPiz,m);
    break;
      }
      case PicPicKK: {
    static double mPic = EvtPDL::getMass( EvtPDL::getId( "pi+" ) );
    //G1 = _g1*_g1*psFactor(mPic,mPic,m);
    G1 = _g1*psFactor(mPic,mPic,m);
    static double mKc = EvtPDL::getMass( EvtPDL::getId( "K+" ) );
    static double mKz = EvtPDL::getMass( EvtPDL::getId( "K0" ) );
    //G2 = _g2*_g2*psFactor(mKc,mKc,mKz,mKz,m);
    G2 = _g2*psFactor(mKc,mKc,mKz,mKz,m);
    break;
      }
      default:
    std::cout << "EvtDalitzReso:evaluate(): PANIC, wrong coupling2 state." << std::endl;
    assert(0);
    break;
      }
      // calculate standard couple BW propagator
      if (_coupling2 != WA76)
    prop = _g1*propBreitWignerRelCoupled(_m0,G1,G2,m);
    } 
  }
  amp *= prop;
 
  // Compute form-factors (Blatt-Weisskopf penetration factor)
  amp *= _vb.formFactor(vb);  
  amp *= _vd.formFactor(vd);  
 
  // Compute numerator (angular distribution)
  amp *= numerator(x,vb,vd);  
 
  // Compute electromagnetic mass mixing factor
  if (_m0_mix>0.) {
    EvtComplex prop_mix;
    if (_typeN==NBW) {
      prop_mix = propBreitWigner(_m0_mix,_g0_mix,m);
    } else {
      assert(_g1<0.); // running width only
      double g_mix = _g0_mix*_vd.widthFactor(vd);
      prop_mix = propBreitWignerRel(_m0_mix,g_mix,m);
    }
    amp *= mixFactor(prop,prop_mix);
  }
 
  return amp;
}

set_fb()

void EvtDalitzReso::set_fb ( double R )

inline

Definition at line 82 of file EvtDalitzReso.hh.

{ _vb.set_f( R ); }

set_fd()

void EvtDalitzReso::set_fd ( double R )

inline

Definition at line 81 of file EvtDalitzReso.hh.

{ _vd.set_f( R ); }

---

The documentation for this class was generated from the following files:

• EvtDalitzReso.hh
• EvtDalitzReso.cc