TgcCoinDataCnv_p3.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonTrigCoinData/TgcCoinData.h"
#include "MuonEventTPCnv/MuonTrigCoinData/TgcCoinDataCnv_p3.h"
#include <algorithm>
#include <atomic>
#include "AthLinks/tools/IdentContIndex.h"
#include "TrkEventTPCnv/helpers/EigenHelpers.h"
 
void TgcCoinDataCnv_p3::
persToTrans( const Muon::TgcCoinData_p3 *persObj, Muon::TgcCoinData *transObj,MsgStream &log ) 
{
   log << MSG::DEBUG << "TgcCoinDataCnv_p3::persToTrans" << endmsg;
   
   Muon::TgcCoinData::CoinDataType type = Muon::TgcCoinData::TYPE_UNKNOWN;
   switch(persObj->m_type){
   case 0: type=Muon::TgcCoinData::TYPE_TRACKLET; break;
   case 1: type=Muon::TgcCoinData::TYPE_HIPT; break;
   case 2: type=Muon::TgcCoinData::TYPE_SL; break;
   case 3: type=Muon::TgcCoinData::TYPE_UNKNOWN; break;
   case 4: type=Muon::TgcCoinData::TYPE_TRACKLET_EIFI; break;
   default: break;
   }
 
   auto posIn = std::unique_ptr<Amg::Vector2D>
     (createTransFromPStore( &m_localPosCnv, persObj->m_posIn, log ));
   auto posOut = std::unique_ptr<Amg::Vector2D>
     (createTransFromPStore( &m_localPosCnv, persObj->m_posOut, log ));
 
   Trk::ErrorMatrix dummy;
   fillTransFromPStore( &m_errorMxCnv, persObj->m_errMat, &dummy, log );
   auto tempMat = std::make_unique<Amg::MatrixX>();
   if (!dummy.values.empty())
     EigenHelpers::vectorToEigenMatrix(dummy.values, *tempMat, "TgcCoinDataCnv_p3");
 
   *transObj = Muon::TgcCoinData (Identifier(persObj->m_channelIdIn),
                                  Identifier(persObj->m_channelIdOut),
                                  persObj->m_collectionIdHash,
                                  transObj->detectorElementIn(),
                                  transObj->detectorElementOut(),
                                  type,
                                  persObj->m_isAside,
                                  persObj->m_phi,
                                  persObj->m_isInner,
                                  persObj->m_isForward,
                                  persObj->m_isStrip,
                                  persObj->m_trackletId,
                                  posIn.release(),
                                  posOut.release(),
                                  persObj->m_widthIn,
                                  persObj->m_widthOut,
                                  persObj->m_delta,
                                  persObj->m_sub,
                                  persObj->m_inner,
                                  persObj->m_roi,
                                  persObj->m_pt,
                                  persObj->m_trackletIdStrip,
                                  tempMat.release(),
                                  persObj->m_isPositiveDeltaR,
                                  persObj->m_veto
                                  );
   transObj->setHashAndIndex(persObj->m_indexAndHash);
}
 
void TgcCoinDataCnv_p3::
transToPers( const Muon::TgcCoinData *transObj, Muon::TgcCoinData_p3 *persObj, MsgStream &log )
{
  log << MSG::DEBUG << "TgcCoinDataCnv_p3::transToPers" << endmsg;
  
  // Tracklet and HiPt should have valid channelIdIn  
  if((transObj->type()==Muon::TgcCoinData::TYPE_TRACKLET || 
      transObj->type()==Muon::TgcCoinData::TYPE_HIPT || 
      transObj->type()==Muon::TgcCoinData::TYPE_TRACKLET_EIFI)&&
     (!transObj->isInner())) {
    if(!transObj->channelIdIn().get_identifier32().is_valid()) {
      log << MSG::ERROR << "TgcCoinDataCnv_p3::transToPers - invalid trans channelIdIn!! " 
          << (transObj->type()==Muon::TgcCoinData::TYPE_TRACKLET ? "(TRACKLET)" : 
              (transObj->type()==Muon::TgcCoinData::TYPE_TRACKLET_EIFI ? "(TRACKLET_EIFI)" : "(HIPT)"))
          << endmsg;
    }
  }
 
  // Tracklet, HiPt and SL should have valid channelIdOut 
  if((transObj->type()==Muon::TgcCoinData::TYPE_TRACKLET ||
      transObj->type()==Muon::TgcCoinData::TYPE_HIPT || 
      transObj->type()==Muon::TgcCoinData::TYPE_SL) &&
     (!transObj->isInner())) {
    if(!transObj->channelIdOut().get_identifier32().is_valid()) {
      log << MSG::ERROR << "TgcCoinDataCnv_p3::transToPers - invalid trans channelIdOut!! " 
          << (transObj->type()==Muon::TgcCoinData::TYPE_TRACKLET ? "(TRACKLET)" : 
              (transObj->type()==Muon::TgcCoinData::TYPE_HIPT ? "(HIPT)" : "(SL)"))
          << endmsg;
    }
  }
 
  persObj->m_channelIdIn = transObj->channelIdIn().get_identifier32().get_compact(); 
  persObj->m_channelIdOut = transObj->channelIdOut().get_identifier32().get_compact(); 
  persObj->m_collectionIdHash = transObj->identifyHash();
  persObj->m_indexAndHash=transObj->getHashAndIndex().hashAndIndex();
  switch(transObj->type()){
  case Muon::TgcCoinData::TYPE_TRACKLET: persObj->m_type=0; break;
  case Muon::TgcCoinData::TYPE_HIPT: persObj->m_type=1;  break;
  case Muon::TgcCoinData::TYPE_SL: persObj->m_type=2; break;
  case Muon::TgcCoinData::TYPE_UNKNOWN: persObj->m_type=3; break;
  case Muon::TgcCoinData::TYPE_TRACKLET_EIFI: persObj->m_type=4; break;
  default: break;
  }  
 
  persObj->m_isAside = transObj->isAside(); 
  persObj->m_phi = transObj->phi(); 
  persObj->m_isInner = transObj->isInner(); 
  persObj->m_isForward = transObj->isForward(); 
  persObj->m_isStrip = transObj->isStrip(); 
  persObj->m_trackletId = transObj->trackletId(); 
  persObj->m_trackletIdStrip = transObj->trackletIdStrip(); 
  persObj->m_widthIn = transObj->widthIn(); 
  persObj->m_widthOut = transObj->widthOut(); 
  persObj->m_posIn = toPersistent( &m_localPosCnv, transObj->posInPtr(), log );
  persObj->m_posOut = toPersistent( &m_localPosCnv, transObj->posOutPtr(), log );
  // persObj->m_errMat = toPersistent( &m_errorMxCnv, transObj->m_errMat, log );
   
  if (transObj->hasErrMat() )
  {
    if (! transObj->errMat().cols() || !transObj->errMat().rows()) {
      static std::atomic<unsigned int> numErrsPrinted=0;
      if (numErrsPrinted<10) {
        if (numErrsPrinted++ < 10) {
          log << MSG::WARNING << "Have been given a TgcCoinData error matrix with size = ["
              <<transObj->errMat().rows()<<","<<transObj->errMat().cols()<<"]. Skipping writing of the matrix."<<endmsg;
        }
      }
    } else {
      Trk::ErrorMatrix pMat;
      EigenHelpers::eigenMatrixToVector(pMat.values, transObj->errMat(), "TgcCoinDataCnv_p3");
      persObj->m_errMat = toPersistent( &m_errorMxCnv, &pMat, log );
    }
  }
   
  persObj->m_delta = transObj->delta(); 
  persObj->m_roi = transObj->roi(); 
  persObj->m_pt = transObj->pt(); 
  persObj->m_veto = transObj->veto(); 
  persObj->m_inner = transObj->inner(); 
 
  persObj->m_sub = transObj->sub(); 
  persObj->m_isPositiveDeltaR = transObj->isPositiveDeltaR(); 
}