CaloScaleCluster.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
 
#include "CaloScaleCluster.h"
#include "CaloClusterCorrection/interpolate.h"
#include <cmath>
 
 
using xAOD::CaloCluster;
using CaloClusterCorr::interpolate;
 

void CaloScaleCluster::makeTheCorrection (const Context& myctx,
                                          CaloCluster* cluster,
                                          const CaloDetDescrElement* /*elt*/,
                                          float /*eta*/,
                                          float adj_eta,
                                          float /*phi*/,
                                          float /*adj_phi*/,
                                          CaloSampling::CaloSample /*samp*/)
  const
{
  const CxxUtils::Array<2> correction = m_correction (myctx);
  if (correction.size() == 0) {
    ATH_MSG_ERROR("Empty correction array");
    return;
  }
  
  float adj_aeta = std::abs (adj_eta);
  float fac;
  if (std::abs (adj_eta) > m_etamax (myctx)) {
    fac = correction[correction.size()-1][1];
  }
  else {
    // The rfac-v5 tables contain a couple spurious zeros
    // (from empty bins)?  If we don't do anything, this correction will blow
    // up in the immediate vicinity of those points.  Work around by replacing
    // the zeros with an average of the adjacent points, by setting
    // fixZero=true in the interpolate call.
    fac = interpolate (correction, adj_aeta, m_degree (myctx),
                       1, CaloRec::Array<1>(), -1,
                       /*fixZero=*/ true);
  }
  cluster->setE (cluster->e() / fac);
}