DerivationFramework::ClusterEnergyPerLayerDecorator Class Reference

Decorate egamma objects with the energy per layer for a rectangular cluster of size neta X nphi built on the fly. More...

#include <ClusterEnergyPerLayerDecorator.h>

Inheritance diagram for DerivationFramework::ClusterEnergyPerLayerDecorator:

Collaboration diagram for DerivationFramework::ClusterEnergyPerLayerDecorator:

Public Member Functions
	ClusterEnergyPerLayerDecorator (const std::string &t, const std::string &n, const IInterface *p)
	~ClusterEnergyPerLayerDecorator ()
StatusCode	initialize ()
StatusCode	finalize ()
virtual StatusCode	addBranches () const

Private Member Functions
std::vector< float >	decorateObject (const EventContext &ctx, const xAOD::Egamma *egamma, const CaloCellContainer *cellCont) const

Private Attributes
int	m_eta_size = 0
int	m_phi_size = 0
std::vector< unsigned int >	m_layers
SG::ReadHandleKey< xAOD::EgammaContainer >	m_SGKey_photons { this, "SGKey_photons", "", "SG key of photon container" }
SG::ReadHandleKey< xAOD::EgammaContainer >	m_SGKey_electrons
SG::ReadHandleKey< CaloCellContainer >	m_SGKey_caloCells
ToolHandle< CaloClusterProcessor >	m_caloFillRectangularTool
SG::WriteDecorHandleKeyArray< xAOD::EgammaContainer >	m_SGKey_photons_decorations
SG::WriteDecorHandleKeyArray< xAOD::EgammaContainer >	m_SGKey_electrons_decorations
const CaloFillRectangularCluster *	m_tool = nullptr

Detailed Description

Decorate egamma objects with the energy per layer for a rectangular cluster of size neta X nphi built on the fly.

Definition at line 30 of file ClusterEnergyPerLayerDecorator.h.

Constructor & Destructor Documentation

ClusterEnergyPerLayerDecorator()

DerivationFramework::ClusterEnergyPerLayerDecorator::ClusterEnergyPerLayerDecorator	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Definition at line 19 of file ClusterEnergyPerLayerDecorator.cxx.

  : base_class(t, n, p)
{
  declareProperty("neta", m_eta_size);
  declareProperty("nphi", m_phi_size);
  declareProperty("layers", m_layers = { 0, 1, 2, 3 });
}

~ClusterEnergyPerLayerDecorator()

DerivationFramework::ClusterEnergyPerLayerDecorator::~ClusterEnergyPerLayerDecorator ( )

default

Member Function Documentation

addBranches()

StatusCode DerivationFramework::ClusterEnergyPerLayerDecorator::addBranches ( ) const

virtual

Definition at line 97 of file ClusterEnergyPerLayerDecorator.cxx.

{
  const EventContext& ctx = Gaudi::Hive::currentContext();
 
  // Retrieve cell container
 
  SG::ReadHandle<CaloCellContainer> cellContainer(m_SGKey_caloCells, ctx);
  const CaloCellContainer* cellCont = cellContainer.ptr();
 
  // 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;
    for (unsigned int i = 0; i < m_layers.size(); i++) {
      decorations.emplace_back(
        
          m_SGKey_photons_decorations[i], ctx);
    }
 
    // Decorate photons
    for (const auto* photon : *importedPhotons) {
      std::vector<float> result = decorateObject(ctx, photon, cellCont);
      for (unsigned int i = 0; i < m_layers.size(); i++) {
        decorations[i](*photon) = result[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;
    for (unsigned int i = 0; i < m_layers.size(); i++) {
      decorations.emplace_back(
        
          m_SGKey_electrons_decorations[i], ctx);
    }
 
    // Decorate electrons
    for (const auto* electron : *importedElectrons) {
      std::vector<float> result = decorateObject(ctx, electron, cellCont);
      for (unsigned int i = 0; i < m_layers.size(); i++) {
        decorations[i](*electron) = result[i];
      }
    }
  }
 
  return StatusCode::SUCCESS;
}

decorateObject()

std::vector< float > DerivationFramework::ClusterEnergyPerLayerDecorator::decorateObject ( const EventContext &  ctx, const xAOD::Egamma *  egamma, const CaloCellContainer *  cellCont  ) const

private

Definition at line 161 of file ClusterEnergyPerLayerDecorator.cxx.

{
  std::vector<float> result;
  result.clear();
  if (not egamma or not egamma->caloCluster())
    return result;
 
  std::unique_ptr<xAOD::CaloCluster> egcClone;
  if (not egamma->author(xAOD::EgammaParameters::AuthorCaloTopo35)) {
    egcClone =
      CaloClusterStoreHelper::makeCluster(cellCont,
                                          egamma->caloCluster()->eta0(),
                                          egamma->caloCluster()->phi0(),
                                          egamma->caloCluster()->clusterSize());
    m_tool->makeCorrection(ctx, egcClone.get());
  }
 
  for (unsigned int layer : m_layers) {
    result.emplace_back(egcClone ? egcClone->energyBE(layer) : 0.);
  }
 
  return result;
}

finalize()

StatusCode DerivationFramework::ClusterEnergyPerLayerDecorator::finalize

(

)

Definition at line 89 of file ClusterEnergyPerLayerDecorator.cxx.

{
 
  return StatusCode::SUCCESS;
}

initialize()

StatusCode DerivationFramework::ClusterEnergyPerLayerDecorator::initialize

(

)

Definition at line 37 of file ClusterEnergyPerLayerDecorator.cxx.

{
  ATH_MSG_VERBOSE("initialize() ...");
 
  if (m_SGKey_photons.empty() and m_SGKey_electrons.empty()) {
    ATH_MSG_FATAL("No e-gamma collection provided for thinning. At least one "
                  "egamma collection (photons/electrons) must be provided!");
    return StatusCode::FAILURE;
  }
 
  ATH_CHECK(m_caloFillRectangularTool.retrieve());
  // how to add a statement to theck that the tool indeed has the size matching
  // the neta and nphi properties? ATH_MSG_DEBUG("CaloFillRectangularCluster
  // print size ...
  m_tool =
    dynamic_cast<const CaloFillRectangularCluster*>(
      &(*m_caloFillRectangularTool));
  if (m_tool==nullptr) {
    ATH_MSG_ERROR("Pointer to CaloFillRectantularCluster tool is invalid");
    return StatusCode::FAILURE;
  }
 
  ATH_CHECK(m_SGKey_caloCells.initialize());
 
  if (!m_SGKey_electrons.key().empty()) {
    ATH_MSG_DEBUG("Using " << m_SGKey_electrons << " for electrons");
    ATH_CHECK(m_SGKey_electrons.initialize());
 
    const char* containerKey = m_SGKey_electrons.key().c_str();
    for (int layer : m_layers) {
      m_SGKey_electrons_decorations.emplace_back(
        Form("%s.E%dx%d_Lr%d", containerKey, m_eta_size, m_phi_size, layer));
    }
    ATH_CHECK(m_SGKey_electrons_decorations.initialize());
  }
 
  if (!m_SGKey_photons.key().empty()) {
    ATH_MSG_DEBUG("Using " << m_SGKey_photons << " for photons");
    ATH_CHECK(m_SGKey_photons.initialize());
 
    const char* containerKey = m_SGKey_photons.key().c_str();
    for (int layer : m_layers) {
      m_SGKey_photons_decorations.emplace_back(
        Form("%s.E%dx%d_Lr%d", containerKey, m_eta_size, m_phi_size, layer));
    }
    ATH_CHECK(m_SGKey_photons_decorations.initialize());
  }
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_caloFillRectangularTool

ToolHandle<CaloClusterProcessor> DerivationFramework::ClusterEnergyPerLayerDecorator::m_caloFillRectangularTool

private

Initial value:

{
    this,
    "CaloFillRectangularClusterTool",
    "",
    "Handle of the CaloFillRectangularClusterTool"
  }

Definition at line 63 of file ClusterEnergyPerLayerDecorator.h.

m_eta_size

int DerivationFramework::ClusterEnergyPerLayerDecorator::m_eta_size = 0

private

Definition at line 42 of file ClusterEnergyPerLayerDecorator.h.

m_layers

std::vector<unsigned int> DerivationFramework::ClusterEnergyPerLayerDecorator::m_layers

private

Definition at line 44 of file ClusterEnergyPerLayerDecorator.h.

m_phi_size

int DerivationFramework::ClusterEnergyPerLayerDecorator::m_phi_size = 0

private

Definition at line 43 of file ClusterEnergyPerLayerDecorator.h.

m_SGKey_caloCells

SG::ReadHandleKey<CaloCellContainer> DerivationFramework::ClusterEnergyPerLayerDecorator::m_SGKey_caloCells

private

Initial value:

{
    this,
    "SGKey_caloCells",
    "AllCalo",
    "SG key of the cell container"
  }

Definition at line 56 of file ClusterEnergyPerLayerDecorator.h.

m_SGKey_electrons

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

private

Initial value:

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

Definition at line 49 of file ClusterEnergyPerLayerDecorator.h.

m_SGKey_electrons_decorations

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

private

Initial value:

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

Definition at line 79 of file ClusterEnergyPerLayerDecorator.h.

m_SGKey_photons

SG::ReadHandleKey<xAOD::EgammaContainer> DerivationFramework::ClusterEnergyPerLayerDecorator::m_SGKey_photons { this, "SGKey_photons", "", "SG key of photon container" }

private

Definition at line 47 of file ClusterEnergyPerLayerDecorator.h.

m_SGKey_photons_decorations

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

private

Initial value:

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

Definition at line 71 of file ClusterEnergyPerLayerDecorator.h.

m_tool

const CaloFillRectangularCluster* DerivationFramework::ClusterEnergyPerLayerDecorator::m_tool = nullptr

private

Definition at line 86 of file ClusterEnergyPerLayerDecorator.h.

---

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

• ClusterEnergyPerLayerDecorator.h
• ClusterEnergyPerLayerDecorator.cxx