EvtDalitzPlot Class Reference

#include <EvtDalitzPlot.hh>

Public Member Functions
	EvtDalitzPlot ()
	EvtDalitzPlot (double mA, double mB, double mC, double bigM, double ldel=0., double rdel=0.)
	EvtDalitzPlot (const EvtDecayMode &mode, double ldel=0., double rdel=0.)
	EvtDalitzPlot (const EvtDalitzPlot &other)
	~EvtDalitzPlot ()
bool	operator== (const EvtDalitzPlot &other) const
const EvtDalitzPlot *	clone () const
double	qAbsMin (EvtCyclic3::Pair i) const
double	qAbsMax (EvtCyclic3::Pair i) const
double	mAbsMin (EvtCyclic3::Pair i) const
double	mAbsMax (EvtCyclic3::Pair i) const
double	qResAbsMin (EvtCyclic3::Pair i) const
double	qResAbsMax (EvtCyclic3::Pair i) const
double	qHelAbsMin (EvtCyclic3::Pair i) const
double	qHelAbsMax (EvtCyclic3::Pair i) const
double	qSumMin () const
double	qSumMax () const
bool	fuzzy () const
double	getArea (int N=1000, EvtCyclic3::Pair i=EvtCyclic3::AB, EvtCyclic3::Pair j=EvtCyclic3::BC) const
double	qMin (EvtCyclic3::Pair i, EvtCyclic3::Pair j, double q) const
double	qMax (EvtCyclic3::Pair i, EvtCyclic3::Pair j, double q) const
double	cosTh (EvtCyclic3::Pair i1, double q1, EvtCyclic3::Pair i2, double q2) const
double	e (EvtCyclic3::Index i, EvtCyclic3::Pair j, double q) const
double	p (EvtCyclic3::Index i, EvtCyclic3::Pair j, double q) const
double	q (EvtCyclic3::Pair i1, double cosTh, EvtCyclic3::Pair i2, double q2) const
double	jacobian (EvtCyclic3::Pair i, double q) const
EvtTwoBodyVertex	vD (EvtCyclic3::Pair iRes, double m0, int L) const
EvtTwoBodyVertex	vB (EvtCyclic3::Pair iRes, double m0, int L) const
double	sum () const
double	bigM () const
double	mA () const
double	mB () const
double	mC () const
double	m (EvtCyclic3::Index i) const
void	print () const
void	sanityCheck () const

Protected Attributes
double	_mA
double	_mB
double	_mC
double	_bigM
double	_ldel
double	_rdel

Detailed Description

Definition at line 25 of file EvtDalitzPlot.hh.

Constructor & Destructor Documentation

EvtDalitzPlot() [1/4]

EvtDalitzPlot::EvtDalitzPlot

(

)

Definition at line 33 of file EvtDalitzPlot.cc.

  : _mA(0.), _mB(0.), _mC(0.), _bigM(0.),
    _ldel(0.), _rdel(0.)
{}

Referenced by clone().

EvtDalitzPlot() [2/4]

EvtDalitzPlot::EvtDalitzPlot	(	double	mA,
		double	mB,
		double	mC,
		double	bigM,
		double	ldel = 0.,
		double	rdel = 0. )

Definition at line 39 of file EvtDalitzPlot.cc.

  : _mA(mA), _mB(mB), _mC(mC), _bigM(bigM),
    _ldel(ldel), _rdel(rdel) 
{
  sanityCheck();
}

EvtDalitzPlot() [3/4]

EvtDalitzPlot::EvtDalitzPlot	(	const EvtDecayMode &	mode,
		double	ldel = 0.,
		double	rdel = 0. )

Definition at line 47 of file EvtDalitzPlot.cc.

{
  _mA = EvtPDL::getMeanMass(EvtPDL::getId(mode.dau(A)));
  _mB = EvtPDL::getMeanMass(EvtPDL::getId(mode.dau(B)));
  _mC = EvtPDL::getMeanMass(EvtPDL::getId(mode.dau(C)));
  _bigM = EvtPDL::getMeanMass(EvtPDL::getId(mode.mother()));
 
  _ldel = ldel;
  _rdel = rdel;
  
  sanityCheck();
}

EvtDalitzPlot() [4/4]

EvtDalitzPlot::EvtDalitzPlot

(

const EvtDalitzPlot &

other

)

Definition at line 61 of file EvtDalitzPlot.cc.

  : _mA(other._mA), _mB(other._mB), _mC(other._mC), _bigM(other._bigM),
    _ldel(other._ldel), _rdel(other._rdel)
{}

~EvtDalitzPlot()

EvtDalitzPlot::~EvtDalitzPlot

(

)

Definition at line 67 of file EvtDalitzPlot.cc.

{}

Member Function Documentation

bigM()

double EvtDalitzPlot::bigM ( ) const

inline

Definition at line 92 of file EvtDalitzPlot.hh.

{ return _bigM; }

Referenced by e(), EvtDalitzReso::evaluate(), EvtDalitzReso::EvtDalitzReso(), EvtDalitzReso::EvtDalitzReso(), EvtDalitzReso::EvtDalitzReso(), print(), qMax(), qMin(), and vB().

clone()

const EvtDalitzPlot * EvtDalitzPlot::clone

(

)

const

Definition at line 83 of file EvtDalitzPlot.cc.

{
  return new EvtDalitzPlot(*this);
}

cosTh()

double EvtDalitzPlot::cosTh	(	EvtCyclic3::Pair	i1,
		double	q1,
		EvtCyclic3::Pair	i2,
		double	q2 ) const

Definition at line 260 of file EvtDalitzPlot.cc.

{
  if(i1 == i2) return 1.;
  
  double qmax = qMax(i1,i2,q2);
  double qmin = qMin(i1,i2,q2);
  
  double cos = (qmax + qmin - 2*q1)/(qmax - qmin);
  
  return cos;
}

Referenced by EvtDalitzPoint::cosTh(), and q().

e()

double EvtDalitzPlot::e	(	EvtCyclic3::Index	i,
		EvtCyclic3::Pair	j,
		double	q ) const

Definition at line 273 of file EvtDalitzPlot.cc.

{
  if(i == other(j)) {
 
    // i does not belong to pair j
 
    return (bigM()*bigM()-q-m(i)*m(i))/2/sqrt(q);
  }
  else {
    
    // i and k make pair j
 
    Index k;
    if(first(j) == i) k = second(j);
    else k = first(j); 
 
    double e = (q + m(i)*m(i) - m(k)*m(k))/2/sqrt(q);          
    return e;
  }
}

Referenced by e(), EvtDalitzPoint::e(), p(), and q().

fuzzy()

bool EvtDalitzPlot::fuzzy ( ) const

inline

Definition at line 55 of file EvtDalitzPlot.hh.

{ return (_rdel - _ldel != 0.); }

getArea()

double EvtDalitzPlot::getArea	(	int	N = 1000,
		EvtCyclic3::Pair	i = EvtCyclic3::AB,
		EvtCyclic3::Pair	j = EvtCyclic3::BC ) const

Definition at line 240 of file EvtDalitzPlot.cc.

{
  // Trapezoidal integral over qi. qj can be calculated.
  // The first and the last point are zero, so they are not counted
 
  double dh = (qAbsMax(i) - qAbsMin(i))/((double) N);
  double sum = 0;
 
  int ii;
  for(ii=1;ii<N;ii++) {
 
    double x = qAbsMin(i) + ii*dh;
    double dy = qMax(j,i,x) - qMin(j,i,x);
    sum += dy;
  }
 
  return sum * dh;
}

Referenced by EvtDalitzFlatPdf::compute_integral().

jacobian()

double EvtDalitzPlot::jacobian	(	EvtCyclic3::Pair	i,
		double	q ) const

Definition at line 319 of file EvtDalitzPlot.cc.

{
  return 2*p(first(i),i,q)*p(other(i),i,q);  // J(BC) = 2pA*pB = 2pA*pC
}

m()

double EvtDalitzPlot::m

(

EvtCyclic3::Index

i

)

const

Definition at line 101 of file EvtDalitzPlot.cc.

                                     {
 
  double m = _mA;
  if(i == B) m = _mB;
  else
    if(i == C) m = _mC;
 
  return m;
}

Referenced by e(), EvtDalitzReso::evaluate(), EvtDalitzReso::EvtDalitzReso(), EvtDalitzReso::EvtDalitzReso(), EvtDalitzReso::EvtDalitzReso(), m(), p(), q(), qAbsMax(), qAbsMin(), qMax(), qMin(), vB(), and vD().

mA()

double EvtDalitzPlot::mA ( ) const

inline

Definition at line 93 of file EvtDalitzPlot.hh.

{ return _mA; } 

mAbsMax()

double EvtDalitzPlot::mAbsMax

(

EvtCyclic3::Pair

i

)

const

Definition at line 165 of file EvtDalitzPlot.cc.

{
  return sqrt(qAbsMax(i));
}

mAbsMin()

double EvtDalitzPlot::mAbsMin

(

EvtCyclic3::Pair

i

)

const

Definition at line 159 of file EvtDalitzPlot.cc.

{
  return sqrt(qAbsMin(i));
}

mB()

double EvtDalitzPlot::mB ( ) const

inline

Definition at line 94 of file EvtDalitzPlot.hh.

{ return _mB; } 

mC()

double EvtDalitzPlot::mC ( ) const

inline

Definition at line 95 of file EvtDalitzPlot.hh.

{ return _mC; } 

operator==()

bool EvtDalitzPlot::operator==

(

const EvtDalitzPlot &

other

)

const

Definition at line 71 of file EvtDalitzPlot.cc.

{
  bool ret = false;
  if(_mA == other._mA && 
     _mB == other._mB &&
     _mC == other._mC &&
     _bigM == other._bigM) ret = true;
 
  return ret;
}

p()

double EvtDalitzPlot::p	(	EvtCyclic3::Index	i,
		EvtCyclic3::Pair	j,
		double	q ) const

Definition at line 295 of file EvtDalitzPlot.cc.

{
  double en = e(i,j,q);
  double p2 = en*en - m(i)*m(i);
  
  if(p2 < 0) {
    printf("Bad value of p2 %f %d %d %f %f\n",p2,i,j,en,m(i));
    assert(0);
  }
 
  return sqrt(p2);  
}

Referenced by jacobian(), EvtDalitzPoint::p(), and q().

print()

void EvtDalitzPlot::print

(

)

const

Definition at line 338 of file EvtDalitzPlot.cc.

{
  // masses
  printf("Mass  M    %f\n",bigM());
  printf("Mass mA    %f\n",_mA);
  printf("Mass mB    %f\n",_mB);
  printf("Mass mC    %f\n",_mC);
  // limits
  printf("Limits qAB %f : %f\n",qAbsMin(AB),qAbsMax(AB));
  printf("Limits qBC %f : %f\n",qAbsMin(BC),qAbsMax(BC));
  printf("Limits qCA %f : %f\n",qAbsMin(CA),qAbsMax(CA));
  printf("Sum q       %f\n",sum());
  printf("Limit qsum  %f : %f\n",qSumMin(),qSumMax());
}

Referenced by EvtDalitzPoint::print().

q()

double EvtDalitzPlot::q	(	EvtCyclic3::Pair	i1,
		double	cosTh,
		EvtCyclic3::Pair	i2,
		double	q2 ) const

Definition at line 309 of file EvtDalitzPlot.cc.

{
  if(i1 == i2) return q2;
 
  EvtCyclic3::Index f = first(i1);
  EvtCyclic3::Index s = second(i1);
  return m(f)*m(f) + m(s)*m(s) + 2*e(f,i2,q2)*e(s,i2,q2) - 2*p(f,i2,q2)*p(s,i2,q2)*cosTh;
}

qAbsMax()

double EvtDalitzPlot::qAbsMax

(

EvtCyclic3::Pair

i

)

const

Definition at line 127 of file EvtDalitzPlot.cc.

{
  Index j = other(i);
  return (_bigM-m(j))*(_bigM-m(j));
}

Referenced by EvtDalitzResPdf::compute_integral(), getArea(), EvtDalitzPoint::isValid(), mAbsMax(), print(), qHelAbsMax(), qHelAbsMin(), qResAbsMax(), EvtDalitzFlatPdf::randomPoint(), and EvtDalitzResPdf::randomPoint().

qAbsMin()

double EvtDalitzPlot::qAbsMin

(

EvtCyclic3::Pair

i

)

const

Definition at line 118 of file EvtDalitzPlot.cc.

{
  Index j = first(i);
  Index k = second(i);
 
  return (m(j) + m(k))*(m(j) + m(k));
}

Referenced by EvtDalitzResPdf::compute_integral(), getArea(), EvtDalitzPoint::isValid(), mAbsMin(), print(), qHelAbsMax(), qHelAbsMin(), qResAbsMin(), EvtDalitzFlatPdf::randomPoint(), and EvtDalitzResPdf::randomPoint().

qHelAbsMax()

double EvtDalitzPlot::qHelAbsMax

(

EvtCyclic3::Pair

i

)

const

Definition at line 151 of file EvtDalitzPlot.cc.

{
  Pair j = next(i);
  Pair k = prev(i);
  return  (qAbsMax(j) - qAbsMin(k))/2.;
}

qHelAbsMin()

double EvtDalitzPlot::qHelAbsMin

(

EvtCyclic3::Pair

i

)

const

Definition at line 144 of file EvtDalitzPlot.cc.

{
  Pair j = next(i);
  Pair k = prev(i);
  return  (qAbsMin(j) - qAbsMax(k))/2.;
}

qMax()

double EvtDalitzPlot::qMax	(	EvtCyclic3::Pair	i,
		EvtCyclic3::Pair	j,
		double	q ) const

Definition at line 208 of file EvtDalitzPlot.cc.

{
 
  if(i == j) return q;
  else {
 
    // Particle pair j defines the rest-frame
    // 0 - particle common to r.f. and angle calculations
    // 1 - particle belonging to r.f. but not angle
    // 2 - particle not belonging to r.f.
    
    Index k0 = common(i,j);
    Index k2 = other(j);
    Index k1 = other(k0,k2); 
 
    // Energy, momentum of particle common to rest-frame and angle
    EvtTwoBodyKine jpair(m(k0),m(k1),sqrt(q)); 
    double ek = jpair.e(EvtTwoBodyKine::A,EvtTwoBodyKine::AB);
    double pk = jpair.p();
 
    // Energy and momentum of the other particle
    EvtTwoBodyKine mother(sqrt(q),m(k2),bigM());
    double ej = mother.e(EvtTwoBodyKine::B,EvtTwoBodyKine::A);
    double pj = mother.p(EvtTwoBodyKine::A);
 
    
    // See PDG 34.4.3.1
    return (ek+ej)*(ek+ej) - (pk-pj)*(pk-pj);
  }
}

Referenced by EvtDalitzResPdf::compute_integral(), cosTh(), getArea(), and EvtDalitzPoint::qMax().

qMin()

double EvtDalitzPlot::qMin	(	EvtCyclic3::Pair	i,
		EvtCyclic3::Pair	j,
		double	q ) const

Definition at line 173 of file EvtDalitzPlot.cc.

{
  if(i == j) return q;
 
  else {
 
    // Particle pair j defines the rest-frame
    // 0 - particle common to r.f. and angle calculations
    // 1 - particle belonging to r.f. but not angle
    // 2 - particle not belonging to r.f.
 
    Index k0 = common(i,j);
    Index k2 = other(j);
    Index k1 = other(k0,k2);
 
    // Energy, momentum of particle common to rest-frame and angle
    EvtTwoBodyKine jpair(m(k0),m(k1),sqrt(q)); 
    double pk = jpair.p();
    double ek = jpair.e(EvtTwoBodyKine::A,EvtTwoBodyKine::AB);
 
 
    // Energy and momentum of the other particle
    EvtTwoBodyKine mother(sqrt(q),m(k2),bigM());
    double ej = mother.e(EvtTwoBodyKine::B,EvtTwoBodyKine::A);
    double pj = mother.p(EvtTwoBodyKine::A);
 
 
    // See PDG 34.4.3.1
    return (ek+ej)*(ek+ej) - (pk+pj)*(pk+pj);
  }
}

Referenced by EvtDalitzResPdf::compute_integral(), cosTh(), getArea(), and EvtDalitzPoint::qMin().

qResAbsMax()

double EvtDalitzPlot::qResAbsMax

(

EvtCyclic3::Pair

i

)

const

Definition at line 139 of file EvtDalitzPlot.cc.

{
  return  qAbsMax(i) - sum()/3.;
}

qResAbsMin()

double EvtDalitzPlot::qResAbsMin

(

EvtCyclic3::Pair

i

)

const

Definition at line 134 of file EvtDalitzPlot.cc.

{
  return qAbsMin(i) - sum()/3.;
}

qSumMax()

double EvtDalitzPlot::qSumMax ( ) const

inline

Definition at line 54 of file EvtDalitzPlot.hh.

{ return sum() + _rdel; }

Referenced by print().

qSumMin()

double EvtDalitzPlot::qSumMin ( ) const

inline

Definition at line 53 of file EvtDalitzPlot.hh.

{ return sum() + _ldel; }

Referenced by print().

sanityCheck()

void EvtDalitzPlot::sanityCheck

(

)

const

Definition at line 89 of file EvtDalitzPlot.cc.

{
  if(_mA < 0 || _mB < 0 || _mC < 0 || _bigM <= 0 || _bigM - _mA - _mB - _mC < 0.) {
    
    printf("Invalid Dalitz plot %f %f %f %f\n",_mA,_mB,_mC,_bigM);
    assert(0);
  }
  assert(_ldel <= 0.);
  assert(_rdel >= 0.);
}

Referenced by EvtDalitzPlot(), and EvtDalitzPlot().

sum()

double EvtDalitzPlot::sum

(

)

const

Definition at line 112 of file EvtDalitzPlot.cc.

{
  return _mA*_mA + _mB*_mB + _mC*_mC + _bigM*_bigM;
}

Referenced by EvtDalitzPoint::EvtDalitzPoint(), getArea(), print(), qResAbsMax(), qResAbsMin(), qSumMax(), and qSumMin().

vB()

EvtTwoBodyVertex EvtDalitzPlot::vB	(	EvtCyclic3::Pair	iRes,
		double	m0,
		int	L ) const

Definition at line 332 of file EvtDalitzPlot.cc.

{
  return EvtTwoBodyVertex(m0,m(other(iRes)),bigM(),L);
}

vD()

EvtTwoBodyVertex EvtDalitzPlot::vD	(	EvtCyclic3::Pair	iRes,
		double	m0,
		int	L ) const

Definition at line 325 of file EvtDalitzPlot.cc.

{
  return EvtTwoBodyVertex(m(first(iRes)),
              m(second(iRes)),m0,L);
}

Member Data Documentation

_bigM

double EvtDalitzPlot::_bigM

protected

Definition at line 110 of file EvtDalitzPlot.hh.

Referenced by bigM(), EvtDalitzPlot(), operator==(), qAbsMax(), sanityCheck(), and sum().

_ldel

double EvtDalitzPlot::_ldel

protected

Definition at line 114 of file EvtDalitzPlot.hh.

Referenced by EvtDalitzPlot(), fuzzy(), qSumMin(), and sanityCheck().

_mA

double EvtDalitzPlot::_mA

protected

Definition at line 107 of file EvtDalitzPlot.hh.

Referenced by EvtDalitzPlot(), m(), mA(), operator==(), print(), sanityCheck(), and sum().

_mB

double EvtDalitzPlot::_mB

protected

Definition at line 108 of file EvtDalitzPlot.hh.

Referenced by EvtDalitzPlot(), m(), mB(), operator==(), print(), sanityCheck(), and sum().

_mC

double EvtDalitzPlot::_mC

protected

Definition at line 109 of file EvtDalitzPlot.hh.

Referenced by EvtDalitzPlot(), m(), mC(), operator==(), print(), sanityCheck(), and sum().

_rdel

double EvtDalitzPlot::_rdel

protected

Definition at line 115 of file EvtDalitzPlot.hh.

Referenced by EvtDalitzPlot(), fuzzy(), qSumMax(), and sanityCheck().

---

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

• EvtDalitzPlot.hh
• EvtDalitzPlot.cc