RecCgemSegmentCnv.cxx

Go to the documentation of this file.

#ifndef RECCGEMSEGMENTCNV_CXX
#define RECCGEMSEGMENTCNV_CXX
 
#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/RecCgemSegment.h"
 
#include "RootEventData/TRecCgemSegment.h"
#include "RootEventData/TRecTrackEvent.h"
 
#include "RootCnvSvc/Rec/RecTrackCnv.h"
#include "RootCnvSvc/Rec/RecCgemSegmentCnv.h"
#include "RootCnvSvc/RootAddress.h"
 
#include <vector>
 
using namespace std;
 
RecCgemSegmentCnv::RecCgemSegmentCnv(ISvcLocator* svc)
    :RootEventBaseCnv(classID(), svc)
{
 
        MsgStream log(msgSvc(),"RecCgemSegmentCnv");
        log << MSG::INFO << " RecCgemSegmentCnv ------ " << endreq;
        m_rootBranchname = "m_recCgemSegmentCol";
        m_adresses.push_back(&m_recCgemSegmentCol);
        m_recCgemSegmentCol = 0;
}
StatusCode RecCgemSegmentCnv::TObjectToDataObject(DataObject*& refpObject)
{
    MsgStream log(msgSvc(),"RecCgemSegmentCnv");
    log << MSG::DEBUG << "RecCgemSegmentCnv::TObjectToDataObject" << endreq;
 
 
    RecCgemSegmentCol* recCgemSegmentCol = new RecCgemSegmentCol;
    refpObject = recCgemSegmentCol;
 
 
    //convert//
    if (!m_recCgemSegmentCol) return StatusCode::SUCCESS;
    TIter cgemSegmentIter(m_recCgemSegmentCol);
    TRecCgemSegment *recCgemSegmentRoot = 0;
    while ((recCgemSegmentRoot = (TRecCgemSegment*)cgemSegmentIter.Next()))
    {
        int clusterid_1 = recCgemSegmentRoot->getclusterid_1();
        int clusterid_2 = recCgemSegmentRoot->getclusterid_2();
        int clusterid_3 = recCgemSegmentRoot->getclusterid_3();
        int eventid = recCgemSegmentRoot->geteventid();
        int match = recCgemSegmentRoot->getmatch();
        int segmentid = recCgemSegmentRoot->getsegmentid();
        double  helix[5];
        double  helix_err[15];
        for(int i = 0; i < 5; i++){
           helix[i]     = (recCgemSegmentRoot)->gethelix(i);
        }
        for(int i = 0; i < 15; i++){
           helix_err[i] = (recCgemSegmentRoot)->gethelix_err(i);
        }
 
        RecCgemSegment* recCgemSegment = new RecCgemSegment();
        m_common.m_rootRecCgemSegmentMap[recCgemSegmentRoot] = recCgemSegment;
 
        recCgemSegment->setclusterid_1(clusterid_1);
        recCgemSegment->setclusterid_2(clusterid_2);
        recCgemSegment->setclusterid_3(clusterid_3);
        recCgemSegment->seteventid(eventid);
        recCgemSegment->setmatch(match);
        recCgemSegment->setsegmentid(segmentid);
        recCgemSegment->sethelix( helix );
        recCgemSegment->sethelix_err( helix_err );
 
        recCgemSegmentCol->push_back(recCgemSegment);
    }
    return StatusCode::SUCCESS;
}
StatusCode RecCgemSegmentCnv::DataObjectToTObject(DataObject* obj,RootAddress* rootaddr)
{
    MsgStream log(msgSvc(), "RecCgemSegmentCnv");
    log << MSG::DEBUG << "RecCgemSegmentCnv::DataObjectToTObject" << endreq;
 
    RecCgemSegmentCol * cgemSegmentColTds=dynamic_cast<RecCgemSegmentCol *> (obj);
    if (!cgemSegmentColTds) {
        log << MSG::ERROR << "Could not downcast to RecCgemSegmentCol" << 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();
//  m_cnvSvc->getRecSegmentCnv()->createRep(evt,addr);
//  TRecTrackEvent *recEvt=m_cnvSvc->getRecSegmentCnv()->getWriteObject();
 
    const TObjArray *m_recCgemSegmentCol = recEvt->getRecCgemSegmentCol();
    if (!m_recCgemSegmentCol) return StatusCode::SUCCESS;
    recEvt->clearRecCgemSegmentCol();
 
    RecCgemSegmentCol::const_iterator recCgemSegment;
    for (recCgemSegment = cgemSegmentColTds->begin(); recCgemSegment != cgemSegmentColTds->end(); recCgemSegment++)
    {
        Int_t eventid = (*recCgemSegment)->geteventid();
        Int_t match = (*recCgemSegment)->getmatch();
        Int_t segmentid = (*recCgemSegment)->getsegmentid();
        Int_t clusterid_1 = (*recCgemSegment)->getclusterid_1();
        Int_t clusterid_2 = (*recCgemSegment)->getclusterid_2();
        Int_t clusterid_3 = (*recCgemSegment)->getclusterid_3();
        Double_t  helix[5];
        Double_t  helix_err[15];
        for(int i = 0; i < 5; i++){
           helix[i]     = (*recCgemSegment)->gethelix(i);
        }
        for(int i = 0; i < 15; i++){
           helix_err[i] = (*recCgemSegment)->gethelix_err(i);
        }
 
        TRecCgemSegment *recCgemSegmentRoot = new TRecCgemSegment();
 
        recCgemSegmentRoot->seteventid(eventid);
        recCgemSegmentRoot->setmatch(match);
        recCgemSegmentRoot->setsegmentid(segmentid);
        recCgemSegmentRoot->setclusterid_1(clusterid_1);
        recCgemSegmentRoot->setclusterid_2(clusterid_2);
        recCgemSegmentRoot->setclusterid_3(clusterid_3);
        recCgemSegmentRoot->sethelix( helix );
        recCgemSegmentRoot->sethelix_err( helix_err );
 
        recEvt->addRecCgemSegment(recCgemSegmentRoot);
    }
    return StatusCode::SUCCESS;
}
#endif