EvtItgAbsIntegrator Class Reference

#include <EvtItgAbsIntegrator.hh>

Inheritance diagram for EvtItgAbsIntegrator:

List of all members.

Public Member Functions
	EvtItgAbsIntegrator (const EvtItgAbsFunction &)
virtual	~EvtItgAbsIntegrator ()
double	evaluate (double lower, double upper) const
double	normalisation () const
Protected Member Functions
double	trapezoid (double lower, double higher, int n, double &result) const
virtual double	evaluateIt (double lower, double higher) const =0
double	myFunction (double x) const
Private Member Functions
void	boundsCheck (double &, double &) const
	EvtItgAbsIntegrator ()
	EvtItgAbsIntegrator (const EvtItgAbsIntegrator &)
EvtItgAbsIntegrator &	operator= (const EvtItgAbsIntegrator &)
Private Attributes
const EvtItgAbsFunction &	_myFunction

---

Detailed Description

Definition at line 30 of file EvtItgAbsIntegrator.hh.

---

Constructor & Destructor Documentation

EvtItgAbsIntegrator::EvtItgAbsIntegrator

(

const EvtItgAbsFunction &

)

Definition at line 39 of file EvtItgAbsIntegrator.cc.

                                                                            :
  _myFunction(theFunction)
{}

EvtItgAbsIntegrator::~EvtItgAbsIntegrator

(

)

[virtual]

Definition at line 43 of file EvtItgAbsIntegrator.cc.

{}

EvtItgAbsIntegrator::EvtItgAbsIntegrator

(

)

[private]

EvtItgAbsIntegrator::EvtItgAbsIntegrator

(

const EvtItgAbsIntegrator &

)

[private]

---

Member Function Documentation

void EvtItgAbsIntegrator::boundsCheck	(	double &	,
		double &
	)			const [private]

Definition at line 89 of file EvtItgAbsIntegrator.cc.

References _myFunction, EvtItgAbsFunction::lowerRange(), report(), EvtItgAbsFunction::upperRange(), and Bes_Common::WARNING.

Referenced by evaluate().

                                                                  {

  if (lower < _myFunction.lowerRange() ) {
    report(WARNING,"EvtGen") << "Warning in EvtItgAbsIntegrator::evaluate.  Lower bound " << lower << " of integral " 
                    << " is less than lower bound " << _myFunction.lowerRange() 
                    << " of function.  No contribution from this range will be counted." << endl;
    lower = _myFunction.lowerRange();
  }

  if (upper > _myFunction.upperRange() ) {
    report(WARNING,"EvtGen") << "Warning in EvtItgAbsIntegrator::evaluate.  Upper bound " << upper << " of integral "
                    << " is greater than upper bound " << _myFunction.upperRange() 
                    << " of function.  No contribution from this range will be counted." << endl;  
    upper = _myFunction.upperRange();
  }

}

double EvtItgAbsIntegrator::evaluate	(	double	lower,
		double	upper
	)			const

Definition at line 52 of file EvtItgAbsIntegrator.cc.

References boundsCheck(), and evaluateIt().

Referenced by EvtBtoXsgammaKagan::computeHadronicMass(), EvtVubNLO::Gamma(), EvtBtoXsgammaRootFinder::GetGaussIntegFcnRoot(), and EvtVubNLO::tripleDiff().

                                                             {

  double newLower(lower), newUpper(upper);

  boundsCheck(newLower, newUpper);

  return evaluateIt(newLower, newUpper);
}

virtual double EvtItgAbsIntegrator::evaluateIt	(	double	lower,
		double	higher
	)			const [protected, pure virtual]

Implemented in EvtItgSimpsonIntegrator.

Referenced by evaluate(), and normalisation().

double EvtItgAbsIntegrator::myFunction

(

double

x

)

const [inline, protected]

Definition at line 49 of file EvtItgAbsIntegrator.hh.

References _myFunction.

{return _myFunction(x);}

double EvtItgAbsIntegrator::normalisation

(

)

const

Definition at line 47 of file EvtItgAbsIntegrator.cc.

References _myFunction, evaluateIt(), EvtItgAbsFunction::lowerRange(), and EvtItgAbsFunction::upperRange().

Referenced by EvtBtoXsgammaKagan::computeHadronicMass().

                                         {
  return evaluateIt(_myFunction.lowerRange(), _myFunction.upperRange());
}

EvtItgAbsIntegrator& EvtItgAbsIntegrator::operator=

(

const EvtItgAbsIntegrator &

)

[private]

double EvtItgAbsIntegrator::trapezoid	(	double	lower,
		double	higher,
		int	n,
		double &	result
	)			const [protected]

Definition at line 62 of file EvtItgAbsIntegrator.cc.

References _myFunction, ganga-rec::j, and x.

Referenced by EvtItgSimpsonIntegrator::evaluateIt().

                                                                                       {

  if (n==1) return 0.5*(higher-lower)*(_myFunction(lower) + _myFunction(higher));
  
  int it, j;
  
  for (it=1, j=1;j<n-1;j++) it <<=1;
  
  double itDouble(it);
  
  double sum(0.0);

  double deltaX((higher - lower)/itDouble);
  
  double x(lower + 0.5* deltaX);
    
  for (j=1;j<=it;j++){
    sum+=_myFunction(x);
    x+=deltaX;
  }
  
  result = 0.5*(result+(higher - lower)*sum/itDouble);

  return result;
}

---

Member Data Documentation

const EvtItgAbsFunction& EvtItgAbsIntegrator::_myFunction [private]

Definition at line 53 of file EvtItgAbsIntegrator.hh.

Referenced by boundsCheck(), myFunction(), normalisation(), and trapezoid().

---