egammaTopoClusterCopier.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef EGAMMAALGS_EGAMMATOPOCLUSTERCOPIER_H
#define EGAMMAALGS_EGAMMATOPOCLUSTERCOPIER_H
 
#include "xAODCaloEvent/CaloClusterContainer.h"
 
#include "AthenaBaseComps/AthReentrantAlgorithm.h"
#include "StoreGate/ReadHandleKey.h"
#include "StoreGate/WriteHandleKey.h"
#include "AthContainers/ConstDataVector.h"
 
#include "GaudiKernel/SystemOfUnits.h"
 
#include <Gaudi/Accumulators.h>

class egammaTopoClusterCopier : public AthReentrantAlgorithm {
 
public:
 
  egammaTopoClusterCopier(const std::string& name, ISvcLocator* pSvcLocator);
 
  virtual StatusCode initialize() override final;
  virtual StatusCode execute(const EventContext& ctx) const override final;
  virtual StatusCode finalize() override final;
 
private:
 
  SG::ReadHandleKey<xAOD::CaloClusterContainer> m_inputTopoCollection {
    this,
    "InputTopoCollection",
    "CaloTopoClusters",
    "input topocluster collection"
  };
 
  SG::WriteHandleKey<xAOD::CaloClusterContainer> m_outputTopoCollectionShallow {
    this,
    "TopoShallow_doNotConfig",
    "",
    "Shallow copy of input collection that allows properties to be modified"
  };
 
  SG::WriteHandleKey<ConstDataVector<xAOD::CaloClusterContainer>> m_outputTopoCollection {
    this,
    "OutputTopoCollection",
    "egammaTopoCluster",
    "View container of selected topoclusters"
  };
 
  SG::WriteHandleKey<ConstDataVector<xAOD::CaloClusterContainer>> m_outputFwdTopoCollection {
    this,
    "OutputFwdTopoCollection",
    "",
    "View container of selected fwd topoclusters"
  };
 
  Gaudi::Property<float> m_etaCut {
    this,
    "EtaCut",
    2.6,
    "maximum |eta| of selected clusters"
  };
 
  Gaudi::Property<double> m_fwdEtaCut {
    this,
    "fwdEtaCut",
    2.5,
    "minimum |eta| of selected fwd clusters"
  };
 
  Gaudi::Property<float> m_ECut {
    this,
    "ECut",
    700,
    "minimum EM energy of selected clusters"
  };
 
  Gaudi::Property<double> m_fwdETCut {
    this,
    "fwdETCut",
    5. * Gaudi::Units::GeV,
    "Fwd ET cut"
  };
 
  Gaudi::Property<float> m_EMFracCut {
    this,
    "EMFracCut",
    0.5,
    "mimimum EM fraction of selected clusters"
  };

  Gaudi::Property<bool> m_hasITk {
    this,
    "hasITk",
    false,
    "Boolean to do track matching"
  };

  bool m_doForwardClusters = false;
 
  mutable Gaudi::Accumulators::Counter<> m_AllClusters {};
  mutable Gaudi::Accumulators::Counter<> m_CentralPassPreSelection {};
  mutable Gaudi::Accumulators::Counter<> m_CentralPassSelection {};
  mutable Gaudi::Accumulators::Counter<> m_FwdPassPreSelection {};
  mutable Gaudi::Accumulators::Counter<> m_FwdPassSelection {};
  mutable Gaudi::Accumulators::Counter<> m_SharedPassPreSelection {};
  mutable Gaudi::Accumulators::Counter<> m_SharedPassSelection {};
};
 
#endif // EGAMMATOOLS_EMCLUSTERTOOL_H