Loading [MathJax]/extensions/tex2jax.js
ATLAS Offline Software
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
Functions
LArFEBMonAlg Namespace Reference

Functions

def LArFEBMonConfig (flags, cellDebug=False, dspDebug=False)
 
def LArFEBMonConfigCore (helper, algoinstance, flags, cellDebug=False, dspDebug=False)
 

Function Documentation

◆ LArFEBMonConfig()

def LArFEBMonAlg.LArFEBMonConfig (   flags,
  cellDebug = False,
  dspDebug = False 
)

Definition at line 4 of file LArFEBMonAlg.py.

4 def LArFEBMonConfig(flags, cellDebug=False, dspDebug=False):
5  from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
6  from AthenaMonitoring import AthMonitorCfgHelper
7  helper = AthMonitorCfgHelper(flags,'LArFEBMonAlgCfg')
8 
9  from AthenaConfiguration.ComponentFactory import CompFactory
10  LArFEBMonConfigCore(helper, CompFactory.LArFEBMonAlg,flags,cellDebug, dspDebug)
11 
12  rv = ComponentAccumulator()
13 
14  # adding LArFebErrorSummary algo
15  from LArROD.LArFebErrorSummaryMakerConfig import LArFebErrorSummaryMakerCfg
16  rv.merge(LArFebErrorSummaryMakerCfg(flags))
17 
18  rv.merge(helper.result())
19 
20  return rv
21 

◆ LArFEBMonConfigCore()

def LArFEBMonAlg.LArFEBMonConfigCore (   helper,
  algoinstance,
  flags,
  cellDebug = False,
  dspDebug = False 
)

Definition at line 22 of file LArFEBMonAlg.py.

22 def LArFEBMonConfigCore(helper,algoinstance,flags, cellDebug=False, dspDebug=False):
23 
24  from LArMonitoring.GlobalVariables import lArDQGlobals
25 
26  larFEBMonAlg = helper.addAlgorithm(algoinstance,'larFEBMonAlg')
27 
28  GroupName="FEBMon"
29  nslots=[]
30  for i in range(0,len(lArDQGlobals.FEB_Slot)):
31  nslots.append(lArDQGlobals.FEB_Slot[lArDQGlobals.Partitions[i]][1])
32 
33  larFEBMonAlg.MonGroup=GroupName
34  larFEBMonAlg.PartitionNames=lArDQGlobals.Partitions
35  larFEBMonAlg.SubDetNames=lArDQGlobals.SubDet
36  larFEBMonAlg.Streams=lArDQGlobals.defaultStreamNames
37 
38  isCOMP200="COMP200" in flags.IOVDb.DatabaseInstance
39 
40  if not isCOMP200:
41  dbString="<db>COOLONL_LAR/CONDBR2</db>"
42  persClass="AthenaAttributeList"
43  fld="/LAR/Configuration/DSPThresholdFlat/Thresholds"
44  havethem=False
45  for c in helper.resobj.getServices():
46  if c.getName()=="IOVDbSvc":
47  iovDbSvc=c
48  condLoader=helper.resobj.getCondAlgo("CondInputLoader")
49  havethem=True
50  break
51 
52  if not havethem:
53  from IOVDbSvc.IOVDbSvcConfig import IOVDbSvcCfg
54  helper.resobj.merge(IOVDbSvcCfg(flags))
55  condLoader=helper.resobj.getCondAlgo("CondInputLoader")
56  iovDbSvc=helper.resobj.getService("IOVDbSvc")
57 
58  iovDbSvc.Folders.append(fld+dbString)
59  condLoader.Load.add((persClass,fld))
60  larFEBMonAlg.Run2DSPThresholdsKey = fld
61  else:
62  fld='/LAR/Configuration/DSPThreshold/Thresholds'
63  db='LAR_ONL'
64  obj='LArDSPThresholdsComplete'
65  from IOVDbSvc.IOVDbSvcConfig import addFolders
66  helper.resobj.merge(addFolders(flags,fld,db,obj))
67  larFEBMonAlg.Run1DSPThresholdsKey = 'LArDSPThresholds'
68 
69 
70  Group = helper.addGroup(
71  larFEBMonAlg,
72  GroupName,
73  '/LAr/'+GroupName+'/'
74  )
75 
76 
77  #Summary histos
78  summary_hist_path='Summary/'
79 
80  #-- TTree for corrupted events timestamp
81  Group.defineTree('timestamp,time_ns,febHwId,febErrorType;LArCorrupted',
82  path=summary_hist_path,
83  title='Timestamps of corrupted LAr events',
84  treedef='timestamp/i:time_ns/i:febHwId/vector<int>:febErrorType/vector<int>')
85 
86  Group.defineHistogram('nbFEB;NbOfReadoutFEBGlobal',
87  title='# of readout FEB/DSP header',
88  type='TH1I',
89  path=summary_hist_path,
90  opt='kAlwaysCreate',
91  xbins=lArDQGlobals.N_FEB+11, xmin=-0.5, xmax=lArDQGlobals.N_FEB+10+0.5)
92  Group.defineHistogram('nbFEBpart,part;NbOfEvts2d',
93  title='# of readout FEB/DSP header;Num. FEBs;Partition',
94  type='TH2I',
95  path=summary_hist_path,
96  opt='kAlwaysCreate',
97  xbins=lArDQGlobals.N_FEB_Parttions_Max, xmin=-0.5, xmax=lArDQGlobals.N_FEB_Parttions_Max-0.5,
98  ybins=lArDQGlobals.N_Partitions, ymin=-0.5, ymax=lArDQGlobals.N_Partitions-0.5,
99  ylabels=lArDQGlobals.Partitions
100  )
101  Group.defineHistogram('febError,part;NbOfLArFEBMonErrors_dE',
102  title='# of data corruption errors',
103  type='TH2I',
104  path=summary_hist_path,
105  opt='kAlwaysCreate',
106  xbins=lArDQGlobals.N_FEBErrors, xmin=0.5, xmax=lArDQGlobals.N_FEBErrors+0.5,
107  ybins=lArDQGlobals.N_Partitions, ymin=-0.5, ymax=lArDQGlobals.N_Partitions-0.5,
108  xlabels=lArDQGlobals.FEBErrors, ylabels=lArDQGlobals.Partitions)
109  Group.defineHistogram('dspThrADC;dspThresholdsADC',
110  title='DSP thresholds to readout samples;Number of cells;Cell threshold in ADC counts',
111  type='TH1I',
112  path=summary_hist_path,
113  xbins=lArDQGlobals.DSPThr_Bins+1, xmin=-0.5, xmax=lArDQGlobals.DSPThr_Bins+0.5,
114  merge='identical')
115  Group.defineHistogram('dspThrQT;dspThresholds_qfactortime',
116  title='DSP thresholds to readout (qfactor+time);Number of cells;Cell threshold in ADC counts',
117  type='TH1I',
118  path=summary_hist_path,
119  xbins=lArDQGlobals.DSPThr_Bins+1, xmin=-0.5, xmax=lArDQGlobals.DSPThr_Bins+0.5,
120  merge='identical')
121  Group.defineHistogram('EvtType;Eventtype',
122  title='Event type (1st readout FEB)',
123  type='TH1I',
124  path=summary_hist_path,
125  xbins=lArDQGlobals.Evt_Bins, xmin=lArDQGlobals.Evt_Min, xmax=lArDQGlobals.Evt_Max,
126  xlabels=lArDQGlobals.Evt_labels)
127  Group.defineHistogram('LVL1Trig;TriggerWord',
128  title='Number of Events per L1 trigger word (8 bits);L1 trigger word',
129  type='TH1I',
130  path=summary_hist_path,
131  xbins=lArDQGlobals.L1Trig_Bins, xmin=lArDQGlobals.L1Trig_Min, xmax=lArDQGlobals.L1Trig_Max)
132  Group.defineHistogram('LVL1TrigAllDSP;TriggerWordAllDSP',
133  title='Number of L1 trigger word per DSP (8 bits);L1 trigger word',
134  type='TH1I',
135  path=summary_hist_path,
136  xbins=lArDQGlobals.L1Trig_Bins, xmin=lArDQGlobals.L1Trig_Min, xmax=lArDQGlobals.L1Trig_Max)
137  Group.defineHistogram('EvtRej;EventsRejected',
138  title='Nb of events rejected (at least one error)',
139  type='TH1I',
140  path=summary_hist_path,
141  xbins=lArDQGlobals.EvtRej_Bins, xmin=lArDQGlobals.EvtRej_Min, xmax=lArDQGlobals.EvtRej_Max,
142  xlabels=lArDQGlobals.EvtRej_labels)
143  Group.defineHistogram('LB0,EvtRejYield;RAW_YieldOfRejectedEventsVsLB',
144  title='Yield of corrupted events (DATACORRUPTED);Luminosity Block;Yield(%)',
145  type='TProfile',
146  path=summary_hist_path,
147  xbins=lArDQGlobals.LB_Bins, xmin=lArDQGlobals.LB_Min, xmax=lArDQGlobals.LB_Max)
148  Group.defineHistogram('LB0,EvtRejYieldOut;RAW_YieldOfRejectedEventsVsLBout',
149  title='Yield of corrupted events (DATACORRUPTED) not vetoed by time window;Luminosity Block;Yield(%)',
150  type='TProfile',
151  path=summary_hist_path,
152  xbins=lArDQGlobals.LB_Bins, xmin=lArDQGlobals.LB_Min, xmax=lArDQGlobals.LB_Max)
153  Group.defineHistogram('LB0,EvtOneErrorYield;RAW_YieldOfOneErrorEventsVsLB',
154  title='Yield of events with >=1 FEB in error;Luminosity Block;Yield(%)',
155  type='TProfile',
156  path=summary_hist_path,
157  xbins=lArDQGlobals.LB_Bins, xmin=lArDQGlobals.LB_Min, xmax=lArDQGlobals.LB_Max)
158  Group.defineHistogram('rejBits;rejectionBits',
159  title='Errors at the origin of event rejection;Bits;Number of (rejected) events',
160  type='TH1I',
161  path=summary_hist_path,
162  xbins=lArDQGlobals.rejBits_Bins, xmin=-0.5, xmax=lArDQGlobals.rejBits_Bins-0.5)
163  Group.defineHistogram('LB0;NbOfEventsVsLB',
164  title='Nb of events per LB;Luminosity Block',
165  type='TH1I',
166  path=summary_hist_path,
167  xbins=lArDQGlobals.LB_Bins, xmin=lArDQGlobals.LB_Min, xmax=lArDQGlobals.LB_Max)
168  Group.defineHistogram('NbOfSweet2;NbOfSw2',
169  title='# of cells with samples readout;Number of cells;Number of events',
170  type='TH1I',
171  path=summary_hist_path,
172  xbins=int(lArDQGlobals.N_Cells/10), xmin=-1000, xmax=lArDQGlobals.N_Cells-1000)
173  Group.defineHistogram('LB0,LArEvSize;eventSizeVsLB',
174  title='LAr event size (w/o ROS headers);Luminosity Block;Megabytes',
175  type='TProfile',
176  path=summary_hist_path,
177  xbins=lArDQGlobals.LB_Bins, xmin=lArDQGlobals.LB_Min, xmax=lArDQGlobals.LB_Max)
178  Group.defineHistogram('NbOfSamp;NbOfSamples',
179  title='# of samples (1st readout FEB);Samples;Number of events',
180  type='TH1I',
181  path=summary_hist_path,
182  xbins=lArDQGlobals.Samples_Bins, xmin=lArDQGlobals.Samples_Min, xmax=lArDQGlobals.Samples_Max)
183 
184  isOnline = flags.DQ.Environment=='online'
185 
186  if isOnline:
187  Group.defineHistogram('LBf,EvtRejYield;RAW_EventsRejectedLB',
188  title='% of events rejected in current LB (online only)',
189  type='TProfile',
190  path=summary_hist_path,
191  opt='kLBNHistoryDepth=10,kAlwaysCreate',
192  xbins=1, xmin=0, xmax=1, xlabels=['% of events'])
193  Group.defineHistogram('LB,streamBin,LArEvSizePart;eventSizeStreamVsLB',
194  title='LAr event size per stream per LB (w/o ROS headers)',
195  type='TProfile2D',
196  path=summary_hist_path,
197  xbins=lArDQGlobals.LB_Bins, xmin=lArDQGlobals.LB_Min, xmax=lArDQGlobals.LB_Max,
198  ybins=len(larFEBMonAlg.Streams),ymin=-0.5, ymax= len(larFEBMonAlg.Streams)-0.5,
199  ylabels=list(larFEBMonAlg.Streams)
200  )
201 
202  # Now per partition histograms
203  for subdet in range(0,lArDQGlobals.N_SubDet):
204  hist_path='/LAr/'+GroupName+'/'+lArDQGlobals.SubDet[subdet]+'/'
205  slot_low = lArDQGlobals.FEB_Slot[lArDQGlobals.Partitions[subdet*2]][0] - 0.5
206  slot_up = lArDQGlobals.FEB_Slot[lArDQGlobals.Partitions[subdet*2]][1] + 0.5
207  slot_n = int(slot_up - slot_low)
208  ft_low = lArDQGlobals.FEB_Feedthrough[lArDQGlobals.Partitions[subdet*2]][0] - 0.5
209  ft_up = lArDQGlobals.FEB_Feedthrough[lArDQGlobals.Partitions[subdet*2]][1] + 0.5
210  ft_n = int(ft_up - ft_low)
211 
212  darray = helper.addArray([lArDQGlobals.Partitions[2*subdet:2*subdet+2]],larFEBMonAlg,lArDQGlobals.SubDet[subdet],topPath='/')
213 
214  darray.defineHistogram('slotPar,FTPar;RAW_Parity',
215  title='Parity error;Slot;FT',
216  type='TH2I',
217  path=hist_path,
218  xbins=slot_n,xmin=slot_low,xmax=slot_up,
219  ybins=ft_n, ymin=ft_low, ymax=ft_up)
220 
221  darray.defineHistogram('slotBcid,FTBcid;RAW_BCID',
222  title='BCID mismatch betw. 2 halves of FEB;Slot;FT',
223  type='TH2I',
224  path=hist_path,
225  xbins=slot_n,xmin=slot_low,xmax=slot_up,
226  ybins=ft_n, ymin=ft_low, ymax=ft_up)
227 
228  darray.defineHistogram('slotRadd,FTRadd;RAW_RADD',
229  title='Sample header mismatch betw. 2 halves of FEB;Slot;FT',
230  type='TH2I',
231  path=hist_path,
232  xbins=slot_n,xmin=slot_low,xmax=slot_up,
233  ybins=ft_n, ymin=ft_low, ymax=ft_up)
234 
235  darray.defineHistogram('slotEvtid,FTEvtid;RAW_EVTID',
236  title='EVTID mismatch betw. 2 halves of FEB;Slot;FT',
237  type='TH2I',
238  path=hist_path,
239  xbins=slot_n,xmin=slot_low,xmax=slot_up,
240  ybins=ft_n, ymin=ft_low, ymax=ft_up)
241 
242  darray.defineHistogram('slotScac,FTScac;RAW_SCACStatus',
243  title='Wrong SCAC status in one half of a FEB;Slot;FT',
244  type='TH2I',
245  path=hist_path,
246  xbins=slot_n,xmin=slot_low,xmax=slot_up,
247  ybins=ft_n, ymin=ft_low, ymax=ft_up)
248 
249  darray.defineHistogram('slotscout,FTscout;RAW_scaOutOfRange',
250  title='Sca out of range;Slot;FT',
251  type='TH2I',
252  path=hist_path,
253  xbins=slot_n,xmin=slot_low,xmax=slot_up,
254  ybins=ft_n, ymin=ft_low, ymax=ft_up)
255 
256  darray.defineHistogram('slotgain,FTgain;RAW_gainMismatch',
257  title='Gain mismatch within time samples;Slot;FT',
258  type='TH2I',
259  path=hist_path,
260  xbins=slot_n,xmin=slot_low,xmax=slot_up,
261  ybins=ft_n, ymin=ft_low, ymax=ft_up)
262 
263  darray.defineHistogram('slottype,FTtype;RAW_typeMismatch',
264  title='Event type mismatch;Slot;FT',
265  type='TH2I',
266  path=hist_path,
267  xbins=slot_n,xmin=slot_low,xmax=slot_up,
268  ybins=ft_n, ymin=ft_low, ymax=ft_up)
269 
270  darray.defineHistogram('slotsmp,FTsmp;RAW_badNbOfSamp',
271  title='Non uniform number of samples;Slot;FT',
272  type='TH2I',
273  path=hist_path,
274  xbins=slot_n,xmin=slot_low,xmax=slot_up,
275  ybins=ft_n, ymin=ft_low, ymax=ft_up)
276 
277  darray.defineHistogram('slotzero,FTzero;RAW_zeroSamp',
278  title='Empty FEB data blocks;Slot;FT',
279  type='TH2I',
280  path=hist_path,
281  xbins=slot_n,xmin=slot_low,xmax=slot_up,
282  ybins=ft_n, ymin=ft_low, ymax=ft_up)
283 
284  darray.defineHistogram('slotsum,FTsum;RAW_checkSum',
285  title='Checksum / DSP block size;Slot;FT',
286  type='TH2I',
287  path=hist_path,
288  opt='kAlwaysCreate',
289  xbins=slot_n,xmin=slot_low,xmax=slot_up,
290  ybins=ft_n, ymin=ft_low, ymax=ft_up)
291 
292  darray.defineHistogram('slotmis,FTmis;RAW_missingHeader',
293  title='Missing header ;Slot;FT',
294  type='TH2I',
295  path=hist_path,
296  xbins=slot_n,xmin=slot_low,xmax=slot_up,
297  ybins=ft_n, ymin=ft_low, ymax=ft_up)
298 
299  darray.defineHistogram('slotgain,FTgain;RAW_badGain',
300  title='Bad gain ;Slot;FT',
301  type='TH2I',
302  path=hist_path,
303  xbins=slot_n,xmin=slot_low,xmax=slot_up,
304  ybins=ft_n, ymin=ft_low, ymax=ft_up)
305 
306  darray.defineHistogram('slotabs,FTabs;RAW_LArFEBMonErrorsAbsolute',
307  title='Nb of events with at least one error ;Slot;FT',
308  type='TH2I',
309  path=hist_path,
310  opt='kAlwaysCreate',
311  xbins=slot_n,xmin=slot_low,xmax=slot_up,
312  ybins=ft_n, ymin=ft_low, ymax=ft_up)
313 
314  darray.defineHistogram('slotmist,FTmist;RAW_missingTriggerType',
315  title='LVL1 trigger type missing or different from event type ;Slot;FT',
316  type='TH2I',
317  path=hist_path,
318  opt='kAlwaysCreate',
319  xbins=slot_n,xmin=slot_low,xmax=slot_up,
320  ybins=ft_n, ymin=ft_low, ymax=ft_up)
321 
322  darray.defineHistogram('slotnb,FTnb;RAW_nbOfEvts',
323  title='Nb of events (DSP header check only) ;Slot;FT',
324  type='TH2I',
325  path=hist_path,
326  opt='kAlwaysCreate',
327  xbins=slot_n,xmin=slot_low,xmax=slot_up,
328  ybins=ft_n, ymin=ft_low, ymax=ft_up)
329 
330  darray.defineHistogram('slotnb,FTnb,weightsweet1;RAW_NbOfSweet1PerFEB',
331  title='Average # of cells with (qfactor+time) readout ;Slot;FT',
332  type='TProfile2D',
333  path=hist_path,
334  opt='kAlwaysCreate',
335  xbins=slot_n,xmin=slot_low,xmax=slot_up,
336  ybins=ft_n, ymin=ft_low, ymax=ft_up)
337 
338  darray.defineHistogram('slotnb,FTnb,weightsweet2;RAW_NbOfSweet2PerFEB',
339  title='Average # of cells with samples readout ;Slot;FT',
340  type='TProfile2D',
341  path=hist_path,
342  opt='kAlwaysCreate',
343  xbins=slot_n,xmin=slot_low,xmax=slot_up,
344  ybins=ft_n, ymin=ft_low, ymax=ft_up)
345 
346  darray.defineHistogram('nbFEBpart;nbOfFebBlocks',
347  title='# of readout FEBs (DSP header check only) ;Slot;FT',
348  type='TH1I',
349  path=hist_path,
350  xbins=lArDQGlobals.N_FEB_Parttions_Max, xmin=-0.5, xmax=lArDQGlobals.N_FEB_Parttions_Max-0.5)
351 
352  darray.defineHistogram('slotMasked,FTMasked;RAW_knownFaultyFEB',
353  title='FEB with known errors (1:err. ignored 2:FEB masked);Slot;FT',
354  type='TH2I',
355  path=hist_path,
356  xbins=slot_n,xmin=slot_low,xmax=slot_up,
357  ybins=ft_n, ymin=ft_low, ymax=ft_up,
358  merge='identical')
359 
360  darray.defineHistogram('LB,LArEvSizePart;eventSizeVsLB',
361  title='LAr event size per LB (w/o ROS headers);Luminosity Block',
362  type='TProfile',
363  path=hist_path,
364  xbins=lArDQGlobals.LB_Bins, xmin=lArDQGlobals.LB_Min, xmax=lArDQGlobals.LB_Max)
365 
366  if isOnline:
367  darray.defineHistogram('LBf,erronl;RAW_EventsRejectedLB',
368  title='% of events rejected in current LB (online only)',
369  type='TProfile',
370  path=hist_path,
371  #opt='kAlwaysCreate',
372  opt='kLBNHistoryDepth=10,kAlwaysCreate',
373  xbins=1, xmin=0, xmax=1, xlabels=['% of events'])
374  darray.defineHistogram('LB,streamBin,LArEvSizePart;eventSizeStreamVsLB',
375  title='LAr event size per stream per LB (w/o ROS headers)',
376  type='TProfile2D',
377  path=hist_path,
378  xbins=lArDQGlobals.LB_Bins, xmin=lArDQGlobals.LB_Min, xmax=lArDQGlobals.LB_Max,
379  ybins=len(larFEBMonAlg.Streams),ymin=-0.5, ymax= len(larFEBMonAlg.Streams)-0.5,
380  ylabels=list(larFEBMonAlg.Streams)
381  )
382  pass
383 
384 
385 
386 
387 
python.JetAnalysisCommon.ComponentAccumulator
ComponentAccumulator
Definition: JetAnalysisCommon.py:302
LArFEBMonAlg.LArFEBMonConfigCore
def LArFEBMonConfigCore(helper, algoinstance, flags, cellDebug=False, dspDebug=False)
Definition: LArFEBMonAlg.py:22
LArFEBMonAlg.LArFEBMonConfig
def LArFEBMonConfig(flags, cellDebug=False, dspDebug=False)
Definition: LArFEBMonAlg.py:4
python.LArFebErrorSummaryMakerConfig.LArFebErrorSummaryMakerCfg
def LArFebErrorSummaryMakerCfg(flags)
Definition: LArFebErrorSummaryMakerConfig.py:7
plotBeamSpotVxVal.range
range
Definition: plotBeamSpotVxVal.py:195
histSizes.list
def list(name, path='/')
Definition: histSizes.py:38
python.IOVDbSvcConfig.addFolders
def addFolders(flags, folderStrings, detDb=None, className=None, extensible=False, tag=None, db=None, modifiers='')
Definition: IOVDbSvcConfig.py:81
python.IOVDbSvcConfig.IOVDbSvcCfg
def IOVDbSvcCfg(flags, **kwargs)
Definition: IOVDbSvcConfig.py:28
python.CaloAddPedShiftConfig.int
int
Definition: CaloAddPedShiftConfig.py:45
Generated on Mon May 5 2025 22:04:25 for ATLAS Offline Software by   doxygen 1.8.18