CaloTopoEMphioff.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef CALOCLUSTERCORRECTION_CALOTOPOEMPHIOFF_H
#define CALOCLUSTERCORRECTION_CALOTOPOEMPHIOFF_H
/********************************************************************
 
NAME:     CaloTopoEMphioff.h
PACKAGE:  offline/Calorimeter/CaloClusterCorrection
 
AUTHORS:  M.Boonekamp, N.Kerschen
CREATED:  March 2005
 
PURPOSE:  correction for the phi offset due to accordion structure
          base class: CaloClusterCorrection (AlgTool)
 
Updated:  March 12, 2005   (MB)
          corrections for the TopoCluster 
********************************************************************/
 
#include "CaloClusterCorrection/CaloClusterCorrectionCommon.h"
#include "CaloConditions/Array.h"
class CaloCluster;
 
class CaloTopoEMphioff : public CaloClusterCorrectionCommon
{
public:
  using CaloClusterCorrectionCommon::CaloClusterCorrectionCommon;
 
 
  // virtual method in CaloClusterCorrection
  virtual void makeTheCorrection(const Context& myctx,
                                 xAOD::CaloCluster* cluster,
                 const CaloDetDescrElement* elt,
                 float eta,
                 float adj_eta,
                 float phi,
                                 float adj_phi,
                 CaloSampling::CaloSample samp) const;
 
 private:
  // Comments on the parametrization :
  // The energy dependence has been fitted in eta bins of 0.1. The function that have been found to best describe this dependence are:
  // f(x) = A*(sqrt(1/x)) + B   for   (0 < eta < 0.8) where x is the energy in GeV
  // f(x) = A*(1/x) + B  for  (0.8 < eta < 2.3) 
  // f(x) = A*x + B
  // For the fits, abs(offset) has been taken, so that the adequate sign should be given to the correction (-1) for (0.8 < eta < 2.3)
  // For the G4 data, the sign should be flipped for eta < 0
  // So : 
  // EdepA is A in 25 eta bins
  // EdepB is B in 25 eta bins
  // Granularity is the Granularity
  // EtaFrontier[i] is 0.8, 2.3, 2.5 for i=0,1,2
  Constant<CxxUtils::Array<1> > m_EdepA        { this, "EdepA", "" };
  Constant<CxxUtils::Array<1> > m_EdepB        { this, "EdepB", "" };
  Constant<CxxUtils::Array<1> > m_EtaFrontier  { this, "EtaFrontier", "" };
 
  Constant<float> m_Granularity    { this, "Granularity",  "" };
  Constant<int> m_FlipPhi          { this, "FlipPhi",      "" };
  Constant<float> m_EndcapOffset   { this, "EndcapOffset", "" };
};
 
#endif