EvtPdf< T > Class Template Reference

#include <EvtPdf.hh>

Inheritance diagram for EvtPdf< T >:

Public Member Functions
	EvtPdf ()
	EvtPdf (const EvtPdf &other)
virtual	~EvtPdf ()
virtual EvtPdf< T > *	clone () const =0
double	evaluate (const T &p) const
EvtPdfMax< T >	findMax (const EvtPdf< T > &pc, int N)
EvtValError	findGenEff (const EvtPdf< T > &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< T > &pc, int N)
EvtPredGen< EvtPdfGen< T >, EvtPdfPred< T > >	accRejGen (const EvtPdf< T > &pc, int nMax, double factor=1.)
virtual T	randomPoint ()

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

Protected Attributes
EvtValError	_itg

Detailed Description

template<class T>
class EvtPdf< T >

Definition at line 57 of file EvtPdf.hh.

Constructor & Destructor Documentation

EvtPdf() [1/2]

template<class T>

EvtPdf< T >::EvtPdf ( )

inline

Definition at line 60 of file EvtPdf.hh.

{}

Referenced by accRejGen(), EvtPdfSum< T >::addOwnedTerm(), EvtPdfSum< T >::addTerm(), clone(), EvtPdfDiv< T >::clone(), EvtPdfSum< T >::clone(), compute_mc_integral(), EvtAmpAmpPdf< T >::EvtAmpAmpPdf(), EvtAmpAmpPdf< T >::EvtAmpAmpPdf(), EvtAmpPdf< T >::EvtAmpPdf(), EvtAmpPdf< T >::EvtAmpPdf(), EvtPdf(), EvtPdfDiv< T >::EvtPdfDiv(), EvtPdfDiv< T >::EvtPdfDiv(), EvtPdfSum< T >::EvtPdfSum(), findGenEff(), findMax(), and EvtPdfSum< T >::getPdf().

EvtPdf() [2/2]

template<class T>

EvtPdf< T >::EvtPdf ( const EvtPdf< T > & other )

inline

Definition at line 61 of file EvtPdf.hh.

: _itg(other._itg) {}

~EvtPdf()

template<class T>

virtual EvtPdf< T >::~EvtPdf ( )

inlinevirtual

Definition at line 62 of file EvtPdf.hh.

{}

Member Function Documentation

accRejGen()

template<class T>

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

Definition at line 301 of file EvtPdf.hh.

{
  EvtPdfGen<T> gen(pc);
  EvtPdfDiv<T> pdfdiv(*this,pc);
  EvtPdfPred<T> pred(pdfdiv);
  pred.compute_max(iter(gen,nMax),iter(gen),factor);
  return EvtPredGen<EvtPdfGen<T>,EvtPdfPred<T> >(gen,pred);
}

Referenced by findGenEff().

clone()

template<class T>

virtual EvtPdf< T > * EvtPdf< T >::clone ( ) const

pure virtual

Implemented in EvtAmpAmpPdf< T >, EvtAmpPdf< T >, EvtBreitWignerPdf, EvtDalitzFlatPdf, EvtDalitzResPdf, EvtIntervalFlatPdf, EvtPdfDiv< T >, EvtPdfSum< T >, and EvtPdfSum< EvtDalitzPoint >.

compute_integral() [1/2]

template<class T>

virtual EvtValError EvtPdf< T >::compute_integral ( ) const

inlinevirtual

Reimplemented in EvtIntegPdf1D, EvtIntervalFlatPdf, EvtPdfSum< T >, and EvtPdfSum< EvtDalitzPoint >.

Definition at line 92 of file EvtPdf.hh.

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

Referenced by compute_integral(), getItg(), and getItg().

compute_integral() [2/2]

template<class T>

virtual EvtValError EvtPdf< T >::compute_integral ( int N ) const

inlinevirtual

Reimplemented in EvtDalitzFlatPdf, EvtDalitzResPdf, EvtPdfSum< T >, and EvtPdfSum< EvtDalitzPoint >.

Definition at line 95 of file EvtPdf.hh.

{ return compute_integral(); }

Referenced by compute_integral().

compute_mc_integral()

template<class T>

EvtValError EvtPdf< T >::compute_mc_integral	(	const EvtPdf< T > &	pc,
		int	N )

Definition at line 253 of file EvtPdf.hh.

{
  assert(N > 0);
 
  EvtValError otherItg = pc.getItg();
  EvtPdfDiv<T> pdfdiv(*this,pc);
  EvtPdfUnary<T> unary(pdfdiv);  
  
  EvtPdfGen<T> gen(pc);    
  EvtStreamInputIterator<T> begin = iter(gen,N);
  EvtStreamInputIterator<T> end;
 
  double sum = 0.;
  double sum2 = 0.;
  while(!(begin == end)) {
    
    double value = pdfdiv.evaluate(*begin++);
    sum += value;
    sum2 += value*value;
  }
  
  EvtValError x;
  if(N > 0) {
    double av = sum/((double) N);
    if(N > 1) {
      double dev2 = (sum2 - av*av*N)/((double) (N - 1));
      // Due to numerical precision dev2 may sometimes be negative
      if(dev2 < 0.) dev2 = 0.;
      double error = sqrt(dev2/((double) N));
      x = EvtValError(av,error);
    }
    else x = EvtValError(av);
  }
  _itg = x * pc.getItg();
  return _itg;
}

evaluate()

template<class T>

double EvtPdf< T >::evaluate ( const T & p ) const

inline

Definition at line 65 of file EvtPdf.hh.

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

Referenced by compute_mc_integral(), EvtBtoKD3P::decay(), EvtIntervalDecayAmp< T >::decay(), EvtRelBreitWignerBarrierFact::getRandMass(), and EvtPdfSum< T >::pdf().

findGenEff()

template<class T>

EvtValError EvtPdf< T >::findGenEff	(	const EvtPdf< T > &	pc,
		int	N,
		int	nFindMax )

Definition at line 241 of file EvtPdf.hh.

{
  assert(N > 0 || nFindMax > 0);
  EvtPredGen<EvtPdfGen<T>,EvtPdfPred<T> > gen = accRejGen(pc,nFindMax);
  int i;
  for(i=0;i<N;i++) gen();
  double eff = double(gen.getPassed())/double(gen.getTried());
  double err = sqrt(double(gen.getPassed()))/double(gen.getTried());
  return EvtValError(eff,err);
}

findMax()

template<class T>

EvtPdfMax< T > EvtPdf< T >::findMax	(	const EvtPdf< T > &	pc,
		int	N )

Definition at line 230 of file EvtPdf.hh.

{
  EvtPdfPred<T> pred(*this);
  EvtPdfGen<T> gen(pc);
  pred.compute_max(iter(gen,N),iter(gen));
  EvtPdfMax<T> p = pred.getMax();
  return p;
}

Referenced by EvtIntervalDecayAmp< T >::initProbMax().

getItg() [1/2]

template<class T>

EvtValError EvtPdf< T >::getItg ( ) const

inline

Definition at line 83 of file EvtPdf.hh.

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

Referenced by EvtPdfSum< T >::compute_integral(), EvtPdfSum< T >::compute_integral(), and compute_mc_integral().

getItg() [2/2]

template<class T>

EvtValError EvtPdf< T >::getItg ( int N ) const

inline

Definition at line 87 of file EvtPdf.hh.

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

pdf()

template<class T>

virtual double EvtPdf< T >::pdf ( const T & ) const

protectedpure virtual

Implemented in EvtAmpAmpPdf< T >, EvtAmpPdf< T >, EvtBreitWignerPdf, EvtDalitzFlatPdf, EvtDalitzResPdf, EvtIntervalFlatPdf, EvtPdfDiv< T >, and EvtPdfSum< T >.

Referenced by evaluate().

randomPoint()

template<class T>

T EvtPdf< T >::randomPoint ( )

virtual

Reimplemented in EvtDalitzFlatPdf, EvtDalitzResPdf, EvtIntegPdf1D, EvtIntervalFlatPdf, EvtPdfSum< T >, and EvtPdfSum< EvtDalitzPoint >.

Definition at line 291 of file EvtPdf.hh.

{
  printf("Function defined for analytic PDFs only\n");
  assert(0);
  T temp;
  return temp;
}

setItg()

template<class T>

void EvtPdf< T >::setItg ( EvtValError itg )

inline

Definition at line 81 of file EvtPdf.hh.

{_itg = itg; }

Member Data Documentation

_itg

template<class T>

EvtValError EvtPdf< T >::_itg

mutableprotected

Definition at line 111 of file EvtPdf.hh.

Referenced by compute_mc_integral(), EvtPdf(), getItg(), getItg(), and setItg().

---

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

• EvtPdf.hh