DalitzMCIntegrator Class Reference

#include <DalitzMCIntegrator.h>

Inheritance diagram for DalitzMCIntegrator:

Classes
class	integrationWeight

Public Member Functions
	DalitzMCIntegrator (const DalitzEventPattern &pattern, MINT::IReturnRealForEvent< IDalitzEvent > *weightFunction=0, MINT::IEventGenerator< IDalitzEvent > *eventGenerator=0, TRandom *rnd=gRandom, double precision=1.e-2)
	DalitzMCIntegrator ()
bool	initialise (const DalitzEventPattern &pattern, MINT::IReturnRealForEvent< IDalitzEvent > *weightFunction=0, MINT::IEventGenerator< IDalitzEvent > *eventGenerator=0, TRandom *rnd=gRandom, double precision=1.e-2)
bool	resetIntegrand (MINT::IReturnRealForEvent< IDalitzEvent > *weightFunction=0)
void	setPrecision (double prec)
bool	initialised () const
double	getVal ()
double	RealVal ()
MINT::MinuitParameterSet *	getMPS ()
const MINT::MinuitParameterSet *	getMPS () const
virtual void	doFinalStats (MINT::Minimiser *mini=0)
DalitzHistoSet	histoSet ()
virtual	~DalitzMCIntegrator ()
Public Member Functions inherited from MINT::IReturnReal
virtual	~IReturnReal ()

Protected Member Functions
int	updateEventSet (int Nevents)
double	evaluateSum ()
int	determineNumEvents ()
int	generateEnoughEvents ()
int	addEvents (int Nevents)
Protected Member Functions inherited from MINT::IReturnReal
	IReturnReal ()

Protected Attributes
MINT::MinuitParameterSet *	_mps
DalitzEventPattern	_pat
MINT::IReturnRealForEvent< IDalitzEvent > *	_w
integrationWeight	_iw
DalitzEventPtrList	_eventPtrList
TRandom *	_rnd
double	_precision
MINT::IEventGenerator< IDalitzEvent > *	_generator

Private Attributes
long int	_Ncalls
double	_sumT
double	_mean
double	_variance
int	_numEvents
double	_weightSum
bool	_initialised

Static Private Attributes
static const int	_minEvents =1000

Detailed Description

Definition at line 28 of file DalitzMCIntegrator.h.

Constructor & Destructor Documentation

DalitzMCIntegrator() [1/2]

DalitzMCIntegrator::DalitzMCIntegrator	(	const DalitzEventPattern &	pattern,
		MINT::IReturnRealForEvent< IDalitzEvent > *	weightFunction = 0,
		MINT::IEventGenerator< IDalitzEvent > *	eventGenerator = 0,
		TRandom *	rnd = gRandom,
		double	precision = 1.e-2
	)

Definition at line 20 of file DalitzMCIntegrator.cpp.

  : _Ncalls(0)
  , _sumT(0)
  , _mean(0)
  , _variance(0)
  , _numEvents(0)
  , _weightSum(0)
  , _initialised(0)
  , _iw(weightFunction)
{
  initialise(pattern, weightFunction, generator, rnd, precision);
}

DalitzMCIntegrator() [2/2]

DalitzMCIntegrator::DalitzMCIntegrator

(

)

Definition at line 37 of file DalitzMCIntegrator.cpp.

  : _Ncalls(0)
  , _sumT(0)
  , _mean(0)
  , _variance(0)
  , _numEvents(0)
  , _weightSum(0)
  , _initialised(0)
  , _iw(0)
{
}

~DalitzMCIntegrator()

virtual DalitzMCIntegrator::~DalitzMCIntegrator ( )

inlinevirtual

Definition at line 116 of file DalitzMCIntegrator.h.

{};

Member Function Documentation

addEvents()

int DalitzMCIntegrator::addEvents ( int  Nevents )

protected

Definition at line 99 of file DalitzMCIntegrator.cpp.

                                            {
  bool dbThis=false;
  if(Nevents <= 0) return _eventPtrList.size();

  if(0 == _generator){
    _eventPtrList.generatePhaseSpaceEvents(Nevents, _pat, _rnd);
    return _eventPtrList.size();
  }
  
  int N_success = 0;
  for(int i=0; i < Nevents*1000 && N_success < Nevents; i++){
    //counted_ptr<IDalitzEvent> ptr(_generator->newEvent());
    counted_ptr<DalitzEvent> ptr(new DalitzEvent(_generator->newEvent().get()));
    if(dbThis) cout << "got event with ptr: " << ptr << endl;
    if(0 != ptr){
      _eventPtrList.push_back(ptr); // change event list to list of counted ptrs
      N_success++;
      int printEvery = 1000;
      if(N_success > 10000) printEvery=10000;
      bool printout = (N_success%printEvery == 0);
      if(printout){
    cout << " DalitzMCIntegrator::addEvents: added " << N_success 
         << " th event" << endl;
      }
    }
  }
  cout << " added " << Nevents << " _eventList.size() "
       << _eventPtrList.size() << endl;
  return _eventPtrList.size();
}

determineNumEvents()

int DalitzMCIntegrator::determineNumEvents ( )

protected

Definition at line 213 of file DalitzMCIntegrator.cpp.

                                          {
  updateEventSet(_minEvents);
  evaluateSum();

  double safetyFactor = 1.15;
  double maxVariance = _precision*_precision * _mean*_mean;
  _numEvents = _eventPtrList.size() * ( ((int)(_variance*safetyFactor/maxVariance)) + 1);
  cout << " DalitzMCIntegrator::determineNumEvents():"
       << "\n    > mean: " << _mean << ", rms " << sqrt(_variance)
       << "\n    > currently, our precision is " << sqrt(_variance)/_mean
       << "\n    > to achieve requested pecision of " << _precision << ","
       << "\n    > I think we'll need " << (_variance*safetyFactor/maxVariance)
       << " times that, i.e. " << _numEvents << " events" << endl;

  return _numEvents;
}

doFinalStats()

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

virtual

Implements IDalitzIntegrator.

Definition at line 265 of file DalitzMCIntegrator.cpp.

{}

evaluateSum()

double DalitzMCIntegrator::evaluateSum ( )

protected

Definition at line 130 of file DalitzMCIntegrator.cpp.

                                      {
  /*
  if(! _initialised){
    cout << " DalitzMCIntegrator::evaluateSum "
     << " need to know pattern first." << endl;
    return 0;
  }

  if(_eventPtrList.empty()){
    cout << "WARNING in DalitzMCIntegrator::evaluateSum()"
     << " no events!" << endl;
    _mean = _variance = 0;
    return 0;
  }
  */

  double sum   = 0;
  double sumsq = 0;
  double ws  =0;
  _weightSum = 0;
  //int printEveryNEvents = (_eventPtrList.size()/4);
  //if(printEveryNEvents < 1) printEveryNEvents = 1;

  //time_t tstart = time(0);

  //#pragma omp parallel for reduction(+:sum, sumsq, ws);
  for(unsigned int N=0; N < _eventPtrList.size(); N++){
    //DalitzEvent thisEvt(_eventPtrList[N]);
    /*  
    double ps = thisEvt.phaseSpace();
    if(ps <= 0.0){
      cout << "WARNING in DalitzMCIntegrator::evaluateSum()"
       << " event with phase space = " << ps << endl;
      continue; // should not happen.
    }
    */
    double weight =  _eventPtrList.getEventRef(N).getWeight() / _eventPtrList.getEventRef(N).getGeneratorPdfRelativeToPhaseSpace();
    double val = _iw.RealVal(_eventPtrList.getEventRef(N)); // _w->RealVal() * weight;

    sum   += val;
    sumsq += val*val;

    ws += weight;// /ps;

    /*  
    if(N%printEveryNEvents == 0){
      cout << "DalitzMCIntegrator::evaluateSum() N = " << N 
       << " val = " << val 
       << " _w->RealVal(thisEvt)*weight " << _w->RealVal(thisEvt)*weight 
       << endl;
    }
    */
  }

  _weightSum = ws;

  double fN       = (double) _eventPtrList.size();
  _mean           = sum   / fN;
  double meanOfSq = sumsq / fN;
  _variance       = (meanOfSq - _mean * _mean)/fN;

  if(_weightSum > 0){
    double sf = fN/_weightSum;
    _mean *= sf;
    _variance *= sf*sf;
  }

  /*
  double delT = difftime(time(0), tstart);
  double sigma = -9999;
  if(_variance >= 0) sigma = sqrt(_variance);
  double actualPrecision = -9999;
  if(_mean != 0 && sigma > 0) actualPrecision = sigma/_mean;
  cout << " DalitzMCIntegrator::evaluateSum() for " << N << " events:"
       << "\n     > getting: " << _mean << " +/- " << sigma
       << "\n     > precision: requested: " << _precision*100
       << "%, actual: " << actualPrecision*100 << "%"
       << "\n     > This took " << delT << " s."
       << endl;
  */
  return _mean;
}

generateEnoughEvents()

int DalitzMCIntegrator::generateEnoughEvents ( )

protected

Definition at line 230 of file DalitzMCIntegrator.cpp.

                                            {
  addEvents(300000);
  //updateEventSet(_minEvents);
  //determineNumEvents();
  //updateEventSet(_numEvents);
  return _eventPtrList.size();
}

getMPS() [1/2]

MINT::MinuitParameterSet* DalitzMCIntegrator::getMPS

(

)

getMPS() [2/2]

const MINT::MinuitParameterSet* DalitzMCIntegrator::getMPS

(

)

const

getVal()

double DalitzMCIntegrator::getVal

(

)

Definition at line 238 of file DalitzMCIntegrator.cpp.

                                 {
  _Ncalls++;

  int                  printEvery =     10;
  if(_Ncalls >    100) printEvery =    100;
  if(_Ncalls >   1000) printEvery =   1000;
  if(_Ncalls >  10000) printEvery =  10000;
  if(_Ncalls > 100000) printEvery = 100000;
  bool printout = _Ncalls < 10 || (0 == _Ncalls%printEvery);

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

  _sumT += delT;

  if(_Ncalls > 0 && printout){
    cout << "DalitzMCIntegrator::getVal, " << _Ncalls << "th call: Returning " << _mean
     << " this call took " << delT << " s,"
     << " average call " << _sumT/_Ncalls << " s, "
     << " for " << _eventPtrList.size() << " events."
     << endl;
  }

  return _mean;
}

histoSet()

DalitzHistoSet DalitzMCIntegrator::histoSet ( )

inline

Definition at line 109 of file DalitzMCIntegrator.h.

                           {
        double den= getVal();
        den *= _weightSum;
        if(den <= 0) den=1;
        return _eventPtrList.reWeightedHistoSet(&_iw)/den;
  }

initialise()

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

Definition at line 49 of file DalitzMCIntegrator.cpp.

                     {
  _Ncalls = 0;
  _sumT   = 0.0;
  _pat    = pattern;
  _w      = weightFunction;
  _rnd    = rnd;
  _precision = precision;
  
  _generator = generator;

  _iw.setWeight(weightFunction);

  _initialised = true;


  generateEnoughEvents();// must come last.

  return _initialised;
}

initialised()

bool DalitzMCIntegrator::initialised ( ) const

inline

Definition at line 94 of file DalitzMCIntegrator.h.

                          {
    return _initialised;
  }

RealVal()

double DalitzMCIntegrator::RealVal ( )

inlinevirtual

Implements MINT::IReturnReal.

Definition at line 100 of file DalitzMCIntegrator.h.

                  {
    return getVal();
  }

resetIntegrand()

bool DalitzMCIntegrator::resetIntegrand

(

MINT::IReturnRealForEvent< IDalitzEvent > *

weightFunction = 0

)

Definition at line 73 of file DalitzMCIntegrator.cpp.

                     {
  _w   = weightFunction;
  _iw.setWeight(weightFunction);
  return true;
}

setPrecision()

void DalitzMCIntegrator::setPrecision ( double  prec )

inline

Definition at line 84 of file DalitzMCIntegrator.h.

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

updateEventSet()

int DalitzMCIntegrator::updateEventSet ( int  Nevents )

protected

Definition at line 81 of file DalitzMCIntegrator.cpp.

                                                 {
  if(! _initialised){
    cout << " DalitzMCIntegrator::updateEventSet "
     << " need to know pattern first." << endl;
    return 0;

  }
  int missingEvents = Nevents - _eventPtrList.size();
  cout << "missing events: " << missingEvents 
       << " = " << Nevents 
       << " - " << _eventPtrList.size()
       << endl;
  if(missingEvents > 0){
    addEvents(missingEvents);
  }
  return _eventPtrList.size();
}

Member Data Documentation

_eventPtrList

DalitzEventPtrList DalitzMCIntegrator::_eventPtrList

protected

Definition at line 53 of file DalitzMCIntegrator.h.

_generator

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

protected

Definition at line 63 of file DalitzMCIntegrator.h.

_initialised

bool DalitzMCIntegrator::_initialised

private

Definition at line 47 of file DalitzMCIntegrator.h.

_iw

integrationWeight DalitzMCIntegrator::_iw

protected

Definition at line 52 of file DalitzMCIntegrator.h.

_mean

double DalitzMCIntegrator::_mean

mutableprivate

Definition at line 42 of file DalitzMCIntegrator.h.

_minEvents

const int DalitzMCIntegrator::_minEvents =1000

staticprivate

Definition at line 43 of file DalitzMCIntegrator.h.

_mps

MINT::MinuitParameterSet* DalitzMCIntegrator::_mps

protected

Definition at line 49 of file DalitzMCIntegrator.h.

_Ncalls

long int DalitzMCIntegrator::_Ncalls

private

Definition at line 39 of file DalitzMCIntegrator.h.

_numEvents

int DalitzMCIntegrator::_numEvents

private

Definition at line 44 of file DalitzMCIntegrator.h.

_pat

DalitzEventPattern DalitzMCIntegrator::_pat

protected

Definition at line 50 of file DalitzMCIntegrator.h.

_precision

double DalitzMCIntegrator::_precision

protected

Definition at line 55 of file DalitzMCIntegrator.h.

_rnd

TRandom* DalitzMCIntegrator::_rnd

protected

Definition at line 54 of file DalitzMCIntegrator.h.

_sumT

double DalitzMCIntegrator::_sumT

private

Definition at line 40 of file DalitzMCIntegrator.h.

_variance

double DalitzMCIntegrator::_variance

mutableprivate

Definition at line 42 of file DalitzMCIntegrator.h.

_w

MINT::IReturnRealForEvent<IDalitzEvent>* DalitzMCIntegrator::_w

protected

Definition at line 51 of file DalitzMCIntegrator.h.

_weightSum

double DalitzMCIntegrator::_weightSum

private

Definition at line 45 of file DalitzMCIntegrator.h.

---

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

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