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 ()
	EvtPdf (const EvtPdf &other)
virtual	~EvtPdf ()
virtual EvtPdf< EvtDalitzPoint > *	clone () const=0
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
EvtValError	getItg (int N) const
virtual EvtValError	compute_integral () const
virtual EvtValError	compute_integral (int N) 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.)
virtual EvtDalitzPoint	randomPoint ()

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

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 20 of file EvtDalitzResPdf.cc.

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

EvtDalitzResPdf() [2/2]

EvtDalitzResPdf::EvtDalitzResPdf

(

const EvtDalitzResPdf &

other

)

Definition at line 27 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 33 of file EvtDalitzResPdf.cc.

{}

Member Function Documentation

clone()

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

inlinevirtual

Implements EvtPdf< EvtDalitzPoint >.

Definition at line 40 of file EvtDalitzResPdf.hh.

{ return new EvtDalitzResPdf(*this); }

compute_integral()

EvtValError EvtDalitzResPdf::compute_integral ( int  N ) const

virtual

Reimplemented from EvtPdf< EvtDalitzPoint >.

Definition at line 36 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 93 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 101 of file EvtDalitzResPdf.cc.

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

randomPoint()

EvtDalitzPoint EvtDalitzResPdf::randomPoint ( )

virtual

Reimplemented from EvtPdf< EvtDalitzPoint >.

Definition at line 63 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");
  assert(0);
}

---

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

• EvtDalitzResPdf.hh
• EvtDalitzResPdf.cc