PhysicsAnalysis/TopPhys/xAOD/TopObjectSelectionTools/Root/Jet.cxx

Go to the documentation of this file.

/*
   Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
 */
 
#include "TopObjectSelectionTools/Jet.h"
 
#include <cmath>
 
#include "TopEvent/EventTools.h"
 
namespace top {
  Jet::Jet(const double ptcut,
           const double etamax,
           const double minmass,
           const bool doJVTCut):
    m_ptcut(ptcut),
    m_etamax(etamax),
    m_masscut(minmass),
    m_applyJVTCut(doJVTCut),
    m_appyMassCut(minmass > 0.),
    m_jvt_tool("JetJvtEfficiencyTool") {
    top::check(m_jvt_tool.retrieve(),
               "Failed to retrieve JVT tool");
  }
 
  // This version of the constructor always perform JVT cut 
  Jet::Jet(const double ptcut,
           const double etamax) : Jet::Jet(ptcut, etamax, -1.0, true) {}
 
  // DEPRECIATED - fwdJetSel string now defunct, keeping blank string input for backwards compatibility
  Jet::Jet(const double ptcut,
           const double etamax,
           const std::string) : Jet::Jet(ptcut, etamax, -1.0, true) {}
 
  // DEPRECATED - only kept for backward compatibility
  Jet::Jet(const double ptcut,
           const double etamax,
           const double) : Jet::Jet(ptcut, etamax) {}
 
  bool Jet::passSelection(const xAOD::Jet& jet) {
    if (m_applyJVTCut) {
      jet.auxdecor<char>("passJVT") = (m_jvt_tool->passesJvtCut(jet) ? 1 : 0);
    }
    //fJVT pass/fail decision already attatched in CP tools
 
    if (jet.pt() < m_ptcut) return false;
 
    if (std::fabs(jet.eta()) > m_etamax) return false;
 
    if (m_appyMassCut) {
      if (jet.m() < m_masscut) return false;
    }
 
    jet.auxdecor<char>("good") = 1;
    jet.auxdecor<char>("closeToLepton") = 0;
 
    return true;
  }
 
  void Jet::print(std::ostream& os) const {
    os << "Jet\n"
       << "    * pT > " << m_ptcut << "\n"
       << "    * |eta| < " << m_etamax << "\n";
    if (m_appyMassCut)
      os << "    * m > " << m_masscut << "\n";
  }
}  // namespace top