ObjectCutFlowHistAlg.h

Go to the documentation of this file.

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

 
 
#ifndef ASG_ANALYSIS_ALGORITHMS__OBJECT_CUT_FLOW_HIST_ALG_H
#define ASG_ANALYSIS_ALGORITHMS__OBJECT_CUT_FLOW_HIST_ALG_H
 
#include <AnaAlgorithm/AnaAlgorithm.h>
#include <AsgTools/PropertyWrapper.h>
#include <SelectionHelpers/SysReadSelectionHandle.h>
#include <SystematicsHandles/SysReadHandle.h>
#include <SystematicsHandles/SysListHandle.h>
#include <xAODBase/IParticleContainer.h>
 
namespace CP
{
 
  class ObjectCutFlowHistAlg final : public EL::AnaAlgorithm
  {
  public:
    using EL::AnaAlgorithm::AnaAlgorithm;
    StatusCode initialize () override;
    StatusCode execute (const EventContext& ctx) override;
 

  private:
    SysListHandle m_systematicsList {this};

  private:
    SysReadHandle<xAOD::IParticleContainer> m_inputHandle {
      this, "input", "", "the input particle container to run on"};

  private:
    SysReadSelectionHandle m_preselection {
      this, "preselection", "", "the preselection to apply"};

  private:
    Gaudi::Property<std::string> m_histPattern {this, "histPattern", "cutflow_%SYS%", "the pattern for histogram names"};

  private:
    Gaudi::Property<std::string> m_histTitle {this, "histTitle", "object cut flow", "title for the created histograms"};

  private:
    SysReadSelectionHandleArray m_selections {
      this, "selections", {}, "the inputs to the object cutflow"};

  private:
    Gaudi::Property<bool> m_forceCutSequence {this, "forceCutSequence", false, "force cut sequence and do not accept if previous cuts failed"};

  private:
    unsigned m_allCutsNum = 0;

  private:
    std::unordered_map<CP::SystematicSet,TH1*> m_hist;

  private:
    std::vector<std::string> m_labels;
  };
}
 
#endif