egammaLargeClusterMaker.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/

 
#ifndef EGAMMATOOLS_EGAMMALARGECLUSTERMAKER_H
#define EGAMMATOOLS_EGAMMALARGECLUSTERMAKER_H
 
#include "AthenaBaseComps/AthAlgTool.h"
 
#include "CaloDetDescr/CaloDetDescrManager.h"
 
#include "xAODCaloEvent/CaloClusterContainer.h"
#include "xAODCaloEvent/CaloClusterFwd.h"
 
#include "CaloUtils/CaloClusterCollectionProcessor.h"
#include "GaudiKernel/ToolHandle.h"
 
#include "StoreGate/ReadCondHandleKey.h"
#include "StoreGate/ReadHandleKey.h"
 
class CaloCellContainer;
class CaloClusterCellLink;
 
class egammaLargeClusterMaker final
  : public AthAlgTool
  , virtual public CaloClusterCollectionProcessor
{
public:
  egammaLargeClusterMaker(const std::string& type,
                          const std::string& name,
                          const IInterface* parent);
 
  using CaloClusterCollectionProcessor::execute;

  virtual StatusCode initialize() override final;
  virtual StatusCode execute(
    const EventContext& ctx,
    xAOD::CaloClusterContainer* collection) const override final;
 
private:
  SG::ReadHandleKey<xAOD::CaloClusterContainer> m_inputClusterCollection{
    this,
    "InputClusterCollection",
    "egammaClusters",
    "The name of the cluster container for electrons and photons"
  };

  SG::ReadHandleKey<CaloCellContainer> m_cellsKey{
    this,
    "CellsName",
    "AllCalo",
    "Names of containers which contain cells"
  };
 
  SG::ReadCondHandleKey<CaloDetDescrManager> m_caloDetDescrMgrKey{
    this,
    "CaloDetDescrManager",
    "CaloDetDescrManager",
    "SG Key for CaloDetDescrManager in the Condition Store"
  };

  Gaudi::Property<bool> m_isFWD{ this,
                                 "doFWDelesurraundingWindows",
                                 false,
                                 "Save FWD electron surraunding windows" };
 
  Gaudi::Property<double> m_neta{
    this,
    "Neta",
    7.0,
    "Number of eta cells in each sampling in which to look for hottest cell"
  };
 
  Gaudi::Property<double> m_nphi{
    this,
    "Nphi",
    7.0,
    "Number of phi cell in each sampling in which to look for hottest cell"
  };

  Gaudi::Property<double> m_centEtThr{
    this,
    "CentralEtThreshold",
    3000.,
    "Value of Et cut on input central cluster to produce large cluster" };
 
  Gaudi::Property<double> m_netaFWD{
    this,
    "NetaFWD",
    6.0,
    "Number of eta cells in each sampling in which to look for hottest cell"
  };
 
  Gaudi::Property<double> m_nphiFWD{
    this,
    "NphiFWD",
    6.0,
    "Number of phi cell in each sampling in which to look for hottest cell"
  };
 
  Gaudi::Property<double> m_drFWD{
    this,
    "deltaR_FCAL",
    0.4,
    "Cone size to collec cell around hottest-cell FCAL"
  };
 
  Gaudi::Property<double> m_drEM{
    this,
    "deltaR_EM",
    0.3,
    "Cone size to eventually collect EME1 cells around hottest-cell"
  };

  Gaudi::Property<double> m_fwdEtThr{
    this,
    "ForwardEtThreshold",
    5000.,
    "Value of Et cut on input forward cluster to produce large cluster" };

  Gaudi::Property<bool> m_addCellsFromOtherSamplings{
    this,
    "AddCellsFromOtherSamplings",
    true,
    "Keep all cells from relevant samplings or just from the max sampling" };
 
};
 
#endif // EGAMMATOOLS_EMCLUSTERTOOL_H