SignalGenerator Class Reference

#include <SignalGenerator.h>

Inheritance diagram for SignalGenerator:

Public Member Functions
	SignalGenerator (const DalitzEventPattern &pat, TRandom *rnd=gRandom)
	SignalGenerator (const DalitzEventPattern &pat, IFastAmplitudeIntegrable *amps, TRandom *rnd=gRandom)
	SignalGenerator (const DalitzEventPattern &pat, double rB, double phase, TRandom *rnd=gRandom)
virtual	~SignalGenerator ()
virtual MINT::counted_ptr< IDalitzEvent >	tryDalitzEvent ()
virtual MINT::counted_ptr< IDalitzEvent >	newDalitzEvent ()
virtual MINT::counted_ptr< IDalitzEvent >	newEvent ()
virtual bool	exhausted () const
virtual bool	ensureFreshEvents ()
bool	am_I_generating_what_I_think_I_am_generating (int Nevents=100000)
bool	compareGenerationMethodsForFullPDF (int Nevents=100000)
Public Member Functions inherited from BaseGenerator
	BaseGenerator (TRandom *rnd=gRandom)
	BaseGenerator (const BaseGenerator &other)
void	setMothers3Momentum (MINT::counted_ptr< MINT::IReturnReal > p, MINT::counted_ptr< MINT::IReturnReal > cosTheta=(MINT::counted_ptr< MINT::IReturnReal >) 0, MINT::counted_ptr< MINT::IReturnReal > phi=(MINT::counted_ptr< MINT::IReturnReal >) 0)
void	unsetMothers3Momentum ()
bool	mothers3MomentumIsSet () const
void	setSaveEvents (const std::string &fname="GeneratorEvents.root", const std::string &opt="RECREATE")
void	dontSaveEvents ()
void	setUnweighted (bool uw=true)
void	setWeighted (bool w=true)
void	FillEventList (DalitzEventList &evtList, int NEvents)
void	FillEventList (DiskResidentEventList &evtList, int NEvents)
void	noPrintout ()
virtual	~BaseGenerator ()
Public Member Functions inherited from MINT::IEventGenerator< IDalitzEvent >
virtual	~IEventGenerator ()

Protected Member Functions
bool	makeBoxes ()
MINT::MinuitParameterSet *	myMPS ()
Protected Member Functions inherited from BaseGenerator
bool	saveEvents ()
TVector3	mothers3Momentum () const

Protected Attributes
DalitzEventPattern	_pat
MINT::MinuitParameterSet	_myOwnPSet
MINT::counted_ptr< IFastAmplitudeIntegrable >	_counted_amps
IFastAmplitudeIntegrable *	_amps
MINT::counted_ptr< MINT::IUnweightedEventGenerator< IDalitzEvent > >	_boxes
Protected Attributes inherited from BaseGenerator
std::string	_fname
bool	_unWeighted
MINT::counted_ptr< DiskResidentEventList >	_evtList
bool	_noPrintout
TRandom *	_rnd
MINT::counted_ptr< MINT::IReturnReal >	_mums_p
MINT::counted_ptr< MINT::IReturnReal >	_mums_cosTheta
MINT::counted_ptr< MINT::IReturnReal >	_mums_phi

Detailed Description

Definition at line 25 of file SignalGenerator.h.

Constructor & Destructor Documentation

SignalGenerator() [1/3]

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

Definition at line 11 of file SignalGenerator.cpp.

  : BaseGenerator(rnd)
  , _pat(pat)
  , _myOwnPSet()
  , _counted_amps(new FitAmpSum(pat, myMPS()))
  , _amps(_counted_amps.get())
  , _boxes(_amps->makeEventGenerator(pat))
{
  _boxes->setRnd(rnd);
}

SignalGenerator() [2/3]

SignalGenerator::SignalGenerator	(	const DalitzEventPattern &	pat,
		IFastAmplitudeIntegrable *	amps,
		TRandom *	rnd = gRandom
	)

Definition at line 55 of file SignalGenerator.cpp.

  : BaseGenerator(rnd)
  , _pat(pat)
  , _myOwnPSet()
  , _counted_amps(0)
  , _amps(amps)
  , _boxes(_amps->makeEventGenerator(pat))
{
  _boxes->setRnd(_rnd);
}

SignalGenerator() [3/3]

SignalGenerator::SignalGenerator	(	const DalitzEventPattern &	pat,
		double	rB,
		double	phase,
		TRandom *	rnd = gRandom
	)

Definition at line 23 of file SignalGenerator.cpp.

  : BaseGenerator(rnd)
  , _pat(pat)
  , _myOwnPSet()
  , _counted_amps(0)
  , _amps(0)
  , _boxes(0)
{
  DalitzEventPattern cpPat(pat);
  cpPat[0].antiThis();

  counted_ptr<FitAmpSum> fs(new FitAmpSum(pat, myMPS()));
  FitAmpSum cpAmps(pat, myMPS());
  cpAmps *= polar(rB, phase);
  fs->add(cpAmps);
  _counted_amps = (counted_ptr<IFastAmplitudeIntegrable>) fs;
  _amps = _counted_amps.get();

  makeBoxes();

  /*
  counted_ptr< IUnweightedEventGenerator<IDalitzEvent> > 
    bpt(fs->makeEventGenerator());

  _boxes = bpt;
  _boxes->setRnd(_rnd);
  */
}

~SignalGenerator()

virtual SignalGenerator::~SignalGenerator ( )

inlinevirtual

Definition at line 56 of file SignalGenerator.h.

{}

Member Function Documentation

am_I_generating_what_I_think_I_am_generating()

bool SignalGenerator::am_I_generating_what_I_think_I_am_generating

(

int

Nevents = 100000

)

Definition at line 134 of file SignalGenerator.cpp.

                                                                             {
  if(0 == _boxes) return false;
  return _boxes->am_I_generating_what_I_think_I_am_generating(Nevents);

}

compareGenerationMethodsForFullPDF()

bool SignalGenerator::compareGenerationMethodsForFullPDF

(

int

Nevents = 100000

)

Definition at line 139 of file SignalGenerator.cpp.

                                                                   {
  if(0 == _boxes) return false;
  return _boxes->compareGenerationMethodsForFullPDF(Nevents);
}

ensureFreshEvents()

bool SignalGenerator::ensureFreshEvents ( )

virtual

Reimplemented from BaseGenerator.

Definition at line 126 of file SignalGenerator.cpp.

                                       {
  _boxes->ensureFreshEvents();
  return BaseGenerator::ensureFreshEvents();
}

exhausted()

virtual bool SignalGenerator::exhausted ( ) const

inlinevirtual

Reimplemented from BaseGenerator.

Definition at line 67 of file SignalGenerator.h.

{return false;}

makeBoxes()

bool SignalGenerator::makeBoxes ( )

protected

Definition at line 77 of file SignalGenerator.cpp.

                               {
  if(0 == _amps) return 0;
  counted_ptr< IUnweightedEventGenerator<IDalitzEvent> > 
    bpt(_amps->makeEventGenerator(_pat));
  _boxes = bpt;
  if(0 == _boxes) return 0;
  _boxes->setRnd(_rnd);
  return true;
}

myMPS()

MINT::MinuitParameterSet * SignalGenerator::myMPS ( )

protected

Definition at line 68 of file SignalGenerator.cpp.

                                              {
  // seems like an unnecessary complication
  // but it's a "hook" to add the option to pass
  // other MinuitParameterSets to SignalGenerator
  // should this become necessary.
  return & _myOwnPSet;
}

newDalitzEvent()

counted_ptr< IDalitzEvent > SignalGenerator::newDalitzEvent ( )

virtual

Definition at line 109 of file SignalGenerator.cpp.

                                                         {
//   bool dbThis=false;
  counted_ptr<IDalitzEvent> evtPtr(0);
  int counter(0);
  int largeNumber(1000000);

  do{
    evtPtr = tryDalitzEvent();
  }while(0 == evtPtr &&  counter++ < largeNumber);
  if(saveEvents()) _evtList->Add(evtPtr);
/*  if(dbThis){
    cout << "SignalGenerator::newDalitzEvent:"
     << " just generated this event:\n" << *evtPtr << endl;
  }*/
  return evtPtr;
}

newEvent()

counted_ptr< IDalitzEvent > SignalGenerator::newEvent ( )

virtual

Implements BaseGenerator.

Definition at line 130 of file SignalGenerator.cpp.

                                                   {
  return counted_ptr<IDalitzEvent>(newDalitzEvent());
}

tryDalitzEvent()

counted_ptr< IDalitzEvent > SignalGenerator::tryDalitzEvent ( )

virtual

Definition at line 87 of file SignalGenerator.cpp.

                                                         {
//   bool dbThis=false;
  if(0 == _boxes) makeBoxes();
  if(_unWeighted){
    counted_ptr<IDalitzEvent> evtPtr(_boxes->newUnweightedEvent());
    if(0 != evtPtr) evtPtr->setMothers3Momentum(mothers3Momentum());
/*    if(dbThis && 0 != evtPtr){
      cout << "SignalGenerator::tryDalitzEvent(): made un-weighted event "
       << "with weight " << evtPtr->getWeight() << endl;
    }*/
    return evtPtr;
  }else{
    counted_ptr<IDalitzEvent> evtPtr(_boxes->newEvent());
    if(0 != evtPtr) evtPtr->setMothers3Momentum(mothers3Momentum());
//     if(dbThis && 0 != evtPtr){
//       cout << "SignalGenerator::tryDalitzEvent(): made weighted event "
//     << "with weight " << evtPtr->getWeight() << endl;
//     }
    return evtPtr;
  }
}

Member Data Documentation

_amps

IFastAmplitudeIntegrable* SignalGenerator::_amps

protected

Definition at line 35 of file SignalGenerator.h.

_boxes

MINT::counted_ptr<MINT::IUnweightedEventGenerator<IDalitzEvent> > SignalGenerator::_boxes

protected

Definition at line 36 of file SignalGenerator.h.

_counted_amps

MINT::counted_ptr<IFastAmplitudeIntegrable> SignalGenerator::_counted_amps

protected

Definition at line 34 of file SignalGenerator.h.

_myOwnPSet

MINT::MinuitParameterSet SignalGenerator::_myOwnPSet

protected

Definition at line 30 of file SignalGenerator.h.

_pat

DalitzEventPattern SignalGenerator::_pat

protected

Definition at line 29 of file SignalGenerator.h.

---

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

• Mint/SignalGenerator.h
• src/Mojito/DalitzEventGeneration/SignalGenerator.cpp