RecCgemKalHelixSegCnv.cxx

Go to the documentation of this file.

#ifndef RecCgemKalHelixSegCnv_CXX
#define RecCgemKalHelixSegCnv_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/RecCgemKalHelixSeg.h"
#include "Identifier/Identifier.h"
 
#include "RootEventData/TRecCgemKalHelixSeg.h" // standard root object
#include "RootEventData/TRecTrackEvent.h"
 
#include "RootCnvSvc/Rec/RecTrackCnv.h"
#include "RootCnvSvc/Rec/RecCgemKalHelixSegCnv.h"
#include "RootCnvSvc/RootAddress.h"
 
#include <vector>
 
using namespace std;
 
 
RecCgemKalHelixSegCnv::RecCgemKalHelixSegCnv(ISvcLocator* svc)
: RootEventBaseCnv(classID(), svc)
{
 
    // Here we associate this converter with the /Event path on the TDS.
    MsgStream log(msgSvc(), "RecCgemKalHelixSegCnv");
    //log << MSG::DEBUG << "Constructor called for " << objType() << endreq;
    //m_rootTreename ="Rec";
    m_rootBranchname ="m_recCgemKalHelixSegCol";
    //declareObject(EventModel::Recon::RecCgemKalHelixSegCol, objType(), m_rootTreename, m_rootBranchname);
    m_adresses.push_back(&m_recCgemKalHelixSegCol);
    m_recCgemKalHelixSegCol=0;
}
 
StatusCode RecCgemKalHelixSegCnv::TObjectToDataObject(DataObject*& refpObject) {
  // creation of TDS object from root object
  MsgStream log(msgSvc(), "RecCgemKalHelixSegCnv");
  log << MSG::DEBUG << "RecCgemKalHelixSegCnv::TObjectToDataObject" << endreq;
 
  // create the TDS location for the RecCgemKalHelixSeg Collection
  RecCgemKalHelixSegCol* recCgemKalHelixSegCol = new RecCgemKalHelixSegCol;
  refpObject = recCgemKalHelixSegCol;
 
  // now convert
  if (!m_recCgemKalHelixSegCol) return StatusCode::SUCCESS;
  TIter cgemHitIter(m_recCgemKalHelixSegCol);
  TRecCgemKalHelixSeg *recCgemKalHelixSegRoot = 0;
  while ((recCgemKalHelixSegRoot = (TRecCgemKalHelixSeg*)cgemHitIter.Next())) {
    int      trackId = recCgemKalHelixSegRoot->getTrackId();
//     unsigned int cgemIdMSB = recCgemKalHelixSegRoot->getCgemIdMSB();
//     unsigned int cgemIdLSB = recCgemKalHelixSegRoot->getCgemIdLSB();
//     Identifier64 cgemId(cgemIdMSB,cgemIdLSB); 
    Identifier cgemId(recCgemKalHelixSegRoot->getCgemId());
    double   tdc = recCgemKalHelixSegRoot->getTdc();
    double   adc = recCgemKalHelixSegRoot->getAdc();
    double   zhit = recCgemKalHelixSegRoot->getZhit();
    double   tof = recCgemKalHelixSegRoot->getTof();
    double   docaIncl = recCgemKalHelixSegRoot->getDocaIncl();
    double   docaExcl = recCgemKalHelixSegRoot->getDocaExcl();
    double   entra = recCgemKalHelixSegRoot->getEntra();
  //  int      layerId = recCgemKalHelixSegRoot->getLayerId();
  //  double   drIncl = recCgemKalHelixSegRoot->getDrIncl();
  //  double   phi0Incl = recCgemKalHelixSegRoot->getFi0Incl();
  //  double   kappaIncl = recCgemKalHelixSegRoot->getCpaIncl();
  //  double   dzIncl = recCgemKalHelixSegRoot->getDzIncl();
  //  double   tanlIncl = recCgemKalHelixSegRoot->getTanlIncl();
  //  double   drExcl = recCgemKalHelixSegRoot->getDrExcl();
  //  double   phi0Excl = recCgemKalHelixSegRoot->getFi0Excl();
  //  double   kappaExcl = recCgemKalHelixSegRoot->getCpaExcl();
  //  double   dzExcl = recCgemKalHelixSegRoot->getDzExcl();
  //  double   tanlExcl = recCgemKalHelixSegRoot->getTanlExcl();
    double helixIncl[5]; // errIncl[15]; 
    double helixExcl[5]; // errExcl[15];
    for(int i=0, k=0; i<5; i++){
        helixIncl[i] = recCgemKalHelixSegRoot->getHelixIncl(i);
        helixExcl[i] = recCgemKalHelixSegRoot->getHelixExcl(i);
       // for(int j=0; j<=i; j++){
       //         errIncl[k] = recCgemKalHelixSegRoot->getErrorIncl(i,j);
       //         errExcl[k] = recCgemKalHelixSegRoot->getErrorExcl(i,j);
       // }
    }
 
    RecCgemKalHelixSeg *recCgemKalHelixSeg= new RecCgemKalHelixSeg();
    m_common.m_rootRecCgemKalHelixSegMap[recCgemKalHelixSegRoot] = recCgemKalHelixSeg;
 
    recCgemKalHelixSeg->setTrackId( trackId );
    recCgemKalHelixSeg->setCgemId( Identifier(cgemId) );
    recCgemKalHelixSeg->setTdc( tdc );
    recCgemKalHelixSeg->setAdc( adc );
    recCgemKalHelixSeg->setZhit( zhit );
    recCgemKalHelixSeg->setTof( tof );
    recCgemKalHelixSeg->setDocaIncl( docaIncl );
    recCgemKalHelixSeg->setDocaExcl( docaExcl );
    recCgemKalHelixSeg->setEntra( entra );
//    recCgemKalHelixSeg->setLayerId( layerId );
//    recCgemKalHelixSeg->setDrIncl( drIncl );
//    recCgemKalHelixSeg->setFi0Incl( phi0Incl );
//    recCgemKalHelixSeg->setCpaIncl( kappaIncl );
//    recCgemKalHelixSeg->setDzIncl( dzIncl );
//    recCgemKalHelixSeg->setTanlIncl( tanlIncl );
//    recCgemKalHelixSeg->setDrExcl( drIncl );
//    recCgemKalHelixSeg->setFi0Excl( phi0Excl );
//    recCgemKalHelixSeg->setCpaExcl( kappaExcl );
//    recCgemKalHelixSeg->setDzExcl( dzExcl );
//    recCgemKalHelixSeg->setTanlExcl( tanlExcl );
    recCgemKalHelixSeg->setHelixIncl( helixIncl );
    recCgemKalHelixSeg->setHelixExcl( helixExcl );
//    recCgemKalHelixSeg->setErrorIncl( errIncl );
//    recCgemKalHelixSeg->setErrorExcl( errExcl );
 
    recCgemKalHelixSegCol->push_back(recCgemKalHelixSeg);
  }
    return StatusCode::SUCCESS;
}
 
StatusCode RecCgemKalHelixSegCnv::DataObjectToTObject(DataObject* obj,RootAddress* rootaddr) {
 
  MsgStream log(msgSvc(), "RecCgemKalHelixSegCnv");
  log << MSG::DEBUG << "RecCgemKalHelixSegCnv::DataObjectToTObject" << endreq;
 
  RecCgemKalHelixSegCol * cgemHitColTds=dynamic_cast<RecCgemKalHelixSegCol *> (obj);
  if (!cgemHitColTds) {
    log << MSG::ERROR << "Could not downcast to RecCgemKalHelixSegCol" << 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 << "RecCgemKalHelixSegCnv:Could not downcast to TDS ReconEvent" << endreq;
  }
 
  IOpaqueAddress *addr;
 
  m_cnvSvc->getRecTrackCnv()->createRep(evt,addr); 
  TRecTrackEvent *recEvt=m_cnvSvc->getRecTrackCnv()->getWriteObject();
 
  const TObjArray *m_recCgemKalHelixSegCol = recEvt->getRecCgemKalHelixSegCol();
  if (!m_recCgemKalHelixSegCol) return StatusCode::SUCCESS;
  recEvt->clearRecCgemKalHelixSegCol(); //necessary in case there is I/O at the same time since array is static
 
  RecCgemKalHelixSegCol::const_iterator recCgemKalHelixSeg;
  for (recCgemKalHelixSeg = cgemHitColTds->begin(); recCgemKalHelixSeg != cgemHitColTds->end(); recCgemKalHelixSeg++) {
     Int_t      trackId =(*recCgemKalHelixSeg)->getTrackId();
//      UInt_t     cgemIdMSB = (*recCgemKalHelixSeg)->getCgemId().get_valueMSB();
//      UInt_t     cgemIdLSB = (*recCgemKalHelixSeg)->getCgemId().get_valueLSB();
     UInt_t     cgemId = (*recCgemKalHelixSeg)->getCgemId().get_value();
     Double_t   tdc = (*recCgemKalHelixSeg)->getTdc();
     Double_t   adc = (*recCgemKalHelixSeg)->getAdc();
     Double_t   zhit = (*recCgemKalHelixSeg)->getZhit();
     Double_t   tof = (*recCgemKalHelixSeg)->getTof();
     Double_t   docaIncl = (*recCgemKalHelixSeg)->getDocaIncl();
     Double_t   docaExcl = (*recCgemKalHelixSeg)->getDocaExcl();
     Double_t   entra = (*recCgemKalHelixSeg)->getEntra();
//     Int_t      layerId = (*recCgemKalHelixSeg)->getLayerId();
//     Double_t   drIncl =  (*recCgemKalHelixSeg)->getDrIncl();
//     Double_t   phi0Incl =(*recCgemKalHelixSeg)->getFi0Incl();
//     Double_t   kappaIncl = (*recCgemKalHelixSeg)->getCpaIncl();
//     Double_t   dzIncl =  (*recCgemKalHelixSeg)->getDzIncl();
//     Double_t   tanlIncl = (*recCgemKalHelixSeg)->getTanlIncl();
//     Double_t   drExcl =  (*recCgemKalHelixSeg)->getDrExcl();
//     Double_t   phi0Excl = (*recCgemKalHelixSeg)->getFi0Excl();
//     Double_t   kappaExcl = (*recCgemKalHelixSeg)->getCpaExcl();
//     Double_t   dzExcl =   (*recCgemKalHelixSeg)->getDzExcl();
//     Double_t   tanlExcl = (*recCgemKalHelixSeg)->getTanlExcl();
     HepVector  helixIncl_vct = (*recCgemKalHelixSeg)->getHelixIncl();
     HepVector  helixExcl_vct = (*recCgemKalHelixSeg)->getHelixExcl();
//     HepSymMatrix  errIncl_mtx = (*recCgemKalHelixSeg)->getErrorIncl();
//     HepSymMatrix  errExcl_mtx = (*recCgemKalHelixSeg)->getErrorExcl();
    
    Double_t helixIncl[5], helixExcl[5]; //errIncl[5][5], errExcl[5][5];
     for(Int_t i=0; i<5; i++){
             helixIncl[i] = helixIncl_vct[i];
             helixExcl[i] = helixExcl_vct[i];
     }
//     for (Int_t i=0; i<5; i++){
//             for(Int_t j=0; j<5; j++){
//               errIncl[i][j] = errIncl_mtx[i][j];
//               errExcl[i][j] = errExcl_mtx[i][j];
//             }
//     }
 
     TRecCgemKalHelixSeg *recCgemKalHelixSegRoot = new TRecCgemKalHelixSeg();
     //m_common.m_recCgemKalHelixSegMap[(*recCgemKalHelixSeg)] = recCgemKalHelixSegRoot;
 
     recCgemKalHelixSegRoot->setTrackId( trackId );
//      recCgemKalHelixSegRoot->setCgemIdMSB( cgemIdMSB );
//      recCgemKalHelixSegRoot->setCgemIdLSB( cgemIdLSB );
     recCgemKalHelixSegRoot->setCgemId( cgemId );
     recCgemKalHelixSegRoot->setTdc( tdc );
     recCgemKalHelixSegRoot->setAdc( adc );
     recCgemKalHelixSegRoot->setZhit( zhit );
     recCgemKalHelixSegRoot->setTof( tof );
     recCgemKalHelixSegRoot->setDocaIncl( docaIncl );
     recCgemKalHelixSegRoot->setDocaExcl( docaExcl );
     recCgemKalHelixSegRoot->setEntra( entra );
//     recCgemKalHelixSegRoot->setLayerId( layerId );
//     recCgemKalHelixSegRoot->setDrIncl( drIncl );
//     recCgemKalHelixSegRoot->setFi0Incl( phi0Incl );
//     recCgemKalHelixSegRoot->setCpaIncl( kappaIncl );
//     recCgemKalHelixSegRoot->setDzIncl( dzIncl );
//     recCgemKalHelixSegRoot->setTanlIncl( tanlIncl );
//     recCgemKalHelixSegRoot->setDrExcl( drIncl );
//     recCgemKalHelixSegRoot->setFi0Excl( phi0Excl );
//     recCgemKalHelixSegRoot->setCpaExcl( kappaExcl );
//     recCgemKalHelixSegRoot->setDzExcl( dzExcl );
//     recCgemKalHelixSegRoot->setTanlExcl( tanlExcl );
     recCgemKalHelixSegRoot->setHelixIncl( helixIncl );
     recCgemKalHelixSegRoot->setHelixExcl( helixExcl );
//     recCgemKalHelixSegRoot->setErrorIncl( errIncl );
//     recCgemKalHelixSegRoot->setErrorExcl( errExcl );
 
     recEvt->addRecCgemKalHelixSeg(recCgemKalHelixSegRoot);
  }
  return StatusCode::SUCCESS;
}
#endif