A small algorithm class using the random number service. More...

#include </builds/gaudi/Gaudi/GaudiTestSuite/src/RandomNumber/RandomNumberAlg.h>

Inheritance diagram for RandomNumberAlg:

[legend]

Collaboration diagram for RandomNumberAlg:

[legend]

Public Member Functions
	RandomNumberAlg (const std::string &name, ISvcLocator *pSvcLocator)
	Constructor: A constructor of this form must be provided.

	~RandomNumberAlg () override
	Standard Destructor.

StatusCode	initialize () override
	Customized initialisation.

StatusCode	finalize () override
	Customized finalisation.

StatusCode	execute () override
	Event callback.

Public Member Functions inherited from Algorithm
	LegacyAlgorithmAdapter (std::string name, ISvcLocator *svcloc, std::string version=PACKAGE_VERSION)

Public Member Functions inherited from Gaudi::details::LegacyAlgorithmAdapter
	LegacyAlgorithmAdapter (std::string name, ISvcLocator *svcloc, std::string version=PACKAGE_VERSION)

const EventContext &	getContext () const

bool	isExecuted () const
	Has this algorithm been executed since the last reset?

void	setExecuted (bool state) const
	Set the executed flag to the specified state.

bool	filterPassed () const
	Did this algorithm pass or fail its filter criterion for the last event?

void	setFilterPassed (bool state) const
	Set the filter passed flag to the specified state.

Public Member Functions inherited from Gaudi::Algorithm
	Algorithm (std::string name, ISvcLocator *svcloc, std::string version=PACKAGE_VERSION)
	Constructor.

StatusCode	sysStart () override
	Reinitialization method invoked by the framework.

StatusCode	sysInitialize () override
	Initialization method invoked by the framework.

StatusCode	sysReinitialize () override
	Reinitialization method invoked by the framework.

StatusCode	sysRestart () override
	Restart method invoked by the framework.

StatusCode	sysExecute (const EventContext &ctx) override
	The actions to be performed by the algorithm on an event.

StatusCode	sysStop () override
	System stop.

StatusCode	sysFinalize () override
	System finalization.

const std::string &	name () const override
	The identifying name of the algorithm object.

const Gaudi::StringKey &	nameKey () const override

const std::string &	type () const override
	The type of the algorithm object.

void	setType (std::string type) override

const std::string &	version () const override

unsigned int	index () const override

StatusCode	configure () override
	Dummy implementation of IStateful::configure() method.

StatusCode	terminate () override
	Dummy implementation of IStateful::terminate() method.

StatusCode	initialize () override
	the default (empty) implementation of IStateful::initialize() method

StatusCode	start () override
	the default (empty) implementation of IStateful::start() method

StatusCode	stop () override
	the default (empty) implementation of IStateful::stop() method

StatusCode	finalize () override
	the default (empty) implementation of IStateful::finalize() method

StatusCode	reinitialize () override
	the default (empty) implementation of IStateful::reinitialize() method

StatusCode	restart () override
	the default (empty) implementation of IStateful::restart() method

Gaudi::StateMachine::State	FSMState () const override
	returns the current state of the algorithm

Gaudi::StateMachine::State	targetFSMState () const override
	returns the state the algorithm will be in after the ongoing transition

bool	isEnabled () const override
	Is this algorithm enabled or disabled?

bool	isSequence () const override
	Are we a Sequence?

unsigned int	errorCount () const
	Get the number of failures of the algorithm.

SmartIF< IService >	service (std::string_view name, const bool createIf=true, const bool quiet=false) const
	Return a pointer to the service identified by name (or "type/name")

template<class T>
SmartIF< T >	service (std::string_view name, bool createIf=true, bool quiet=false) const

SmartIF< IAuditorSvc > &	auditorSvc () const
	The standard auditor service.May not be invoked before sysInitialize() has been invoked.

SmartIF< IChronoStatSvc > &	chronoSvc () const
	The standard Chrono & Stat service, Return a pointer to the service if present.

SmartIF< IDataProviderSvc > &	detSvc () const
	The standard detector data service.

SmartIF< IConversionSvc > &	detCnvSvc () const
	The standard detector data persistency conversion service.

SmartIF< IDataProviderSvc > &	eventSvc () const
	The standard event data service.

SmartIF< IDataProviderSvc > &	evtSvc () const
	shortcut for method eventSvc

SmartIF< IConversionSvc > &	eventCnvSvc () const
	The standard event data persistency conversion service.

SmartIF< IHistogramSvc > &	histoSvc () const
	The standard histogram service.

SmartIF< INTupleSvc > &	ntupleSvc () const
	The standard N tuple service.

SmartIF< IRndmGenSvc > &	randSvc () const
	The standard RandomGen service, Return a pointer to the service if present.

SmartIF< IToolSvc > &	toolSvc () const
	The standard ToolSvc service, Return a pointer to the service if present.

SmartIF< IExceptionSvc > &	exceptionSvc () const
	Get the exception Service.

SmartIF< IAlgContextSvc > &	contextSvc () const
	get Algorithm Context Service

SmartIF< ITimelineSvc > &	timelineSvc () const

SmartIF< ISvcLocator > &	serviceLocator () const override
	The standard service locator.

SmartIF< ISvcLocator > &	svcLoc () const
	shortcut for method serviceLocator

SmartIF< IHiveWhiteBoard > &	whiteboard () const

SmartIF< IAlgExecStateSvc > &	algExecStateSvc () const

bool	registerContext () const
	register for Algorithm Context Service?

template<class T>
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, ToolHandle< T > &hndl, const std::string &doc="none")

template<class T>
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, ToolHandleArray< T > &hndlArr, const std::string &doc="none")

SmartIF< IMonitorSvc > &	monitorSvc () const
	Access the monitor service.

template<class T>
void	declareInfo (const std::string &name, const T &var, const std::string &desc) const
	Declare monitoring information.

void	declareInfo (const std::string &name, const std::string &format, const void *var, int size, const std::string &desc) const
	Declare monitoring information (special case)

void	acceptDHVisitor (IDataHandleVisitor *) const override

void	registerTool (IAlgTool *tool) const

void	deregisterTool (IAlgTool *tool) const

template<class T>
StatusCode	declareTool (ToolHandle< T > &handle, bool createIf=true)

template<class T>
StatusCode	declareTool (ToolHandle< T > &handle, const std::string &toolTypeAndName, bool createIf=true)

template<class T>
void	addToolsArray (ToolHandleArray< T > &hndlArr)

const std::vector< IAlgTool * > &	tools () const

bool	isAsynchronous () const

void	setAsynchronous (bool value)

unsigned int	cardinality () const override
	Return the cardinality.

const std::vector< std::string > &	neededResources () const override

AlgExecStateRef	execState (const EventContext &ctx) const override
	get the AlgExecStateRef of current algorithm Actually a small wrapper around it, thus the plain object return type

std::ostream &	toControlFlowExpression (std::ostream &os) const override
	Produce string represention of the control flow expression.

Public Member Functions inherited from DataHandleHolderBase< PropertyHolder< CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > > > >
std::vector< Gaudi::DataHandle * >	inputHandles () const override

std::vector< Gaudi::DataHandle * >	outputHandles () const override

virtual const DataObjIDColl &	extraInputDeps () const override

virtual const DataObjIDColl &	extraOutputDeps () const override

void	declare (Gaudi::DataHandle &handle) override

void	renounce (Gaudi::DataHandle &handle) override

bool	renounceInput (const DataObjID &id) override

const DataObjIDColl &	inputDataObjs () const override

const DataObjIDColl &	outputDataObjs () const override

void	addDependency (const DataObjID &id, const Gaudi::DataHandle::Mode &mode) override

Public Member Functions inherited from extends< PropertyHolder< CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > > >, IDataHandleHolder >
void const *	i_cast (const InterfaceID &tid) const override
	Implementation of IInterface::i_cast.

StatusCode	queryInterface (const InterfaceID &ti, void **pp) override
	Implementation of IInterface::queryInterface.

std::vector< std::string >	getInterfaceNames () const override
	Implementation of IInterface::getInterfaceNames.

Public Member Functions inherited from PropertyHolder< CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > > >
StatusCode	setProperty (const Gaudi::Details::PropertyBase &p)
	Set the property from a property.

StatusCode	setProperty (const std::string &name, const char *v)
	Special case for string literals.

StatusCode	setProperty (const std::string &name, const std::string &v)
	Special case for std::string.

StatusCode	setProperty (const std::string &name, const TYPE &value)
	set the property form the value

	PropertyHolder ()=default

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Details::PropertyBase &prop)
	Declare a property.

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, TYPE &value, const std::string &doc="none")
	Helper to wrap a regular data member and use it as a regular property.

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, Gaudi::Property< TYPE, VERIFIER, HANDLERS > &prop, const std::string &doc="none")
	Declare a PropertyBase instance setting name and documentation.

Gaudi::Details::PropertyBase *	declareRemoteProperty (const std::string &name, IProperty *rsvc, const std::string &rname="")
	Declare a remote property.

StatusCode	setProperty (const std::string &name, const Gaudi::Details::PropertyBase &p) override
	set the property from another property with a different name

StatusCode	setProperty (const std::string &s) override
	set the property from the formatted string

StatusCode	setProperty (const Gaudi::Details::PropertyBase &p)
	Set the property from a property.

StatusCode	setProperty (const std::string &name, const char *v)
	Special case for string literals.

StatusCode	setProperty (const std::string &name, const std::string &v)
	Special case for std::string.

StatusCode	setProperty (const std::string &name, const TYPE &value)
	set the property form the value

StatusCode	setPropertyRepr (const std::string &n, const std::string &r) override
	set the property from name and value string representation

StatusCode	getProperty (Gaudi::Details::PropertyBase *p) const override
	get the property

const Gaudi::Details::PropertyBase &	getProperty (std::string_view name) const override
	get the property by name

StatusCode	getProperty (std::string_view n, std::string &v) const override
	convert the property to the string

const std::vector< Gaudi::Details::PropertyBase * > &	getProperties () const override
	get all properties

bool	hasProperty (std::string_view name) const override
	Return true if we have a property with the given name.

Gaudi::Details::PropertyBase *	property (std::string_view name) const
	\fixme property and bindPropertiesTo should be protected

void	bindPropertiesTo (Gaudi::Interfaces::IOptionsSvc &optsSvc)

	PropertyHolder (const PropertyHolder &)=delete

PropertyHolder &	operator= (const PropertyHolder &)=delete

Public Member Functions inherited from CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > >
MSG::Level	msgLevel () const
	get the cached level (originally extracted from the embedded MsgStream)

bool	msgLevel (MSG::Level lvl) const
	get the output level from the embedded MsgStream

Public Member Functions inherited from CommonMessagingBase
virtual	~CommonMessagingBase ()=default
	Virtual destructor.

const SmartIF< IMessageSvc > &	msgSvc () const
	The standard message service.

MsgStream &	msgStream () const
	Return an uninitialized MsgStream.

MsgStream &	msgStream (const MSG::Level level) const
	Predefined configurable message stream for the efficient printouts.

MsgStream &	always () const
	shortcut for the method msgStream(MSG::ALWAYS)

MsgStream &	fatal () const
	shortcut for the method msgStream(MSG::FATAL)

MsgStream &	err () const
	shortcut for the method msgStream(MSG::ERROR)

MsgStream &	error () const
	shortcut for the method msgStream(MSG::ERROR)

MsgStream &	warning () const
	shortcut for the method msgStream(MSG::WARNING)

MsgStream &	info () const
	shortcut for the method msgStream(MSG::INFO)

MsgStream &	debug () const
	shortcut for the method msgStream(MSG::DEBUG)

MsgStream &	verbose () const
	shortcut for the method msgStream(MSG::VERBOSE)

MsgStream &	msg () const
	shortcut for the method msgStream(MSG::INFO)

Protected Attributes
Rndm::Numbers	m_numbers
	Allocate wrapper for random number generator.

uint32_t	m_initial { 0 }
	Initial seed to fill deterministic random numbers.

NTuple::Tuple *	m_ntuple { nullptr }
	Pointer to N-tuple.

NTuple::Item< int >	m_int
	N-tuple items.

NTuple::Item< uint32_t >	m_deter

NTuple::Item< float >	m_gauss

NTuple::Item< float >	m_exponential

NTuple::Item< float >	m_poisson

Protected Attributes inherited from Gaudi::Algorithm
std::unique_ptr< IDataHandleVisitor >	m_updateDataHandles
	Hook for for derived classes to provide a custom visitor for data handles.

Protected Attributes inherited from DataHandleHolderBase< PropertyHolder< CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > > > >
DataObjIDColl	m_inputDataObjs

DataObjIDColl	m_outputDataObjs

Additional Inherited Members
Public Types inherited from Gaudi::Algorithm
typedef Gaudi::PluginService::Factory< IAlgorithm (const std::string &, ISvcLocator )>	Factory

Public Types inherited from extends< PropertyHolder< CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > > >, IDataHandleHolder >
using	base_class
	Typedef to this class.

using	extend_interfaces_base
	Typedef to the base of this class.

Public Types inherited from PropertyHolder< CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > > >
using	PropertyHolderImpl
	Typedef used to refer to this class from derived classes, as in.

Public Types inherited from CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > >
using	base_class

Public Types inherited from extend_interfaces< Interfaces... >
using	ext_iids
	take union of the ext_iids of all Interfaces...

Public Attributes inherited from Gaudi::Algorithm
friend	AlgorithmManager

Protected Member Functions inherited from Gaudi::details::LegacyAlgorithmAdapter
bool	isReEntrant () const override

Protected Member Functions inherited from Gaudi::Algorithm
std::vector< IAlgTool * > &	tools ()

bool	isInitialized () const override
	Has the Algorithm already been initialized?

bool	isFinalized () const override
	Has the Algorithm already been finalized?

void	setIndex (const unsigned int &idx) override
	set instantiation index of Alg

bool	isReEntrant () const override

Protected Member Functions inherited from DataHandleHolderBase< PropertyHolder< CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > > > >
void	initDataHandleHolder ()
	initializes all handles - called by the sysInitialize method of any descendant of this

Protected Member Functions inherited from CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > >
MSG::Level	setUpMessaging () const
	Set up local caches.

MSG::Level	resetMessaging ()
	Reinitialize internal states.

void	updateMsgStreamOutputLevel (int level)
	Update the output level of the cached MsgStream.

Detailed Description

A small algorithm class using the random number service.

Author: M.Frank

Date: 1999

Definition at line 24 of file RandomNumberAlg.h.

Constructor & Destructor Documentation

◆ RandomNumberAlg()

RandomNumberAlg::RandomNumberAlg	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Constructor: A constructor of this form must be provided.

Algorithm parameters which can be set at run time must be declared.

This should be done in the constructor.

Definition at line 40 of file RandomNumberAlg.cpp.

41 : Algorithm( name, pSvcLocator ) {}

Gaudi::Algorithm::Algorithm

Algorithm(std::string name, ISvcLocator *svcloc, std::string version=PACKAGE_VERSION)

Constructor.

Definition Algorithm.h:98

Gaudi::Algorithm::name

const std::string & name() const override

The identifying name of the algorithm object.

Definition Algorithm.cpp:529

◆ ~RandomNumberAlg()

RandomNumberAlg::~RandomNumberAlg ( )

override

Standard Destructor.

Definition at line 44 of file RandomNumberAlg.cpp.

                                  {
  // do not print messages if we are created in genconf
  const std::string cmd = System::cmdLineArgs()[0];
  if ( cmd.find( "genconf" ) != std::string::npos ) return;
 
  std::cout << "Destructor Called for " << name() << std::endl;
}

Member Function Documentation

◆ execute()

StatusCode RandomNumberAlg::execute ( )

overridevirtual

Event callback.

Implements Algorithm.

Definition at line 125 of file RandomNumberAlg.cpp.

                                    {
  StatusCode status;
  static int count = 0;
 
  Rndm::Numbers gauss( randSvc(), Rndm::Gauss( 0.5, 0.2 ) );
  Rndm::Numbers exponential( randSvc(), Rndm::Exponential( 0.2 ) );
  Rndm::Numbers poisson( randSvc(), Rndm::Poisson( 0.3 ) );
 
  // Return integer in interval [0, size) from random integer in interval [0, MAX_INT]
  auto scale = []( uint32_t x, uint32_t size ) {
    const uint32_t denom = boost::integer_traits<uint32_t>::const_max / size;
    return x / denom;
  };
 
  m_int         = ++count;
  m_deter       = scale( QuasiRandom::mix32( m_initial, m_int ), 100 );
  m_gauss       = (float)gauss();
  m_exponential = (float)exponential();
  m_poisson     = (float)poisson();
 
  status = m_ntuple->write();
  if ( !status.isSuccess() ) { error() << "Cannot fill NTuple" << endmsg; }
  return StatusCode::SUCCESS;
}

◆ finalize()

StatusCode RandomNumberAlg::finalize ( )

override

Customized finalisation.

Definition at line 150 of file RandomNumberAlg.cpp.

                                     {
  m_numbers.finalize().ignore();
  return StatusCode::SUCCESS;
}

◆ initialize()

StatusCode RandomNumberAlg::initialize ( )

override

Customized initialisation.

Definition at line 54 of file RandomNumberAlg.cpp.

                                       {
  StatusCode status = Algorithm::initialize();
  if ( !status ) return status;
 
  //
  // The first example is for purists:
  // Every step is done by hand....tends to become complicated,
  // but shows the usage of the raw interfaces
  //
  // Get random number generator:
  auto gen = randSvc()->generator( Rndm::Gauss( 0.5, 0.2 ) );
  if ( gen ) {
    std::vector<double> numbers;
    gen->shootArray( numbers, 5000 ).ignore();
    IHistogram1D* his = histoSvc()->book( "1", "Gauss", 40, 0., 3. );
    for ( unsigned int i = 0; i < numbers.size(); i++ ) his->fill( numbers[i], 1.0 );
 
    for ( int j = 0; j < 5000; j++ ) his->fill( gen->shoot(), 1.0 );
  }
 
  //
  // Now we do it as proposed for LOCAL usage of the wrapper
  //
  Rndm::Numbers exponential( randSvc(), Rndm::Exponential( 0.2 ) );
  if ( exponential ) {
    IHistogram1D* his = histoSvc()->book( "2", "Exponential", 40, 0., 3. );
    for ( long j = 0; j < 5000; j++ ) his->fill( exponential(), 1.0 );
  } else {
    return StatusCode::FAILURE;
  }
 
  //
  // Now we do it as proposed for GLOBAL usage of the wrapper
  // - Initialize the wrapper allocated in the header file
  //
  status = m_numbers.initialize( randSvc(), Rndm::Poisson( 0.3 ) );
  if ( !status.isSuccess() ) { return status; }
 
  // The GLOBAL wrapper is now initialized and ready for use.
  // The code below could go anywhere. It is only for simplicity
  // in the "initialize" method!
  {
    IHistogram1D* hispoisson = histoSvc()->book( "3", "Poisson", 40, 0., 3. );
    for ( long j = 0; j < 5000; j++ ) hispoisson->fill( m_numbers(), 1.0 );
  }
  //
  // Test Gaussian Tail distribution
  //
  Rndm::Numbers gaussiantail( randSvc(), Rndm::GaussianTail( 20., 10. ) );
  if ( gaussiantail ) {
    IHistogram1D* his = histoSvc()->book( "4", "GaussianTail", 50, 0., 50. );
    for ( long j = 0; j < 50009; j++ ) his->fill( gaussiantail(), 1.0 );
  } else {
    return StatusCode::FAILURE;
  }
 
  // Initial randomness for deterministic random numbers
  m_initial = QuasiRandom::mixString( name().size(), name() );
 
  // Book N-tuple
  m_ntuple = ntupleSvc()->book( "/NTUPLES/FILE1/100", CLID_RowWiseTuple, "Hello World" );
  if ( m_ntuple ) {
    status = m_ntuple->addItem( "Event#", m_int );
    status = m_ntuple->addItem( "DeterInt", m_deter );
    status = m_ntuple->addItem( "Gauss", m_gauss );
    status = m_ntuple->addItem( "Exp", m_exponential );
    status = m_ntuple->addItem( "Poisson", m_poisson );
  }
  return status;
}

Member Data Documentation

◆ m_deter

NTuple::Item<uint32_t> RandomNumberAlg::m_deter

protected

Definition at line 34 of file RandomNumberAlg.h.

◆ m_exponential

NTuple::Item<float> RandomNumberAlg::m_exponential

protected

Definition at line 36 of file RandomNumberAlg.h.

◆ m_gauss

NTuple::Item<float> RandomNumberAlg::m_gauss

protected

Definition at line 35 of file RandomNumberAlg.h.

◆ m_initial

uint32_t RandomNumberAlg::m_initial { 0 }

protected

Initial seed to fill deterministic random numbers.

Definition at line 29 of file RandomNumberAlg.h.

29{ 0 };

◆ m_int

NTuple::Item<int> RandomNumberAlg::m_int

protected

N-tuple items.

Definition at line 33 of file RandomNumberAlg.h.

◆ m_ntuple

NTuple::Tuple* RandomNumberAlg::m_ntuple { nullptr }

protected

Pointer to N-tuple.

Definition at line 31 of file RandomNumberAlg.h.

31{ nullptr };

◆ m_numbers

Rndm::Numbers RandomNumberAlg::m_numbers

protected

Allocate wrapper for random number generator.

Definition at line 27 of file RandomNumberAlg.h.

◆ m_poisson

NTuple::Item<float> RandomNumberAlg::m_poisson

protected

Definition at line 37 of file RandomNumberAlg.h.

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

GaudiTestSuite/src/RandomNumber/RandomNumberAlg.h
GaudiTestSuite/src/RandomNumber/RandomNumberAlg.cpp

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ RandomNumberAlg()

◆ ~RandomNumberAlg()

Member Function Documentation

◆ execute()

◆ finalize()

◆ initialize()

Member Data Documentation

◆ m_deter

◆ m_exponential

◆ m_gauss

◆ m_initial

◆ m_int

◆ m_ntuple

◆ m_numbers

◆ m_poisson