digiRootReaderAlg Class Reference

Reads Digitization data from a persistent ROOT file and stores the the data in the TDS. Based on digiRootReaderAlg of Glast. More...

#include <digiRootReaderAlg.h>

Inheritance diagram for digiRootReaderAlg:

Public Member Functions
	digiRootReaderAlg (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
	Handles setup by opening ROOT file in read mode and creating a new TTree.
StatusCode	execute ()
	Orchastrates reading from ROOT file and storing the data on the TDS for each event.
StatusCode	finalize ()
	Closes the ROOT file and cleans up.

Detailed Description

Reads Digitization data from a persistent ROOT file and stores the the data in the TDS. Based on digiRootReaderAlg of Glast.

Definition at line 35 of file digiRootReaderAlg.h.

Constructor & Destructor Documentation

digiRootReaderAlg()

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

Definition at line 87 of file digiRootReaderAlg.h.

                                                                                    : 
Algorithm(name, pSvcLocator)
{
    // Input pararmeters that may be set via the jobOptions file
    // Input ROOT file name
    declareProperty("digiRootFile",m_fileName="");
    StringArrayProperty initList;
    std::vector<std::string> initVec;
    initVec.push_back("digicopy.root");
    initList.setValue(initVec);
    declareProperty("digiRootFileList",m_fileList=initList);
    // Input TTree name
    initVec.clear();
    declareProperty("digiTreeName", m_treeName="Rec");// wensp midify for test 2005/05/14
 
}

Member Function Documentation

execute()

StatusCode digiRootReaderAlg::execute

(

)

Orchastrates reading from ROOT file and storing the data on the TDS for each event.

Definition at line 177 of file digiRootReaderAlg.h.

{
    // Purpose and Method:  Called once per event.  This method calls
    //   the appropriate methods to read data from the ROOT file and store
    //   data on the TDS.
 
    MsgStream log(msgSvc(), name());
 
    StatusCode sc = StatusCode::SUCCESS;
    
    if (m_digiEvt) m_digiEvt->Clear();
 
    static Int_t evtId = 0;
    int readInd, numBytes;
    std::pair<int,int> runEventPair = (m_rootIoSvc) ? m_rootIoSvc->runEventPair() : std::pair<int,int>(-1,-1);
    
    if ((m_rootIoSvc) && (m_rootIoSvc->index() >= 0)) {
        readInd = m_rootIoSvc->index();
    } else if ((m_rootIoSvc) && (runEventPair.first != -1) && (runEventPair.second != -1)) {
        int run = runEventPair.first;
        int evt = runEventPair.second;
        readInd = m_digiTree->GetEntryNumberWithIndex(run, evt);
    } else {
        readInd = evtId;
    }
 
    if (readInd >= m_numEvents) {
        log << MSG::WARNING << "Requested index is out of bounds - no digi data loaded" << endreq;
        return StatusCode::SUCCESS;
    }
 
    numBytes = m_digiTree->GetEvent(readInd);
    
    if ((numBytes <= 0) || (!m_digiEvt)) {
            log << MSG::WARNING << "Failed to load digi event" << endreq;
            return StatusCode::SUCCESS;
    }
 
 
    sc = readDigiEvent();
    if (sc.isFailure()) {
        log << MSG::ERROR << "Failed to read top level DigiEvent" << endreq;
        return sc;
    }
 
    sc = readMdcDigi();
    if (sc.isFailure()) {
        log << MSG::ERROR << "Failed to load MdcDigi" << endreq;
        return sc;
    }
 
    sc = readCgemDigi();
    if (sc.isFailure()) {
        log << MSG::ERROR << "Failed to load CgemDigi" << endreq;
        return sc;
    }
 
    evtId = readInd+1;
    return sc;
}

finalize()

StatusCode digiRootReaderAlg::finalize

(

)

Closes the ROOT file and cleans up.

Definition at line 343 of file digiRootReaderAlg.h.

{
    close();
    
    StatusCode sc = StatusCode::SUCCESS;
    return sc;
}

initialize()

StatusCode digiRootReaderAlg::initialize

(

)

Handles setup by opening ROOT file in read mode and creating a new TTree.

Definition at line 104 of file digiRootReaderAlg.h.

{
    // Purpose and Method:  Called once before the run begins.  This method
    //    opens a new ROOT file and prepares for reading.
 
    StatusCode sc = StatusCode::SUCCESS;
    MsgStream log(msgSvc(), name());
    
    // Use the Job options service to set the Algorithm's parameters
    // This will retrieve parameters set in the job options file
    setProperties();
   
    if ( service("RootIoSvc", m_rootIoSvc, true).isFailure() ){
        log << MSG::INFO << "Couldn't find the RootIoSvc!" << endreq;
        log << MSG::INFO << "Event loop will not terminate gracefully" << endreq;
        m_rootIoSvc = 0;
        //return StatusCode::FAILURE;
    } 
 
    facilities::Util::expandEnvVar(&m_fileName);
    
    // Save the current directory for the ntuple writer service
    TDirectory *saveDir = gDirectory;   
 
    m_digiTree = new TChain(m_treeName.c_str());
 
    std::string emptyStr("");
    if (m_fileName.compare(emptyStr) != 0) {
      TFile f(m_fileName.c_str());
      if (!f.IsOpen()) {
        log << MSG::ERROR << "ROOT file " << m_fileName.c_str()
            << " could not be opened for reading." << endreq;
        return StatusCode::FAILURE;
      }
      f.Close();
      m_digiTree->Add(m_fileName.c_str());
          log << MSG::INFO << "Opened file: " << m_fileName.c_str() << endreq;
    } else {
      const std::vector<std::string> fileList = m_fileList.value( );
      std::vector<std::string>::const_iterator it;
      std::vector<std::string>::const_iterator itend = fileList.end( );
      for (it = fileList.begin(); it != itend; it++) {
        std::string theFile = (*it);
        TFile f(theFile.c_str());
        if (!f.IsOpen()) {
          log << MSG::ERROR << "ROOT file " << theFile.c_str()
              << " could not be opened for reading." << endreq;
          return StatusCode::FAILURE;
        }
      f.Close();
      m_digiTree->Add(theFile.c_str());
          log << MSG::INFO << "Opened file: " << theFile.c_str() << endreq;
      }
    }
 
 
    m_digiEvt = 0;
    m_digiTree->SetBranchAddress("DigiEvent", &m_digiEvt);
    //m_common.m_digiEvt = m_digiEvt;
    m_numEvents = m_digiTree->GetEntries();
    
    if (m_rootIoSvc) {
      m_rootIoSvc->setRootEvtMax(m_numEvents);
      if (!m_digiTree->GetIndex()) m_digiTree->BuildIndex("m_runId", "m_eventId");
      m_rootIoSvc->registerRootTree(m_digiTree);
    }
 
 
    saveDir->cd();
    return sc;
    
}

---

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

• digiRootReaderAlg.h