This is a tool to reconstruct the L1 Jet trigger sums for a particular RoI location from the stored JetInputs. More...

#include <L1JetElementTools.h>

Inheritance diagram for LVL1::L1JetElementTools:

Collaboration diagram for LVL1::L1JetElementTools:

Public Member Functions
	L1JetElementTools (const std::string &, const std::string &, const IInterface *)
	Constructor.
virtual	~L1JetElementTools ()
	default destructor
virtual StatusCode	initialize ()
	standard Athena-Algorithm method
virtual StatusCode	finalize ()
	standard Athena-Algorithm method
virtual void	makeJetElements (const xAOD::TriggerTowerContainer tts, xAOD::JetElementContainer jetElements) const
	Fill DataVector of JetElements from user-supplied TriggerTowers.
virtual void	mapJetElements (const xAOD::JetElementContainer jetElements, xAOD::JetElementMap_t map) const
	Fill map of JetElements from user-supplied vector of TriggerTowers.
virtual void	makeJetElements (const DataVector< TriggerTower > tts, DataVector< JetElement > jetElements) const
	Fill DataVector of JetElements from user-supplied vector of TriggerTowers.
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

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

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

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

Private Attributes
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

Static Private Attributes
static const int	m_ttSat = 0xFF
	Trigger Tower and Jet Element properties.
static const int	m_layerSat = 0x1FF

Detailed Description

This is a tool to reconstruct the L1 Jet trigger sums for a particular RoI location from the stored JetInputs.

Used for offline monitoring and trigger reconstruction.

Useage: L1JetElementTools->mapTowers() // set up event for analysis L1JetElementTools->formSums(eta, phi) // for each RoI L1JetElementTools->ET8x8() // etc. Return values

Author: Alan Watson Alan..nosp@m.Wats.nosp@m.on@ce.nosp@m.rn.c.nosp@m.h

Definition at line 38 of file L1JetElementTools.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ L1JetElementTools()

LVL1::L1JetElementTools::L1JetElementTools	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p )

Constructor.

Definition at line 24 of file L1JetElementTools.cxx.

  :
  AthAlgTool(t,n,p)
{
  declareInterface<IL1JetElementTools>(this);
 
}

◆ ~L1JetElementTools()

LVL1::L1JetElementTools::~L1JetElementTools ( )

virtual

default destructor

Destructor.

Definition at line 36 of file L1JetElementTools.cxx.

37{

38}

Member Function Documentation

◆ declareGaudiProperty()

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareProperty()

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

95{ return m_detStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ evtStore()

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

85{ return m_evtStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

◆ finalize()

StatusCode LVL1::L1JetElementTools::finalize ( )

virtual

standard Athena-Algorithm method

Finalisation.

Definition at line 50 of file L1JetElementTools.cxx.

{
  return StatusCode::SUCCESS;
}

◆ initialize()

StatusCode LVL1::L1JetElementTools::initialize ( )

virtual

standard Athena-Algorithm method

Initialisation.

Definition at line 43 of file L1JetElementTools.cxx.

{
  return StatusCode::SUCCESS;
}

◆ inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

◆ interfaceID()

const InterfaceID & LVL1::IL1JetElementTools::interfaceID ( )

inlinestaticinherited

Definition at line 50 of file IL1JetElementTools.h.

    { 
      return IID_IL1JetElementTools;
    }

◆ makeJetElements() [1/2]

void LVL1::L1JetElementTools::makeJetElements	(	const DataVector< TriggerTower > *	tts,
		DataVector< JetElement > *	jetElements ) const

virtual

Fill DataVector of JetElements from user-supplied vector of TriggerTowers.

Need to get pointers to the TriggerTowers corresponding to each JetElement
Store these in a vector for each JetElement
Keep vectors of pointers in a std::map so can easily locate and add more
towers to correct element

Each entry in the "Sums" map should now be a std::vector of pointers
to all of the towers in the collection making up each JetElement.
Now we need to go through them, add their energies, and form the
actual JetElement objects.
The complication is that the EM and Had vectors may not always have
the same lengths. In this case, we need to be careful how we combine
them.

At this point we have 2 vectors of ET. If their lengths and peak
positions are equal the rest is easy

If not, we align peaks, then create vectors using only those
slices before and after peaks where both layers present

Implements LVL1::IL1JetElementTools.

Definition at line 154 of file L1JetElementTools.cxx.

                                                                                                                      {
  
  // Clear collection before filling
  jetElements->clear();
  
  std::map< int, TriggerTowerVector_t> Sums;
 
  // Step over all TriggerTowers, and put into map
  DataVector<TriggerTower>::const_iterator it ;
  JetElementKey testKey(0.0, 0.0);
 
  for( it = tts->begin(); it != tts->end(); ++it ) {
    // Find JetElementKey for this TriggerTower
    int key = testKey.jeKey((*it)->phi(),(*it)->eta());
    // Does the map already contain an entry for this JetElement?
    std::map< int, TriggerTowerVector_t>::iterator mapIt=Sums.find(key);
    if (mapIt != Sums.end()) {
      // Add pointer to this tower to the list
      (mapIt->second).push_back((*it));
    }
    else {
      // New entry in map.
      TriggerTowerVector_t vec;
      vec.push_back((*it));
      Sums.insert(std::map< int, TriggerTowerVector_t >::value_type(key,vec));
    }
  } // end of loop over towers
 
    for (std::map< int, TriggerTowerVector_t >::iterator mapIt = Sums.begin();
         mapIt != Sums.end(); ++mapIt) {
 
      // Get first TT for this JE
      TriggerTowerVector_t::iterator it = (mapIt->second).begin();
      if (it != (mapIt->second).end()) {
        // Get JE eta, phi using first tower in vector (any tower in JE should do)
        unsigned int key = testKey.jeKey((*it)->phi(),(*it)->eta());
        double phi = testKey.phi();
        double eta = testKey.eta();
        // get peak positions in ET vectors
        int emPeak = (*it)->emPeak();
        int hadPeak = (*it)->hadPeak();
        // create empty ET vectors
        std::vector<int> emET;
        std::vector<int> hadET;
        // now loop through all TT present and add their ET values
        // if tower saturated, saturate corresponding sum also
        for (; it != (mapIt->second).end(); ++it) {
          std::vector<int> tempEM = (*it)->emLUT();
          if (tempEM.size() > emET.size()) emET.resize(tempEM.size(),0);
          for (size_t iem = 0; iem < tempEM.size(); ++iem)
            emET[iem] += (tempEM[iem] < m_ttSat ? tempEM[iem] : m_layerSat);
          std::vector<int> tempHad = (*it)->hadLUT();
          if (tempHad.size() > hadET.size()) hadET.resize(tempHad.size(),0);
          for (size_t ihd = 0; ihd < tempHad.size(); ++ihd)
            hadET[ihd] += (tempHad[ihd] < m_ttSat ? tempHad[ihd] : m_layerSat);
        }
        if (emET.size() == hadET.size() && emPeak == hadPeak) {
          std::vector<int> Dummy(emET.size());
          JetElement* jetElement = new JetElement(phi, eta, emET, hadET, key, Dummy, Dummy, Dummy, emPeak);
          jetElements->push_back(jetElement);
        }
        else {
          std::vector<int> emSums;
          std::vector<int> hadSums;
          int Peak = (emPeak < hadPeak ? emPeak : hadPeak);
          int Tail = (emET.size() - emPeak < hadET.size() - hadPeak ?
                      emET.size() - emPeak : hadET.size() - hadPeak) - 1;
          for (int j = emPeak - Peak; j <= emPeak + Tail; ++j)   emSums.push_back(emET[j]);
          for (int j = hadPeak - Peak; j <= hadPeak + Tail; ++j) hadSums.push_back(hadET[j]);
 
          std::vector<int> Dummy(emSums.size());
          JetElement* jetElement = new JetElement(phi, eta, emSums, hadSums, key, Dummy, Dummy, Dummy, Peak);
          jetElements->push_back(jetElement);
        } // unequal EM and Had readout
        
      } // end of check that first element of vector present
      
    } // end of loop through Sums map
     
  return;
}

◆ makeJetElements() [2/2]

void LVL1::L1JetElementTools::makeJetElements	(	const xAOD::TriggerTowerContainer *	tts,
		xAOD::JetElementContainer *	jetElements ) const

virtual

Fill DataVector of JetElements from user-supplied TriggerTowers.

Fill DataVector of JetElements from user-supplied vector of TriggerTowers.

Need to get pointers to the TriggerTowers corresponding to each JetElement
Store these in a vector for each JetElement
Keep vectors of pointers in a std::map so can easily locate and add more
towers to correct element

Each entry in the "Sums" map should now be a std::vector of pointers
to all of the towers in the collection making up each JetElement.
Now we need to go through them, add their energies, and form the
actual JetElement objects.
The complication is that the EM and Had vectors may not always have
the same lengths. In this case, we need to be careful how we combine
them.

Create JetElement

Implements LVL1::IL1JetElementTools.

Definition at line 56 of file L1JetElementTools.cxx.

                                                                                                                        {
 
  // Clear collection before filling
  jetElements->clear();
  
  std::map< int, std::vector<const xAOD::TriggerTower*> > Sums;
 
  // Step over all TriggerTowers, and put into map
  xAOD::TriggerTowerContainer::const_iterator it ;
  JetElementKey testKey(0.0, 0.0);
 
  for ( it = tts->begin(); it != tts->end(); ++it ) {
    // Find JetElementKey for this TriggerTower
    int key = testKey.jeKey((*it)->phi(),(*it)->eta());
    // Does the map already contain an entry for this JetElement?
    std::map< int, std::vector<const xAOD::TriggerTower*> >::iterator mapIt=Sums.find(key);
    if (mapIt != Sums.end()) {
      // Add pointer to this tower to the list
      (mapIt->second).push_back((*it));
    }
    else {
      // New entry in map.
      std::vector<const xAOD::TriggerTower*> vec;
      vec.push_back((*it));
      Sums.insert(std::map< int, std::vector<const xAOD::TriggerTower*> >::value_type(key,vec));
    }
  } // end of loop over towers
 
    for (std::map< int, std::vector<const xAOD::TriggerTower*> >::iterator mapIt = Sums.begin();
         mapIt != Sums.end(); ++mapIt) {
 
      // Get first TT for this JE
      std::vector<const xAOD::TriggerTower*>::iterator it = (mapIt->second).begin();
      if (it != (mapIt->second).end()) {
        // Get JE eta, phi using first tower in vector (any tower in JE should do)
        unsigned int key = testKey.jeKey((*it)->phi(),(*it)->eta());
        float phi = testKey.phi();
        float eta = testKey.eta();
        // peak position in ET vectors
        int peak = 0;
        // create empty ET vectors with 1 element each
        std::vector<uint16_t> emET(1);
        std::vector<uint16_t> hadET(1);
        // now loop through all TT present and add their ET values
        // if tower saturated, saturate corresponding sum also
        // Note that this version will only create 1 time sample per JE
        for (; it != (mapIt->second).end(); ++it) {
          int layer = (*it)->layer();
          int et = (*it)->jepET();
          if (layer == 0) emET[peak] += ( et < m_ttSat ? et : m_layerSat );
          else           hadET[peak] += ( et < m_ttSat ? et : m_layerSat );
        }        
        std::vector<uint32_t> Dummy(emET.size());
        xAOD::JetElement* jetElement = new xAOD::JetElement();
        jetElements->push_back(jetElement);
        jetElement->initialize(eta, phi, key, emET, hadET, Dummy, Dummy, Dummy, peak);
 
      } // end of check that first element of vector present
      
    } // end of loop through Sums map
     
  return;
}

◆ mapJetElements()

void LVL1::L1JetElementTools::mapJetElements	(	const xAOD::JetElementContainer *	jetElements,
		xAOD::JetElementMap_t *	map ) const

virtual

Fill map of JetElements from user-supplied vector of TriggerTowers.

Implements LVL1::IL1JetElementTools.

Definition at line 136 of file L1JetElementTools.cxx.

                                                                                                                 {
 
  // Clear map before filling
  map->clear();
  
  // Iterate through the DataVector and insert into the map
 
  for (xAOD::JetElementContainer::const_iterator it = jetElements->begin();
       it != jetElements->end(); ++it) {
    // As long as JetElement contains the key value, lazy (& quick) option is to use it
    int key = (*it)->key();
    map->insert(std::make_pair(key,(*it)));
  }
       
  return;
}

◆ msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

◆ outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

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

◆ renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

◆ sysInitialize()

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

overridevirtualinherited

Perform system initialization for an algorithm.

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

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

◆ sysStart()

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

overridevirtualinherited

Handle START transition.

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

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_layerSat

const int LVL1::L1JetElementTools::m_layerSat = 0x1FF

staticprivate

Definition at line 63 of file L1JetElementTools.h.

◆ m_ttSat

const int LVL1::L1JetElementTools::m_ttSat = 0xFF

staticprivate

Trigger Tower and Jet Element properties.

Definition at line 62 of file L1JetElementTools.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ L1JetElementTools()

◆ ~L1JetElementTools()

Member Function Documentation

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ evtStore()

◆ extraDeps_update_handler()

◆ finalize()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ makeJetElements() [1/2]

◆ makeJetElements() [2/2]

◆ mapJetElements()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_detStore

◆ m_evtStore

◆ m_layerSat

◆ m_ttSat

◆ m_varHandleArraysDeclared

◆ m_vhka