 |
ATLAS Offline Software
|
#include <LArRawChannelBuilderSCAlg.h>
|
StatusCode | initialize () override |
|
StatusCode | execute (const EventContext &ctx) const override |
|
StatusCode | finalize () override |
|
virtual StatusCode | sysInitialize () override |
| Override sysInitialize. More...
|
|
virtual bool | isClonable () const override |
| Specify if the algorithm is clonable. More...
|
|
virtual unsigned int | cardinality () const override |
| Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...
|
|
virtual StatusCode | sysExecute (const EventContext &ctx) override |
| Execute an algorithm. More...
|
|
virtual const DataObjIDColl & | extraOutputDeps () const override |
| Return the list of extra output dependencies. More...
|
|
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 . 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 | 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 |
|
|
SG::ReadHandleKey< LArDigitContainer > | m_digitKey |
|
SG::WriteHandleKey< CaloCellContainer > | m_cellKey |
|
SG::ReadCondHandleKey< ILArPedestal > | m_pedestalKey {this,"PedestalKey","LArPedestalSC","SG Key of Pedestal conditions object"} |
|
SG::ReadCondHandleKey< LArADC2MeV > | m_adc2MeVKey {this,"ADC2MeVKey","LArADC2MeVSC","SG Key of ADC2MeV conditions object"} |
|
SG::ReadCondHandleKey< ILArOFC > | m_ofcKey {this,"OFCKey","LArOFCSC","SG Key of OFC conditions object"} |
|
SG::ReadCondHandleKey< ILArShape > | m_shapeKey {this,"ShapeKey","LArShapeSC","SG Key of Shape conditions object"} |
|
SG::ReadCondHandleKey< LArOnOffIdMapping > | m_cablingKey {this,"CablingKey","LArOnOffIdMapSC","SG Key of LArOnOffIdMapping object"} |
|
SG::ReadCondHandleKey< LArBadChannelCont > | m_bcContKey {this, "LArBadChannelKey", "LArBadChannelSC", "Key of the LArBadChannelCont SC" } |
| Bad Channel masking for Super-Cells. More...
|
|
Gaudi::Property< float > | m_eCutFortQ {this,"ECutFortQ",256.0,"Time and Quality will be computed only for channels with E above this value"} |
|
Gaudi::Property< bool > | m_absECutFortQ {this,"absECut",true,"Cut on fabs(E) for Q and t computation"} |
|
Gaudi::Property< bool > | m_useShapeDer {this,"useShapeDer",true,"Use shape derivative in Q-factor computation"} |
|
Gaudi::Property< int > | m_firstSample {this,"firstSample",0,"first of the 32 sampels of the MC shape to be used"} |
|
const LArOnline_SuperCellID * | m_onlineId = nullptr |
|
SG::ReadCondHandleKey< CaloSuperCellDetDescrManager > | m_caloSuperCellMgrKey |
|
DataObjIDColl | m_extendedExtraObjects |
| Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks. More...
|
|
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 |
|
Definition at line 31 of file LArRawChannelBuilderSCAlg.h.
◆ StoreGateSvc_t
◆ cardinality()
Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
Override this to return 0 for reentrant algorithms.
Override this to return 0 for reentrant algorithms.
Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
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.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
◆ declareProperty() [2/6]
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.
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
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.
◆ declareProperty() [5/6]
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.
◆ declareProperty() [6/6]
◆ detStore()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode LArRawChannelBuilderSCAlg::execute |
( |
const EventContext & |
ctx | ) |
const |
|
override |
Definition at line 36 of file LArRawChannelBuilderSCAlg.cxx.
42 auto outputContainerCellPtr = std::make_unique<CaloCellContainer>(
SG::VIEW_ELEMENTS);
46 if (dataPool.allocated()==0){
47 dataPool.reserve (hash_max);
49 outputContainerCellPtr->reserve( hash_max );
73 badchannel = *larBadChan;
80 if (!(*cabling)->isOnlineConnected(
id))
continue;
84 const std::vector<short>& samples=
digit->samples();
85 const size_t nSamples=samples.size();
96 return StatusCode::FAILURE;
101 return StatusCode::FAILURE;
105 std::vector<float> samp_no_ped(
nSamples,0.0);
107 samp_no_ped[
i]=samples[
i]-
p;
114 bool passBCIDmax=
false;
116 size_t nOFC=ofca.size();
119 <<
". Found " <<
nSamples <<
", expect at least " << nOFC+2 <<
".");
122 for (
size_t i=0;
i<nOFC;++
i) {
123 A+=
static_cast<double>(samp_no_ped[
i+1])*ofca[
i];
127 for (
size_t i=0;
i<nOFC;++
i) {
128 Abefore+=
static_cast<double>(samp_no_ped[
i])*ofca[
i];
132 for (
size_t i=0;
i<nOFC;++
i) {
133 Aafter+=
static_cast<double>(samp_no_ped[
i+2])*ofca[
i];
136 if ( (
A>Abefore) && (
A>Aafter) ) passBCIDmax=
true;
140 const float E=adc2mev[0]+
A*adc2mev[1];
153 for (
size_t i=0;
i<nOFC;++
i) {
154 At+=
static_cast<double>(samp_no_ped[
i+1])*ofcb[
i];
157 tau=(std::fabs(
A)>0.1) ? At/
A : 0.0;
158 const auto& fullShape=shapes->
Shape(
id,
gain);
171 <<
" gain " <<
gain);
174 nOFC=fullShape.size()-firstSample;
177 const float* shape=&*fullShape.begin()+firstSample;
185 <<
" gain " <<
gain <<
". Will not use shape derivative.");
189 const float* shapeDer=&*fullshapeDer.begin()+firstSample;
190 for (
size_t i=0;
i<nOFC;++
i) {
191 q +=
std::pow((
A*(shape[
i]-tau*shapeDer[
i])-(samp_no_ped[
i+1])),2);
197 for (
size_t i=0;
i<nOFC;++
i) {
202 int iqua =
static_cast<int>(
q);
203 if (iqua > 0xFFFF) iqua=0xFFFF;
204 iquaShort =
static_cast<uint16_t>(iqua & 0xFFFF);
217 float et =
ss->et()*1
e-3;
231 ss->setProvenance(prov);
233 ss->setQuality(iquaShort);
234 outputContainerCellPtr->push_back(
ss);
239 ATH_CHECK(outputContainer.record(std::move(outputContainerCellPtr) ) );
241 return StatusCode::SUCCESS;
◆ extraDeps_update_handler()
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()
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.
◆ filterPassed()
◆ finalize()
StatusCode LArRawChannelBuilderSCAlg::finalize |
( |
| ) |
|
|
override |
◆ initialize()
StatusCode LArRawChannelBuilderSCAlg::initialize |
( |
| ) |
|
|
override |
◆ inputHandles()
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() [1/2]
◆ msg() [2/2]
◆ msgLvl()
◆ outputHandles()
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()
◆ renounceArray()
◆ setFilterPassed()
◆ sysExecute()
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.
◆ sysInitialize()
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 InputMakerBase, and HypoBase.
Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.
110 if (
sc.isFailure()) {
118 if ( cs.retrieve().isFailure() ) {
120 return StatusCode::SUCCESS;
122 if (cs->regHandle(
this,*
h).isFailure()) {
123 sc = StatusCode::FAILURE;
124 ATH_MSG_ERROR(
"unable to register WriteCondHandle " <<
h->fullKey()
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ updateVHKA()
◆ m_absECutFortQ
Gaudi::Property<bool> LArRawChannelBuilderSCAlg::m_absECutFortQ {this,"absECut",true,"Cut on fabs(E) for Q and t computation"} |
|
private |
◆ m_adc2MeVKey
◆ m_bcContKey
◆ m_cablingKey
◆ m_caloSuperCellMgrKey
Initial value:{
this,"CaloSuperCellDetDescrManager","CaloSuperCellDetDescrManager","SG key of the resulting CaloSuperCellDetDescrManager"}
Definition at line 75 of file LArRawChannelBuilderSCAlg.h.
◆ m_cellKey
Initial value:{this,"CaloCellKey","SCellnoBCID",
"SG key of the output CaloCellContainer"}
Definition at line 46 of file LArRawChannelBuilderSCAlg.h.
◆ m_detStore
◆ m_digitKey
◆ m_eCutFortQ
Gaudi::Property<float> LArRawChannelBuilderSCAlg::m_eCutFortQ {this,"ECutFortQ",256.0,"Time and Quality will be computed only for channels with E above this value"} |
|
private |
◆ m_evtStore
◆ m_extendedExtraObjects
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_firstSample
Gaudi::Property<int> LArRawChannelBuilderSCAlg::m_firstSample {this,"firstSample",0,"first of the 32 sampels of the MC shape to be used"} |
|
private |
◆ m_ofcKey
◆ m_onlineId
◆ m_pedestalKey
◆ m_shapeKey
◆ m_useShapeDer
Gaudi::Property<bool> LArRawChannelBuilderSCAlg::m_useShapeDer {this,"useShapeDer",true,"Use shape derivative in Q-factor computation"} |
|
private |
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
def retrieve(aClass, aKey=None)
virtual float pedestal(const HWIdentifier &id, int gain) const =0
const LArVectorProxy ADC2MEV(const HWIdentifier &id, int gain) const
Extra patterns decribing particle interation process.
bool good() const
Returns true if no problems at all (all bits at zero)
SG::ReadCondHandleKey< ILArPedestal > m_pedestalKey
Gaudi::Property< int > m_firstSample
@ VIEW_ELEMENTS
this data object is a view, it does not own its elmts
Gaudi::Property< bool > m_absECutFortQ
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Conditions-Data class holding LAr Bad Channel or Bad Feb information.
SG::ReadCondHandleKey< LArOnOffIdMapping > m_cablingKey
This class groups all DetDescr information related to a CaloCell. Provides a generic interface for al...
virtual OFCRef_t OFC_b(const HWIdentifier &id, int gain, int tbin=0) const =0
const CaloDetDescrElement * get_element(const Identifier &cellId) const
get element by its identifier
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)
virtual ShapeRef_t ShapeDer(const HWIdentifier &id, int gain, int tbin=0, int mode=0) const =0
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
#define ATH_MSG_VERBOSE(x)
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual void setOwner(IDataHandleHolder *o)=0
virtual const DataObjIDColl & extraOutputDeps() const override
Return the list of extra output dependencies.
SG::ReadHandleKey< LArDigitContainer > m_digitKey
virtual float timeOffset(const HWIdentifier &CellID, int gain) const =0
bool isHECchannel(const HWIdentifier id) const override final
bool statusBad(PosType pb) const
Returns true if corresponding status bit its set
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
Liquid Argon digit base class.
::StatusCode StatusCode
StatusCode definition for legacy code.
virtual StatusCode sysExecute(const EventContext &ctx) override
Execute an algorithm.
SG::ReadCondHandleKey< CaloSuperCellDetDescrManager > m_caloSuperCellMgrKey
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Gaudi::Property< float > m_eCutFortQ
const LArOnline_SuperCellID * m_onlineId
SG::WriteHandleKey< CaloCellContainer > m_cellKey
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
SG::ReadCondHandleKey< LArADC2MeV > m_adc2MeVKey
Gaudi::Property< bool > m_useShapeDer
size_type channelHashMax(void) const
Define channel hash tables max size.
SG::ReadCondHandleKey< ILArOFC > m_ofcKey
StatusCode initialize(bool used=true)
SG::ReadCondHandleKey< ILArShape > m_shapeKey
Data object for each calorimeter readout cell.
#define ATH_MSG_WARNING(x)
a typed memory pool that saves time spent allocation small object. This is typically used by containe...
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
SG::ReadCondHandleKey< LArBadChannelCont > m_bcContKey
Bad Channel masking for Super-Cells.
virtual OFCRef_t OFC_a(const HWIdentifier &id, int gain, int tbin=0) const =0
access to OFCs by online ID, gain, and tbin (!=0 for testbeam)
virtual StatusCode sysInitialize() override
Override sysInitialize.
std::string channel_name(const HWIdentifier id) const
Return a string corresponding to a feedthrough name given an identifier.
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
constexpr int pow(int base, int exp) noexcept
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
virtual ShapeRef_t Shape(const HWIdentifier &id, int gain, int tbin=0, int mode=0) const =0