DerivationFramework::EGammaClusterCoreCellRecovery Class Reference

#include <EGammaClusterCoreCellRecovery.h>

Inheritance diagram for DerivationFramework::EGammaClusterCoreCellRecovery:

Collaboration diagram for DerivationFramework::EGammaClusterCoreCellRecovery:

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

Private Member Functions
IegammaCellRecoveryTool::Info	decorateObject (const xAOD::Egamma *&egamma) const
StatusCode	findMaxECell (const xAOD::CaloCluster *clus, double &etamax, double &phimax) const

Private Attributes
SG::ReadHandleKey< xAOD::EgammaContainer >	m_SGKey_photons { this, "SGKey_photons", "Photons", "SG key of photon container" }
SG::ReadHandleKey< xAOD::EgammaContainer >	m_SGKey_electrons { this, "SGKey_electrons", "Electrons", "SG key of electron container" }
SG::ReadHandleKey< CaloCellContainer >	m_SGKey_CaloCells { this, "SGKey_CaloCells", "AllCalo", "SG key of calo cell container" }
SG::WriteDecorHandleKeyArray< xAOD::EgammaContainer >	m_SGKey_photons_decorations
SG::WriteDecorHandleKeyArray< xAOD::EgammaContainer >	m_SGKey_electrons_decorations
Gaudi::Property< bool >	m_UseWeightForMaxCell
ToolHandle< IegammaCellRecoveryTool >	m_egammaCellRecoveryTool
	Pointer to the egammaCellRecoveryTool. More...

Detailed Description

Definition at line 25 of file EGammaClusterCoreCellRecovery.h.

Member Function Documentation

addBranches()

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

finaloverridevirtual

Definition at line 64 of file EGammaClusterCoreCellRecovery.cxx.

{
 
  std::vector<SG::WriteDecorHandle<xAOD::EgammaContainer, char>> decon;
  std::vector<SG::WriteDecorHandle<xAOD::EgammaContainer, float>> decoE;
  decon.reserve(2);
  decoE.reserve(2);
 
  // Photon decorations
  if (!m_SGKey_photons.key().empty()) {
 
    for (int i = 0; i < 2; i++) {
      decon.emplace_back(m_SGKey_photons_decorations[i * 2], ctx);
      decoE.emplace_back(m_SGKey_photons_decorations[i * 2 + 1], ctx);
    }
 
    // Retrieve photon container
    SG::ReadHandle<xAOD::EgammaContainer> photonContainer(m_SGKey_photons, ctx);
 
    // Decorate photons
    for (const auto* photon : *photonContainer.ptr()) {
      IegammaCellRecoveryTool::Info res = decorateObject(photon);
      for (int i = 0; i < 2; i++) {
        decon[i](*photon) = res.nCells[i];
        decoE[i](*photon) = res.eCells[i];
      }
    }
  }
 
  // Electron decorations
  if (!m_SGKey_electrons.key().empty()) {
 
    decon.clear(); decon.reserve(2);
    decoE.clear(); decoE.reserve(2);
 
    for (int i = 0; i < 2; i++) {
      decon.emplace_back(m_SGKey_electrons_decorations[i * 2], ctx);
      decoE.emplace_back(m_SGKey_electrons_decorations[i * 2 + 1], ctx);
    }
 
    // Retrieve electron container
    SG::ReadHandle<xAOD::EgammaContainer> electronContainer(m_SGKey_electrons,
                                                            ctx);
 
    // Decorate electrons
    for (const auto* electron : *electronContainer.ptr()) {
      IegammaCellRecoveryTool::Info res = decorateObject(electron);
      for (int i = 0; i < 2; i++) {
        decon[i](*electron) = res.nCells[i];
        decoE[i](*electron) = res.eCells[i];
      }
    }
  }
 
  return StatusCode::SUCCESS;
}

decorateObject()

IegammaCellRecoveryTool::Info DerivationFramework::EGammaClusterCoreCellRecovery::decorateObject ( const xAOD::Egamma *&  egamma ) const

private

Definition at line 122 of file EGammaClusterCoreCellRecovery.cxx.

{
  IegammaCellRecoveryTool::Info info{};
 
  ATH_MSG_DEBUG("Trying to recover cell for object of type " << egamma->type()
                << " pT = " << egamma->pt()
                << " eta = " << egamma->eta()
                << " phi = " << egamma->phi());
 
  const xAOD::CaloCluster *clus = egamma->caloCluster();
  if (!clus) {
    ATH_MSG_WARNING("No associated egamma cluster. Do nothing");
    return info;
  }
 
  // Find max energy cell in layer 2
  double etamax = -999., phimax = -999.;
  if (findMaxECell(clus,etamax,phimax).isFailure()) {
    ATH_MSG_WARNING("Problem in finding maximum energy cell in layer 2");
    return info;
  }
 
  info.etamax = etamax;
  info.phimax = phimax;
  if (m_egammaCellRecoveryTool->execute(*clus,info).isFailure()) {
    ATH_MSG_WARNING("Issue trying to recover cells");
  }
  return info;
}

findMaxECell()

StatusCode DerivationFramework::EGammaClusterCoreCellRecovery::findMaxECell ( const xAOD::CaloCluster *  clus, double &  etamax, double &  phimax  ) const

private

Definition at line 154 of file EGammaClusterCoreCellRecovery.cxx.

{
  const CaloClusterCellLink* cellLinks = clus->getCellLinks();
  if (!cellLinks) {
    ATH_MSG_WARNING("No cell link for cluster. Do nothing");
    return StatusCode::FAILURE;
  }
 
  // First check :
  // if we run in MT, this would be nullptr without initializing a ReadHandleKey
  const CaloCellContainer *caloCells = cellLinks->getCellContainer();
  if (!caloCells) {
    ATH_MSG_WARNING("No cells for cluster. Do nothing");
    return StatusCode::FAILURE;
  }
 
  // Second check :
  // one might have used a custom cell container for the linked cells
  if (cellLinks->getCellContainerLink().dataID() != m_SGKey_CaloCells.key()) {
    ATH_MSG_ERROR("Wrong key for the calo cells");
    return StatusCode::FAILURE;
  }
 
  CaloClusterCellLink::const_iterator it_cell = cellLinks->begin(),
    it_cell_e = cellLinks->end();
 
  // find maximum cell energy in layer 2
  double emax = 0.;
  std::pair<const CaloCell*,double> maxcell{nullptr,0}; //just for debug
  for(; it_cell != it_cell_e; ++it_cell) {
    const CaloCell* cell = (*it_cell);
    if (cell) {
      if (!cell->caloDDE()) {
        ATH_MSG_WARNING("Calo cell without detector element ?? eta = "
                        << cell->eta() << " phi = " << cell->phi());
        continue;
      }
      int layer = cell->caloDDE()->getSampling();
      if (layer == CaloSampling::EMB2 || layer == CaloSampling::EME2) {
        double w     = it_cell.weight();
        double eCell = cell->energy();
        if (m_UseWeightForMaxCell) eCell *= w;
        if (eCell > emax) {
          emax           = eCell;
          maxcell.first  = cell;
          maxcell.second = w;
        }
      }
    }
  }
 
  if (emax > 0) {
    etamax = maxcell.first->caloDDE()->eta_raw();
    phimax = maxcell.first->caloDDE()->phi_raw();
    if (msgLvl(MSG::DEBUG)) {
      CaloSampling::CaloSample sam = maxcell.first->caloDDE()->getSampling();
      double etaAmax = clus->etamax(sam);
      double phiAmax = clus->phimax(sam);
      double vemax   = clus->energy_max(sam);
      ATH_MSG_DEBUG("Cluster energy in sampling 2 = " << clus->energyBE(2)
                    << " maximum layer 2 energy cell, E = " << maxcell.first->energy()
                    << " check E = " << vemax
                    << " w = " << maxcell.second << "\n"
                    << " in calo  frame, eta = " << etamax << " phi = " << phimax << "\n"
                    << " in ATLAS frame, eta = " << etaAmax << " phi = " << phiAmax);
    }
  } else {
    ATH_MSG_WARNING("No layer 2 cell with positive energy ! Should never happen");
    return StatusCode::FAILURE;
  }
  return StatusCode::SUCCESS;
}

initialize()

StatusCode DerivationFramework::EGammaClusterCoreCellRecovery::initialize ( )

finaloverridevirtual

Definition at line 20 of file EGammaClusterCoreCellRecovery.cxx.

{
  ATH_MSG_VERBOSE("initialize() ...");
 
  // The main tool
  ATH_CHECK(m_egammaCellRecoveryTool.retrieve());
 
  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;
  }
 
  ATH_CHECK(m_SGKey_electrons.initialize(SG::AllowEmpty));
  ATH_CHECK(m_SGKey_electrons_decorations.initialize(!m_SGKey_electrons.key().empty()));
  if (!m_SGKey_electrons.key().empty()) {
    ATH_MSG_DEBUG("Using " << m_SGKey_electrons.key() << " for electrons");
    if (msgLvl(MSG::DEBUG)) {
      ATH_MSG_DEBUG("Decorations for " << m_SGKey_electrons.key());
      for (const auto& s : m_SGKey_electrons_decorations)
        { ATH_MSG_DEBUG(s.key()); }
    }
  }
 
  // This tool needs the calo cells linked to the clusters...
  ATH_CHECK(m_SGKey_CaloCells.initialize());
 
  ATH_CHECK(m_SGKey_photons.initialize(SG::AllowEmpty));
  ATH_CHECK(m_SGKey_photons_decorations.initialize(!m_SGKey_photons.key().empty()));
  if (!m_SGKey_photons.key().empty()) {
    ATH_MSG_DEBUG("Using " << m_SGKey_photons.key() << " for photons");
    if (msgLvl(MSG::DEBUG)) {
      ATH_MSG_DEBUG("Decorations for " << m_SGKey_photons.key());
      for (const auto& s : m_SGKey_photons_decorations)
        { ATH_MSG_DEBUG(s.key()); }
    }
  }
 
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_egammaCellRecoveryTool

ToolHandle<IegammaCellRecoveryTool> DerivationFramework::EGammaClusterCoreCellRecovery::m_egammaCellRecoveryTool

private

Initial value:

{
      this,
      "egammaCellRecoveryTool",
      "egammaCellRecoveryTool/egammaCellRecoveryTool",
      "Optional tool that adds cells in L2 or L3 "
      "that could have been rejected by timing cut"
    }

Pointer to the egammaCellRecoveryTool.

Definition at line 73 of file EGammaClusterCoreCellRecovery.h.

m_SGKey_CaloCells

SG::ReadHandleKey<CaloCellContainer> DerivationFramework::EGammaClusterCoreCellRecovery::m_SGKey_CaloCells { this, "SGKey_CaloCells", "AllCalo", "SG key of calo cell container" }

private

Definition at line 41 of file EGammaClusterCoreCellRecovery.h.

m_SGKey_electrons

SG::ReadHandleKey<xAOD::EgammaContainer> DerivationFramework::EGammaClusterCoreCellRecovery::m_SGKey_electrons { this, "SGKey_electrons", "Electrons", "SG key of electron container" }

private

Definition at line 38 of file EGammaClusterCoreCellRecovery.h.

m_SGKey_electrons_decorations

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

private

Initial value:

{
      this,
      "SGKey_electrons_decorations",
      m_SGKey_electrons, {"nadded_Lr2", "Eadded_Lr2", "nadded_Lr3", "Eadded_Lr3"},
      "SG keys for electrons"
    }

Definition at line 53 of file EGammaClusterCoreCellRecovery.h.

m_SGKey_photons

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

private

Definition at line 35 of file EGammaClusterCoreCellRecovery.h.

m_SGKey_photons_decorations

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

private

Initial value:

{
      this,
      "SGKey_photons_decorations",
      m_SGKey_photons, {"nadded_Lr2", "Eadded_Lr2", "nadded_Lr3", "Eadded_Lr3"},
      "SG keys for photon decorations"
    }

Definition at line 45 of file EGammaClusterCoreCellRecovery.h.

m_UseWeightForMaxCell

Gaudi::Property<bool> DerivationFramework::EGammaClusterCoreCellRecovery::m_UseWeightForMaxCell

private

Initial value:

{
      this,
      "UseWeightForMaxCell",
      false,
      "Use the cell weights when finding the L2 max energy cell"
    }

Definition at line 60 of file EGammaClusterCoreCellRecovery.h.

---

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

• EGammaClusterCoreCellRecovery.h
• EGammaClusterCoreCellRecovery.cxx