#include <TrigTauRecMerged.h>

Inheritance diagram for TrigTauRecMerged:

Collaboration 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. More...

virtual bool	isClonable () const override
	Specify if the algorithm is clonable. More...

virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...

virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm. More...

virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...

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`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

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

Private Types
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 T , class U , class V >
StatusCode	deepCopy (T &containerOut, U &containerStoreOut, const V dummyContainerType, const T &oldContainer)

template<class W , class V , class T >
StatusCode	deepCopy (W &writeHandle, const V dummyContainerType, const T &oldContainer) const

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
const ToolHandleArray< ITauToolBase >	m_commonTools {this, "ComTools", {}, "List of ITauToolBase tools"}
	internal tool store More...

const ToolHandleArray< ITauToolBase >	m_vertexFinderTools {this, "VFTools", {}, "Vertex Finder tools"}

const ToolHandleArray< ITauToolBase >	m_trackFinderTools {this, "TFTools", {}, "Track Finder tools"}

const ToolHandleArray< ITauToolBase >	m_vertexVarsTools {this, "VVTools", {}, "Vertex Variables tools"}

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

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

SG::ReadHandleKey< TrigRoiDescriptorCollection >	m_roIInputKey { this,"RoIInputKey","InputRoI","Input RoI name"}

SG::ReadHandleKey< xAOD::CaloClusterContainer >	m_clustersKey { this, "clustersKey", "CaloClusters", "caloclusters in view" }

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

SG::ReadHandleKey< xAOD::VertexContainer >	m_vertexKey { this, "Key_vertexInputContainer", "HLT_IDVertex_Tau", "input vertex container key"}

SG::ReadHandleKey< xAOD::TauJetContainer >	m_trigTauJetKey { this, "Key_trigTauJetInputContainer", "HLT_taujet", "input taujet container" }

SG::ReadHandleKey< xAOD::TauTrackContainer >	m_trigTauTrackInKey { this, "Key_trigTauTrackInputContainer", "HLT_tautrack_input", "input tautrack container" }

SG::WriteHandleKey< xAOD::JetContainer >	m_trigtauSeedOutKey { this,"Key_trigJetSeedOutputKey","HLT_jet_seed","Key for output jets which are seed for tau jets"}

SG::WriteHandleKey< xAOD::TauJetContainer >	m_trigtauRecOutKey {this,"Key_trigTauJetOutputContainer","HLT_taujet","Output taujet container"}

SG::WriteHandleKey< xAOD::TauTrackContainer >	m_trigtauTrkOutKey {this,"Key_trigTauTrackOutputContainer","HLT_tautrack","Output tautrack container"}

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

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 24 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

◆ TAUEFCALOMON

enum TrigTauRecMerged::TAUEFCALOMON

private

Enumerator
NoROIDescr
NoCellCont
EmptyCellCont
NoClustCont
NoClustKey
EmptyClustCont
NoJetAttach
NoHLTtauAttach
NoHLTtauDetAttach
NoHLTtauXdetAttach

Definition at line 41 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 54 of file TrigTauRecMerged.h.

                   {
     NoTrkCont=0,
     NoVtxCont=1
   };

Constructor & Destructor Documentation

◆ TrigTauRecMerged()

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

Definition at line 34 of file TrigTauRecMerged.cxx.

   :AthReentrantAlgorithm(name, pSvcLocator)
 {
 }

Member Function Documentation

◆ cardinality()

unsigned int AthReentrantAlgorithm::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 55 of file AthReentrantAlgorithm.cxx.

 {
   return 0;
 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyArrayType &
	)

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

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

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleType &
	)

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T > &	t,
		const SG::NotHandleType &
	)

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty	(	const std::string &	name,
		SG::VarHandleBase &	hndl,
		const std::string &	doc,
		const SG::VarHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKey &	hndl,
		const std::string &	doc,
		const SG::VarHandleKeyType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKeyArray &	hndArr,
		const std::string &	doc,
		const SG::VarHandleKeyArrayType &
	)

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc,
		const SG::NotHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(name, property, doc);
   }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc = `"none"`
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

   {
     typedef typename SG::HandleClassifier<T>::type htype;
     return declareProperty (name, property, doc, htype());
   }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

◆ deepCopy() [1/2]

template<class T , class U , class V >

StatusCode TrigTauRecMerged::deepCopy	(	T *&	containerOut,
		U *&	containerStoreOut,
		const V *	dummyContainerType,
		const T *&	oldContainer
	)

private

◆ deepCopy() [2/2]

template<class W , class V , class T >

StatusCode TrigTauRecMerged::deepCopy	(	W &	writeHandle,
		const V *	dummyContainerType,
		const T *&	oldContainer
	)		const

private

Definition at line 87 of file TrigTauRecMerged.h.

                                                                         {
    if(!writeHandle.isValid()){
       ATH_MSG_FATAL("Provided with an invalid write handle ");
       return StatusCode::FAILURE;
    }
    if(oldContainer != nullptr){
      for( const V* v : *oldContainer ){
        V* newV = new V();
        // Put objects into new container
        writeHandle->push_back(newV);
        // Copy across aux store
        *newV = *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.

95 { return m_detStore; }

◆ evtStore() [1/2]

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.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

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

overridevirtual

Definition at line 86 of file TrigTauRecMerged.cxx.

 {
   ATH_MSG_DEBUG("Execution");
  
   // variables to initialize and keep values for monitoring variables
   std::vector<unsigned char> calo_errors(0);
   std::vector<unsigned char> track_errors(0);
   //std::vector<const xAOD::CaloCluster *> clusters;
   std::vector<float> cluster_et_log(0);
   std::vector<float> cluster_dEta(0); 
   std::vector<float> cluster_dPhi(0);
   std::vector<float> cluster_log_SECOND_R(0);
   std::vector<float> cluster_SECOND_LAMBDA(0);
   std::vector<float> cluster_CENTER_LAMBDA(0);
   std::vector<float> track_pt_log(0);
   std::vector<float> track_dEta(0);
   std::vector<float> track_dPhi(0);
   std::vector<float> track_z0sinthetaTJVA_abs_log(0);
   std::vector<float> track_d0_abs_log(0);
   std::vector<float> track_nIBLHitsAndExp(0);
   std::vector<float> track_nPixelHitsPlusDeadSensors(0);
   std::vector<float> track_nSCTHitsPlusDeadSensors(0);
  
   auto nCells             = Monitored::Scalar<int>("nRoI_EFTauCells",    0);
   auto nTracks            = Monitored::Scalar<int>("nRoI_EFTauTracks", -10);
   auto dEta               = Monitored::Scalar<float>("dEtaEFTau_RoI",    -10.);
   auto dPhi               = Monitored::Scalar<float>("dPhiEFTau_RoI",    -10.);
   auto emRadius           = Monitored::Scalar<float>("EMRadius", -0.099);
   auto hadRadius          = Monitored::Scalar<float>("HadRadius", -0.099);
   auto EtFinal            = Monitored::Scalar<float>("EtFinal", 0.);
   auto Et                 = Monitored::Scalar<float>("Et", 0.);
   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 nWideTrk           = Monitored::Scalar<int>("nWideTrk",-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 numTrack           = Monitored::Scalar<int>("NTrk", -10);
   auto trkAvgDist         = Monitored::Scalar<float>("TrkAvgDist",-1.0);
   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 SumPtTrkFrac       = Monitored::Scalar<float>("SumPtTrkFrac",-999.9);
   auto innerTrkAvgDist    = Monitored::Scalar<float>("innerTrkAvgDist",-1.0);
   auto Ncand              = Monitored::Scalar<int>("nCand",0);
   auto EtaL1              = Monitored::Scalar<float>("EtaL1",-99.9);
   auto PhiL1              = Monitored::Scalar<float>("PhiL1",-99.9);
   auto EtaEF              = Monitored::Scalar<float>("EtaEF",-99.9);
   auto PhiEF              = Monitored::Scalar<float>("PhiEF",-99.9);
   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 ptDetectorAxis     = Monitored::Scalar<float>("ptDetectorAxis",-99.9);
   auto etaDetectorAxis     = Monitored::Scalar<float>("etaDetectorAxis",-99.9);
   auto ptDetectorAxis_log  = Monitored::Scalar<float>("ptDetectorAxis_log",-99.9);
   auto RNN_tracknumber    = Monitored::Scalar<int>("RNN_tracknumber",0);
   auto RNN_clusternumber  = Monitored::Scalar<int>("RNN_clusternumber",0); 
   auto RNNJetScore_0p         = Monitored::Scalar<float>("RNNJetScore_0p",-999);
   auto RNNJetScoreSigTrans_0p = Monitored::Scalar<float>("RNNJetScoreSigTrans_0p",-999);  
   auto RNNJetScore_1p         = Monitored::Scalar<float>("RNNJetScore_1p",-999);
   auto RNNJetScoreSigTrans_1p = Monitored::Scalar<float>("RNNJetScoreSigTrans_1p",-999);
   auto RNNJetScore_mp         = Monitored::Scalar<float>("RNNJetScore_mp",-999);
   auto RNNJetScoreSigTrans_mp = Monitored::Scalar<float>("RNNJetScoreSigTrans_mp",-999);
   auto EF_vertex_x      = Monitored::Scalar<float>("vertex_x", -999.9);
   auto EF_vertex_y      = Monitored::Scalar<float>("vertex_y", -999.9); 
   auto EF_vertex_z      = Monitored::Scalar<float>("vertex_z", -999.9);
  
   auto EF_calo_errors     = Monitored::Collection("calo_errors",calo_errors);
   auto EF_track_errors    = Monitored::Collection("track_errors",track_errors);
  
   auto Cluster_et_log     = Monitored::Collection("cluster_et_log", cluster_et_log);
   auto Cluster_dEta       = Monitored::Collection("cluster_dEta", cluster_dEta);
   auto Cluster_dPhi       = Monitored::Collection("cluster_dPhi", cluster_dPhi);
   auto Cluster_log_SECOND_R = Monitored::Collection("cluster_log_SECOND_R", cluster_log_SECOND_R);
   auto Cluster_SECOND_LAMBDA = Monitored::Collection("cluster_SECOND_LAMBDA", cluster_SECOND_LAMBDA);
   auto Cluster_CENTER_LAMBDA = Monitored::Collection("cluster_CENTER_LAMBDA", cluster_CENTER_LAMBDA);
  
   auto Track_pt_log = Monitored::Collection("track_pt_log", track_pt_log);
   auto Track_dEta = Monitored::Collection("track_dEta", track_dEta);
   auto Track_dPhi = Monitored::Collection("track_dPhi", track_dPhi);
   auto Track_z0sinthetaTJVA_abs_log = Monitored::Collection("track_z0sinthetaTJVA_abs_log", track_z0sinthetaTJVA_abs_log); 
   auto Track_d0_abs_log = Monitored::Collection("track_d0_abs_log", track_d0_abs_log);
   auto Track_nIBLHitsAndExp = Monitored::Collection("track_nIBLHitsAndExp", track_nIBLHitsAndExp);
   auto Track_nPixelHitsPlusDeadSensors = Monitored::Collection("track_nPixelHitsPlusDeadSensors", track_nPixelHitsPlusDeadSensors);
   auto Track_nSCTHitsPlusDeadSensors = Monitored::Collection("track_nSCTHitsPlusDeadSensors", track_nSCTHitsPlusDeadSensors);
  
   auto clustersMeanCenterLambda  = Monitored::Scalar<float>("clustersMeanCenterLambda",0.0);
   auto clustersMeanFirstEngDens  = Monitored::Scalar<float>("clustersMeanFirstEngDens",0.0);
   auto clustersMeanEMProbability  = Monitored::Scalar<float>("clustersMeanEMProbability",0.0);
   auto clustersMeanSecondLambda  = Monitored::Scalar<float>("clustersMeanSecondLambda",0.0);
   auto clustersMeanPresamplerFrac  = Monitored::Scalar<float>("clustersMeanPresamplerFrac",0.0);
  
   auto monitorIt = Monitored::Group( m_monTool, nCells, nTracks, dEta, dPhi, emRadius, hadRadius,
                    EtFinal, Et, EtHad, EtEm, EMFrac, IsoFrac, centFrac, nWideTrk, ipSigLeadTrk, trFlightPathSig, massTrkSys,
                    dRmax, numTrack, trkAvgDist, etovPtLead, PSSFraction, EMPOverTrkSysP, ChPiEMEOverCaloEME, SumPtTrkFrac,
                    innerTrkAvgDist, Ncand, EtaL1, PhiL1, EtaEF, PhiEF, mEflowApprox, ptRatioEflowApprox, pt_jetseed_log, 
                    ptDetectorAxis, etaDetectorAxis, ptDetectorAxis_log, RNN_clusternumber, RNNJetScore_0p, RNNJetScoreSigTrans_0p, 
                    RNNJetScore_1p, RNNJetScoreSigTrans_1p,RNNJetScore_mp,RNNJetScoreSigTrans_mp,Cluster_et_log, Cluster_dEta, Cluster_dPhi, Cluster_log_SECOND_R,
                    Cluster_SECOND_LAMBDA, Cluster_CENTER_LAMBDA, RNN_tracknumber, EF_vertex_x, EF_vertex_y, EF_vertex_z, EF_calo_errors, EF_track_errors, Track_pt_log, Track_dEta, Track_dPhi, Track_z0sinthetaTJVA_abs_log, Track_d0_abs_log, Track_nIBLHitsAndExp,
                    Track_nPixelHitsPlusDeadSensors, Track_nSCTHitsPlusDeadSensors, clustersMeanCenterLambda, clustersMeanFirstEngDens, clustersMeanEMProbability, clustersMeanSecondLambda, clustersMeanPresamplerFrac); 
  
  
   // Retrieve store.
   ATH_MSG_DEBUG("Executing TrigTauRecMerged");
  
   // Get RoiDescriptor
   SG::ReadHandle< TrigRoiDescriptorCollection > roisHandle = SG::makeHandle( m_roIInputKey, ctx );
   if ( not roisHandle.isValid() ) {
     ATH_MSG_ERROR("No roisHandle found");
     return StatusCode::FAILURE;
   }
  
   if(roisHandle->empty()){
     ATH_MSG_ERROR("RoIHandle size = Zero");
     return StatusCode::FAILURE;
   }
   const TrigRoiDescriptor *roiDescriptor = roisHandle->at(0);
  
   if ( roiDescriptor != nullptr ) {
     ATH_MSG_DEBUG(" RoI " << *roiDescriptor);
   }
   else {
     ATH_MSG_ERROR("Failed to find RoiDescriptor ");
     calo_errors.push_back(NoROIDescr);
     return StatusCode::FAILURE;
   }
  
   // get TauJetContainer from SG
   const xAOD::TauJetContainer *pTauContainer = nullptr;
   const xAOD::TauTrackContainer *pTauTrackContainer = nullptr;
  
   if (!m_trigTauJetKey.key().empty() && m_clustersKey.key().empty()) {
     SG::ReadHandle<xAOD::TauJetContainer> tauInputHandle(m_trigTauJetKey, ctx);
     pTauContainer = tauInputHandle.cptr();
     ATH_MSG_DEBUG("Input TauJet Container size: " << pTauContainer->size());
   }
  
   if (!m_trigTauTrackInKey.key().empty() && m_clustersKey.key().empty()) {
     SG::ReadHandle<xAOD::TauTrackContainer> tauTrackInputHandle(m_trigTauTrackInKey, ctx);
     pTauTrackContainer = tauTrackInputHandle.cptr();
     ATH_MSG_DEBUG("Tau Track Container Size " << pTauTrackContainer->size());
   }
  
   // Make new container which is deep copy of that
   std::unique_ptr<xAOD::TauJetContainer>    pContainer    = std::make_unique<xAOD::TauJetContainer>();
   std::unique_ptr<xAOD::TauJetAuxContainer> pAuxContainer = std::make_unique<xAOD::TauJetAuxContainer>();
   pContainer->setStore(pAuxContainer.get());
  
   // Write final taujets container
   SG::WriteHandle<xAOD::TauJetContainer> outputTauHandle(m_trigtauRecOutKey, ctx);
   ATH_CHECK(outputTauHandle.record(std::move(pContainer),std::move(pAuxContainer)));
  
   xAOD::TauJet* tau(nullptr);
   xAOD::TauJet* p_tau(nullptr);
  
   ATH_CHECK(deepCopy(outputTauHandle, tau, pTauContainer));
   if(!outputTauHandle->empty()) p_tau = outputTauHandle->back();
  
   std::unique_ptr<xAOD::TauTrackContainer>    pTrackContainer    = std::make_unique<xAOD::TauTrackContainer>();
   std::unique_ptr<xAOD::TauTrackAuxContainer> pTrackAuxContainer = std::make_unique<xAOD::TauTrackAuxContainer>();
   pTrackContainer->setStore(pTrackAuxContainer.get());
  
   SG::WriteHandle<xAOD::TauTrackContainer> tauTrackHandle(m_trigtauTrkOutKey, ctx);
   ATH_MSG_DEBUG("  write: " << tauTrackHandle.key() << " = " << "..." );
   ATH_CHECK(tauTrackHandle.record(std::move(pTrackContainer), std::move(pTrackAuxContainer)));
  
   xAOD::TauTrack* tautrack(nullptr);
  
   ATH_CHECK(deepCopy(tauTrackHandle, tautrack, pTauTrackContainer));
  
   if(!m_trigTauTrackInKey.key().empty() && m_clustersKey.key().empty()){
     p_tau->clearTauTrackLinks();
   }
  
   const xAOD::CaloClusterContainer *RoICaloClusterContainer = nullptr;
  
   if(!m_clustersKey.key().empty()){
     // Retrieve Calocluster container
     SG::ReadHandle< xAOD::CaloClusterContainer > CCContainerHandle = SG::makeHandle( m_clustersKey,ctx );
     CHECK( CCContainerHandle.isValid() );
  
     RoICaloClusterContainer = CCContainerHandle.get();
  
     if(RoICaloClusterContainer != nullptr) {
       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;
     }
  
     if(p_tau==nullptr){
       p_tau = new xAOD::TauJet();
       outputTauHandle->push_back(p_tau);
       p_tau->setROIWord(roiDescriptor->roiWord());
     }
  
     std::unique_ptr<xAOD::JetContainer> theJetContainer{std::make_unique<xAOD::JetContainer>()};
     std::unique_ptr<xAOD::JetAuxContainer> theJetAuxContainer{std::make_unique<xAOD::JetAuxContainer>()};
  
     theJetContainer->setStore(theJetAuxContainer.get());
  
     xAOD::Jet *aJet = new xAOD::Jet;
  
     theJetContainer->push_back(aJet);
     
     // Build the jet, also keep track of the kinematics by hand
     // We are using calibrated clusters, we need to keep track of this
     aJet->setConstituentsSignalState(xAOD::JetConstitScale::CalibratedJetConstituent);
     
     // Make a minimal effort to speed things up ;)
     // Eventually, want to use FastJet here?
     TLorentzVector myCluster;
     TLorentzVector TauBarycenter(0., 0., 0., 0.);
  
     xAOD::CaloClusterContainer::const_iterator clusterIt;
     for (clusterIt=RoICaloClusterContainer->begin(); clusterIt != RoICaloClusterContainer->end(); ++clusterIt) {
       ATH_MSG_DEBUG(" Cluster (e, eta, phi) : ("<< (*clusterIt)->e() << " , " <<(*clusterIt)->eta()<<" , "<<(*clusterIt)->phi()<< " )");
     
       if((*clusterIt)->e() < 0)
         {
       ATH_MSG_DEBUG(" Negative energy cluster is rejected");
       continue;
         }
  
       myCluster.SetPtEtaPhiE((*clusterIt)->pt(), (*clusterIt)->eta(), (*clusterIt)->phi(), (*clusterIt)->e());
       aJet->addConstituent(*clusterIt);
  
       TauBarycenter += myCluster;
     }
      
     aJet->setJetP4(xAOD::JetFourMom_t(TauBarycenter.Pt(), TauBarycenter.Eta(), TauBarycenter.Phi(), TauBarycenter.M() ) ); 
  
     ATH_MSG_DEBUG("jet formed"<< aJet->eta() <<" , " << aJet->phi() <<" , " << aJet->pt() << " , "<< aJet->e() );
  
     //-------------------------------------------------------------------------
     // using Jet collection
     // setup tau candidate structure
     //-------------------------------------------------------------------------
  
     p_tau->setJet(theJetContainer.get(), aJet);
  
     if(aJet->e()<=0) {
       ATH_MSG_DEBUG( "Roi: changing eta due to energy " << aJet->e() );
       p_tau->setP4(p_tau->pt(), roiDescriptor->eta(), roiDescriptor->phi(), p_tau->m());        
       ATH_MSG_DEBUG( "Roi: " << roiDescriptor->roiId() << " Tau eta: " << p_tau->eta() << " Tau phi: " << p_tau->phi() );
     }
  
     // Save Outputs
     SG::WriteHandle< xAOD::JetContainer > outTauSeedHandle = SG::makeHandle( m_trigtauSeedOutKey,ctx );
     CHECK( outTauSeedHandle.record( std::move( theJetContainer ), std::move( theJetAuxContainer ) ) );
   }
  
   //Check if jetLink is valid for all taus
   CHECK(p_tau->jetLink().isValid());
  
   // get TrackContainer
   if(!m_tracksKey.key().empty()){
     SG::ReadHandle< xAOD::TrackParticleContainer > TPContainerHandle = SG::makeHandle( m_tracksKey,ctx );
  
     const xAOD::TrackParticleContainer *RoITrackParticleContainer = nullptr;
  
     if ( !TPContainerHandle.isValid() ) {
       ATH_MSG_DEBUG(" No Track container found.");
       track_errors.push_back(NoTrkCont);
     }
     else {
       RoITrackParticleContainer = TPContainerHandle.get();
       ATH_MSG_DEBUG(" Size of vector Track container is " << RoITrackParticleContainer->size());
       if(RoITrackParticleContainer != nullptr) nTracks = RoITrackParticleContainer->size();
     }
   }
  
   const xAOD::VertexContainer* RoIVxContainer = nullptr;
  
   // get Vertex Container
   if(!m_vertexKey.key().empty()){
     SG::ReadHandle< xAOD::VertexContainer > VertexContainerHandle = SG::makeHandle( m_vertexKey,ctx );
  
     if( !VertexContainerHandle.isValid() ) {
       ATH_MSG_DEBUG(" No VxContainers retrieved for the trigger element");
       track_errors.push_back(NoVtxCont);
     }
     else {
       RoIVxContainer = VertexContainerHandle.get();
       ATH_MSG_DEBUG(" Size of vector Vertex  container " << RoIVxContainer->size());
     }
   }
  
   ATH_MSG_DEBUG(" roidescriptor roiword " << roiDescriptor->roiWord() << " saved " << p_tau->ROIWord() );
  
   //-------------------------------------------------------------------------
   // loop over booked tau tools
   //-------------------------------------------------------------------------
  
   StatusCode processStatus = StatusCode::SUCCESS;
   
   // dummy container passed to TauVertexVariables, not used in trigger though
   xAOD::VertexContainer dummyVxCont; 
  
   for (const auto& tool: m_vertexFinderTools){
     ATH_MSG_DEBUG("Starting Tool: " <<  tool->name() );
  
     processStatus = tool->executeVertexFinder(*p_tau,RoIVxContainer);
  
     if ( !processStatus.isFailure() ) {
       ATH_MSG_DEBUG(" "<< tool->name() << " executed successfully ");
     } else {
       ATH_MSG_DEBUG(" "<< tool->name() << " execution failed ");
       break;
     }
  
   }
  
   if ( !processStatus.isFailure() ) {
      ATH_MSG_DEBUG("Starting Tool: " <<  m_commonTools[0]->name() );
      processStatus = m_commonTools[0]->execute( *p_tau);
   }
  
   for (const auto& tool: m_trackFinderTools){
  
     if( !processStatus.isFailure() ) ATH_MSG_DEBUG("Starting Tool: " <<  tool->name() );
     else break;
  
     processStatus = tool->executeTrackFinder(*p_tau, *tauTrackHandle);
  
     if ( !processStatus.isFailure() ) {
       ATH_MSG_DEBUG(" "<< tool->name() << " executed successfully ");
     } else {
       ATH_MSG_DEBUG(" "<< tool->name() << " execution failed ");
       break;
     }
  
   }
  
   for (unsigned int i=1; i< m_commonTools.size(); i++) {
  
     if( !processStatus.isFailure() ) ATH_MSG_DEBUG("Starting Tool: " <<  m_commonTools[i]->name() );
     else break;
  
     processStatus = m_commonTools[i]->execute(*p_tau);
  
     if ( !processStatus.isFailure() ) {
       ATH_MSG_DEBUG(" "<< m_commonTools[i]->name() << " executed successfully ");
     } else {
       ATH_MSG_DEBUG(" "<< m_commonTools[i]->name() << " execution failed ");
       break;
     }
  
   }
  
   for (const auto& tool: m_vertexVarsTools){
  
     if( !processStatus.isFailure() ) ATH_MSG_DEBUG("Starting Tool: " <<  tool->name() );
     else break;
  
     processStatus = tool->executeVertexVariables(*p_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(*p_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 " << p_tau->allTracks() << " tracks linked");
  
   //check status
   if ( !processStatus.isSuccess() )  {   // some problem
     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 {
  
     float fJetEnergy = (*p_tau->jetLink())->e();
     ATH_MSG_DEBUG("Roi: jet e "<< fJetEnergy);
       
     if( fJetEnergy < 0.00001 ) {
       ATH_MSG_DEBUG("Roi: changing eta phi to L1 ones due to energy negative (PxPyPzE flips eta and phi)");
       ATH_MSG_DEBUG("Roi: this is probably not needed anymore, method PxPyPzE has been corrected");
         
       p_tau->setP4(p_tau->pt(), roiDescriptor->eta(), roiDescriptor->phi(), p_tau->m());
         
       ATH_MSG_DEBUG("Roi: " << roiDescriptor->roiId()
             << " Tau eta: " << p_tau->eta()
             << " Tau phi: " << p_tau->phi()
             << " Tau pT : "<< p_tau->pt());
     }
  
     // get tau detail variables for Monitoring
  
     numTrack = p_tau->nTracks();
     nWideTrk = p_tau->nTracksIsolation();
     p_tau->detail(xAOD::TauJetParameters::numCells, nCells);
     p_tau->detail(xAOD::TauJetParameters::nChargedTracks, nTracks);
     p_tau->detail(xAOD::TauJetParameters::trkAvgDist, trkAvgDist);
     p_tau->detail(xAOD::TauJetParameters::etOverPtLeadTrk, etovPtLead);
     p_tau->detail(xAOD::TauJetParameters::EMRadius, emRadius);
     p_tau->detail(xAOD::TauJetParameters::hadRadius, hadRadius);
     p_tau->detail(xAOD::TauJetParameters::isolFrac, IsoFrac);
     p_tau->detail(xAOD::TauJetParameters::centFrac, centFrac);
     p_tau->detail(xAOD::TauJetParameters::ipSigLeadTrk, ipSigLeadTrk);
     p_tau->detail(xAOD::TauJetParameters::trFlightPathSig, trFlightPathSig);
     p_tau->detail(xAOD::TauJetParameters::dRmax, dRmax);
     p_tau->detail(xAOD::TauJetParameters::massTrkSys, massTrkSys);
       
     p_tau->detail(xAOD::TauJetParameters::PSSFraction, PSSFraction);
     p_tau->detail(xAOD::TauJetParameters::EMPOverTrkSysP, EMPOverTrkSysP);
     p_tau->detail(xAOD::TauJetParameters::ChPiEMEOverCaloEME, ChPiEMEOverCaloEME);
     p_tau->detail(xAOD::TauJetParameters::innerTrkAvgDist, innerTrkAvgDist);     
     p_tau->detail(xAOD::TauJetParameters::SumPtTrkFrac, SumPtTrkFrac);
  
     massTrkSys /= Gaudi::Units::GeV;
     p_tau->detail(xAOD::TauJetParameters::etEMAtEMScale, EtEm);
     EtEm /= Gaudi::Units::GeV;
     p_tau->detail(xAOD::TauJetParameters::etHadAtEMScale, EtHad);
     EtHad /= Gaudi::Units::GeV;
     Et            = EtEm + EtHad;
     EtFinal       = p_tau->pt()/Gaudi::Units::GeV;
     
     EtaEF = p_tau->eta();
     PhiEF = p_tau->phi();
       
     if( Et !=0) EMFrac =  EtEm/ Et ;
  
     EtaL1 = roiDescriptor->eta();
     PhiL1 = roiDescriptor->phi();     
     dEta =  EtaEF - roiDescriptor->eta();
     dPhi =  PhiEF - roiDescriptor->phi();
     if(dPhi<-M_PI) dPhi += 2.0*M_PI;
     if(dPhi>M_PI)  dPhi -= 2.0*M_PI;
  
    
     float pre_mEflowApprox;
     p_tau->detail(xAOD::TauJetParameters::mEflowApprox, pre_mEflowApprox);  
     mEflowApprox = std::log10(std::max(pre_mEflowApprox, 140.0f));
  
     float pre_ptRatioEflowApprox;
     p_tau->detail(xAOD::TauJetParameters::ptRatioEflowApprox, pre_ptRatioEflowApprox);
     ptRatioEflowApprox = std::min(pre_ptRatioEflowApprox, 4.0f);
     
     pt_jetseed_log  = std::log10(p_tau->ptJetSeed());
  
     ptDetectorAxis  =  std::min(p_tau->ptDetectorAxis() / 1000.0, 10000.0);
     ptDetectorAxis_log  =  std::log10(std::min(p_tau->ptDetectorAxis() / 1000.0, 10000.0));
    
     etaDetectorAxis  = p_tau->etaDetectorAxis();   
  
     if( p_tau->nTracks() > 0 ) {
       ipSigLeadTrk = std::abs(p_tau->track(0)->d0SigTJVA());
     }
  
     // track variables monitoring 
     for( auto track : p_tau->allTracks()){
     
         RNN_tracknumber += 1;
         track_pt_log.push_back(std::log10( track->pt()));
         track_dEta.push_back(track->eta()- p_tau->eta()); 
         track_dPhi.push_back(track->p4().DeltaPhi(p_tau->p4()));
         track_z0sinthetaTJVA_abs_log.push_back(track->z0sinthetaTJVA());
         track_d0_abs_log.push_back(std::log10( std::abs(track->track()->d0()) + 1e-6));
  
         uint8_t inner_pixel_hits, inner_pixel_exp;                    
         const auto success1_innerPixel_hits = track->track()->summaryValue(inner_pixel_hits, xAOD::numberOfInnermostPixelLayerHits);                        
         const auto success2_innerPixel_exp = track->track()->summaryValue(inner_pixel_exp, xAOD::expectInnermostPixelLayerHit);                                       
         float nIBLHitsAndExp = -999;                                              
         if (success1_innerPixel_hits && success2_innerPixel_exp) {nIBLHitsAndExp=inner_pixel_exp ? inner_pixel_hits : 1.;};        
         track_nIBLHitsAndExp.push_back(nIBLHitsAndExp);
  
         uint8_t pixel_hits, pixel_dead;                                    
         const auto success1_pixel_hits = track->track()->summaryValue(pixel_hits, xAOD::numberOfPixelHits);          
         const auto 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 auto success1_sct_hits = track->track()->summaryValue(sct_hits, xAOD::numberOfSCTHits);                   
         const auto 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);
     }
  
     RNN_clusternumber = p_tau->clusters().size();
  
     // cluster variables monitoring
     for ( auto cluster : p_tau->clusters()){
  
         auto cls = dynamic_cast<const xAOD::CaloCluster*>(cluster);
  
         cluster_et_log.push_back(std::log10( cls->et()));
         cluster_dEta.push_back( cls->eta()- p_tau->eta());
         cluster_dPhi.push_back( cls ->p4().DeltaPhi(p_tau->p4()));
      
         double log_second_R = -999.;
         const auto 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 auto 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 auto 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);     
     }
  
     // monitoring tau vertex
     if( p_tau->vertexLink().isValid() && p_tau->vertex() && p_tau->vertex()->vertexType() != xAOD::VxType::NoVtx ){
         EF_vertex_x = p_tau->vertex()->x();
         EF_vertex_y = p_tau->vertex()->y();       
         EF_vertex_z = p_tau->vertex()->z();
     }
  
     ATH_MSG_DEBUG(" Roi: " << roiDescriptor->roiId()
           << " Tau being saved eta: " << EtaEF << " Tau phi: " << PhiEF
           << " wrt L1 dEta "<< dEta<<" dPhi "<<dPhi
           << " Tau Et (GeV): "<< EtFinal);
       
     // monitor RNN score
     if(p_tau->hasDiscriminant(xAOD::TauJetParameters::RNNJetScore)){
         if(p_tau->nTracks() == 0){
            RNNJetScore_0p = p_tau->discriminant(xAOD::TauJetParameters::RNNJetScore);
         } else if (p_tau->nTracks() == 1) {
            RNNJetScore_1p = p_tau->discriminant(xAOD::TauJetParameters::RNNJetScore);
         } else {
            RNNJetScore_mp = p_tau->discriminant(xAOD::TauJetParameters::RNNJetScore);
         }
     }
    
     if(p_tau->hasDiscriminant(xAOD::TauJetParameters::RNNJetScoreSigTrans)){
         if(p_tau->nTracks() == 0){
            RNNJetScoreSigTrans_0p = p_tau->discriminant(xAOD::TauJetParameters::RNNJetScoreSigTrans);
         } else if (p_tau->nTracks() == 1){
            RNNJetScoreSigTrans_1p = p_tau->discriminant(xAOD::TauJetParameters::RNNJetScoreSigTrans);
         } else {
            RNNJetScoreSigTrans_mp = p_tau->discriminant(xAOD::TauJetParameters::RNNJetScoreSigTrans);
         }
     }
  
     // monitor BRT variables
  
     // use dummy variable to temporarily store results retrieved from tau EDM before final assignment
     float avariable = 0.;
  
     bool test = p_tau->detail(xAOD::TauJetParameters::ClustersMeanCenterLambda, avariable);
     if(test) clustersMeanCenterLambda = avariable;
  
     test = p_tau->detail(xAOD::TauJetParameters::ClustersMeanFirstEngDens, avariable);
     if(test) clustersMeanFirstEngDens = avariable;
  
     test = p_tau->detail(xAOD::TauJetParameters::ClustersMeanEMProbability, avariable);
     if(test) clustersMeanEMProbability = avariable;
  
     test = p_tau->detail(xAOD::TauJetParameters::ClustersMeanSecondLambda, avariable);
     if(test) clustersMeanSecondLambda = avariable;
  
     test = p_tau->detail(xAOD::TauJetParameters::ClustersMeanPresamplerFrac, avariable);
     if(test) clustersMeanPresamplerFrac = avariable;
  
     ++Ncand;
   }
  
   //-------------------------------------------------------------------------
   // all done, register the tau Container in TDS.
   //-------------------------------------------------------------------------
   
   ATH_MSG_DEBUG("Output TauJetContainer size:"<< outputTauHandle->size());
   ATH_MSG_DEBUG("Output TauJetTrackContainer size:"<< tauTrackHandle->size());
   
   
   ATH_MSG_DEBUG("Recorded a tau container: HLT_TrigTauRecMerged");
   ATH_MSG_DEBUG("the tau object has been registered in the tau container");
   
   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 & AthReentrantAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 79 of file AthReentrantAlgorithm.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 Algorithm::extraOutputDeps();
 }

◆ filterPassed()

virtual bool AthReentrantAlgorithm::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 135 of file AthReentrantAlgorithm.h.

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

◆ initialize()

StatusCode TrigTauRecMerged::initialize ( )

overridevirtual

Definition at line 39 of file TrigTauRecMerged.cxx.

 {
   ATH_MSG_DEBUG("TrigTauRecMerged::initialize()");
  
   if ( m_commonTools.begin() == m_commonTools.end() ) {
     ATH_MSG_ERROR(" no tools given for this algorithm.");
     return StatusCode::FAILURE;
   }
  
   for(const auto& tool : m_commonTools) {
     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 ( not m_monTool.name().empty() ) {
     ATH_CHECK( m_monTool.retrieve() );
   }
   
   ATH_MSG_DEBUG("Initialising Handle Keys");
   ATH_CHECK(m_roIInputKey.initialize());
   ATH_CHECK(m_clustersKey.initialize(SG::AllowEmpty));
   ATH_CHECK(m_tracksKey.initialize(SG::AllowEmpty));
   ATH_CHECK(m_vertexKey.initialize(SG::AllowEmpty));
   ATH_CHECK(m_trigTauJetKey.initialize(SG::AllowEmpty));
   ATH_CHECK(m_trigTauTrackInKey.initialize(SG::AllowEmpty));
   ATH_CHECK(m_trigtauSeedOutKey.initialize());
   ATH_CHECK(m_trigtauRecOutKey.initialize());
   ATH_CHECK(m_trigtauTrkOutKey.initialize());
  
   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 AthReentrantAlgorithm::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

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

Definition at line 44 of file AthReentrantAlgorithm.cxx.

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

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ 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 AthReentrantAlgorithm::setFilterPassed	(	bool	state,
		const EventContext &	ctx
	)		const

inlinevirtualinherited

Definition at line 139 of file AthReentrantAlgorithm.h.

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

◆ sysExecute()

StatusCode AthReentrantAlgorithm::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 67 of file AthReentrantAlgorithm.cxx.

 {
   return Gaudi::Algorithm::sysExecute (ctx);
 }

◆ sysInitialize()

StatusCode AthReentrantAlgorithm::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 InputMakerBase, and HypoBase.

Definition at line 96 of file AthReentrantAlgorithm.cxx.

                                                 {
   StatusCode sc=AthCommonDataStore<AthCommonMsg<Gaudi::Algorithm>>::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_clustersKey

SG::ReadHandleKey< xAOD::CaloClusterContainer > TrigTauRecMerged::m_clustersKey { this, "clustersKey", "CaloClusters", "caloclusters in view" }

private

Definition at line 73 of file TrigTauRecMerged.h.

◆ m_commonTools

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

private

internal tool store

Definition at line 60 of file TrigTauRecMerged.h.

◆ 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 AthReentrantAlgorithm::m_extendedExtraObjects

privateinherited

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

Empty if no symlinks were found.

Definition at line 153 of file AthReentrantAlgorithm.h.

◆ m_idTools

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

private

Definition at line 64 of file TrigTauRecMerged.h.

◆ m_monTool

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

private

Definition at line 69 of file TrigTauRecMerged.h.

◆ m_roIInputKey

SG::ReadHandleKey< TrigRoiDescriptorCollection > TrigTauRecMerged::m_roIInputKey { this,"RoIInputKey","InputRoI","Input RoI name"}

private

Definition at line 72 of file TrigTauRecMerged.h.

◆ m_trackFinderTools

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

private

Definition at line 62 of file TrigTauRecMerged.h.

◆ m_tracksKey

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

private

Definition at line 74 of file TrigTauRecMerged.h.

◆ m_trigTauJetKey

SG::ReadHandleKey< xAOD::TauJetContainer> TrigTauRecMerged::m_trigTauJetKey { this, "Key_trigTauJetInputContainer", "HLT_taujet", "input taujet container" }

private

Definition at line 76 of file TrigTauRecMerged.h.

◆ m_trigtauRecOutKey

SG::WriteHandleKey< xAOD::TauJetContainer > TrigTauRecMerged::m_trigtauRecOutKey {this,"Key_trigTauJetOutputContainer","HLT_taujet","Output taujet container"}

private

Definition at line 80 of file TrigTauRecMerged.h.

◆ m_trigtauSeedOutKey

SG::WriteHandleKey< xAOD::JetContainer > TrigTauRecMerged::m_trigtauSeedOutKey { this,"Key_trigJetSeedOutputKey","HLT_jet_seed","Key for output jets which are seed for tau jets"}

private

Definition at line 79 of file TrigTauRecMerged.h.

◆ m_trigTauTrackInKey

SG::ReadHandleKey< xAOD::TauTrackContainer> TrigTauRecMerged::m_trigTauTrackInKey { this, "Key_trigTauTrackInputContainer", "HLT_tautrack_input", "input tautrack container" }

private

Definition at line 77 of file TrigTauRecMerged.h.

◆ m_trigtauTrkOutKey

SG::WriteHandleKey< xAOD::TauTrackContainer > TrigTauRecMerged::m_trigtauTrkOutKey {this,"Key_trigTauTrackOutputContainer","HLT_tautrack","Output tautrack container"}

private

Definition at line 81 of file TrigTauRecMerged.h.

◆ 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, "VFTools", {}, "Vertex Finder tools"}

private

Definition at line 61 of file TrigTauRecMerged.h.

◆ m_vertexKey

SG::ReadHandleKey< xAOD::VertexContainer> TrigTauRecMerged::m_vertexKey { this, "Key_vertexInputContainer", "HLT_IDVertex_Tau", "input vertex container key"}

private

Definition at line 75 of file TrigTauRecMerged.h.

◆ m_vertexVarsTools

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

private

Definition at line 63 of file TrigTauRecMerged.h.

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

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Member Enumeration Documentation

◆ TAUEFCALOMON

◆ TAUEFTRKMON

Constructor & Destructor Documentation

◆ TrigTauRecMerged()

Member Function Documentation

◆ cardinality()

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ deepCopy() [1/2]

◆ deepCopy() [2/2]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ filterPassed()

◆ initialize()

◆ inputHandles()

◆ isClonable()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ setFilterPassed()

◆ sysExecute()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_clustersKey

◆ m_commonTools

◆ m_detStore

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_idTools

◆ m_monTool

◆ m_roIInputKey

◆ m_trackFinderTools

◆ m_tracksKey

◆ m_trigTauJetKey

◆ m_trigtauRecOutKey

◆ m_trigtauSeedOutKey

◆ m_trigTauTrackInKey

◆ m_trigtauTrkOutKey

◆ m_varHandleArraysDeclared

◆ m_vertexFinderTools

◆ m_vertexKey

◆ m_vertexVarsTools

◆ m_vhka