/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Event/RootCnvSvc/RootCnvSvc-02-01-12/src/Mc/MdcMcHitCnv.cxx

Go to the documentation of this file.

#ifndef MdcMcHitCnv_CXX
#define MdcMcHitCnv_CXX  1

#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/ObjectVector.h"

#include "TClonesArray.h"

#include "EventModel/EventModel.h"
#include "McTruth/McEvent.h"   //TDS object
#include "McTruth/MdcMcHit.h"   //TDS object
#include "RootEventData/TMdcMc.h" // standard root object
#include "RootEventData/TMcEvent.h"
#include "RootCnvSvc/Mc/McCnv.h"
#include "RootCnvSvc/Mc/MdcMcHitCnv.h"
#include "RootCnvSvc/RootAddress.h"


// Instantiation of a static factory class used by clients to create
// instances of this service
//static CnvFactory<MdcMcHitCnv> s_factory;
//const ICnvFactory& MdcMcHitCnvFactory = s_factory;

MdcMcHitCnv::MdcMcHitCnv(ISvcLocator* svc)
: RootEventBaseCnv(classID(), svc)
{
    // Here we associate this converter with the /Event path on the TDS.
    MsgStream log(msgSvc(), "MdcMcHitCnv");
    //log << MSG::DEBUG << "Constructor called for " << objType() << endreq;
    m_rootBranchname ="m_mdcMcHitCol";
    //declareObject(EventModel::MC::MdcMcHitCol, objType(), m_rootTreename, m_rootBranchname);
    m_adresses.push_back(&m_mdcMcHitCol);
    m_mdcMcHitCol=0;
}

StatusCode MdcMcHitCnv::TObjectToDataObject(DataObject*& refpObject) {
  // creation of TDS object from root object

    MsgStream log(msgSvc(), "MdcMcHitCnv");
    log << MSG::DEBUG << "MdcMcHitCnv::TObjectToDataObject" << endreq;
    StatusCode sc=StatusCode::SUCCESS;

    // create the TDS location for the MdcMc Collection
    MdcMcHitCol* mdcMcTdsCol = new MdcMcHitCol;
    refpObject=mdcMcTdsCol;


    // now convert
    if (!m_mdcMcHitCol) return sc;
    TIter mdcMcIter(m_mdcMcHitCol);
    TMdcMc *mdcMcRoot = 0;
    while ((mdcMcRoot = (TMdcMc*)mdcMcIter.Next())) {
        unsigned int id = mdcMcRoot ->getId();
                unsigned int trackIndex = mdcMcRoot ->getTrackIndex();
                double xPosition =   mdcMcRoot ->getPositionX() ;
                double yPosition =   mdcMcRoot ->getPositionY() ;
                double zPosition =   mdcMcRoot ->getPositionZ() ;
                double  driftDistance= mdcMcRoot ->getDriftDistance() ;
                double depositEnergy = mdcMcRoot ->getDepositEnergy() ;
          int positionFlag = mdcMcRoot->getPositionFlag();
      
                MdcMcHit *mdcMcTds = new MdcMcHit;
                m_common.m_rootMdcMcHitMap[mdcMcRoot] = mdcMcTds;
        
                mdcMcTds->setIdentifier(id);
                mdcMcTds->setTrackIndex(trackIndex);
                mdcMcTds->setPositionX(xPosition);
                mdcMcTds->setPositionY(yPosition);
                mdcMcTds->setPositionZ(zPosition);
                mdcMcTds->setDriftDistance(driftDistance);
                mdcMcTds->setDepositEnergy(depositEnergy);
    mdcMcTds->setPositionFlag(positionFlag);
    
        mdcMcTdsCol->push_back(mdcMcTds);
     }
    //m_mdcMcHitCol->Delete();  // wensp add 2005/12/30
    delete m_mdcMcHitCol;
    m_mdcMcHitCol = 0;
   return StatusCode::SUCCESS;
}

StatusCode MdcMcHitCnv::DataObjectToTObject(DataObject* obj,RootAddress* rootaddr) {

  MsgStream log(msgSvc(), "MdcMcHitCnv");
  log << MSG::DEBUG << "MdcMcHitCnv::DataObjectToTObject" << endreq;
  StatusCode sc=StatusCode::SUCCESS;
 
 MdcMcHitCol * mdcMcHitCnvTds=dynamic_cast<MdcMcHitCol *> (obj);
  if (!mdcMcHitCnvTds) {
    log << MSG::ERROR << "Could not downcast to MdcMcHitCol" << endreq;
    return StatusCode::FAILURE;
  }
 
  DataObject *evt;
  m_eds->findObject(EventModel::MC::Event,evt);
  if (evt==NULL) {
    log << MSG::ERROR << "Could not get McEvent in TDS "  << endreq;
    return StatusCode::FAILURE;
  }
  McEvent * devtTds=dynamic_cast<McEvent *> (evt);
  if (!devtTds) {
    log << MSG::ERROR << "MdcMcHitCnv:Could not downcast to TDS McEvent" << endreq;
  }
  IOpaqueAddress *addr;

  m_cnvSvc->getMcCnv()->createRep(evt,addr); 
  TMcEvent *McEvt=m_cnvSvc->getMcCnv()->getWriteObject();

  const TObjArray *m_mdcMcHitCol = McEvt->getMdcMcHitCol();
  if (!m_mdcMcHitCol) return sc;
  McEvt->clearMdcMcHitCol(); //necessary in case there is I/O at the same time since array is static
  MdcMcHitCol::const_iterator mdcMcTds;

  for (mdcMcTds = mdcMcHitCnvTds->begin(); mdcMcTds != mdcMcHitCnvTds->end(); mdcMcTds++) {
                UInt_t id  = ((*mdcMcTds)->identify()).get_value() ;
                UInt_t trackIndex = (*mdcMcTds) ->getTrackIndex();
                Double_t xPosition =   (*mdcMcTds) ->getPositionX() ;
                Double_t yPosition =   (*mdcMcTds) ->getPositionY() ;
                Double_t zPosition =   (*mdcMcTds) ->getPositionZ() ;
                Double_t driftDistance=   (*mdcMcTds) ->getDriftDistance();
                Double_t depositEnergy = (*mdcMcTds) ->getDepositEnergy() ;
    Int_t positionFlag = (*mdcMcTds) ->getPositionFlag();
    
                TMdcMc *mdcMcRoot = new TMdcMc();
                //m_common.m_mdcMcHitMap[(*mdcMcTds)] = mdcMcRoot;

                mdcMcRoot->setId(id);
                mdcMcRoot->setTrackIndex(trackIndex);
                mdcMcRoot->setPositionX(xPosition);
                mdcMcRoot->setPositionY(yPosition);
                mdcMcRoot->setPositionZ(zPosition);
                mdcMcRoot->setDriftDistance(driftDistance);
                mdcMcRoot->setDepositEnergy(depositEnergy);
    mdcMcRoot->setPositionFlag(positionFlag);

                McEvt->addMdcMc(mdcMcRoot);
  }

  return StatusCode::SUCCESS;
}
#endif

---