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 45 of file McEventSelector.cxx.

Constructor & Destructor Documentation

McContext() [1/2]

McContext::McContext

(

const McEventSelector *

pSelector

)

Standard constructor.

Definition at line 124 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();
  IRealizationSvc *tmpReal;
  StatusCode status = svcLocator->service("RealizationSvc", tmpReal);
  if (!status.isSuccess())
  {
     std::cout << " Could not initialize Realization Service" << std::endl;
  } else {
    m_RealizationSvc=dynamic_cast<RealizationSvc*>(tmpReal);
  }
  //get jobSvc
  IDataInfoSvc* tmpInfoSvc;
  DataInfoSvc* jobInfoSvc;
  status = svcLocator->service("DataInfoSvc",tmpInfoSvc);
  if (status.isSuccess()) {
    std::cout << "get the DataInfoSvc" << std::endl;
    jobInfoSvc=dynamic_cast<DataInfoSvc *>(tmpInfoSvc);
  }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 242 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 252 of file McEventSelector.cxx.

                      {
}

Member Function Documentation

eventNumber()

unsigned int McContext::eventNumber ( ) const

inline

Definition at line 117 of file McEventSelector.cxx.

{ return m_eventNo; }

Referenced by McEventSelector::createAddress().

identifier()

virtual void * McContext::identifier ( ) const

inlinevirtual

Definition at line 75 of file McEventSelector.cxx.

                                   {
    return (void*)m_pSelector;
  }

next()

void McContext::next ( )

inline

Definition at line 92 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 112 of file McEventSelector.cxx.

                  {
    m_eventNo--;
  }

Referenced by McEventSelector::previous().

rewind()

void McContext::rewind ( )

inline

Definition at line 82 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 116 of file McEventSelector.cxx.

{ return m_runNo; }

Referenced by McEventSelector::createAddress().

setCriteria()

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

inline

Definition at line 78 of file McEventSelector.cxx.

                                          {
    m_criteria = crit;
  }

setEventNumber()

void McContext::setEventNumber ( unsigned int  eventNo )

inline

Definition at line 120 of file McEventSelector.cxx.

{ m_eventNo = eventNo ; }

setRunNumber()

void McContext::setRunNumber ( int  runNo )

inline

Definition at line 119 of file McEventSelector.cxx.

{ m_runNo = runNo; }

---

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

• McEventSelector.cxx