GenericTruthStrategy.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef ISF_TOOLS_GENERICTRUTHSTRATEGY_H
#define ISF_TOOLS_GENERICTRUTHSTRATEGY_H 1
 
// stl includes
#include <set>
#include <vector>
 
// Athena includes
#include "AthenaBaseComps/AthAlgTool.h"
#include "AtlasDetDescr/AtlasRegion.h"
 
// ISF includes
#include "ISF_HepMC_Interfaces/ITruthStrategy.h"
 
namespace ISF {
  typedef std::vector<int>              VertexTypesVector;
  typedef std::set<int>                 VertexTypesSet;
  typedef std::vector<int>              PDGCodesVector;
  typedef std::set<int>                 PDGCodesSet;

  class GenericTruthStrategy final : public extends<AthAlgTool, ITruthStrategy> {
 
  public:
    GenericTruthStrategy( const std::string& t, const std::string& n, const IInterface* p );

    ~GenericTruthStrategy() = default;
 
    // Athena algtool's Hooks
    virtual StatusCode  initialize() override final;

    virtual bool pass( ITruthIncident& incident) const override final;
 
    virtual bool appliesToRegion(unsigned short geoID) const override final;
  private:
    // provide either a pT or Ekin cut for the parent and child particles respectively.
    // if none are given for either type, it will not use pT or Ekin cuts
    // (the Pt2 variables get squared in the initialize() method)

    bool m_useParentPt{true};         
    double m_parentPt2{-1.};
    Gaudi::Property<double> m_parentPt{this, "ParentMinPt", -1.};           
    Gaudi::Property<double> m_parentEkin{this, "ParentMinEkin", -1.};          

    bool m_useChildPt{true};          
    double m_childPt2{-1.};
    Gaudi::Property<double> m_childPt{this, "ChildMinPt", -1.};            
    Gaudi::Property<double> m_childEkin{this, "ChildMinEkin", -1.};           
    Gaudi::Property<bool> m_allowChildrenOrParentPass{this, "AllowChildrenOrParentPassKineticCuts", false}; 
    // if set to true, kinetic cuts are passed even if only child particles pass them
    // (used for special cases such as de-excitation)

    Gaudi::Property<VertexTypesVector> m_vertexTypesVector{this, "VertexTypes", 0};  
    VertexTypesSet m_vertexTypes{};        
    bool m_doVertexRangeCheck{false};
    Gaudi::Property<int> m_vertexTypeRangeLow{this, "VertexTypeRangeLow", 0};
    Gaudi::Property<int> m_vertexTypeRangeHigh{this, "VertexTypeRangeHigh", 0};
    unsigned m_vertexTypeRangeLength{0};

    Gaudi::Property<PDGCodesVector> m_parentPdgCodesVector{this, "ParentPDGCodes", 0};  
    PDGCodesSet m_parentPdgCodes{};        
 
    IntegerArrayProperty            m_regionListProperty{this, "Regions", {}};
  };
 
}
 
 
#endif //> !ISF_TOOLS_GENERICTRUTHSTRATEGY_H