EvtItgAbsFunction.cc

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Copyright Information: See EvtGen/COPYRIGHT
//
// Environment:
//      This software is part of the EvtGen package developed jointly
//      for the BaBar and CLEO collaborations.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Module: EvtItgAbsFunction.hh
//
// Description:
//      Abstraction of a generic function for use in integration methods elsewhere
//      in this package. (Stolen and modified from the BaBar IntegrationUtils package
//      - author: Phil Strother).
//
// Modification history:
//
//    Jane Tinslay                March 21, 2001       Module adapted for use in
//                                                     EvtGen
//
//------------------------------------------------------------------------
#include "../EvtGenBase/EvtPatches.hh"
 
#include "EvtItgAbsFunction.hh"
 
//-------------
// C Headers --
//-------------
extern "C" {}
#include "../EvtGenBase/EvtReport.hh"
#include "assert.h"
using std::endl;
 
EvtItgAbsFunction::EvtItgAbsFunction( double lowerRange, double upperRange )
    : _upperRange( upperRange ), _lowerRange( lowerRange ) {}
 
EvtItgAbsFunction::~EvtItgAbsFunction() {}
 
double EvtItgAbsFunction::value( double x ) const {
  if ( x >= _lowerRange && x <= _upperRange ) return myFunction( x );
  report( ERROR, "EvtGen" ) << "Error in EvtItgAbsFunction::value.  Given co-ordinate " << x
                            << " is outside of allowed range [" << _lowerRange << ", "
                            << _upperRange << "].  Returning 0.0" << endl;
  return 0.0; // Never get here
}
 
double EvtItgAbsFunction::operator()( double x ) const { return myFunction( x ); }