EntryLayerFilter.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
// class header include
#include "EntryLayerFilter.h"
 
// ISF includes
#include "ISF_Event/ISFParticle.h"
 
// ISF interfaces
#include "ISF_Interfaces/IParticleFilter.h"
 
ISF::EntryLayerFilter::EntryLayerFilter(const std::string& t, const std::string& n, const IInterface* p) :
  base_class(t,n,p),
  m_ekinCharged(-1),
  m_ekinNeutral(-1)
{
  declareProperty( "MinEkinCharged",
                   m_ekinCharged=-1.,
                   "Minimum kinetic energy cut for charged particles (set -1 to filter out all)");
  declareProperty( "MinEkinNeutral",
                   m_ekinNeutral=-1.,
                   "Minimum kinetic energy cut for neutral particles (set -1 to filter out all)");
}
 
 
bool ISF::EntryLayerFilter::passFilter(const ISFParticle &p) const {
  double        ekin = p.ekin();
  bool     isCharged = fabs(p.charge()) > 1e-5;
 
  double  curEkinCut = isCharged ? m_ekinCharged : m_ekinNeutral;
 
  // pass if cut is activated and energy is above threshold
  bool pass = ( curEkinCut > 0. ) && ( ekin > curEkinCut );
 
  return pass;
}