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

The aim of this algorithm is the simulation of the LAr analogue Super-Cell sums. More...

#include <LArCalibInject_timeSh.h>

Inheritance diagram for LArCalibInject_timeSh:
Collaboration diagram for LArCalibInject_timeSh:

Public Member Functions

 LArCalibInject_timeSh (const std::string &name, ISvcLocator *pSvcLocator)
 constructor
 ~LArCalibInject_timeSh ()
 destructor
StatusCode initialize ()
 Read ascii files for auxiliary data (puslse shapes, noise, etc...)
StatusCode execute (const EventContext &context) const
 Create LAr object save in TES (2 containers: 1 EM, 1 hadronic)
StatusCode finalize ()
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
virtual unsigned int cardinality () const override
 Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
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.
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

template<class T>
const T * retrieve (const EventContext &context, const SG::ReadCondHandleKey< T > &handleKey) const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

SG::ReadCondHandleKey< ILArfSamplm_fracSKey
 Property: Fraction of Energy Sampled (conditions input).
SG::ReadHandleKey< LArHitEMapm_hitMapKey
 hit map
SG::WriteHandleKey< LArHitEMapm_hitMapOutKey
 output Lar Digits Cell container
SG::ReadCondHandleKey< LArOnOffIdMappingm_cablingKey
 seedless option
const CaloCell_IDm_calocell_id
SG::ReadCondHandleKey< CaloDetDescrManagerm_caloMgrKey
bool m_seedLess
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
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

The aim of this algorithm is the simulation of the LAr analogue Super-Cell sums.

It includes correct allocation of cells to Super-Cell, pulse profile as a function of time, saturation, appropriate noise addition, pile-up.
The resulting signals are an input to the Preprocessor (which in turn performs digitization, filtering, bunch crossing id., noise suppression,...)
Since it needs hits, the simulation only takes "simul" datasets as input, NOT digitized datasets.

Warning
although the output is not digitized
Author
Denis O. Damazio (BNL)

Definition at line 48 of file LArCalibInject_timeSh.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ LArCalibInject_timeSh()

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

constructor

Definition at line 36 of file LArCalibInject_timeSh.cxx.

38 : AthReentrantAlgorithm(name, pSvcLocator), m_calocell_id(nullptr) {
39 //
40 // ........ default values of private data
41 //
42
43 m_seedLess = false;
44 return;
45}
const CaloCell_ID * m_calocell_id

◆ ~LArCalibInject_timeSh()

LArCalibInject_timeSh::~LArCalibInject_timeSh ( )

destructor

CalibInject_timeSh Maker destructor

Definition at line 47 of file LArCalibInject_timeSh.cxx.

47 {
49
50 return;
51}

Member Function Documentation

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::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 75 of file AthCommonReentrantAlgorithm.cxx.

64{
65 return 0;
66}

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::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< Gaudi::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< 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.

◆ evtStore()

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.

◆ execute()

StatusCode LArCalibInject_timeSh::execute ( const EventContext & context) const

Create LAr object save in TES (2 containers: 1 EM, 1 hadronic)

Definition at line 75 of file LArCalibInject_timeSh.cxx.

75 {
76
77 ATH_MSG_DEBUG("Begining of execution");
78
79 //
80 // ....... fill the LArHitEMap
81 //
82
83 /* will be used in the near future for better energy to rec energy matching
84 auto fracS = this->retrieve(context,m_fracSKey);
85 */
86
87 SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{m_caloMgrKey, context};
88 ATH_CHECK(caloMgrHandle.isValid());
89
90 const auto* cabling = this->retrieve(context, m_cablingKey);
91
92 SG::WriteHandle<LArHitEMap> hitEMapOutHandle(m_hitMapOutKey, context);
93 auto hitEMapOut = std::make_unique<LArHitEMap>(cabling, m_calocell_id,
94 *caloMgrHandle, false);
95
96 //
97 // ... initialise vectors for sums of energy in each TT
98 //
99
100 int it = 0;
101 int it_end(200);
102
103 // m_nSamples=5;
104
105 // Copy what exist
106 if (!m_hitMapKey.empty()) {
107 SG::ReadHandle<LArHitEMap> hitmap(m_hitMapKey, context);
108 const auto* hitmapPtr = hitmap.cptr();
109 it_end = hitmapPtr->GetNbCells();
110
111 for (; it != it_end; ++it) {
112 const LArHitList& hitlist = hitmapPtr->GetCell(it);
113 const std::vector<std::pair<float, float> >& timeE = hitlist.getData();
114 if (!timeE.empty()) {
115 const IdentifierHash idhash(it);
116 for (size_t i = 0; i < timeE.size(); i++) {
117 std::pair<float, float> energy_time = timeE[i];
118 hitEMapOut->AddEnergy(idhash, energy_time.first, energy_time.second);
119 }
120 } // end of if timeE.size()
121 } // it end
122 } // end of empty Key check
123 // Now, become creative
124 it = 0;
125 unsigned seed = std::chrono::system_clock::now().time_since_epoch().count();
126 std::default_random_engine generator(seed);
127 std::uniform_real_distribution<double> distribution(-12.5, 12.5);
128 it_end = 480;
129 for (; it != it_end; it += 16) {
130 const IdentifierHash idhash(it);
131 float time_inject = (float)distribution(generator);
132 hitEMapOut->AddEnergy(idhash, 10500.0, time_inject);
133 }
134 ATH_CHECK(hitEMapOutHandle.record(std::move(hitEMapOut)));
135
136 return StatusCode::SUCCESS;
137}
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_DEBUG(x)
SG::ReadCondHandleKey< LArOnOffIdMapping > m_cablingKey
seedless option
SG::WriteHandleKey< LArHitEMap > m_hitMapOutKey
output Lar Digits Cell container
SG::ReadCondHandleKey< CaloDetDescrManager > m_caloMgrKey
const T * retrieve(const EventContext &context, const SG::ReadCondHandleKey< T > &handleKey) const
SG::ReadHandleKey< LArHitEMap > m_hitMapKey
hit map
const LARLIST & getData() const
Definition LArHitList.h:25
generator
Configure Herwig7 These are the commands corresponding to what would go into the regular Herwig infil...

◆ 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 & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

90{
91 // If we didn't find any symlinks to add, just return the collection
92 // from the base class. Otherwise, return the extended collection.
93 if (!m_extendedExtraObjects.empty()) {
95 }
97}
An algorithm that can be simultaneously executed in multiple threads.

◆ filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx) const
inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

96 {
97 return execState( ctx ).filterPassed();
98 }
virtual bool filterPassed(const EventContext &ctx) const

◆ finalize()

StatusCode LArCalibInject_timeSh::finalize ( )

Definition at line 139 of file LArCalibInject_timeSh.cxx.

139 {
140
141 ATH_MSG_INFO(" LArCalibInject_timeSh finalize completed successfully");
142
143 return StatusCode::SUCCESS;
144}
#define ATH_MSG_INFO(x)

◆ initialize()

StatusCode LArCalibInject_timeSh::initialize ( )

Read ascii files for auxiliary data (puslse shapes, noise, etc...)

Definition at line 53 of file LArCalibInject_timeSh.cxx.

53 {
54 //
55 // ......... declaration
56 //
57
58 ATH_MSG_INFO("***********************************************");
59 ATH_MSG_INFO("* Steering options for LArCalibInject_timeSh algorithm *");
60 ATH_MSG_INFO("***********************************************");
61
62 CHECK(m_fracSKey.initialize());
63 CHECK(m_cablingKey.initialize());
64 CHECK(m_hitMapKey.initialize(SG::AllowEmpty));
65 CHECK(m_hitMapOutKey.initialize());
66 ATH_CHECK(detStore()->retrieve(m_calocell_id, "CaloCell_ID"));
67 ATH_CHECK(m_caloMgrKey.initialize());
68
69 // Incident Service:
70 ATH_MSG_DEBUG("Initialization completed successfully");
71
72 return StatusCode::SUCCESS;
73}
#define CHECK(...)
Evaluate an expression and check for errors.
const ServiceHandle< StoreGateSvc > & detStore() const
SG::ReadCondHandleKey< ILArfSampl > m_fracSKey
Property: Fraction of Energy Sampled (conditions input).

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

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msgLvl()

bool AthCommonMsg< Gaudi::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< 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.

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< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ retrieve()

template<class T>
const T * LArCalibInject_timeSh::retrieve ( const EventContext & context,
const SG::ReadCondHandleKey< T > & handleKey ) const
inlineprivate

Definition at line 76 of file LArCalibInject_timeSh.h.

77 {
78 SG::ReadCondHandle<T> handle(handleKey, context);
79 if (not handle.isValid()) {
80 ATH_MSG_ERROR("could not retrieve : " << handle.key());
81 return nullptr;
82 } else
83 return handle.cptr();
84 }
#define ATH_MSG_ERROR(x)

◆ setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state,
const EventContext & ctx ) const
inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

100 {
102 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const

◆ sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::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 85 of file AthCommonReentrantAlgorithm.cxx.

77{
78 return BaseAlg::sysExecute (ctx);
79}

◆ sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::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 HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

107 {
109
110 if (sc.isFailure()) {
111 return sc;
112 }
113
114 ServiceHandle<ICondSvc> cs("CondSvc",name());
115 for (auto h : outputHandles()) {
116 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
117 // do this inside the loop so we don't create the CondSvc until needed
118 if ( cs.retrieve().isFailure() ) {
119 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
120 return StatusCode::SUCCESS;
121 }
122 if (cs->regHandle(this,*h).isFailure()) {
124 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
125 << " with CondSvc");
126 }
127 }
128 }
129 return sc;
130}
#define ATH_MSG_WARNING(x)
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override

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

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 }

Member Data Documentation

◆ m_cablingKey

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

seedless option

Definition at line 98 of file LArCalibInject_timeSh.h.

98 {
99 this, "CablingKey", "LArOnOffIdMap",
100 "SG Key of LArOnOffIdMapping object"};

◆ m_calocell_id

const CaloCell_ID* LArCalibInject_timeSh::m_calocell_id
private

Definition at line 102 of file LArCalibInject_timeSh.h.

◆ m_caloMgrKey

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

Definition at line 103 of file LArCalibInject_timeSh.h.

103 {
104 this, "CaloDetDescrManager", "CaloDetDescrManager",
105 "SG Key for CaloDetDescrManager in the Condition Store"};

◆ 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 AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

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

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_fracSKey

SG::ReadCondHandleKey<ILArfSampl> LArCalibInject_timeSh::m_fracSKey
private
Initial value:
{
this, "FracSKey", "LArfSampl", "SG Key of fSamplS conditions object"}

Property: Fraction of Energy Sampled (conditions input).

Definition at line 72 of file LArCalibInject_timeSh.h.

72 {
73 this, "FracSKey", "LArfSampl", "SG Key of fSamplS conditions object"};

◆ m_hitMapKey

SG::ReadHandleKey<LArHitEMap> LArCalibInject_timeSh::m_hitMapKey
private
Initial value:
{this, "LArHitEMapKey", "LArHitEMap",
"Input HitMap from cell hits"}

hit map

Definition at line 91 of file LArCalibInject_timeSh.h.

91 {this, "LArHitEMapKey", "LArHitEMap",
92 "Input HitMap from cell hits"};

◆ m_hitMapOutKey

SG::WriteHandleKey<LArHitEMap> LArCalibInject_timeSh::m_hitMapOutKey
private
Initial value:
{
this, "LArHitEMapOutKey", "LArHitEMapOut", "Output LArDigit container"}

output Lar Digits Cell container

Definition at line 94 of file LArCalibInject_timeSh.h.

94 {
95 this, "LArHitEMapOutKey", "LArHitEMapOut", "Output LArDigit container"};

◆ m_seedLess

bool LArCalibInject_timeSh::m_seedLess
private

Definition at line 107 of file LArCalibInject_timeSh.h.

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.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: