RecCgemDedxCnv.cxx

Go to the documentation of this file.

#ifndef RecCgemDedxCnv_CXX
#define RecCgemDedxCnv_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/RecCgemDedx.h"
 
#include "RootEventData/TRecCgemDedx.h" // standard root object
#include "RootEventData/TRecTrackEvent.h"
 
#include "RootCnvSvc/Rec/RecTrackCnv.h"
#include "RootCnvSvc/Rec/RecCgemDedxCnv.h"
#include "RootCnvSvc/RootAddress.h"
 
#include <vector>
 
using namespace std;
 
// Instantiation of a static factory class used by clients to create
// instances of this service
//static CnvFactory<RecCgemDedxCnv> s_factory;
//const ICnvFactory& RecCgemDedxCnvFactory = s_factory;
 
RecCgemDedxCnv::RecCgemDedxCnv(ISvcLocator* svc)
: RootEventBaseCnv(classID(), svc)
{
 
    // Here we associate this converter with the /Event path on the TDS.
    MsgStream log(msgSvc(), "RecCgemDedxCnv");
    //log << MSG::DEBUG << "Constructor called for " << objType() << endreq;
    //m_rootTreename ="Rec";
    m_rootBranchname ="m_recCgemDedxCol";
    // declareObject(EventModel::Recon::RecCgemDedxCol, objType(), m_rootTreename, m_rootBranchname);
    m_adresses.push_back(&m_recCgemDedxCol);
    m_recCgemDedxCol=0;
}
 
StatusCode RecCgemDedxCnv::TObjectToDataObject(DataObject*& refpObject) {
  // creation of TDS object from root object
    MsgStream log(msgSvc(), "RecCgemDedxCnv");
    log << MSG::DEBUG << "RecCgemDedxCnv::TObjectToDataObject" << endreq;
    StatusCode sc=StatusCode::SUCCESS;
 
    // create the TDS location for the Dedx Collection
    RecCgemDedxCol* recCgemDedxCol = new RecCgemDedxCol;
    refpObject=recCgemDedxCol;
 
 
    // now convert
    if (!m_recCgemDedxCol) return sc;
    TIter dedxIter(m_recCgemDedxCol);
    TRecCgemDedx *recCgemDedxRoot = 0;
    while ((recCgemDedxRoot = (TRecCgemDedx*)dedxIter.Next())) {
     double dedxHit = recCgemDedxRoot->dedxHit();
     double dedxEsat = recCgemDedxRoot->dedxEsat();
     double dedxNoRun = recCgemDedxRoot->dedxNoRun();
     double dedxMoment = recCgemDedxRoot->dedxMoment();        
 
     int trackId =  recCgemDedxRoot->trackId();     
     int particleId = recCgemDedxRoot->particleId(); 
     int status = recCgemDedxRoot->status();
     int truncAlg = recCgemDedxRoot->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] = recCgemDedxRoot->chiE();
     chi[1] = recCgemDedxRoot->chiMu();
     chi[2] = recCgemDedxRoot->chiPi();
     chi[3] = recCgemDedxRoot->chiK();
     chi[4] = recCgemDedxRoot->chiP();
     int numGoodHits = recCgemDedxRoot->numGoodHits(); 
     int numTotalHits = recCgemDedxRoot->numTotalHits();
 
     double probPH = recCgemDedxRoot->probPH();      
     double normPH = recCgemDedxRoot->normPH();
     double errorPH = recCgemDedxRoot->errorPH();
     double twentyPH = recCgemDedxRoot->twentyPH(); 
     double dedxExpect[5],sigmaDedx[5],pidProb[5];
     for (int i=0; i<5; i++){
     dedxExpect[i] = recCgemDedxRoot->dedxExpect(i);
     sigmaDedx[i] = recCgemDedxRoot->sigmaDedx(i);
     pidProb[i] = recCgemDedxRoot->pidProb(i);}
     
     RecCgemDedx *recCgemDedx = new RecCgemDedx();
     m_common.m_rootRecCgemDedxMap[recCgemDedxRoot] = recCgemDedx;
 
     recCgemDedx->setTrackId(trackId);
     recCgemDedx->setParticleId(particleId);
     recCgemDedx->setStatus (status);
     recCgemDedx->setTruncAlg(truncAlg);
     // dedxTds->setProb(pb); 
     // dedxTds->setNumSigmaE(numSigmaE);       
     // dedxTds->setNumSigmaMu(numSigmaMu);     
     // dedxTds->setNumSigmaPi(numSigmaPi);     
     // dedxTds->setNumSigmaK(numSigmaK);       
     // dedxTds->setNumSigmaP(numSigmaP);       
     recCgemDedx->setChi(chi);
     recCgemDedx->setNumGoodHits( numGoodHits);     
     recCgemDedx->setNumTotalHits( numTotalHits); 
 
     recCgemDedx->setProbPH(probPH);             
     recCgemDedx->setNormPH(normPH);
     recCgemDedx->setErrorPH(errorPH);
     recCgemDedx->setTwentyPH(twentyPH);
     //for (int i=0; i<5; i++){
     recCgemDedx->setDedxExpect(dedxExpect);
     recCgemDedx->setSigmaDedx(sigmaDedx);
     recCgemDedx->setPidProb(pidProb);
     
     recCgemDedx->setDedxHit(dedxHit);   
     recCgemDedx->setDedxEsat(dedxEsat);  
     recCgemDedx->setDedxNoRun(dedxNoRun); 
     recCgemDedx->setDedxMoment(dedxMoment);     
 
     //}
              
             
     
     recCgemDedxCol->push_back(recCgemDedx); 
     //delete dedxTds;
     // dedxTds = NULL;
     }
    //m_dedxCol->Delete();
    delete m_recCgemDedxCol;
    m_recCgemDedxCol = 0;
 
    return StatusCode::SUCCESS;
}
 
 
StatusCode RecCgemDedxCnv::DataObjectToTObject(DataObject* obj,RootAddress* rootaddr) {
  MsgStream log(msgSvc(), "RecCgemDedxCnv");
  log << MSG::DEBUG << "RecCgemDedxCnv::DataObjectToTObject" << endreq;
  StatusCode sc=StatusCode::SUCCESS;
 
  RecCgemDedxCol * recCgemDedxCol=dynamic_cast<RecCgemDedxCol *> (obj);
  if (!recCgemDedxCol) {
    log << MSG::ERROR << "Could not downcast to RecCgemDedxCol" << endreq;
    return StatusCode::FAILURE;
  }
 
  DataObject *evt;
  m_eds->findObject(EventModel::Recon::Event,evt);
  if (evt==NULL) {
    log << MSG::ERROR << "Could not get RecEvent in TDS "  << endreq;
    return StatusCode::FAILURE;
  }
  ReconEvent * devtTds=dynamic_cast<ReconEvent *> (evt);
  if (!devtTds) {
    log << MSG::ERROR << "RecCgemDedxCnv:Could not downcast to TDS RecEvent" << endreq;
  }
  IOpaqueAddress *addr;
 
  m_cnvSvc->getRecTrackCnv()->createRep(evt,addr); 
  TRecTrackEvent *recEvt=m_cnvSvc->getRecTrackCnv()->getWriteObject();
 
  const TObjArray *m_recCgemDedxCol = recEvt->getRecCgemDedxCol();
  if (!m_recCgemDedxCol) return sc;
  recEvt->clearRecCgemDedxCol(); //necessary in case there is I/O at the same time since array is static
  RecCgemDedxCol::const_iterator recCgemDedx;
 
  for (recCgemDedx = recCgemDedxCol->begin(); recCgemDedx != recCgemDedxCol->end(); recCgemDedx++) {
    Int_t trackId = (*recCgemDedx)->trackId();
    Int_t particleId = (*recCgemDedx)->particleId(); 
    Int_t status = (*recCgemDedx)->status();     
    Int_t truncAlg = (*recCgemDedx)->truncAlg();
    // Double_t pb[5];
    // for (Int_t i = 0; i < 5; i++) 
    // pb[i] = (*dedxTds)->prob(i);
    Double_t dedxHit    = (*recCgemDedx)->getDedxHit();
    Double_t dedxEsat   = (*recCgemDedx)->getDedxEsat();
    Double_t dedxNoRun  = (*recCgemDedx)->getDedxNoRun();
    Double_t dedxMoment = (*recCgemDedx)->getDedxMoment();    
 
                 
    Double_t chiE  = (*recCgemDedx)->chi(0);   
    Double_t chiMu = (*recCgemDedx)->chi(1);  
    Double_t chiPi = (*recCgemDedx)->chi(2);  
    Double_t chiK = (*recCgemDedx)->chi(3);   
    Double_t chiP = (*recCgemDedx)->chi(4);   
  
    Int_t numGoodHits = (*recCgemDedx)->numGoodHits(); 
    Int_t numTotalHits = (*recCgemDedx)->numTotalHits();
 
    Double_t probPH = (*recCgemDedx)->probPH();
    Double_t normPH = (*recCgemDedx)->normPH();
    Double_t errorPH = (*recCgemDedx)->errorPH();
    Double_t twentyPH = (*recCgemDedx)->twentyPH();
    //Double_t fracErrPH = (*dedxTds)-> fracErrPH();  
    //Double_t minIronPH = (*dedxTds)->minIronPH();   
    //Double_t corrPH = (*dedxTds)->corrPH();      
    double dedxExpect[5],sigmaDedx[5],pidProb[5],chi[5];
   for (int i=0; i<5; i++){
      chi[i] = (*recCgemDedx)->chi(i);
      dedxExpect[i] = (*recCgemDedx)->getDedxExpect(i);
      sigmaDedx[i] = (*recCgemDedx)->getSigmaDedx(i);
      pidProb[i] = (*recCgemDedx)->getPidProb(i);
  }
                    
       
    TRecCgemDedx *recCgemDedxRoot = new TRecCgemDedx();
    //m_common.m_recCgemDedxMap[(*recCgemDedx)] = recCgemDedxRoot;
    //std::cout<<"check write to Reconstrunction root particle Id is "<<dedxRoot->particleId()<<endl; 
    recCgemDedxRoot->setTrackId(trackId);
    recCgemDedxRoot->setParticleId(particleId);  
    recCgemDedxRoot->setStatus (status); 
    recCgemDedxRoot->setTruncAlg(truncAlg);
    log << MSG::INFO<<"check Reconstrunction root"<<" particle Id is : "<<particleId
        <<"track id is : "<<trackId
        <<" and status is : "<<status<<endreq; 
    //dedxRoot->setProb(pb); 
    recCgemDedxRoot->setChiE(chiE);       
    recCgemDedxRoot->setChiMu(chiMu);     
    recCgemDedxRoot->setChiPi(chiPi);     
    recCgemDedxRoot->setChiK(chiK);       
    recCgemDedxRoot->setChiP(chiP);       
 
    recCgemDedxRoot->setNumGoodHits( numGoodHits);     
    recCgemDedxRoot->setNumTotalHits(   numTotalHits); 
 
    recCgemDedxRoot->setProbPH(probPH);             
    recCgemDedxRoot->setNormPH(normPH);
    recCgemDedxRoot->setErrorPH(errorPH);
    recCgemDedxRoot->setTwentyPH(twentyPH);
    //  for (int i=0; i<5; i++){
    recCgemDedxRoot->setChi(chi);
    recCgemDedxRoot->setDedxExpect(dedxExpect);
    recCgemDedxRoot->setSigmaDedx(sigmaDedx);
    recCgemDedxRoot->setPidProb(pidProb);
  
    recCgemDedxRoot->setDedxHit(dedxHit);   
    recCgemDedxRoot->setDedxEsat(dedxEsat);
    recCgemDedxRoot->setDedxNoRun(dedxNoRun);
    recCgemDedxRoot->setDedxMoment(dedxMoment);    
    
    //}
    log << MSG::INFO<<"check Reconstrunction root"<<" particle Id is : "<<particleId
              <<"track id is : "<<trackId
          <<" and status is : "<<status<<endreq;
 
    // dedxRoot->setFracErrPH(fracErrPH);       
    //dedxRoot->setMinIronPH(minIronPH);       
    //dedxRoot->setCorrPH(corrPH);             
    recEvt->addRecCgemDedx(recCgemDedxRoot);
  }
 
  return StatusCode::SUCCESS;
}
 
#endif