GenericTrackSelector.cxx

Go to the documentation of this file.

 
#include "xAODEgamma/ElectronContainer.h"
 
#include "GenericTrackSelector.h"
 
#include <cmath>
 
 
IDTPM2::GenericTrackSelector::GenericTrackSelector( const std::string& name )
  : asg::AsgTool( name ) { }
 
 
 
StatusCode IDTPM2::GenericTrackSelector::initialize() {
  if ( m_initialized ) return StatusCode::SUCCESS;
  ATH_CHECK( asg::AsgTool::initialize() );  
  m_initialized = true;
  return StatusCode::SUCCESS;
}
 
 
  
bool IDTPM2::GenericTrackSelector::selectTrack( const TrackView& t ) const {


  
  if ( t.pt() < m_minPt )                return false;  
  if ( t.pt() > m_maxPt )                return false;

 
  if ( t.hasParent() ) {
    if ( t.parent()->ET() < m_minET     ) return false;
    if ( std::fabs(t.parent()->ET()/t.pt()) < m_minETovpt ) return false;
  }
  else { 
    if ( t.ET() < m_minET     ) return false;
    if ( std::fabs(t.ET()/t.pt()) < m_minETovpt ) return false;
  }
    
  if ( t.eta() > m_maxEta )              return false;
  if ( t.eta() < m_minEta )              return false; 
  if ( std::fabs(t.eta()) < m_minAbsEta )       return false; 
 
  if ( t.d0() < m_minD0 )                return false; 
  if ( t.d0() > m_maxD0 )                return false; 
  if ( t.z0() < m_minZ0 )                return false; 
  if ( t.z0() > m_maxZ0 )                return false; 
 
  if ( std::fabs(t.d0()) < m_minAbsD0 )         return false; 
  if ( std::fabs(t.d0()) > m_maxAbsD0 )         return false; 
  if ( std::fabs(t.z0()) < m_minAbsZ0 )         return false; 
  if ( std::fabs(t.z0()) > m_maxAbsZ0 )         return false; 
 
  if ( t.phi() < m_minPhi )              return false; 
  if ( t.phi() > m_maxPhi )              return false;  

  //  if ( std::fabs(t.phi()) < m_minAbsPhi )       return false; 
  //  if ( std::fabs(t.phi()) > m_maxAbsPhi )       return false;

  
  double qoverpt = t.qOverP()/std::sin(t.theta()); 
 
  if ( qoverpt < m_minQoPT )         return false; 
  if ( qoverpt > m_maxQoPT )         return false; 
  
  if ( std::fabs(qoverpt) < m_minAbsQoPT )  return false; 
  if ( std::fabs(qoverpt) > m_maxAbsQoPT )  return false;
  
  //}
 
  return true;
}