DalitzBoxSet Class Reference

#include <DalitzBoxSet.h>

Inheritance diagram for DalitzBoxSet:

Public Member Functions
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 >	tryEvent ()
virtual MINT::counted_ptr< DalitzEvent >	tryWeightedEvent ()
virtual MINT::counted_ptr< DalitzEvent >	generateEventForOwner ()
virtual DalitzEvent	generateEvent ()
virtual MINT::counted_ptr< IDalitzEvent >	newEvent ()
virtual void	print (std::ostream &os=std::cout) const
virtual bool	exhausted () 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 ()

Protected Member Functions
double	VolumeSum () const
double	AreaSum () const
void	makeVolumeProbIntervals ()
int	pickRandomVolume ()
void	setBoxPDFs ()
void	setBoxRnds ()
void	getBoxesReady ()
virtual void	getReady ()

Protected Attributes
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 17 of file DalitzBoxSet.h.

Constructor & Destructor Documentation

DalitzBoxSet() [1/2]

DalitzBoxSet::DalitzBoxSet

(

TRandom *

rnd = gRandom

)

Definition at line 68 of file DalitzBoxSet.cpp.

  : _ready(false)
{
  _rnd = rnd;
}

DalitzBoxSet() [2/2]

DalitzBoxSet::DalitzBoxSet

(

const DalitzBoxSet &

other

)

Definition at line 74 of file DalitzBoxSet.cpp.

  : IEventGenerator<IDalitzEvent>()
  , MINT::PolymorphVector<DalitzBox>(other)
  , _ready(other._ready)
  , _volumeProbs(other._volumeProbs)
  , _amps(other._amps)
  , _rnd(other._rnd)
{
  if(_ready) makeRelations();
}

~DalitzBoxSet()

virtual DalitzBoxSet::~DalitzBoxSet ( )

inlinevirtual

Definition at line 73 of file DalitzBoxSet.h.

{}

Member Function Documentation

add() [1/3]

void DalitzBoxSet::add

(

const DalitzBox &

box

)

Definition at line 218 of file DalitzBoxSet.cpp.

                                          {
  _ready = false;
  cout << " DalitzBoxSet::addBox: adding box "  
       << " with " << box.eventList().size()
       << " events in it."
       << endl;
  this->push_back(box);
  cout << " done that." << endl;
}

add() [2/3]

void DalitzBoxSet::add

(

const MINT::PolymorphVector< DalitzBox > &

boxes

)

Definition at line 228 of file DalitzBoxSet.cpp.

                                                                 {
  _ready = false;
  for(unsigned int i=0; i<boxes.size(); i++){
    this->push_back(boxes[i]);
  }
}

add() [3/3]

void DalitzBoxSet::add

(

const DalitzBoxSet &

boxes

)

Definition at line 235 of file DalitzBoxSet.cpp.

                                               {
  _ready = false;
  for(unsigned int i=0; i<boxes.size(); i++){
    this->push_back(boxes[i]);
  }
}

AreaSum()

double DalitzBoxSet::AreaSum ( ) const

protected

Definition at line 170 of file DalitzBoxSet.cpp.

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

callSortYourselfForDebug()

void DalitzBoxSet::callSortYourselfForDebug

(

)

Definition at line 265 of file DalitzBoxSet.cpp.

                                           {
  sortYourself();
}

ensureFreshEvents()

bool DalitzBoxSet::ensureFreshEvents ( )

virtual

Implements MINT::IEventGenerator< IDalitzEvent >.

Definition at line 391 of file DalitzBoxSet.cpp.

                                    {
  _rnd->SetSeed(time(0)*4);
  setRnd(_rnd);
  return true;
}

exhausted()

virtual bool DalitzBoxSet::exhausted ( ) const

inlinevirtual

Implements MINT::IEventGenerator< IDalitzEvent >.

Reimplemented in DalitzBoxSet_Method2.

Definition at line 70 of file DalitzBoxSet.h.

{return false;}

generateEvent()

DalitzEvent DalitzBoxSet::generateEvent ( )

virtual

Definition at line 369 of file DalitzBoxSet.cpp.

                                       {
  counted_ptr<DalitzEvent> evtPtr = generateEventForOwner();
  if(0==evtPtr){
    cout << "WARNING: DalitzBoxSet::generateEvent()"
     << " failed to generate event despite LARGE"
     << " number of tries." << endl;
    DalitzEvent event;
  }
  DalitzEvent evt(*evtPtr);
  //  delete evtPtr;
  return evt;
}

generateEventForOwner()

counted_ptr< DalitzEvent > DalitzBoxSet::generateEventForOwner ( )

virtual

Definition at line 358 of file DalitzBoxSet.cpp.

                                                            {
  counted_ptr<DalitzEvent> newEvent(0);
  unsigned int counter=0;
  unsigned int veryLargeNumber = 1000000000;
  do{
    newEvent = tryEvent();
    counter++;
  }while(0==newEvent && counter < veryLargeNumber);
  
  return newEvent;
}

getBoxesReady()

void DalitzBoxSet::getBoxesReady ( )

protected

Definition at line 257 of file DalitzBoxSet.cpp.

                                {
  for(unsigned int i=0; i < this->size(); i++){
    cout << " getting box " << i << " ready." << endl;
    (*this)[i].getReady();
    cout << " got it." << endl;
  }
}

getReady()

void DalitzBoxSet::getReady ( )

protectedvirtual

Reimplemented in DalitzBoxSet_Method2.

Definition at line 268 of file DalitzBoxSet.cpp.

                           {
  time_t startTime = time(0);

  cout << "DalitzBoxSet::getReady()" << endl;
  cout << "before removing dupliates: " 
       << this->size() << " boxes." << endl;
  removeDuplicates();
  cout << "after removing dupliates: " 
       << this->size() << " boxes." << endl;
  sortYourself();
  cout << " sorted" << endl;
  makeRelations();
  cout << " made relations" << endl;
  cout << "I have the following " << this->size() << " boxes: " << endl;
  cout << *this << endl;

  setBoxPDFs();
  cout << " set the pdfs " << endl;
  setBoxRnds();
  cout << " set the rnds" << endl;
  getBoxesReady();
  cout << " got Boxes ready" << endl;
  makeVolumeProbIntervals();
  cout << " made the intervalse " << endl;

  double deltaT = difftime(time(0), startTime);
  cout << "DalitzBoxSet: after only:" << deltaT/60 << " min" << endl;
  cout << "DalitzBoxSet: I am ready." << endl;
  _ready = true;
}

makeRelations()

void DalitzBoxSet::makeRelations

(

)

Definition at line 147 of file DalitzBoxSet.cpp.

                                {
  if(this->empty()) return;
  (*this)[0].setDaddy(0);

  if(this->size() <=1) return;
  //  sortYourself();
  for(unsigned int i=1; i< this->size(); i++){
    (*this)[i].setDaddy(&((*this)[i-1]));
  }
}

makeVolumeProbIntervals()

void DalitzBoxSet::makeVolumeProbIntervals ( )

protected

Definition at line 178 of file DalitzBoxSet.cpp.

                                          {
  bool dbThis=false;
  if(this->empty()) return;
  // intervals each with a length proportional
  // to the volume of the corresponding box.
  // All put together add up to one.
  // For random number generation.

  _volumeProbs.clear();
  _volumeProbs.resize(this->size());

  double totalVolume = VolumeSum();
  double sum=0;
  for(unsigned int i=0; i < this->size(); i++){
    sum += (*this)[i].volume()/totalVolume;
    if(dbThis) cout << " volume probs [" << i <<"] = " << sum << endl;
    _volumeProbs[i] = sum;
  }
}

newEvent()

counted_ptr< IDalitzEvent > DalitzBoxSet::newEvent ( )

virtual

Implements MINT::IEventGenerator< IDalitzEvent >.

Reimplemented in DalitzBoxSet_Method2.

Definition at line 345 of file DalitzBoxSet.cpp.

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

pickRandomVolume()

int DalitzBoxSet::pickRandomVolume ( )

protected

Definition at line 198 of file DalitzBoxSet.cpp.

                                  {
  if(! _ready) getReady();

  if(_volumeProbs.size() != this->size()){
    makeVolumeProbIntervals();
  }

  double rndNumber = _rnd->Rndm();
  for(unsigned int i=0; i< _volumeProbs.size(); i++){
    if(_volumeProbs[i] > rndNumber) return i;
  }
  cout << "WARNING in DalitzBoxSet::pickRandomVolume():"
       << " How odd - should never have gotten here."
       << " rndNumber = " << rndNumber
       << ", _volumeProbs[this->size()-1] = "
       << _volumeProbs[this->size()-1]
       << endl;
  return this->size()-1;
}

print()

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

virtual

Definition at line 396 of file DalitzBoxSet.cpp.

                                            {
  os << " DalitzBoxSet with " << this->size() << " boxes:"
     << "\n======================================";
  for(unsigned int i=0; i< this->size(); i++){
    os << "\n box " << i << "\n" << (*this)[i];
    os << "\n--------------------------------------";
  }
}

removeDuplicates()

void DalitzBoxSet::removeDuplicates

(

)

Definition at line 85 of file DalitzBoxSet.cpp.

                                   {

  //  lessBoxBySimilarMappedArea lesser;
  lesserBoxByMappedArea lesserArea;
  sameBoxByMappedArea same;
  moreBoxByGuessedHeight moreHeight;
  //sameBoxBySimilarMappedArea same;

  sort(this->begin(), this->end(), moreHeight);
  stable_sort(this->begin(), this->end(), lesserArea);

  std::vector<DalitzBox>::iterator 
    uniqueEnd = unique(this->begin(), this->end(), same);


  cout << "removeDuplicates: after sorting: " << endl;
  for(unsigned int i=0; i< this->size(); i++){
    bool equalToPrevious = false;
    if(i >= 1) equalToPrevious = same((*this)[i], (*this)[i-1]);
    cout << " " << i << ") "
     << (*this)[i]
     << "\n is equal to previous? "
     << equalToPrevious
     << '\n' << endl;
  }

  //  sameBoxBySimilarMappedArea same;
  /*
  this->erase( (unique(this->begin(), this->end(), same))
           , this->end());
  */
  this->erase(uniqueEnd, this->end());

  cout << "removeDuplicates: after removing: " << endl;
  for(unsigned int i=0; i< this->size(); i++){
    cout << " " << i << ") "
     << (*this)[i]
     << '\n' << endl;
  }
}

setBoxPDFs()

void DalitzBoxSet::setBoxPDFs ( )

protected

Definition at line 247 of file DalitzBoxSet.cpp.

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

setBoxRnds()

void DalitzBoxSet::setBoxRnds ( )

protected

Definition at line 252 of file DalitzBoxSet.cpp.

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

setPDF()

void DalitzBoxSet::setPDF

(

MINT::IReturnRealForEvent< IDalitzEvent > *

amps

)

Definition at line 242 of file DalitzBoxSet.cpp.

                                                                {
  _ready = false;
  _amps = amps;
}

setRnd()

bool DalitzBoxSet::setRnd

(

TRandom *

rnd

)

Definition at line 382 of file DalitzBoxSet.cpp.

                                     {
  _ready = false;
  if(0 == rnd) _rnd = gRandom;
  else _rnd = rnd; // this *should* do it, but for safety:
  for(unsigned int i=0; this->size(); i++){
    (*this)[i].setRnd(_rnd);
  }
  return true;
}

sortYourself()

void DalitzBoxSet::sortYourself

(

)

Definition at line 125 of file DalitzBoxSet.cpp.

                               {
  cout << " DalitzBoxSet::sortYourself() " << endl;
    
  if(this->size() <= 1) return;
  cout << " creating sba" << endl;
  //sortBoxByAreaSize sba;
  moreBoxByGuessedHeight moreHeight;

  /*
  if(is_sorted(this->begin(), this->end(), sba)){
    return;
  }
  */
  stable_sort(this->begin(), this->end(), moreHeight);
  //  stable_sort(this->begin(), this->end(), sba);

  for(unsigned int i=0; i< this->size(); i++){
    (*this)[i].setNumber(i);
  }
  cout << " done sorting " << endl;
}

split()

DalitzBoxSet DalitzBoxSet::split

(

unsigned int

nWays

)

const

Definition at line 300 of file DalitzBoxSet.cpp.

                                                        {
  if(nWays == 1) return (*this);
 
  DalitzBoxSet newSet;
  if(nWays == 0) return newSet;

  cout << " DalitzBoxSet: before splitting I have "
       << this->size() << " boxes"
       << endl;
  for(unsigned int i=0; i< this->size(); i++){
    newSet.add((*this)[i].split(nWays));
  }
  cout << "after, it's " << newSet.size() << endl;
  return newSet;
}

splitIfWiderThan()

DalitzBoxSet DalitzBoxSet::splitIfWiderThan

(

double

maxWidth

)

const

Definition at line 316 of file DalitzBoxSet.cpp.

                                                                {
 
  DalitzBoxSet newSet;

  cout << " DalitzBoxSet: before splitting I have "
       << this->size() << " boxes"
       << endl;
  for(unsigned int i=0; i< this->size(); i++){
    newSet.add((*this)[i].splitIfWiderThan(maxWidth));
  }
  cout << "after, it's " << newSet.size() << endl;
  return newSet;
}

tryEvent()

counted_ptr< DalitzEvent > DalitzBoxSet::tryEvent ( )

virtual

Reimplemented in DalitzBoxSet_Method2.

Definition at line 331 of file DalitzBoxSet.cpp.

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

tryWeightedEvent()

counted_ptr< DalitzEvent > DalitzBoxSet::tryWeightedEvent ( )

virtual

Reimplemented in DalitzBoxSet_Method2.

Definition at line 338 of file DalitzBoxSet.cpp.

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

VolumeSum()

double DalitzBoxSet::VolumeSum ( ) const

protected

Definition at line 158 of file DalitzBoxSet.cpp.

                                    {
  double sum = 0;
  for(unsigned int i=0; i< this->size(); i++){
    cout << " volume number " << i << ") " 
     << (*this)[i].volume()
     << endl;
    sum += (*this)[i].volume();
  }
  cout << " Volume sum: " << sum;
  return sum;
}

Member Data Documentation

_amps

MINT::IReturnRealForEvent<IDalitzEvent>* DalitzBoxSet::_amps

protected

Definition at line 25 of file DalitzBoxSet.h.

_ready

bool DalitzBoxSet::_ready

protected

Definition at line 23 of file DalitzBoxSet.h.

_rnd

TRandom* DalitzBoxSet::_rnd

protected

Definition at line 27 of file DalitzBoxSet.h.

_volumeProbs

MINT::PolymorphVector<double> DalitzBoxSet::_volumeProbs

protected

Definition at line 24 of file DalitzBoxSet.h.

---

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

• Mint/DalitzBoxSet.h
• src/Mojito/BetterMC/DalitzBoxSet.cpp