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

#include <LArOFCCondAlg.h>

Inheritance diagram for LArOFCCondAlg:
Collaboration diagram for LArOFCCondAlg:

Public Member Functions

 LArOFCCondAlg (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~LArOFCCondAlg () override
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

unsigned int findTheNSamples (ILArShape::ShapeRef_t Shape, unsigned int nsamples_AC_OFC, unsigned int nsamples_shape) 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< LArOnOffIdMappingm_LArOnOffIdMappingObjKey
SG::ReadCondHandleKey< ILArShapem_LArShapeObjKey
SG::ReadCondHandleKey< ILArNoisem_LArNoiseObjKey
SG::ReadCondHandleKey< ILArPedestalm_LArPedestalObjKey
SG::ReadCondHandleKey< LArAutoCorrTotalm_LArAutoCorrTotalObjKey
SG::WriteCondHandleKey< LArOFCm_LArOFCObjKey
float m_Nminbias
bool m_isMC
bool m_isSuperCell
unsigned int m_firstSample
bool m_useHighestGainAutoCorr
bool m_Dump
size_t m_nGains = 0UL
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

AlgoTool to compute OFC on fly
  • OFC are computed consistently with autocorrelation (by default (NMinBias=-1), these OFC are retrieved)
  • possibility to retrieve OFC for a wanted luminosity (set by NMinBias, the number of minimum bias events per bunch-crossing) but they won't be consistent with autocorrelation anymore. slower since re-do calculation
  • flag FromDB switches between "take OFC stored in DataBase" and "take OFC stored in datamembers (computed consistently from LArAutoCorrTotal) or recomputed for a wanted luminosity"

Definition at line 42 of file LArOFCCondAlg.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

◆ LArOFCCondAlg()

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

Definition at line 22 of file LArOFCCondAlg.cxx.

24 : ::AthAlgorithm(name, pSvcLocator),
25 m_LArOnOffIdMappingObjKey("LArOnOffIdMap"),
26 m_LArShapeObjKey("LArShapeSym"),
27 m_LArNoiseObjKey("LArNoiseSym"),
28 m_LArPedestalObjKey("LArPedestal"),
29 m_LArAutoCorrTotalObjKey("LArAutoCorrTotal"),
30 m_LArOFCObjKey("LArOFC"),
31 m_Nminbias(0), m_isMC(true),
32 m_isSuperCell(false), m_firstSample(0),
33 m_useHighestGainAutoCorr(false), m_Dump(false) {
34 declareProperty("LArOnOffIdMappingObjKey", m_LArOnOffIdMappingObjKey,
35 "Key to read LArOnOffIdMapping object");
36 declareProperty("LArShapeObjKey", m_LArShapeObjKey,
37 "Key to read LArShape object");
38 declareProperty("LArNoiseObjKey", m_LArNoiseObjKey,
39 "Key to read LArNoise object");
40 declareProperty("LArPedestalObjKey", m_LArPedestalObjKey,
41 "Key to read LArPedestal object");
42 declareProperty("LArAutoCorrTotalObjKey", m_LArAutoCorrTotalObjKey,
43 "Key to read LArAutoCorrTotal object");
44 declareProperty("LArOFCObjKey", m_LArOFCObjKey,
45 "Key to write LArOFC object");
46 declareProperty("Nminbias", m_Nminbias);
47 declareProperty("isMC", m_isMC);
48 declareProperty("isSuperCell", m_isSuperCell);
50 "firstSample", m_firstSample,
51 "First sample to use for in-time event on the full pulse shape");
52 declareProperty("useHighestGainAutoCorr",m_useHighestGainAutoCorr);
53 declareProperty("DumpOFCCondAlg",m_Dump);
54 }
AthAlgorithm()
Default constructor:
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
unsigned int m_firstSample
SG::WriteCondHandleKey< LArOFC > m_LArOFCObjKey
SG::ReadCondHandleKey< LArOnOffIdMapping > m_LArOnOffIdMappingObjKey
bool m_useHighestGainAutoCorr
SG::ReadCondHandleKey< ILArShape > m_LArShapeObjKey
SG::ReadCondHandleKey< LArAutoCorrTotal > m_LArAutoCorrTotalObjKey
SG::ReadCondHandleKey< ILArNoise > m_LArNoiseObjKey
SG::ReadCondHandleKey< ILArPedestal > m_LArPedestalObjKey

◆ ~LArOFCCondAlg()

LArOFCCondAlg::~LArOFCCondAlg ( )
overridevirtual

Definition at line 56 of file LArOFCCondAlg.cxx.

56{}

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

StatusCode LArOFCCondAlg::execute ( )
overridevirtual

Definition at line 84 of file LArOFCCondAlg.cxx.

84 {
85
86 // WriteHandle setup
87 SG::WriteCondHandle<LArOFC> writeHandle(m_LArOFCObjKey);
88 // So the following should not be called usually?!
89 if (writeHandle.isValid()) {
91 "CondHandle "
92 << writeHandle.fullKey() << " is already valid.");
93 return StatusCode::SUCCESS;
94 }
95
96 // Identifier helper
97 // Copied from LArADC2MeVCondAlg.cxx
98 const LArOnlineID_Base *larOnlineID = nullptr;
99 if (m_isSuperCell) {
100 const LArOnline_SuperCellID *scidhelper;
101 ATH_CHECK(detStore()->retrieve(scidhelper, "LArOnline_SuperCellID"));
102 larOnlineID = scidhelper; // cast to base-class
103 } else { // regular cells
104 const LArOnlineID *idhelper;
105 ATH_CHECK(detStore()->retrieve(idhelper, "LArOnlineID"));
106 larOnlineID = idhelper; // cast to base-class
107 }
108
109 // Mapping helper
110 const LArOnOffIdMapping *larOnOffIdMapping = nullptr;
111 SG::ReadCondHandle<LArOnOffIdMapping> larOnOffIdMappingHdl{
113 };
114 larOnOffIdMapping = *larOnOffIdMappingHdl;
115 if (larOnOffIdMapping == nullptr) {
116 ATH_MSG_ERROR("Failed to retrieve LArOnOffIdMapping object");
117 return StatusCode::FAILURE;
118 }
119 writeHandle.addDependency(larOnOffIdMappingHdl);
120
121 // Get pointers to inputs
122 // Retrieve validity ranges and determine their intersection
123 SG::ReadCondHandle<ILArShape> ShapeHdl{ m_LArShapeObjKey };
124 // FIXME: should check if handle is properly created and/or check if handle is
125 // properly retrieved
126 // operator star of a ReadCondHandle returns a const pointer to type T
127 const ILArShape *larShape{ *ShapeHdl };
128 if (larShape == nullptr) {
129 ATH_MSG_ERROR("Failed to retrieve LArShape object");
130 return StatusCode::FAILURE;
131 }
132 writeHandle.addDependency(ShapeHdl);
133
134 SG::ReadCondHandle<LArAutoCorrTotal> AutoCorrTotalHdl{ m_LArAutoCorrTotalObjKey };
135 const LArAutoCorrTotal *larAutoCorrTotal = nullptr;
136 larAutoCorrTotal= *AutoCorrTotalHdl;
137 if (larAutoCorrTotal == nullptr) {
138 ATH_MSG_ERROR("Failed to retrieve LArADC2MeV object");
139 return StatusCode::FAILURE;
140 }
141 writeHandle.addDependency(AutoCorrTotalHdl);
142
143 // Consider the determinstic objects
144 const ILArNoise *larNoise = nullptr;
145 const ILArPedestal *larPedestal = nullptr;
146
147 if (m_isMC) {
148 SG::ReadCondHandle<ILArNoise> NoiseHdl{ m_LArNoiseObjKey };
149 larNoise = *NoiseHdl;
150 if (larNoise == nullptr) {
151 ATH_MSG_ERROR("Failed to retrieve object LArNoise");
152 return StatusCode::FAILURE;
153 }
154 writeHandle.addDependency(NoiseHdl);
155 } else {
156 SG::ReadCondHandle<ILArPedestal> PedestalHdl{ m_LArPedestalObjKey };
157 larPedestal = *PedestalHdl;
158 if (larPedestal == nullptr) {
159 ATH_MSG_ERROR("Failed to retrieve object LArPedestal");
160 return StatusCode::FAILURE;
161 }
162 writeHandle.addDependency(PedestalHdl);
163 }
164
165 ATH_MSG_INFO("IOV found from intersection for LArOFCCondObj object: "
166 << writeHandle.getRange());
167
168 // make output object
169 // dimensions: number of gains x number of channel IDs x elements of
170 // OFC
171 std::unique_ptr<LArOFC> larOFC =
172 std::make_unique<LArOFC>(larOnlineID, larOnOffIdMapping, m_nGains);
173
174 std::vector<float> OFCa_tmp, OFCb_tmp;
175
176 for (const HWIdentifier chid : larOnlineID->channel_range()) {
177 const IdentifierHash hid = larOnlineID->channel_Hash(chid);
178
179 //if (!(larOnOffIdMapping->isOnlineConnected(chid))) continue;
180 if (larOnOffIdMapping->isOnlineConnected(chid)) {
181 for (size_t igain = 0; igain < m_nGains; igain++) {
182
183
184 //:::::::::::::::::::::::::::::::
185 //retrieve the data
186 //:::::::::::::::::::::::::::::::
187 ILArShape::ShapeRef_t Shape = larShape->Shape(chid,igain);
188 unsigned int nsamples_shape = Shape.size();
189 ILArShape::ShapeRef_t ShapeDer = larShape->ShapeDer(chid,igain);
190 //:::::::::::::::::::::::::::::::
191
192 // get Noise autocorrelation for gain
193 int igain_autocorr = 0;
194 // to use only Autocorr fro highest gain in optimization: HEC/FCAL=> medium gain EM=>high gain
196 &&(larOnlineID->isHECchannel(chid) || larOnlineID->isFCALchannel(chid)) ) igain_autocorr=1;
197
198 const std::vector<double> AutoCorr =
199 larAutoCorrTotal->autoCorrTotal(chid,igain_autocorr,m_Nminbias);
200 //unsigned int nsamples_AC_OFC=AutoCorr.size()+1;
201 unsigned int nsamples_AC_OFC = (1+((int)(sqrt(1+8*AutoCorr.size()))))/2;
202
203 const std::vector<double>& rmsSampl =
204 larAutoCorrTotal->samplRMS(chid,igain_autocorr,m_Nminbias);
205 unsigned int nsamples2 = rmsSampl.size();
206 if (nsamples2 != nsamples_AC_OFC) {
207 ATH_MSG_WARNING( " bad size for rmsSampl " );
208 //return (m_OFCtmp); // return empty vector
209 }
210 //:::::::::::::::::::::::::::::::
211 //unsigned int iBeginOfNSamples=findTheNSamples(Shape,
212 // nsamples_AC_OFC,
213 // nsamples_shape);
214 unsigned int firstSample = m_firstSample;
215 if(larOnlineID->isHECchannel(chid) && m_firstSample == 0 && nsamples_AC_OFC==4){
216 firstSample=1;
217 }
218 unsigned int iBeginOfNSamples = firstSample;
219 if(nsamples_AC_OFC + iBeginOfNSamples > nsamples_shape)
220 iBeginOfNSamples=0;
221 //:::::::::::::::::::::::::::::::
222
223 if(m_isMC) {
224 }
225 else
226 {
227 float RMSpedestal = larPedestal->pedestalRMS(chid,igain);
228 if(RMSpedestal>= (1.0+LArElecCalib::ERRORCODE))
229 ;
230 else
231 {
232 ATH_MSG_WARNING(" PedestalRMS vector empty for "
233 <<chid<<" at gain "<<igain );
234 }
235 }
236 //:::::::::::::::::::::::::::::::
237 //protection against missing data
238 //:::::::::::::::::::::::::::::::
239 if(Shape.size()==0 || ShapeDer.size()==0 || AutoCorr.empty())
240 {
241 ATH_MSG_WARNING("Some data are missing -> OFC will be empty for "
242 <<chid<<" at gain "<<igain );
243 //return (m_OFCtmp);
244 //returns an empty vector
245 }
246 if (Shape.size()!=ShapeDer.size()) {
247 ATH_MSG_ERROR("Got invalid shape object: Size=" << Shape.size() << ", DerSize=" << ShapeDer.size());
248 return StatusCode::SUCCESS;
249 }
250 //:::::::::::::::::::::::::::::::
251 unsigned int l,c,i;
252 //:::::::::::::::::::::::::::::::
253 //calculations
254 //:::::::::::::::::::::::::::::::
255 // fill and inverrt AC matrix
256 //HepMatrix AC(nsamples_AC_OFC,nsamples_AC_OFC),
257 //ACinv(nsamples_AC_OFC,nsamples_AC_OFC);
258 Eigen::MatrixXf AC = Eigen::MatrixXf::Zero(nsamples_AC_OFC,nsamples_AC_OFC);
259 Eigen::MatrixXf ACinv = Eigen::MatrixXf::Zero(nsamples_AC_OFC,nsamples_AC_OFC);
260 for(l=0;l<nsamples_AC_OFC;++l) { // l=line c=column
261 for(c=0;c<nsamples_AC_OFC;++c) {
262 if (l==c) {
263 AC(l,c)=1.;
264 }
265 else {
266 int i1=std::min(l,c);
267 int i2=std::max(l,c);
268 int index = i1*nsamples_AC_OFC - i1*(i1+1)/2 -(i1+1) + i2;
269 AC(l,c)=AutoCorr[index];
270 }
271 AC(l,c) = AC(l,c)*rmsSampl[l]*rmsSampl[c];
272 }
273 }
274 ACinv=AC.inverse();
275 //:::::::::::::::::::::::::::::::
276
277
278 float ACinv_PS[32];//ACinv_PS
279 float ACinv_PSD[32]; //ACinv_PSD
280 //Q1 Q2 Q3 DELTA
281 float Q1=0.;
282 float Q2=0.;
283 float Q3=0.;
284
285 for(l=0;l<nsamples_AC_OFC;++l)
286 {
287 ACinv_PS[l]=0.; ACinv_PSD[l]=0.;
288 for(c=0;c<nsamples_AC_OFC;++c){
289 ACinv_PS[l]+=ACinv(l,c)*Shape[c+iBeginOfNSamples];
290 ACinv_PSD[l]+=ACinv(l,c)*ShapeDer[c+iBeginOfNSamples];
291 }
292 Q1+=Shape[l+iBeginOfNSamples]*ACinv_PS[l];
293 Q2+=ShapeDer[l+iBeginOfNSamples]*ACinv_PSD[l];
294 Q3+=ShapeDer[l+iBeginOfNSamples]*ACinv_PS[l];
295 }
296 float DELTA=Q1*Q2-Q3*Q3;
297 //:::::::::::::::::::::::::::::::
298 //OFCa
299 OFCa_tmp.resize(nsamples_AC_OFC);
300 for(i=0;i<nsamples_AC_OFC;++i)
301 OFCa_tmp[i]=(ACinv_PS[i]*Q2-ACinv_PSD[i]*Q3)/DELTA;
302 //OFCb
303 OFCb_tmp.resize(nsamples_AC_OFC);
304 for(i=0;i<nsamples_AC_OFC;++i)
305 OFCb_tmp[i]=(ACinv_PS[i]*Q3-ACinv_PSD[i]*Q1)/DELTA;
306
307 //for debugging only
308 if(m_Dump)
309 {
310 std::cout<<larOnlineID
311 ->show_to_string(larOnOffIdMapping->cnvToIdentifier(chid))
312 <<" gain="<<igain<<" Nminbias="<<m_Nminbias<<std::endl;
313 std::cout<<"Shape: ";
314 for(c=0;c<nsamples_shape;++c)
315 std::cout<<Shape[c]<<" ";
316 std::cout<<std::endl;
317 std::cout<<"ShapeDer: ";
318 for(c=0;c<nsamples_shape;++c)
319 std::cout<<ShapeDer[c]<<" ";
320 std::cout<<std::endl;
321 for(c=0;c<nsamples_AC_OFC;++c)
322 std::cout<<Shape[c+iBeginOfNSamples]<<" ";
323 std::cout<<" <- "<<iBeginOfNSamples<<std::endl;
324 for(i=0;i<nsamples_AC_OFC;++i) std::cout<<ACinv_PS[i]<<" ";
325 std::cout<<std::endl;
326 for(i=0;i<nsamples_AC_OFC;++i) std::cout<<ACinv_PSD[i]<<" ";
327 std::cout<<std::endl;
328 std::cout<<" Q1="<<Q1<<" Q2="<<Q2<<" Q3="<<Q3
329 <<" DELTA="<<DELTA<<std::endl;
330 std::cout << " OFCa: ";
331 for(i=0;i<nsamples_AC_OFC;++i)
332 std::cout<<(ACinv_PS[i]*Q2-ACinv_PSD[i]*Q3)/DELTA<<" ";
333 std::cout<<std::endl;
334 }
335 bool stat = larOFC->setOFC(hid, igain, std::make_pair(OFCa_tmp, OFCb_tmp));
336 if (!stat) {
337 msg(MSG::ERROR) << "LArOFC::setOFC fails for gain " << igain << ", hash " << hid << endmsg;
338 }
339 } // end loop over gains
340 } else { // end loop over connected channels -- now, set empty for disc. chnanels
341 for (unsigned igain=0;igain<m_nGains;++igain) {
342 std::vector<float> empty;
343 bool stat = larOFC->setOFC(hid,igain, std::make_pair(empty, empty));
344 if (!stat) {
345 msg(MSG::ERROR) << "LArOFC::setOFC fails for gain " << igain << ", hash " << hid << endmsg;
346 }
347 } // end loop over gains of disconnected channels
348 } // end loop over disconnected channels
349
350 } // end loop over all channels
351
352
353 ATH_CHECK(writeHandle.record(std::move(larOFC)));
354 ATH_MSG_INFO("Wrote LArOFC obj to CondStore");
355 return StatusCode::SUCCESS;
356}
#define endmsg
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_ERROR(x)
#define ATH_MSG_INFO(x)
#define ATH_MSG_WARNING(x)
#define ATH_MSG_DEBUG(x)
static const double DELTA
static const Attributes_t empty
const ServiceHandle< StoreGateSvc > & detStore() const
MsgStream & msg() const
std::string show_to_string(Identifier id, const IdContext *context=0, char sep='.') const
or provide the printout in string form
virtual float pedestalRMS(const HWIdentifier &id, int gain) const =0
access to RMS of Pedestal index by Identifier, and gain setting
LArVectorProxy ShapeRef_t
This class defines the interface for accessing Shape (Nsample variable, Dt = 25 ns fixed) @stereotype...
Definition ILArShape.h:26
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
const std::vector< double > autoCorrTotal(const IdentifierHash &hid, int gain, float Nminbias) const
const std::vector< double > samplRMS(const IdentifierHash &hid, int gain, float Nminbias) const
Identifier cnvToIdentifier(const HWIdentifier &sid) const
create an Identifier from a HWIdentifier (inline)
bool isOnlineConnected(const HWIdentifier &sid) const
Test whether a HWIdentifier is connected of not (inline)
id_range channel_range() const
IdentifierHash channel_Hash(HWIdentifier channelId) const
Create channel_hash from channel_Id.
bool isFCALchannel(const HWIdentifier id) const
virtual bool isHECchannel(const HWIdentifier id) const =0
str index
Definition DeMoScan.py:362
l
Printing final latex table to .tex output file.
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

◆ findTheNSamples()

unsigned int LArOFCCondAlg::findTheNSamples ( ILArShape::ShapeRef_t Shape,
unsigned int nsamples_AC_OFC,
unsigned int nsamples_shape ) const
private

Definition at line 358 of file LArOFCCondAlg.cxx.

361{
362 unsigned int i_ShapeMax=0;
363 double ShapeMax=0;
364 for(unsigned int i=0;i<nsamples_shape;++i)
365 {
366 double value=Shape[i];
367 if(value>ShapeMax) { ShapeMax=value; i_ShapeMax=i; }
368 else if(value<0 && i>3) break;//after the peak
369 }
370
371 unsigned int tmp=int(nsamples_AC_OFC/2.);
372 if(tmp>i_ShapeMax) return 0;
373 else return i_ShapeMax-tmp;
374}

◆ initialize()

StatusCode LArOFCCondAlg::initialize ( )
overridevirtual

Definition at line 58 of file LArOFCCondAlg.cxx.

58 {
59 ATH_MSG_DEBUG("initialize " << name());
60
61 // ReadCondHandle initialization
62 ATH_CHECK(m_LArShapeObjKey.initialize());
63
66
67 //WriteHandle initialization
68 ATH_CHECK(m_LArOFCObjKey.initialize());
69
72
73 // Number of gains (does this have to be in initialize now b/c of AthenaMT?)
74 // Copied from LArADC2MeVCondAlg.cxx
75 if (m_isSuperCell) {
76 m_nGains = 1;
77 } else {
78 m_nGains = 3;
79 }
80
81 return StatusCode::SUCCESS;
82}

◆ 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}
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_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 LArOFCCondAlg::m_Dump
private

Definition at line 69 of file LArOFCCondAlg.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_firstSample

unsigned int LArOFCCondAlg::m_firstSample
private

Definition at line 62 of file LArOFCCondAlg.h.

◆ m_isMC

bool LArOFCCondAlg::m_isMC
private

Definition at line 60 of file LArOFCCondAlg.h.

◆ m_isSuperCell

bool LArOFCCondAlg::m_isSuperCell
private

Definition at line 61 of file LArOFCCondAlg.h.

◆ m_LArAutoCorrTotalObjKey

SG::ReadCondHandleKey<LArAutoCorrTotal> LArOFCCondAlg::m_LArAutoCorrTotalObjKey
private

Definition at line 55 of file LArOFCCondAlg.h.

◆ m_LArNoiseObjKey

SG::ReadCondHandleKey<ILArNoise> LArOFCCondAlg::m_LArNoiseObjKey
private

Definition at line 53 of file LArOFCCondAlg.h.

◆ m_LArOFCObjKey

SG::WriteCondHandleKey<LArOFC> LArOFCCondAlg::m_LArOFCObjKey
private

Definition at line 57 of file LArOFCCondAlg.h.

◆ m_LArOnOffIdMappingObjKey

SG::ReadCondHandleKey<LArOnOffIdMapping> LArOFCCondAlg::m_LArOnOffIdMappingObjKey
private

Definition at line 51 of file LArOFCCondAlg.h.

◆ m_LArPedestalObjKey

SG::ReadCondHandleKey<ILArPedestal> LArOFCCondAlg::m_LArPedestalObjKey
private

Definition at line 54 of file LArOFCCondAlg.h.

◆ m_LArShapeObjKey

SG::ReadCondHandleKey<ILArShape> LArOFCCondAlg::m_LArShapeObjKey
private

Definition at line 52 of file LArOFCCondAlg.h.

◆ m_nGains

size_t LArOFCCondAlg::m_nGains = 0UL
private

Definition at line 71 of file LArOFCCondAlg.h.

◆ m_Nminbias

float LArOFCCondAlg::m_Nminbias
private

Definition at line 59 of file LArOFCCondAlg.h.

◆ m_useHighestGainAutoCorr

bool LArOFCCondAlg::m_useHighestGainAutoCorr
private

Definition at line 63 of file LArOFCCondAlg.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: