HIJetDiscriminatorTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#include "HIJetDiscriminatorTool.h"
 
namespace{
  struct  Discrim{
    Discrim(float maxOverMean, float ETMin) :  _maxOverMean(maxOverMean),
                           _ETMin(ETMin){};
 
    float _maxOverMean;
    float _ETMin;
 
    //return true if jet should be removed
    bool operator()(const xAOD::Jet* j1)
    {
      float max_ET=j1->getAttribute<float>("MaxConstituentET");
      float disc=j1->getAttribute<float>("MaxOverMean");
      if( max_ET < _ETMin ) return true;
      if( disc < _maxOverMean ) return true;
      return false;
    }
  };
} //annonymous namespace
 
HIJetDiscriminatorTool::HIJetDiscriminatorTool(const std::string& t) : JetModifierBase(t)
{
}
 
StatusCode HIJetDiscriminatorTool::modify(xAOD::JetContainer& jets) const
{
  xAOD::JetContainer::iterator itB = jets.begin();
  xAOD::JetContainer::iterator itE = jets.end();
  xAOD::JetContainer::iterator lastFiltered=std::remove_if(itB, itE, Discrim(m_MaxOverMeanCut,m_ETminCut));
  jets.erase( lastFiltered, itE );
  return StatusCode::SUCCESS;
}