RecRpcEventCnv.cc

Go to the documentation of this file.

#include "RecRpcEventCnv.h"
#include "PerBaseEvent/HeaderObjectCnv.h"
#include "Conventions/Detectors.h"
#include <algorithm>

using namespace DayaBay;

RecRpcEventCnv::RecRpcEventCnv(ISvcLocator* svc)
    : RootIOTypedCnv<PerRecRpcHeader,RecRpcHeader>("PerRecRpcHeader",
                                             classID(),svc)
{
}

RecRpcEventCnv::~RecRpcEventCnv()
{
}


StatusCode RecRpcEventCnv::PerToTran(const PerRecRpcHeader& perobj,
                                   RecRpcHeader& tranobj)
{
    MsgStream log(msgSvc(), "RecRpcEventCnv::PerToTran");
    StatusCode sc = HeaderObjectCnv::toTran(perobj, tranobj);
    if(sc.isFailure()) return sc;

    RecRpcTrigger& recTr = tranobj.recTrigger();
    recTr.setHeader(&tranobj);
    Detector det((Site::Site_t)perobj.site,
                 (DetectorId::DetectorId_t)perobj.detector);
    recTr.setDetector(det);
    recTr.setTriggerNumber( perobj.triggerNumber );
    recTr.setTriggerType(Trigger::TriggerType_t (perobj.triggerType));

    timespec value = {perobj.triggerTimeSec, perobj.triggerTimeNanoSec};
    TimeStamp tstamp(value);
    recTr.setTriggerTime(tstamp);

    recTr.setEventType(perobj.eventType);

    for ( int i = 0; i < perobj.nClusters; ++i ) {
        const PerRecRpcHeader::RpcCluster& per = perobj.clusters.at(i);
        RecRpcCluster* tran = new RecRpcCluster(
                CLHEP::HepLorentzVector(per.x, per.y, per.z, per.t),
                ReconStatus::ReconStatus_t(per.status),
                per.quality,
                RecRpcCluster::ClusterType(per.type));
        tran->setMaxHitsInLayer(per.maxHitsInLayer);
        tran->setIs4of4(per.is4of4);
        recTr.addCluster(tran);
    }

    const std::vector<const RecRpcCluster*>& clusters = recTr.clusters();
    for ( int i = 0; i < perobj.nTracks; ++i ) {
        const PerRecRpcHeader::RpcTrack& per = perobj.tracks.at(i);
        recTr.addTrack(RecRpcTrack(clusters[per.start], clusters[per.end]));
    }

    return StatusCode::SUCCESS;
}

StatusCode RecRpcEventCnv::TranToPer(const RecRpcHeader& tranobj,
                                   PerRecRpcHeader& perobj)
{
    MsgStream log(msgSvc(), "RecRpcEventCnv::TranToPer");
    StatusCode sc = HeaderObjectCnv::toPer(tranobj,perobj);
    if (sc.isFailure()) return sc;
    const RecRpcTrigger& recTr = tranobj.recTrigger();

    perobj.site = recTr.detector().site();
    perobj.detector = recTr.detector().detectorId();
    perobj.triggerNumber = recTr.triggerNumber();
    perobj.triggerType = recTr.triggerType();
    perobj.triggerTimeSec = recTr.triggerTime().GetSec();
    perobj.triggerTimeNanoSec = recTr.triggerTime().GetNanoSec();

    perobj.combinationType = recTr.combinationType();
    perobj.eventType = recTr.eventType();

    perobj.reset(recTr.numCluster(), recTr.numTrack());

    const std::vector<const RecRpcCluster*>& clusters = recTr.clusters();
    for ( int i = 0; i < perobj.nClusters; ++i ) {
        const RecRpcCluster* iCluster = clusters.at(i);
        perobj.clusters[i].status = iCluster->status();
        perobj.clusters[i].type = iCluster->type();
        perobj.clusters[i].maxHitsInLayer = iCluster->maxHitsInLayer();
        perobj.clusters[i].is4of4 = iCluster->is4of4();
        perobj.clusters[i].x    = iCluster->position().x();
        perobj.clusters[i].y    = iCluster->position().y();
        perobj.clusters[i].z    = iCluster->position().z();
        perobj.clusters[i].t    = iCluster->position().t();
        perobj.clusters[i].quality = iCluster->quality();
    }

    std::vector<const RecRpcCluster*>::const_iterator begin = clusters.begin();
    std::vector<const RecRpcCluster*>::const_iterator end = clusters.end();
    const std::vector<RecRpcTrack>& tracks = recTr.tracks();
    for ( int i = 0; i < perobj.nTracks; ++i ) {
        perobj.tracks[i].start = std::find(begin, end, tracks[i].start()) - begin;
        perobj.tracks[i].end = std::find(begin, end, tracks[i].end()) - begin;
    }

    return StatusCode::SUCCESS;
}

StatusCode RecRpcEventCnv::fillRepRefs(IOpaqueAddress*, DataObject* tobj) {
  
  MsgStream log(msgSvc(), "RecRpcEventCnv::fillRepRefs");
  RecRpcHeader* rh = dynamic_cast<RecRpcHeader*>(tobj);

  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 RecRpcHeader part...
  return sc;
}

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

  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 RecRpcHeader part...
  return sc;
}