---

ReadDigi.cpp

Go to the documentation of this file.

#define GlastDigi_ReadDigi_CPP 

// Include files
#include "ReadDigi.h"

#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/SmartDataPtr.h"

#include "GlastSvc/GlastDetSvc/IGlastDetSvc.h"
#include "GlastEvent/TopLevel/EventModel.h"
#include "GlastEvent/TopLevel/Event.h"
#include "GaudiKernel/ObjectVector.h"

#include "GlastEvent/TopLevel/IrfEvent.h"
#include "GlastEvent/Irf/IrfAcdHit.h"

#include "GlastEvent/Digi/AcdDigi.h"


static const AlgFactory<ReadDigi>  Factory;
const IAlgFactory& ReadDigiFactory = Factory;

//------------------------------------------------------------------------------
ReadDigi::ReadDigi(const std::string& name, ISvcLocator* pSvcLocator) :
Algorithm(name, pSvcLocator) {

}


//------------------------------------------------------------------------------
StatusCode ReadDigi::initialize() {
    
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "initialize" << endreq;
    
    // Use the Job options service to set the Algorithm's parameters
    setProperties();

    return StatusCode::SUCCESS;
}


//------------------------------------------------------------------------------
StatusCode ReadDigi::execute() {
    
    StatusCode  sc = StatusCode::SUCCESS;
    MsgStream   log( msgSvc(), name() );
    log << MSG::INFO << "execute" << endreq;

    SmartDataPtr<AcdDigiVector> DigiAcd(eventSvc(),"/Event/Digi/AcdDigis");

    if (DigiAcd == 0) {
        log << MSG::INFO << "no acd data is available" << endreq;
        return sc;
    }

    // Otherwise there was ACD data for this event

    // Iterate over all hit ACD tiles
    for (ObjectVector<AcdDigi>::const_iterator it = DigiAcd->begin(); it != DigiAcd->end(); it++) {
        const unsigned int myId = int((*it)->ID());
        log << MSG::INFO << " ACD Digi Id = " 
            << (*it)->ID() << " "
            << " PHA = " << (*it)->PulseHeight() 
            << " LowThres = " << (*it)->LowThreshold()
            << " VetoThres = " << (*it)->Veto()
            << " HighThres = " << (*it)->HighThreshold()
            << endreq;
    }


    return sc;
}


//------------------------------------------------------------------------------
/*
Finalize is called once at the end of event processing.
Any cleanup to occur - will occur here.
*/
StatusCode ReadDigi::finalize() {
    
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "finalize" << endreq;
    
    return StatusCode::SUCCESS;
}

---