eg_resolution.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef EG_RESOLUTION_H
#define EG_RESOLUTION_H
 
#include <cmath>
#include <cstdlib>
#include <memory>
#include <array>
#include "TArrayD.h"
#include "TFile.h"
#include "TH1.h"
#include "xAODEgamma/Egamma.h"
 
class eg_resolution
{
 
public:
  static constexpr size_t samplings = 4;
  static constexpr size_t resolTypes = 3;
  eg_resolution(const std::string& configuration);
  ~eg_resolution() = default;

  double getResolution(int particle_type,
                       double energy,
                       double eta,
                       int resol_type = 2) const;
  double getResolution(const xAOD::Egamma& particle, int resol_type = 2) const;
 
private:
  // histograms to store resolution parameters
  std::array<std::array<std::unique_ptr<TH1>, resolTypes>, samplings> m_hSampling;
  std::array<std::array<std::unique_ptr<TH1>, resolTypes>, samplings> m_hNoise;
  std::array<std::array<std::unique_ptr<TH1>, resolTypes>, samplings> m_hConst;
  const TArrayD* m_etaBins; // Not owning pointer
};
 
#endif