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

#include <LArGeoWeightsFill.h>

Inheritance diagram for LArGeoWeightsFill:
Collaboration diagram for LArGeoWeightsFill:

Public Member Functions

 LArGeoWeightsFill (const std::string &name, ISvcLocator *pSvcLocator)
 ~LArGeoWeightsFill ()
virtual StatusCode initialize () override
virtual StatusCode execute () override
virtual StatusCode stop () 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

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_cablingKey {this, "OnOffMap", "LArOnOffIdMap", "SG key for mapping object"}
SG::ReadCondHandleKey< CaloDetDescrManagerm_caloMgrKey {this,"CaloDetDescrManager","CaloDetDescrManager","SG Key for CaloDetDescrManager in the Condition Store" }
const LArOnlineIDm_onlineID
std::string m_nameOfSet
std::string m_key
bool m_fill
bool m_dump
std::string m_outFileName
ToolHandle< CaloTriggerTowerServicem_ttService
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 17 of file LArGeoWeightsFill.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

◆ LArGeoWeightsFill()

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

Definition at line 18 of file LArGeoWeightsFill.cxx.

18 :
19 AthAlgorithm(name,pSvcLocator),
20 m_onlineID(0),
21 m_ttService("CaloTriggerTowerService")
22{
23 declareProperty("Key",m_key="GeoWeights");
24 declareProperty("Dump",m_dump=false);
25 declareProperty("OutFile",m_outFileName="out.txt");
26 //declareProperty("InFile",m_inFileName="");
27
28 declareProperty("Fill",m_fill=true);
29}
AthAlgorithm()
Default constructor:
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
ToolHandle< CaloTriggerTowerService > m_ttService
const LArOnlineID * m_onlineID
std::string m_outFileName

◆ ~LArGeoWeightsFill()

LArGeoWeightsFill::~LArGeoWeightsFill ( )

Definition at line 31 of file LArGeoWeightsFill.cxx.

31{}

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 }

◆ 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()

virtual StatusCode LArGeoWeightsFill::execute ( )
inlineoverridevirtual

Definition at line 22 of file LArGeoWeightsFill.h.

22{return StatusCode::SUCCESS;}

◆ 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 LArGeoWeightsFill::initialize ( )
overridevirtual

Definition at line 34 of file LArGeoWeightsFill.cxx.

34 {
35
36 ATH_MSG_DEBUG ( "start initialize()" );
37 ATH_CHECK( detStore()->retrieve(m_onlineID,"LArOnlineID") );
38 ATH_CHECK( m_cablingKey.initialize() );
39 ATH_CHECK( m_caloMgrKey.initialize() );
40 ATH_CHECK( m_ttService.retrieve() );
41 return StatusCode::SUCCESS;
42}
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_DEBUG(x)
const ServiceHandle< StoreGateSvc > & detStore() const
SG::ReadCondHandleKey< CaloDetDescrManager > m_caloMgrKey
SG::ReadCondHandleKey< LArOnOffIdMapping > m_cablingKey
retrieve(aClass, aKey=None)
Definition PyKernel.py:110

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

◆ stop()

StatusCode LArGeoWeightsFill::stop ( )
overridevirtual

Definition at line 44 of file LArGeoWeightsFill.cxx.

44 {
45
46 ATH_MSG_DEBUG ( "start stop()" );
47
48 const unsigned hashMax=m_onlineID->channelHashMax();
49
50 if (m_fill) {
51
52 ATH_MSG_INFO ( "Filling database" );
53
54 coral::AttributeListSpecification* spec = new coral::AttributeListSpecification();
55 //cool::RecordSpecification* spec = new cool::RecordSpecification();
56 spec->extend("costheta", "blob");
57 spec->extend("sinthetacosphi", "blob");
58 spec->extend("sinthetasinphi", "blob");
59 spec->extend("offlineTTid","blob");
60
61 AthenaAttributeList* attr=new AthenaAttributeList(*spec);
62 coral::Blob& costhetaBlob = (*attr)["costheta"].data<coral::Blob>();
63 coral::Blob& sinthetacosphiBlob = (*attr)["sinthetacosphi"].data<coral::Blob>();
64 coral::Blob& sinthetasinphiBlob = (*attr)["sinthetasinphi"].data<coral::Blob>();
65 coral::Blob& offlineTTidBlob = (*attr)["offlineTTid"].data<coral::Blob>();
66
67
68 costhetaBlob.resize(hashMax*sizeof(float));
69 sinthetacosphiBlob.resize(hashMax*sizeof(float));
70 sinthetasinphiBlob.resize(hashMax*sizeof(float));
71 offlineTTidBlob.resize(hashMax*sizeof(uint32_t));
72
73 float* pcostheta=static_cast<float*>(costhetaBlob.startingAddress());
74 float* psinthetacosphi=static_cast<float*>(sinthetacosphiBlob.startingAddress());
75 float* psinthetasinphi=static_cast<float*>(sinthetasinphiBlob.startingAddress());
76 uint32_t* pofflineTTid=static_cast<uint32_t*>(offlineTTidBlob.startingAddress());
77
78 ATH_CHECK( detStore()->record(attr,m_key) );
79
80 SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{m_caloMgrKey};
81 ATH_CHECK(caloMgrHandle.isValid());
82 const CaloDetDescrManager *theCaloDDM = *caloMgrHandle;
83 ATH_MSG_INFO ( "theCaloDDM retrieved" );
84 SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKey};
85 const LArOnOffIdMapping* cabling{*cablingHdl};
86 if(!cabling) {
87 ATH_MSG_ERROR("Do not have mapping object " << m_cablingKey.key());
88 return StatusCode::FAILURE;
89 }
90
91 for (unsigned hs=0;hs<hashMax;++hs) {
92 const HWIdentifier chid=m_onlineID->channel_Id(hs);
93
94 if(!cabling->isOnlineConnected(chid)){
95 ATH_MSG_DEBUG ( "cell chid: " << chid.get_compact() << " not connected channel, skip " );
96 //Set values for disconnected cells to 0
97 pcostheta[hs]=0.0;
98 psinthetacosphi[hs]=0.0;
99 psinthetasinphi[hs]=0.0;
100 pofflineTTid[hs]=0;
101 continue;
102 }
103 const Identifier id=cabling->cnvToIdentifier(chid);
104 const CaloDetDescrElement *caloDDE = theCaloDDM->get_element(id);
105 if(!caloDDE){
106 ATH_MSG_ERROR ( "Failed to return CaloDetDescrElement" );
107 return StatusCode::FAILURE;
108 }
109 ATH_MSG_DEBUG ( "hash, eta, phi: " << caloDDE->calo_hash() << ", " << caloDDE->eta() << ", " << caloDDE->phi() );
110
111 const float eta = caloDDE->eta();
112 const float phi = caloDDE->phi();
113 const float v = 1./cosh(eta); // sintheta
114 const float costheta=tanh(eta);
115 const float sinthetacosphi=v*cos(phi);
116 const float sinthetasinphi=v*sin(phi);
117
118 pcostheta[hs]=costheta;
119 psinthetacosphi[hs]=sinthetacosphi;
120 psinthetasinphi[hs]=sinthetasinphi;
121 pofflineTTid[hs]=m_ttService->whichTTID(id).get_identifier32().get_compact();
122 } // end loop over hash ids
123 }//end if FILL
124
125
126 if (m_dump) {
127
128 const AthenaAttributeList* attr=0;
129
130 ATH_CHECK( detStore()->retrieve(attr,m_key) );
131
132 const coral::Blob& costhetaBlob = (*attr)["costheta"].data<coral::Blob>();
133 const coral::Blob& sinthetacosphiBlob = (*attr)["sinthetacosphi"].data<coral::Blob>();
134 const coral::Blob& sinthetasinphiBlob = (*attr)["sinthetasinphi"].data<coral::Blob>();
135 const coral::Blob& offlineTTidBlob = (*attr)["offlineTTid"].data<coral::Blob>();
136
137 const float* pcostheta=static_cast<const float*>(costhetaBlob.startingAddress());
138 const float* psinthetacosphi=static_cast<const float*>(sinthetacosphiBlob.startingAddress());
139 const float* psinthetasinphi=static_cast<const float*>(sinthetasinphiBlob.startingAddress());
140 const uint32_t* pofflineTTid=static_cast<const uint32_t*>(offlineTTidBlob.startingAddress());
141
142
143 std::ostream *out = &(std::cout);
144 std::ofstream outfile;
145 if (m_outFileName.size()) {
146 outfile.open(m_outFileName.c_str(),std::ios::out);
147 if (outfile.is_open()) {
148 ATH_MSG_INFO ( "Writing to file " << m_outFileName );
149 out = &outfile;
150 }
151 else
152 ATH_MSG_ERROR ( "Failed to open file " << m_outFileName );
153 }
154
155 for (unsigned hs=0;hs<hashMax;++hs) {
156 const HWIdentifier chid=m_onlineID->channel_Id(hs);
157 (*out) << std::hex << chid.get_identifier32().get_compact() << std::dec << " "
158 << pcostheta[hs] << " " << psinthetacosphi[hs] << " " << psinthetasinphi[hs]
159 << " 0x"<< std::hex << pofflineTTid[hs] << std::dec << std::endl;
160 }
161
162
163 if (outfile.is_open())
164 outfile.close();
165 }// end if m_dump
166
167 return StatusCode::SUCCESS;
168}
Scalar eta() const
pseudorapidity method
Scalar phi() const
phi method
#define ATH_MSG_ERROR(x)
#define ATH_MSG_INFO(x)
const CaloDetDescrElement * get_element(const Identifier &cellId) const
get element by its identifier
value_type get_compact() const
Get the compact id.
value_type get_compact() const
Get the compact id.
Identifier32 get_identifier32() const
Get the 32-bit version Identifier, will be invalid if >32 bits needed.
setEventNumber uint32_t

◆ 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> LArGeoWeightsFill::m_cablingKey {this, "OnOffMap", "LArOnOffIdMap", "SG key for mapping object"}
private

Definition at line 26 of file LArGeoWeightsFill.h.

26{this, "OnOffMap", "LArOnOffIdMap", "SG key for mapping object"};

◆ m_caloMgrKey

SG::ReadCondHandleKey<CaloDetDescrManager> LArGeoWeightsFill::m_caloMgrKey {this,"CaloDetDescrManager","CaloDetDescrManager","SG Key for CaloDetDescrManager in the Condition Store" }
private

Definition at line 27 of file LArGeoWeightsFill.h.

27{this,"CaloDetDescrManager","CaloDetDescrManager","SG Key for CaloDetDescrManager in the Condition Store" };

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

bool LArGeoWeightsFill::m_dump
private

Definition at line 34 of file LArGeoWeightsFill.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_fill

bool LArGeoWeightsFill::m_fill
private

Definition at line 34 of file LArGeoWeightsFill.h.

◆ m_key

std::string LArGeoWeightsFill::m_key
private

Definition at line 32 of file LArGeoWeightsFill.h.

◆ m_nameOfSet

std::string LArGeoWeightsFill::m_nameOfSet
private

Definition at line 31 of file LArGeoWeightsFill.h.

◆ m_onlineID

const LArOnlineID* LArGeoWeightsFill::m_onlineID
private

Definition at line 29 of file LArGeoWeightsFill.h.

◆ m_outFileName

std::string LArGeoWeightsFill::m_outFileName
private

Definition at line 35 of file LArGeoWeightsFill.h.

◆ m_ttService

ToolHandle< CaloTriggerTowerService > LArGeoWeightsFill::m_ttService
private

Definition at line 37 of file LArGeoWeightsFill.h.

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