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 32 of file EvtDalitzPlot.cc.

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

Referenced by clone(), EvtDalitzPlot(), and operator==().

EvtDalitzPlot() [2/4]

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

Definition at line 35 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 41 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 53 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 61 of file EvtDalitzPlot.cc.

{}

Member Function Documentation

bigM()

double EvtDalitzPlot::bigM ( ) const

inline

Definition at line 90 of file EvtDalitzPlot.hh.

{ return _bigM; }

Referenced by e(), EvtDalitzPlot(), print(), qMax(), qMin(), and vB().

clone()

const EvtDalitzPlot * EvtDalitzPlot::clone

(

)

const

Definition at line 71 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 214 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 226 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 54 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 195 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;
}

jacobian()

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

Definition at line 271 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 84 of file EvtDalitzPlot.cc.

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

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

mA()

double EvtDalitzPlot::mA ( ) const

inline

Definition at line 91 of file EvtDalitzPlot.hh.

{ return _mA; }

Referenced by EvtDalitzPlot().

mAbsMax()

double EvtDalitzPlot::mAbsMax

(

EvtCyclic3::Pair

i

)

const

Definition at line 129 of file EvtDalitzPlot.cc.

{ return sqrt( qAbsMax( i ) ); }

mAbsMin()

double EvtDalitzPlot::mAbsMin

(

EvtCyclic3::Pair

i

)

const

Definition at line 127 of file EvtDalitzPlot.cc.

{ return sqrt( qAbsMin( i ) ); }

mB()

double EvtDalitzPlot::mB ( ) const

inline

Definition at line 92 of file EvtDalitzPlot.hh.

{ return _mB; }

Referenced by EvtDalitzPlot().

mC()

double EvtDalitzPlot::mC ( ) const

inline

Definition at line 93 of file EvtDalitzPlot.hh.

{ return _mC; }

Referenced by EvtDalitzPlot().

operator==()

bool EvtDalitzPlot::operator==

(

const EvtDalitzPlot &

other

)

const

Definition at line 63 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 248 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 283 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 261 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;
}

Referenced by e(), jacobian(), p(), qMax(), and qMin().

qAbsMax()

double EvtDalitzPlot::qAbsMax

(

EvtCyclic3::Pair

i

)

const

Definition at line 102 of file EvtDalitzPlot.cc.

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

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

qAbsMin()

double EvtDalitzPlot::qAbsMin

(

EvtCyclic3::Pair

i

)

const

Definition at line 95 of file EvtDalitzPlot.cc.

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

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

qHelAbsMax()

double EvtDalitzPlot::qHelAbsMax

(

EvtCyclic3::Pair

i

)

const

Definition at line 121 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 115 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 165 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 cosTh(), getArea(), and EvtDalitzPoint::qMax().

qMin()

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

Definition at line 133 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 cosTh(), getArea(), and EvtDalitzPoint::qMin().

qResAbsMax()

double EvtDalitzPlot::qResAbsMax

(

EvtCyclic3::Pair

i

)

const

Definition at line 111 of file EvtDalitzPlot.cc.

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

qResAbsMin()

double EvtDalitzPlot::qResAbsMin

(

EvtCyclic3::Pair

i

)

const

Definition at line 107 of file EvtDalitzPlot.cc.

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

qSumMax()

double EvtDalitzPlot::qSumMax ( ) const

inline

Definition at line 53 of file EvtDalitzPlot.hh.

{ return sum() + _rdel; }

Referenced by print().

qSumMin()

double EvtDalitzPlot::qSumMin ( ) const

inline

Definition at line 52 of file EvtDalitzPlot.hh.

{ return sum() + _ldel; }

Referenced by print().

sanityCheck()

void EvtDalitzPlot::sanityCheck

(

)

const

Definition at line 73 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 93 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 279 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 275 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 106 of file EvtDalitzPlot.hh.

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

_ldel

double EvtDalitzPlot::_ldel

protected

Definition at line 110 of file EvtDalitzPlot.hh.

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

_mA

double EvtDalitzPlot::_mA

protected

Definition at line 103 of file EvtDalitzPlot.hh.

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

_mB

double EvtDalitzPlot::_mB

protected

Definition at line 104 of file EvtDalitzPlot.hh.

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

_mC

double EvtDalitzPlot::_mC

protected

Definition at line 105 of file EvtDalitzPlot.hh.

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

_rdel

double EvtDalitzPlot::_rdel

protected

Definition at line 111 of file EvtDalitzPlot.hh.

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

---

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

• EvtDalitzPlot.hh
• EvtDalitzPlot.cc