EvtDalitzFlatPdf Class Reference

#include <EvtDalitzFlatPdf.hh>

Inheritance diagram for EvtDalitzFlatPdf:

List of all members.

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 ()
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 Attributes
EvtDalitzPlot	_dp
EvtValError	_itg

---

Detailed Description

Definition at line 22 of file EvtDalitzFlatPdf.hh.

---

Constructor & Destructor Documentation

EvtDalitzFlatPdf::EvtDalitzFlatPdf

(

const EvtDalitzPlot &

dp

)

Definition at line 14 of file EvtDalitzFlatPdf.cc.

Referenced by clone().

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

EvtDalitzFlatPdf::EvtDalitzFlatPdf

(

const EvtDalitzFlatPdf &

other

)

Definition at line 18 of file EvtDalitzFlatPdf.cc.

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

EvtDalitzFlatPdf::~EvtDalitzFlatPdf

(

)

[virtual]

Definition at line 22 of file EvtDalitzFlatPdf.cc.

{}

---

Member Function Documentation

EvtPredGen<EvtPdfGen<EvtDalitzPoint >,EvtPdfPred<EvtDalitzPoint > > EvtPdf< EvtDalitzPoint >::accRejGen	(	const EvtPdf< EvtDalitzPoint > &	pc,
		int	nMax,
		double	factor = 1.
	)			[inherited]

EvtPdf< EvtDalitzPoint > * EvtDalitzFlatPdf::clone

(

)

const [virtual]

Implements EvtPdf< EvtDalitzPoint >.

Definition at line 25 of file EvtDalitzFlatPdf.cc.

References EvtDalitzFlatPdf().

{
  return new EvtDalitzFlatPdf(*this);
}

virtual EvtValError EvtPdf< EvtDalitzPoint >::compute_integral

(

)

const [inline, virtual, inherited]

Reimplemented in EvtPdfSum< EvtDalitzPoint >.

Definition at line 92 of file EvtPdf.hh.

  { printf("Analytic integration of PDF is not defined\n"); assert(0); return compute_integral();}

EvtValError EvtDalitzFlatPdf::compute_integral

(

int

N

)

const [virtual]

Reimplemented from EvtPdf< EvtDalitzPoint >.

Definition at line 35 of file EvtDalitzFlatPdf.cc.

References _dp, and EvtDalitzPlot::getArea().

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

EvtValError EvtPdf< EvtDalitzPoint >::compute_mc_integral	(	const EvtPdf< EvtDalitzPoint > &	pc,
		int	N
	)			[inherited]

double EvtPdf< EvtDalitzPoint >::evaluate

(

const EvtDalitzPoint &

p

)

const [inline, inherited]

Definition at line 65 of file EvtPdf.hh.

                                    { 
    if(p.isValid()) return pdf(p); 
    else return 0.;
  }

EvtValError EvtPdf< EvtDalitzPoint >::findGenEff	(	const EvtPdf< EvtDalitzPoint > &	pc,
		int	N,
		int	nFindMax
	)			[inherited]

EvtPdfMax<EvtDalitzPoint > EvtPdf< EvtDalitzPoint >::findMax	(	const EvtPdf< EvtDalitzPoint > &	pc,
		int	N
	)			[inherited]

EvtValError EvtPdf< EvtDalitzPoint >::getItg

(

int

N

)

const [inline, inherited]

Definition at line 87 of file EvtPdf.hh.

                                  {
    if(!_itg.valueKnown()) _itg = compute_integral(N);
    return _itg;
  }

EvtValError EvtPdf< EvtDalitzPoint >::getItg

(

)

const [inline, inherited]

Definition at line 83 of file EvtPdf.hh.

                             {
    if(!_itg.valueKnown()) _itg = compute_integral();
    return _itg;
  }

double EvtDalitzFlatPdf::pdf

(

const EvtDalitzPoint &

)

const [protected, virtual]

Implements EvtPdf< EvtDalitzPoint >.

Definition at line 30 of file EvtDalitzFlatPdf.cc.

{
  return 1.;
}

EvtDalitzPoint EvtDalitzFlatPdf::randomPoint

(

)

[virtual]

Reimplemented from EvtPdf< EvtDalitzPoint >.

Definition at line 40 of file EvtDalitzFlatPdf.cc.

References _dp, EvtCyclic3::BC, EvtCyclic3::CA, EvtRandom::Flat(), EvtDalitzPlot::qAbsMax(), EvtDalitzPlot::qAbsMin(), and x.

{
  // 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);
}

void EvtPdf< EvtDalitzPoint >::setItg

(

EvtValError

itg

)

[inline, inherited]

Definition at line 81 of file EvtPdf.hh.

{_itg = itg; }

---

Member Data Documentation

EvtDalitzPlot EvtDalitzFlatPdf::_dp [protected]

Definition at line 37 of file EvtDalitzFlatPdf.hh.

Referenced by compute_integral(), and randomPoint().

EvtValError EvtPdf< EvtDalitzPoint >::_itg [mutable, protected, inherited]

Reimplemented in EvtPdfSum< EvtDalitzPoint >.

Definition at line 111 of file EvtPdf.hh.

---