#include <DalitzBoxSet_Method2.h>

Inheritance diagram for DalitzBoxSet_Method2:

Public Member Functions
	DalitzBoxSet_Method2 (const DalitzEventPattern &pat, TRandom *rnd=gRandom)

	DalitzBoxSet_Method2 (const DalitzBoxSet &boxSet, const DalitzEventPattern &pat)

	DalitzBoxSet_Method2 (const DalitzBoxSet_Method2 &other)

double	maxHeight () const

virtual MINT::counted_ptr< DalitzEvent >	tryEvent ()

virtual MINT::counted_ptr< DalitzEvent >	tryWeightedEvent ()

virtual MINT::counted_ptr< IDalitzEvent >	newEvent ()

virtual bool	exhausted () const

virtual	~DalitzBoxSet_Method2 ()

Public Member Functions inherited from DalitzBoxSet
void	makeRelations ()

void	removeDuplicates ()

void	sortYourself ()

void	callSortYourselfForDebug ()

	DalitzBoxSet (TRandom *rnd=gRandom)

	DalitzBoxSet (const DalitzBoxSet &other)

void	add (const DalitzBox &box)

void	add (const MINT::PolymorphVector< DalitzBox > &boxes)

void	add (const DalitzBoxSet &boxes)

void	setPDF (MINT::IReturnRealForEvent< IDalitzEvent > *amps)

DalitzBoxSet	split (unsigned int nWays) const

DalitzBoxSet	splitIfWiderThan (double maxWidth) const

virtual MINT::counted_ptr< DalitzEvent >	generateEventForOwner ()

virtual DalitzEvent	generateEvent ()

virtual void	print (std::ostream &os=std::cout) const

bool	setRnd (TRandom *rnd)

virtual bool	ensureFreshEvents ()

virtual	~DalitzBoxSet ()

Public Member Functions inherited from MINT::PolymorphVector< DalitzBox >
	PolymorphVector ()

	PolymorphVector (unsigned int N)

	PolymorphVector (unsigned int N, const DalitzBox &c)

	PolymorphVector (const PolymorphVector &other)

	PolymorphVector (const typename std::vector< DalitzBox > &other)

virtual	~PolymorphVector ()

std::vector< DalitzBox > &	theVector ()

const std::vector< DalitzBox > &	theVector () const

DalitzBox &	operator[] (unsigned int i)

const DalitzBox &	operator[] (unsigned int i) const

DalitzBox &	at (unsigned int i)

const DalitzBox &	at (unsigned int i) const

std::vector< DalitzBox >::iterator	begin ()

std::vector< DalitzBox >::const_iterator	begin () const

std::vector< DalitzBox >::iterator	end ()

std::vector< DalitzBox >::const_iterator	end () const

std::vector< DalitzBox >::iterator	find (const DalitzBox &c)

std::vector< DalitzBox >::const_iterator	find (const DalitzBox &c) const

DalitzBox &	front ()

const DalitzBox &	front () const

DalitzBox &	back ()

const DalitzBox &	back () const

unsigned int	size () const

bool	empty () const

void	push_back (const DalitzBox &c)

void	pop_back ()

void	erase (typename std::vector< DalitzBox >::iterator pos)

void	erase (typename std::vector< DalitzBox >::iterator first, typename std::vector< DalitzBox >::iterator last)

PolymorphVector< DalitzBox > &	operator= (const PolymorphVector< DalitzBox > &other)

void	clear ()

void	resize (unsigned int N)

void	resize (unsigned int N, const DalitzBox &c)

	operator const typename std::vector< DalitzBox > & () const

	operator typename std::vector< DalitzBox > & ()

bool	operator== (const MINT::PolymorphVector< DalitzBox > &v2) const

bool	operator!= (const MINT::PolymorphVector< DalitzBox > &v2) const

bool	operator< (const MINT::PolymorphVector< DalitzBox > &v2) const

bool	operator> (const MINT::PolymorphVector< DalitzBox > &v2) const

Public Member Functions inherited from MINT::IEventGenerator< IDalitzEvent >
virtual	~IEventGenerator ()

Private Member Functions
void	setBoxesToFlat ()

void	getMaxHeight ()

void	setEmptyBoxesHeight ()

virtual void	getReady ()

double	eventsPDF (DalitzEvent &evt)

Private Attributes
DalitzEventPattern	_pat

double	_maxHeight

Additional Inherited Members
Protected Member Functions inherited from DalitzBoxSet
double	VolumeSum () const

double	AreaSum () const

void	makeVolumeProbIntervals ()

int	pickRandomVolume ()

void	setBoxPDFs ()

void	setBoxRnds ()

void	getBoxesReady ()

Protected Attributes inherited from DalitzBoxSet
bool	_ready

MINT::PolymorphVector< double >	_volumeProbs

MINT::IReturnRealForEvent< IDalitzEvent > *	_amps

TRandom *	_rnd

Protected Attributes inherited from MINT::PolymorphVector< DalitzBox >
std::vector< DalitzBox >	_vec

Detailed Description

Definition at line 12 of file DalitzBoxSet_Method2.h.

Constructor & Destructor Documentation

◆ DalitzBoxSet_Method2() [1/3]

DalitzBoxSet_Method2::DalitzBoxSet_Method2	(	const DalitzEventPattern &	pat,
		TRandom *	rnd = `gRandom`
	)

Definition at line 16 of file DalitzBoxSet_Method2.cpp.

   : DalitzBoxSet(rnd)
   , _pat(pat)
 {
 }

◆ DalitzBoxSet_Method2() [2/3]

DalitzBoxSet_Method2::DalitzBoxSet_Method2	(	const DalitzBoxSet &	boxSet,
		const DalitzEventPattern &	pat
	)

Definition at line 23 of file DalitzBoxSet_Method2.cpp.

   : DalitzBoxSet(boxSet)
   , _pat(pat)
 {
   cout << " DalitzBoxSet_Method2 from DalitzBoxSet, and pattern " 
        << pat
        << ", ready? " << _ready
        << "; this is me:\n" << (*this)
        << endl;
   _ready = false;
 }

◆ DalitzBoxSet_Method2() [3/3]

DalitzBoxSet_Method2::DalitzBoxSet_Method2 ( const DalitzBoxSet_Method2 & other )

Definition at line 9 of file DalitzBoxSet_Method2.cpp.

   : IEventGenerator<IDalitzEvent>()
   , DalitzBoxSet(other)
   , _pat(other._pat)
 {
 }

◆ ~DalitzBoxSet_Method2()

virtual DalitzBoxSet_Method2::~DalitzBoxSet_Method2 ( )

inlinevirtual

Definition at line 45 of file DalitzBoxSet_Method2.h.

45 {}

Member Function Documentation

◆ eventsPDF()

double DalitzBoxSet_Method2::eventsPDF ( DalitzEvent & evt )

private

Definition at line 74 of file DalitzBoxSet_Method2.cpp.

                                                       {
 
   if(0 == _amps ) return 0;
   double val = _amps->RealVal(evt);
   return val;
 }

◆ exhausted()

virtual bool DalitzBoxSet_Method2::exhausted ( ) const

inlinevirtual

Reimplemented from DalitzBoxSet.

Definition at line 44 of file DalitzBoxSet_Method2.h.

44 {return false;}

◆ getMaxHeight()

void DalitzBoxSet_Method2::getMaxHeight ( )

private

Definition at line 44 of file DalitzBoxSet_Method2.cpp.

                                        {
   _maxHeight=-1;
   for(unsigned int i=0; i < this->size(); i++){
     double thisHeight = (*this)[i].height();
     if(thisHeight > _maxHeight) _maxHeight = thisHeight;
   }
 }

◆ getReady()

void DalitzBoxSet_Method2::getReady ( )

privatevirtual

Reimplemented from DalitzBoxSet.

Definition at line 65 of file DalitzBoxSet_Method2.cpp.

                                    {
   cout << " getting ready " << endl;
   setBoxesToFlat();
   DalitzBoxSet::getReady();
   getMaxHeight();
   setEmptyBoxesHeight();
   cout << "DalitzBoxSet_Method2: I am ready.\n" << (*this) << endl;
 }

◆ maxHeight()

double DalitzBoxSet_Method2::maxHeight ( ) const

Definition at line 61 of file DalitzBoxSet_Method2.cpp.

                                             {
   return _maxHeight;
 }

◆ newEvent()

counted_ptr< IDalitzEvent > DalitzBoxSet_Method2::newEvent ( )

virtual

Reimplemented from DalitzBoxSet.

Definition at line 88 of file DalitzBoxSet_Method2.cpp.

                                                         {
   bool dbThis=false;
   if(dbThis)cout << "_ready = " << _ready << endl;
   if(dbThis)cout << " maxHeight " << maxHeight() << endl;
   if(maxHeight() < -9998 || ! _ready) getReady();
   if(dbThis)cout << " DalitzBoxSet_Method2::newEvent "
          << " getting event " << endl;
   counted_ptr<DalitzEvent> ptr(0);
   while(0 == ptr){
     ptr = tryWeightedEvent();
   }
   if(dbThis) cout << "got event, returning " << ptr << endl;
   return ptr;
 }

◆ setBoxesToFlat()

void DalitzBoxSet_Method2::setBoxesToFlat ( )

private

Definition at line 38 of file DalitzBoxSet_Method2.cpp.

                                          {
   for(unsigned int i=0; i < this->size(); i++){
     (*this)[i].setFlatPhaseSpace();
   }
 }

◆ setEmptyBoxesHeight()

void DalitzBoxSet_Method2::setEmptyBoxesHeight ( )

private

Definition at line 52 of file DalitzBoxSet_Method2.cpp.

                                               {
   for(unsigned int i=0; i < this->size(); i++){
     double thisHeight = (*this)[i].height();
     if(thisHeight < 0){
        (*this)[i].setHeight(_maxHeight);
     }
   }
 }

◆ tryEvent()

counted_ptr< DalitzEvent > DalitzBoxSet_Method2::tryEvent ( )

virtual

Reimplemented from DalitzBoxSet.

Definition at line 103 of file DalitzBoxSet_Method2.cpp.

                                                        {
   if(! _ready) getReady();
   counted_ptr<DalitzEvent> evtPtr(new DalitzEvent(_pat, _rnd));
   int boxIndex=-1;
   /*
   cout << "Hello from DalitzBoxSet_Method2::tryEvent()"
        << endl;
   cout << "I've got an event to try " << *evtPtr << endl;
   */
   for(unsigned int i=0; i< this->size(); i++){
     bool isIn = (*this)[i].insideArea(*evtPtr);
     /*
       cout << "Is event in box " << i << " ? " << endl;
       cout << "\t" << (isIn? "yes" : "no")
        << endl;
     */
     if(isIn){
       boxIndex=i;
       break;
     }
   }
   if(boxIndex < 0 || boxIndex >= (int) this->size()){
     cout << "ERROR in DalitzBoxSet_Method2::tryEvent()"
      << " boxIndex = " << boxIndex
      << ", outside [" << 0 << ", " << this->size() << "]"
      << endl;
     return counted_ptr<DalitzEvent>(0);
   }
   
   double thisHeight = (*this)[boxIndex].height();
   if(_rnd->Rndm()*maxHeight() >= thisHeight){
     return counted_ptr<DalitzEvent>(0);
   }
   // OK, so it passed the box stage, now internal:
   double pdf = eventsPDF(*evtPtr);
   if(pdf > thisHeight){
     cout << "ERROR, underestimated box height [" << boxIndex << "]"
      << "\n\t\t thisHeight = " << thisHeight 
      << "\n\t\t      < pdf = " << pdf
      << "\n\t\t (maxHeight = " << maxHeight() << ")" << endl;
     (*this)[boxIndex].setHeight(3*pdf);
     getMaxHeight();
   }
   if(_rnd->Rndm()* thisHeight >= pdf) {
     return counted_ptr<DalitzEvent>(0);
   }
 
   return evtPtr;
 }

◆ tryWeightedEvent()

counted_ptr< DalitzEvent > DalitzBoxSet_Method2::tryWeightedEvent ( )

virtual

Reimplemented from DalitzBoxSet.

Definition at line 81 of file DalitzBoxSet_Method2.cpp.

                                                                {
   if(! _ready) getReady();
   int boxIndex = pickRandomVolume();
   //cout << "DalitzBoxSet::tryEvent() in volume " << boxIndex << endl;
   return (*this)[boxIndex].tryWeightedEventForOwner();
 }

Member Data Documentation

◆ _maxHeight

double DalitzBoxSet_Method2::_maxHeight

private

Definition at line 19 of file DalitzBoxSet_Method2.h.

◆ _pat

DalitzEventPattern DalitzBoxSet_Method2::_pat

private

Definition at line 17 of file DalitzBoxSet_Method2.h.

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

Mint/DalitzBoxSet_Method2.h
src/Mojito/BetterMC/DalitzBoxSet_Method2.cpp

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DalitzBoxSet_Method2() [1/3]

◆ DalitzBoxSet_Method2() [2/3]

◆ DalitzBoxSet_Method2() [3/3]

◆ ~DalitzBoxSet_Method2()

Member Function Documentation

◆ eventsPDF()

◆ exhausted()

◆ getMaxHeight()

◆ getReady()

◆ maxHeight()

◆ newEvent()

◆ setBoxesToFlat()

◆ setEmptyBoxesHeight()

◆ tryEvent()

◆ tryWeightedEvent()

Member Data Documentation

◆ _maxHeight

◆ _pat