xAODPhotonAuxContainerCnv_v1.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
 
// System include(s):
#include <stdexcept>
 
// Gaudi/Athena include(s):
#include "GaudiKernel/MsgStream.h"
 
// Core EDM include(s):
#include "AthContainers/tools/copyAuxStoreThinned.h"
 
// Local include(s):
#include "xAODPhotonAuxContainerCnv_v1.h"
#include "xAODEgamma/PhotonContainer.h"
#include "xAODEgamma/versions/PhotonContainer_v1.h"
 
 
xAODPhotonAuxContainerCnv_v1::xAODPhotonAuxContainerCnv_v1()
{
}
 
void xAODPhotonAuxContainerCnv_v1::
persToTrans( const xAOD::PhotonAuxContainer_v1* oldObj,
             xAOD::PhotonAuxContainer* newObj,
             MsgStream& /*log*/ ) const {
 
   // Clear the transient object:
   newObj->resize( 0 );
 
   // Copy the payload of the v1 object into the latest one by misusing
   // the thinning code a bit...
   SG::copyAuxStoreThinned( *oldObj, *newObj, nullptr );
 
   // Set up interface containers on top of them:
 
   //The old  uses v_
   xAOD::PhotonContainer_v1 oldInt;
   for( size_t i = 0; i < oldObj->size(); ++i ) { 
     oldInt.push_back( new xAOD::Photon_v1() ); 
   }
   oldInt.setStore( oldObj );
   
   xAOD::PhotonContainer newInt;
   for( size_t i = 0; i < newObj->size(); ++i ) { 
     newInt.push_back( new xAOD::Photon() ); 
   }
   newInt.setStore( newObj );
 
   for( size_t i = 0; i < oldInt.size(); ++i ) {
     float e237 = oldInt[ i ]->showerShapeValue(xAOD::EgammaParameters::e237);
     float e277 = oldInt[ i ]->showerShapeValue(xAOD::EgammaParameters::e277);
     float Reta= ( e277 != 0 ? e237/e277 : 0);
     newInt[ i ]->setShowerShapeValue(Reta, xAOD::EgammaParameters::Reta);
     //
     float e233 = oldInt[ i ]->showerShapeValue(xAOD::EgammaParameters::e233);
     float Rphi= ( e237 != 0 ? e233/e237 : 0);
     newInt[ i ]->setShowerShapeValue(Rphi, xAOD::EgammaParameters::Rphi);
     //
     float emax = oldInt[ i ]->showerShapeValue(xAOD::EgammaParameters::emaxs1);
     float emax2 = oldInt[ i ]->showerShapeValue(xAOD::EgammaParameters::e2tsts1);
     float Eratio = fabs(emax+emax2)>0. ? (emax-emax2)/(emax+emax2) : 0.;
     newInt[ i ]->setShowerShapeValue(Eratio, xAOD::EgammaParameters::Eratio);
     //
     float emin = oldInt[ i ]->showerShapeValue(xAOD::EgammaParameters::emins1);
     float DeltaE = emax2-emin;
     newInt[ i ]->setShowerShapeValue(DeltaE, xAOD::EgammaParameters::DeltaE);
     //
     const xAOD::CaloCluster* cluster  = oldInt[ i ]->caloCluster();
     if (cluster) {
       const float eta2   = fabsf(cluster->etaBE(2));
       // transverse energy in calorimeter (using eta position in second sampling)
       const double energy =  cluster->e();
       double et = 0.;
       if (eta2<999.) {
         const double cosheta = cosh(eta2);
         et = (cosheta != 0.) ? energy /cosheta : 0.;
       }
       float ethad = oldInt[ i ]->showerShapeValue(xAOD::EgammaParameters::ethad);
       float Rhad =  fabs(et) > 0. ? ethad/et : 0.;
       newInt[ i ]->setShowerShapeValue(Rhad, xAOD::EgammaParameters::Rhad);
       //
       float ethad1 = oldInt[ i ]->showerShapeValue(xAOD::EgammaParameters::ethad1);
       float Rhad1 = fabs(et) > 0. ? ethad1/et : 0.;
       newInt[ i ]->setShowerShapeValue(Rhad1, xAOD::EgammaParameters::Rhad1); 
     }
   }
 
   return;
}
 
void xAODPhotonAuxContainerCnv_v1::
transToPers( const xAOD::PhotonAuxContainer*,
             xAOD::PhotonAuxContainer_v1*,
             MsgStream& log ) const {
 
   log << MSG::ERROR
       << "Somebody called xAODPhotonAuxContainerCnv_v1::transToPers"
       << endmsg;
   throw std::runtime_error( "Somebody called xAODPhotonAuxContainerCnv_v1::"
                             "transToPers" );
 
   return;
}