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 59 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 60 of file EvtPdf.hh.

: _itg( other._itg ) {}

~EvtPdf()

template<class T>

virtual EvtPdf< T >::~EvtPdf ( )

inlinevirtual

Definition at line 61 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 281 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 91 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 EvtPdfSum< T >.

Definition at line 98 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 234 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 64 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 224 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 215 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 82 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 86 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 273 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 80 of file EvtPdf.hh.

{ _itg = itg; }

Member Data Documentation

_itg

template<class T>

EvtValError EvtPdf< T >::_itg

mutableprotected

Definition at line 114 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