DerivationFramework::GainDecorator Class Reference

#include <GainDecorator.h>

Inheritance diagram for DerivationFramework::GainDecorator:

Collaboration diagram for DerivationFramework::GainDecorator:

Classes
struct	calculation

Public Member Functions
virtual StatusCode	initialize () override final
virtual StatusCode	addBranches (const EventContext &ctx) const override final

Static Public Member Functions
static int	getLayer (const CaloCell *cell)

Private Member Functions
calculation	decorateObject (const xAOD::Egamma *&egamma) const

Private Attributes
SG::ReadHandleKey< xAOD::EgammaContainer >	m_SGKey_photons
SG::ReadHandleKey< xAOD::EgammaContainer >	m_SGKey_electrons
Gaudi::Property< std::map< int, std::string > >	m_gainNames
Gaudi::Property< std::vector< unsigned int > >	m_layers {this, "layers", { 0, 1, 2, 3 }}
std::vector< std::pair< int, int > >	m_names_E
SG::WriteDecorHandleKeyArray< xAOD::EgammaContainer >	m_SGKey_photons_decorations
SG::WriteDecorHandleKeyArray< xAOD::EgammaContainer >	m_SGKey_electrons_decorations

Detailed Description

Definition at line 26 of file GainDecorator.h.

Member Function Documentation

addBranches()

StatusCode DerivationFramework::GainDecorator::addBranches ( const EventContext & ctx ) const

finaloverridevirtual

Definition at line 51 of file GainDecorator.cxx.

{
  ATH_MSG_DEBUG("m_gainNames.size() = " << m_gainNames.size() << ", m_layers.size() = "<< m_layers.size() );
  const size_t limit{m_gainNames.size()*m_layers.size()};
  if (limit == 0) { return StatusCode::FAILURE; } // FIXME Add ERROR message here
  // Photon decorations
 
  if (!m_SGKey_photons.key().empty()) {
 
    // Retrieve photon container
    SG::ReadHandle<xAOD::EgammaContainer> photonContainer(m_SGKey_photons, ctx);
    const xAOD::EgammaContainer* importedPhotons = photonContainer.ptr();
 
    // Setup vectors of photon decorations
    std::vector<SG::WriteDecorHandle<xAOD::EgammaContainer, float>>
      decorations_E;
    std::vector<SG::WriteDecorHandle<xAOD::EgammaContainer, float>>
      decorations_rnoW;
    std::vector<SG::WriteDecorHandle<xAOD::EgammaContainer, char>>
      decorations_nCells;
    for (unsigned int i = 0; i < limit; ++i) {
      decorations_E.emplace_back(
                                 m_SGKey_photons_decorations[i * 3], ctx);
      decorations_rnoW.emplace_back(
                                    m_SGKey_photons_decorations[i * 3 + 1], ctx);
      decorations_nCells.emplace_back(
                                      m_SGKey_photons_decorations[i * 3 + 2], ctx);
    }
    // Decorate photons
    int i(0);
    for (const auto* photon : *importedPhotons) {
      DerivationFramework::GainDecorator::calculation res =
        decorateObject(photon);
      i = 0;
      for (const auto& kv : m_gainNames) {
        for (const auto layer : m_layers) {
          std::pair<int, int> key(kv.first, layer);
          decorations_E[i](*photon) = res.E[key];
          decorations_rnoW[i](*photon) =
            res.EnoW[key] != 0 ? res.E[key]/res.EnoW[key] : 1;
          decorations_nCells[i](*photon) = res.nCells[key];
          i++;
        }
      }
    }
  }
 
  // Electron decorations
 
  if (!m_SGKey_electrons.key().empty()) {
 
    // Retrieve electron container
    SG::ReadHandle<xAOD::EgammaContainer> electronContainer(m_SGKey_electrons,
                                                            ctx);
    const xAOD::EgammaContainer* importedElectrons = electronContainer.ptr();
 
    // Setup vectors of electron decorations
    std::vector<SG::WriteDecorHandle<xAOD::EgammaContainer, float>>
      decorations_E;
    std::vector<SG::WriteDecorHandle<xAOD::EgammaContainer, float>>
      decorations_rnoW;
    std::vector<SG::WriteDecorHandle<xAOD::EgammaContainer, char>>
      decorations_nCells;
    for (unsigned int i = 0; i < limit; ++i) {
      decorations_E.emplace_back(
                                 m_SGKey_electrons_decorations[i * 3], ctx);
      decorations_rnoW.emplace_back(
                                    m_SGKey_electrons_decorations[i * 3 + 1], ctx);
      decorations_nCells.emplace_back(
                                      m_SGKey_electrons_decorations[i * 3 + 2], ctx);
    }
    // Decorate electrons
    int i(0);
    for (const auto* electron : *importedElectrons) {
      DerivationFramework::GainDecorator::calculation res =
        decorateObject(electron);
      i = 0;
      for (const auto& kv : m_gainNames) {
        for (const auto layer : m_layers) {
          std::pair<int, int> key(kv.first, layer);
          decorations_E[i](*electron) = res.E[key];
          decorations_rnoW[i](*electron) =
            res.EnoW[key] != 0 ? res.E[key]/res.EnoW[key] : 1;
          decorations_nCells[i](*electron) = res.nCells[key];
          i++;
        }
      }
    }
  }
 
  return StatusCode::SUCCESS;
}

decorateObject()

DerivationFramework::GainDecorator::calculation DerivationFramework::GainDecorator::decorateObject ( const xAOD::Egamma *& egamma ) const

private

Definition at line 145 of file GainDecorator.cxx.

{
 
  // Compute energy and number of cells per gain per layer
  // Set the initial values to 0 (needed?)
  DerivationFramework::GainDecorator::calculation result;
  for (const auto& kv : m_names_E) {
    result.E[kv] = 0.;
    result.EnoW[kv] = 0.;
    result.nCells[kv] = 0;
  }
 
  // Skip the computation for missing cell links (like topo-seeded photons)
  // but decorate anyway
  const CaloClusterCellLink* cellLinks =
    egamma->caloCluster() ? egamma->caloCluster()->getCellLinks() : nullptr;
  if (cellLinks) {
    CaloClusterCellLink::const_iterator it = cellLinks->begin(),
      itE = cellLinks->end();
    for (; it != itE; ++it) {
      const CaloCell *cell = *it;
      if (!cell)
        continue;
      std::pair<int, int> key(static_cast<int>(cell->gain()), getLayer(cell));
      // Increment the corresponding entry (not important if it is not
      // initialised)
      double weight = it.weight();
      result.E[key] += cell->energy()*weight;
      result.EnoW[key] += cell->energy();
      result.nCells[key]++;
    }
  }
 
  return result;
}

getLayer()

int DerivationFramework::GainDecorator::getLayer ( const CaloCell * cell )

static

Definition at line 183 of file GainDecorator.cxx.

{
  int sampling =
    (cell && cell->caloDDE() ? cell->caloDDE()->getSampling() : -1);
  if (sampling == CaloSampling::PreSamplerB ||
      sampling == CaloSampling::PreSamplerE)
    return 0;
  if (sampling == CaloSampling::EMB1 || sampling == CaloSampling::EME1)
    return 1;
  if (sampling == CaloSampling::EMB2 || sampling == CaloSampling::EME2)
    return 2;
  if (sampling == CaloSampling::EMB3 || sampling == CaloSampling::EME3)
    return 3;
  return -1;
}

initialize()

StatusCode DerivationFramework::GainDecorator::initialize ( )

finaloverridevirtual

Definition at line 18 of file GainDecorator.cxx.

{
  ATH_MSG_VERBOSE("initialize() ...");
 
  if (m_SGKey_photons.key().empty() && m_SGKey_electrons.key().empty()) {
    ATH_MSG_FATAL("No e-gamma collection provided for thinning. At least one "
                  "egamma collection (photons/electrons) must be provided!");
    return StatusCode::FAILURE;
  }
 
  for (const auto& kv : m_gainNames) {
    for (const auto layer : m_layers) {
      m_names_E.emplace_back(kv.first, layer);
    }
  }
 
  ATH_CHECK(m_SGKey_electrons.initialize(SG::AllowEmpty));
  if (!m_SGKey_electrons.key().empty()) {
    ATH_MSG_DEBUG("Using " << m_SGKey_electrons << " for electrons");
  }
  ATH_CHECK(m_SGKey_electrons_decorations.initialize(!m_SGKey_electrons.key().empty()));
 
  ATH_CHECK(m_SGKey_photons.initialize(SG::AllowEmpty));
  if (!m_SGKey_photons.key().empty()) {
    ATH_MSG_DEBUG("Using " << m_SGKey_photons << " for photons");
  }
  ATH_CHECK(m_SGKey_photons_decorations.initialize(!m_SGKey_photons.key().empty()));
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_gainNames

Gaudi::Property< std::map<int, std::string> > DerivationFramework::GainDecorator::m_gainNames

private

Initial value:

{this, "gain_names",
                                                              { { CaloGain::LARHIGHGAIN, "Hi" },
                                                                { CaloGain::LARMEDIUMGAIN, "Med" },
                                                                { CaloGain::LARLOWGAIN, "Low" } } }

Definition at line 50 of file GainDecorator.h.

                                                       {this, "gain_names",
                                                              { { CaloGain::LARHIGHGAIN, "Hi" },
                                                                { CaloGain::LARMEDIUMGAIN, "Med" },
                                                                { CaloGain::LARLOWGAIN, "Low" } } };

m_layers

Gaudi::Property< std::vector<unsigned int> > DerivationFramework::GainDecorator::m_layers {this, "layers", { 0, 1, 2, 3 }}

private

Definition at line 54 of file GainDecorator.h.

{this, "layers", { 0, 1, 2, 3 }};

m_names_E

std::vector<std::pair<int, int> > DerivationFramework::GainDecorator::m_names_E

private

Definition at line 57 of file GainDecorator.h.

m_SGKey_electrons

SG::ReadHandleKey<xAOD::EgammaContainer> DerivationFramework::GainDecorator::m_SGKey_electrons

private

Initial value:

{
      this, "SGKey_electrons", "", "SG key of electron container" }

Definition at line 47 of file GainDecorator.h.

                                                          {
      this, "SGKey_electrons", "", "SG key of electron container" };

m_SGKey_electrons_decorations

SG::WriteDecorHandleKeyArray<xAOD::EgammaContainer> DerivationFramework::GainDecorator::m_SGKey_electrons_decorations

private

Initial value:

{
      this,
      "SGKey_electrons_decorations",
      m_SGKey_electrons, {},
      "SG keys for electrons decorations not really configurable"
    }

Definition at line 68 of file GainDecorator.h.

                                 {
      this,
      "SGKey_electrons_decorations",
      m_SGKey_electrons, {},
      "SG keys for electrons decorations not really configurable"
    };

m_SGKey_photons

SG::ReadHandleKey<xAOD::EgammaContainer> DerivationFramework::GainDecorator::m_SGKey_photons

private

Initial value:

{
      this, "SGKey_photons", "", "SG key of photon container" }

Definition at line 44 of file GainDecorator.h.

                                                        {
      this, "SGKey_photons", "", "SG key of photon container" };

m_SGKey_photons_decorations

SG::WriteDecorHandleKeyArray<xAOD::EgammaContainer> DerivationFramework::GainDecorator::m_SGKey_photons_decorations

private

Initial value:

{
      this,
      "SGKey_photons_decorations",
      m_SGKey_photons, {},
      "SG keys for photon decorations not really configurable"
    }

Definition at line 60 of file GainDecorator.h.

                               {
      this,
      "SGKey_photons_decorations",
      m_SGKey_photons, {},
      "SG keys for photon decorations not really configurable"
    };

---

The documentation for this class was generated from the following files:

• GainDecorator.h
• GainDecorator.cxx