CgemAliEvent.cxx

Go to the documentation of this file.

#include "CgemAlignAlg/CgemAliEvent.h"
 
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/IMessageSvc.h"
#include "GaudiKernel/StatusCode.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/PropertyMgr.h"
 
#include "BesTimerSvc/IBesTimerSvc.h"
#include "BesTimerSvc/BesTimerSvc.h"
 
#include "EventModel/Event.h"
#include "RawEvent/RawDataUtil.h"
#include "EventModel/Event.h"
 
#include "CgemCalibFunSvc/ICgemCalibFunSvc.h"
#include "CgemCalibFunSvc/CgemCalibFunSvc.h"
 
#include "EvTimeEvent/RecEsTime.h"
#include "MdcRecEvent/RecMdcTrack.h"
#include "MdcRecEvent/RecMdcKalTrack.h"
#include "MdcRecEvent/RecMdcHit.h"
#include "ReconEvent/ReconEvent.h"
 
#include "Identifier/Identifier.h"
#include "Identifier/CgemID.h"
 
#include <iostream>
#include <math.h>
 
using namespace Event;
 
StatusCode CgemAliEvent::setRecEvent(){
     IMessageSvc *msgSvc;
     Gaudi::svcLocator()->service("MessageSvc", msgSvc);
     MsgStream log(msgSvc, "CgemAliEvent");
     log << MSG::INFO << "CgemAliEvent::setRecEvent" << endreq;
 
     IDataProviderSvc* eventSvc = NULL;
     Gaudi::svcLocator()->service("EventDataSvc", eventSvc);
 
     SmartDataPtr<RecMdcTrackCol> newtrkCol(eventSvc, "/Event/Recon/RecMdcTrackCol");
     if(!newtrkCol){
      //log << MSG::ERROR << "Could not find RecMdcTrackCol" << endreq;
      return ( StatusCode::FAILURE );
     }
 
     int i = 0;
     CgemAliRecTrk* rectrk;
     RecMdcTrackCol::iterator it_trk = newtrkCol->begin();
     for(; it_trk != newtrkCol->end(); it_trk++){
      rectrk = new CgemAliRecTrk();
      rectrk -> setRecTrk(it_trk);
      m_rectrk.push_back(rectrk);
 
      i++;
     }
     log << MSG::WARNING << "RecTrack ntrk = " << i << endreq;
 
     return StatusCode::SUCCESS;
}
 
StatusCode CgemAliEvent::setKalEvent(){
     IMessageSvc *msgSvc;
     Gaudi::svcLocator()->service("MessageSvc", msgSvc);
     MsgStream log(msgSvc, "CgemAliEvent");
     log << MSG::INFO << "CgemAliEvent::setKalEvent" << endreq;
 
     IDataProviderSvc* eventSvc = NULL;
     Gaudi::svcLocator()->service("EventDataSvc", eventSvc);
 
 
     SmartDataPtr<RecMdcKalTrackCol> kaltrkCol(eventSvc,"/Event/Recon/RecMdcKalTrackCol");
     if (!kaltrkCol) {
      log << MSG::FATAL << "Could not find RecMdcKalTrackCol" << endreq;
      return StatusCode::FAILURE;
     }
 
     int i = 0;
     CgemAliRecTrk* rectrk;
     RecMdcKalTrackCol::iterator iter_trk = kaltrkCol->begin();
     for(; iter_trk != kaltrkCol->end(); iter_trk++) {
      rectrk = new CgemAliRecTrk();
      rectrk -> setKalTrk(iter_trk);
      m_rectrk.push_back(rectrk);
 
      i++;
     }
     log << MSG::INFO << "KalTrack ntrk = " << i << endreq;
 
     return StatusCode::SUCCESS;
}
 
void CgemAliEvent::clear(){
     unsigned int i;
     for(i=0; i<m_rectrk.size(); i++){
      delete m_rectrk[i];
     }
     m_rectrk.clear();
}