Associate tracks to the tau candidate. More...

#include <TauTrackFinder.h>

Inheritance diagram for TauTrackFinder:

Collaboration diagram for TauTrackFinder:

Public Types
enum	TauTrackType { TauTrackCore = 0 , TauTrackWide = 1 , TauTrackOther = 2 , NotTauTrack = 3 }
	Enumerator defining type of tau track. More...

Public Member Functions
	TauTrackFinder (const std::string &name)
	Constructor and Destructor.
	ASG_TOOL_CLASS2 (TauTrackFinder, TauRecToolBase, ITauToolBase)
	~TauTrackFinder ()
virtual StatusCode	initialize () override
	Algorithm functions.
virtual StatusCode	executeTrackFinder (xAOD::TauJet &pTau, xAOD::TauTrackContainer &tauTrackCon) const override
virtual StatusCode	eventInitialize () override
	Event initializer - called at the beginning of each event.
virtual StatusCode	execute (xAOD::TauJet &pTau) const override
	Execute - called for each tau candidate.
virtual StatusCode	executePi0CreateROI (xAOD::TauJet &pTau, CaloConstCellContainer &caloCellContainer, boost::dynamic_bitset<> &map) const override
virtual StatusCode	executeVertexFinder (xAOD::TauJet &pTau, const xAOD::VertexContainer *vertexContainer=nullptr) const override
virtual StatusCode	executeTrackClassifier (xAOD::TauJet &pTau, xAOD::TauTrackContainer &tauTrackContainer) const override
virtual StatusCode	executeShotFinder (xAOD::TauJet &pTau, xAOD::CaloClusterContainer &shotClusterContainer, xAOD::PFOContainer &PFOContainer) const override
virtual StatusCode	executePi0ClusterCreator (xAOD::TauJet &pTau, xAOD::PFOContainer &neutralPFOContainer, xAOD::PFOContainer &hadronicPFOContainer, const xAOD::CaloClusterContainer &pCaloClusterContainer) const override
virtual StatusCode	executeVertexVariables (xAOD::TauJet &pTau, xAOD::VertexContainer &vertexContainer) const override
virtual StatusCode	executePi0ClusterScaler (xAOD::TauJet &pTau, xAOD::PFOContainer &neutralPFOContainer, xAOD::PFOContainer &chargedPFOContainer) const override
virtual StatusCode	executePi0nPFO (xAOD::TauJet &pTau, xAOD::PFOContainer &neutralPFOContainer) const override
virtual StatusCode	executePanTau (xAOD::TauJet &pTau, xAOD::ParticleContainer &particleContainer, xAOD::PFOContainer &neutralPFOContainer) const override
virtual StatusCode	eventFinalize () override
	Event finalizer - called at the end of each event.
virtual StatusCode	finalize () override
	Finalizer.
std::string	find_file (const std::string &fname) const
virtual void	print () const
	Print the state of the tool.
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
Additional helper functions, not directly mimicking Athena
template<class T>
const T *	getProperty (const std::string &name) const
	Get one of the tool's properties.
const std::string &	msg_level_name () const __attribute__((deprecated))
	A deprecated function for getting the message level's name.
const std::string &	getName (const void *ptr) const
	Get the name of an object that is / should be in the event store.
SG::sgkey_t	getKey (const void *ptr) const
	Get the (hashed) key of an object that is in the event store.

Protected Member Functions
bool	inTrigger () const
bool	inAOD () const
bool	inEleRM () const
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.

Protected Attributes
Gaudi::Property< bool >	m_in_trigger {this, "inTrigger", false, "Indicate if the tool is running on trigger"}
Gaudi::Property< bool >	m_in_AOD {this, "inAOD", false, "Indicate if the tool is running on AOD"}
Gaudi::Property< bool >	m_in_EleRM {this, "inEleRM", false, "Indicate if the tool is running on EleRM routine"}
Gaudi::Property< std::string >	m_tauRecToolsTag {this, "calibFolder", "tauRecTools/R22_preprod", "CVMFS path to the tau calibration folder"}

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	extrapolateToCaloSurface (const EventContext &ctx, xAOD::TauJet &pTau, xAOD::TauTrackContainer &tauTrackCon) const
	Extrapolate track eta and phi to the calorimeter middle surface.
TauTrackType	tauTrackType (const xAOD::TauJet &tauJet, const xAOD::TrackParticle &trackParticle, const xAOD::Vertex *primaryVertex) const
void	getTauTracksFromPV (const xAOD::TauJet &tauJet, const std::vector< const xAOD::TrackParticle * > &vecTrackParticles, const xAOD::Vertex primaryVertex, const bool &useGhostTracks, const xAOD::JetContainer jetContainer, std::vector< const xAOD::TrackParticle * > &tauTracks, std::vector< const xAOD::TrackParticle * > &wideTracks, std::vector< const xAOD::TrackParticle * > &otherTracks) const
void	removeOffsideTracksWrtLeadTrk (const EventContext &ctx, std::vector< const xAOD::TrackParticle * > &tauTracks, std::vector< const xAOD::TrackParticle * > &wideTracks, std::vector< const xAOD::TrackParticle * > &otherTracks, const xAOD::Vertex *tauOrigin, double maxDeltaZ0) const
float	getZ0 (const EventContext &ctx, const xAOD::TrackParticle track, const xAOD::Vertex vertex) const
	Some internally used functions.
bool	isLargeD0Track (const xAOD::TrackParticle *track) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ToolHandle< Trk::IParticleCaloExtensionTool >	m_caloExtensionTool {this, "ParticleCaloExtensionTool", "Trk::ParticleCaloExtensionTool/ParticleCaloExtensionTool", "Tool for the extrapolation of charged tracks"}
	tools
ToolHandle< Trk::ITrackSelectorTool >	m_trackSelectorTool_tau {this, "TrackSelectorToolTau", "", "Tool for track selection"}
ToolHandle< Reco::ITrackToVertex >	m_trackToVertexTool {this, "TrackToVertexTool", "Reco::TrackToVertex"}
ToolHandle< Trk::ITrackToVertexIPEstimator >	m_trackToVertexIPEstimator {this, "TrackToVertexIPEstimator", ""}
Gaudi::Property< double >	m_maxJetDr_tau {this, "MaxJetDrTau", 0.2}
Gaudi::Property< double >	m_maxJetDr_wide {this, "MaxJetDrWide", 0.4}
Gaudi::Property< bool >	m_applyZ0cut {this, "removeTracksOutsideZ0wrtLeadTrk", false}
Gaudi::Property< float >	m_z0maxDelta {this, "maxDeltaZ0wrtLeadTrk", 1000.}
Gaudi::Property< bool >	m_removeDuplicateCoreTracks {this, "removeDuplicateCoreTracks", true}
Gaudi::Property< bool >	m_bypassExtrapolator {this, "BypassExtrapolator", false}
Gaudi::Property< bool >	m_useGhostTracks {this, "useGhostTracks", false}
Gaudi::Property< double >	m_ghostTrackDR {this, "ghostTrackDR", 0.25}
SG::ReadHandleKey< xAOD::TrackParticleContainer >	m_trackPartInputContainer {this,"Key_trackPartInputContainer", "InDetTrackParticles", "input track particle container key"}
SG::ReadHandleKey< xAOD::TrackParticleContainer >	m_largeD0TracksInputContainer {this,"Key_LargeD0TrackInputContainer", "", "input LRT particle container key"}
SG::ReadHandleKey< xAOD::JetContainer >	m_jetContainer {this,"Key_jetContainer", "", "Name of the seed jet container, when using ghost matching"}
SG::ReadHandleKey< CaloExtensionCollection >	m_ParticleCacheKey {this,"tauParticleCache", "ParticleCaloExtension", "Name of the particle measurement extrapolation cache for TauTrackFinder"}
SG::ReadCondHandleKey< InDet::BeamSpotData >	m_beamSpotKey { this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot" }
std::set< CaloSampling::CaloSample >	m_EMSamplings
std::set< CaloSampling::CaloSample >	m_HadSamplings
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

Associate tracks to the tau candidate.

The tracks have to pass dedicated quality criteria and have to match to a primary vertex consistent with the tau origin.

Author: KG Tan Kong..nosp@m.Guan.nosp@m..Tan@.nosp@m.cern.nosp@m..ch; Felix Friedrich

Definition at line 39 of file TauTrackFinder.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

◆ TauTrackType

enum TauTrackFinder::TauTrackType

Enumerator defining type of tau track.

Enumerator
TauTrackCore
TauTrackWide
TauTrackOther
NotTauTrack

Definition at line 51 of file TauTrackFinder.h.

    {
        TauTrackCore  = 0,
        TauTrackWide  = 1,
        TauTrackOther = 2,
        NotTauTrack   = 3
    };

Constructor & Destructor Documentation

◆ TauTrackFinder()

TauTrackFinder::TauTrackFinder ( const std::string & name )

Constructor and Destructor.

Definition at line 13 of file TauTrackFinder.cxx.

                                                    :
  TauRecToolBase(name),
  m_EMSamplings {CaloSampling::EME1, CaloSampling::EMB1},
  m_HadSamplings {CaloSampling::TileBar1, CaloSampling::HEC1, CaloSampling::TileExt1}
{}

◆ ~TauTrackFinder()

TauTrackFinder::~TauTrackFinder ( )

Definition at line 19 of file TauTrackFinder.cxx.

19{}

Member Function Documentation

◆ ASG_TOOL_CLASS2()

TauTrackFinder::ASG_TOOL_CLASS2	(	TauTrackFinder	,
		TauRecToolBase	,
		ITauToolBase	)

◆ 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

◆ eventFinalize()

StatusCode TauRecToolBase::eventFinalize ( )

overridevirtualinherited

Event finalizer - called at the end of each event.

Implements ITauToolBase.

Definition at line 105 of file TauRecToolBase.cxx.

                                         {
  return StatusCode::SUCCESS;
}

◆ eventInitialize()

StatusCode TauRecToolBase::eventInitialize ( )

overridevirtualinherited

Event initializer - called at the beginning of each event.

Implements ITauToolBase.

Definition at line 36 of file TauRecToolBase.cxx.

                                          {
  return StatusCode::SUCCESS;
}

◆ 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

◆ execute()

StatusCode TauRecToolBase::execute ( xAOD::TauJet & pTau ) const

overridevirtualinherited

Execute - called for each tau candidate.

Implements ITauToolBase.

Reimplemented in JetSeedBuilder, MvaTESEvaluator, MvaTESVariableDecorator, TauAODMuonRemovalTool, TauAODSelector, TauAxisSetter, TauCalibrateLC, TauCellVariables, TauClusterFinder, TauCombinedTES, TauCommonCalcVars, TauDecayModeNNClassifier, TauElectronVetoVariables, TauEleOverlapChecker, TauGNNEvaluator, TauIDVarCalculator, TauJetRNNEvaluator, TauSubstructureVariables, TauVertexedClusterDecorator, and TauWPDecorator.

Definition at line 41 of file TauRecToolBase.cxx.

                                                    {
  ATH_MSG_ERROR("function not implemented");
  return StatusCode::FAILURE;
}

◆ executePanTau()

StatusCode TauRecToolBase::executePanTau	(	xAOD::TauJet &	pTau,
		xAOD::ParticleContainer &	particleContainer,
		xAOD::PFOContainer &	neutralPFOContainer ) const

overridevirtualinherited

Implements ITauToolBase.

Reimplemented in PanTau::PanTauProcessor.

Definition at line 100 of file TauRecToolBase.cxx.

                                                                                                           {
  ATH_MSG_ERROR("function not implemented");
  return StatusCode::FAILURE;
}

◆ executePi0ClusterCreator()

StatusCode TauRecToolBase::executePi0ClusterCreator	(	xAOD::TauJet &	pTau,
		xAOD::PFOContainer &	neutralPFOContainer,
		xAOD::PFOContainer &	hadronicPFOContainer,
		const xAOD::CaloClusterContainer &	pCaloClusterContainer ) const

overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauPi0ClusterCreator.

Definition at line 78 of file TauRecToolBase.cxx.

                                                                    {
  ATH_MSG_ERROR("function not implemented");
  return StatusCode::FAILURE;
}

◆ executePi0ClusterScaler()

StatusCode TauRecToolBase::executePi0ClusterScaler	(	xAOD::TauJet &	pTau,
		xAOD::PFOContainer &	neutralPFOContainer,
		xAOD::PFOContainer &	chargedPFOContainer ) const

overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauPi0ClusterScaler.

Definition at line 90 of file TauRecToolBase.cxx.

                                                                                                                 {
  ATH_MSG_ERROR("function not implemented");
  return StatusCode::FAILURE;
} 

◆ executePi0CreateROI()

StatusCode TauRecToolBase::executePi0CreateROI	(	xAOD::TauJet &	pTau,
		CaloConstCellContainer &	caloCellContainer,
		boost::dynamic_bitset<> &	map ) const

overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauPi0CreateROI.

Definition at line 52 of file TauRecToolBase.cxx.

                                                                                                                        {
  ATH_MSG_ERROR("function not implemented");
  return StatusCode::FAILURE;
}

◆ executePi0nPFO()

StatusCode TauRecToolBase::executePi0nPFO	(	xAOD::TauJet &	pTau,
		xAOD::PFOContainer &	neutralPFOContainer ) const

overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauPi0ScoreCalculator, and TauPi0Selector.

Definition at line 95 of file TauRecToolBase.cxx.

                                                                                  {
  ATH_MSG_ERROR("function not implemented");
  return StatusCode::FAILURE;
}

◆ executeShotFinder()

StatusCode TauRecToolBase::executeShotFinder	(	xAOD::TauJet &	pTau,
		xAOD::CaloClusterContainer &	shotClusterContainer,
		xAOD::PFOContainer &	PFOContainer ) const

overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauShotFinder.

Definition at line 73 of file TauRecToolBase.cxx.

                                                                                                                   {
  ATH_MSG_ERROR("function not implemented");
  return StatusCode::FAILURE;
}

◆ executeTrackClassifier()

StatusCode TauRecToolBase::executeTrackClassifier	(	xAOD::TauJet &	pTau,
		xAOD::TauTrackContainer &	tauTrackContainer ) const

overridevirtualinherited

Implements ITauToolBase.

Reimplemented in tauRecTools::TauTrackRNNClassifier.

Definition at line 68 of file TauRecToolBase.cxx.

                                                                                           {
  ATH_MSG_ERROR("function not implemented");
  return StatusCode::FAILURE;
}

◆ executeTrackFinder()

StatusCode TauTrackFinder::executeTrackFinder	(	xAOD::TauJet &	pTau,
		xAOD::TauTrackContainer &	tauTrackCon ) const

overridevirtual

Reimplemented from TauRecToolBase.

Definition at line 53 of file TauTrackFinder.cxx.

                                                                                                        {
  const EventContext &ctx = Gaudi::Hive::currentContext();
  std::vector<const xAOD::TrackParticle*> tauTracks;
  std::vector<const xAOD::TrackParticle*> wideTracks;
  std::vector<const xAOD::TrackParticle*> otherTracks;
 
  //Retrieve standard track container
  const xAOD::TrackParticleContainer* trackParticleCont = nullptr; 
  
  SG::ReadHandle<xAOD::TrackParticleContainer> trackPartInHandle( m_trackPartInputContainer, ctx );
  if (!trackPartInHandle.isValid()) {
    ATH_MSG_ERROR ("Could not retrieve HiveDataObj with key " << trackPartInHandle.key());
    return StatusCode::FAILURE;
  }
  trackParticleCont = trackPartInHandle.cptr();
 
  //Retrieve LRT container
  const xAOD::TrackParticleContainer* largeD0TracksParticleCont = nullptr; 
  std::vector<const xAOD::TrackParticle*> vecTrksLargeD0;
  if (! m_largeD0TracksInputContainer.empty()) { 
    SG::ReadHandle<xAOD::TrackParticleContainer> trackPartInHandle( m_largeD0TracksInputContainer, ctx );
    if (!trackPartInHandle.isValid()) {
      ATH_MSG_VERBOSE ("Could not retrieve HiveDataObj with key " << trackPartInHandle.key());
      ATH_MSG_VERBOSE ("LRT container " << trackPartInHandle.key()<<" is not being used for tau tracks");
    }
    else { 
      largeD0TracksParticleCont = trackPartInHandle.cptr();
      vecTrksLargeD0 = std::vector<const xAOD::TrackParticle*>(largeD0TracksParticleCont->begin(), largeD0TracksParticleCont->end());
    }  
  }
 
  
 
  // retrieve the seed jet container when using ghost-matching
  const xAOD::JetContainer* jetContainer = nullptr; 
  if (! m_jetContainer.empty()) {
    SG::ReadHandle<xAOD::JetContainer> jetContHandle( m_jetContainer, ctx );
    if (!jetContHandle.isValid()) {
      ATH_MSG_ERROR ("Could not retrieve HiveDataObj with key " << jetContHandle.key());
      return StatusCode::FAILURE;
    }
    jetContainer = jetContHandle.cptr();
  }
  // in EleRM reco, we need the original track particles
  std::vector<const xAOD::TrackParticle*> vecTrks;
  vecTrks.reserve( trackParticleCont->size() );
  for (auto trk : *trackParticleCont){
    if (!inEleRM()) { vecTrks.push_back(trk); }
    else{
      static const SG::ConstAccessor<ElementLink<xAOD::TrackParticleContainer>> acc_originalObject("ERMOriginalTrack");
      auto original_id_track_link = acc_originalObject(*trk);
      if (!original_id_track_link.isValid()) {
          ATH_MSG_ERROR("Original track link is not valid");
          continue;
      }
      vecTrks.push_back(*original_id_track_link);
    }
  }
 
  // get the primary vertex
  const xAOD::Vertex* pVertex = pTau.vertex();
 
  // retrieve tracks wrt a vertex                                                                                                                              
  // as a vertex is used: tau origin / PV / beamspot / 0,0,0 (in this order, depending on availability)                                                        
 
  getTauTracksFromPV(pTau, vecTrks, pVertex, m_useGhostTracks, jetContainer, tauTracks, wideTracks, otherTracks);
 
  bool foundLRTCont = bool (largeD0TracksParticleCont != nullptr);
  // additional LRT with vertex association added to tracks
  if (foundLRTCont){
    // for now, use cone association for LRTs, not ghost association
    getTauTracksFromPV(pTau, vecTrksLargeD0, pVertex, false, nullptr, tauTracks, wideTracks, otherTracks);
  }
 
  // remove core and wide tracks outside a maximal delta z0 wrt lead core track                                                                                
  if (m_applyZ0cut) {
    this->removeOffsideTracksWrtLeadTrk(ctx, tauTracks, wideTracks, otherTracks, pVertex, m_z0maxDelta);
  }
 
  if(m_removeDuplicateCoreTracks){
    bool alreadyUsed = false;
    for (std::vector<const xAOD::TrackParticle*>::iterator track_it = tauTracks.begin(); track_it != tauTracks.end() ;)
      {
    alreadyUsed = false;
    //loop over all up-to-now core tracks   
    for( const xAOD::TauTrack* tau_trk : tauTrackCon ) {
      //originally it was coreTrack&passTrkSelector
      if(! tau_trk->flagWithMask( (1<<xAOD::TauJetParameters::TauTrackFlag::coreTrack) | (1<<xAOD::TauJetParameters::TauTrackFlag::passTrkSelector))) continue; 
      if( (*track_it) == tau_trk->track()) alreadyUsed = true;
    }
    //if this track has already been used by another tau, don't associate it to this new one                                                               
    if(alreadyUsed) {
          ATH_MSG_INFO( "Found Already Used track new, now removing: " << *track_it );
      track_it = tauTracks.erase(track_it);
        } else ++track_it;
      }
  }
 
  // associate track to tau candidate and calculate charge                                                                                                    
  float charge = 0.;  
  for (unsigned int i = 0; i < tauTracks.size(); ++i) {
    const xAOD::TrackParticle* trackParticle = tauTracks.at(i);
 
    ATH_MSG_VERBOSE(name() << " adding core track nr: " << i
            << " eta " << trackParticle->eta()
            << " phi " << trackParticle->phi());
 
    charge += trackParticle->charge();
 
    xAOD::TauTrack* track = new xAOD::TauTrack();
    tauTrackCon.push_back(track);
 
    ElementLink<xAOD::TrackParticleContainer> linkToTrackParticle;
    linkToTrackParticle.toContainedElement(*static_cast<const xAOD::TrackParticleContainer*>(trackParticle->container()), trackParticle);
    if (foundLRTCont && isLargeD0Track(trackParticle)){//Check LRT track and link to container
      track->setFlag(xAOD::TauJetParameters::TauTrackFlag::LargeRadiusTrack, true);
    }
    else {
      track->setFlag(xAOD::TauJetParameters::TauTrackFlag::LargeRadiusTrack, false);
    }
    track->addTrackLink(linkToTrackParticle);
 
    track->setP4(trackParticle->pt(), trackParticle->eta(), trackParticle->phi(), trackParticle->m());
    track->setFlag(xAOD::TauJetParameters::TauTrackFlag::coreTrack, true);
    track->setFlag(xAOD::TauJetParameters::TauTrackFlag::passTrkSelector, true);
    // in case TrackClassifier is not run, still get sensible results                                                                                        
    track->setFlag(xAOD::TauJetParameters::TauTrackFlag::classifiedCharged, true);
    track->setFlag(xAOD::TauJetParameters::TauTrackFlag::unclassified, true);
    
    ElementLink<xAOD::TauTrackContainer> linkToTauTrack;
    linkToTauTrack.toContainedElement(tauTrackCon, track);
    pTau.addTauTrackLink(linkToTauTrack);
 
    ATH_MSG_VERBOSE(name() << " added core track nr: " << i
            << " eta " << pTau.track(i)->eta()
            << " phi " << pTau.track(i)->phi());
  }
  // set the charge, which is defined by the core tau tracks only                                                                                              
  pTau.setCharge(charge);
 
  for (unsigned int i = 0; i < wideTracks.size(); ++i) {
    const xAOD::TrackParticle* trackParticle = wideTracks.at(i);
 
    ATH_MSG_VERBOSE(name() << " adding wide track nr: " << i
            << " eta " << trackParticle->eta()
            << " phi " << trackParticle->phi());
 
    xAOD::TauTrack* track = new xAOD::TauTrack();
    tauTrackCon.push_back(track);
 
    ElementLink<xAOD::TrackParticleContainer> linkToTrackParticle;
    linkToTrackParticle.toContainedElement(*static_cast<const xAOD::TrackParticleContainer*>(trackParticle->container()), trackParticle);
    if (foundLRTCont && isLargeD0Track(trackParticle)){//Check LRT track and link to container
      track->setFlag(xAOD::TauJetParameters::TauTrackFlag::LargeRadiusTrack, true);
    }
    else {
      track->setFlag(xAOD::TauJetParameters::TauTrackFlag::LargeRadiusTrack, false);
    }
    track->addTrackLink(linkToTrackParticle);
 
    track->setP4(trackParticle->pt(), trackParticle->eta(), trackParticle->phi(), trackParticle->m());
    track->setFlag(xAOD::TauJetParameters::TauTrackFlag::wideTrack, true);
    track->setFlag(xAOD::TauJetParameters::TauTrackFlag::passTrkSelector, true);
    // in case TrackClassifier is not run, still get sensible results                                                                                        
    track->setFlag(xAOD::TauJetParameters::TauTrackFlag::classifiedIsolation, true); // for sake of trigger, reset in TauTrackClassifier
    track->setFlag(xAOD::TauJetParameters::TauTrackFlag::modifiedIsolationTrack, true); // for sake of trigger, reset in TauTrackClassifier
    track->setFlag(xAOD::TauJetParameters::TauTrackFlag::unclassified, true);
 
    ElementLink<xAOD::TauTrackContainer> linkToTauTrack;
    linkToTauTrack.toContainedElement(tauTrackCon, track);
    pTau.addTauTrackLink(linkToTauTrack);
  }
 
  //These are set again in TauTrackClassifier                                                                                                                  
  pTau.setDetail(xAOD::TauJetParameters::nChargedTracks, static_cast<int>(pTau.nTracks()));
  pTau.setDetail(xAOD::TauJetParameters::nIsolatedTracks, static_cast<int>(pTau.nTracks(xAOD::TauJetParameters::classifiedIsolation)));
 
  for (unsigned int i = 0; i < otherTracks.size(); ++i) {
    const xAOD::TrackParticle* trackParticle = otherTracks.at(i);
 
    ATH_MSG_VERBOSE(name()     << " adding other track nr: " << i
            << " eta " << trackParticle->eta()
            << " phi " << trackParticle->phi());
 
    xAOD::TauTrack* track = new xAOD::TauTrack();
    tauTrackCon.push_back(track);
 
    ElementLink<xAOD::TrackParticleContainer> linkToTrackParticle;
    if (foundLRTCont && isLargeD0Track(trackParticle)){//Check LRT track and link to container
      linkToTrackParticle.toContainedElement(*static_cast<const xAOD::TrackParticleContainer*>(trackParticle->container()), trackParticle);
      track->setFlag(xAOD::TauJetParameters::TauTrackFlag::LargeRadiusTrack, true); 
    }
    else {
      linkToTrackParticle.toContainedElement(*static_cast<const xAOD::TrackParticleContainer*>(trackParticle->container()), trackParticle);
      track->setFlag(xAOD::TauJetParameters::TauTrackFlag::LargeRadiusTrack, false);
    }
    track->addTrackLink(linkToTrackParticle);
 
    track->setP4(trackParticle->pt(), trackParticle->eta(), trackParticle->phi(), trackParticle->m());
    float dR = track->p4().DeltaR(pTau.p4());
    if(dR<=m_maxJetDr_tau) track->setFlag(xAOD::TauJetParameters::TauTrackFlag::coreTrack, true);
    else track->setFlag(xAOD::TauJetParameters::TauTrackFlag::wideTrack, true);
    track->setFlag(xAOD::TauJetParameters::TauTrackFlag::unclassified, true);
 
    ElementLink<xAOD::TauTrackContainer> linkToTauTrack;
    linkToTauTrack.toContainedElement(tauTrackCon, track);
    pTau.addTauTrackLink(linkToTauTrack);
  }
 
  pTau.setDetail(xAOD::TauJetParameters::nLargeRadiusTracks, static_cast<int>(pTau.nTracks(xAOD::TauJetParameters::LargeRadiusTrack)));
  // keep track of total number of associated tracks, in case of tau track thinning
  pTau.setDetail(xAOD::TauJetParameters::nAllTracks, static_cast<int>(pTau.nAllTracks()));
 
  ATH_MSG_DEBUG("numTrack: " << "/" << pTau.nTracks());
  ATH_MSG_DEBUG("charge: " << "/" << pTau.charge());
 
  // impact parameter variables w.r.t. tau vertex 
  const xAOD::Vertex* vxcand = nullptr;
 
  xAOD::Vertex vxbkp;
  vxbkp.makePrivateStore();
 
  // FIXME: don't we want to use the beamspot in the offline reconstruction too, when there is no reconstructed primary vertex?
  if (pTau.vertex()!=nullptr && pTau.vertex()->vertexType() != xAOD::VxType::NoVtx) {
    vxcand = pTau.vertex();
  }
  else if (inTrigger()) { // online: use vertex with x-y coordinates from the beamspot and the z from the leading track
    vxbkp.setX(0); vxbkp.setY(0); vxbkp.setZ(0);
 
    SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle { m_beamSpotKey, ctx };
    if(beamSpotHandle.isValid()) {
      vxbkp.setPosition(beamSpotHandle->beamPos());
      const auto& cov = beamSpotHandle->beamVtx().covariancePosition();
      vxbkp.setCovariancePosition(cov);
 
      if(!tauTracks.empty()) {
         vxbkp.setZ(tauTracks.at(0)->z0());
      }
    }
    else {
      ATH_MSG_DEBUG("No Beamspot object in tau candidate");
    }
    vxcand = & vxbkp;
  }
 
 
  // this could be replaced with TauTrack::setDetail
  static const SG::Accessor<float> dec_d0TJVA("d0TJVA");
  static const SG::Accessor<float> dec_z0sinthetaTJVA("z0sinthetaTJVA");
  static const SG::Accessor<float> dec_d0SigTJVA("d0SigTJVA");
  static const SG::Accessor<float> dec_z0sinthetaSigTJVA("z0sinthetaSigTJVA");
 
  for(const ElementLink<xAOD::TauTrackContainer>& trackLink : pTau.allTauTrackLinks())
  {
    assert (trackLink.getStorableObjectPointer() == &tauTrackCon);
    xAOD::TauTrack* track = tauTrackCon[trackLink.index()];
    dec_d0TJVA(*track) = track->track()->d0();
    dec_z0sinthetaTJVA(*track) = track->z0sinThetaTJVA(pTau);
    dec_d0SigTJVA(*track) = -999.;
    dec_z0sinthetaSigTJVA(*track) = -999.;
 
    // in the trigger, z0sintheta and corresponding significance are meaningless if we use the beamspot
    if(vxcand) {
      std::unique_ptr<const Trk::ImpactParametersAndSigma> myIPandSigma 
    = std::unique_ptr<const Trk::ImpactParametersAndSigma>(m_trackToVertexIPEstimator->estimate(ctx, track->track(), vxcand));
      
      if(myIPandSigma) {
    dec_d0TJVA(*track) = myIPandSigma->IPd0;
    dec_z0sinthetaTJVA(*track) = myIPandSigma->IPz0SinTheta;
    dec_d0SigTJVA(*track) = (myIPandSigma->sigmad0 != 0.) ? static_cast<float>( myIPandSigma->IPd0 / myIPandSigma->sigmad0 ) : -999.f;
    dec_z0sinthetaSigTJVA(*track) = (myIPandSigma->sigmaz0SinTheta != 0.) ? static_cast<float>( myIPandSigma->IPz0SinTheta / myIPandSigma->sigmaz0SinTheta ) : -999.f;
      }
    }
  }
 
  // extrapolate core tracks to calorimeter surface
  // store information only in ExtraDetailsContainer
  if(!m_bypassExtrapolator)
    {
      StatusCode sc = extrapolateToCaloSurface(ctx, pTau, tauTrackCon);
      if (sc.isFailure() && !sc.isRecoverable()) {
    ATH_MSG_ERROR("couldn't extrapolate tracks to calo surface");
    return StatusCode::FAILURE;
      }
    }
  
  return StatusCode::SUCCESS;
}

◆ executeVertexFinder()

StatusCode TauRecToolBase::executeVertexFinder	(	xAOD::TauJet &	pTau,
		const xAOD::VertexContainer *	vertexContainer = nullptr ) const

overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauVertexFinder.

Definition at line 58 of file TauRecToolBase.cxx.

                                                                                            {
  ATH_MSG_ERROR("function not implemented");
  return StatusCode::FAILURE;
}

◆ executeVertexVariables()

StatusCode TauRecToolBase::executeVertexVariables	(	xAOD::TauJet &	pTau,
		xAOD::VertexContainer &	vertexContainer ) const

overridevirtualinherited

Implements ITauToolBase.

Reimplemented in TauVertexVariables.

Definition at line 85 of file TauRecToolBase.cxx.

                                                                                              {
  ATH_MSG_ERROR("function not implemented");
  return StatusCode::FAILURE;
}

◆ 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

◆ extrapolateToCaloSurface()

StatusCode TauTrackFinder::extrapolateToCaloSurface	(	const EventContext &	ctx,
		xAOD::TauJet &	pTau,
		xAOD::TauTrackContainer &	tauTrackCon ) const

private

Extrapolate track eta and phi to the calorimeter middle surface.

Definition at line 437 of file TauTrackFinder.cxx.

{
  Trk::TrackParametersIdHelper parsIdHelper;
 
  int trackIndex = -1;
  const Trk::CaloExtension * caloExtension = nullptr;
  std::unique_ptr<Trk::CaloExtension> uniqueExtension;
  for(const ElementLink<xAOD::TauTrackContainer>& trackLink : pTau.allTauTrackLinks())
  {
    assert (trackLink.getStorableObjectPointer() == &tauTrackCon);
    xAOD::TauTrack* tauTrack = tauTrackCon[trackLink.index()];
    const xAOD::TrackParticle *orgTrack = tauTrack->track();
    if( !orgTrack ) continue;
    trackIndex = orgTrack->index();   
 
    // set default values
    float etaEM = -10.0;
    float phiEM = -10.0;
    float etaHad = -10.0;
    float phiHad = -10.0;
 
    // get the extrapolation into the calo
    ATH_MSG_DEBUG( "Try extrapolation of track with pt = " << orgTrack->pt() 
           << ", eta " << orgTrack->eta() 
           << ", phi" << orgTrack->phi() );
 
    if(!m_ParticleCacheKey.key().empty()){
      /*get the CaloExtension object*/
      ATH_MSG_VERBOSE("Using the CaloExtensionBuilder Cache");
      SG::ReadHandle<CaloExtensionCollection>  particleCache {m_ParticleCacheKey, ctx};
      caloExtension = (*particleCache)[trackIndex];
      ATH_MSG_VERBOSE("Getting element " << trackIndex << " from the particleCache");
      if( not caloExtension ){
        ATH_MSG_VERBOSE("Cache does not contain a calo extension -> "
                        "Calculating with the a CaloExtensionTool");
        uniqueExtension = m_caloExtensionTool->caloExtension(ctx, *orgTrack);
        caloExtension = uniqueExtension.get();
      }
    }
    else {
      /* If CaloExtensionBuilder is unavailable, use the calo extension tool */
      ATH_MSG_VERBOSE("Using the CaloExtensionTool");
      uniqueExtension = m_caloExtensionTool->caloExtension(ctx, *orgTrack);
      caloExtension = uniqueExtension.get();
    }
 
    if (!caloExtension)
      { 
    ATH_MSG_DEBUG("Track extrapolation failed");
      }
    else {
      const std::vector<Trk::CurvilinearParameters>& clParametersVector = caloExtension->caloLayerIntersections();
      if (clParametersVector.empty()) {
    ATH_MSG_DEBUG("Track extrapolation failed");
      }
 
      ATH_MSG_DEBUG("Scanning samplings");
      bool validECal = false;
      bool validHCal = false;
      for( const Trk::CurvilinearParameters& cur : clParametersVector ){
    ATH_MSG_DEBUG("Sampling " << parsIdHelper.caloSample(cur.cIdentifier()) );
                
    // only use entry layer
    if( not parsIdHelper.isEntryToVolume(cur.cIdentifier()) ) continue;
 
    CaloSampling::CaloSample sample = parsIdHelper.caloSample(cur.cIdentifier());
 
    // ECal
    if( not validECal and m_EMSamplings.count(sample))
      {
        validECal = true;
        etaEM = cur.position().eta();
        phiEM = cur.position().phi();
        ATH_MSG_DEBUG("Extrapolated to ECal layer " << sample);
      }
 
    // HCal
    if( not validHCal and m_HadSamplings.count(sample))
      {
        validHCal = true;
        etaHad = cur.position().eta();
        phiHad = cur.position().phi();
        ATH_MSG_DEBUG("Extrapolated to HCal layer " << sample);
      }
    if( validECal and validHCal ) break;
      }
      // EM failure warn if within acceptance 
      if( not validECal and std::abs(orgTrack->pt()) < 2.48 ){
    ATH_MSG_DEBUG("Failed extrapolation to ECal");
      }
      // Had failure warn if enough pt to reach HCal
      if( not validHCal and orgTrack->pt() > 2000. ){
    ATH_MSG_DEBUG("Failed extrapolation to HCal");
      }
 
      ATH_MSG_DEBUG( "Extrapolated track with eta=" << orgTrack->eta()
             << " phi="<<orgTrack->phi()
             << " to ECal eta=" << etaEM 
             << " phi="<< phiEM
             << " HCal eta=" << etaHad 
             << " phi="<< phiHad
             );
    }
    tauTrack->setDetail(xAOD::TauJetParameters::CaloSamplingEtaEM, etaEM);
    tauTrack->setDetail(xAOD::TauJetParameters::CaloSamplingPhiEM, phiEM);
    tauTrack->setDetail(xAOD::TauJetParameters::CaloSamplingEtaHad, etaHad);
    tauTrack->setDetail(xAOD::TauJetParameters::CaloSamplingPhiHad, phiHad);
  }
 
  return StatusCode::SUCCESS;
 
}

◆ finalize()

StatusCode TauRecToolBase::finalize ( )

overridevirtualinherited

Finalizer.

Implements ITauToolBase.

Reimplemented in PanTau::PanTauProcessor.

Definition at line 109 of file TauRecToolBase.cxx.

                                   {
  return StatusCode::SUCCESS;
}

◆ find_file()

std::string TauRecToolBase::find_file ( const std::string & fname ) const

inherited

Definition at line 19 of file TauRecToolBase.cxx.

                                                                {
  std::string full_path;
  //offline calib files are in GroupData
  //online calib files are in release
  full_path = PathResolverFindCalibFile(m_tauRecToolsTag+"/"+fname);
  if(full_path.empty()) full_path = PathResolverFindCalibFile(fname);
  return full_path;
}

◆ getKey()

SG::sgkey_t asg::AsgTool::getKey ( const void * ptr ) const

inherited

Get the (hashed) key of an object that is in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::Event both provide ways for getting the SG::sgkey_t key for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also: asg::AsgTool::getName

Parameters

ptr	The bare pointer to the object that the event store should know about

Returns: The hashed key of the object in the store. If not found, an invalid (zero) key.

Definition at line 119 of file AsgTool.cxx.

                                                    {
 
#ifdef XAOD_STANDALONE
      // In case we use @c xAOD::Event, we have a direct function call
      // for this.
      return evtStore()->event()->getKey( ptr );
#else
      const SG::DataProxy* proxy = evtStore()->proxy( ptr );
      return ( proxy == nullptr ? 0 : proxy->sgkey() );
#endif // XAOD_STANDALONE
   }

◆ getName()

const std::string & asg::AsgTool::getName ( const void * ptr ) const

inherited

Get the name of an object that is / should be in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::Event both provide ways for getting the std::string name for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also: asg::AsgTool::getKey

Parameters

ptr	The bare pointer to the object that the event store should know about

Returns: The string name of the object in the store. If not found, an empty string.

Definition at line 106 of file AsgTool.cxx.

                                                            {
 
#ifdef XAOD_STANDALONE
      // In case we use @c xAOD::Event, we have a direct function call
      // for this.
      return evtStore()->event()->getName( ptr );
#else
      const SG::DataProxy* proxy = evtStore()->proxy( ptr );
      static const std::string dummy = "";
      return ( proxy == nullptr ? dummy : proxy->name() );
#endif // XAOD_STANDALONE
   }

◆ getProperty()

template<class T>

const T * asg::AsgTool::getProperty ( const std::string & name ) const

inherited

Get one of the tool's properties.

◆ getTauTracksFromPV()

void TauTrackFinder::getTauTracksFromPV	(	const xAOD::TauJet &	tauJet,
		const std::vector< const xAOD::TrackParticle * > &	vecTrackParticles,
		const xAOD::Vertex *	primaryVertex,
		const bool &	useGhostTracks,
		const xAOD::JetContainer *	jetContainer,
		std::vector< const xAOD::TrackParticle * > &	tauTracks,
		std::vector< const xAOD::TrackParticle * > &	wideTracks,
		std::vector< const xAOD::TrackParticle * > &	otherTracks ) const

private

Definition at line 364 of file TauTrackFinder.cxx.

{
  std::vector<const xAOD::TrackParticle*> ghostTracks;
  if(useGhostTracks) {
    const xAOD::Jet* seedJet = pTau.jet();
    if(seedJet) {
      if(!seedJet->getAssociatedObjects("GhostTrack", ghostTracks)) {
    ATH_MSG_WARNING("Could not retrieve GhostTrack from seed jet.");
      }
    }
  }
  // in EleRM reco, we need the original track particles
  if (inEleRM()){
    static const SG::ConstAccessor<ElementLink<xAOD::TrackParticleContainer>> acc_originalTrack("ERMOriginalTrack");
    for (uint i = 0; i < ghostTracks.size(); i++){
      auto original_id_track_link = acc_originalTrack(*(ghostTracks[i]));
      if (!original_id_track_link.isValid()) {
          ATH_MSG_ERROR("Original track link is not valid");
          continue;
      }
      ghostTracks[i] = *original_id_track_link;
    }
  }
  
  for (const xAOD::TrackParticle *trackParticle : vecTrackParticles) {
    TauTrackType type = tauTrackType(pTau, *trackParticle, primaryVertex);
    if(type == NotTauTrack) continue;
 
    if(useGhostTracks) {
      // require that tracks are ghost-matched with the seed jet at large dR(tau,track), to avoid using tracks from another tau
      double dR = pTau.p4().DeltaR(trackParticle->p4());
      if (dR > m_ghostTrackDR) {
    if (std::find(ghostTracks.begin(), ghostTracks.end(), trackParticle) == ghostTracks.end()) {
      // check whether the jet closest to the track is the current seed jet
      // if so, recover the track even if not ghost-matched, to improve tau-track association efficiency at low pt (esp. for 3p)
      double dRmin = 999.;
      bool isSeedClosest = false;
      for (const xAOD::Jet* jet : *jetContainer) {
        xAOD::JetFourMom_t jetP4 = jet->jetP4(xAOD::JetScale::JetConstitScaleMomentum);
        TLorentzVector jetLV;
        jetLV.SetPtEtaPhiM(jetP4.Pt(), jetP4.Eta(), jetP4.Phi(), jetP4.M());
        double dRjet = trackParticle->p4().DeltaR(jetLV);
        if(dRjet < dRmin) {
          dRmin = dRjet;
          isSeedClosest = (jet == pTau.jet());
        }
      }
      if(!isSeedClosest) continue;
    }
      }
    }
    
    if (type == TauTrackCore)
      tauTracks.push_back(trackParticle);
    else if (type == TauTrackWide)
      wideTracks.push_back(trackParticle);
    else if (type == TauTrackOther)
      otherTracks.push_back(trackParticle);
  }
  std::sort(tauTracks.begin(), tauTracks.end(), TrackSort());
  std::sort(wideTracks.begin(), wideTracks.end(), TrackSort());
  std::sort(otherTracks.begin(), otherTracks.end(), TrackSort());
}

◆ getZ0()

float TauTrackFinder::getZ0	(	const EventContext &	ctx,
		const xAOD::TrackParticle *	track,
		const xAOD::Vertex *	vertex ) const

private

Some internally used functions.

Definition at line 612 of file TauTrackFinder.cxx.

{
  float MAX=1e5;
 
  if (!track) return MAX;
 
  std::unique_ptr<Trk::Perigee> perigee;
  if (vertex) perigee = m_trackToVertexTool->perigeeAtVertex(ctx, *track, vertex->position());
  else        perigee = m_trackToVertexTool->perigeeAtVertex(ctx, *track); //will use beamspot or 0,0,0 instead
 
  if (!perigee) {
    ATH_MSG_WARNING("Bad track; can't find perigee at vertex.");
    return MAX;
  }
 
  return static_cast<float>(perigee->parameters()[Trk::z0]);
}

◆ inAOD()

bool TauRecToolBase::inAOD ( ) const

inlineprotectedinherited

Definition at line 87 of file TauRecToolBase.h.

87{ return m_in_AOD; }

TauRecToolBase::m_in_AOD

Gaudi::Property< bool > m_in_AOD

Definition TauRecToolBase.h:76

◆ inEleRM()

bool TauRecToolBase::inEleRM ( ) const

inlineprotectedinherited

Definition at line 88 of file TauRecToolBase.h.

88{ return m_in_EleRM; }

TauRecToolBase::m_in_EleRM

Gaudi::Property< bool > m_in_EleRM

Definition TauRecToolBase.h:77

◆ initialize()

StatusCode TauTrackFinder::initialize ( void )

overridevirtual

Algorithm functions.

Reimplemented from TauRecToolBase.

Definition at line 22 of file TauTrackFinder.cxx.

                                      {
 
  // retrieve tools
  ATH_CHECK( m_trackSelectorTool_tau.retrieve() );
  ATH_CHECK( m_trackToVertexTool.retrieve() );
  ATH_CHECK( m_caloExtensionTool.retrieve() );
  ATH_CHECK( m_trackToVertexIPEstimator.retrieve() );
 
  // initialize ReadHandleKey
  ATH_CHECK( m_trackPartInputContainer.initialize() );
  // use CaloExtensionTool when key is empty 
  ATH_CHECK( m_ParticleCacheKey.initialize(SG::AllowEmpty) );
  // allow empty for LRT
  ATH_CHECK( m_largeD0TracksInputContainer.initialize(SG::AllowEmpty) );
 
  if(m_useGhostTracks) {
    if(inTrigger()) {
      ATH_MSG_ERROR ("Ghost matching is not a valid tau-track association scheme for trigger, use cone association. Aborting.");
      return StatusCode::FAILURE;
    }
    ATH_MSG_INFO ("Using ghost matching for tau-track association" << m_ghostTrackDR );
    // allow empty for trigger 
    ATH_CHECK( m_jetContainer.initialize(SG::AllowEmpty) );
  }
 
  ATH_CHECK( m_beamSpotKey.initialize(inTrigger()) );
 
  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.

◆ inTrigger()

bool TauRecToolBase::inTrigger ( ) const

inlineprotectedinherited

Definition at line 86 of file TauRecToolBase.h.

86{ return m_in_trigger; }

TauRecToolBase::m_in_trigger

Gaudi::Property< bool > m_in_trigger

Definition TauRecToolBase.h:75

◆ isLargeD0Track()

bool TauTrackFinder::isLargeD0Track ( const xAOD::TrackParticle * track ) const

private

Definition at line 630 of file TauTrackFinder.cxx.

{
  const std::bitset<xAOD::NumberOfTrackRecoInfo> patternReco = track->patternRecoInfo();
  if (patternReco.test(xAOD::TrackPatternRecoInfo::SiSpacePointsSeedMaker_LargeD0)) {
    ATH_MSG_DEBUG("LargeD0Track found");
    return true;
  }
 
  return false;
}

◆ msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msg_level_name()

const std::string & asg::AsgTool::msg_level_name ( ) const

inherited

A deprecated function for getting the message level's name.

Instead of using this, weirdly named function, user code should get the string name of the current minimum message level (in case they really need it...), with:

MSG::name( msg().level() )

This function's name doesn't follow the ATLAS coding rules, and as such will be removed in the not too distant future.

Returns: The string name of the current minimum message level that's printed

Definition at line 101 of file AsgTool.cxx.

                                                  {
 
      return MSG::name( msg().level() );
   }

◆ 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.

◆ print()

void asg::AsgTool::print ( ) const

virtualinherited

Print the state of the tool.

Implements asg::IAsgTool.

Reimplemented in AsgHelloTool, HI::HIPileupTool, JetBottomUpSoftDrop, JetConstituentsRetriever, JetDumper, JetFinder, JetFromPseudojet, JetModifiedMassDrop, JetPileupLabelingTool, JetPruner, JetPseudojetRetriever, JetReclusterer, JetReclusteringTool, JetRecTool, JetRecursiveSoftDrop, JetSoftDrop, JetSplitter, JetSubStructureMomentToolsBase, JetToolRunner, JetTrimmer, JetTruthLabelingTool, and KtDeltaRTool.

Definition at line 131 of file AsgTool.cxx.

                             {
 
      ATH_MSG_INFO( "AsgTool " << name() << " @ " << this );
      return;
   }

◆ removeOffsideTracksWrtLeadTrk()

void TauTrackFinder::removeOffsideTracksWrtLeadTrk	(	const EventContext &	ctx,
		std::vector< const xAOD::TrackParticle * > &	tauTracks,
		std::vector< const xAOD::TrackParticle * > &	wideTracks,
		std::vector< const xAOD::TrackParticle * > &	otherTracks,
		const xAOD::Vertex *	tauOrigin,
		double	maxDeltaZ0 ) const

private

Definition at line 553 of file TauTrackFinder.cxx.

{
  float MAX=1e5;
 
  // need at least one core track to have a leading trk to compare with
  if (tauTracks.empty()) return;
 
  // get lead trk parameters
  const xAOD::TrackParticle *leadTrack = tauTracks.at(0);
  float z0_leadTrk = getZ0(ctx, leadTrack, tauOrigin);
 
  if (z0_leadTrk > MAX-1) return; // bad lead trk -> do nothing
 
  ATH_MSG_VERBOSE("before z0 cut: #coreTracks=" << tauTracks.size() << ", #wideTracks=" << wideTracks.size() << ", #otherTracks=" << otherTracks.size());
 
  std::vector<const xAOD::TrackParticle*>::iterator itr;
    
  // skip leading track, because it is the reference
  itr = tauTracks.begin()+1;
  while (itr!=tauTracks.end()) {
    float z0 = getZ0(ctx, *itr, tauOrigin);
    float deltaZ0=z0 - z0_leadTrk;
    ATH_MSG_VERBOSE("core Trks: deltaZ0= " << deltaZ0);
 
    if ( std::abs(deltaZ0) < maxDeltaZ0 ) {++itr;}
    else {
      otherTracks.push_back(*itr);
      itr = tauTracks.erase(itr); //remove from core track collection
    }
  }
 
  // check wide tracks
  itr = wideTracks.begin();
  while (itr!=wideTracks.end()) {
    float z0 = getZ0(ctx, *itr, tauOrigin);
    float deltaZ0=z0 - z0_leadTrk;
    ATH_MSG_VERBOSE("wide Trks: deltaZ0= " << deltaZ0);
 
    if ( std::abs(deltaZ0) < maxDeltaZ0 ) { ++itr; }
    else {
      otherTracks.push_back(*itr);
      itr = wideTracks.erase(itr); //remove from wide track collection
    }
  }
 
  ATH_MSG_VERBOSE("after z0 cut: #coreTracks=" << tauTracks.size() << ", #wideTracks=" << wideTracks.size() << ", #otherTracks=" << otherTracks.size());
 
  // sort again
  std::sort(tauTracks.begin(), tauTracks.end(), TrackSort());
  std::sort(wideTracks.begin(), wideTracks.end(), TrackSort());
  std::sort(otherTracks.begin(), otherTracks.end(), TrackSort());
}

◆ 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();
  }

◆ 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 >, and AthCheckedComponent<::AthAlgTool >.

◆ 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.

◆ tauTrackType()

TauTrackFinder::TauTrackType TauTrackFinder::tauTrackType	(	const xAOD::TauJet &	tauJet,
		const xAOD::TrackParticle &	trackParticle,
		const xAOD::Vertex *	primaryVertex ) const

private

Definition at line 344 of file TauTrackFinder.cxx.

{
  double dR = pTau.p4().DeltaR(trackParticle.p4());
 
  if (dR > m_maxJetDr_wide) return NotTauTrack;
 
  bool goodTrack = m_trackSelectorTool_tau->decision(trackParticle, primaryVertex);
    
  if (goodTrack) {
    if (dR > m_maxJetDr_tau)
      return TauTrackWide;
    else
      return TauTrackCore;
  } else
    return TauTrackOther;
}

◆ 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_applyZ0cut

Gaudi::Property<bool> TauTrackFinder::m_applyZ0cut {this, "removeTracksOutsideZ0wrtLeadTrk", false}

private

Definition at line 110 of file TauTrackFinder.h.

110{this, "removeTracksOutsideZ0wrtLeadTrk", false};

◆ m_beamSpotKey

SG::ReadCondHandleKey<InDet::BeamSpotData> TauTrackFinder::m_beamSpotKey { this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot" }

private

Definition at line 122 of file TauTrackFinder.h.

122{ this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot" };

◆ m_bypassExtrapolator

Gaudi::Property<bool> TauTrackFinder::m_bypassExtrapolator {this, "BypassExtrapolator", false}

private

Definition at line 113 of file TauTrackFinder.h.

113{this, "BypassExtrapolator", false};

◆ m_caloExtensionTool

ToolHandle<Trk::IParticleCaloExtensionTool> TauTrackFinder::m_caloExtensionTool {this, "ParticleCaloExtensionTool", "Trk::ParticleCaloExtensionTool/ParticleCaloExtensionTool", "Tool for the extrapolation of charged tracks"}

private

tools

Definition at line 103 of file TauTrackFinder.h.

103{this, "ParticleCaloExtensionTool", "Trk::ParticleCaloExtensionTool/ParticleCaloExtensionTool", "Tool for the extrapolation of charged tracks"};

◆ 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_EMSamplings

std::set<CaloSampling::CaloSample> TauTrackFinder::m_EMSamplings

private

Definition at line 124 of file TauTrackFinder.h.

◆ 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_ghostTrackDR

Gaudi::Property<double> TauTrackFinder::m_ghostTrackDR {this, "ghostTrackDR", 0.25}

private

Definition at line 115 of file TauTrackFinder.h.

115{this, "ghostTrackDR", 0.25};

◆ m_HadSamplings

std::set<CaloSampling::CaloSample> TauTrackFinder::m_HadSamplings

private

Definition at line 125 of file TauTrackFinder.h.

◆ m_in_AOD

Gaudi::Property<bool> TauRecToolBase::m_in_AOD {this, "inAOD", false, "Indicate if the tool is running on AOD"}

protectedinherited

Definition at line 76 of file TauRecToolBase.h.

76{this, "inAOD", false, "Indicate if the tool is running on AOD"};

◆ m_in_EleRM

Gaudi::Property<bool> TauRecToolBase::m_in_EleRM {this, "inEleRM", false, "Indicate if the tool is running on EleRM routine"}

protectedinherited

Definition at line 77 of file TauRecToolBase.h.

77{this, "inEleRM", false, "Indicate if the tool is running on EleRM routine"};

◆ m_in_trigger

Gaudi::Property<bool> TauRecToolBase::m_in_trigger {this, "inTrigger", false, "Indicate if the tool is running on trigger"}

protectedinherited

Definition at line 75 of file TauRecToolBase.h.

75{this, "inTrigger", false, "Indicate if the tool is running on trigger"};

◆ m_jetContainer

SG::ReadHandleKey<xAOD::JetContainer> TauTrackFinder::m_jetContainer {this,"Key_jetContainer", "", "Name of the seed jet container, when using ghost matching"}

private

Definition at line 119 of file TauTrackFinder.h.

119{this,"Key_jetContainer", "", "Name of the seed jet container, when using ghost matching"};

◆ m_largeD0TracksInputContainer

SG::ReadHandleKey<xAOD::TrackParticleContainer> TauTrackFinder::m_largeD0TracksInputContainer {this,"Key_LargeD0TrackInputContainer", "", "input LRT particle container key"}

private

Definition at line 118 of file TauTrackFinder.h.

118{this,"Key_LargeD0TrackInputContainer", "", "input LRT particle container key"}; //Expecting InDetLargeD0TrackParticles (offline tracks) if using LRT used

◆ m_maxJetDr_tau

Gaudi::Property<double> TauTrackFinder::m_maxJetDr_tau {this, "MaxJetDrTau", 0.2}

private

Definition at line 108 of file TauTrackFinder.h.

108{this, "MaxJetDrTau", 0.2};

◆ m_maxJetDr_wide

Gaudi::Property<double> TauTrackFinder::m_maxJetDr_wide {this, "MaxJetDrWide", 0.4}

private

Definition at line 109 of file TauTrackFinder.h.

109{this, "MaxJetDrWide", 0.4};

◆ m_ParticleCacheKey

SG::ReadHandleKey<CaloExtensionCollection> TauTrackFinder::m_ParticleCacheKey {this,"tauParticleCache", "ParticleCaloExtension", "Name of the particle measurement extrapolation cache for TauTrackFinder"}

private

Definition at line 120 of file TauTrackFinder.h.

120{this,"tauParticleCache", "ParticleCaloExtension", "Name of the particle measurement extrapolation cache for TauTrackFinder"};

◆ m_removeDuplicateCoreTracks

Gaudi::Property<bool> TauTrackFinder::m_removeDuplicateCoreTracks {this, "removeDuplicateCoreTracks", true}

private

Definition at line 112 of file TauTrackFinder.h.

112{this, "removeDuplicateCoreTracks", true};

◆ m_tauRecToolsTag

Gaudi::Property<std::string> TauRecToolBase::m_tauRecToolsTag {this, "calibFolder", "tauRecTools/R22_preprod", "CVMFS path to the tau calibration folder"}

protectedinherited

Definition at line 78 of file TauRecToolBase.h.

78{this, "calibFolder", "tauRecTools/R22_preprod", "CVMFS path to the tau calibration folder"};

◆ m_trackPartInputContainer

SG::ReadHandleKey<xAOD::TrackParticleContainer> TauTrackFinder::m_trackPartInputContainer {this,"Key_trackPartInputContainer", "InDetTrackParticles", "input track particle container key"}

private

Definition at line 117 of file TauTrackFinder.h.

117{this,"Key_trackPartInputContainer", "InDetTrackParticles", "input track particle container key"};

◆ m_trackSelectorTool_tau

ToolHandle<Trk::ITrackSelectorTool> TauTrackFinder::m_trackSelectorTool_tau {this, "TrackSelectorToolTau", "", "Tool for track selection"}

private

Definition at line 104 of file TauTrackFinder.h.

104{this, "TrackSelectorToolTau", "", "Tool for track selection"};

◆ m_trackToVertexIPEstimator

ToolHandle<Trk::ITrackToVertexIPEstimator> TauTrackFinder::m_trackToVertexIPEstimator {this, "TrackToVertexIPEstimator", ""}

private

Definition at line 106 of file TauTrackFinder.h.

106{this, "TrackToVertexIPEstimator", ""};

◆ m_trackToVertexTool

ToolHandle<Reco::ITrackToVertex> TauTrackFinder::m_trackToVertexTool {this, "TrackToVertexTool", "Reco::TrackToVertex"}

private

Definition at line 105 of file TauTrackFinder.h.

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

◆ m_useGhostTracks

Gaudi::Property<bool> TauTrackFinder::m_useGhostTracks {this, "useGhostTracks", false}

private

Definition at line 114 of file TauTrackFinder.h.

114{this, "useGhostTracks", false};

◆ 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_z0maxDelta

Gaudi::Property<float> TauTrackFinder::m_z0maxDelta {this, "maxDeltaZ0wrtLeadTrk", 1000.}

private

Definition at line 111 of file TauTrackFinder.h.

111{this, "maxDeltaZ0wrtLeadTrk", 1000.};

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

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Member Enumeration Documentation

◆ TauTrackType

Constructor & Destructor Documentation

◆ TauTrackFinder()

◆ ~TauTrackFinder()

Member Function Documentation

◆ ASG_TOOL_CLASS2()

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ eventFinalize()

◆ eventInitialize()

◆ evtStore()

◆ execute()

◆ executePanTau()

◆ executePi0ClusterCreator()

◆ executePi0ClusterScaler()

◆ executePi0CreateROI()

◆ executePi0nPFO()

◆ executeShotFinder()

◆ executeTrackClassifier()

◆ executeTrackFinder()

◆ executeVertexFinder()

◆ executeVertexVariables()

◆ extraDeps_update_handler()

◆ extrapolateToCaloSurface()

◆ finalize()

◆ find_file()

◆ getKey()

◆ getName()

◆ getProperty()

◆ getTauTracksFromPV()

◆ getZ0()

◆ inAOD()

◆ inEleRM()

◆ initialize()

◆ inputHandles()

◆ inTrigger()

◆ isLargeD0Track()

◆ msg()

◆ msg_level_name()

◆ msgLvl()

◆ outputHandles()

◆ print()

◆ removeOffsideTracksWrtLeadTrk()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ tauTrackType()

◆ updateVHKA()

Member Data Documentation

◆ m_applyZ0cut

◆ m_beamSpotKey

◆ m_bypassExtrapolator

◆ m_caloExtensionTool

◆ m_detStore

◆ m_EMSamplings

◆ m_evtStore

◆ m_ghostTrackDR

◆ m_HadSamplings

◆ m_in_AOD

◆ m_in_EleRM

◆ m_in_trigger

◆ m_jetContainer

◆ m_largeD0TracksInputContainer

◆ m_maxJetDr_tau

◆ m_maxJetDr_wide

◆ m_ParticleCacheKey

◆ m_removeDuplicateCoreTracks