---

AcdDigiAlg.cpp

Go to the documentation of this file.

#define GlastDigi_AcdDigiAlg_CPP 

// Include files
#include "AcdDigiAlg.h"

#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/SmartDataPtr.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"

#include <algorithm>
#include <math.h>

#include "xml/IFile.h"

static const AlgFactory<AcdDigiAlg>  Factory;
const IAlgFactory& AcdDigiAlgFactory = Factory;

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

    declareProperty ("xmlFile", m_xmlFile);
}


//------------------------------------------------------------------------------
StatusCode AcdDigiAlg::initialize() {
    
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "initialize" << endreq;
    
    // Use the Job options service to set the Algorithm's parameters
    setProperties();
    
    xml::IFile m_ifile(m_xmlFile.c_str());
    double lowThreshold = m_ifile.getDouble("acd", "lowThreshold");
    double vetoThreshold = m_ifile.getDouble("acd", "vetoThreshold");
    double highThreshold = m_ifile.getDouble("acd", "highThreshold");
    double adcChannelsPerMeV = m_ifile.getDouble("acd", "adcChannelsPerMeV");

    // Convert from MeV to GeV
    m_lowThresholdGeV = lowThreshold / 1000.f;
    m_vetoThresholdGeV = vetoThreshold / 1000.f;
    m_highThresholdGeV = highThreshold / 1000.f;
    m_adcChannelsPerGeV = adcChannelsPerMeV * 1000.0f;

    return StatusCode::SUCCESS;
}


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

    SmartDataPtr<IrfAcdHitVector> IrfAcd(eventSvc(),"/Event/Irf/IrfAcdHits");

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

    
    AcdDigiVector *DigiAcd = new AcdDigiVector;
    for (ObjectVector<IrfAcdHit>::const_iterator it = IrfAcd->begin(); it != IrfAcd->end(); it++) {
        
        bool overLow = false;
        bool overVeto = false;
        bool overHigh = false;

        double energyDeposited = (*it)->energy();  // GeV

        if (energyDeposited < m_lowThresholdGeV) continue;
        overLow = true;
        if (energyDeposited > m_vetoThresholdGeV) overVeto = true;
        if (energyDeposited > m_highThresholdGeV) overHigh = true;
        
        int pha = static_cast<int>(std::min((float)floor(energyDeposited * m_adcChannelsPerGeV), 4095.0f)) ;

        DigiAcd->push_back(new AcdDigi(idents::AcdId((*it)->id()), pha, overVeto, overLow, overHigh) );
    }

    sc = eventSvc()->registerObject("/Event/Digi/AcdDigis", DigiAcd);

    if (sc != StatusCode::SUCCESS) log << MSG::INFO << "Failed to register ACD digi vector" << endreq;


    return sc;
}


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

---