MissingETCnv_p2.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2019 CERN for the benefit of the ATLAS collaboration
*/
 
/********************************************************************
 
NAME:     MissingETCnv_p2.cxx
PACKAGE:  offline/Reconstruction/RecTPCnv
 
PURPOSE:  Transient/Persisten converter for MissingET class
********************************************************************/
 
// AthenaPoolCnvSvc includes
#include "AthenaPoolCnvSvc/T_AthenaPoolTPConverter.h"
 
// MissingETEvent includes
#include "MissingETEvent/MissingET.h"
#include "MissingETEvent/MissingEtRegions.h"
 
// RecTPCnv includes
#include "RecTPCnv/MissingETCnv_p2.h"
#include "RecTPCnv/MissingEtRegionsCnv_p2.h"
 
 
#include <stdlib.h>
 
// region converter
static const MissingEtRegionsCnv_p2 regCnv;
 
 
 
void MissingETCnv_p2::persToTrans(  const MissingET_p2* pers, MissingET* trans, MsgStream& /* msg */ ) const {
    std::vector<float>::const_iterator it=pers->m_AllTheData.begin();
    persToTrans(trans, it);
  return;
}
 
void MissingETCnv_p2::transToPers(  const MissingET* trans, MissingET_p2* pers, MsgStream& /* msg */ ) const {
    transToPers(trans, pers->m_AllTheData); 
 return;
}
 
 
void MissingETCnv_p2::persToTrans( MissingET* trans, std::vector<float>::const_iterator i) const {
    union conv{  unsigned int i;   float f;  } c;
    c.f=(*i); ++i;
        MissingET::Source source = static_cast<MissingET::Source>(c.i>>1);
    double ex     = (*i);++i;
    double ey     = (*i);++i;
    double etSum  = (*i);++i;
        std::unique_ptr<MissingEtRegions> regions;
    if( c.i & 1) {
          regions = std::make_unique<MissingEtRegions>();
          regCnv.persToTrans( regions.get(), i);
    }
 
        *trans = MissingET (source,
                            std::move (regions),
                            ex, ey, etSum);
    // std::cout<<"READING: source: "<<trans->m_source<<"\t et: "<<trans->m_etSum;
    // if (trans->m_regions != 0) std::cout<<"\t has Regions"<<std::endl; else std::cout<<"\t NO Regions"<<std::endl;
    return;
}
 
void  MissingETCnv_p2::transToPers(  const MissingET* trans,  std::vector<float> &all ) const {
        unsigned int tmp=static_cast<unsigned int>(trans->getSource());
    tmp <<= 1;
    if( trans->getRegions() != 0) tmp|=1;
    union conv{  unsigned int i;   float f;  } c;
    c.i=tmp;
    all.push_back(c.f);
    all.push_back(trans->etx());
    all.push_back(trans->ety());
    all.push_back(trans->sumet());
    if( trans->getRegions() != 0) regCnv.transToPers(trans->getRegions(), all);
    // std::cout<<"WRITING: source: "<<trans->m_source<<"\t et: "<<trans->m_etSum;
    // if (trans->m_regions != 0) std::cout<<"\t has Regions"<<std::endl; else std::cout<<"\t NO Regions"<<std::endl;
    return;
}