This class selects tracks for b-taggers. More...

#include <TrackSelector.h>

Inheritance diagram for Analysis::TrackSelector:

Collaboration diagram for Analysis::TrackSelector:

Public Member Functions
	TrackSelector (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~TrackSelector ()
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
bool	selectTrack (const Amg::Vector3D &pv, const xAOD::TrackParticle *track, double refPt=0) const
	Returns true if the argument track fulfills the selection.
bool	selectTrack (const Amg::Vector3D &pv, const xAOD::TrackParticle *track, std::bitset< 17 > &failedCuts, double refPt=0) const
	Also return individual cuts.
const std::bitset< 17 > &	currentTrackpassedCuts () const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
enum	m_Cuts { pTMin , d0Max , z0Max , sigd0Max , sigz0Max , etaMax , nHitBLayer , deadBLayer , nHitPix , nHitSct , nHitSi , nHitTrt , nHitTrtHighE , fitChi2 , fitProb , fitChi2OnNdfMax , trackingTightDef , numCuts }
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
bool	m_usePerigeeParameters {}
	Properties for V0 finding: not yet implemented.
double	m_pTMin {}
	if true use perigee parameters instead of parameters w.r.t. primary vertex
bool	m_usepTDepTrackSel {}
	min. pT: \|pT\|>pTMin
double	m_pTMinOffset {}
double	m_pTMinSlope {}
double	m_d0Max {}
double	m_z0Max {}
	max. d0: \|d0\|<d0Max
double	m_sigd0Max {}
	max. z0:
double	m_sigz0Max {}
double	m_etaMax {}
bool	m_useTrackSummaryInfo {}
	max. pseudo-rapidity
bool	m_useBLayerHitPrediction {}
	if false the following cuts are ignored
int	m_nHitBLayer {}
int	m_nHitPix {}
	at least n hits in Blayer
int	m_nHitSct {}
	at least n hits in pixels
int	m_nHitSi {}
	at least n hits in SCT
int	m_nHitTrt {}
	at least n hits in pixels+SCT
int	m_nHitTrtHighE {}
	at least n hits in TRT
bool	m_useDeadPixInfo {}
	at least n transition hits in TRT
bool	m_useDeadSctInfo {}
	if true uses dead pixel sensors from conddb (except b-layer) to compute nPix
bool	m_useTrackQualityInfo {}
	if true uses dead SCT sensors to compute nSct
double	m_fitChi2 {}
	if false the following cuts are ignored
double	m_fitProb {}
	max. fit chi2
double	m_fitChi2OnNdfMax {}
	min. fit chi2 probability
std::string	m_inputTrackCollection
std::string	m_outputTrackCollection
	location of inputTracks in StoreGate
std::string	m_sharedHitMapLocation
	location of outputTracks in StoreGate
bool	m_useAntiPileUpCuts {}
	location of shared map in StoreGate
double	m_antiPileUpSigD0Cut {}
	activate antiPU cut d0z0
double	m_antiPileUpSigZ0Cut {}
int	m_antiPileUpNHitSiCut {}
int	m_antiPileUpNHolePixCut {}
bool	m_useTrackingTightDefinition {}
int	m_ntri {}
int	m_ntrf {}
std::atomic< int >	m_ntrc [numCuts] {}
ToolHandle< Reco::ITrackToVertex >	m_trackToVertexTool {this, "trackToVertexTool", "Reco::TrackToVertex"}
std::bitset< numCuts >	m_passedCuts
bool	m_useEtaDependentCuts {}
ServiceHandle< InDet::IInDetEtaDependentCutsSvc >	m_etaDependentCutsSvc {this, "InDetEtaDependentCutsSvc", ""}
	service to get cut values depending on different variable
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

This class selects tracks for b-taggers.

Author: Laurent Vacavant (vacav.nosp@m.ant@.nosp@m.in2p3.nosp@m..fr)

Definition at line 29 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ m_Cuts

enum Analysis::TrackSelector::m_Cuts

private

Enumerator
pTMin
d0Max
z0Max
sigd0Max
sigz0Max
etaMax
nHitBLayer
deadBLayer
nHitPix
nHitSct
nHitSi
nHitTrt
nHitTrtHighE
fitChi2
fitProb
fitChi2OnNdfMax
trackingTightDef
numCuts

Definition at line 96 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

              { pTMin, d0Max, z0Max, sigd0Max, sigz0Max, etaMax,
                nHitBLayer, deadBLayer, nHitPix, nHitSct, nHitSi, nHitTrt, nHitTrtHighE,
                fitChi2, fitProb,fitChi2OnNdfMax, trackingTightDef,
                numCuts };

Constructor & Destructor Documentation

◆ TrackSelector()

Analysis::TrackSelector::TrackSelector	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Definition at line 23 of file TrackSelector.cxx.

                                                                    :
    AthAlgTool(type, name, parent),
    m_ntri(0),
    m_ntrf(0)
  {
 
    declareInterface<TrackSelector>(this);
    declareProperty("useBLayerHitPrediction", m_useBLayerHitPrediction = false);
    declareProperty("usePerigeeParameters", m_usePerigeeParameters = false);
    declareProperty("pTMin", m_pTMin = 1.*Gaudi::Units::GeV);
    declareProperty("usepTDepTrackSel", m_usepTDepTrackSel = false);
    declareProperty("pTMinOffset", m_pTMinOffset = 0.);
    declareProperty("pTMinSlope", m_pTMinSlope = 0.);
    declareProperty("d0Max", m_d0Max = 1.*Gaudi::Units::mm);
    declareProperty("z0Max", m_z0Max = 1.5*Gaudi::Units::mm);
    declareProperty("sigd0Max",m_sigd0Max = 999.*Gaudi::Units::mm);
    declareProperty("sigz0Max",m_sigz0Max = 999.*Gaudi::Units::mm);
    declareProperty("etaMax", m_etaMax = 9999.);
    declareProperty("useTrackSummaryInfo", m_useTrackSummaryInfo = true);
    declareProperty("nHitBLayer", m_nHitBLayer = 1);
    declareProperty("nHitPix", m_nHitPix = 1);
    declareProperty("nHitSct", m_nHitSct = 0);
    declareProperty("nHitSi", m_nHitSi = 7);
    declareProperty("nHitTrt", m_nHitTrt = 0);
    declareProperty("nHitTrtHighE", m_nHitTrtHighE = 0);
    declareProperty("useDeadPixInfo", m_useDeadPixInfo = true);
    declareProperty("useDeadSctInfo", m_useDeadSctInfo = true);
    declareProperty("useTrackQualityInfo", m_useTrackQualityInfo = true);
    declareProperty("fitChi2", m_fitChi2 = 99999.);
    declareProperty("fitProb", m_fitProb = -1.);
    declareProperty("fitChi2OnNdfMax",m_fitChi2OnNdfMax=999.);
    declareProperty("inputTrackCollection", m_inputTrackCollection);
    declareProperty("outputTrackCollection", m_outputTrackCollection);
    declareProperty("useAntiPileUpCuts", m_useAntiPileUpCuts = false);
    declareProperty("antiPileUpSigD0Cut", m_antiPileUpSigD0Cut = 3.0);
    declareProperty("antiPileUpSigZ0Cut", m_antiPileUpSigZ0Cut = 3.8);
    declareProperty("antiPileUpNHitSiCut", m_antiPileUpNHitSiCut = 9);
    declareProperty("antiPileUpNHolePixCut", m_antiPileUpNHolePixCut = 9);
    declareProperty("useTrackingTightDefinition", m_useTrackingTightDefinition = false); 
  }

◆ ~TrackSelector()

Analysis::TrackSelector::~TrackSelector ( )

virtual

Definition at line 65 of file TrackSelector.cxx.

65 {

66 }

Member Function Documentation

◆ currentTrackpassedCuts()

const std::bitset< 17 > & Analysis::TrackSelector::currentTrackpassedCuts ( ) const

inline

Definition at line 54 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

54{ return m_passedCuts; }

Analysis::TrackSelector::m_passedCuts

std::bitset< numCuts > m_passedCuts

Definition PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h:107

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	hndl,
		const SG::VarHandleKeyType &	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95{ return m_detStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85{ return m_evtStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ finalize()

StatusCode Analysis::TrackSelector::finalize ( )

overridevirtual

Definition at line 161 of file TrackSelector.cxx.

                                     {
    ATH_MSG_VERBOSE("#BTAG#  tracks selected: In= " << m_ntri);
    for(int i=0;i<numCuts;i++) ATH_MSG_VERBOSE("#BTAG#  cut" << i << "= " << m_ntrc[i]);
    ATH_MSG_VERBOSE("#BTAG#  Out= " << m_ntrf);
    return StatusCode::SUCCESS;
  }

◆ initialize()

StatusCode Analysis::TrackSelector::initialize ( )

overridevirtual

retrieving TrackToVertex:

dump cuts:

Definition at line 72 of file TrackSelector.cxx.

                                       {
    for(int i=0;i<numCuts;i++) m_ntrc[i]=0;
    if ( m_trackToVertexTool.retrieve().isFailure() ) {
      ATH_MSG_FATAL("#BTAG# Failed to retrieve tool " << m_trackToVertexTool);
      return StatusCode::FAILURE;
    } else {
      ATH_MSG_DEBUG("#BTAG# Retrieved tool " << m_trackToVertexTool);
    }
 
    m_useEtaDependentCuts = !(m_etaDependentCutsSvc.name().empty());
    if (m_useEtaDependentCuts){
       ATH_CHECK(m_etaDependentCutsSvc.retrieve());
       ATH_MSG_INFO("Using InDetEtaDependentCutsSvc. Track selections from config not used");
     }
     else{
       ATH_MSG_INFO("Using track selections from config");
     }
 
    if (m_useTrackingTightDefinition){
      m_d0Max = 2*Gaudi::Units::mm;
      m_z0Max = 3*Gaudi::Units::mm;
      m_pTMin = 0.4*Gaudi::Units::GeV;
    }
 
    if (msgLvl(MSG::DEBUG)) {
      msg(MSG::DEBUG) << "#BTAG# TrackSelector " << name() << " cuts: " << endmsg;
      if(!m_useEtaDependentCuts){
    msg(MSG::DEBUG) << "#BTAG#     - pT >= " << m_pTMin << endmsg;
    msg(MSG::DEBUG) << "#BTAG#     - |eta| <= " << m_etaMax << endmsg;
    msg(MSG::DEBUG) << "#BTAG#     - |d0| <= " << m_d0Max << endmsg;
    msg(MSG::DEBUG) << "#BTAG#     - |z0| <= " << m_z0Max << endmsg;
      }
 
      msg(MSG::DEBUG) << "#BTAG#     - |sigd0| <= " << m_sigd0Max << endmsg;
      msg(MSG::DEBUG) << "#BTAG#     - |sigz0| <= " << m_sigz0Max << endmsg;
      if(m_useAntiPileUpCuts) {
        msg(MSG::DEBUG) << "#BTAG#     - antiPUcut: reject tracks with |sigz0| > " << m_antiPileUpSigZ0Cut 
                        << " when |sigd0| < " << m_antiPileUpSigD0Cut << endmsg;
      }
 
      if(m_useTrackSummaryInfo) {
 
    if(!m_useEtaDependentCuts){
      msg(MSG::DEBUG) << "#BTAG#     - nbHitsBLayer >= " << m_nHitBLayer << endmsg;
      if(m_useDeadPixInfo) {
        msg(MSG::DEBUG) << "#BTAG#     - nbHitsPix+nbDeadPix >= " << m_nHitPix << endmsg;
      } else {
        msg(MSG::DEBUG) << "#BTAG#     - nbHitsPix >= " << m_nHitPix << endmsg;
      }
      if(m_useBLayerHitPrediction)
        msg(MSG::DEBUG) << "#BTAG#     using conddb for b-layer hit requirements " << endmsg;
      
      if(m_useDeadSctInfo) {
        msg(MSG::DEBUG) << "#BTAG#     - nbHitsSct+nbDeadSct >= " << m_nHitSct << endmsg;
      } else {
        msg(MSG::DEBUG) << "#BTAG#     - nbHitsSct >= " << m_nHitSct << endmsg;
      }
      int nhsi = m_nHitSi;
      if(m_useAntiPileUpCuts) nhsi = m_antiPileUpNHitSiCut;
      if(m_useDeadPixInfo) {
        if(m_useDeadSctInfo) {
          msg(MSG::DEBUG) << "#BTAG#     - nbHitsSi+nbDeadPix+nbDeadSct >= " << nhsi << endmsg;
        } else {
          msg(MSG::DEBUG) << "#BTAG#     - nbHitsSi+nbDeadPix >= " << nhsi << endmsg;
        }
      } else {
        if(m_useDeadSctInfo) {
          msg(MSG::DEBUG) << "#BTAG#     - nbHitsSi+nbDeadSct >= " << nhsi << endmsg;
        } else {
          msg(MSG::DEBUG) << "#BTAG#     - nbHitsSi >= " << nhsi << endmsg;
        }
      }
    }
 
    msg(MSG::DEBUG) << "#BTAG#     - nbHitsTrt >= " << m_nHitTrt << endmsg;
    msg(MSG::DEBUG) << "#BTAG#     - nbHitsTrtHighE >= " << m_nHitTrtHighE << endmsg;
      }
 
      msg(MSG::DEBUG) << "#BTAG#     - fit chi2 <= " << m_fitChi2 << endmsg;
      msg(MSG::DEBUG) << "#BTAG#     - fit proba >= " << m_fitProb << endmsg;
      msg(MSG::DEBUG) << "#BTAG#     - fit chi2 / ndf <= " << m_fitChi2OnNdfMax << endmsg;
    }
 
    return StatusCode::SUCCESS;
  }

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

const InterfaceID & Analysis::TrackSelector::interfaceID ( )

static

Definition at line 68 of file TrackSelector.cxx.

                                                {
    return IID_ITrackSelector;
  }

◆ msg()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

◆ selectTrack() [1/2]

bool Analysis::TrackSelector::selectTrack	(	const Amg::Vector3D &	pv,
		const xAOD::TrackParticle *	track,
		double	refPt = 0 ) const

Returns true if the argument track fulfills the selection.

Definition at line 169 of file TrackSelector.cxx.

{
  std::bitset<17> failedCuts;
  return selectTrack (pv, track, failedCuts, refPt);
}

◆ selectTrack() [2/2]

bool Analysis::TrackSelector::selectTrack	(	const Amg::Vector3D &	pv,
		const xAOD::TrackParticle *	track,
		std::bitset< 17 > &	failedCuts,
		double	refPt = 0 ) const

Also return individual cuts.

apply cuts:

std::string showPassedCuts = ~failedCuts.to_string(); // not available !

Definition at line 177 of file TrackSelector.cxx.

{
    double trackD0;
    double trackZ0;
    double tracksigD0;
    double tracksigZ0;
    if(m_usePerigeeParameters) {
      trackD0 = track->d0();
      trackZ0 = track->z0() - pv.z();
      tracksigD0 = std::sqrt(track->definingParametersCovMatrix()(0,0));
      tracksigZ0 = std::sqrt(track->definingParametersCovMatrix()(1,1));
    } else {
      // extrapolate with the TrackToVertex tool:
      auto ctx = Gaudi::Hive::currentContext();
      auto perigee = m_trackToVertexTool->perigeeAtVertex(ctx, *track, pv);
      if (perigee==nullptr) {
        ATH_MSG_WARNING("#BTAG#  Extrapolation failed. Rejecting track... ");
        return false;
      }
      trackD0 = perigee->parameters()[Trk::d0];
      trackZ0 = perigee->parameters()[Trk::z0];
      tracksigD0 = std::sqrt((*perigee->covariance())(Trk::d0,Trk::d0));
      tracksigZ0 = std::sqrt((*perigee->covariance())(Trk::z0,Trk::z0));
    }
 
    ATH_MSG_VERBOSE( "#BTAG#  Track "
             << " Eta= " << track->eta() << " Phi= " << track->phi() << " pT= " <<track->pt()
             << " d0= " << trackD0
             << " z0= " << trackZ0 << " sigd0= " << tracksigD0 << " sigz0: " << tracksigZ0 );
    double eta = track->eta();
    bool pass = true;
    double pTMin_cut = m_useEtaDependentCuts ? m_etaDependentCutsSvc->getMinPtAtEta(eta) : m_pTMin;
    if(track->pt() < pTMin_cut) {
      pass = false;
      failedCuts.set(pTMin);
    } else if (refPt > 0 && m_usepTDepTrackSel && track->pt() < m_pTMinOffset + m_pTMinSlope*refPt) {
      pass = false;
      failedCuts.set(pTMin);
    }
 
    double d0_cut = m_useEtaDependentCuts ? m_etaDependentCutsSvc->getMaxPrimaryImpactAtEta(eta) : m_d0Max;
    if(std::abs(trackD0)>d0_cut) {
      pass = false;
      failedCuts.set(d0Max);
    }
 
    double z0_cut = m_useEtaDependentCuts ? m_etaDependentCutsSvc->getMaxZImpactAtEta(eta) : m_z0Max;
    if(std::abs(trackZ0*sin(track->theta()))>z0_cut) {
      pass = false;
      failedCuts.set(z0Max);
    }
    if (tracksigD0>m_sigd0Max) {
      pass = false;
      failedCuts.set(sigd0Max);
    }
    if (tracksigZ0>m_sigz0Max) {
      pass = false;
      failedCuts.set(sigz0Max);
    }
    if(m_useAntiPileUpCuts) {
      if(fabs(trackZ0/tracksigZ0)>m_antiPileUpSigZ0Cut && fabs(trackD0/tracksigD0)<m_antiPileUpSigD0Cut) {
        pass = false;
        failedCuts.set(sigz0Max);
        failedCuts.set(sigd0Max);
      }
    }
 
    double eta_cut = m_useEtaDependentCuts ? m_etaDependentCutsSvc->getMaxEta() : m_etaMax;
    if(std::abs(track->eta())>eta_cut) {
      pass = false;
      failedCuts.set(etaMax);
    }
    if(m_useTrackSummaryInfo) {
      uint8_t nb=0;
      track->summaryValue(nb, xAOD::numberOfInnermostPixelLayerHits);
      int InPix_cut = m_useEtaDependentCuts ? m_etaDependentCutsSvc->getMinInnermostPixelHitsAtEta(eta) : m_nHitBLayer;
      if(nb < InPix_cut) {
    failedCuts.set(nHitBLayer);
    if(!m_useBLayerHitPrediction) { 
      pass = false;
      failedCuts.set(deadBLayer);
    } else {
      uint8_t ehib=1;
      if (!track->summaryValue(ehib,xAOD::expectInnermostPixelLayerHit)) {
        ATH_MSG_WARNING("#BTAG# expectInnermostPixelLayerHit not computed in  TrackSummary: assuming true");
        ehib=1;
      }
      if(ehib) {  // check if module was alive
        pass = false;
        failedCuts.set(deadBLayer);
      }
    }
      }
      uint8_t nhp=0;
      track->summaryValue(nhp, xAOD::numberOfPixelHoles);
      if(m_useAntiPileUpCuts) {
    int PixHole_cut = m_useEtaDependentCuts ? m_etaDependentCutsSvc->getMaxPixelHolesAtEta(eta) : m_antiPileUpNHolePixCut;
    if(nhp>=PixHole_cut) {
      pass = false;
    }
      }
      uint8_t np=0;
      track->summaryValue(np, xAOD::numberOfPixelHits);
      if(m_useDeadPixInfo) 
      {
    uint8_t ndead;
    track->summaryValue(ndead, xAOD::numberOfPixelDeadSensors);
    np += std::max((int)ndead, 0);
      }
      int Pix_cut = m_useEtaDependentCuts ? m_etaDependentCutsSvc->getMinPixelHitsAtEta(eta) : m_nHitPix;
      if(np < Pix_cut) {
    pass = false;
    failedCuts.set(nHitPix);
      }
 
      uint8_t ns=0;
      track->summaryValue(ns, xAOD::numberOfSCTHits);
      if(m_useDeadSctInfo)
      {
    uint8_t ndead;
        track->summaryValue(ndead, xAOD::numberOfSCTDeadSensors);
    ns += std::max((int)ndead,0);
      }
      int SCTStrip_cut = m_useEtaDependentCuts ? m_etaDependentCutsSvc->getMinStripHitsAtEta(eta) : m_nHitSct;
      if(ns < SCTStrip_cut) {
    pass = false;
    failedCuts.set(nHitSct);
      }
 
      int Si_cut = m_useEtaDependentCuts ? m_etaDependentCutsSvc->getMinSiHitsAtEta(eta) : m_nHitSi;
      if((np+ns) < Si_cut) {
    pass = false;
    failedCuts.set(nHitSi);
      }
      if(m_useAntiPileUpCuts) {
    if((np+ns) < m_antiPileUpNHitSiCut) {
      pass = false;
      failedCuts.set(nHitSi);
    }
      }
      uint8_t nh=0;
      track->summaryValue(nh, xAOD::numberOfTRTHits);//ms
      if(nh < m_nHitTrt) {
    pass = false;
    failedCuts.set(nHitTrt);
      }
      uint8_t nhe=0;
      track->summaryValue(nhe, xAOD::numberOfTRTHighThresholdHits);//ms
      if(nhe < m_nHitTrtHighE) {
    pass = false;
    failedCuts.set(nHitTrtHighE);
      }
      // Tracking CP group Tight definition
      if(m_useTrackingTightDefinition) {
    uint8_t nibl, nnibl;
    track->summaryValue(nibl , xAOD::numberOfPixelHits);
    track->summaryValue(nnibl, xAOD::numberOfPixelHits);
    bool innerHitsCrit = ((nibl+nnibl)>0);
    bool lowetaCrit  = fabs(track->eta())> 1.65 && (ns+np)>=11;
    bool highetaCrit = fabs(track->eta())<=1.65 && (ns+np)>=9 ;
    bool pixholeCrit = (nhp==0) ;
    bool isTight = innerHitsCrit && pixholeCrit && (lowetaCrit || highetaCrit);
    if (!isTight){
      pass=false;
      failedCuts.set(trackingTightDef);
    }
      
      }
    }
    
    // Now the fit quality
    double chi2 = track->chiSquared();
    int ndf = track->numberDoF();
    double proba = 1.;
    if(ndf>0 && chi2>=0.) {
      Genfun::CumulativeChiSquare myCumulativeChiSquare(ndf);
      proba = 1.-myCumulativeChiSquare(chi2);
    }
    if(chi2>m_fitChi2) {
      pass = false;
      failedCuts.set(fitChi2);
    }
    if(proba<m_fitProb) {
      pass = false;
      failedCuts.set(fitProb);
    }
    if(ndf>0) {
      if(chi2/double(ndf)>m_fitChi2OnNdfMax) {
        pass = false;
        failedCuts.set(fitChi2OnNdfMax);
      }
    } else {
      pass = false;
      failedCuts.set(fitChi2OnNdfMax);
    }
    if( msgLvl(MSG::VERBOSE) ){
      ATH_MSG_VERBOSE("#BTAG#  passedCuts for track ");
      for(int i=0;i<numCuts;i++) {
    int passl = ~failedCuts[(m_Cuts)i];
    if(passl) m_ntrc[i]++;
    msg(MSG::VERBOSE) << passl;
      } 
      msg(MSG::VERBOSE) << endmsg;
    }
 
    return pass;
  }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

◆ m_antiPileUpNHitSiCut

int Analysis::TrackSelector::m_antiPileUpNHitSiCut {}

private

Definition at line 92 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

92{};

◆ m_antiPileUpNHolePixCut

int Analysis::TrackSelector::m_antiPileUpNHolePixCut {}

private

Definition at line 93 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

93{};

◆ m_antiPileUpSigD0Cut

double Analysis::TrackSelector::m_antiPileUpSigD0Cut {}

private

activate antiPU cut d0z0

Definition at line 90 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

90{};

◆ m_antiPileUpSigZ0Cut

double Analysis::TrackSelector::m_antiPileUpSigZ0Cut {}

private

Definition at line 91 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

91{};

◆ m_d0Max

double Analysis::TrackSelector::m_d0Max {}

private

Definition at line 67 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

67{};

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_etaDependentCutsSvc

ServiceHandle<InDet::IInDetEtaDependentCutsSvc> Analysis::TrackSelector::m_etaDependentCutsSvc {this, "InDetEtaDependentCutsSvc", ""}

private

service to get cut values depending on different variable

Definition at line 111 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

111{this, "InDetEtaDependentCutsSvc", ""};

◆ m_etaMax

double Analysis::TrackSelector::m_etaMax {}

private

Definition at line 71 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

71{};

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_fitChi2

double Analysis::TrackSelector::m_fitChi2 {}

private

if false the following cuts are ignored

Definition at line 83 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

83{};

◆ m_fitChi2OnNdfMax

double Analysis::TrackSelector::m_fitChi2OnNdfMax {}

private

min. fit chi2 probability

Definition at line 85 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

85{}; //max. fitchi2/ndf

◆ m_fitProb

double Analysis::TrackSelector::m_fitProb {}

private

max. fit chi2

Definition at line 84 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

84{};

◆ m_inputTrackCollection

std::string Analysis::TrackSelector::m_inputTrackCollection

private

Definition at line 86 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

◆ m_nHitBLayer

int Analysis::TrackSelector::m_nHitBLayer {}

private

Definition at line 74 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

74{};

◆ m_nHitPix

int Analysis::TrackSelector::m_nHitPix {}

private

at least n hits in Blayer

Definition at line 75 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

75{};

◆ m_nHitSct

int Analysis::TrackSelector::m_nHitSct {}

private

at least n hits in pixels

Definition at line 76 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

76{};

◆ m_nHitSi

int Analysis::TrackSelector::m_nHitSi {}

private

at least n hits in SCT

Definition at line 77 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

77{};

◆ m_nHitTrt

int Analysis::TrackSelector::m_nHitTrt {}

private

at least n hits in pixels+SCT

Definition at line 78 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

78{};

◆ m_nHitTrtHighE

int Analysis::TrackSelector::m_nHitTrtHighE {}

private

at least n hits in TRT

Definition at line 79 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

79{};

◆ m_ntrc

std::atomic<int> Analysis::TrackSelector::m_ntrc[numCuts] {}

mutableprivate

Definition at line 103 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

103{};

◆ m_ntrf

int Analysis::TrackSelector::m_ntrf {}

private

Definition at line 102 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

102{};

◆ m_ntri

int Analysis::TrackSelector::m_ntri {}

private

Definition at line 101 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

101{};

◆ m_outputTrackCollection

std::string Analysis::TrackSelector::m_outputTrackCollection

private

location of inputTracks in StoreGate

Definition at line 87 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

◆ m_passedCuts

std::bitset<numCuts> Analysis::TrackSelector::m_passedCuts

private

Definition at line 107 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

◆ m_pTMin

double Analysis::TrackSelector::m_pTMin {}

private

if true use perigee parameters instead of parameters w.r.t. primary vertex

Definition at line 63 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

63{};

◆ m_pTMinOffset

double Analysis::TrackSelector::m_pTMinOffset {}

private

Definition at line 65 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

65{};

◆ m_pTMinSlope

double Analysis::TrackSelector::m_pTMinSlope {}

private

Definition at line 66 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

66{};

◆ m_sharedHitMapLocation

std::string Analysis::TrackSelector::m_sharedHitMapLocation

private

location of outputTracks in StoreGate

Definition at line 88 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

◆ m_sigd0Max

double Analysis::TrackSelector::m_sigd0Max {}

private

max. z0:

Definition at line 69 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

69{}; // max sig d0:

◆ m_sigz0Max

double Analysis::TrackSelector::m_sigz0Max {}

private

Definition at line 70 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

70{}; // max sig z0:

◆ m_trackToVertexTool

ToolHandle< Reco::ITrackToVertex > Analysis::TrackSelector::m_trackToVertexTool {this, "trackToVertexTool", "Reco::TrackToVertex"}

private

Definition at line 105 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

105{this, "trackToVertexTool", "Reco::TrackToVertex"};

◆ m_useAntiPileUpCuts

bool Analysis::TrackSelector::m_useAntiPileUpCuts {}

private

location of shared map in StoreGate

Definition at line 89 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

89{};

◆ m_useBLayerHitPrediction

bool Analysis::TrackSelector::m_useBLayerHitPrediction {}

private

if false the following cuts are ignored

Definition at line 73 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

73{}; // use b-layer tool tied to conddb

◆ m_useDeadPixInfo

bool Analysis::TrackSelector::m_useDeadPixInfo {}

private

at least n transition hits in TRT

Definition at line 80 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

80{};

◆ m_useDeadSctInfo

bool Analysis::TrackSelector::m_useDeadSctInfo {}

private

if true uses dead pixel sensors from conddb (except b-layer) to compute nPix

Definition at line 81 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

81{};

◆ m_useEtaDependentCuts

bool Analysis::TrackSelector::m_useEtaDependentCuts {}

private

Definition at line 109 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

109{};

◆ m_usePerigeeParameters

bool Analysis::TrackSelector::m_usePerigeeParameters {}

private

Properties for V0 finding: not yet implemented.

Properties for track selection: all cuts are ANDed

Definition at line 62 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

62{};

◆ m_usepTDepTrackSel

bool Analysis::TrackSelector::m_usepTDepTrackSel {}

private

min. pT: |pT|>pTMin

Definition at line 64 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

64{};

◆ m_useTrackingTightDefinition

bool Analysis::TrackSelector::m_useTrackingTightDefinition {}

private

Definition at line 94 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

94{};

◆ m_useTrackQualityInfo

bool Analysis::TrackSelector::m_useTrackQualityInfo {}

private

if true uses dead SCT sensors to compute nSct

Definition at line 82 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

82{};

◆ m_useTrackSummaryInfo

bool Analysis::TrackSelector::m_useTrackSummaryInfo {}

private

max. pseudo-rapidity

Definition at line 72 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

72{};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_z0Max

double Analysis::TrackSelector::m_z0Max {}

private

max. d0: |d0|<d0Max

Definition at line 68 of file PhysicsAnalysis/JetTagging/JetTagTools/JetTagTools/TrackSelector.h.

68{};

The documentation for this class was generated from the following files:

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Member Enumeration Documentation

◆ m_Cuts

Constructor & Destructor Documentation

◆ TrackSelector()

◆ ~TrackSelector()

Member Function Documentation

◆ currentTrackpassedCuts()

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ evtStore()

◆ extraDeps_update_handler()

◆ finalize()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ selectTrack() [1/2]

◆ selectTrack() [2/2]

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_antiPileUpNHitSiCut

◆ m_antiPileUpNHolePixCut

◆ m_antiPileUpSigD0Cut

◆ m_antiPileUpSigZ0Cut

◆ m_d0Max

◆ m_detStore

◆ m_etaDependentCutsSvc

◆ m_etaMax

◆ m_evtStore

◆ m_fitChi2

◆ m_fitChi2OnNdfMax

◆ m_fitProb

◆ m_inputTrackCollection

◆ m_nHitBLayer

◆ m_nHitPix

◆ m_nHitSct

◆ m_nHitSi

◆ m_nHitTrt

◆ m_nHitTrtHighE

◆ m_ntrc

◆ m_ntrf

◆ m_ntri

◆ m_outputTrackCollection

◆ m_passedCuts

◆ m_pTMin

◆ m_pTMinOffset

◆ m_pTMinSlope

◆ m_sharedHitMapLocation

◆ m_sigd0Max

◆ m_sigz0Max

◆ m_trackToVertexTool

◆ m_useAntiPileUpCuts

◆ m_useBLayerHitPrediction

◆ m_useDeadPixInfo

◆ m_useDeadSctInfo

◆ m_useEtaDependentCuts

◆ m_usePerigeeParameters

◆ m_usepTDepTrackSel

◆ m_useTrackingTightDefinition

◆ m_useTrackQualityInfo

◆ m_useTrackSummaryInfo

◆ m_varHandleArraysDeclared

◆ m_vhka

◆ m_z0Max