ScaleFactorRetriever.h

Go to the documentation of this file.

/*
   Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
 */
 
#ifndef SCALEFACTORRETRIEVER_H_
#define SCALEFACTORRETRIEVER_H_
 
#include <memory>
#include <string>
#include <vector>
 
#include "AsgTools/AsgTool.h"
#include "xAODEgamma/ElectronContainer.h"
#include "xAODEgamma/PhotonContainer.h"
#include "xAODMuon/MuonContainer.h"
#include "xAODTau/TauJetContainer.h"
#include "AthContainers/AuxElement.h"
 
// Forward declare:
namespace top {
  class Event;
  class TopConfig;
}
 
 
namespace top {
  enum topSFSyst {
    nominal = 0,
    // Electron Trigger SFs
    EL_SF_Trigger_UP, EL_SF_Trigger_DOWN,
    // Electron reconstruction SFs
    EL_SF_Reco_UP, EL_SF_Reco_DOWN,
    // Electron ID SFs
    EL_SF_ID_UP, EL_SF_ID_DOWN,
    // Electron isolation SFs
    EL_SF_Isol_UP, EL_SF_Isol_DOWN,
    // Electron Charge flip SFs
    EL_SF_ChargeID_UP, EL_SF_ChargeID_DOWN,
    EL_SF_ChargeMisID_STAT_UP, EL_SF_ChargeMisID_STAT_DOWN,
    EL_SF_ChargeMisID_SYST_UP, EL_SF_ChargeMisID_SYST_DOWN,
    // Muon Trigger SFs
    MU_SF_Trigger_UP, MU_SF_Trigger_DOWN,
    MU_SF_Trigger_STAT_UP, MU_SF_Trigger_STAT_DOWN,
    MU_SF_Trigger_SYST_UP, MU_SF_Trigger_SYST_DOWN,
    // Muon ID SFs (regular)
    MU_SF_ID_STAT_UP, MU_SF_ID_STAT_DOWN,
    MU_SF_ID_SYST_UP, MU_SF_ID_SYST_DOWN,
    // Muon ID SFs (low pT)
    MU_SF_ID_STAT_LOWPT_UP, MU_SF_ID_STAT_LOWPT_DOWN,
    MU_SF_ID_SYST_LOWPT_UP, MU_SF_ID_SYST_LOWPT_DOWN,
    // Muon isolation SFs
    MU_SF_Isol_STAT_UP, MU_SF_Isol_STAT_DOWN,
    MU_SF_Isol_SYST_UP, MU_SF_Isol_SYST_DOWN,
    // Muon TTVA SFs
    MU_SF_TTVA_STAT_UP, MU_SF_TTVA_STAT_DOWN,
    MU_SF_TTVA_SYST_UP, MU_SF_TTVA_SYST_DOWN,
    // For more complex systematics
    MU_SF_ID_BKG_FRACTION_UP, MU_SF_ID_BKG_FRACTION_DOWN,
    MU_SF_ID_FIT_MODEL_LOWPT_UP, MU_SF_ID_FIT_MODEL_LOWPT_DOWN,
    MU_SF_ID_LUMI_UNCERT_UP, MU_SF_ID_LUMI_UNCERT_DOWN,
    MU_SF_ID_MATCHING_UP, MU_SF_ID_MATCHING_DOWN,
    MU_SF_ID_MATCHING_LOWPT_UP, MU_SF_ID_MATCHING_LOWPT_DOWN,
    MU_SF_ID_MC_XSEC_UP, MU_SF_ID_MC_XSEC_DOWN,
    MU_SF_ID_PT_DEPENDENCY_UP, MU_SF_ID_PT_DEPENDENCY_DOWN,
    MU_SF_ID_QCD_TEMPLATE_UP, MU_SF_ID_QCD_TEMPLATE_DOWN,
    MU_SF_ID_SUPRESSION_SCALE_UP, MU_SF_ID_SUPRESSION_SCALE_DOWN,
    MU_SF_ID_TRUTH_UP, MU_SF_ID_TRUTH_DOWN,
    MU_SF_ID_TRUTH_LOWPT_UP, MU_SF_ID_TRUTH_LOWPT_DOWN,
    MU_SF_ID_BAD_MUON_VETO_UP, MU_SF_ID_BAD_MUON_VETO_DOWN,
    MU_SF_ID_CR1_UP, MU_SF_ID_CR1_DOWN,
    MU_SF_ID_CR2_UP, MU_SF_ID_CR2_DOWN,
    MU_SF_ID_CR3_UP, MU_SF_ID_CR3_DOWN,
    MU_SF_ID_HIGHETA_PROBEIP_UP, MU_SF_ID_HIGHETA_PROBEIP_DOWN,
    MU_SF_ID_HIGHETA_PROBEISO_UP, MU_SF_ID_HIGHETA_PROBEISO_DOWN,
    MU_SF_ID_TAGPT_UP, MU_SF_ID_TAGPT_DOWN,
    MU_SF_ID_EXTRAPOLATION_UP, MU_SF_ID_EXTRAPOLATION_DOWN,
    MU_SF_Isol_BKG_FRACTION_UP, MU_SF_Isol_BKG_FRACTION_DOWN,
    MU_SF_Isol_DRMUJ_UP, MU_SF_Isol_DRMUJ_DOWN,
    MU_SF_Isol_LUMI_UNCERT_UP, MU_SF_Isol_LUMI_UNCERT_DOWN,
    MU_SF_Isol_MC_XSEC_UP, MU_SF_Isol_MC_XSEC_DOWN,
    MU_SF_Isol_MLLWINDOW_UP, MU_SF_Isol_MLLWINDOW_DOWN,
    MU_SF_Isol_QCD_TEMPLATE_UP, MU_SF_Isol_QCD_TEMPLATE_DOWN,
    MU_SF_Isol_SHERPA_POWHEG_UP, MU_SF_Isol_SHERPA_POWHEG_DOWN,
    MU_SF_Isol_SUPRESSION_SCALE_UP, MU_SF_Isol_SUPRESSION_SCALE_DOWN,
    MU_SF_Isol_EXTRAPOLATION_UP, MU_SF_Isol_EXTRAPOLATION_DOWN,
    MU_SF_TTVA_BKG_FRACTION_UP, MU_SF_TTVA_BKG_FRACTION_DOWN,
    MU_SF_TTVA_LUMI_UNCERT_UP, MU_SF_TTVA_LUMI_UNCERT_DOWN,
    MU_SF_TTVA_MC_XSEC_UP, MU_SF_TTVA_MC_XSEC_DOWN,
    MU_SF_TTVA_QCD_TEMPLATE_UP, MU_SF_TTVA_QCD_TEMPLATE_DOWN,
    MU_SF_TTVA_SUPRESSION_SCALE_UP, MU_SF_TTVA_SUPRESSION_SCALE_DOWN,
    MU_SF_TTVA_EXTRAPOLATION_UP, MU_SF_TTVA_EXTRAPOLATION_DOWN,
 
    TAU_SF_NOMINAL,
    // Tau EleOLR SFs
    TAUS_TRUEHADTAU_EFF_ELEOLR_TOTAL__1up, TAUS_TRUEHADTAU_EFF_ELEOLR_TOTAL__1down,
    TAUS_TRUEELECTRON_EFF_ELEOLR_STAT__1up, TAUS_TRUEELECTRON_EFF_ELEOLR_STAT__1down,
    TAUS_TRUEELECTRON_EFF_ELEOLR_SYST__1up, TAUS_TRUEELECTRON_EFF_ELEOLR_SYST__1down,
    // Tau RNNID SFs
    TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2025__1up, TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2025__1down,
    TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2530__1up, TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2530__1down,
    TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT3040__1up, TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT3040__1down,
    TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPTGE40__1up, TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPTGE40__1down,
      TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2025__1up, TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2025__1down,
      TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2530__1up, TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2530__1down,
      TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT3040__1up, TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT3040__1down,
      TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPTGE40__1up, TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPTGE40__1down,
    TAUS_TRUEHADTAU_EFF_RNNID_HIGHPT__1up, TAUS_TRUEHADTAU_EFF_RNNID_HIGHPT__1down,
    TAUS_TRUEHADTAU_EFF_RNNID_SYST__1up, TAUS_TRUEHADTAU_EFF_RNNID_SYST__1down,
    // Tau Reconstruction SFs
    TAUS_TRUEHADTAU_EFF_RECO_TOTAL__1up, TAUS_TRUEHADTAU_EFF_RECO_TOTAL__1down,
   TAU_SF_END,
    // Photon SFs
    PHOTON_IDSF_UP, PHOTON_IDSF_DOWN,
    PHOTON_EFF_ISO_UP, PHOTON_EFF_ISO_DOWN,
    PHOTON_EFF_TRIGGER_UNCERTAINTY_UP,
    PHOTON_EFF_TRIGGER_UNCERTAINTY_DOWN,
    // B-tagging SFs
    BTAG_SF_EIGEN_B, BTAG_SF_EIGEN_C,
    BTAG_SF_EIGEN_LIGHT,
    BTAG_SF_NAMED_UP, BTAG_SF_NAMED_DOWN,
    // JVT SFs
    JVT_UP, JVT_DOWN,
    // fJVT SFs
    FJVT_UP, FJVT_DOWN,
    // FWD Electron ID SFs
    FWDEL_SF_ID_UP, FWDEL_SF_ID_DOWN,
 
    lastsys
  };
  const std::map<topSFSyst, TString> tauSF_name = {
    {TAU_SF_NOMINAL, "TAU_SF_NOMINAL"},
    {TAUS_TRUEHADTAU_EFF_RECO_TOTAL__1up, "TAUS_TRUEHADTAU_EFF_RECO_TOTAL__1up"},
    {TAUS_TRUEHADTAU_EFF_RECO_TOTAL__1down, "TAUS_TRUEHADTAU_EFF_RECO_TOTAL__1down"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2025__1up,
     "TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2025__1up"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2025__1down,
     "TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2025__1down"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2530__1up,
     "TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2530__1up"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2530__1down,
     "TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2530__1down"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT3040__1up,
     "TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT3040__1up"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT3040__1down,
      "TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT3040__1down"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPTGE40__1up,
       "TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPTGE40__1up"},
      {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPTGE40__1down,
       "TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPTGE40__1down"},
    {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2025__1up,
     "TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2025__1up"},
    {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2025__1down,
     "TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2025__1down"},
    {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2530__1up,
     "TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2530__1up"},
    {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2530__1down,
     "TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2530__1down"},
    {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT3040__1up,
     "TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT3040__1up"},
    {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT3040__1down,
      "TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT3040__1down"},
    {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPTGE40__1up,
       "TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPTGE40__1up"},
      {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPTGE40__1down,
       "TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPTGE40__1down"},
    {TAUS_TRUEHADTAU_EFF_RNNID_HIGHPT__1up, "TAUS_TRUEHADTAU_EFF_RNNID_HIGHPT__1up"},
    {TAUS_TRUEHADTAU_EFF_RNNID_HIGHPT__1down, "TAUS_TRUEHADTAU_EFF_RNNID_HIGHPT__1down"},
    {TAUS_TRUEHADTAU_EFF_RNNID_SYST__1up, "TAUS_TRUEHADTAU_EFF_RNNID_SYST__1up"},
    {TAUS_TRUEHADTAU_EFF_RNNID_SYST__1down, "TAUS_TRUEHADTAU_EFF_RNNID_SYST__1down"},
    {TAUS_TRUEHADTAU_EFF_ELEOLR_TOTAL__1up, "TAUS_TRUEHADTAU_EFF_ELEOLR_TOTAL__1up"},
    {TAUS_TRUEHADTAU_EFF_ELEOLR_TOTAL__1down, "TAUS_TRUEHADTAU_EFF_ELEOLR_TOTAL__1down"},
    {TAUS_TRUEELECTRON_EFF_ELEOLR_STAT__1up, "TAUS_TRUEELECTRON_EFF_ELEOLR_STAT__1up"},
    {TAUS_TRUEELECTRON_EFF_ELEOLR_STAT__1down, "TAUS_TRUEELECTRON_EFF_ELEOLR_STAT__1down"},
    {TAUS_TRUEELECTRON_EFF_ELEOLR_SYST__1up, "TAUS_TRUEELECTRON_EFF_ELEOLR_SYST__1up"},
    {TAUS_TRUEELECTRON_EFF_ELEOLR_SYST__1down, "TAUS_TRUEELECTRON_EFF_ELEOLR_SYST__1down"},
  };
 
  const std::map<topSFSyst, TString> tauSF_alias = {
    {TAU_SF_NOMINAL, "TAU_SF_NOMINAL"},
    {TAUS_TRUEHADTAU_EFF_RECO_TOTAL__1up, "EFF_RECO_TOTAL_UP"},
    {TAUS_TRUEHADTAU_EFF_RECO_TOTAL__1down, "EFF_RECO_TOTAL_DOWN"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2025__1up, "EFF_RNNID_1PRONGSTATSYSTPT2025_UP"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2025__1down, "EFF_RNNID_1PRONGSTATSYSTPT2025_DOWN"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2530__1up, "EFF_RNNID_1PRONGSTATSYSTPT2530_UP"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT2530__1down, "EFF_RNNID_1PRONGSTATSYSTPT2530_DOWN"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT3040__1up, "EFF_RNNID_1PRONGSTATSYSTPT3040_UP"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPT3040__1down, "EFF_RNNID_1PRONGSTATSYSTPT3040_DOWN"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPTGE40__1up, "EFF_RNNID_1PRONGSTATSYSTPTGE40_UP"},
    {TAUS_TRUEHADTAU_EFF_RNNID_1PRONGSTATSYSTPTGE40__1down, "EFF_RNNID_1PRONGSTATSYSTPTGE40_DOWN"},
      {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2025__1up, "EFF_RNNID_3PRONGSTATSYSTPT2025_UP"},
      {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2025__1down, "EFF_RNNID_3PRONGSTATSYSTPT2025_DOWN"},
      {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2530__1up, "EFF_RNNID_3PRONGSTATSYSTPT2530_UP"},
      {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT2530__1down, "EFF_RNNID_3PRONGSTATSYSTPT2530_DOWN"},
      {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT3040__1up, "EFF_RNNID_3PRONGSTATSYSTPT3040_UP"},
      {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPT3040__1down, "EFF_RNNID_3PRONGSTATSYSTPT3040_DOWN"},
      {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPTGE40__1up, "EFF_RNNID_3PRONGSTATSYSTPTGE40_UP"},
      {TAUS_TRUEHADTAU_EFF_RNNID_3PRONGSTATSYSTPTGE40__1down, "EFF_RNNID_3PRONGSTATSYSTPTGE40_DOWN"},
    {TAUS_TRUEHADTAU_EFF_RNNID_HIGHPT__1up, "EFF_RNNID_HIGHPT_UP"},
    {TAUS_TRUEHADTAU_EFF_RNNID_HIGHPT__1down, "EFF_RNNID_HIGHPT_DOWN"},
    {TAUS_TRUEHADTAU_EFF_RNNID_SYST__1up, "EFF_RNNID_SYST_UP"},
    {TAUS_TRUEHADTAU_EFF_RNNID_SYST__1down, "EFF_RNNID_SYST_DOWN"},
    {TAUS_TRUEHADTAU_EFF_ELEOLR_TOTAL__1up, "EFF_ELEOLR_TOTAL_UP"},
    {TAUS_TRUEHADTAU_EFF_ELEOLR_TOTAL__1down, "EFF_ELEOLR_TOTAL_DOWN"},
    {TAUS_TRUEELECTRON_EFF_ELEOLR_STAT__1up, "ELE_EFF_ELEOLR_STAT_UP"},
    {TAUS_TRUEELECTRON_EFF_ELEOLR_STAT__1down, "ELE_EFF_ELEOLR_STAT_DOWN"},
    {TAUS_TRUEELECTRON_EFF_ELEOLR_SYST__1up, "ELE_EFF_ELEOLR_SYST_UP"},
    {TAUS_TRUEELECTRON_EFF_ELEOLR_SYST__1down, "ELE_EFF_ELEOLR_SYST_DOWN"},
  };
 
 
  enum topSFComp {
    ALL = 0, TRIGGER, RECO, ID, ISOLATION, CHARGEID, CHARGEMISID, TTVA
  };
 
  class ScaleFactorRetriever final: public asg::AsgTool {
    // Specify that we will overload print with top::Event later on
    using asg::AsgTool::print;
  public:
    explicit ScaleFactorRetriever(const std::string& name);
    virtual ~ScaleFactorRetriever() {}
 
    // Delete Standard constructors
    ScaleFactorRetriever(const ScaleFactorRetriever& rhs) = delete;
    ScaleFactorRetriever(ScaleFactorRetriever&& rhs) = delete;
    ScaleFactorRetriever& operator = (const ScaleFactorRetriever& rhs) = delete;
 
    StatusCode initialize();
 
    // Pile up SF
    static bool hasPileupSF(const top::Event& event);
    static float pileupSF(const top::Event& event, int var = 0);
 
    // Obtain the lepton SF
    float leptonSF(const top::Event& event, const top::topSFSyst SFSyst) const;
 
    // Obtain the trigger SF
    float triggerSF(const top::Event& event, const top::topSFSyst SFSyst) const;
 
    // Obtain trigger SF via global tool
    float globalTriggerSF(const top::Event& event, const top::topSFSyst SFSyst) const;
 
    // Obtain the trigger SF via approximate formula
    float oldTriggerSF(const top::Event& event, const top::topSFSyst SFSyst) const;
 
    // Obtain the electron SF
    float electronSF(const top::Event& event,
                     const top::topSFSyst SFSyst,
                     const top::topSFComp SFComp) const;
 
    //get vector of systematic variations of electron SFs, var=1 for up, var=-1 for down
    std::vector<float> electronSFSystVariationVector(const top::Event& event, const top::topSFComp SFComp, int var = 1) const;
 
    //get vector of systematic variations of electron SFs, var=1 for up, var=-1 for down and size of systematic variations
    std::vector<float> electronSFSystVariationVector(const top::Event& event, const top::topSFComp SFComp, int var = 1, int size=50) const;
 
    // Obtain the fwd electron SF
    float fwdElectronSF(const top::Event& event,
                        const top::topSFSyst SFSyst,
                        const top::topSFComp SFComp = top::topSFComp::ALL) const;
 
    // Obtain the muon SF
    float muonSF(const top::Event& event,
                 const top::topSFSyst SFSyst,
                 const top::topSFComp SFComp) const;
 
    // Obtain the soft muon SF
    float softmuonSF(const top::Event& event,
                     const top::topSFSyst SFSyst,
                     const top::topSFComp SFComp) const;
 
    float tauSF(const top::Event& event,
                const top::topSFSyst SFSyst) const;
 
    float photonSF(const top::Event& event,
                   const top::topSFSyst SFSyst) const;
 
    float btagSF(const top::Event& event,
                 const top::topSFSyst SFSyst = top::topSFSyst::nominal,
                 std::string WP = "FixedCutBEff_77",
                 bool do_trackjets = false,
                 std::string uncert_name = "") const;
 
    void btagSF_eigen_vars(const top::Event& event,
                           const top::topSFSyst SFSyst,
                           std::vector<float>& btagSF_up,
                           std::vector<float>& btagSF_down,
                           std::string WP = "FixedCutBEff_77",
                           bool do_trackjets = false) const;
 
    float jvtSF(const top::Event& event,
                const top::topSFSyst SFsyst) const;
 
    float fjvtSF(const top::Event& event,
                const top::topSFSyst SFsyst) const;
 
    void print(const top::Event& event);
 
    float electronSF_Reco(const xAOD::Electron& x,
                          const top::topSFSyst SFSyst
              ) const;
 
    float electronSF_Trigger(const xAOD::Electron& x, const top::topSFSyst SFSyst, bool useLooseDef) const;
 
    float electronEff_Trigger(const xAOD::Electron& x, const top::topSFSyst SFSyst, bool useLooseDef) const;
 
    float electronSF_ID(const xAOD::Electron& x,
                        const top::topSFSyst SFSyst,
                        bool useLooseDef) const;
 
    float electronSF_Isol(const xAOD::Electron& x,
                          const top::topSFSyst SFSyst,
                          bool useLooseDef) const;
 
    float electronSF_ChargeID(const xAOD::Electron& x,
                              const top::topSFSyst SFSyst,
                              bool useLooseDef) const;
 
    float electronSF_ChargeMisID(const xAOD::Electron& x,
                                 const top::topSFSyst SFSyst,
                                 bool useLooseDef) const;
 
    float fwdElectronSF_ID(const xAOD::Electron& x,
                           const top::topSFSyst SFSyst,
                           bool useLooseDef) const;
 
    float muonSF_Trigger(const xAOD::Muon& x,
                         const top::topSFSyst SFSyst,
                         bool useLooseDef) const;
 
    float muonEff_Trigger(const xAOD::Muon& x,
                          const top::topSFSyst SFSyst,
                          bool useLooseDef) const;
 
    float muonSF_ID(const xAOD::Muon& x,
                    const top::topSFSyst SFSyst,
                    bool useLooseDef) const;
 
    float muonSF_Isol(const xAOD::Muon& x,
                      const top::topSFSyst SFSyst,
                      bool useLooseDef) const;
 
    float muonSF_TTVA(const xAOD::Muon& x,
                      const top::topSFSyst SFSyst) const;
 
    float softmuonSF_ID(const xAOD::Muon& x,
                        const top::topSFSyst SFSyst) const;
 
    float tauSF(const xAOD::TauJet& x,
                const top::topSFSyst SFSyst,
                bool isLoose) const;
 
    float photonSF_Reco(const xAOD::Photon& x,
            const top::topSFSyst SFSyst) const;
 
    float photonSF_Isol(const xAOD::Photon& x,
                        const top::topSFSyst SFSyst,
                        bool isLoose) const;
  private:
    float electronSF_Trigger(const xAOD::Electron& x,
                             const std::string& id,
                             const top::topSFSyst SFSyst,
                             bool useLooseDef) const;
 
    float electronEff_Trigger(const xAOD::Electron& x,
                              const std::string& id,
                              const top::topSFSyst SFSyst,
                              bool useLooseDef) const;
 
    float electronSF_ID(const xAOD::Electron& x,
                        const std::string& id,
                        const top::topSFSyst SFSyst,
                        bool useLooseDef) const;
 
    float electronSF_Isol(const xAOD::Electron& x,
                          const std::string& id,
                          const top::topSFSyst SFSyst,
                          bool useLooseDef) const;
 
    float electronSF_ChargeID(const xAOD::Electron& x,
                              const std::string& id, const std::string& iso,
                              const top::topSFSyst SFSyst,
                              bool useLooseDef) const;
 
    float electronSF_ChargeMisID(const xAOD::Electron& x,
                                 const std::string& id, const std::string& iso,
                                 const top::topSFSyst SFSyst,
                                 bool useLooseDef) const;
 
    float fwdElectronSF_ID(const xAOD::Electron& x,
                           const std::string& id,
                           const top::topSFSyst SFSyst,
                           bool useLooseDef) const;
 
    float photonSF_Isol(const xAOD::Photon& x,
                        const std::string& iso,
                        const top::topSFSyst SFSyst,
                        bool isLoose) const;
 
    float muonSF_Trigger(const xAOD::Muon& x,
                         const std::string& id,
                         const top::topSFSyst SFSyst,
                         bool useLooseDef) const;
    float muonEff_Trigger(const xAOD::Muon& x,
                          const std::string& id,
                          const top::topSFSyst SFSyst,
                          bool useLooseDef) const;
 
    float muonSF_ID(const xAOD::Muon& x,
                    const std::string& id,
                    const top::topSFSyst SFSyst,
                    bool useLooseDef) const;
 
    float softmuonSF_ID(const xAOD::Muon& x,
                        const std::string& id,
                        const top::topSFSyst SFSyst) const;
 
    float muonSF_Isol(const xAOD::Muon& x,
                      const std::string& iso,
                      const top::topSFSyst SFSyst,
                      bool useLooseDef) const;
 
    // List of triggers to 'or' together for each event.
    // If any one passes, the event passes
    std::vector<std::pair<std::string, int> > m_electronTriggers_Tight;
    std::vector<std::pair<std::string, int> > m_muonTriggers_Tight;
    std::vector<std::pair<std::string, int> > m_photonTriggers_Tight;
    std::vector<std::pair<std::string, int> > m_electronTriggers_Loose;
    std::vector<std::pair<std::string, int> > m_muonTriggers_Loose;
    std::vector<std::pair<std::string, int> > m_photonTriggers_Loose;
    // Do we need to add tau triggers?
 
    // Configuration
    std::shared_ptr<top::TopConfig> m_config;
    bool m_preferGlobalTriggerSF;
 
    static size_t s_warn_counter; // print WARNING messages limited times to prevent excessive log spam
  };
}  // namespace top
 
#endif  // SCALEFACTORRETRIEVER_H_