/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Event/RootCnvSvc/RootCnvSvc-02-01-12/src/Rec/RecMdcHitCnv.cxx

Go to the documentation of this file.

#ifndef RecMdcHitCnv_CXX
#define RecMdcHitCnv_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 "MdcRecEvent/RecMdcHit.h"
#include "Identifier/Identifier.h"

#include "RootEventData/TRecMdcHit.h" // standard root object
#include "RootEventData/TRecTrackEvent.h"

#include "RootCnvSvc/Rec/RecTrackCnv.h"
#include "RootCnvSvc/Rec/RecMdcHitCnv.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() << endreq;
    //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" << endreq;

  // 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" << endreq;
 
  RecMdcHitCol * mdcHitColTds=dynamic_cast<RecMdcHitCol *> (obj);
  if (!mdcHitColTds) {
    log << MSG::ERROR << "Could not downcast to RecMdcHitCol" << 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 << "RecMdcHitCnv:Could not downcast to TDS ReconEvent" << endreq;
  }

  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()
                   <<endreq;
     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

---