McContext Class Reference

Inheritance diagram for McContext:

Public Member Functions
	McContext (const McEventSelector *pSelector)
	Standard constructor.
	McContext (const McEventSelector *pSelector, const int &initRunNo, const unsigned int &initEventNo, const unsigned int &evPR)
virtual	~McContext ()
	Standard destructor.
virtual void *	identifier () const
void	setCriteria (const std::string &crit)
void	rewind ()
void	next ()
void	previous ()
int	runNumber () const
unsigned int	eventNumber () const
void	setRunNumber (int runNo)
void	setEventNumber (unsigned int eventNo)

Detailed Description

Definition at line 39 of file McEventSelector.cxx.

Constructor & Destructor Documentation

McContext() [1/2]

McContext::McContext

(

const McEventSelector *

pSelector

)

Standard constructor.

Definition at line 111 of file McEventSelector.cxx.

    : m_pSelector( pSelector ), m_firstEvent( true ), m_lumiBlockNo( 0 ), m_initEventNo( 0 ) {
  // Obtain the IProperty of the ApplicationMgr
  // SmartIF<IProperty> prpMgr(IProperty::interfaceID(), Gaudi::svcLocator());
  SmartIF<IProperty> prpMgr( Gaudi::svcLocator() );
  if ( !prpMgr.isValid() )
  { std::cout << "IProperty interface not found in ApplicationMgr." << std::endl; }
  else { m_appMgrProperty = prpMgr; }
  // get IRealizationSvc
  ISvcLocator* svcLocator = Gaudi::svcLocator();
  StatusCode   status     = svcLocator->service( "RealizationSvc", m_RealizationSvc );
  if ( !status.isSuccess() )
  { std::cout << " Could not initialize Realization Service" << std::endl; }
 
  // get jobSvc
  IDataInfoSvc* jobInfoSvc;
  status = svcLocator->service( "DataInfoSvc", jobInfoSvc );
  if ( status.isSuccess() ) { std::cout << "get the DataInfoSvc" << std::endl; }
  else { std::cout << "could not get the DataInfoSvc." << std::endl; }
 
  // set initial event ID
  m_initEventNo = m_RealizationSvc->getInitEvtID();
 
  std::string sMax;
  status     = m_appMgrProperty->getProperty( "EvtMax", sMax );
  int EvtMax = std::atoi( sMax.c_str() );
  // for realization
  if ( m_RealizationSvc->UseDBFlag() == true )
  {
    std::vector<int> totEvtNo;
    totEvtNo.clear();
    std::vector<float> lumi;
    std::vector<int>   tmp_runList = m_RealizationSvc->getRunList();
    float              totLumi     = 0;
    for ( unsigned int i = 0; i < tmp_runList.size(); i++ )
    {
      float lumi_value = m_RealizationSvc->getLuminosity( tmp_runList[i] );
      lumi.push_back( lumi_value );
      totLumi += lumi_value;
    }
 
    m_runList.clear();
 
    int totSimEvt = 0;
    int evtSubNo  = 0;
    if ( totLumi != 0 )
    {
      for ( unsigned int i = 0; i < lumi.size(); i++ )
      {
        // if(i == 0) m_evtNoList.push_back(0);
        double ratio = lumi[i] / totLumi * EvtMax;
        evtSubNo     = int( ratio );
        if ( ( ratio - evtSubNo ) >= 0.5 ) evtSubNo = evtSubNo + 1;
        totSimEvt += evtSubNo;
        if ( evtSubNo == 0 )
        {
          std::cout << "The run " << tmp_runList[i]
                    << " is not simulated, due to the luminosity is too small!" << std::endl;
        }
        else
        {
          // m_evtNoList.push_back(totSimEvt);
          m_evtNoList.push_back( evtSubNo );
          m_runList.push_back( tmp_runList[i] );
          totEvtNo.push_back( tmp_runList[i] );
          totEvtNo.push_back( evtSubNo );
        }
        std::cout << "Total " << evtSubNo << " events need to be simulated in run "
                  << tmp_runList[i] << std::endl;
      }
 
      if ( ( EvtMax - totSimEvt ) != 0 )
      {
        unsigned int effRunSize = m_evtNoList.size();
        if ( effRunSize > 0 )
        {
          m_evtNoList[effRunSize - 1] = m_evtNoList[effRunSize - 1] + ( EvtMax - totSimEvt );
          effRunSize                  = totEvtNo.size();
          totEvtNo[effRunSize - 1]    = totEvtNo[effRunSize - 1] + ( EvtMax - totSimEvt );
          std::cout << "Additional " << EvtMax - totSimEvt
                    << " events need to be simulated in run "
                    << m_runList[m_runList.size() - 1] << std::endl;
        }
        else
        {
          bool found = false;
          for ( unsigned int i = 0; i < lumi.size(); i++ )
          {
            if ( lumi[i] > 0 )
            {
              m_evtNoList.push_back( EvtMax - totSimEvt );
              m_runList.push_back( tmp_runList[i] );
              totEvtNo.push_back( tmp_runList[i] );
              totEvtNo.push_back( EvtMax - totSimEvt );
              found = true;
              std::cout << "The max event number maybe too small, all " << evtSubNo
                        << " events need to be simulated in run " << tmp_runList[i]
                        << std::endl;
              break;
            }
          }
          if ( !found )
          {
            std::cerr << "ERORR: "
                      << "Total luminosity is ZERO, please check your run list. " << std::endl;
            std::exit( 0 );
          }
        }
      }
    }
    else
    {
      std::cerr << "ERORR: "
                << "Total luminosity is ZERO!!! Please check your run list. " << std::endl;
      std::exit( 0 );
    }
    m_RealizationSvc->setRunId( -std::abs( m_runList[0] ) );
    m_RealizationSvc->setRunEvtNum( m_evtNoList[0] );
    m_initRunNo = -std::abs( m_runList[0] );
    jobInfoSvc->setTotEvtNo( totEvtNo );
  }
  else { m_initRunNo = -std::abs( ( m_RealizationSvc->getRunList() )[0] ); }
}

McContext() [2/2]

McContext::McContext	(	const McEventSelector *	pSelector,
		const int &	initRunNo,
		const unsigned int &	initEventNo,
		const unsigned int &	evPR )

Definition at line 235 of file McEventSelector.cxx.

    : m_pSelector( pSelector )
    , m_initRunNo( initRunNo )
    , m_initEventNo( initEventNo )
    , m_eventsP( evPR )
    , m_firstEvent( true )
    , m_lumiBlockNo( 0 ) {}

~McContext()

McContext::~McContext ( )

virtual

Standard destructor.

Definition at line 244 of file McEventSelector.cxx.

{}

Member Function Documentation

eventNumber()

unsigned int McContext::eventNumber ( ) const

inline

Definition at line 105 of file McEventSelector.cxx.

{ return m_eventNo; }

Referenced by McEventSelector::createAddress().

identifier()

virtual void * McContext::identifier ( ) const

inlinevirtual

Definition at line 67 of file McEventSelector.cxx.

{ return (void*)m_pSelector; }

next()

void McContext::next ( )

inline

Definition at line 80 of file McEventSelector.cxx.

              {
    if ( m_firstEvent )
    {
      m_firstEvent = false;
 
      m_eventNo = m_initEventNo;
      m_runNo   = m_initRunNo;
 
      return;
    }
    m_eventNo++;
    if ( m_RealizationSvc->UseDBFlag() == false ) return;
    if ( m_eventNo >= ( m_initEventNo + m_evtNoList[m_lumiBlockNo] ) )
    {
      m_lumiBlockNo++;
      m_runNo = -std::abs( m_runList[m_lumiBlockNo] );
      m_RealizationSvc->setRunId( m_runNo );
      m_RealizationSvc->setRunEvtNum( m_evtNoList[m_lumiBlockNo] );
      m_eventNo = m_initEventNo;
    }
  }

Referenced by McEventSelector::next().

previous()

void McContext::previous ( )

inline

Definition at line 102 of file McEventSelector.cxx.

{ m_eventNo--; }

Referenced by McEventSelector::previous().

rewind()

void McContext::rewind ( )

inline

Definition at line 70 of file McEventSelector.cxx.

                {
 
    // Question: should this rewind to before the first event, or to
    // it? ie, will next() be called right after rewind()?
    // if so, then should set m_firstEvent to true;
 
    m_eventNo = m_initEventNo;
    m_runNo   = m_initRunNo;
  }

Referenced by McEventSelector::rewind().

runNumber()

int McContext::runNumber ( ) const

inline

Definition at line 104 of file McEventSelector.cxx.

{ return m_runNo; }

Referenced by McEventSelector::createAddress().

setCriteria()

void McContext::setCriteria ( const std::string & crit )

inline

Definition at line 68 of file McEventSelector.cxx.

{ m_criteria = crit; }

setEventNumber()

void McContext::setEventNumber ( unsigned int eventNo )

inline

Definition at line 108 of file McEventSelector.cxx.

{ m_eventNo = eventNo; }

setRunNumber()

void McContext::setRunNumber ( int runNo )

inline

Definition at line 107 of file McEventSelector.cxx.

{ m_runNo = runNo; }

---

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

• McEventSelector.cxx