BesRndmGenSvc Class Reference

#include <BesRndmGenSvc.h>

Inheritance diagram for BesRndmGenSvc:

Public Types
typedef std::map< std::string, CLHEP::RanluxEngine * >	engineMap
	CLHEP engines typedefs:
typedef engineMap::iterator	engineIter
typedef engineMap::const_iterator	engineConstIter
typedef engineMap::value_type	engineValType

Public Member Functions
CLHEP::HepRandomEngine *	GetEngine (const std::string &StreamName)
	Interface to the CLHEP engine.
void	CreateStream (long seed, const std::string &StreamName)
engineConstIter	begin (void) const
engineConstIter	end (void) const
unsigned int	number_of_streams (void) const
void	setGenseed (long i)
long	getGenseed ()
void	print (const std::string &StreamName)
	Print methods.
void	print (void)
CLHEP::HepRandomEngine *	setOnDefinedSeeds (int EventNumber, const std::string &StreamName)
	set the seeds for an engine. First param will usually be the event number
StatusCode	initialize ()
	Gaudi Service Implementation.
StatusCode	finalize ()
void	handle (const Incident &)
	IIncidentListener implementation. Handles EndEvent incident.
	BesRndmGenSvc (const std::string &name, ISvcLocator *svc)
	Standard Constructor.
virtual	~BesRndmGenSvc ()
	Standard Destructor.

Friends
class	SvcFactory< BesRndmGenSvc >

Detailed Description

Definition at line 33 of file BesRndmGenSvc.h.

Member Typedef Documentation

engineConstIter

typedef engineMap::const_iterator BesRndmGenSvc::engineConstIter

Definition at line 48 of file BesRndmGenSvc.h.

engineIter

typedef engineMap::iterator BesRndmGenSvc::engineIter

Definition at line 47 of file BesRndmGenSvc.h.

engineMap

typedef std::map<std::string, CLHEP::RanluxEngine*> BesRndmGenSvc::engineMap

CLHEP engines typedefs:

Definition at line 46 of file BesRndmGenSvc.h.

engineValType

typedef engineMap::value_type BesRndmGenSvc::engineValType

Definition at line 49 of file BesRndmGenSvc.h.

Constructor & Destructor Documentation

BesRndmGenSvc()

BesRndmGenSvc::BesRndmGenSvc	(	const std::string &	name,
		ISvcLocator *	svc )

Standard Constructor.

Definition at line 43 of file BesRndmGenSvc.cxx.

    : base_class( name, svc ) {
  // Property Default values
  m_read_from_file = false;
  m_file_to_read   = "BesRndmGenSvc.out";
  ofile            = "BesRndmGenSvc.out";
 
  // Set Default values
  m_default_seed         = 11000;
  m_EVTGEN_default_seed  = 35910;
  m_PYTHIA_default_seed  = 93531;
  m_HERWIG_default_seed  = 35391;
  m_LUNDCRM_default_seed = 12456;
  m_SIM_default_seed     = 23569;
  m_MIX_default_seed     = 76543;
  m_KKMC_default_seed    = 26877;
  // Get user's input
  declareProperty( "Seeds", m_streams_seeds );
  declareProperty( "StreamName", m_StreamName = "EVTGEN" );
  declareProperty( "ReadFromFile", m_read_from_file );
  declareProperty( "FileToRead", m_file_to_read );
  declareProperty( "RdmOutFile", ofile );
  declareProperty( "RndmSeed", m_default_seed );
 
  Stream_EVTGEN  = "EVTGEN";
  Stream_PYTHIA  = "PYTHIA";
  Stream_HERWIG  = "HERWIG";
  Stream_LUNDCRM = "LUNDCRM";
  Stream_SIM     = "SIM";
  Stream_MIX     = "MIX";
  Stream_KKMC    = "KKMC";
 
  m_engines = new engineMap();
}

Referenced by SvcFactory< BesRndmGenSvc >.

~BesRndmGenSvc()

BesRndmGenSvc::~BesRndmGenSvc ( )

virtual

Standard Destructor.

Definition at line 79 of file BesRndmGenSvc.cxx.

{ delete m_engines; }

Member Function Documentation

begin()

BesRndmGenSvc::engineConstIter BesRndmGenSvc::begin ( void ) const

inline

Definition at line 122 of file BesRndmGenSvc.h.

                                                                     {
  return m_engines->begin();
}

Referenced by handle(), and initialize().

CreateStream()

void BesRndmGenSvc::CreateStream	(	long	seed,
		const std::string &	StreamName )

Definition at line 259 of file BesRndmGenSvc.cxx.

                                                                         {
  engineConstIter citer = m_engines->find( StreamName );
  if ( citer == m_engines->end() )
    m_engines->insert( engineValType( StreamName, new CLHEP::RanluxEngine() ) );
  engineIter iter = m_engines->find( StreamName );
  ( ( *iter ).second )->setSeed( seed, 3 );
}

Referenced by initialize().

end()

BesRndmGenSvc::engineConstIter BesRndmGenSvc::end ( void ) const

inline

Definition at line 126 of file BesRndmGenSvc.h.

                                                                   {
  return m_engines->end();
}

Referenced by handle(), and initialize().

finalize()

StatusCode BesRndmGenSvc::finalize

(

)

Definition at line 241 of file BesRndmGenSvc.cxx.

                                   {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << " FINALISING " << endmsg;
  return Service::finalize();
}

GetEngine()

CLHEP::HepRandomEngine * BesRndmGenSvc::GetEngine

(

const std::string &

StreamName

)

Interface to the CLHEP engine.

Definition at line 247 of file BesRndmGenSvc.cxx.

                                                                            {
  engineConstIter citer = m_engines->find( StreamName );
  if ( citer == m_engines->end() )
  {
    m_engines->insert( engineValType( StreamName, new CLHEP::RanluxEngine() ) );
    SetStreamSeeds( StreamName );
  }
 
  engineIter iter = m_engines->find( StreamName );
  return (CLHEP::HepRandomEngine*)( *iter ).second;
}

Referenced by handle(), and initialize().

getGenseed()

long BesRndmGenSvc::getGenseed ( )

inline

Definition at line 57 of file BesRndmGenSvc.h.

{ return genSeed; }

handle()

void BesRndmGenSvc::handle

(

const Incident &

inc

)

IIncidentListener implementation. Handles EndEvent incident.

Definition at line 200 of file BesRndmGenSvc.cxx.

                                                {
  MsgStream log( msgSvc(), name() );
  log << MSG::DEBUG << " Handle EndEvent " << endmsg;
 
  if ( inc.type() == "EndEvent" )
  {
 
    m_engines_copy.clear();
    for ( engineConstIter i = begin(); i != end(); ++i )
    {
      CLHEP::HepRandomEngine*    engine = GetEngine( ( *i ).first );
      std::vector<unsigned long> s      = engine->put();
      std::vector<long int>      tseeds;
      for ( unsigned int j = 0; j < s.size(); j++ ) { tseeds.push_back( s[j] ); }
      m_engines_copy.insert(
          std::map<std::string, std::vector<long int>>::value_type( ( *i ).first, tseeds ) );
    }
    // Write the status of the Service into a file
    /*
        ofstream    outfile( ofile.c_str() );
        if ( !outfile )
        {
          log <<MSG::ERROR << "error: unable to open: " << ofile << endmsg;
        }
        else
        {
            for (std::map<std::string, std::vector<long int> >::const_iterator i =
       m_engines_copy.begin(); i != m_engines_copy.end();
                 ++i)
            {
                outfile << (*i).first;
                for(unsigned int j=0; j<((*i).second).size(); j++)
                  outfile << " " << (*i).second[j];
                outfile <<"\n";
            }
        }
    */
    //    print();
  }
}

initialize()

StatusCode BesRndmGenSvc::initialize

(

)

Gaudi Service Implementation.

Incident Service

Definition at line 100 of file BesRndmGenSvc.cxx.

                                     {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "Initializing " << name() << " - package version "
      << "FIXME: PACKAGE_VERSION"
      << "\n INITIALISING RANDOM NUMBER STREAMS. " << endmsg;
 
  if ( !( Service::initialize().isSuccess() ) ) { return StatusCode::FAILURE; }
 
  // set up different stream seeds by default seed
  m_EVTGEN_default_seed  = m_default_seed + 35910;
  m_PYTHIA_default_seed  = m_default_seed + 93531;
  m_HERWIG_default_seed  = m_default_seed + 35391;
  m_LUNDCRM_default_seed = m_default_seed + 12456;
  m_SIM_default_seed     = m_default_seed + 23569;
  m_MIX_default_seed     = m_default_seed + 76543;
  m_KKMC_default_seed    = m_default_seed + 26877;
 
  /// Incident Service
  IIncidentSvc* pIncSvc( 0 );
 
  // set up the incident service:
  if ( !( service( "IncidentSvc", pIncSvc, true ) ).isSuccess() )
  {
    log << MSG::ERROR << "Could not locate IncidentSvc " << endmsg;
    return StatusCode::FAILURE;
  }
 
  // start listening to "EndEvent"
  static const int PRIORITY = 100;
  pIncSvc->addListener( this, "EndEvent", PRIORITY );
 
  if ( m_read_from_file )
  {
    // Read from a file
    ifstream infile( m_file_to_read.c_str() );
    if ( !infile )
    {
      log << MSG::ERROR << " Unable to open: " << m_file_to_read << endmsg;
      return StatusCode::FAILURE;
    }
    else
    {
      std::string buffer;
      while ( std::getline( infile, buffer ) )
      {
        string                     stream;
        std::vector<unsigned long> v;
        // split the space-separated string in 3 words:
        if ( interpretSeeds( buffer, stream, v ) )
        {
          log << MSG::DEBUG << " INITIALISING " << stream << " stream, "
              << " read from file " << m_file_to_read << endmsg;
          CLHEP::HepRandomEngine* engine = GetEngine( stream );
          engine->get( v );
        }
        else
        {
          log << MSG::ERROR << "bad line\n"
              << buffer << "\n in input file " << m_file_to_read << endmsg;
          return StatusCode::FAILURE;
        }
      }
    }
  }
 
  // Create the various streams according to user's request
  for ( VStrings::const_iterator i = m_streams_seeds.begin(); i != m_streams_seeds.end(); ++i )
  {
    string                     stream;
    std::vector<unsigned long> v;
    // split the space-separated string in 3 words:
    if ( interpretSeeds( *i, stream, v ) )
    { log << MSG::VERBOSE << "Seeds property: stream " << stream << endmsg; }
    else
    {
      log << MSG::ERROR << "bad Seeds property\n" << *i << endmsg;
      return StatusCode::FAILURE;
    }
 
    // Check if stream already generated (e.g. from reading a file)
    bool not_found = true;
    if ( number_of_streams() != 0 )
    {
      BesRndmGenSvc::engineConstIter sf = begin();
      do {
        if ( ( *sf ).first == stream ) not_found = false;
        ++sf;
      } while ( sf != end() && not_found );
    }
 
    if ( not_found )
    {
      log << MSG::DEBUG << " INITIALISING " << stream << endmsg;
      CreateStream( v[0], stream );
    }
  }
 
  return StatusCode::SUCCESS;
}

number_of_streams()

unsigned int BesRndmGenSvc::number_of_streams ( void ) const

inline

Definition at line 130 of file BesRndmGenSvc.h.

                                                                 {
  return m_engines->size();
}

Referenced by initialize().

print() [1/2]

void BesRndmGenSvc::print

(

const std::string &

StreamName

)

Print methods.

Definition at line 288 of file BesRndmGenSvc.cxx.

                                                       {
  MsgStream       log( msgSvc(), name() );
  engineConstIter citer = m_engines->find( StreamName );
  if ( citer == m_engines->end() )
  { log << MSG::WARNING << " Stream =  " << StreamName << " NOT FOUND" << endmsg; }
  else
  {
    const long s = ( ( *citer ).second )->getSeed();
    log << MSG::INFO << " Stream =  " << StreamName << ", Seed =  " << s << endmsg;
  }
}

print() [2/2]

void BesRndmGenSvc::print

(

void

)

Definition at line 300 of file BesRndmGenSvc.cxx.

{ print( m_StreamName ); }

Referenced by print().

setGenseed()

void BesRndmGenSvc::setGenseed ( long i )

inline

Definition at line 56 of file BesRndmGenSvc.h.

{ genSeed = i; }

setOnDefinedSeeds()

CLHEP::HepRandomEngine * BesRndmGenSvc::setOnDefinedSeeds	(	int	EventNumber,
		const std::string &	StreamName )

set the seeds for an engine. First param will usually be the event number

Definition at line 302 of file BesRndmGenSvc.cxx.

                                                                                        {
  engineConstIter citer = m_engines->find( StreamName );
  if ( citer == m_engines->end() )
    m_engines->insert( engineValType( StreamName, new CLHEP::RanluxEngine() ) );
  engineIter iter = m_engines->find( StreamName );
  // (SG::simpleStringHash(StreamName));
  long theseed = 10000 * EventNumber;
  ( ( *iter ).second )->setSeed( theseed, 3 );
  return (CLHEP::HepRandomEngine*)( *iter ).second;
}

Friends And Related Symbol Documentation

SvcFactory< BesRndmGenSvc >

friend class SvcFactory< BesRndmGenSvc >

friend

Definition at line 110 of file BesRndmGenSvc.h.

---

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

• BesRndmGenSvc.h
• BesRndmGenSvc.cxx