RecCgemHitCnv.cxx

Go to the documentation of this file.

#ifndef RecCgemHitCnv_CXX
#define RecCgemHitCnv_CXX  1
 
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/ObjectVector.h"
 
#include "TClonesArray.h"
 
#include "EventModel/EventModel.h"
 
#include "ReconEvent/ReconEvent.h"
#include "CgemRecEvent/RecCgemHit.h"
#include "Identifier/Identifier.h"
 
#include "RootEventData/TRecCgemHit.h" // standard root object
#include "RootEventData/TRecTrackEvent.h"
 
#include "RootCnvSvc/Rec/RecTrackCnv.h"
#include "RootCnvSvc/Rec/RecCgemHitCnv.h"
#include "RootCnvSvc/RootAddress.h"
 
#include <vector>
 
using namespace std;
 
 
RecCgemHitCnv::RecCgemHitCnv(ISvcLocator* svc)
: RootEventBaseCnv(classID(), svc)
{
 
    // Here we associate this converter with the /Event path on the TDS.
    MsgStream log(msgSvc(), "RecCgemHitCnv");
    //log << MSG::DEBUG << "Constructor called for " << objType() << endreq;
    //m_rootTreename ="Rec";
    m_rootBranchname ="m_recCgemHitCol";
    //declareObject(EventModel::Recon::RecCgemHitCol, objType(), m_rootTreename, m_rootBranchname);
    m_adresses.push_back(&m_recCgemHitCol);
    m_recCgemHitCol=0;
}
 
StatusCode RecCgemHitCnv::TObjectToDataObject(DataObject*& refpObject) {
  // creation of TDS object from root object
  MsgStream log(msgSvc(), "RecCgemHitCnv");
  log << MSG::DEBUG << "RecCgemHitCnv::TObjectToDataObject" << endreq;
 
  // create the TDS location for the RecCgemHit Collection
  RecCgemHitCol* recCgemHitCol = new RecCgemHitCol;
  refpObject = recCgemHitCol;
 
  // now convert
  if (!m_recCgemHitCol) return StatusCode::SUCCESS;
  TIter cgemHitIter(m_recCgemHitCol);
  TRecCgemHit *recCgemHitRoot = 0;
  while ((recCgemHitRoot = (TRecCgemHit*)cgemHitIter.Next())) {
    int      id = recCgemHitRoot->getId();
    int      trackId = recCgemHitRoot->getTrkId();
    double   chisqAdd = recCgemHitRoot->getChisqAdd();
    int      stat = recCgemHitRoot->getStat(); 
//     unsigned int cgemIdMSB=recCgemHitRoot->getCgemIdMSB(); 
//     unsigned int cgemIdLSB=recCgemHitRoot->getCgemIdLSB(); 
    unsigned int cgemId=recCgemHitRoot->getCgemId(); 
    double   tdc = recCgemHitRoot->getTdc();
    double   adc = recCgemHitRoot->getAdc();
    double   doca = recCgemHitRoot->getDoca();
    double   entra = recCgemHitRoot->getEntra();
    double   zhit = recCgemHitRoot->getZhit();
    double   fltLen = recCgemHitRoot->getFltLen();
 
    RecCgemHit *recCgemHit= new RecCgemHit();
    m_common.m_rootRecCgemHitMap[recCgemHitRoot] = recCgemHit;
 
    recCgemHit->setId( id );
    recCgemHit->setTrkId( trackId );
    recCgemHit->setChisqAdd( chisqAdd );
    recCgemHit->setStat( stat );
//     recCgemHit->setCgemId( Identifier64(cgemIdMSB,cgemIdLSB) );
    recCgemHit->setCgemId( Identifier(cgemId) );
    recCgemHit->setTdc( tdc );
    recCgemHit->setAdc( adc );
    recCgemHit->setDoca( doca );
    recCgemHit->setEntra( entra );
    recCgemHit->setZhit( zhit );
    recCgemHit->setFltLen( fltLen );
 
    recCgemHitCol->push_back(recCgemHit);
  }
    return StatusCode::SUCCESS;
}
 
StatusCode RecCgemHitCnv::DataObjectToTObject(DataObject* obj,RootAddress* rootaddr) {
 
  MsgStream log(msgSvc(), "RecCgemHitCnv");
  log << MSG::DEBUG << "RecCgemHitCnv::DataObjectToTObject" << endreq;
 
  RecCgemHitCol * cgemHitColTds=dynamic_cast<RecCgemHitCol *> (obj);
  if (!cgemHitColTds) {
    log << MSG::ERROR << "Could not downcast to RecCgemHitCol" << endreq;
    return StatusCode::FAILURE;
  }
 
  DataObject *evt;
  m_eds->findObject(EventModel::Recon::Event,evt);
  if (evt==NULL) {
    log << MSG::ERROR << "Could not get ReconEvent in TDS "  << endreq;
    return StatusCode::FAILURE;
  }
 
  ReconEvent * devtTds=dynamic_cast<ReconEvent *> (evt);
  if (!devtTds) {
    log << MSG::ERROR << "RecCgemHitCnv:Could not downcast to TDS ReconEvent" << endreq;
  }
 
  IOpaqueAddress *addr;
 
  m_cnvSvc->getRecTrackCnv()->createRep(evt,addr); 
  TRecTrackEvent *recEvt=m_cnvSvc->getRecTrackCnv()->getWriteObject();
 
  const TObjArray *m_recCgemHitCol = recEvt->getRecCgemHitCol();
  if (!m_recCgemHitCol) return StatusCode::SUCCESS;
  recEvt->clearRecCgemHitCol(); //necessary in case there is I/O at the same time since array is static
 
  RecCgemHitCol::const_iterator recCgemHit;
  for (recCgemHit = cgemHitColTds->begin(); recCgemHit != cgemHitColTds->end(); recCgemHit++) {
     Bool_t     isGrouped =(*recCgemHit)->isGrouped();
     Int_t      id =(*recCgemHit)->getId();
     Int_t      trackId =(*recCgemHit)->getTrkId();
     Double_t   chisqAdd =(*recCgemHit)->getChisqAdd();
     Int_t      stat =(*recCgemHit)->getStat(); 
//      UInt_t     cgemIdMSB =(*recCgemHit)->getCgemId().get_valueMSB(); 
//      UInt_t     cgemIdLSB =(*recCgemHit)->getCgemId().get_valueLSB(); 
     UInt_t     cgemId =(*recCgemHit)->getCgemId().get_value(); 
     Double_t   tdc =(*recCgemHit)->getTdc();
     Double_t   adc =(*recCgemHit)->getAdc();
     Double_t   doca =(*recCgemHit)->getDoca();
     Double_t   entra =(*recCgemHit)->getEntra();
     Double_t   zhit =(*recCgemHit)->getZhit();
     Double_t   fltLen =(*recCgemHit)->getFltLen();
 
     TRecCgemHit *recCgemHitRoot = new TRecCgemHit();
     //m_common.m_recCgemHitMap[(*recCgemHit)] = recCgemHitRoot;
 
     recCgemHitRoot->setIsGrouped( isGrouped );
     recCgemHitRoot->setId( id );
     recCgemHitRoot->setTrkId( trackId );
     recCgemHitRoot->setChisqAdd( chisqAdd );
     recCgemHitRoot->setStat( stat );
//      recCgemHitRoot->setCgemIdMSB( cgemIdMSB );
//      recCgemHitRoot->setCgemIdLSB( cgemIdLSB );
     recCgemHitRoot->setCgemId( cgemId );
     recCgemHitRoot->setTdc( tdc );
     recCgemHitRoot->setAdc( adc );
     recCgemHitRoot->setDoca( doca );
     recCgemHitRoot->setEntra( entra );
     recCgemHitRoot->setZhit( zhit );
     recCgemHitRoot->setFltLen( fltLen );
 
     recEvt->addRecCgemHit(recCgemHitRoot);
  }
  return StatusCode::SUCCESS;
}
#endif