RecEventCnv.cc

Go to the documentation of this file.

#include "RecEventCnv.h"
#include "PerBaseEvent/HeaderObjectCnv.h"
#include "Conventions/Electronics.h"

using namespace DayaBay;
using namespace std;

RecHeaderCnv::RecHeaderCnv(ISvcLocator* svc)
    : RootIOTypedCnv<PerRecHeader,RecHeader>("PerRecHeader",
                                             classID(),svc)
{
}

RecHeaderCnv::~RecHeaderCnv()
{
}


StatusCode RecHeaderCnv::PerToTran(const PerRecHeader& perobj,
                                   DayaBay::RecHeader& tranobj)
{
    MsgStream log(msgSvc(), "RecHeaderCnv::PerToTran");
    StatusCode sc = HeaderObjectCnv::toTran(perobj, tranobj);
    if(sc.isFailure()) return sc;

    DayaBay::RecTrigger& recTr = tranobj.recTrigger();
    recTr.setHeader(&tranobj);
    Detector det((Site::Site_t)perobj.site,
                 (DetectorId::DetectorId_t)perobj.detector);
    recTr.setDetector(det);
    recTr.setTriggerType(Trigger::TriggerType_t (perobj.triggerType));
    recTr.setTriggerNumber( perobj.triggerNumber );
    
    timespec value = {perobj.triggerTimeSec, perobj.triggerTimeNanoSec};
    TimeStamp tstamp(value);
    recTr.setTriggerTime(tstamp);
    
    recTr.setEnergy( perobj.energy );
    recTr.setRawEvis( perobj.rawEvis );
    recTr.setEnrec( perobj.enrec );
    recTr.setEprec( perobj.eprec );
    recTr.setPosition( CLHEP::HepLorentzVector(perobj.x, perobj.y, 
                                               perobj.z, perobj.t) );
    recTr.setDirection( CLHEP::HepLorentzVector(perobj.dx, perobj.dy, 
                                                perobj.dz, perobj.dt) );
    recTr.setEnergyStatus( perobj.energyStatus );
    recTr.setPositionStatus( perobj.positionStatus );
    recTr.setDirectionStatus( perobj.directionStatus );
    recTr.setEnergyQuality( perobj.energyQuality );
    recTr.setPositionQuality( perobj.positionQuality );
    recTr.setDirectionQuality( perobj.directionQuality );
    if(perobj.errorMatrixDim>0){
      CLHEP::HepMatrix errMatrix(perobj.errorMatrixDim,
                                 perobj.errorMatrixDim,
                                 0);
      for(int row=1; row<perobj.errorMatrixDim; row++){
        for(int col=1; col<perobj.errorMatrixDim; col++){
          int rcIndex = (row-1)*perobj.errorMatrixDim + (col-1);
          errMatrix(row,col) = perobj.errorMatrix[rcIndex];
        }
      }
      recTr.setErrorMatrix( errMatrix );
    }else{
      recTr.setErrorMatrix( CLHEP::HepMatrix() );
    }
    return StatusCode::SUCCESS;
}

StatusCode RecHeaderCnv::TranToPer(const DayaBay::RecHeader& tranobj,
                                   PerRecHeader& perobj)
{
    MsgStream log(msgSvc(), "RecHeaderCnv::TranToPer");
    StatusCode sc = HeaderObjectCnv::toPer(tranobj,perobj);
    if (sc.isFailure()) return sc;
    const DayaBay::RecTrigger& recTr = tranobj.recTrigger();

    perobj.site = (recTr.detector()).site();
    perobj.detector = (recTr.detector()).detectorId();
    perobj.triggerType = recTr.triggerType();
    perobj.triggerNumber = recTr.triggerNumber();
    perobj.triggerTimeSec = recTr.triggerTime().GetSec();
    perobj.triggerTimeNanoSec = recTr.triggerTime().GetNanoSec();
    
    perobj.energy = recTr.energy();
    perobj.rawEvis = recTr.rawEvis();
    perobj.enrec = recTr.enrec();
    perobj.eprec = recTr.eprec();
    perobj.x = recTr.position().x();
    perobj.y = recTr.position().y();
    perobj.z = recTr.position().z();
    perobj.t = recTr.position().t();
    perobj.dx = recTr.direction().x();
    perobj.dy = recTr.direction().y();
    perobj.dz = recTr.direction().z();
    perobj.dt = recTr.direction().t();
    perobj.energyStatus = recTr.energyStatus();
    perobj.positionStatus = recTr.positionStatus();
    perobj.directionStatus = recTr.directionStatus();
    perobj.energyQuality = recTr.energyQuality();
    perobj.positionQuality = recTr.positionQuality();
    perobj.directionQuality = recTr.directionQuality();
    perobj.errorMatrixDim = recTr.errorMatrix().num_row();
    perobj.errorMatrix.resize(perobj.errorMatrixDim*perobj.errorMatrixDim);
    if(perobj.errorMatrixDim>0){
      for(int row=1; row<perobj.errorMatrixDim; row++){
        for(int col=1; col<perobj.errorMatrixDim; col++){
          int rcIndex = (row-1)*perobj.errorMatrixDim + (col-1);
          perobj.errorMatrix[rcIndex] = recTr.errorMatrix()(row,col);
        }
      }
    }

    return StatusCode::SUCCESS;
}

StatusCode RecHeaderCnv::fillRepRefs(IOpaqueAddress*, DataObject* dobj) {
  
  MsgStream log(msgSvc(), "RecHeaderCnv::fillRepRefs");
  //MsgStream log(msgSvc(), "fillRepRefs");
  RecHeader* rh = dynamic_cast<RecHeader*>(dobj);

  log << MSG::DEBUG
      << "Saving links to " << rh->inputHeaders().size() 
      << " input headers" << endreq;
  
  StatusCode sc = HeaderObjectCnv::fillPer(m_rioSvc,*rh,*m_perOutObj);

  if (sc.isFailure()) {
    log << MSG::ERROR << "Failed to fill HeaderObject part" << endreq;
    return sc;
  }
  
  // ... fill RecHeader part...
  return sc;
}

StatusCode RecHeaderCnv::fillObjRefs(IOpaqueAddress*, DataObject* dobj) {
  //MsgStream log(msgSvc(), "RecHeaderCnv::fillObjRefs");
  MsgStream log(msgSvc(), "fillObjRefs");
  HeaderObject* hobj = dynamic_cast<HeaderObject*>(dobj);

  StatusCode sc = HeaderObjectCnv::fillTran(m_rioSvc,*m_perInObj,*hobj);
  if (sc.isFailure()) {
    log << MSG::ERROR << "Failed to fill HeaderObject part" << endreq;
    return sc;
  }
  
  log << MSG::DEBUG
      << "Restored links to " << hobj->inputHeaders().size() 
      << " input headers" << endreq;
  
  // ... fill RecHeader part...
  return sc;
}