CaloSwPhioff_v2.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
 
#include "CaloSwPhioff_v2.h"
#include "CaloClusterCorrection/interpolate.h"
#include "CaloDetDescr/CaloDetDescrManager.h"
#include "CaloGeoHelpers/CaloPhiRange.h"
#include <cmath> 
 
 
using xAOD::CaloCluster;
using CaloClusterCorr::interpolate;
 

void CaloSwPhioff_v2::makeTheCorrection (const Context& myctx,
                                         CaloCluster* cluster,
                                         const CaloDetDescrElement* elt,
                                         float /*eta*/,
                                         float adj_eta,
                                         float phi,
                                         float /*adj_phi*/,
                                         CaloSampling::CaloSample samp) const
{
  // Calculate the correction.
  float offs = m_correction_coef(myctx) * energy_interpolation
    (cluster->e(),
     Builder (m_correction (myctx),
              m_interp_barriers (myctx),
              m_degree (myctx),
              std::abs (adj_eta)),
     m_energies(myctx),
     m_energy_degree(myctx));
 
  // Flip the sign, if needed.
  if (m_flip_phi(myctx) && elt->eta_raw() < 0)
    offs = -offs;
 
  // Apply the correction.
  cluster->setPhi (samp, CaloPhiRange::fix (phi + offs));
}
 

CaloSwPhioff_v2::Builder::Builder (const CxxUtils::Array<3>& correction,
                                   const CxxUtils::Array<1>& interp_barriers,
                                   int degree,
                                   float aeta)
  : m_correction (correction),
    m_interp_barriers (interp_barriers),
    m_degree (degree),
    m_aeta (aeta)
{
}
 

float CaloSwPhioff_v2::Builder::calculate (int energy_ndx, bool& good) const
{
  good = true;
  return interpolate (m_correction[energy_ndx],
                      m_aeta,
                      m_degree,
                      1,
                      m_interp_barriers);
}