FastAmplitudeIntegrator Class Reference

#include <FastAmplitudeIntegrator.h>

Inheritance diagram for FastAmplitudeIntegrator:

Public Member Functions
	FastAmplitudeIntegrator ()
	FastAmplitudeIntegrator (const DalitzEventPattern &pattern, IFastAmplitudeIntegrable *amps=0, MINT::IEventGenerator< IDalitzEvent > *eventGenerator=0, TRandom *rnd=0, double precision=1.e-2)
	FastAmplitudeIntegrator (const DalitzEventPattern &pattern, IFastAmplitudeIntegrable *amps, const std::string &fname)
virtual	~FastAmplitudeIntegrator ()
bool	initialise (const DalitzEventPattern &pattern, IFastAmplitudeIntegrable *amps=0, MINT::IEventGenerator< IDalitzEvent > *eventGenerator=0, TRandom *rnd=0, double precision=1.e-2)
bool	initialise (const std::string &commaSeparatedList, const DalitzEventPattern &pattern, IFastAmplitudeIntegrable *amps=0, MINT::IEventGenerator< IDalitzEvent > *eventGenerator=0, TRandom *rnd=0, double precision=1.e-2)
bool	initialiseFromFile (const DalitzEventPattern &pattern, IFastAmplitudeIntegrable *amps, const std::string &commaSeparatedList)
virtual bool	add (const FastAmplitudeIntegrator &other)
bool	setValues (const DalitzEventPattern &pattern, IFastAmplitudeIntegrable *amps=0, MINT::IEventGenerator< IDalitzEvent > *eventGenerator=0, TRandom *rnd=0, double precision=1.e-2)
void	setPrecision (double prec)
bool	initialised () const
double	getVal ()
void	Gradient (MINT::MinuitParameterSet *mps, std::vector< double > &grad)
void	GradientForLasso (MINT::MinuitParameterSet *mps, std::vector< double > &grad)
double	RealVal ()
double	sumOfFitFractions ()
double	absSumOfSqrtInterferenceFractions ()
double	absSumOfInterferenceFractions ()
double	sumOfSqrtFitFractions ()
int	numberOfFitFractionsLargerThanThreshold (double threshold)
double	variance () const
DalitzHistoSet	histoSet () const
void	saveEachAmpsHistograms (const std::string &prefix) const
std::vector< DalitzHistoSet >	GetEachAmpsHistograms ()
DalitzHistoSet	interferenceHistoSet () const
void	saveInterferenceHistograms (const std::string &prefix) const
std::vector< DalitzHistoSet >	GetInterferenceHistograms ()
virtual void	doFinalStats (MINT::Minimiser *mini=0)
double	getFractionChi2 () const
bool	ensureFreshEvents ()
bool	save (const std::string &fname) const
Public Member Functions inherited from MINT::IReturnReal
virtual	~IReturnReal ()

Protected Member Functions
double	weight (IDalitzEvent *evtPtr)
int	updateEventSet (long int Nevents)
double	evaluateSum ()
int	determineNumEvents ()
int	generateEnoughEvents ()
int	addEvents (long int Nevents)
TRandom *	makeNewRnd (int seed=-9999)
Protected Member Functions inherited from MINT::IReturnReal
	IReturnReal ()

Protected Attributes
DalitzEventPattern	_pat
IFastAmplitudeIntegrable *	_amps
MINT::counted_ptr< IntegCalculator >	_integCalc
MINT::counted_ptr< IntegCalculator >	_integCalc_copyForDebug
TRandom *	_rnd
MINT::counted_ptr< TRandom >	_localRnd
double	_precision
MINT::IEventGenerator< IDalitzEvent > *	_generator

Private Attributes
double	_mean
int	_numEvents
int	_Ncalls
bool	_initialised
bool	_db

Static Private Attributes
static const long int	_minEvents =100

Detailed Description

Definition at line 32 of file FastAmplitudeIntegrator.h.

Constructor & Destructor Documentation

FastAmplitudeIntegrator() [1/3]

FastAmplitudeIntegrator::FastAmplitudeIntegrator

(

)

Definition at line 45 of file FastAmplitudeIntegrator.cpp.

  : _initialised(0), _db(false)
{
}

FastAmplitudeIntegrator() [2/3]

FastAmplitudeIntegrator::FastAmplitudeIntegrator	(	const DalitzEventPattern &	pattern,
		IFastAmplitudeIntegrable *	amps = 0,
		MINT::IEventGenerator< IDalitzEvent > *	eventGenerator = 0,
		TRandom *	rnd = 0,
		double	precision = 1.e-2
	)

Definition at line 22 of file FastAmplitudeIntegrator.cpp.

  : _db(false)
{
  //  initialise(pattern, amps, weight, generator, rnd, precision);
  initialise(pattern, amps, generator, rnd, precision);
}

FastAmplitudeIntegrator() [3/3]

FastAmplitudeIntegrator::FastAmplitudeIntegrator	(	const DalitzEventPattern &	pattern,
		IFastAmplitudeIntegrable *	amps,
		const std::string &	fname
	)

Definition at line 35 of file FastAmplitudeIntegrator.cpp.

  : _db(false)
{
  //  initialise(pattern, amps, weight, generator, rnd, precision);
  initialiseFromFile(pattern, amps, fname);
}

~FastAmplitudeIntegrator()

virtual FastAmplitudeIntegrator::~FastAmplitudeIntegrator ( )

inlinevirtual

Definition at line 83 of file FastAmplitudeIntegrator.h.

{}

Member Function Documentation

absSumOfInterferenceFractions()

double FastAmplitudeIntegrator::absSumOfInterferenceFractions ( )

inline

Definition at line 149 of file FastAmplitudeIntegrator.h.

                                         {
        return _integCalc->absSumOfInterferenceFractions();
  }     

absSumOfSqrtInterferenceFractions()

double FastAmplitudeIntegrator::absSumOfSqrtInterferenceFractions ( )

inline

Definition at line 145 of file FastAmplitudeIntegrator.h.

                                             {
        return _integCalc->absSumOfSqrtInterferenceFractions();
  }     

add()

bool FastAmplitudeIntegrator::add ( const FastAmplitudeIntegrator &  other )

virtual

Definition at line 114 of file FastAmplitudeIntegrator.cpp.

                                                                     {
  int totalEvents = _numEvents + other._numEvents;
  if(totalEvents > 0){
    _mean = (_numEvents * _mean + other._numEvents * other._mean)/totalEvents;
  }
  _numEvents = totalEvents;
  _Ncalls   += other._Ncalls;
  if(0 != this->_integCalc){
    _integCalc->add(other._integCalc);
  }else{
    if(0 != other._integCalc){
      counted_ptr<IntegCalculator> nL(new IntegCalculator(*(other._integCalc)));
      _integCalc = nL;
    }
  }
  if(0 != this->_integCalc_copyForDebug){
    _integCalc_copyForDebug->add(other._integCalc_copyForDebug);
  }else{
    if(0 != other._integCalc_copyForDebug){
      counted_ptr<IntegCalculator> 
    nLc(new IntegCalculator(*(other._integCalc_copyForDebug)));
      _integCalc_copyForDebug = nLc;
    }
  }
  _initialised |= other._initialised;
  _db |= other._db;

  // other quantities we take over if they aren't defined, here:
  if(0 == _amps) _amps = other._amps;
  if(0 == _rnd) _rnd = other._rnd;
  if(0 == _localRnd) _localRnd = other._localRnd;
  if(_precision <= 0.0) _precision = other._precision;

  return true;
}

addEvents()

int FastAmplitudeIntegrator::addEvents ( long int  Nevents )

protected

Definition at line 242 of file FastAmplitudeIntegrator.cpp.

                                                      {
  bool dbThis=false;
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::addEvents "
     << " need to get initialised first." << endl;
    return 0;
  }
  if(dbThis) cout << "FastAmplitudeIntegrator::addEvents "
          << " got called " << endl;
  if(Nevents <= 0) return _integCalc->numEvents();

  if(0 == _generator){
    cout << "FastAmplitudeIntegrator::addEvents "
     << " adding flat events. " << endl;
    for(int i=0; i < Nevents; i++){
      DalitzEvent evt(_pat, _rnd);
      _integCalc->addEvent(&evt, weight(&evt));
    }
  }else{
    int N_success = 0;
    unsigned long int maxTries = (unsigned long int) Nevents * 1000;
    // dealing with possible overflow:
    if(maxTries < (unsigned long int) Nevents) maxTries = Nevents * 100;
    if(maxTries < (unsigned long int) Nevents) maxTries = Nevents * 10;
    if(maxTries < (unsigned long int) Nevents) maxTries = Nevents;

    time_t startTime = time(0);
    for(unsigned long int i=0; i < maxTries && N_success < Nevents; i++){
      counted_ptr<IDalitzEvent> ptr(_generator->newEvent());
      if(dbThis) cout << "got event with ptr: " << ptr << endl;
      if(0 == ptr){
    if(_generator->exhausted()){
      cout << "ERROR in FastAmplitudeIntegrator::addEvents"
           << "\n\t Generator exhausted, cannot generate more events."
           << endl;
      break;
    }
      }else{
    _integCalc->addEvent(ptr, weight(ptr.get()));
    N_success++;
    int                    printEvery =   1000;
    if(N_success >  1000)  printEvery =  50000;
    if(N_success > 100000) printEvery = 100000;
    bool printout = (N_success%printEvery == 0);
    if(dbThis) printout |= (N_success < 20);

    if(printout){
      cout << " FastAmplitudeIntegrator::addEvents: added " << N_success 
           << " th event";
      evaluateSum();
      double v= variance();
      double sigma = -9999;
      if(v > 0) sigma = sqrt(v);
      cout << "\t integ= " << _mean << " +/- " << sigma;
      cout << "\t("
           << MINT::stringtime(difftime(time(0), startTime))
           << ")";
      cout << endl;
    }
      }
    }
  }
  if(_integCalc->numEvents() <  Nevents){
    cout << "ERROR in FastAmplitudeIntegrator::addEvents"
     << "\n\t> I was supposed to add " << Nevents << " events,"
     << "\n\t> but now the total list size is only " 
     << _integCalc->numEvents()
     << endl;
  }else{
    cout << " added " << Nevents << " _integCalc->numEvents "
     << _integCalc->numEvents() << endl;

  }
  return _integCalc->numEvents();
}

determineNumEvents()

int FastAmplitudeIntegrator::determineNumEvents ( )

protected

Definition at line 373 of file FastAmplitudeIntegrator.cpp.

                                               {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::determineNumEvents "
     << " need to get initialised first." << endl;
    return 0;
  }

  updateEventSet(_minEvents);
  evaluateSum();
  double v = variance();
  double safetyFactor = 1.15;
  double maxVariance = _precision*_precision * _mean*_mean;
  _numEvents = _integCalc->numEvents() * 
    (v*safetyFactor/maxVariance);
  // round to the nearest 100:
  _numEvents = (_numEvents/100 + 1)*100;

  cout << " FastAmplitudeIntegrator::determineNumEvents():"
       << "\n\t> mean: " << _mean << ", rms " << sqrt(v)
       << "\n\t> currently, our relative precision is " 
       << sqrt(v)/_mean
       << "\n\t> to achieve the requested pecision of " << _precision << ","
       << "\n\t> (with a safety factor) I think we'll need " 
       << (v*safetyFactor/maxVariance)
       << " times that, i.e. " << _numEvents << " events" << endl;

  return _numEvents;
}

doFinalStats()

void FastAmplitudeIntegrator::doFinalStats ( MINT::Minimiser *  mini = 0 )

virtual

Implements IDalitzIntegrator.

Definition at line 481 of file FastAmplitudeIntegrator.cpp.

                                                         {
  bool dbThis=true;
  if(dbThis){
    cout << "FastAmplitudeIntegrator::doFinalStats() called" << endl;
    cout << "Am I initialised? " << _initialised << endl;
    cout << "Got called with these pointers: mini " << mini
     << " and my own _integCalc " << _integCalc << endl;
  }
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::doFinalStats "
     << " need to get initialised first." << endl;
    cout << " So I won't do anything." << endl;
    return;
  }
    
  _integCalc->doFinalStats(mini);
}

ensureFreshEvents()

bool FastAmplitudeIntegrator::ensureFreshEvents

(

)

Definition at line 150 of file FastAmplitudeIntegrator.cpp.

                                               {
  _rnd = makeNewRnd(time(0));
  if(0 != _generator){
    _generator->ensureFreshEvents();
  }
  return true;
}

evaluateSum()

double FastAmplitudeIntegrator::evaluateSum ( )

protected

Definition at line 326 of file FastAmplitudeIntegrator.cpp.

                                           {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::evaluateSum "
     << " need to get initialised first." << endl;
    return 0;
  }

  time_t tstart = time(0);
  _Ncalls++;

  _mean    = _integCalc->integral();
  double v = variance();
  double sigma = -9999;
  if(v >= 0) sigma = sqrt(v);
  double actualPrecision = -9999;
  if(_mean != 0 && sigma > 0) actualPrecision = sigma/_mean;

  double delT = difftime(time(0), tstart);

  int                  printEvery =    100;
  if(_Ncalls >   1000) printEvery =   1000;
  if(_Ncalls >  10000) printEvery =  10000;
  if(_Ncalls > 100000) printEvery = 100000;
  bool printout = _Ncalls < 3 || (0 == _Ncalls%printEvery);
 
  if(printout || _db){
    cout << " FastAmplitudeIntegrator::evaluateSum() for "
     << _integCalc->numEvents() << " events, "
     << _Ncalls << " th call:"
     << "\n\t> getting: " << _mean << " +/- " << sigma
     << "\n\t> precision: requested: " << _precision*100
     << "%, actual: " << actualPrecision*100 << "%"
     << "\n\t> This took " << delT << " s."
     << endl;
    if(_db){
      cout << "checking save and read-back:" << endl;
      _integCalc->save("checkSaveAndReadBack");
      _integCalc_copyForDebug->retrieve("checkSaveAndReadBack");
      cout << " copy values:"
       << "\n\t     mean: " << _integCalc_copyForDebug->integral()
       << "\n\t variance: " << _integCalc_copyForDebug->variance()
       << endl;
    }
  }
  return _mean;
}

generateEnoughEvents()

int FastAmplitudeIntegrator::generateEnoughEvents ( )

protected

Definition at line 402 of file FastAmplitudeIntegrator.cpp.

                                                 {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::generateEnoughEvents "
     << " need to get initialised first." << endl;
    return 0;
  }

  updateEventSet(_minEvents);
  int counter(0);
  while(determineNumEvents() > _integCalc->numEvents() &&
    counter < 10 && _numEvents > 0){
    updateEventSet(_numEvents);
    counter++;
  }
  return _integCalc->numEvents();
}

GetEachAmpsHistograms()

std::vector< DalitzHistoSet > FastAmplitudeIntegrator::GetEachAmpsHistograms

(

)

Definition at line 442 of file FastAmplitudeIntegrator.cpp.

                                                                        {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::GetEachAmpsHistograms "
     << " need to get initialised first." << endl;
    std::vector<DalitzHistoSet> dummy;
    return dummy;
  }
  return _integCalc->GetEachAmpsHistograms();
}

getFractionChi2()

double FastAmplitudeIntegrator::getFractionChi2

(

)

const

Definition at line 498 of file FastAmplitudeIntegrator.cpp.

                                                    {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::getFractionChi2() "
     << " need to get initialised first." << endl;
    return -9999;
  }
  return _integCalc->getFractionChi2();
}

GetInterferenceHistograms()

std::vector< DalitzHistoSet > FastAmplitudeIntegrator::GetInterferenceHistograms

(

)

Definition at line 470 of file FastAmplitudeIntegrator.cpp.

                                                                            {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::GetInterferenceHistograms "
     << " need to get initialised first." << endl;
    std::vector<DalitzHistoSet> dummy;
    return dummy;
  }

  return _integCalc->GetInterferenceHistograms();
}

getVal()

double FastAmplitudeIntegrator::getVal

(

)

Definition at line 419 of file FastAmplitudeIntegrator.cpp.

                                      {
  evaluateSum();
  return _mean;
}

Gradient()

void FastAmplitudeIntegrator::Gradient ( MINT::MinuitParameterSet *  mps, std::vector< double > &  grad  )

inline

Definition at line 130 of file FastAmplitudeIntegrator.h.

                                                                     {
        _integCalc->Gradient(mps,grad);
  }

GradientForLasso()

void FastAmplitudeIntegrator::GradientForLasso ( MINT::MinuitParameterSet *  mps, std::vector< double > &  grad  )

inline

Definition at line 133 of file FastAmplitudeIntegrator.h.

                                                                             {
        _integCalc->GradientForLasso(mps,grad);
  }  

histoSet()

DalitzHistoSet FastAmplitudeIntegrator::histoSet

(

)

const

Definition at line 424 of file FastAmplitudeIntegrator.cpp.

                                                      {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::histoSet "
     << " need to get initialised first." << endl;
    throw "bummer";
  }

  return _integCalc->histoSet();
}

initialise() [1/2]

bool FastAmplitudeIntegrator::initialise	(	const DalitzEventPattern &	pattern,
		IFastAmplitudeIntegrable *	amps = 0,
		MINT::IEventGenerator< IDalitzEvent > *	eventGenerator = 0,
		TRandom *	rnd = 0,
		double	precision = 1.e-2
	)

Definition at line 52 of file FastAmplitudeIntegrator.cpp.

          {
  bool sc=true;
  sc &= setValues(pattern, amps, generator, rnd, precision);
  sc &= generateEnoughEvents();
  return sc;
}

initialise() [2/2]

bool FastAmplitudeIntegrator::initialise	(	const std::string &	commaSeparatedList,
		const DalitzEventPattern &	pattern,
		IFastAmplitudeIntegrable *	amps = 0,
		MINT::IEventGenerator< IDalitzEvent > *	eventGenerator = 0,
		TRandom *	rnd = 0,
		double	precision = 1.e-2
	)

Definition at line 84 of file FastAmplitudeIntegrator.cpp.

          {
  // with file and top up
  bool sc=true;
  sc &= setValues(pattern, amps, generator, rnd, precision);
  if(! sc) return false;
  std::string firstDir;
  std::stringstream is(commaSeparatedList);
  getline(is, firstDir, ',');

  struct stat buf;
  //  if(firstDir != "" && access(firstDir.c_str(), 0) == 0){
  if(firstDir != "" && stat(firstDir.c_str(), &buf) == 0){
    // it exists, we can retrieve it
    sc &= _integCalc->retrieve(commaSeparatedList);
    // if we add events to this, we need to avoid
    // duplication. The below will refuse to re-use
    // events read from a file, and will randomise
    // the random seed:
    this->ensureFreshEvents();
  }
  sc &= generateEnoughEvents();
  return sc;
}

initialised()

bool FastAmplitudeIntegrator::initialised ( ) const

inline

Definition at line 125 of file FastAmplitudeIntegrator.h.

                          {
    return _initialised;
  }

initialiseFromFile()

bool FastAmplitudeIntegrator::initialiseFromFile	(	const DalitzEventPattern &	pattern,
		IFastAmplitudeIntegrable *	amps,
		const std::string &	commaSeparatedList
	)

Definition at line 66 of file FastAmplitudeIntegrator.cpp.

              {
  bool sc = setValues(pattern, amps, 0, 0, 1.e-10);
  if(! sc) return false;
  sc &= _integCalc->retrieve(commaSeparatedList);
  return sc;
}

interferenceHistoSet()

DalitzHistoSet FastAmplitudeIntegrator::interferenceHistoSet

(

)

const

Definition at line 452 of file FastAmplitudeIntegrator.cpp.

                                                                  {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::interferenceHistoSet "
     << " need to get initialised first." << endl;
    DalitzHistoSet dummy;
    return dummy;
  }
  return _integCalc->interferenceHistoSet();
}

makeNewRnd()

TRandom * FastAmplitudeIntegrator::makeNewRnd ( int  seed = -9999 )

protected

Definition at line 158 of file FastAmplitudeIntegrator.cpp.

                                                    {
  if(seed < -9998) seed = time(0);
  cout << "FastAmplitudeIntegrator::makeNewRnd with seed " << seed << endl;
  counted_ptr<TRandom> lr(new TRandom3(seed));
  _localRnd = lr;
  _rnd = lr.get();
  return _rnd;
}

numberOfFitFractionsLargerThanThreshold()

int FastAmplitudeIntegrator::numberOfFitFractionsLargerThanThreshold ( double  threshold )

inline

Definition at line 157 of file FastAmplitudeIntegrator.h.

                                                               {
        return _integCalc->numberOfFitFractionsLargerThanThreshold(threshold);
  } 

RealVal()

double FastAmplitudeIntegrator::RealVal ( )

inlinevirtual

Implements MINT::IReturnReal.

Definition at line 137 of file FastAmplitudeIntegrator.h.

                  {
    return getVal();
  }

save()

bool FastAmplitudeIntegrator::save

(

const std::string &

fname

)

const

Definition at line 507 of file FastAmplitudeIntegrator.cpp.

                                                             {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::save( " << fname << " ):"
     << " need to get initialised first." << endl;
    cout << " So I won't do anything." << endl;
    return false;
  }
  return _integCalc->save(fname);
}

saveEachAmpsHistograms()

void FastAmplitudeIntegrator::saveEachAmpsHistograms

(

const std::string &

prefix

)

const

Definition at line 433 of file FastAmplitudeIntegrator.cpp.

                                                                                 {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::saveEachAmpsHistograms "
     << " need to get initialised first." << endl;
    cout << "So I won't do anything" << endl;
    return;
  }
  return _integCalc->saveEachAmpsHistograms(prefix);
}

saveInterferenceHistograms()

void FastAmplitudeIntegrator::saveInterferenceHistograms

(

const std::string &

prefix

)

const

Definition at line 461 of file FastAmplitudeIntegrator.cpp.

                                                                                     {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::saveInterferenceAmpsHistograms "
     << " need to get initialised first." << endl;
    cout << "So I won't do anything" << endl;
    return;
  }
  return _integCalc->saveInterferenceHistograms(prefix);
}

setPrecision()

void FastAmplitudeIntegrator::setPrecision ( double  prec )

inline

Definition at line 115 of file FastAmplitudeIntegrator.h.

                                {
    _precision=prec;
    if(initialised()){
      generateEnoughEvents();
    }
  }

setValues()

bool FastAmplitudeIntegrator::setValues	(	const DalitzEventPattern &	pattern,
		IFastAmplitudeIntegrable *	amps = 0,
		MINT::IEventGenerator< IDalitzEvent > *	eventGenerator = 0,
		TRandom *	rnd = 0,
		double	precision = 1.e-2
	)

Definition at line 167 of file FastAmplitudeIntegrator.cpp.

         {
  bool dbThis=false;
  _pat  = pattern;
  _amps = amps;
  //_weight = weight;
  if(0 != rnd){
    _rnd = rnd;
  }else{
    _rnd = makeNewRnd(time(0));
  }
  _precision = precision;
  
  _generator = generator;

  _integCalc = amps->makeIntegCalculator();
  if(dbThis){
    cout << "FastAmplitudeIntegrator::setValues: Have _integCalc ptr " << _integCalc
     << ", now printing it." << endl;
    _integCalc->print();
  }

  if(_db) _integCalc_copyForDebug = amps->makeIntegCalculator();

  if(0 == _integCalc){
    cout << "ERROR in FastAmplitudeIntegrator::initialise"
     << " amps->makeIntegCalculator()"
     << " return NULL pointer!" << endl;
    cout << "\t This will go wrong." << endl;
  }
  _Ncalls=0;

  _initialised = true;
  return _initialised;
}

sumOfFitFractions()

double FastAmplitudeIntegrator::sumOfFitFractions ( )

inline

Definition at line 141 of file FastAmplitudeIntegrator.h.

                             {
        return _integCalc->sumOfFitFractions();
  }   

sumOfSqrtFitFractions()

double FastAmplitudeIntegrator::sumOfSqrtFitFractions ( )

inline

Definition at line 153 of file FastAmplitudeIntegrator.h.

                                 {
        return _integCalc->sumOfSqrtFitFractions();
  }     

updateEventSet()

int FastAmplitudeIntegrator::updateEventSet ( long int  Nevents )

protected

Definition at line 208 of file FastAmplitudeIntegrator.cpp.

                                                           {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::updateEventSet "
     << " need to know pattern first." << endl;
    return 0;
  }
  int missingEvents = Nevents - _integCalc->numEvents();

  if(missingEvents > 0){
    cout << "missing events: " << missingEvents 
     << " = " << Nevents 
     << " - " << _integCalc->numEvents()
     << endl;
    addEvents(missingEvents);
  }else{
    cout << "want " << Nevents
     << ", have " << _integCalc->numEvents()
     << ", so no more events missing." << endl;
  }
  return _integCalc->numEvents();
}

variance()

double FastAmplitudeIntegrator::variance

(

)

const

Definition at line 318 of file FastAmplitudeIntegrator.cpp.

                                              {
  if(! _initialised){
    cout << " FastAmplitudeIntegrator::variance "
     << " need to get initialised first." << endl;
    return 0;
  }
  return _integCalc->variance();
}

weight()

double FastAmplitudeIntegrator::weight ( IDalitzEvent *  evtPtr )

protected

Definition at line 230 of file FastAmplitudeIntegrator.cpp.

                                                    {
  return 1;
  // (not in use, but left as a hook should the need arise)
  /*
  if(0 == _weight) return 1;

  _weight->setEvent(evtPtr);
  double wght = _weight->RealVal();
  _weight->resetEventRecord();
  return wght;
  */
}

Member Data Documentation

_amps

IFastAmplitudeIntegrable* FastAmplitudeIntegrator::_amps

protected

Definition at line 44 of file FastAmplitudeIntegrator.h.

_db

bool FastAmplitudeIntegrator::_db

private

Definition at line 40 of file FastAmplitudeIntegrator.h.

_generator

MINT::IEventGenerator<IDalitzEvent>* FastAmplitudeIntegrator::_generator

protected

Definition at line 65 of file FastAmplitudeIntegrator.h.

_initialised

bool FastAmplitudeIntegrator::_initialised

private

Definition at line 38 of file FastAmplitudeIntegrator.h.

_integCalc

MINT::counted_ptr<IntegCalculator> FastAmplitudeIntegrator::_integCalc

protected

Definition at line 49 of file FastAmplitudeIntegrator.h.

_integCalc_copyForDebug

MINT::counted_ptr<IntegCalculator> FastAmplitudeIntegrator::_integCalc_copyForDebug

protected

Definition at line 50 of file FastAmplitudeIntegrator.h.

_localRnd

MINT::counted_ptr<TRandom> FastAmplitudeIntegrator::_localRnd

protected

Definition at line 53 of file FastAmplitudeIntegrator.h.

_mean

double FastAmplitudeIntegrator::_mean

mutableprivate

Definition at line 33 of file FastAmplitudeIntegrator.h.

_minEvents

const long int FastAmplitudeIntegrator::_minEvents =100

staticprivate

Definition at line 34 of file FastAmplitudeIntegrator.h.

_Ncalls

int FastAmplitudeIntegrator::_Ncalls

private

Definition at line 36 of file FastAmplitudeIntegrator.h.

_numEvents

int FastAmplitudeIntegrator::_numEvents

private

Definition at line 35 of file FastAmplitudeIntegrator.h.

_pat

DalitzEventPattern FastAmplitudeIntegrator::_pat

protected

Definition at line 43 of file FastAmplitudeIntegrator.h.

_precision

double FastAmplitudeIntegrator::_precision

protected

Definition at line 54 of file FastAmplitudeIntegrator.h.

_rnd

TRandom* FastAmplitudeIntegrator::_rnd

protected

Definition at line 52 of file FastAmplitudeIntegrator.h.

---

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

• Mint/FastAmplitudeIntegrator.h
• src/Mojito/DalitzIntegrator/FastAmplitudeIntegrator.cpp