EvtDalitzResPdf Class Reference

#include <EvtDalitzResPdf.hh>

Inheritance diagram for EvtDalitzResPdf:

Public Member Functions
	EvtDalitzResPdf (const EvtDalitzPlot &dp, double m0, double g0, EvtCyclic3::Pair pairRes)
	EvtDalitzResPdf (const EvtDalitzResPdf &other)
virtual	~EvtDalitzResPdf ()
EvtPdf< EvtDalitzPoint > *	clone () const
virtual EvtValError	compute_integral (int N) const
virtual EvtDalitzPoint	randomPoint ()
double	pdfMaxValue () const
Public Member Functions inherited from EvtPdf< EvtDalitzPoint >
	EvtPdf ()
virtual	~EvtPdf ()
double	evaluate (const EvtDalitzPoint &p) const
EvtPdfMax< EvtDalitzPoint >	findMax (const EvtPdf< EvtDalitzPoint > &pc, int N)
EvtValError	findGenEff (const EvtPdf< EvtDalitzPoint > &pc, int N, int nFindMax)
void	setItg (EvtValError itg)
EvtValError	getItg () const
virtual EvtValError	compute_integral () const
EvtValError	compute_mc_integral (const EvtPdf< EvtDalitzPoint > &pc, int N)
EvtPredGen< EvtPdfGen< EvtDalitzPoint >, EvtPdfPred< EvtDalitzPoint > >	accRejGen (const EvtPdf< EvtDalitzPoint > &pc, int nMax, double factor=1.)

Protected Member Functions
virtual double	pdf (const EvtDalitzPoint &) const

Additional Inherited Members
Protected Attributes inherited from EvtPdf< EvtDalitzPoint >
EvtValError	_itg

Detailed Description

Definition at line 31 of file EvtDalitzResPdf.hh.

Constructor & Destructor Documentation

EvtDalitzResPdf() [1/2]

EvtDalitzResPdf::EvtDalitzResPdf	(	const EvtDalitzPlot &	dp,
		double	m0,
		double	g0,
		EvtCyclic3::Pair	pairRes )

Definition at line 21 of file EvtDalitzResPdf.cc.

    : EvtPdf<EvtDalitzPoint>(), _dp( dp ), _m0( _m0 ), _g0( _g0 ), _pair( pair ) {}

Referenced by clone(), and EvtDalitzResPdf().

EvtDalitzResPdf() [2/2]

EvtDalitzResPdf::EvtDalitzResPdf

(

const EvtDalitzResPdf &

other

)

Definition at line 25 of file EvtDalitzResPdf.cc.

    : EvtPdf<EvtDalitzPoint>( other )
    , _dp( other._dp )
    , _m0( other._m0 )
    , _g0( other._g0 )
    , _pair( other._pair ) {}

~EvtDalitzResPdf()

EvtDalitzResPdf::~EvtDalitzResPdf ( )

virtual

Definition at line 32 of file EvtDalitzResPdf.cc.

{}

Member Function Documentation

clone()

EvtPdf< EvtDalitzPoint > * EvtDalitzResPdf::clone ( ) const

inlinevirtual

Implements EvtPdf< EvtDalitzPoint >.

Definition at line 38 of file EvtDalitzResPdf.hh.

{ return new EvtDalitzResPdf( *this ); }

compute_integral()

EvtValError EvtDalitzResPdf::compute_integral ( int N ) const

virtual

Definition at line 34 of file EvtDalitzResPdf.cc.

                                                           {
  assert( N != 0 );
 
  EvtCyclic3::Pair i = _pair;
  EvtCyclic3::Pair j = EvtCyclic3::next( i );
 
  // Trapezoidal integral
 
  double dh  = ( _dp.qAbsMax( j ) - _dp.qAbsMin( j ) ) / ( (double)N );
  double sum = 0;
 
  int ii;
  for ( ii = 1; ii < N; ii++ )
  {
 
    double x   = _dp.qAbsMin( j ) + ii * dh;
    double min = ( _dp.qMin( i, j, x ) - _m0 * _m0 ) / _m0 / _g0;
    double max = ( _dp.qMax( i, j, x ) - _m0 * _m0 ) / _m0 / _g0;
    double itg = 1 / EvtConst::pi * ( atan( max ) - atan( min ) );
    sum += itg;
  }
  EvtValError ret( sum * dh, 0. );
 
  return ret;
}

pdf()

double EvtDalitzResPdf::pdf ( const EvtDalitzPoint & x ) const

protectedvirtual

Implements EvtPdf< EvtDalitzPoint >.

Definition at line 88 of file EvtDalitzResPdf.cc.

                                                           {
  EvtCyclic3::Pair i  = _pair;
  double           dq = x.q( i ) - _m0 * _m0;
  return 1 / EvtConst::pi * _g0 * _m0 / ( dq * dq + _g0 * _g0 * _m0 * _m0 );
}

pdfMaxValue()

double EvtDalitzResPdf::pdfMaxValue

(

)

const

Definition at line 94 of file EvtDalitzResPdf.cc.

{ return 1 / ( EvtConst::pi * _g0 * _m0 ); }

randomPoint()

EvtDalitzPoint EvtDalitzResPdf::randomPoint ( )

virtual

Reimplemented from EvtPdf< EvtDalitzPoint >.

Definition at line 60 of file EvtDalitzResPdf.cc.

                                            {
  // Random point generation must be done in a box encompassing the
  // Dalitz plot
 
  EvtCyclic3::Pair i = _pair;
  EvtCyclic3::Pair j = EvtCyclic3::next( i );
  double min         = 1 / EvtConst::pi * atan( ( _dp.qAbsMin( i ) - _m0 * _m0 ) / _m0 / _g0 );
  double max         = 1 / EvtConst::pi * atan( ( _dp.qAbsMax( i ) - _m0 * _m0 ) / _m0 / _g0 );
 
  int n = 0;
  while ( n++ < 1000 )
  {
 
    double         qj = EvtRandom::Flat( _dp.qAbsMin( j ), _dp.qAbsMax( j ) );
    double         r  = EvtRandom::Flat( min, max );
    double         qi = tan( EvtConst::pi * r ) * _g0 * _m0 + _m0 * _m0;
    EvtDalitzCoord x( i, qi, j, qj );
    EvtDalitzPoint ret( _dp, x );
    if ( ret.isValid() ) return ret;
  }
 
  // All generated points turned out to be outside of the Dalitz plot
  // (in the outer box)
 
  printf( "No point generated for dalitz plot after 1000 tries\n" );
  abort();
}

---

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

• EvtDalitzResPdf.hh
• EvtDalitzResPdf.cc