#include <L1CPMTowerTools.h>

Inheritance diagram for LVL1::L1CPMTowerTools:

Collaboration diagram for LVL1::L1CPMTowerTools:

Public Types
typedef DataVector< TriggerTower >	TTCollection

typedef DataVector< xAOD::TriggerTower >	xAODTTCollection

Public Member Functions
	L1CPMTowerTools (const std::string &, const std::string &, const IInterface *)
	Constructor. More...

virtual	~L1CPMTowerTools ()
	default destructor More...

virtual StatusCode	initialize ()
	standard Athena-Algorithm method More...

virtual StatusCode	finalize ()
	standard Athena-Algorithm method More...

virtual void	makeCPMTowers (const DataVector< TriggerTower > triggerTowers, DataVector< CPMTower > cpmTowers, bool zeroSuppress=true)
	Fill DataVector of CPMTowers from user-supplied TriggerTowers. More...

virtual void	makeCPMTowers (const DataVector< xAOD::TriggerTower > triggerTowers, DataVector< CPMTower > cpmTowers, bool zeroSuppress=true)
	Fill DataVector of CPMTowers from user-supplied xAOD::TriggerTowers. More...

virtual void	makeCPMTowers (const DataVector< xAOD::TriggerTower > triggerTowers, DataVector< xAOD::CPMTower > cpmTowers, bool zeroSuppress=true)
	Fill DataVector of CPMTowers from user-supplied xAOD::TriggerTowers. More...

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	sysInitialize () override
	Perform system initialization for an algorithm. 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

Static Public Member Functions
static const InterfaceID &	interfaceID ()

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
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
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
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

This is a tool to form CPMTowers from TriggerTowers or to copy CPMTowers into a standard map, using keys from TrigT1Utils::TriggerTowerKey

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

Definition at line 34 of file L1CPMTowerTools.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

◆ TTCollection

typedef DataVector<TriggerTower> LVL1::L1CPMTowerTools::TTCollection

Definition at line 37 of file L1CPMTowerTools.h.

◆ xAODTTCollection

typedef DataVector<xAOD::TriggerTower> LVL1::L1CPMTowerTools::xAODTTCollection

Definition at line 38 of file L1CPMTowerTools.h.

Constructor & Destructor Documentation

◆ L1CPMTowerTools()

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

Constructor.

Definition at line 21 of file L1CPMTowerTools.cxx.

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

◆ ~L1CPMTowerTools()

LVL1::L1CPMTowerTools::~L1CPMTowerTools ( )

virtual

default destructor

Destructor.

Definition at line 33 of file L1CPMTowerTools.cxx.

34 {

35 }

Member Function Documentation

◆ declareGaudiProperty() [1/4]

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

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

◆ evtStore() [1/2]

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; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::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; }

◆ 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::L1CPMTowerTools::finalize ( )

virtual

standard Athena-Algorithm method

Finalisation.

Definition at line 47 of file L1CPMTowerTools.cxx.

 {
   return StatusCode::SUCCESS;
 }

◆ initialize()

StatusCode LVL1::L1CPMTowerTools::initialize ( )

virtual

standard Athena-Algorithm method

Initialisation.

Definition at line 40 of file L1CPMTowerTools.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::IL1CPMTowerTools::interfaceID ( )

inlinestaticinherited

Definition at line 38 of file IL1CPMTowerTools.h.

     { 
       return IID_IL1CPMTowerTools;
     }

◆ makeCPMTowers() [1/3]

void LVL1::L1CPMTowerTools::makeCPMTowers	(	const DataVector< TriggerTower > *	triggerTowers,
		DataVector< CPMTower > *	cpmTowers,
		bool	zeroSuppress = `true`
	)

virtual

Fill DataVector of CPMTowers from user-supplied TriggerTowers.

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

Need to get pointers to the TriggerTowers corresponding to each CPMTower
Store these in a vector for each CPMTower
Keep vectors of pointers in a std::map so can easily locate and add more
towers to correct element Right now this is redundant, but when TT becomes a single-layer object
we'll need this bit

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 CPMTower.
Now we need to go through them and add their energies to the appropriate
layer to form the actual CPMTower 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.

May not have signals in both layers. If not, fill empty layer with null data. TODO when have single-layer TriggerTower as input

Now we can create the tower object

Implements LVL1::IL1CPMTowerTools.

Definition at line 53 of file L1CPMTowerTools.cxx.

                                                                                                                                     {
  
   // Clear collection before filling
   cpmTowers->clear();
   
   std::map< int, TriggerTowerVector_t > Sums;
  
   // Step over all TriggerTowers, and put into map
   TTCollection::const_iterator it ;
   TriggerTowerKey testKey(0.0, 0.0);
  
   for( it = triggerTowers->begin(); it != triggerTowers->end(); ++it ) {
     // Check within CPM tower coverage
     if (fabs((*it)->eta()) > 2.5) continue;
     
     // Find TriggerTowerKey for this TriggerTower
     int key = testKey.ttKey((*it)->phi(),(*it)->eta());
     // Does the map already contain an entry for this CPMTower?
     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 CPMT
       TriggerTowerVector_t::iterator it = (mapIt->second).begin();
       if (it != (mapIt->second).end()) {
         // Get CPMT eta, phi using first tower in vector (either tower in CPMT should do)
         double phi = (*it)->phi();
         double eta = (*it)->eta();
         // create empty result vectors 
         std::vector<int> emET;
         std::vector<int> hadET;
         std::vector<int> emError;
         std::vector<int> hadError;
         // Both should contain same number of samples, and same peak position
         int Peak = 0;
         // now loop through all TT present and add their ET values to CPMT
         for (; it != (mapIt->second).end(); ++it) {
           // Right now there should only be one TT. Add logic for dual layers later.
           Peak = (*it)->emPeak();
           emET = (*it)->emLUT();
           hadET = (*it)->hadLUT();
           emError.assign(emET.size(),(*it)->emError());
           hadError.assign(hadET.size(),(*it)->hadError());
         }
         
         if (!zeroSuppress || emET[Peak] > 0 || hadET[Peak] > 0) {
           CPMTower* cpmTower = new CPMTower(phi, eta, emET, emError, hadET, hadError, Peak);
           cpmTowers->push_back(cpmTower);
         }
       } // end of check that first element of vector present
       
     } // end of loop through Sums map
      
   return;
 }

◆ makeCPMTowers() [2/3]

void LVL1::L1CPMTowerTools::makeCPMTowers	(	const DataVector< xAOD::TriggerTower > *	triggerTowers,
		DataVector< CPMTower > *	cpmTowers,
		bool	zeroSuppress = `true`
	)

virtual

Fill DataVector of CPMTowers from user-supplied xAOD::TriggerTowers.

Fill DataVector of CPMTowers from user-supplied vector of xAOD::TriggerTowers.

Need to get pointers to the TriggerTowers corresponding to each CPMTower
Store these in a vector for each CPMTower
Keep vectors of pointers in a std::map so can easily locate and add more
towers to correct element Right now this is redundant, but when TT becomes a single-layer object
we'll need this bit

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 CPMTower.
Now we need to go through them and add their energies to the appropriate
layer to form the actual CPMTower 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.

Now we can create the tower object

Implements LVL1::IL1CPMTowerTools.

Definition at line 141 of file L1CPMTowerTools.cxx.

                                                                                                                                         {
  
   // Clear collection before filling
   cpmTowers->clear();
   
   std::map< int, xAODTriggerTowerVector_t > Sums;
  
   // Step over all TriggerTowers, and put into map
   xAODTTCollection::const_iterator it ;
   TriggerTowerKey testKey(0.0, 0.0);
  
   for( it = triggerTowers->begin(); it != triggerTowers->end(); ++it ) {
     // Check within CPM tower coverage
     if (fabs((*it)->eta()) > 2.5) continue;
     
     // Find TriggerTowerKey for this TriggerTower
     int key = testKey.ttKey((*it)->phi(),(*it)->eta());
     // Does the map already contain an entry for this CPMTower?
     std::map< int, xAODTriggerTowerVector_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.
       xAODTriggerTowerVector_t vec;
       vec.push_back((*it));
       Sums.insert(std::map< int, xAODTriggerTowerVector_t >::value_type(key,vec));
     }
   } // end of loop over towers
  
     for (std::map< int, xAODTriggerTowerVector_t >::iterator mapIt = Sums.begin();
          mapIt != Sums.end(); ++mapIt) {
  
       // create empty result vectors containing a single element
       // That way if one layer is missing from TT collection we have a 0 entered already
       std::vector<int> emET(1);
       std::vector<int> hadET(1);
       std::vector<int> emError(1);
       std::vector<int> hadError(1);
       // Both should contain same number of samples, and same peak position
       int Peak = 0;
       
       // Get first TT for this CPMT
       xAODTriggerTowerVector_t::iterator it = (mapIt->second).begin();
       if (it != (mapIt->second).end()) {
         // Get CPMT eta, phi using first tower in vector (either tower in CPMT should do)
         double phi = (*it)->phi();
         double eta = (*it)->eta();
         // now loop through all TT present and add their ET values to CPMT
         for (; it != (mapIt->second).end(); ++it) {
           // Going to play safe and fill just one entry into each vector
           // Avoids potential problem of different vector lengths in EM and Had
           if ((*it)->layer() == 0) {
             emET[Peak] = (*it)->cpET();
             emError[Peak] = (*it)->errorWord();
           }
           else {
             hadET[Peak] = (*it)->cpET();
             hadError[Peak] = (*it)->errorWord();
           }
           
         } // Loop through TT for this CPMT
        
         if (!zeroSuppress || emET[Peak] > 0 || hadET[Peak] > 0) {
           CPMTower* cpmTower = new CPMTower(phi, eta, emET, emError, hadET, hadError, Peak);
           cpmTowers->push_back(cpmTower);
         }
         
       } // Check vector has non-zero length
     } // end of loop through Sums map
     
   return;
 }

◆ makeCPMTowers() [3/3]

void LVL1::L1CPMTowerTools::makeCPMTowers	(	const DataVector< xAOD::TriggerTower > *	triggerTowers,
		DataVector< xAOD::CPMTower > *	cpmTowers,
		bool	zeroSuppress = `true`
	)

virtual

Fill DataVector of CPMTowers from user-supplied xAOD::TriggerTowers.

Fill DataVector of xAOD::CPMTowers from user-supplied vector of xAOD::TriggerTowers.

Need to get pointers to the TriggerTowers corresponding to each CPMTower
Store these in a vector for each CPMTower
Keep vectors of pointers in a std::map so can easily locate and add more
towers to correct element Right now this is redundant, but when TT becomes a single-layer object
we'll need this bit

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 CPMTower.
Now we need to go through them and add their energies to the appropriate
layer to form the actual CPMTower 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.

Now we can create the tower object

Implements LVL1::IL1CPMTowerTools.

Definition at line 234 of file L1CPMTowerTools.cxx.

                                                                                                                                             {
  
   // Clear collection before filling
   cpmTowers->clear();
   
   std::map< int, xAODTriggerTowerVector_t > Sums;
  
   // Step over all TriggerTowers, and put into map
   xAODTTCollection::const_iterator it ;
   TriggerTowerKey testKey(0.0, 0.0);
  
   for( it = triggerTowers->begin(); it != triggerTowers->end(); ++it ) {
     // Check within CPM tower coverage
     if (fabs((*it)->eta()) > 2.5) continue;
     
     // Find TriggerTowerKey for this TriggerTower
     int key = testKey.ttKey((*it)->phi(),(*it)->eta());
     // Does the map already contain an entry for this CPMTower?
     std::map< int, xAODTriggerTowerVector_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.
       xAODTriggerTowerVector_t vec;
       vec.push_back((*it));
       Sums.insert(std::map< int, xAODTriggerTowerVector_t >::value_type(key,vec));
     }
   } // end of loop over towers
  
     for (std::map< int, xAODTriggerTowerVector_t >::iterator mapIt = Sums.begin();
          mapIt != Sums.end(); ++mapIt) {
  
       // create empty result vectors containing a single element
       // That way if one layer is missing from TT collection we have a 0 entered already
       std::vector<uint8_t> emET(1);
       std::vector<uint8_t> hadET(1);
       std::vector<uint32_t> emError(1);
       std::vector<uint32_t> hadError(1);
       // Both should contain same number of samples, and same peak position
       int Peak = 0;
       
       // Get first TT for this CPMT
       xAODTriggerTowerVector_t::iterator it = (mapIt->second).begin();
       if (it != (mapIt->second).end()) {
         // Get CPMT eta, phi using first tower in vector (either tower in CPMT should do)
         double phi = (*it)->phi();
         double eta = (*it)->eta();
         // now loop through all TT present and add their ET values to CPMT
         for (; it != (mapIt->second).end(); ++it) {
           // Going to play safe and fill just one entry into each vector
           // Avoids potential problem of different vector lengths in EM and Had
           if ((*it)->layer() == 0) {
             emET[Peak] = int((*it)->cpET());
             emError[Peak] = (*it)->errorWord();
           }
           else {
             hadET[Peak] = int((*it)->cpET());
             hadError[Peak] = (*it)->errorWord();
           }
           
         } // Loop through TT for this CPMT
        
         if (!zeroSuppress || emET[Peak] > 0 || hadET[Peak] > 0) {
           xAOD::CPMTower* cpmTower = new xAOD::CPMTower();
           cpmTowers->push_back(cpmTower);
           cpmTower->initialize(eta, phi, emET, hadET, emError, hadError, Peak);
         }
         
       } // Check vector has non-zero length
     } // end of loop through Sums map
     
   return;
 }

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ 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 DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ 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_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 Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

◆ TTCollection

◆ xAODTTCollection

Constructor & Destructor Documentation

◆ L1CPMTowerTools()

◆ ~L1CPMTowerTools()

Member Function Documentation

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

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ finalize()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ makeCPMTowers() [1/3]

◆ makeCPMTowers() [2/3]

◆ makeCPMTowers() [3/3]

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_detStore

◆ m_evtStore

◆ m_varHandleArraysDeclared

◆ m_vhka