EvtDalitzFlatPdf Class Reference

#include <EvtDalitzFlatPdf.hh>

Inheritance diagram for EvtDalitzFlatPdf:

Public Member Functions
	EvtDalitzFlatPdf (const EvtDalitzPlot &dp)
	EvtDalitzFlatPdf (const EvtDalitzFlatPdf &other)
virtual	~EvtDalitzFlatPdf ()
virtual EvtPdf< EvtDalitzPoint > *	clone () const
virtual EvtValError	compute_integral (int N) const
virtual EvtDalitzPoint	randomPoint ()
Public Member Functions inherited from EvtPdf< EvtDalitzPoint >
	EvtPdf ()
	EvtPdf (const EvtPdf &other)
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
EvtValError	getItg (int N) 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
Protected Member Functions inherited from EvtPdf< EvtDalitzPoint >

Protected Attributes
EvtDalitzPlot	_dp
Protected Attributes inherited from EvtPdf< EvtDalitzPoint >
EvtValError	_itg

Detailed Description

Definition at line 22 of file EvtDalitzFlatPdf.hh.

Constructor & Destructor Documentation

EvtDalitzFlatPdf() [1/2]

EvtDalitzFlatPdf::EvtDalitzFlatPdf

(

const EvtDalitzPlot &

dp

)

Definition at line 14 of file EvtDalitzFlatPdf.cc.

  : EvtPdf<EvtDalitzPoint>(), _dp(dp)
{}

Referenced by clone().

EvtDalitzFlatPdf() [2/2]

EvtDalitzFlatPdf::EvtDalitzFlatPdf

(

const EvtDalitzFlatPdf &

other

)

Definition at line 18 of file EvtDalitzFlatPdf.cc.

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

~EvtDalitzFlatPdf()

EvtDalitzFlatPdf::~EvtDalitzFlatPdf ( )

virtual

Definition at line 22 of file EvtDalitzFlatPdf.cc.

{}

Member Function Documentation

clone()

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

virtual

Implements EvtPdf< EvtDalitzPoint >.

Definition at line 25 of file EvtDalitzFlatPdf.cc.

{
  return new EvtDalitzFlatPdf(*this);
}

compute_integral()

EvtValError EvtDalitzFlatPdf::compute_integral ( int N ) const

virtual

Reimplemented from EvtPdf< EvtDalitzPoint >.

Definition at line 35 of file EvtDalitzFlatPdf.cc.

{
  return EvtValError(_dp.getArea(N),0.);
}

pdf()

double EvtDalitzFlatPdf::pdf ( const EvtDalitzPoint & ) const

protectedvirtual

Implements EvtPdf< EvtDalitzPoint >.

Definition at line 30 of file EvtDalitzFlatPdf.cc.

{
  return 1.;
}

randomPoint()

EvtDalitzPoint EvtDalitzFlatPdf::randomPoint ( )

virtual

Reimplemented from EvtPdf< EvtDalitzPoint >.

Definition at line 40 of file EvtDalitzFlatPdf.cc.

{
  // To obtain a uniform distribution generate 
  // in terms of q's. Generate in a box that circumscribes the 
  // Dalitz plot. Accept points inside. If there are two 
  // many unsuccessful attempts it's a hint that the Dalitz plot
  // area is tiny compared to the box. It's a pathological
  // case. Abort.
  
  EvtCyclic3::Pair pair1 = EvtCyclic3::BC;
  EvtCyclic3::Pair pair2 = EvtCyclic3::CA;
  
  int n = 0;
  int maxTries = 1000;
  while(n++ < maxTries) {
    
    double q1 = EvtRandom::Flat(_dp.qAbsMin(pair1),_dp.qAbsMax(pair2));
    double q2 = EvtRandom::Flat(_dp.qAbsMin(pair2),_dp.qAbsMax(pair2)); 
    
    EvtDalitzCoord point(pair1,q1,pair2,q2);
    EvtDalitzPoint x(_dp,point);
    
    if(x.isValid()) return x;
  }
  
  printf("No point generated for dalitz plot after %d tries\n",maxTries);
  assert(0);
}

Member Data Documentation

_dp

EvtDalitzPlot EvtDalitzFlatPdf::_dp

protected

Definition at line 37 of file EvtDalitzFlatPdf.hh.

Referenced by compute_integral(), and randomPoint().

---

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

• EvtDalitzFlatPdf.hh
• EvtDalitzFlatPdf.cc