RecMdcHitCnv.cxx

Go to the documentation of this file.

#ifndef RecMdcHitCnv_CXX
#define RecMdcHitCnv_CXX 1
 
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/ObjectVector.h"
 
#include "TClonesArray.h"
 
#include "EventModel/EventModel.h"
 
#include "Identifier/Identifier.h"
#include "MdcRecEvent/RecMdcHit.h"
#include "ReconEvent/ReconEvent.h"
 
#include "RootEventData/TRecMdcHit.h" // standard root object
#include "RootEventData/TRecTrackEvent.h"
 
#include "RootCnvSvc/Rec/RecMdcHitCnv.h"
#include "RootCnvSvc/Rec/RecTrackCnv.h"
#include "RootCnvSvc/RootAddress.h"
 
#include <vector>
 
using namespace std;
 
RecMdcHitCnv::RecMdcHitCnv( ISvcLocator* svc ) : RootEventBaseCnv( classID(), svc ) {
 
  // Here we associate this converter with the /Event path on the TDS.
  MsgStream log( msgSvc(), "RecMdcHitCnv" );
  // log << MSG::DEBUG << "Constructor called for " << objType() << endmsg;
  // m_rootTreename ="Rec";
  m_rootBranchname = "m_recMdcHitCol";
  // declareObject(EventModel::Recon::RecMdcHitCol, objType(), m_rootTreename,
  // m_rootBranchname);
  m_adresses.push_back( &m_recMdcHitCol );
  m_recMdcHitCol = 0;
}
 
StatusCode RecMdcHitCnv::TObjectToDataObject( DataObject*& refpObject ) {
  // creation of TDS object from root object
  MsgStream log( msgSvc(), "RecMdcHitCnv" );
  log << MSG::DEBUG << "RecMdcHitCnv::TObjectToDataObject" << endmsg;
 
  // create the TDS location for the RecMdcHit Collection
  RecMdcHitCol* recMdcHitCol = new RecMdcHitCol;
  refpObject                 = recMdcHitCol;
 
  // now convert
  if ( !m_recMdcHitCol ) return StatusCode::SUCCESS;
  TIter       mdcHitIter( m_recMdcHitCol );
  TRecMdcHit* recMdcHitRoot = 0;
  while ( ( recMdcHitRoot = (TRecMdcHit*)mdcHitIter.Next() ) )
  {
    int      id                = recMdcHitRoot->getId();
    int      trackId           = recMdcHitRoot->getTrkId();
    double   driftDistLeft     = recMdcHitRoot->getDriftDistLeft();
    double   driftDistRight    = recMdcHitRoot->getDriftDistRight();
    double   errDriftDistLeft  = recMdcHitRoot->getErrDriftDistLeft();
    double   errDriftDistRight = recMdcHitRoot->getErrDriftDistRight();
    double   chisqAdd          = recMdcHitRoot->getChisqAdd();
    int      flagLR            = recMdcHitRoot->getFlagLR();
    int      stat              = recMdcHitRoot->getStat();
    long int mdcId( recMdcHitRoot->getMdcId() );
    double   tdc    = recMdcHitRoot->getTdc();
    double   adc    = recMdcHitRoot->getAdc();
    double   driftT = recMdcHitRoot->getDriftT();
    double   doca   = recMdcHitRoot->getDoca();
    double   entra  = recMdcHitRoot->getEntra();
    double   zhit   = recMdcHitRoot->getZhit();
    double   fltLen = recMdcHitRoot->getFltLen();
 
    RecMdcHit* recMdcHit                       = new RecMdcHit();
    m_common.m_rootRecMdcHitMap[recMdcHitRoot] = recMdcHit;
 
    recMdcHit->setId( id );
    recMdcHit->setTrkId( trackId );
    recMdcHit->setDriftDistLeft( driftDistLeft );
    recMdcHit->setDriftDistRight( driftDistRight );
    recMdcHit->setErrDriftDistLeft( errDriftDistLeft );
    recMdcHit->setErrDriftDistRight( errDriftDistRight );
    recMdcHit->setChisqAdd( chisqAdd );
    recMdcHit->setFlagLR( flagLR );
    recMdcHit->setStat( stat );
    recMdcHit->setMdcId( Identifier( mdcId ) );
    recMdcHit->setTdc( tdc );
    recMdcHit->setAdc( adc );
    recMdcHit->setDriftT( driftT );
    recMdcHit->setDoca( doca );
    recMdcHit->setEntra( entra );
    recMdcHit->setZhit( zhit );
    recMdcHit->setFltLen( fltLen );
 
    recMdcHitCol->push_back( recMdcHit );
  }
  return StatusCode::SUCCESS;
}
 
StatusCode RecMdcHitCnv::DataObjectToTObject( DataObject* obj, RootAddress* rootaddr ) {
 
  MsgStream log( msgSvc(), "RecMdcHitCnv" );
  log << MSG::DEBUG << "RecMdcHitCnv::DataObjectToTObject" << endmsg;
 
  RecMdcHitCol* mdcHitColTds = dynamic_cast<RecMdcHitCol*>( obj );
  if ( !mdcHitColTds )
  {
    log << MSG::ERROR << "Could not downcast to RecMdcHitCol" << endmsg;
    return StatusCode::FAILURE;
  }
 
  DataObject* evt;
  m_eds->findObject( EventModel::Recon::Event, evt );
  if ( evt == NULL )
  {
    log << MSG::ERROR << "Could not get ReconEvent in TDS " << endmsg;
    return StatusCode::FAILURE;
  }
 
  ReconEvent* devtTds = dynamic_cast<ReconEvent*>( evt );
  if ( !devtTds )
  { log << MSG::ERROR << "RecMdcHitCnv:Could not downcast to TDS ReconEvent" << endmsg; }
 
  IOpaqueAddress* addr;
 
  m_cnvSvc->getRecTrackCnv()->createRep( evt, addr );
  TRecTrackEvent* recEvt = m_cnvSvc->getRecTrackCnv()->getWriteObject();
 
  const TObjArray* m_recMdcHitCol = recEvt->getRecMdcHitCol();
  if ( !m_recMdcHitCol ) return StatusCode::SUCCESS;
  recEvt->clearRecMdcHitCol(); // necessary in case there is I/O at the same time since array
                               // is static
 
  RecMdcHitCol::const_iterator recMdcHit;
  for ( recMdcHit = mdcHitColTds->begin(); recMdcHit != mdcHitColTds->end(); recMdcHit++ )
  {
    Bool_t   isGrouped         = ( *recMdcHit )->isGrouped();
    Int_t    id                = ( *recMdcHit )->getId();
    Int_t    trackId           = ( *recMdcHit )->getTrkId();
    Double_t driftDistLeft     = ( *recMdcHit )->getDriftDistLeft();
    Double_t driftDistRight    = ( *recMdcHit )->getDriftDistRight();
    Double_t errDriftDistLeft  = ( *recMdcHit )->getErrDriftDistLeft();
    Double_t errDriftDistRight = ( *recMdcHit )->getErrDriftDistRight();
    Double_t chisqAdd          = ( *recMdcHit )->getChisqAdd();
    Int_t    flagLR            = ( *recMdcHit )->getFlagLR();
    Int_t    stat              = ( *recMdcHit )->getStat();
    UInt_t   mdcId             = ( *recMdcHit )->getMdcId().get_value();
    Double_t tdc               = ( *recMdcHit )->getTdc();
    Double_t adc               = ( *recMdcHit )->getAdc();
    Double_t driftT            = ( *recMdcHit )->getDriftT();
    Double_t doca              = ( *recMdcHit )->getDoca();
    Double_t entra             = ( *recMdcHit )->getEntra();
    Double_t zhit              = ( *recMdcHit )->getZhit();
    Double_t fltLen            = ( *recMdcHit )->getFltLen();
 
    TRecMdcHit* recMdcHitRoot = new TRecMdcHit();
    // m_common.m_recMdcHitMap[(*recMdcHit)] = recMdcHitRoot;
 
    recMdcHitRoot->setIsGrouped( isGrouped );
    recMdcHitRoot->setId( id );
    recMdcHitRoot->setTrkId( trackId );
    recMdcHitRoot->setDriftDistLeft( driftDistLeft );
    recMdcHitRoot->setDriftDistRight( driftDistRight );
    recMdcHitRoot->setErrDriftDistLeft( errDriftDistLeft );
    recMdcHitRoot->setErrDriftDistRight( errDriftDistRight );
    log << MSG::INFO << " recMdcHitRoot.ddl: " << recMdcHitRoot->getDriftDistLeft()
        << " recMdcHitRoot.erddl: " << recMdcHitRoot->getErrDriftDistLeft() << endmsg;
    recMdcHitRoot->setChisqAdd( chisqAdd );
    recMdcHitRoot->setFlagLR( flagLR );
    recMdcHitRoot->setStat( stat );
    recMdcHitRoot->setMdcId( mdcId );
    recMdcHitRoot->setTdc( tdc );
    recMdcHitRoot->setAdc( adc );
    recMdcHitRoot->setDriftT( driftT );
    recMdcHitRoot->setDoca( doca );
    recMdcHitRoot->setEntra( entra );
    recMdcHitRoot->setZhit( zhit );
    recMdcHitRoot->setFltLen( fltLen );
 
    recEvt->addRecMdcHit( recMdcHitRoot );
  }
  return StatusCode::SUCCESS;
}
#endif