---

CreateEvent.cpp

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/GlastSvc/src/test/CreateEvent.cpp,v 1.17 2001/04/19 01:32:31 igable Exp $
#define GlastApps_CreateEvent_CPP 


//#define GlastApps_CreateEvent_CPP 


// Include files
#include "CreateEvent.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/Irf/IrfCalHit.h"
#include "GlastEvent/Irf/IrfTkrLayer.h"

#include "GlastEvent/MonteCarlo/McIntegratingHit.h"
                                
static const AlgFactory<CreateEvent>  Factory;
const IAlgFactory& CreateEventFactory = Factory;

//extern void GlastSvc_load();

//------------------------------------------------------------------------------
CreateEvent::CreateEvent(const std::string& name, ISvcLocator* pSvcLocator) :
Algorithm(name, pSvcLocator), m_detSvc(0), m_irfLoadSvc(0) {
    //GlastSvc_load();
}


//------------------------------------------------------------------------------
StatusCode CreateEvent::initialize() {
    
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "initialize" << endreq;
    
    // Use the Job options service to set the Algorithm's parameters
    setProperties();
    
    // now try to find the GlastDevSvc service
    IGlastDetSvc* detSvc = 0;
    const IID&  IID_IGlastDetSvc  =  401; // wired it for now!

    StatusCode sc = serviceLocator()->getService ("GlastDetSvc",
        IID_IGlastDetSvc, reinterpret_cast<IInterface*&>( detSvc ));
    
    if (sc.isSuccess ())
        log << MSG::INFO << "Succeeded in accessing the GlastDetSvc!" << endreq;
    else {
        log << MSG::ERROR << "Couldn't find the GlastDetSvc!" << endreq;
        return StatusCode::FAILURE;
    }
    m_detSvc = detSvc;


    return StatusCode::SUCCESS;
}


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

    DataObject* pObject;

    sc = eventSvc()->retrieveObject("/Event/Irf/IrfAcdHits", pObject);

    if( sc.isFailure() ){
        log << MSG::INFO << "Failed to get IrfAcdHits from TDS!" << endreq;
        return sc; 
    }
    
        log << MSG::INFO << "Successfully retrieved ACD Container!!!" << endreq;

    ObjectVector<IrfAcdHit>* acdList;
    try {
        acdList  = dynamic_cast<ObjectVector<IrfAcdHit>*>(pObject);
    } catch(...) {
        log << MSG::INFO << "Failed to convert object to IrfAcdHitVector" << endreq;
        return StatusCode::FAILURE;
    }

    
        // Decrease the amount of output
#if 0

    for (ObjectVector<IrfAcdHit>::const_iterator it = acdList->begin(); it != acdList->end(); it++) {
        log << MSG::INFO << " ACD Tile Hit " 
            << (*it)->id() << " "
            << (*it)->energy() << endreq;
    }
#endif

    sc = eventSvc()->retrieveObject("/Event/Irf/IrfCalHits", pObject);
    if( sc.isFailure() ) {
        log << MSG::INFO << "Failed to get IrfCalHits" << endreq;
        return sc;                                                             
    } 
    
    log << MSG::INFO << "Retrieved IrfCalHit Vector!" << endreq;

    ObjectVector<IrfCalHit>* calList;
    try {
        calList  = dynamic_cast<ObjectVector<IrfCalHit>*>(pObject);
    } catch(...) {
        log << MSG::INFO << "Failed to convert object to IrfCalHitVector" << endreq;
        return StatusCode::FAILURE;
    }

        // Decrease amount of output
#if 0

    for (ObjectVector<IrfCalHit>::const_iterator xtal = calList->begin(); xtal != calList->end(); xtal++) {
        log << MSG::INFO << "Csi Log Hit: " << endreq;
        log <<MSG::INFO << "       " 
                        << (*xtal)->minusResponse() << " " 
                        << (*xtal)->plusResponse() << " "
                        << (*xtal)->energy() << " " << endreq;
    }
#endif

    sc = eventSvc()->retrieveObject("/Event/Irf/IrfTkrHits", pObject);
    if( sc.isFailure() ) return sc;                                                             
    
        log << MSG::INFO << "Successfully retrieved TKR Container!!!" << endreq;

    ObjectVector<IrfTkrLayer>* tkrList;
    try {
        tkrList  = dynamic_cast<ObjectVector<IrfTkrLayer>*>(pObject);
    } catch(...) {
        log << MSG::INFO << "Failed to convert object to IrfTkrLayerVector" << endreq;
        return StatusCode::FAILURE;
    }
    
        // Decrease amount of output
#if 0

    for (ObjectVector<IrfTkrLayer>::const_iterator silayer = tkrList->begin(); silayer != tkrList->end(); silayer++) {
        log << MSG::INFO << "TkrLayer " << (*silayer)->id() << " "
            << (*silayer)->MaxEnergy() << endreq;

        for (ObjectVector<IrfTkrHit>::const_iterator hit = ((*silayer)->getHits())->begin(); hit != ((*silayer)->getHits())->end(); hit++) {
            log << MSG::INFO << "Hit SSD Strip " 
                << (*hit)->id() << " "
                << (*hit)->energy() << " "
                << (*hit)->noise() << endreq;
        }
      }
#endif // End of decreased output comment.



    return sc;


    //
}


//------------------------------------------------------------------------------
StatusCode CreateEvent::finalize() {
    
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "finalize" << endreq;
    
    return StatusCode::SUCCESS;
}
;

StatusCode CreateEvent::testMcClass(){
    McIntegratingHit* integratingHit = new McIntegratingHit();
    McParticle* mcParticle =  new McParticle();

    
    return StatusCode::SUCCESS;
}

---