ATLAS Offline Software
Loading...
Searching...
No Matches
LArShapeToSCShape Class Reference

#include <LArShapeToSCShape.h>

Inheritance diagram for LArShapeToSCShape:
Collaboration diagram for LArShapeToSCShape:

Public Member Functions

 LArShapeToSCShape (const std::string &name, ISvcLocator *pSvcLocator)
 Constructor with parameters:
virtual ~LArShapeToSCShape ()
 Destructor:
virtual StatusCode initialize () override
virtual StatusCode execute () override
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Protected Member Functions

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

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

 LArShapeToSCShape ()
 Default constructor:
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

SG::ReadCondHandleKey< LArOnOffIdMappingm_cablingKeySC {this,"SCCablingKey","LArOnOffIdMapSC","SG Key of SC LArOnOffIdMapping object"}
SG::ReadCondHandleKey< LArOnOffIdMappingm_cablingKey {this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"}
ToolHandle< ICaloSuperCellIDToolm_scidTool {this, "CaloSuperCellIDTool", "CaloSuperCellIDTool"}
DataObjIDColl m_extendedExtraObjects
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default)
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

Definition at line 18 of file LArShapeToSCShape.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ LArShapeToSCShape() [1/2]

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

Constructor with parameters:

Definition at line 30 of file LArShapeToSCShape.cxx.

31 :
32 ::AthAlgorithm( name, pSvcLocator )
33{
34}
AthAlgorithm()
Default constructor:

◆ ~LArShapeToSCShape()

LArShapeToSCShape::~LArShapeToSCShape ( )
virtual

Destructor:

Definition at line 38 of file LArShapeToSCShape.cxx.

39{}

◆ LArShapeToSCShape() [2/2]

LArShapeToSCShape::LArShapeToSCShape ( )
private

Default constructor:

Member Function Documentation

◆ declareGaudiProperty()

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

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

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

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

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode LArShapeToSCShape::execute ( )
overridevirtual

Definition at line 57 of file LArShapeToSCShape.cxx.

58{
59 ATH_MSG_DEBUG ("Executing " << name() << "...");
60
61
62 //Retrieve shape for regular cells
63 const ILArShape* cellShape = nullptr;
64 CHECK(detStore()->retrieve(cellShape,"LArShape"));
65
66 //Retrieve SuperCell id
67 const CaloCell_SuperCell_ID* caloSCID = nullptr;
68 CHECK(detStore()->retrieve(caloSCID));
69
70 //Retrieve SuperCell online id
71 const LArOnline_SuperCellID* onlSCID = nullptr;
72 CHECK(detStore()->retrieve(onlSCID));
73
74 const EventContext& ctx = Gaudi::Hive::currentContext();
75
76 SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKey, ctx};
77 const LArOnOffIdMapping* cabling{*cablingHdl};
78 if(!cabling) {
79 ATH_MSG_ERROR( "Do not have cabling mapping from key " << m_cablingKey.key() );
80 return StatusCode::FAILURE;
81 }
82
83 SG::ReadCondHandle<LArOnOffIdMapping> cablingHdlSC{m_cablingKeySC, ctx};
84 const LArOnOffIdMapping* cablingSC{*cablingHdlSC};
85 if(!cablingSC) {
86 ATH_MSG_ERROR( "Do not have cabling mapping from key " << m_cablingKeySC.key() );
87 return StatusCode::FAILURE;
88 }
89
90 const unsigned hashMax=onlSCID->channelHashMax();
91 ATH_MSG_DEBUG("SuperCell online hash max: " << hashMax);
92
93 //Set up AttributeListCollection
94 coral::AttributeListSpecification* spec = new coral::AttributeListSpecification();
95 spec->extend("Shape", "blob");
96 spec->extend("ShapeDer", "blob");
97 spec->extend("TimeOffset","blob");
98 spec->extend<unsigned>("nSamples");
99 spec->extend<unsigned>("version");
100
101 unsigned nSamples=32; //Assumption....
102
103 const unsigned blobSize=hashMax*nSamples;
104 ATH_MSG_DEBUG("Shape blob has size " << blobSize);
105
106 coral::AttributeList attrList(*spec);
107 attrList["version"].setValue(0U);
108 coral::Blob& shapeBlob=attrList["Shape"].data<coral::Blob>();
109 coral::Blob& shapeDerBlob=attrList["ShapeDer"].data<coral::Blob>();
110 coral::Blob& toBlob=attrList["TimeOffset"].data<coral::Blob>();
111
112 spec->release();
113 // cppcheck-suppress memleak
114 spec = nullptr;
115
116 attrList["nSamples"].setValue(nSamples);
117 shapeBlob.resize(blobSize*sizeof(float));
118 shapeDerBlob.resize(blobSize*sizeof(float));
119 toBlob.resize(hashMax*sizeof(float));
120
121 float *pShape=static_cast<float*>(shapeBlob.startingAddress());
122 float *pShapeDer=static_cast<float*>(shapeDerBlob.startingAddress());
123 float *pTimeOffset=static_cast<float*>(toBlob.startingAddress());
124
125
126 //Initialize blobs to ERRORCODE
127 for (unsigned i=0;i<blobSize;++i) {
129 pShapeDer[i]=LArElecCalib::ERRORCODE;
130 }
131
132 //Initialize time-offset blob to 0.0
133 for (unsigned i=0;i<hashMax;++i) {
134 pTimeOffset[i]=0.0;
135 }
136
137 unsigned nTileIds=0;
138 unsigned nTotalIds=0;
139 auto scIt=caloSCID->cell_begin();
140 auto scIt_e=caloSCID->cell_end();
141 for (;scIt!=scIt_e;++scIt) {
142 ++nTotalIds;
143 const Identifier scId=*scIt;
144 if (caloSCID->is_tile(scId)) {
145 ++nTileIds;
146 continue;
147 }
148 //const IdentifierHash scHash=caloSCID->calo_cell_hash(scId);
149 const std::vector<Identifier> &cellIds=m_scidTool->superCellToOfflineID(scId);
150 if (cellIds.empty()) {
151 ATH_MSG_ERROR("Got empty vector of cell ids for super cell id 0x"
152 << std::hex << scId.get_compact());
153 return StatusCode::FAILURE;
154 }
155 const Identifier cellId=cellIds[0];
156
157 //Get online hash (via online identifier):
158 const HWIdentifier scOnlId=cablingSC->createSignalChannelID(scId);
159 const IdentifierHash scOnlHash=onlSCID->channel_Hash(scOnlId);
160
161 auto shape=cellShape->Shape(cabling->createSignalChannelID(cellId),0);
162 auto shapeDer=cellShape->ShapeDer(cabling->createSignalChannelID(cellId),0);
163 //consistency check
164 if (nSamples!=shape.size() || nSamples!=shapeDer.size()) {
165 ATH_MSG_ERROR("Inconsistent number of samples!");
166 ATH_MSG_ERROR("Expected" << nSamples << ", Shape:" << shape.size() << ", ShapeDer:" << shapeDer.size());
167 return StatusCode::FAILURE;
168 }
169 for (unsigned i=0;i<nSamples;++i) {
170 pShape[scOnlHash*nSamples+i]=shape[i];
171 pShapeDer[scOnlHash*nSamples+i]=shapeDer[i];
172
173 }
174
175 }//end loop over super-cells
176
177 //Add to collection
178 CondAttrListCollection* coll=new CondAttrListCollection(true);
179 CHECK(detStore()->record(coll,"/LAR/ElecCalibSC/Shape"));
180 coll->add(0,attrList);
181
182 ATH_MSG_INFO("Total number of SuperCells:" << nTotalIds << ", Tile:" << nTileIds
183 << ", LAr: " << nTotalIds-nTileIds);
184 LArShapeSC* scShape=new LArShapeSC(coll);
185 CHECK(detStore()->record(scShape,"ShapeSC"));
186
187 return StatusCode::SUCCESS;
188}
#define ATH_MSG_ERROR(x)
#define ATH_MSG_INFO(x)
#define ATH_MSG_DEBUG(x)
#define CHECK(...)
Evaluate an expression and check for errors.
const ServiceHandle< StoreGateSvc > & detStore() const
id_iterator cell_begin() const
begin iterator over full set of Identifiers (LAr + Tiles)
id_iterator cell_end() const
end iterator over full set of Identifiers (LAr + Tiles)
bool is_tile(const Identifier id) const
test if the id belongs to the Tiles
bool add(ChanNum chanNum, const AttributeList &attributeList)
Adding in chan/attrList pairs.
virtual ShapeRef_t Shape(const HWIdentifier &id, int gain, int tbin=0, int mode=0) const =0
virtual ShapeRef_t ShapeDer(const HWIdentifier &id, int gain, int tbin=0, int mode=0) const =0
value_type get_compact() const
Get the compact id.
HWIdentifier createSignalChannelID(const Identifier &id) const
create a HWIdentifier from an Identifier (not inline)
IdentifierHash channel_Hash(HWIdentifier channelId) const
Create channel_hash from channel_Id.
size_type channelHashMax() const
Define channel hash tables max size.
SG::ReadCondHandleKey< LArOnOffIdMapping > m_cablingKeySC
SG::ReadCondHandleKey< LArOnOffIdMapping > m_cablingKey
ToolHandle< ICaloSuperCellIDTool > m_scidTool
retrieve(aClass, aKey=None)
Definition PyKernel.py:110

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< 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 & AthAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 50 of file AthAlgorithm.cxx.

51{
52 // If we didn't find any symlinks to add, just return the collection
53 // from the base class. Otherwise, return the extended collection.
54 if (!m_extendedExtraObjects.empty()) {
56 }
57 return Algorithm::extraOutputDeps();
58}
DataObjIDColl m_extendedExtraObjects

◆ initialize()

StatusCode LArShapeToSCShape::initialize ( )
overridevirtual

Definition at line 43 of file LArShapeToSCShape.cxx.

44{
45
46 ATH_CHECK(m_scidTool.retrieve());
47 ATH_CHECK(m_cablingKeySC.initialize());
48 ATH_CHECK(m_cablingKey.initialize());
49 ATH_MSG_DEBUG ("Initializing " << name() << "...");
50
51
52
53 return StatusCode::SUCCESS;
54}
#define ATH_CHECK
Evaluate an expression and check for errors.

◆ inputHandles()

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

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

24 {
25 return this->msgStream();
26 }

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< 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< Algorithm > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
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)

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ sysInitialize()

StatusCode AthAlgorithm::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< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, AthHistogramAlgorithm, and PyAthena::Alg.

Definition at line 66 of file AthAlgorithm.cxx.

66 {
68
69 if (sc.isFailure()) {
70 return sc;
71 }
72 ServiceHandle<ICondSvc> cs("CondSvc",name());
73 for (auto h : outputHandles()) {
74 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
75 // do this inside the loop so we don't create the CondSvc until needed
76 if ( cs.retrieve().isFailure() ) {
77 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
78 return StatusCode::SUCCESS;
79 }
80 if (cs->regHandle(this,*h).isFailure()) {
81 sc = StatusCode::FAILURE;
82 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
83 << " with CondSvc");
84 }
85 }
86 }
87 return sc;
88}
#define ATH_MSG_WARNING(x)
static Double_t sc
virtual StatusCode sysInitialize() override
Override sysInitialize.
AthCommonDataStore(const std::string &name, T... args)
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
::StatusCode StatusCode
StatusCode definition for legacy code.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< 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< Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }
std::vector< SG::VarHandleKeyArray * > m_vhka

Member Data Documentation

◆ m_cablingKey

SG::ReadCondHandleKey<LArOnOffIdMapping> LArShapeToSCShape::m_cablingKey {this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"}
private

Definition at line 43 of file LArShapeToSCShape.h.

43{this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"};

◆ m_cablingKeySC

SG::ReadCondHandleKey<LArOnOffIdMapping> LArShapeToSCShape::m_cablingKeySC {this,"SCCablingKey","LArOnOffIdMapSC","SG Key of SC LArOnOffIdMapping object"}
private

Definition at line 42 of file LArShapeToSCShape.h.

42{this,"SCCablingKey","LArOnOffIdMapSC","SG Key of SC LArOnOffIdMapping object"};

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< 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< Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_scidTool

ToolHandle<ICaloSuperCellIDTool> LArShapeToSCShape::m_scidTool {this, "CaloSuperCellIDTool", "CaloSuperCellIDTool"}
private

Definition at line 47 of file LArShapeToSCShape.h.

47{this, "CaloSuperCellIDTool", "CaloSuperCellIDTool"};

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.


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