PseudoJetAlgorithm.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
// PseudoJetAlgorithm.h 

 
#ifndef PseudoJetAlgorithm_H
#define PseudoJetAlgorithm_H
 
#include <memory>
#include "fastjet/PseudoJet.hh"
#include "JetRec/PseudoJetContainer.h"
 
#include "xAODTracking/VertexContainer.h"
#include "xAODTracking/Vertex.h"
 
#include "AsgDataHandles/ReadHandleKey.h"
#include "AsgDataHandles/WriteHandleKey.h"
#include "AsgTools/PropertyWrapper.h"
#include "AnaAlgorithm/AnaReentrantAlgorithm.h"
 
class PseudoJetAlgorithm : public EL::AnaReentrantAlgorithm { 
 
public: 
 
  // Can't use "using ctor" because of incompatiblity with pyroot in AnalysisBase
  PseudoJetAlgorithm(const std::string & n, ISvcLocator* l) : EL::AnaReentrantAlgorithm(n,l) {}

  virtual StatusCode  initialize() override final;
 
  // Standard execute, forwards to createAndRecord
  virtual StatusCode  execute(const EventContext& ctx) const override final;
 
private: 

  virtual std::unique_ptr<PseudoJetContainer> createPJContainer(const xAOD::IParticleContainer& cont, const EventContext& ctx) const;

  virtual void print() const;
 
  std::vector<fastjet::PseudoJet> 
  createPseudoJets(const xAOD::IParticleContainer&) const;
 
  std::vector<fastjet::PseudoJet> 
  createPseudoJets(const xAOD::IParticleContainer&, const xAOD::VertexContainer* vertices) const;
 
private:

  SG::ReadHandleKey<xAOD::IParticleContainer> m_incoll{this, "InputContainer", "", "The input IParticleContainer name"};

  SG::WriteHandleKey<PseudoJetContainer> m_outcoll{this, "OutputContainer", "", "The output PseudoJetContainer name"};

  Gaudi::Property<std::string> m_label{this, "Label", "", "String label identifying the pseudojet type"};

  Gaudi::Property<bool> m_skipNegativeEnergy{this, "SkipNegativeEnergy", false, "Whether to skip negative energy inputs"};

  Gaudi::Property<bool> m_negEnergyAsGhosts{this, "TreatNegativeEnergyAsGhost", false, "Whether to convert negative energy inputs into ghosts"};

  Gaudi::Property<bool> m_useCharged{this, "UseCharged", true, "Whether to use charged PFOs/FEs"};

  Gaudi::Property<bool> m_useNeutral{this, "UseNeutral", true, "Whether to use neutral PFOs/FEs"};

  Gaudi::Property<bool> m_useChargedPV{this, "UseChargedPV", true, "Whether to use charged PFOs/FEs matched to the PV"};

  Gaudi::Property<bool> m_useChargedPUsideband{this, "UseChargedPUsideband", false, "Whether to use charged PU sideband only"};

  Gaudi::Property<bool> m_byVertex{this, "DoByVertex", false, "True if jets should be reconstructed by vertex"};
 
  SG::ReadHandleKey<xAOD::VertexContainer> m_vertexContainer_key{this, "VertexContainer", "PrimaryVertices", "Vertex container (for by-vertex reconstruction)"};

  bool m_isGhost{false}; 
  bool m_emtopo{false};  
  bool m_pflow{false};   
  bool m_ufo{false};     
}; 
 
#endif