RecCgemDedxHitCnv.cxx

Go to the documentation of this file.

#ifndef RecCgemDedxHitCnv_CXX
#define RecCgemDedxHitCnv_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/RecCgemDedxHit.h"
#include "Identifier/Identifier.h"
 
#include "RootEventData/TRecCgemDedxHit.h" // standard root object
#include "RootEventData/TRecTrackEvent.h"
 
#include "RootCnvSvc/Rec/RecTrackCnv.h"
#include "RootCnvSvc/Rec/RecCgemDedxHitCnv.h"
#include "RootCnvSvc/RootAddress.h"
 
#include <vector>
 
using namespace std;
 
// Instantiation of a static factory class used by clients to create
// instances of this service
//static CnvFactory<RecCgemDedxCnv> s_factory;
//const ICnvFactory& RecCgemDedxCnvFactory = s_factory;
 
RecCgemDedxHitCnv::RecCgemDedxHitCnv(ISvcLocator* svc)
: RootEventBaseCnv(classID(), svc)
{
 
    // Here we associate this converter with the /Event path on the TDS.
    MsgStream log(msgSvc(), "RecCgemDedxHitCnv");
    //log << MSG::DEBUG << "Constructor called for " << objType() << endreq;
    //m_rootTreename ="Rec";
    m_rootBranchname ="m_recCgemDedxHitCol";
    // declareObject(EventModel::Recon::RecCgemDedxCol, objType(), m_rootTreename, m_rootBranchname);
    m_adresses.push_back(&m_recCgemDedxHitCol);
    m_recCgemDedxHitCol=0;
}
 
StatusCode RecCgemDedxHitCnv::TObjectToDataObject(DataObject*& refpObject) {
  // creation of TDS object from root object
    MsgStream log(msgSvc(), "RecCgemDedxHitCnv");
    log << MSG::DEBUG << "RecCgemDedxHitCnv::TObjectToDataObject" << endreq;
    StatusCode sc=StatusCode::SUCCESS;
    // create the TDS location for the Dedx Collection
    RecCgemDedxHitCol* recCgemDedxHitCol = new RecCgemDedxHitCol;
    refpObject=recCgemDedxHitCol;
 
    // now convert
    if (!m_recCgemDedxHitCol) return sc;
    TIter dedxIter(m_recCgemDedxHitCol);
    TRecCgemDedxHit *recCgemDedxHitRoot = 0;
    while ((recCgemDedxHitRoot = (TRecCgemDedxHit*)dedxIter.Next())) {
      int trackId = recCgemDedxHitRoot->trkId();
//       unsigned int cgemIdMSB(recCgemDedxHitRoot->cgemIdMSB());
//       unsigned int cgemIdLSB(recCgemDedxHitRoot->cgemIdLSB());
      unsigned int cgemId(recCgemDedxHitRoot->cgemId());
      double pathlength =  recCgemDedxHitRoot->pathLength();
     log << MSG::DEBUG<<"check Reconstrunction of dE/dx Hit root"<<" pathlength is : "<<pathlength<<endreq;
     RecCgemDedxHit *recCgemDedxHit = new RecCgemDedxHit();
     m_common.m_rootRecCgemDedxHitMap[recCgemDedxHitRoot] = recCgemDedxHit;
     recCgemDedxHit->setTrkId(trackId);
//      recCgemDedxHit->setCgemId( Identifier64(cgemIdMSB,cgemIdLSB) );
     recCgemDedxHit->setCgemId( Identifier(cgemId) );
     recCgemDedxHit->setPathLength( pathlength);
  //   log << MSG::DEBUG<<"check Reconstrunction of dE/dx Hit root"<<" pathlength 2 is : "<<(*recdedxhit)->getPathLength()<<endreq; 
     recCgemDedxHitCol->push_back(recCgemDedxHit); 
     //delete dedxTds;
     // dedxTds = NULL;
     }
    //m_dedxCol->Delete();
    delete m_recCgemDedxHitCol;
    m_recCgemDedxHitCol = 0;
 
    return StatusCode::SUCCESS;
}
 
 
StatusCode RecCgemDedxHitCnv::DataObjectToTObject(DataObject* obj,RootAddress* rootaddr) {
  MsgStream log(msgSvc(), "RecCgemDedxHitCnv");
  log << MSG::DEBUG << "RecCgemDedxHitCnv::DataObjectToTObject" << endreq;
  StatusCode sc=StatusCode::SUCCESS;
  RecCgemDedxHitCol * recCgemDedxHitCol=dynamic_cast<RecCgemDedxHitCol *> (obj);
  if (!recCgemDedxHitCol) {
    log << MSG::ERROR << "Could not downcast to RecCgemDedxHitCol" << endreq;
    return StatusCode::FAILURE;
  }
  DataObject *evt;
  m_eds->findObject(EventModel::Recon::Event,evt);
  if (evt==NULL) {
    log << MSG::ERROR << "Could not get RecEvent in TDS "  << endreq;
    return StatusCode::FAILURE;
  }
  ReconEvent * devtTds=dynamic_cast<ReconEvent *> (evt);
  if (!devtTds) {
    log << MSG::ERROR << "RecCgemDedxHitCnv:Could not downcast to TDS RecEvent" << endreq;
  }
  IOpaqueAddress *addr;
  m_cnvSvc->getRecTrackCnv()->createRep(evt,addr); 
  TRecTrackEvent *recEvt=m_cnvSvc->getRecTrackCnv()->getWriteObject();
  const TObjArray *m_recCgemDedxHitCol = recEvt->getRecCgemDedxHitCol();
  if (!m_recCgemDedxHitCol) return sc;
  recEvt->clearRecCgemDedxHitCol(); //necessary in case there is I/O at the same time since array is static
  RecCgemDedxHitCol::const_iterator recCgemDedxHit;
  for (recCgemDedxHit = recCgemDedxHitCol->begin(); recCgemDedxHit != recCgemDedxHitCol->end(); recCgemDedxHit++) {
    int trackId = (*recCgemDedxHit)->getTrkId();
//     UInt_t  cgemIdMSB =(*recCgemDedxHit)->getCgemId().get_valueMSB();
//     UInt_t  cgemIdLSB =(*recCgemDedxHit)->getCgemId().get_valueLSB();
    UInt_t  cgemId =(*recCgemDedxHit)->getCgemId().get_value();
    double pathlength = (*recCgemDedxHit)->getPathLength();
     log << MSG::DEBUG<<"check Reconstrunction of dE/dx Hit root"<<" pathlength is : "<<pathlength<<endreq;
     TRecCgemDedxHit *recCgemDedxHitRoot = new TRecCgemDedxHit();
     //m_common.m_recCgemDedxHitMap[(*recCgemDedxHit)] = recCgemDedxHitRoot;
     recCgemDedxHitRoot->setTrkId(trackId);
//      recCgemDedxHitRoot->setCgemIdMSB( cgemIdMSB );
//      recCgemDedxHitRoot->setCgemIdLSB( cgemIdLSB );
     recCgemDedxHitRoot->setCgemId( cgemId );
     
     recCgemDedxHitRoot->setPathLength( pathlength);
     log << MSG::DEBUG<<"check Reconstrunction of dE/dx Hit root"<<" pathlength 2 is : "<<(*recCgemDedxHit)->getPathLength()<<endreq;
     recEvt->addRecCgemDedxHit(recCgemDedxHitRoot);
  }
 
  return StatusCode::SUCCESS;
}
 
#endif