JetUncertaintiesAlg.h

Go to the documentation of this file.

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

 
 
#ifndef JET_ANALYSIS_ALGORITHMS__JET_UNCERTAINTIES_ALG_H
#define JET_ANALYSIS_ALGORITHMS__JET_UNCERTAINTIES_ALG_H
 
#include <AnaAlgorithm/AnaAlgorithm.h>
#include <CxxUtils/checker_macros.h>
#include <JetCPInterfaces/ICPJetUncertaintiesTool.h>
#include <SelectionHelpers/OutOfValidityHelper.h>
#include <SelectionHelpers/SysReadSelectionHandle.h>
#include <SystematicsHandles/SysCopyHandle.h>
#include <SystematicsHandles/SysListHandle.h>
 
#include <optional>
 
namespace CP
{
 
  class ATLAS_NOT_THREAD_SAFE JetUncertaintiesAlg final
    : public EL::AnaAlgorithm
  {
  public:
    using EL::AnaAlgorithm::AnaAlgorithm;
    StatusCode initialize () override;
    StatusCode execute (const EventContext& ctx) override;

  private:
    StatusCode applyTool (ToolHandle<ICPJetUncertaintiesTool>& tool,
                          const std::vector<CP::SystematicSet>& systematics,
                          bool decorateBtag, const EventContext& ctx);
 

  private:
    ToolHandle<ICPJetUncertaintiesTool> m_uncertaintiesTool {this, "uncertaintiesTool", "JetUncertaintiesTool", "the uncertainties tool we apply"};

  private:
    ToolHandle<ICPJetUncertaintiesTool> m_uncertaintiesToolPD {this, "uncertaintiesToolPD", "", "the uncertainties tool we apply specifically for the 'Full'/'All' JER systematic models"};

  private:
    SysListHandle m_systematicsList {this};

  private:
    SysCopyHandle<xAOD::JetContainer> m_jetHandle {
      this, "jets", "AntiKt4EMTopoJets", "the jet collection to run on"};

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

  private:
    OutOfValidityHelper m_outOfValidity {this};

  private:
    std::vector<CP::SystematicSet> m_systematicsVector;

  private:
    std::vector<CP::SystematicSet> m_systematicsVectorOnlyJERPseudoData;
 
  private:
    Gaudi::Property<std::string> m_isJESbtag {this, "isJESbtagLabel", "IsBjet", "The label to apply to truth b-tagged jets, for JES flavour uncertainties"};
    std::optional<SG::Decorator<char>> m_decIsJESbtag;
    std::optional<SG::ConstAccessor<int>> m_accTruthLabel;
 
  };
}
 
#endif