CaloSwEtaoff_v3.h

Go to the documentation of this file.

// This file's extension implies that it's C, but it's really -*- C++ -*-.
/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef CALOCLUSTERCORRECTION_CALOSWETAOFF_V3_H
#define CALOCLUSTERCORRECTION_CALOSWETAOFF_V3_H
 
#include "CaloClusterCorrection/CaloClusterCorrectionCommon.h"
#include "CaloConditions/Array.h"
 
 
class CaloDetDescrElement;
 

class CaloSwEtaoff_v3
  : public CaloClusterCorrectionCommon
{
public:
  using CaloClusterCorrectionCommon::CaloClusterCorrectionCommon;
  

  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 override;
 
private:
  class Builder
    : public TableBuilder
  {
  public:
    Builder (const CxxUtils::Array<4>& correction,
             const CxxUtils::Array<2>& regions,
             const CxxUtils::Array<2>& forms,
             float aeta,
             float u,
             int region_ndx);

    virtual float calculate (int energy_ndx, bool& good) const;
 
  private:
    static float calc0 (float aeta, float u, const CaloRec::Array<2>& coef) ;

    static float calc3 (float aeta, float u, const CaloRec::Array<2>& coef) ;

    static float calc4 (float aeta, float u, const CaloRec::Array<2>& coef) ;

    static float calc5 (float aeta, float u, const CaloRec::Array<2>& coef) ;

    static float calc10 (float aeta, float u, const CaloRec::Array<2>& coef) ;

    static float calc11 (float aeta, float u, const CaloRec::Array<2>& coef) ;

    static float calc13 (float aeta, float u, const CaloRec::Array<2>& coef, float xlo, float xhi) ;
 
    // Tabulated arrays of function parameters.
    const CxxUtils::Array<4>& m_correction;

    const CxxUtils::Array<2>& m_regions;

    const CxxUtils::Array<2>& m_forms;

    int m_region_ndx;

    float m_aeta;

    float m_u;

    int m_form;
  };
  friend class Builder;

  static int find_region (const CxxUtils::Array<2>& regions,
                   float aeta) ;

  Constant<CxxUtils::Array<4> > m_correction
  { this, "correction", "Tabulated arrays of function parameters." };

  Constant<CxxUtils::Array<2> > m_regions
  { this, "regions", "Table of regions." };
 
  enum {
    REG_LO = 0,
    REG_HI = 1,
    REG_FORM = 2,
    REG_CELLSIZE = 3,
    REG_PHASE = 4,
    REG_ENTRIES = 5 // Number of entries expected.
  };

  Constant<CxxUtils::Array<1> > m_energies
  { this, "energies", "Table of energies at which the correction was tabulated." };

  Constant<int> m_energy_degree
  { this, "energy_degree", "Degree of the polynomial interpolation in energy." };

  Constant<CxxUtils::Array<2> > m_forms
  { this, "forms", "Functional form to use per region per energy." };
};
 
 
#endif // not CALOCLUSTERCORRECTION_CALOSWETAOFF_V3_H