CgemDedxCnv.cxx

Go to the documentation of this file.

#ifndef DedxCnv_CXX
#define DedxCnv_CXX  1
 
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/ObjectVector.h"
 
#include "TClonesArray.h"
 
#include "EventModel/EventModel.h"
#include "DstEvent/DstEvent.h"   //TDS object
#include "DstEvent/DstCgemDedx.h"   //TDS object
#include "RootEventData/TCgemDedx.h" // standard root object
#include "RootEventData/TDstEvent.h"
 
#include "RootCnvSvc/Dst/DstCnv.h"
#include "RootCnvSvc/Dst/CgemDedxCnv.h"
#include "RootCnvSvc/RootAddress.h"
 
 
// Instantiation of a static factory class used by clients to create
// instances of this service
//static CnvFactory<DedxCnv> s_factory;
//const ICnvFactory& DedxCnvFactory = s_factory;
 
CgemDedxCnv::CgemDedxCnv(ISvcLocator* svc)
: RootEventBaseCnv(classID(), svc)
{
    // Here we associate this converter with the /Event path on the TDS.
    MsgStream log(msgSvc(), "CgemDedxCnv");
    //log << MSG::DEBUG << "Constructor called for " << objType() << endreq;
    m_rootBranchname ="m_cgemDedxCol";
    //declareObject(EventModel::Dst::DstDedxCol, objType(), m_rootTreename, m_rootBranchname);
    m_adresses.push_back(&m_cgemDedxCol);
    m_cgemDedxCol=0;
}
 
StatusCode CgemDedxCnv::TObjectToDataObject(DataObject*& refpObject) {
  // creation of TDS object from root object
 
    MsgStream log(msgSvc(), "CgemDedxCnv");
    log << MSG::DEBUG << "CgemDedxCnv::TObjectToDataObject" << endreq;
    StatusCode sc=StatusCode::SUCCESS;
 
    // create the TDS location for the Dedx Collection
    DstCgemDedxCol* cgemDedxTdsCol = new DstCgemDedxCol;
    refpObject=cgemDedxTdsCol;
 
 
    // now convert
    if (!m_cgemDedxCol) return sc;
    TIter cgemDedxIter(m_cgemDedxCol);
    TCgemDedx *cgemDedxRoot = 0;
    while ((cgemDedxRoot = (TCgemDedx*)cgemDedxIter.Next())) {
     
     int trackId =  cgemDedxRoot->trackId();        
     int particleId = cgemDedxRoot->particleId(); 
     int status = cgemDedxRoot->status();     
     int truncAlg = cgemDedxRoot->truncAlg();
     /*
     double pb[5];
     for ( int i = 0; i < 5; i++) 
          pb[i] = dedxRoot->prob(i);
     
     double numSigmaE  = dedxRoot->numSigmaE();   
     double numSigmaMu = dedxRoot->numSigmaMu();  
     double numSigmaPi = dedxRoot->numSigmaPi();  
     double numSigmaK = dedxRoot->numSigmaK();   
     double numSigmaP = dedxRoot->numSigmaP();   
     */
     double chi[5];
     chi[0] = cgemDedxRoot->chiE();
     chi[1] = cgemDedxRoot->chiMu();
     chi[2] = cgemDedxRoot->chiPi();
     chi[3] = cgemDedxRoot->chiK();
     chi[4] = cgemDedxRoot->chiP();
     int numGoodHits = cgemDedxRoot->numGoodHits(); 
     int numTotalHits = cgemDedxRoot->numTotalHits();
 
     double probPH = cgemDedxRoot->probPH();      
     double normPH = cgemDedxRoot->normPH();
     double errorPH = cgemDedxRoot->errorPH();
     double twentyPH = cgemDedxRoot->twentyPH();        
     //double fracErrPH = dedxRoot-> fracErrPH();  
     //double minIronPH = dedxRoot->minIronPH();   
     //double corrPH = dedxRoot->corrPH();   
    
     DstCgemDedx *cgemDedxTds = new DstCgemDedx();
     m_common.m_rootCgemDedxMap[cgemDedxRoot] = cgemDedxTds;
 
     cgemDedxTds->setTrackId(trackId);
     cgemDedxTds->setParticleId(particleId);
     cgemDedxTds->setStatus (status);
     cgemDedxTds->setTruncAlg(truncAlg);
    /*
     dedxTds->setProb(pb); 
     dedxTds->setNumSigmaE(numSigmaE);       
     dedxTds->setNumSigmaMu(numSigmaMu);     
     dedxTds->setNumSigmaPi(numSigmaPi);     
     dedxTds->setNumSigmaK(numSigmaK);       
     dedxTds->setNumSigmaP(numSigmaP);       
     */
     cgemDedxTds->setChi(chi);
     
     cgemDedxTds->setNumGoodHits( numGoodHits);     
     cgemDedxTds->setNumTotalHits( numTotalHits); 
 
     cgemDedxTds->setProbPH(probPH);             
     cgemDedxTds->setNormPH(normPH);
     cgemDedxTds->setErrorPH(errorPH);
     cgemDedxTds->setTwentyPH(twentyPH);  
     //dedxTds->setFracErrPH(fracErrPH);       
     //dedxTds->setMinIronPH(minIronPH);       
     //dedxTds->setCorrPH(corrPH); 
     
     cgemDedxTdsCol->push_back(cgemDedxTds); 
     //delete dedxTds;
     // dedxTds = NULL;
     }
    //m_dedxCol->Delete();
    delete m_cgemDedxCol;
    m_cgemDedxCol = 0;
   return StatusCode::SUCCESS;
}
 
StatusCode CgemDedxCnv::DataObjectToTObject(DataObject* obj,RootAddress* rootaddr) {
 
  MsgStream log(msgSvc(), "CgemDedxCnv");
  log << MSG::DEBUG << "CgemDedxCnv::DataObjectToTObject" << endreq;
  StatusCode sc=StatusCode::SUCCESS;
 
  DstCgemDedxCol * cgemDedxColTds=dynamic_cast<DstCgemDedxCol *> (obj);
  if (!cgemDedxColTds) {
    log << MSG::ERROR << "Could not downcast to CgemDedxCol" << endreq;
    return StatusCode::FAILURE;
  }
 
  DataObject *evt;
  m_eds->findObject(EventModel::Dst::Event,evt);
  if (evt==NULL) {
    log << MSG::ERROR << "Could not get DstEvent in TDS "  << endreq;
    return StatusCode::FAILURE;
  }
  DstEvent * devtTds=dynamic_cast<DstEvent *> (evt);
  if (!devtTds) {
    log << MSG::ERROR << "CgemDedxCnv:Could not downcast to TDS DstEvent" << endreq;
  }
  IOpaqueAddress *addr;
 
  m_cnvSvc->getDstCnv()->createRep(evt,addr); 
  TDstEvent *recEvt=m_cnvSvc->getDstCnv()->getWriteObject();
 
  const TObjArray *m_cgemDedxCol = recEvt->getCgemDedxCol();
  if (!m_cgemDedxCol) return sc;
  recEvt->clearCgemDedxCol(); //necessary in case there is I/O at the same time since array is static
  DstCgemDedxCol::const_iterator cgemDedxTds;
 
  for (cgemDedxTds = cgemDedxColTds->begin(); cgemDedxTds != cgemDedxColTds->end(); cgemDedxTds++) {
    Int_t trackId = (*cgemDedxTds)->trackId();
    Int_t particleId = (*cgemDedxTds)->particleId(); 
    Int_t status = (*cgemDedxTds)->status();     
    Int_t truncAlg =(*cgemDedxTds)->truncAlg();  
    /*
    Double_t pb[5];
    for (Int_t i = 0; i < 5; i++) 
         pb[i] = (*dedxTds)->prob(i);
    */
    Double_t chiE  = (*cgemDedxTds)->chi(0);   
    Double_t chiMu = (*cgemDedxTds)->chi(1);  
    Double_t chiPi = (*cgemDedxTds)->chi(2);  
    Double_t chiK = (*cgemDedxTds)->chi(3);   
    Double_t chiP = (*cgemDedxTds)->chi(4);   
  
    Int_t numGoodHits = (*cgemDedxTds)->numGoodHits(); 
    Int_t numTotalHits = (*cgemDedxTds)->numTotalHits();
 
    Double_t probPH = (*cgemDedxTds)->probPH();      
    Double_t normPH = (*cgemDedxTds)->normPH();
    Double_t errorPH = (*cgemDedxTds)->errorPH();
    Double_t twentyPH =(*cgemDedxTds)->twentyPH();
    //Double_t fracErrPH = (*dedxTds)-> fracErrPH();  
    //Double_t minIronPH = (*dedxTds)->minIronPH();   
    //Double_t corrPH = (*dedxTds)->corrPH();      
    log << MSG::INFO <<"check Reconstrunction root"<<" particle Id is : "<<particleId
         <<"track id is : "<<trackId
         <<" and status is : "<<status<<endreq;
 
       
    TCgemDedx *cgemDedxRoot = new TCgemDedx();
    //m_common.m_cgemDedxMap[(*cgemDedxTds)] = cgemDedxRoot;
  
    cgemDedxRoot->setTrackId(trackId);
    cgemDedxRoot->setParticleId(particleId);  
    cgemDedxRoot->setStatus (status); 
    cgemDedxRoot->setTruncAlg(truncAlg);
    //dedxRoot->setProb(pb); 
    cgemDedxRoot->setChiE(chiE);       
    cgemDedxRoot->setChiMu(chiMu);     
    cgemDedxRoot->setChiPi(chiPi);     
    cgemDedxRoot->setChiK(chiK);       
    cgemDedxRoot->setChiP(chiP);       
 
    cgemDedxRoot->setNumGoodHits( numGoodHits);     
    cgemDedxRoot->setNumTotalHits( numTotalHits); 
 
    cgemDedxRoot->setProbPH(probPH);             
    cgemDedxRoot->setNormPH(normPH);
    cgemDedxRoot->setErrorPH(errorPH);
    cgemDedxRoot->setTwentyPH(twentyPH);  
    // dedxRoot->setFracErrPH(fracErrPH);       
    //dedxRoot->setMinIronPH(minIronPH);       
    //dedxRoot->setCorrPH(corrPH);             
    log << MSG::INFO << "check Reconstrunction root"<<" particle Id is : "<<particleId
        <<"track id is : "<<trackId
        <<" and status is : "<<status<<endreq;
    
    
    
    recEvt->addCgemDedx(cgemDedxRoot);
  }
 
  return StatusCode::SUCCESS;
}
#endif