TrigTauRecMerged Class Reference

#include <TrigTauRecMerged.h>

Inheritance diagram for TrigTauRecMerged:

Public Member Functions
	TrigTauRecMerged (const std::string &name, ISvcLocator *pSvcLocator)
virtual StatusCode	initialize () override
virtual StatusCode	execute (const EventContext &ctx) const override
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
virtual void	setFilterPassed (bool state, const EventContext &ctx) 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	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

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	Mode { FromClusters =0 , FromTauJet =1 }
enum	TAUEFCALOMON {   NoROIDescr =0 , NoCellCont =1 , EmptyCellCont =2 , NoClustCont =3 ,   NoClustKey =4 , EmptyClustCont =5 , NoJetAttach =6 , NoHLTtauAttach =7 ,   NoHLTtauDetAttach =8 , NoHLTtauXdetAttach =9 }
enum	TAUEFTRKMON { NoTrkCont =0 , NoVtxCont =1 }
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
template<class V>
StatusCode	deepCopy (SG::WriteHandle< DataVector< V > > &writeHandle, const DataVector< V > *oldContainer) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
const ToolHandleArray< ITauToolBase >	m_commonTools {this, "CommonTools", {}, "List of ITauToolBase common tools"}
const ToolHandleArray< ITauToolBase >	m_commonToolsBeforeTF {this, "CommonToolsBeforeTF", {}, "List of ITauToolBase common tools to execute before the Track Finder tools"}
const ToolHandleArray< ITauToolBase >	m_vertexFinderTools {this, "VertexFinderTools", {}, "Vertex Finder tools"}
const ToolHandleArray< ITauToolBase >	m_trackFinderTools {this, "TrackFinderTools", {}, "Track Finder tools"}
const ToolHandleArray< ITauToolBase >	m_vertexVarsTools {this, "VertexVarsTools", {}, "Vertex Variables tools"}
const ToolHandleArray< ITauToolBase >	m_idTools {this, "IDTools", {}, "Vertex Variables tools"}
const ToolHandle< GenericMonitoringTool >	m_monTool {this, "MonTool", "", "Monitoring tool"}
Gaudi::Property< std::map< std::string, std::pair< std::string, std::string > > >	m_monitoredIdScores {this, "MonitoredIDScores", {}, "Pairs of the TauID score and signal-transformed scores for each TauID algorithm to be monitored"}
Gaudi::Property< std::map< std::string, std::pair< std::string, std::string > > >	m_monitoredHitZRegressions {this, "MonitoredHitZRegressions", {}, "Pairs of the z and sigma regression output variables for each HitZ algorithm to be monitored"}
std::map< std::string, std::pair< SG::ConstAccessor< float >, SG::ConstAccessor< float > > >	m_monitoredInferenceAccessors
SG::ReadHandleKey< TrigRoiDescriptorCollection >	m_roiInputKey {this, "InputRoIs", "", "Input RoI name"}
SG::ReadHandleKey< xAOD::CaloClusterContainer >	m_clustersInputKey {this, "InputCaloClusterContainer", "", "Caloclusters in view"}
SG::ReadHandleKey< xAOD::VertexContainer >	m_vertexInputKey {this, "InputVertexContainer", "", "Input vertex container"}
SG::ReadHandleKey< xAOD::TauJetContainer >	m_tauJetInputKey {this, "InputTauJetContainer", "", "Input TauJet container"}
SG::ReadDecorHandleKeyArray< xAOD::TauJetContainer >	m_tauJetInputDecorKeysArray {this, "InputTauJetCopyDecorKeys", {}, "Array of input TauJet decoration keys to copy on the output TauJet"}
SG::ReadHandleKey< xAOD::TauTrackContainer >	m_tauTrackInputKey {this, "InputTauTrackContainer", "", "Input TauTrack container" }
SG::ReadDecorHandleKey< xAOD::TauJetContainer >	m_hitsInputDecorKey {this, "InputTauJetHitsKey", "", "Input TauJet hits decoration key"}
SG::WriteHandleKey< xAOD::JetContainer >	m_tauSeedOutputKey {this, "OutputJetSeed", "", "Output jets which are seeds for tau jets"}
SG::WriteHandleKey< xAOD::TauJetContainer >	m_tauJetOutputKey {this, "OutputTauJetContainer", "", "Output TauJet container"}
SG::WriteDecorHandleKeyArray< xAOD::TauJetContainer >	m_tauJetOutputDecorKeysArray
SG::WriteHandleKey< xAOD::TauTrackContainer >	m_tauTrackOutputKey {this, "OutputTauTrackContainer", "", "Output TauTrack container"}
SG::WriteDecorHandleKey< xAOD::TauJetContainer >	m_hitsOutputDecorKey {this, "OutputTauJetHitsKey", "", "Output TauJet hits decoration key"}
Mode	m_reco_mode = FromClusters
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
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

Definition at line 31 of file TrigTauRecMerged.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

Mode

enum TrigTauRecMerged::Mode

private

Enumerator
FromClusters
FromTauJet

Definition at line 42 of file TrigTauRecMerged.h.

              {
        FromClusters=0, // Calo reconstruction from calo clusters: CaloMVA step
        FromTauJet=1, // Full tau reconstruction from existing tau jets + tracks: CaloHits and Precision steps
    };

TAUEFCALOMON

enum TrigTauRecMerged::TAUEFCALOMON

private

Enumerator
NoROIDescr
NoCellCont
EmptyCellCont
NoClustCont
NoClustKey
EmptyClustCont
NoJetAttach
NoHLTtauAttach
NoHLTtauDetAttach
NoHLTtauXdetAttach

Definition at line 48 of file TrigTauRecMerged.h.

                      {
        NoROIDescr=0,
        NoCellCont=1,
        EmptyCellCont=2,
        NoClustCont=3,
        NoClustKey=4,
        EmptyClustCont=5,
        NoJetAttach=6,
        NoHLTtauAttach=7,
        NoHLTtauDetAttach=8,
        NoHLTtauXdetAttach=9
    };

TAUEFTRKMON

enum TrigTauRecMerged::TAUEFTRKMON

private

Enumerator
NoTrkCont
NoVtxCont

Definition at line 62 of file TrigTauRecMerged.h.

                     {
        NoTrkCont=0,
        NoVtxCont=1
    };

Constructor & Destructor Documentation

TrigTauRecMerged()

TrigTauRecMerged::TrigTauRecMerged	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Definition at line 40 of file TrigTauRecMerged.cxx.

    : AthReentrantAlgorithm(name, pSvcLocator)
{
    
}

Member Function Documentation

cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const

overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

{
  return 0;
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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());
  }

deepCopy()

template<class V>

StatusCode TrigTauRecMerged::deepCopy ( SG::WriteHandle< DataVector< V > > & writeHandle, const DataVector< V > * oldContainer ) const

private

Definition at line 105 of file TrigTauRecMerged.h.

                                                                                                                                          {
   if(!writeHandle.isValid()) {
        ATH_MSG_FATAL("Provided with an invalid write handle");
        return StatusCode::FAILURE;
   }
 
   if(oldContainer){
       for(const V* v : *oldContainer) {
           // Build a new object in the new container
           writeHandle->push_back(std::make_unique<V>());
           // Copy across the aux store
           *writeHandle->back() = *v;
       }
   }
   return StatusCode::SUCCESS;
}

detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::evtStore ( )

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode TrigTauRecMerged::execute ( const EventContext & ctx ) const

overridevirtual

Definition at line 134 of file TrigTauRecMerged.cxx.

{
    //===============================================================================
    // Initialize monitoring variables
    //===============================================================================
 
    // Common CaloOnly and Precision monitored variables:
    auto n_taus                     = Monitored::Scalar<int>("NTauCandidates", 0);
 
    auto pT                         = Monitored::Scalar<float>("Pt", 0);
    auto eta                        = Monitored::Scalar<float>("Eta", -99.9);
    auto phi                        = Monitored::Scalar<float>("Phi", -99.9);
 
    auto etaRoI                     = Monitored::Scalar<float>("EtaRoI", -99.9);
    auto phiRoI                     = Monitored::Scalar<float>("PhiRoI", -99.9);
    auto dEta_RoI                   = Monitored::Scalar<float>("dEtaTau_RoI", -10);
    auto dPhi_RoI                   = Monitored::Scalar<float>("dPhiTau_RoI", -10);
 
    auto mEflowApprox               = Monitored::Scalar<float>("mEflowApprox", -99.9);
    auto ptRatioEflowApprox         = Monitored::Scalar<float>("ptRatioEflowApprox", -99.9);
    auto pt_jetseed_log             = Monitored::Scalar<float>("pt_jetseed_log", -99.9);
    auto etaDetectorAxis            = Monitored::Scalar<float>("etaDetectorAxis", -99.9);
    auto ptDetectorAxis             = Monitored::Scalar<float>("ptDetectorAxis", -99.9);
    auto ptDetectorAxis_log         = Monitored::Scalar<float>("ptDetectorAxis_log", -99.9);
 
    auto n_cells                    = Monitored::Scalar<int>("NCaloCells", 0);
    auto EMRadius                   = Monitored::Scalar<float>("EMRadius", -0.099);
    auto HadRadius                  = Monitored::Scalar<float>("HadRadius", -0.099);
    auto EtHad                      = Monitored::Scalar<float>("EtHad", -10);
    auto EtEm                       = Monitored::Scalar<float>("EtEm", -10);
    auto EMFrac                     = Monitored::Scalar<float>("EMFrac", -10);
    auto IsoFrac                    = Monitored::Scalar<float>("IsoFrac", -1);
    auto CentFrac                   = Monitored::Scalar<float>("CentFrac", -10);
 
    auto clustersMeanCenterLambda   = Monitored::Scalar<float>("clustersMeanCenterLambda", 0);
    auto clustersMeanFirstEngDens   = Monitored::Scalar<float>("clustersMeanFirstEngDens", 0);
    auto clustersMeanEMProbability  = Monitored::Scalar<float>("clustersMeanEMProbability", 0);
    auto clustersMeanSecondLambda   = Monitored::Scalar<float>("clustersMeanSecondLambda", 0);
    auto clustersMeanPresamplerFrac = Monitored::Scalar<float>("clustersMeanPresamplerFrac", 0);
 
    auto n_clusters                     = Monitored::Scalar<int>("NClusters", 0); 
    std::vector<float> cluster_et_log, cluster_dEta, cluster_dPhi;
    std::vector<float> cluster_log_SECOND_R, cluster_SECOND_LAMBDA, cluster_CENTER_LAMBDA;
    auto mon_cluster_et_log             = Monitored::Collection("cluster_et_log", cluster_et_log);
    auto mon_cluster_dEta               = Monitored::Collection("cluster_dEta", cluster_dEta);
    auto mon_cluster_dPhi               = Monitored::Collection("cluster_dPhi", cluster_dPhi);
    auto mon_cluster_log_SECOND_R       = Monitored::Collection("cluster_log_SECOND_R", cluster_log_SECOND_R);
    auto mon_cluster_SECOND_LAMBDA      = Monitored::Collection("cluster_SECOND_LAMBDA", cluster_SECOND_LAMBDA);
    auto mon_cluster_CENTER_LAMBDA      = Monitored::Collection("cluster_CENTER_LAMBDA", cluster_CENTER_LAMBDA);
    std::vector<unsigned char> calo_errors;
    auto mon_calo_errors                = Monitored::Collection("calo_errors", calo_errors);
 
    // Precision monitored variables
    auto n_tracks                   = Monitored::Scalar<int>("NTracks", -10);
    auto n_iso_tracks               = Monitored::Scalar<int>("NIsoTracks", -10);
 
    auto ipSigLeadTrk               = Monitored::Scalar<float>("ipSigLeadTrk", -1000);
    auto trFlightPathSig            = Monitored::Scalar<float>("trFlightPathSig", -10);
    auto massTrkSys                 = Monitored::Scalar<float>("massTrkSys", -10);
    auto dRmax                      = Monitored::Scalar<float>("dRmax", -10);
    auto trkAvgDist                 = Monitored::Scalar<float>("TrkAvgDist", -1);
    auto innerTrkAvgDist            = Monitored::Scalar<float>("innerTrkAvgDist", -1);
    auto etovPtLead                 = Monitored::Scalar<float>("EtovPtLead", -10);
    auto PSSFraction                = Monitored::Scalar<float>("PSSFraction", -999.9);
    auto EMPOverTrkSysP             = Monitored::Scalar<float>("EMPOverTrkSysP", -999.9);
    auto ChPiEMEOverCaloEME         = Monitored::Scalar<float>("ChPiEMEOverCaloEME", -999.9);
    auto vertex_x                   = Monitored::Scalar<float>("vertex_x", -999.9);
    auto vertex_y                   = Monitored::Scalar<float>("vertex_y", -999.9); 
    auto vertex_z                   = Monitored::Scalar<float>("vertex_z", -999.9);
 
    auto n_all_tracks                           = Monitored::Scalar<int>("NAllTracks", 0);
    std::vector<float> track_pt_log, track_dEta, track_dPhi;
    std::vector<float> track_d0_abs_log, track_z0sinthetaTJVA_abs_log;
    std::vector<float> track_nPixelHitsPlusDeadSensors, track_nSCTHitsPlusDeadSensors;
    auto mon_track_pt_log                       = Monitored::Collection("track_pt_log", track_pt_log);
    auto mon_track_dEta                         = Monitored::Collection("track_dEta", track_dEta);
    auto mon_track_dPhi                         = Monitored::Collection("track_dPhi", track_dPhi);
    auto mon_track_d0_abs_log                   = Monitored::Collection("track_d0_abs_log", track_d0_abs_log);
    auto mon_track_z0sinthetaTJVA_abs_log       = Monitored::Collection("track_z0sinthetaTJVA_abs_log", track_z0sinthetaTJVA_abs_log); 
    auto mon_track_nPixelHitsPlusDeadSensors    = Monitored::Collection("track_nPixelHitsPlusDeadSensors", track_nPixelHitsPlusDeadSensors);
    auto mon_track_nSCTHitsPlusDeadSensors      = Monitored::Collection("track_nSCTHitsPlusDeadSensors", track_nSCTHitsPlusDeadSensors);
    std::vector<unsigned char> track_errors;
    auto mon_track_errors                       = Monitored::Collection("track_errors", track_errors);
 
    std::map<std::string, Monitored::Scalar<float>> monitoredIdVariables;
    for(const auto& [key, p] : m_monitoredIdScores) {
        monitoredIdVariables.emplace(key + "_TauJetScore_0p", Monitored::Scalar<float>(key + "_TauJetScore_0p", -1));
        monitoredIdVariables.emplace(key + "_TauJetScoreTrans_0p", Monitored::Scalar<float>(key + "_TauJetScoreTrans_0p", -1));
        monitoredIdVariables.emplace(key + "_TauJetScore_1p", Monitored::Scalar<float>(key + "_TauJetScore_1p", -1));
        monitoredIdVariables.emplace(key + "_TauJetScoreTrans_1p", Monitored::Scalar<float>(key + "_TauJetScoreTrans_1p", -1));
        monitoredIdVariables.emplace(key + "_TauJetScore_mp", Monitored::Scalar<float>(key + "_TauJetScore_mp", -1));
        monitoredIdVariables.emplace(key + "_TauJetScoreTrans_mp", Monitored::Scalar<float>(key + "_TauJetScoreTrans_mp", -1));
    }
 
    // CaloHits monitored variables
    auto n_hits             = Monitored::Scalar<int>("NHits", 0);
    std::vector<float> hit_z, hit_dPhi;
    auto mon_hit_z          = Monitored::Collection("hit_z", hit_z);
    auto mon_hit_dPhi       = Monitored::Collection("hit_dPhi", hit_dPhi);
 
    std::map<std::string, Monitored::Scalar<float>> monitoredHitZVariables;
    for(const auto& [key, p] : m_monitoredHitZRegressions) {
        monitoredHitZVariables.emplace(key + "_z0", Monitored::Scalar<float>(key + "_z0", -999));
        monitoredHitZVariables.emplace(key + "_z0_sigma", Monitored::Scalar<float>(key + "_z0_sigma", -999));
    }
 
    std::vector<std::reference_wrapper<Monitored::IMonitoredVariable>> monVars = {
        std::ref(n_taus),
        std::ref(pT), std::ref(eta), std::ref(phi),
        std::ref(etaRoI), std::ref(phiRoI), std::ref(dEta_RoI), std::ref(dPhi_RoI),
        std::ref(mEflowApprox), std::ref(ptRatioEflowApprox), std::ref(pt_jetseed_log),
        std::ref(etaDetectorAxis), std::ref(ptDetectorAxis), std::ref(ptDetectorAxis_log),
        std::ref(n_cells),
        std::ref(EMRadius), std::ref(HadRadius), std::ref(EtHad), std::ref(EtEm), std::ref(EMFrac), std::ref(IsoFrac), std::ref(CentFrac),
        std::ref(clustersMeanCenterLambda), std::ref(clustersMeanFirstEngDens), std::ref(clustersMeanEMProbability),
        std::ref(clustersMeanSecondLambda), std::ref(clustersMeanPresamplerFrac),
        std::ref(n_clusters), std::ref(mon_cluster_et_log), std::ref(mon_cluster_dEta), std::ref(mon_cluster_dPhi),
        std::ref(mon_cluster_log_SECOND_R), std::ref(mon_cluster_SECOND_LAMBDA), std::ref(mon_cluster_CENTER_LAMBDA),
        std::ref(mon_calo_errors),
        std::ref(n_tracks), std::ref(n_iso_tracks),
        std::ref(ipSigLeadTrk), std::ref(trFlightPathSig), std::ref(massTrkSys), std::ref(dRmax), std::ref(trkAvgDist), std::ref(innerTrkAvgDist),
        std::ref(etovPtLead), std::ref(PSSFraction), std::ref(EMPOverTrkSysP), std::ref(ChPiEMEOverCaloEME),
        std::ref(vertex_x), std::ref(vertex_y), std::ref(vertex_z),
        std::ref(n_all_tracks), std::ref(mon_track_pt_log), std::ref(mon_track_dEta), std::ref(mon_track_dPhi),
        std::ref(mon_track_d0_abs_log), std::ref(mon_track_z0sinthetaTJVA_abs_log),
        std::ref(mon_track_nPixelHitsPlusDeadSensors), std::ref(mon_track_nSCTHitsPlusDeadSensors),
        std::ref(mon_track_errors),
        std::ref(n_hits), std::ref(mon_hit_z), std::ref(mon_hit_dPhi)
    };
    for(auto& [key, var] : monitoredIdVariables) monVars.push_back(std::ref(var));
    for(auto& [key, var] : monitoredHitZVariables) monVars.push_back(std::ref(var));
    auto monitorIt = Monitored::Group(m_monTool, monVars);
 
 
    ATH_MSG_DEBUG("Executing TrigTauRecMerged");
 
    //===============================================================================
    // Main TauJet object and collections
    //===============================================================================
    xAOD::TauJet* tau = nullptr;
 
    // We store the position of hits assigned to the tau for later monitoring
    std::vector<ROOT::Math::XYZVector> tauHits;
 
 
    //===============================================================================
    // Retrieve RoI
    //===============================================================================
    SG::ReadHandle<TrigRoiDescriptorCollection> roisHandle = SG::makeHandle(m_roiInputKey, ctx);
    if(!roisHandle.isValid()) {
        ATH_MSG_ERROR("No RoIHandle found");
        return StatusCode::FAILURE;
    }
 
    if(roisHandle->empty()) {
        ATH_MSG_ERROR("Empty RoIHandle");
        return StatusCode::FAILURE;
    }
 
    const TrigRoiDescriptor *roiDescriptor = roisHandle->at(0);
    if(roiDescriptor) {
        ATH_MSG_DEBUG("RoI: " << *roiDescriptor);
    } else {
        ATH_MSG_ERROR("Failed to find TrigRoiDescriptor");
        calo_errors.push_back(NoROIDescr);
        return StatusCode::FAILURE;
    }
 
 
    //===============================================================================
    // Prepare Output TauJet and TauTrack containers
    //===============================================================================
    
    // Create and register the output TauJetContainer
    std::unique_ptr<xAOD::TauJetContainer> outputContainer = std::make_unique<xAOD::TauJetContainer>();
    std::unique_ptr<xAOD::TauJetAuxContainer> outputAuxContainer = std::make_unique<xAOD::TauJetAuxContainer>();
    outputContainer->setStore(outputAuxContainer.get());
 
    SG::WriteHandle<xAOD::TauJetContainer> outputTauHandle(m_tauJetOutputKey, ctx);
    ATH_CHECK(outputTauHandle.record(std::move(outputContainer), std::move(outputAuxContainer)));
 
    // Create and register the output TauTrackContainer
    std::unique_ptr<xAOD::TauTrackContainer> outputTrackContainer = std::make_unique<xAOD::TauTrackContainer>();
    std::unique_ptr<xAOD::TauTrackAuxContainer> outputTrackAuxContainer = std::make_unique<xAOD::TauTrackAuxContainer>();
    outputTrackContainer->setStore(outputTrackAuxContainer.get());
 
    SG::WriteHandle<xAOD::TauTrackContainer> tauTrackHandle(m_tauTrackOutputKey, ctx);
    ATH_CHECK(tauTrackHandle.record(std::move(outputTrackContainer), std::move(outputTrackAuxContainer)));
    
 
    //===============================================================================
    // Initial TauJet calo-reco / input-retrieval
    //===============================================================================
    // We now have two options for the TauJet reconstruction:
    //  1) Reconstruct the TauJet from scratch, using the bare calo-clusters (1st trigger step)
    //  2) Fetch a preceding TauJet (and [dummy] TauTracks), and use them to seed the current reconstruction
    
    //-------------------------------------------------------------------------------
    // Option 1: Calorimeter-only reconstruction from clusters (CaloMVA step)
    //-------------------------------------------------------------------------------
    if(m_reco_mode == Mode::FromClusters) {
        // Retrieve Calocluster container
        SG::ReadHandle<xAOD::CaloClusterContainer> CCContainerHandle = SG::makeHandle(m_clustersInputKey, ctx);
        ATH_CHECK(CCContainerHandle.isValid());
 
        const xAOD::CaloClusterContainer *RoICaloClusterContainer = CCContainerHandle.get();
        if(RoICaloClusterContainer) {
            ATH_MSG_DEBUG("CaloCluster container found of size: " << RoICaloClusterContainer->size());
            
            // If size is zero, don't stop just continue to produce empty TauJetCollection
            if(RoICaloClusterContainer->empty()) calo_errors.push_back(NoClustCont);
        } else {
            ATH_MSG_ERROR("No CaloCluster container found");
            calo_errors.push_back(NoClustCont);
            return StatusCode::FAILURE;
        }
 
        // Also create the seed-jet containers
        std::unique_ptr<xAOD::JetContainer> jetContainer{std::make_unique<xAOD::JetContainer>()};
        std::unique_ptr<xAOD::JetAuxContainer> jetAuxContainer{std::make_unique<xAOD::JetAuxContainer>()};
        jetContainer->setStore(jetAuxContainer.get());
 
        SG::WriteHandle<xAOD::JetContainer> outputTauSeedJetHandle = SG::makeHandle(m_tauSeedOutputKey, ctx);
        ATH_CHECK(outputTauSeedJetHandle.record(std::move(jetContainer), std::move(jetAuxContainer)));
 
        // And create the seed-jet object itself
        outputTauSeedJetHandle->push_back(std::make_unique<xAOD::Jet>());
        xAOD::Jet *jet = outputTauSeedJetHandle->back();
        
 
        // Build the jet, also keep track of the kinematics by hand
        // Eventually, want to use FastJet here?
        // We are using calibrated clusters, we need to keep track of this
        jet->setConstituentsSignalState(xAOD::JetConstitScale::CalibratedJetConstituent);
        TLorentzVector cluster_p4, barycenter;
        for(const xAOD::CaloCluster* cluster : *RoICaloClusterContainer) {
            ATH_MSG_DEBUG("Cluster (e, eta, phi): (" << cluster->e() << ", " << cluster->eta() << ", " << cluster->phi() << ")");
        
            if(cluster->e() < 0) {
                ATH_MSG_DEBUG("Negative energy cluster is rejected");
                continue;
            }
 
            cluster_p4.SetPtEtaPhiE(cluster->pt(), cluster->eta(), cluster->phi(), cluster->e());
            jet->addConstituent(cluster);
 
            barycenter += cluster_p4;
        }
         
        jet->setJetP4(xAOD::JetFourMom_t(barycenter.Pt(), barycenter.Eta(), barycenter.Phi(), barycenter.M())); 
        ATH_MSG_DEBUG("Built jet with eta: " << jet->eta() << ", phi: " << jet->phi() << ", pT: " << jet->pt() << ", E: "<< jet->e() );
 
        
        // If we're running calo-clustering, that means we just started the HLT reco,
        // and there's no input TauJet container to this step. Create one instead!
        outputTauHandle->push_back(std::make_unique<xAOD::TauJet>());
        tau = outputTauHandle->back();
        tau->setROIWord(roiDescriptor->roiWord());
 
        // Using the new Jet collection, setup the tau candidate structure
        tau->setJet(outputTauSeedJetHandle.ptr(), jet);
 
        // Fix eta, phi in case the jet's energy is negative
        if(jet->e() <= 0) {
            ATH_MSG_DEBUG("Changing (eta, phi) back to the RoI center due to negative energy: " << jet->e());
            tau->setP4(tau->pt(), roiDescriptor->eta(), roiDescriptor->phi(), tau->m());      
            ATH_MSG_DEBUG("Roi: " << roiDescriptor->roiId() << ", tau eta: " << tau->eta() << ", tau phi: " << tau->phi() );
        }
    }
 
 
    //-------------------------------------------------------------------------------
    // Option 2: Use input TauJet (and TauTracks) as seeds (non calo-only reco)
    //-------------------------------------------------------------------------------
    else if(m_reco_mode == Mode::FromTauJet) {
        // Retrieve input TauJet container
        SG::ReadHandle<xAOD::TauJetContainer> tauInputHandle(m_tauJetInputKey, ctx);
        const xAOD::TauJetContainer* inputTauContainer = tauInputHandle.cptr();
        ATH_MSG_DEBUG("Input TauJet Container size: " << inputTauContainer->size());
 
        // Copy the input TauJets to the output container
        ATH_CHECK(deepCopy(outputTauHandle, inputTauContainer));
 
        // Copy associated hits
        // Store associated hits for later monitoring
        if(!m_hitsInputDecorKey.empty()) {
            // Initialize hits decorators
            SG::ReadDecorHandle<xAOD::TauJetContainer, std::vector<ElementLink<xAOD::TrackMeasurementValidationContainer>>> hitsInputDecorHandle(m_hitsInputDecorKey, ctx);
            ATH_CHECK(hitsInputDecorHandle.isValid());
            SG::WriteDecorHandle<xAOD::TauJetContainer, std::vector<ElementLink<xAOD::TrackMeasurementValidationContainer>>> hitsOutputDecorHandle(m_hitsOutputDecorKey, ctx);
            ATH_CHECK(hitsOutputDecorHandle.isValid());
 
            // We use ConstAccessors instead of ReadDecorHandles to read the hit positions
            // because hits might come from different original containers if we're including
            // both Pixel and Strips. It is ok, because the pre-processing pipeline required to
            // attach the hit ELs to the TauJet includes calls to HitDecoratorAlg, so the
            // decorations will always be there, no matter how the scheduler organizes the execution.
            static const SG::AuxElement::ConstAccessor<float> x("HitsXRelToBeamspot");
            static const SG::AuxElement::ConstAccessor<float> y("HitsYRelToBeamspot");
            static const SG::AuxElement::ConstAccessor<float> z("HitsZRelToBeamspot");
 
            for(size_t i = 0; i < inputTauContainer->size(); ++i) {
                const xAOD::TauJet* inputTau = inputTauContainer->at(i);
                const auto& inputHits = hitsInputDecorHandle(*inputTau);
 
                for(const ElementLink<xAOD::TrackMeasurementValidationContainer>& el : inputHits) {
                    if(el.isValid()) {
                        tauHits.push_back(ROOT::Math::XYZVector(x(**el), y(**el), z(**el)));
                    }
                }
 
                const xAOD::TauJet* outputTau = outputTauHandle->at(i);
                hitsOutputDecorHandle(*outputTau) = inputHits;
            }
        }
 
        
        // Copy input TauJet decorations to output TauJets
        for(size_t i = 0; i < m_tauJetInputDecorKeysArray.size(); ++i) {
            SG::ReadDecorHandle<xAOD::TauJetContainer, float> inputDecorHandle(m_tauJetInputDecorKeysArray[i], ctx);
            ATH_CHECK(inputDecorHandle.isValid());
            SG::WriteDecorHandle<xAOD::TauJetContainer, float> outputDecorHandle(m_tauJetOutputDecorKeysArray[i], ctx);
            ATH_CHECK(outputDecorHandle.isValid());
 
            for(size_t j = 0; j < inputTauContainer->size(); ++j) {
                outputDecorHandle(*outputTauHandle->at(j)) = inputDecorHandle(*inputTauContainer->at(j));
            }
        }
 
        // Retrieve input TauTrack container
        if(!m_tauTrackInputKey.key().empty()) {
            SG::ReadHandle<xAOD::TauTrackContainer> tauTrackInputHandle(m_tauTrackInputKey, ctx);
            const xAOD::TauTrackContainer* inputTauTrackContainer = tauTrackInputHandle.cptr();
            ATH_MSG_DEBUG("Tau Track Container Size " << inputTauTrackContainer->size());
 
            // Copy the input TauTracks to the output container
            ATH_CHECK(deepCopy(tauTrackHandle, inputTauTrackContainer));
        }
 
 
        // Now retrieve the main TauJet object, if available (the recently created copy of the input TauJet)
        if(!outputTauHandle->empty()) {
            tau = outputTauHandle->back();
            
            // Check if the tau has a valid jetLink
            ATH_CHECK(tau->jetLink().isValid());
 
            // Clear all previous TauTrack links (if any), since we will run the
            // TauTrackFinder tool (and the associated InDet helpers) in the
            // Presel and Precision steps, refilling the tracks
            if(!m_tauTrackInputKey.key().empty()) tau->clearTauTrackLinks();
        }
    }
   
 
    //===============================================================================
    // Get Vertex Container (optional)
    //===============================================================================
    const xAOD::VertexContainer* RoIVxContainer = nullptr;
    if(!m_vertexInputKey.key().empty()){
        SG::ReadHandle<xAOD::VertexContainer> VertexContainerHandle = SG::makeHandle(m_vertexInputKey, ctx);
 
        if(!VertexContainerHandle.isValid()) {
            ATH_MSG_DEBUG("No VertexContainer retrieved for the trigger element");
            track_errors.push_back(NoVtxCont);
        } else {
            RoIVxContainer = VertexContainerHandle.get();
            ATH_MSG_DEBUG("Size of VertexContainer: " << RoIVxContainer->size());
        }
    }
 
 
    ATH_MSG_DEBUG("roiDescriptor roiWord: " << roiDescriptor->roiWord() << ", saved in TauJet: " << tau->ROIWord());
 
 
    //===============================================================================
    // Loop over all booked tau tools:
    //===============================================================================
    StatusCode processStatus = StatusCode::SUCCESS;
 
    // Sequence: VertexFinderTools -> CommonToolsBeforeTF -> TrackFinderTools -> CommonTools -> VertexVarsTools -> IDTools
  
    for(const auto& tool : m_vertexFinderTools) {
        ATH_MSG_DEBUG("Starting Tool: " << tool->name());
 
        processStatus = tool->executeVertexFinder(*tau, RoIVxContainer);
 
        if(!processStatus.isFailure()) {
            ATH_MSG_DEBUG(" " << tool->name() << " executed successfully");
        } else {
            ATH_MSG_DEBUG(" " << tool->name() << " execution failed");
            break;
        }
    }
 
    for(const auto& tool : m_commonToolsBeforeTF) {
        if(!processStatus.isFailure()) ATH_MSG_DEBUG("Starting Tool: " << tool->name());
        else break;
 
        processStatus = tool->execute(*tau);
 
        if(!processStatus.isFailure()) {
            ATH_MSG_DEBUG(" " << tool->name() << " executed successfully");
        } else {
            ATH_MSG_DEBUG(" " << tool->name() << " execution failed");
            break;
        }
    }
 
    for(const auto& tool : m_trackFinderTools) {
        if(!processStatus.isFailure()) ATH_MSG_DEBUG("Starting Tool: " << tool->name());
        else break;
 
        processStatus = tool->executeTrackFinder(*tau, *tauTrackHandle);
 
        if(!processStatus.isFailure()) {
            ATH_MSG_DEBUG(" " << tool->name() << " executed successfully");
        } else {
            ATH_MSG_DEBUG(" " << tool->name() << " execution failed");
            break;
        }
    }
 
    for(const auto& tool : m_commonTools) {
        if(!processStatus.isFailure()) ATH_MSG_DEBUG("Starting Tool: " << tool->name());
        else break;
 
        processStatus = tool->execute(*tau);
 
        if(!processStatus.isFailure()) {
            ATH_MSG_DEBUG(" " << tool->name() << " executed successfully");
        } else {
            ATH_MSG_DEBUG(" " << tool->name() << " execution failed");
            break;
        }
    }
 
    // Dummy container passed to TauVertexVariables, not used in trigger though
    xAOD::VertexContainer dummyVxCont; 
    for(const auto& tool : m_vertexVarsTools) {
        if(!processStatus.isFailure()) ATH_MSG_DEBUG("Starting Tool: " << tool->name());
        else break;
 
        processStatus = tool->executeVertexVariables(*tau, dummyVxCont);
 
        if(!processStatus.isFailure()) {
            ATH_MSG_DEBUG(" " << tool->name() << " executed successfully");
        } else {
            ATH_MSG_DEBUG(" " << tool->name() << " execution failed");
            break;
        }
    }
 
    for(const auto& tool : m_idTools) {
        if(!processStatus.isFailure()) ATH_MSG_DEBUG("Starting Tool: " << tool->name());
        else break;
 
        processStatus = tool->execute(*tau);
 
        if(!processStatus.isFailure()) {
            ATH_MSG_DEBUG(" " << tool->name() << " executed successfully");
        } else {
            ATH_MSG_DEBUG(" " << tool->name() << " execution failed");
            break;
        }
    }
 
    ATH_MSG_DEBUG("This tau has " << tau->allTracks() << " tracks linked");
 
 
    // Cleanup in case any of the tools failed (rejected Tau)
    if(!processStatus.isSuccess()) {
        ATH_MSG_DEBUG("The tau object has NOT been registered in the tau container");
 
        xAOD::TauJet* bad_tau = outputTauHandle->back();
        ATH_MSG_DEBUG("Deleting " << bad_tau->nAllTracks() << " tracks associated with tau");
        tauTrackHandle->erase(tauTrackHandle->end() - bad_tau->nAllTracks(), tauTrackHandle->end());
 
        outputTauHandle->pop_back();
 
        ATH_MSG_DEBUG("Clean up done after jet seed");  
 
    } else {
        // Check that the seed-jet energy is positive
        // Otherwise, try to salvage it by replacing the tau position with the RoI's
        float fJetEnergy = (*tau->jetLink())->e();
        ATH_MSG_DEBUG("Seed jet E: " << fJetEnergy);
          
        if(fJetEnergy < 0.00001) {
            ATH_MSG_DEBUG("Changing tau's (eta,phi) to RoI ones due to negative energy (PxPyPzE flips eta and phi)");
            ATH_MSG_DEBUG("This is probably not needed anymore, method PxPyPzE has been corrected");
            // TODO: Do we still need this??
            
            tau->setP4(tau->pt(), roiDescriptor->eta(), roiDescriptor->phi(), tau->m());
            
            ATH_MSG_DEBUG("Roi: " << roiDescriptor->roiId() << ", Tau eta: " << tau->eta() << ", phi: " << tau->phi() << ", pT: " << tau->pt());
        }
 
 
        //===============================================================================
        // Monitor tau variables
        //===============================================================================
        
        pT = tau->pt()/Gaudi::Units::GeV;
        eta = tau->eta();
        phi = tau->phi();
 
 
        etaRoI = roiDescriptor->eta();
        phiRoI = roiDescriptor->phi();     
        dEta_RoI = eta - roiDescriptor->eta();
        dPhi_RoI = phi - roiDescriptor->phi();
        if(dPhi_RoI < -M_PI) dPhi_RoI += 2.0*M_PI;
        if(dPhi_RoI > M_PI) dPhi_RoI -= 2.0*M_PI;
 
 
        float pre_mEflowApprox;
        tau->detail(xAOD::TauJetParameters::mEflowApprox, pre_mEflowApprox);  
        mEflowApprox = std::log10(std::max(pre_mEflowApprox, 140.0f));
 
        float pre_ptRatioEflowApprox;
        tau->detail(xAOD::TauJetParameters::ptRatioEflowApprox, pre_ptRatioEflowApprox);
        ptRatioEflowApprox = std::min(pre_ptRatioEflowApprox, 4.0f);
        
        pt_jetseed_log = std::log10(tau->ptJetSeed());
        etaDetectorAxis = tau->etaDetectorAxis();   
        ptDetectorAxis = std::min(tau->ptDetectorAxis()/Gaudi::Units::GeV, 10000.0);
        ptDetectorAxis_log = std::log10(std::min(tau->ptDetectorAxis()/Gaudi::Units::GeV, 10000.0));
 
 
        tau->detail(xAOD::TauJetParameters::numCells, n_cells);
        tau->detail(xAOD::TauJetParameters::EMRadius, EMRadius);
        tau->detail(xAOD::TauJetParameters::hadRadius, HadRadius);
        tau->detail(xAOD::TauJetParameters::etHadAtEMScale, EtHad);
        EtHad /= Gaudi::Units::GeV;
        tau->detail(xAOD::TauJetParameters::etEMAtEMScale, EtEm);
        EtEm /= Gaudi::Units::GeV;
 
        float Et_raw = EtEm + EtHad;
        if(Et_raw != 0) EMFrac = EtEm / Et_raw;
 
        tau->detail(xAOD::TauJetParameters::isolFrac, IsoFrac);
        tau->detail(xAOD::TauJetParameters::centFrac, CentFrac);
 
 
        // Monitor BRT variables
        float tmp = 0;
        bool test = tau->detail(xAOD::TauJetParameters::ClustersMeanCenterLambda, tmp);
        if(test) clustersMeanCenterLambda = tmp;
        test = tau->detail(xAOD::TauJetParameters::ClustersMeanFirstEngDens, tmp);
        if(test) clustersMeanFirstEngDens = tmp;
        test = tau->detail(xAOD::TauJetParameters::ClustersMeanEMProbability, tmp);
        if(test) clustersMeanEMProbability = tmp;
        test = tau->detail(xAOD::TauJetParameters::ClustersMeanSecondLambda, tmp);
        if(test) clustersMeanSecondLambda = tmp;
        test = tau->detail(xAOD::TauJetParameters::ClustersMeanPresamplerFrac, tmp);
        if(test) clustersMeanPresamplerFrac = tmp;
 
        
        // Cluster variables monitoring
        n_clusters = tau->clusters().size();
        for(const auto& cluster : tau->clusters()) {
            const xAOD::CaloCluster* cls = dynamic_cast<const xAOD::CaloCluster*>(cluster);
 
            cluster_et_log.push_back(std::log10(cls->et()));
            cluster_dEta.push_back(cls->eta() - tau->eta());
            cluster_dPhi.push_back(cls->p4().DeltaPhi(tau->p4()));
         
            double log_second_R = -999;
            const bool success_SECOND_R = cls->retrieveMoment(xAOD::CaloCluster::MomentType::SECOND_R, log_second_R);
            if(success_SECOND_R) log_second_R = std::log10(log_second_R + 0.1);
            cluster_log_SECOND_R.push_back(log_second_R);
 
            double second_lambda = -999;
            const bool success_SECOND_LAMBDA = cls->retrieveMoment(xAOD::CaloCluster::MomentType::SECOND_LAMBDA, second_lambda);
            if(success_SECOND_LAMBDA) second_lambda = std::log10(second_lambda + 0.1);
            cluster_SECOND_LAMBDA.push_back(second_lambda);
 
            double center_lambda = -999;
            const bool success_CENTER_LAMBDA = cls->retrieveMoment(xAOD::CaloCluster::MomentType::CENTER_LAMBDA, center_lambda);
            if(success_CENTER_LAMBDA) center_lambda = std::log10(center_lambda + 1e-6);
            cluster_CENTER_LAMBDA.push_back(center_lambda);     
        }
 
 
        // Tracks summary monitoring
        n_tracks = tau->nTracks();
        n_iso_tracks = tau->nTracksIsolation();
        tau->detail(xAOD::TauJetParameters::ipSigLeadTrk, ipSigLeadTrk);
        if(tau->nTracks() > 0) ipSigLeadTrk = std::abs(tau->track(0)->d0SigTJVA()); // TODO: Is this needed?
        tau->detail(xAOD::TauJetParameters::trFlightPathSig, trFlightPathSig);
        tau->detail(xAOD::TauJetParameters::massTrkSys, massTrkSys);
        massTrkSys /= Gaudi::Units::GeV;
        tau->detail(xAOD::TauJetParameters::dRmax, dRmax);
        tau->detail(xAOD::TauJetParameters::trkAvgDist, trkAvgDist);
        tau->detail(xAOD::TauJetParameters::innerTrkAvgDist, innerTrkAvgDist);  
        tau->detail(xAOD::TauJetParameters::etOverPtLeadTrk, etovPtLead);
        tau->detail(xAOD::TauJetParameters::PSSFraction, PSSFraction);
        tau->detail(xAOD::TauJetParameters::EMPOverTrkSysP, EMPOverTrkSysP);
        tau->detail(xAOD::TauJetParameters::ChPiEMEOverCaloEME, ChPiEMEOverCaloEME);
 
        if(tau->vertexLink().isValid() && tau->vertex() && tau->vertex()->vertexType() != xAOD::VxType::NoVtx) {
            vertex_x = tau->vertex()->x();
            vertex_y = tau->vertex()->y();       
            vertex_z = tau->vertex()->z();
        }
 
 
        // Track variables monitoring 
        n_all_tracks = tau->allTracks().size();
        for(const xAOD::TauTrack* track : tau->allTracks()) {
            track_pt_log.push_back(std::log10(track->pt()));
            track_dEta.push_back(track->eta() - tau->eta()); 
            track_dPhi.push_back(track->p4().DeltaPhi(tau->p4()));
            track_d0_abs_log.push_back(std::log10(std::abs(track->track()->d0()) + 1e-6));
            track_z0sinthetaTJVA_abs_log.push_back(track->z0sinthetaTJVA());
 
            uint8_t pixel_hits, pixel_dead;
            const bool success1_pixel_hits = track->track()->summaryValue(pixel_hits, xAOD::numberOfPixelHits);
            const bool success2_pixel_dead = track->track()->summaryValue(pixel_dead, xAOD::numberOfPixelDeadSensors);
            float nPixelHitsPlusDeadSensor = -999;
            if(success1_pixel_hits && success2_pixel_dead) nPixelHitsPlusDeadSensor = pixel_hits + pixel_dead;
            track_nPixelHitsPlusDeadSensors.push_back(nPixelHitsPlusDeadSensor);
 
            uint8_t sct_hits, sct_dead;
            const bool success1_sct_hits = track->track()->summaryValue(sct_hits, xAOD::numberOfSCTHits);
            const bool success2_sct_dead = track->track()->summaryValue(sct_dead, xAOD::numberOfSCTDeadSensors);
            float nSCTHitsPlusDeadSensors = -999;
            if(success1_sct_hits && success2_sct_dead) nSCTHitsPlusDeadSensors = sct_hits + sct_dead;
            track_nSCTHitsPlusDeadSensors.push_back(nSCTHitsPlusDeadSensors);
        }
 
 
        // Hits summary monitoring
        if(!m_hitsInputDecorKey.empty()) {
            n_hits = tauHits.size();
            for(const ROOT::Math::XYZVector& hit : tauHits) {
                hit_z.push_back(hit.z());
                hit_dPhi.push_back(CxxUtils::wrapToPi(hit.phi() - tau->phi()));
            }
        }
 
 
        // TauID Score monitoring
        for(const auto& [key, p] : m_monitoredIdScores) {
            const auto& [score, scoreTrans] = m_monitoredInferenceAccessors.at(key);
 
            if(score.isAvailable(*tau)) {
                if(tau->nTracks() == 0) {
                    monitoredIdVariables.at(key + "_TauJetScore_0p") = score(*tau);
                } else if(tau->nTracks() == 1) {
                    monitoredIdVariables.at(key + "_TauJetScore_1p") = score(*tau);
                } else { // MP tau
                    monitoredIdVariables.at(key + "_TauJetScore_mp") = score(*tau);
                }
            }
 
            if(scoreTrans.isAvailable(*tau)) {
                if(tau->nTracks() == 0) {
                    monitoredIdVariables.at(key + "_TauJetScoreTrans_0p") = scoreTrans(*tau);
                } else if(tau->nTracks() == 1) {
                    monitoredIdVariables.at(key + "_TauJetScoreTrans_1p") = scoreTrans(*tau);
                } else { // MP tau
                    monitoredIdVariables.at(key + "_TauJetScoreTrans_mp") = scoreTrans(*tau);
                }
            }
        }
 
        
        // HitZ inference monitoring
        for(const auto& [key, p] : m_monitoredHitZRegressions) {
            const auto& [z0, sigma] = m_monitoredInferenceAccessors.at(key);
 
            if(z0.isAvailable(*tau)) {
                monitoredHitZVariables.at(key + "_z0") = z0(*tau);
            }
            if(sigma.isAvailable(*tau)) {
                monitoredHitZVariables.at(key + "_z0_sigma") = sigma(*tau);
            }
        }
 
 
        ++n_taus;
 
 
        ATH_MSG_DEBUG("RoI: " << roiDescriptor->roiId()
              << ", Tau pT (GeV): " << pT << ", Tau eta: " << eta << ", Tau phi: " << phi
              << ", wrt RoI dEta: " << dEta_RoI << ", dPhi: " << dPhi_RoI);
    }
 
    //-------------------------------------------------------------------------------
    // All done!
    //-------------------------------------------------------------------------------
    
    ATH_MSG_DEBUG("Output TauJetContainer size: " << outputTauHandle->size());
    ATH_MSG_DEBUG("Output TauJetTrackContainer size: " << tauTrackHandle->size());
    
    return StatusCode::SUCCESS;
}

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

extraOutputDeps()

const DataObjIDColl & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

initialize()

StatusCode TrigTauRecMerged::initialize ( )

overridevirtual

Definition at line 47 of file TrigTauRecMerged.cxx.

{
    ATH_MSG_DEBUG("Initialize");
 
    for(const auto& tool : m_commonTools) ATH_CHECK(tool.retrieve());
    for(const auto& tool : m_commonToolsBeforeTF) ATH_CHECK(tool.retrieve());
    for(const auto& tool : m_vertexFinderTools) ATH_CHECK(tool.retrieve());
    for(const auto& tool : m_trackFinderTools) ATH_CHECK(tool.retrieve());
    for(const auto& tool : m_vertexVarsTools) ATH_CHECK(tool.retrieve());
    for(const auto& tool : m_idTools) ATH_CHECK(tool.retrieve());
 
    if(!m_monTool.name().empty()) ATH_CHECK(m_monTool.retrieve());
  
    ATH_MSG_DEBUG("Initialising handle keys");
    ATH_CHECK(m_roiInputKey.initialize());
    ATH_CHECK(m_clustersInputKey.initialize(SG::AllowEmpty));
    ATH_CHECK(m_vertexInputKey.initialize(SG::AllowEmpty));
    ATH_CHECK(m_tauJetInputKey.initialize(SG::AllowEmpty));
    ATH_CHECK(m_tauTrackInputKey.initialize(SG::AllowEmpty));
 
    ATH_CHECK(m_tauSeedOutputKey.initialize(SG::AllowEmpty));
    ATH_CHECK(m_tauJetOutputKey.initialize());
    ATH_CHECK(m_tauTrackOutputKey.initialize());
 
    // Hits decorations
    if(!m_tauJetInputKey.empty() && !m_hitsInputDecorKey.empty()) {
        // Initialize output decoration keys with the same names
        // We manually copy the decorations because we want to make the scheduler aware of the hits dependencies
        m_hitsOutputDecorKey = m_tauJetOutputKey.key() + "."
            + (m_hitsOutputDecorKey.empty() ? m_hitsInputDecorKey.key() : m_hitsOutputDecorKey.key());
        ATH_CHECK(m_hitsOutputDecorKey.initialize());
 
        // Initialize input decorator keys
        m_hitsInputDecorKey = m_tauJetInputKey.key() + "." + m_hitsInputDecorKey.key();
        ATH_CHECK(m_hitsInputDecorKey.initialize());
    }
 
 
    // Propagate decorated quantities from the input to the output tau collection
    if(!m_tauJetInputKey.empty() && !m_tauJetInputDecorKeysArray.empty()) {
        for(auto& decor_key : m_tauJetInputDecorKeysArray) {
            m_tauJetOutputDecorKeysArray.push_back(m_tauJetOutputKey.key() + "." + decor_key.key());
            ATH_CHECK(m_tauJetOutputDecorKeysArray.back().initialize());
 
            decor_key = m_tauJetInputKey.key() + "." + decor_key.key();
            ATH_CHECK(decor_key.initialize());
        }
 
    } else if(m_tauJetInputKey.empty() && !m_tauJetInputDecorKeysArray.empty()) {
        ATH_MSG_ERROR("Cannot specify input TauJet decoration keys without an input TauJet container key!");
        return StatusCode::FAILURE;
    }
 
 
    // Determine the reconstruction mode
    if(m_tauJetInputKey.empty() && !m_clustersInputKey.empty()) {
        m_reco_mode = FromClusters;
        ATH_MSG_INFO("Running in Calo reconstruction mode (from clusters)");
    } else if(!m_tauJetInputKey.empty()) {
        m_reco_mode = FromTauJet;
        ATH_MSG_INFO("Running in Full reconstruction mode (from input TauJets + TauTracks)");
    } else {
        ATH_MSG_ERROR("Invalid configuration!");
        return StatusCode::FAILURE;
    }
 
 
    // Set up the monitoring accessors
    for(const auto& [key, p] : m_monitoredIdScores) {
        m_monitoredInferenceAccessors.emplace(
            key,
            std::make_pair(SG::ConstAccessor<float>(p.first), SG::ConstAccessor<float>(p.second))
        );
    }
 
    for(const auto& [key, p] : m_monitoredHitZRegressions) {
        m_monitoredInferenceAccessors.emplace(
            key,
            std::make_pair(SG::ConstAccessor<float>(p.first), SG::ConstAccessor<float>(p.second))
        );
    }
 
 
    return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

isClonable()

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Reimplemented in HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ITkPixelCablingAlg, ITkStripCablingAlg, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, and Simulation::BeamEffectsAlg.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

{
  // Reentrant algorithms are clonable.
  return true;
}

msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state, const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

                                                               {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }      
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;  
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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_clustersInputKey

SG::ReadHandleKey<xAOD::CaloClusterContainer> TrigTauRecMerged::m_clustersInputKey {this, "InputCaloClusterContainer", "", "Caloclusters in view"}

private

Definition at line 86 of file TrigTauRecMerged.h.

{this, "InputCaloClusterContainer", "", "Caloclusters in view"};

m_commonTools

const ToolHandleArray<ITauToolBase> TrigTauRecMerged::m_commonTools {this, "CommonTools", {}, "List of ITauToolBase common tools"}

private

Definition at line 70 of file TrigTauRecMerged.h.

{this, "CommonTools", {}, "List of ITauToolBase common tools"};

m_commonToolsBeforeTF

const ToolHandleArray<ITauToolBase> TrigTauRecMerged::m_commonToolsBeforeTF {this, "CommonToolsBeforeTF", {}, "List of ITauToolBase common tools to execute before the Track Finder tools"}

private

Definition at line 71 of file TrigTauRecMerged.h.

{this, "CommonToolsBeforeTF", {}, "List of ITauToolBase common tools to execute before the Track Finder tools"};

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

m_hitsInputDecorKey

SG::ReadDecorHandleKey<xAOD::TauJetContainer> TrigTauRecMerged::m_hitsInputDecorKey {this, "InputTauJetHitsKey", "", "Input TauJet hits decoration key"}

private

Definition at line 91 of file TrigTauRecMerged.h.

{this, "InputTauJetHitsKey", "", "Input TauJet hits decoration key"};

m_hitsOutputDecorKey

SG::WriteDecorHandleKey<xAOD::TauJetContainer> TrigTauRecMerged::m_hitsOutputDecorKey {this, "OutputTauJetHitsKey", "", "Output TauJet hits decoration key"}

private

Definition at line 98 of file TrigTauRecMerged.h.

{this, "OutputTauJetHitsKey", "", "Output TauJet hits decoration key"};

m_idTools

const ToolHandleArray<ITauToolBase> TrigTauRecMerged::m_idTools {this, "IDTools", {}, "Vertex Variables tools"}

private

Definition at line 75 of file TrigTauRecMerged.h.

{this, "IDTools", {}, "Vertex Variables tools"};

m_monitoredHitZRegressions

Gaudi::Property<std::map<std::string, std::pair<std::string, std::string> > > TrigTauRecMerged::m_monitoredHitZRegressions {this, "MonitoredHitZRegressions", {}, "Pairs of the z and sigma regression output variables for each HitZ algorithm to be monitored"}

private

Definition at line 80 of file TrigTauRecMerged.h.

{this, "MonitoredHitZRegressions", {}, "Pairs of the z and sigma regression output variables for each HitZ algorithm to be monitored"};

m_monitoredIdScores

Gaudi::Property<std::map<std::string, std::pair<std::string, std::string> > > TrigTauRecMerged::m_monitoredIdScores {this, "MonitoredIDScores", {}, "Pairs of the TauID score and signal-transformed scores for each TauID algorithm to be monitored"}

private

Definition at line 79 of file TrigTauRecMerged.h.

{this, "MonitoredIDScores", {}, "Pairs of the TauID score and signal-transformed scores for each TauID algorithm to be monitored"};

m_monitoredInferenceAccessors

std::map<std::string, std::pair<SG::ConstAccessor<float>, SG::ConstAccessor<float> > > TrigTauRecMerged::m_monitoredInferenceAccessors

private

Definition at line 81 of file TrigTauRecMerged.h.

m_monTool

const ToolHandle<GenericMonitoringTool> TrigTauRecMerged::m_monTool {this, "MonTool", "", "Monitoring tool"}

private

Definition at line 78 of file TrigTauRecMerged.h.

{this, "MonTool", "", "Monitoring tool"};

m_reco_mode

Mode TrigTauRecMerged::m_reco_mode = FromClusters

private

Definition at line 101 of file TrigTauRecMerged.h.

m_roiInputKey

SG::ReadHandleKey<TrigRoiDescriptorCollection> TrigTauRecMerged::m_roiInputKey {this, "InputRoIs", "", "Input RoI name"}

private

Definition at line 85 of file TrigTauRecMerged.h.

{this, "InputRoIs", "", "Input RoI name"};

m_tauJetInputDecorKeysArray

SG::ReadDecorHandleKeyArray<xAOD::TauJetContainer> TrigTauRecMerged::m_tauJetInputDecorKeysArray {this, "InputTauJetCopyDecorKeys", {}, "Array of input TauJet decoration keys to copy on the output TauJet"}

private

Definition at line 89 of file TrigTauRecMerged.h.

{this, "InputTauJetCopyDecorKeys", {}, "Array of input TauJet decoration keys to copy on the output TauJet"};

m_tauJetInputKey

SG::ReadHandleKey<xAOD::TauJetContainer> TrigTauRecMerged::m_tauJetInputKey {this, "InputTauJetContainer", "", "Input TauJet container"}

private

Definition at line 88 of file TrigTauRecMerged.h.

{this, "InputTauJetContainer", "", "Input TauJet container"};

m_tauJetOutputDecorKeysArray

SG::WriteDecorHandleKeyArray<xAOD::TauJetContainer> TrigTauRecMerged::m_tauJetOutputDecorKeysArray

private

Definition at line 96 of file TrigTauRecMerged.h.

m_tauJetOutputKey

SG::WriteHandleKey<xAOD::TauJetContainer> TrigTauRecMerged::m_tauJetOutputKey {this, "OutputTauJetContainer", "", "Output TauJet container"}

private

Definition at line 95 of file TrigTauRecMerged.h.

{this, "OutputTauJetContainer", "", "Output TauJet container"};

m_tauSeedOutputKey

SG::WriteHandleKey<xAOD::JetContainer> TrigTauRecMerged::m_tauSeedOutputKey {this, "OutputJetSeed", "", "Output jets which are seeds for tau jets"}

private

Definition at line 94 of file TrigTauRecMerged.h.

{this, "OutputJetSeed", "", "Output jets which are seeds for tau jets"};

m_tauTrackInputKey

SG::ReadHandleKey<xAOD::TauTrackContainer> TrigTauRecMerged::m_tauTrackInputKey {this, "InputTauTrackContainer", "", "Input TauTrack container" }

private

Definition at line 90 of file TrigTauRecMerged.h.

{this, "InputTauTrackContainer", "", "Input TauTrack container" };

m_tauTrackOutputKey

SG::WriteHandleKey<xAOD::TauTrackContainer> TrigTauRecMerged::m_tauTrackOutputKey {this, "OutputTauTrackContainer", "", "Output TauTrack container"}

private

Definition at line 97 of file TrigTauRecMerged.h.

{this, "OutputTauTrackContainer", "", "Output TauTrack container"};

m_trackFinderTools

const ToolHandleArray<ITauToolBase> TrigTauRecMerged::m_trackFinderTools {this, "TrackFinderTools", {}, "Track Finder tools"}

private

Definition at line 73 of file TrigTauRecMerged.h.

{this, "TrackFinderTools", {}, "Track Finder tools"};

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vertexFinderTools

const ToolHandleArray<ITauToolBase> TrigTauRecMerged::m_vertexFinderTools {this, "VertexFinderTools", {}, "Vertex Finder tools"}

private

Definition at line 72 of file TrigTauRecMerged.h.

{this, "VertexFinderTools", {}, "Vertex Finder tools"};

m_vertexInputKey

SG::ReadHandleKey<xAOD::VertexContainer> TrigTauRecMerged::m_vertexInputKey {this, "InputVertexContainer", "", "Input vertex container"}

private

Definition at line 87 of file TrigTauRecMerged.h.

{this, "InputVertexContainer", "", "Input vertex container"};

m_vertexVarsTools

const ToolHandleArray<ITauToolBase> TrigTauRecMerged::m_vertexVarsTools {this, "VertexVarsTools", {}, "Vertex Variables tools"}

private

Definition at line 74 of file TrigTauRecMerged.h.

{this, "VertexVarsTools", {}, "Vertex Variables tools"};

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• TrigTauRecMerged.h
• TrigTauRecMerged.cxx