ATLAS Offline Software
Loading...
Searching...
No Matches
Functions
TileCellMonitorAlgorithm Namespace Reference

Functions

 TileCellMonitoringConfig (flags, **kwargs)

Detailed Description

@file TileCellMonitorAlgorithm.py
@brief Python configuration of TileCellMonitorAlgorithm algorithm for the Run III

Function Documentation

◆ TileCellMonitoringConfig()

TileCellMonitorAlgorithm.TileCellMonitoringConfig ( flags,
** kwargs ) 
Function to configure TileCellMonitorAlgorithm algorithm in the monitoring system.

Definition at line 8 of file TileCellMonitorAlgorithm.py.

8def TileCellMonitoringConfig(flags, **kwargs):
9
10 ''' Function to configure TileCellMonitorAlgorithm algorithm in the monitoring system.'''
11
12 # Define one top-level monitoring algorithm. The new configuration
13 # framework uses a component accumulator.
14 from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
15 result = ComponentAccumulator()
16
17 from TileRecUtils.TileDQstatusConfig import TileDQstatusAlgCfg
18 result.merge( TileDQstatusAlgCfg(flags) )
19
20 from TileGeoModel.TileGMConfig import TileGMCfg
21 result.merge(TileGMCfg(flags))
22
23 from LArGeoAlgsNV.LArGMConfig import LArGMCfg
24 result.merge(LArGMCfg(flags))
25
26 from TileConditions.TileCablingSvcConfig import TileCablingSvcCfg
27 result.merge( TileCablingSvcCfg(flags) )
28
29 from TileConditions.TileBadChannelsConfig import TileBadChannelsCondAlgCfg
30 result.merge( TileBadChannelsCondAlgCfg(flags) )
31
32 # The following class will make a sequence, configure algorithms, and link
33 # them to GenericMonitoringTools
34 from AthenaMonitoring import AthMonitorCfgHelper
35 helper = AthMonitorCfgHelper(flags,'TileMonitoring')
36
37 # Adding an TileCellMonitorAlgorithm algorithm to the helper
38 from AthenaConfiguration.ComponentFactory import CompFactory
39 tileCellMonAlg = helper.addAlgorithm(CompFactory.TileCellMonitorAlgorithm, 'TileCellMonAlg')
40
41 tileCellMonAlg.TriggerChain = ''
42
43 from AthenaCommon.SystemOfUnits import MeV, GeV, ns
44 kwargs.setdefault('EnergyThreshold', 300.0 * MeV)
45 kwargs.setdefault('EnergyThresholdForGain', [10 * GeV, 300.0 * MeV])
46 kwargs.setdefault('NegativeEnergyThreshold', -2000.0 * MeV)
47 kwargs.setdefault('EnergyBalanceThreshold', 3)
48 kwargs.setdefault('TimeBalanceThreshold', 25 * ns)
49 kwargs.setdefault('EnergyThresholdForGapScintilator', 0 * MeV)
50 kwargs.setdefault('fillChannelTimeHistograms', False)
51 kwargs.setdefault('fillTimeAndEnergyDiffHistograms', flags.Common.isOnline)
52 kwargs.setdefault('fillGapScintilatorHistograms', False)
53 kwargs.setdefault('EnergyLimitForTime', 750.0 * GeV)
54
55 from AthenaConfiguration.Enums import BeamType
56 if flags.Beam.Type in [BeamType.Cosmics, BeamType.SingleBeam]:
57 kwargs.setdefault('fillTimeHistograms', True)
58 kwargs.setdefault('EnergyThresholdForTime', 150.0 * MeV)
59 else:
60 kwargs.setdefault('fillTimeHistograms', False)
61 kwargs.setdefault('EnergyThresholdForTime', 500.0 * MeV)
62
63 # L1Trigger Type Bits:
64 # bit0_RNDM, bit1_ZeroBias, bit2_L1Cal, bit3_Muon,
65 # bit4_RPC, bit5_FTK, bit6_CTP, bit7_Calib, AnyPhysTrig
66 kwargs.setdefault('fillHistogramsForL1Triggers', ['AnyPhysTrig', 'bit7_Calib'])
67 l1Triggers = kwargs['fillHistogramsForL1Triggers']
68
69 for k, v in kwargs.items():
70 setattr(tileCellMonAlg, k, v)
71
72 run = str(flags.Input.RunNumbers[0])
73
74 # 1) Configure histogram with TileCellMonAlg algorithm execution time
75 executeTimeGroup = helper.addGroup(tileCellMonAlg, 'TileCellMonExecuteTime', 'Tile/')
76 executeTimeGroup.defineHistogram('TIME_execute', path = 'Cell', type='TH1F',
77 title = 'Time for execute TileCellMonAlg algorithm;time [#mus]',
78 xbins = 100, xmin = 0, xmax = 100000)
79
80
81 from TileCalibBlobObjs.Classes import TileCalibUtils as Tile
82 from TileMonitoring.TileMonitoringCfgHelper import getPartitionName, getGainName
83
84 # 2) Configure histogram with average number of Tile bad cells in partition
85 labelsPartitions = [getPartitionName(ros) for ros in range(1, Tile.MAX_ROS)]
86 badCellGroup = helper.addGroup(tileCellMonAlg, 'TileBadCell', 'Tile/')
87 badCellGroup.defineHistogram('Partition,nBadCells;TileBadCell', path = 'Cell', type='TProfile',
88 xlabels = labelsPartitions, title = ('Run %s: Average number of Bad Tile Cells' % run),
89 xbins = Tile.MAX_ROS - 1, xmin = -0.5, xmax = Tile.MAX_ROS - 1.5)
90
91
92 from TileMonitoring.TileMonitoringCfgHelper import addTileModulePartitionMapsArray
93
94 # 3) Configure histograms with energy ratio > threshold vs module and partion
95 eneBalModPartTitle = ('Tile Cell Energy Ratio > %s' % kwargs['EnergyBalanceThreshold'])
96 addTileModulePartitionMapsArray(helper, tileCellMonAlg, name = 'TileCellEneBalModPart',
97 title = eneBalModPartTitle, path = 'Tile/Cell', type='TH2D',
98 run = run, triggers = l1Triggers, separator = '_')
99
100 # 4) Configure histograms with Tile cell time difference over threshold vs module and partition
101 timeBalModPartTitle = ('Tile Cell Time difference > %s ns. E_{ch} > %s [MeV]')
102 timeBalModPartTitle = (timeBalModPartTitle % (kwargs['TimeBalanceThreshold'], kwargs['EnergyThresholdForTime']))
103 addTileModulePartitionMapsArray(helper, tileCellMonAlg, name = 'TileCellTimeBalModPart',
104 title = timeBalModPartTitle, path = 'Tile/Cell', type='TH2D',
105 run = run, triggers = l1Triggers, separator = '_')
106
107
108 from TileMonitoring.TileMonitoringCfgHelper import addTileModuleChannelMapsArray, addTileModuleDigitizerMapsArray
109
110 # 5) Configure histograms with everagy Tile channel time per partition
111 chanTimeTitle = ('TileCal Average Channel Time [ns]. E_{ch} > %s MeV' % kwargs['EnergyThresholdForTime'])
112 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileChanPartTime', type='TProfile2D',
113 value = 'time', title = chanTimeTitle, path = 'Tile/Cell',
114 subDirectory = True, run = run, triggers = l1Triggers, separator = '_')
115
116 # 6) Configure histograms with everagy Tile digitizer time per partition
117 digiTimeTitle = ('TileCal Average Digitizer Time [ns]. E_{ch} > %s MeV' % kwargs['EnergyThresholdForTime'])
118 addTileModuleDigitizerMapsArray(helper, tileCellMonAlg, name = 'TileDigiPartTime', type='TProfile2D',
119 value = 'time', title = digiTimeTitle, path = 'Tile/Cell',
120 subDirectory = True, run = run, triggers = l1Triggers, separator = '_')
121
122
123 # 7) Configure histograms with average energy maps per partition
124 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileCellDetailOccMap', type='TProfile2D',
125 value = 'energy', title = 'Occupancy Map [MeV]', path = 'Tile/Cell',
126 subDirectory = True, run = run, triggers = l1Triggers, separator = '_')
127
128
129 # 7.2) Configure histograms with average energy maps per partition and gain
130 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileCellDetailOccMapGain', type='TProfile2D',
131 value = 'energy', title = 'Occupancy Map [MeV]', path = 'Tile/Cell',
132 subDirectory = True, run = run, triggers = l1Triggers, perGain = True, separator = '_')
133
134
135 # 8) Configure histograms with occupancy maps over threshold per partition
136 titleMapOvThr = ('Occupancy Map Over Threshod %s MeV' % kwargs['EnergyThreshold'])
137 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileCellDetailOccMapOvThr',
138 weight = 'weight', title = titleMapOvThr, path = 'Tile/Cell',
139 subDirectory = True, run = run, triggers = l1Triggers, separator = '_')
140
141
142 # 9) Configure histograms with occupancy maps over threshold 30GeV per partition
143 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileCellDetailOccMapOvThr30GeV',
144 title = 'Occupancy Map Over Threshod 30 GeV', path = 'Tile/Cell',
145 subDirectory = True, run = run, triggers = l1Triggers, separator = '_')
146
147
148 # 10) Configure histograms with occupancy maps over threshold 300GeV per partition
149 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileCellDetailOccMapOvThr300GeV',
150 title = 'Occupancy Map Over Threshod 300 GeV', path = 'Tile/Cell',
151 subDirectory = True, run = run, triggers = l1Triggers, separator = '_')
152
153 # ) Configure histograms with Tile cell energy difference maps per partition
154 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileCellEneDiffChanMod',
155 title = 'Tile Cell energy difference between PMTs [MeV]', path = 'Tile/Cell',
156 subDirectory = True, type = 'TProfile2D', value = 'energyDiff',
157 run = run, triggers = l1Triggers, separator = '_')
158
159 # 11) Configure histograms with occupancy maps over threshold per partition
160 titleMapOvThrGain = {}
161 for gain in range(0, Tile.MAX_GAIN):
162 gainName = getGainName(gain)
163 energyThresholdForGain = kwargs['EnergyThresholdForGain'][gain]
164 energyThreshold = f'{energyThresholdForGain} MeV' if energyThresholdForGain < 1000.0 else f'{energyThresholdForGain/GeV} GeV'
165 titleMapOvThrGain[gainName] = f'Occupancy Map Over Threshod {energyThreshold}'
166 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileCellDetailOccMapOvThrGain',
167 weight = 'weight', title = titleMapOvThrGain, path = 'Tile/Cell', subDirectory = True,
168 run = run, triggers = l1Triggers, perGain = True, separator = '_')
169
170
171 # 12) Configure histograms with status of Tile channels in DB per partition
172 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileCellStatusInDB',
173 title = 'Cell channel status in DB', path = 'Tile/Cell',
174 run = run, perGain = True, separator = '_')
175
176 # 13) Configure histograms with Tile channels masked on the fly per partition
177 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileCellStatusOnFly',
178 title = 'Channels masked on the fly', path = 'Tile/Cell',
179 run = run, perGain = False, separator = '_')
180
181
182 # 14) Configure histograms with occupancy map below negative threshold per partition
183 titleNegOcc = 'Occupancy Map Below %s GeV' % (kwargs['NegativeEnergyThreshold'] / GeV)
184 addTileModuleChannelMapsArray(helper, tileCellMonAlg, name = 'TileCellDetailNegOccMap',
185 title = titleNegOcc, path = 'Tile/Cell', run = run, separator = '_')
186
187 # 15) Configure histograms with Tile module correlation per partition
188 from TileMonitoring.TileMonitoringCfgHelper import addTileModuleCorrelionMapsArray
189 addTileModuleCorrelionMapsArray(helper, tileCellMonAlg, name = 'TileCellModuleCorrelation',
190 title = 'Tile Cell Module correlation', path = 'Tile/Cell',
191 weight = 'weight', subDirectory = True, run = run,
192 triggers = l1Triggers, allPartitions = True, separator = '_')
193
194
195 from TileMonitoring.TileMonitoringCfgHelper import addTile1DHistogramsArray
196
197 # 16) Configure histograms with number of masked Tile channels on the fly vs lumi block per partition
198 titleMaskOnFlyLB = 'Number of masked channels on the fly'
199 titleMaskOnFlyLB += ';LumiBlock;Number of masked channels'
200 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileMaskChannelOnFlyLB', path = 'Tile/Cell',
201 xvalue = 'lumiBlock', value= 'nMaskedChannelsOnFly', title = titleMaskOnFlyLB,
202 xbins = 1000, xmin = -0.5, xmax = 999.5, type = 'TProfile', run = run, triggers = [],
203 perPartition = True, perSample = False, perGain = False, subDirectory = False,
204 opt = 'kAddBinsDynamically', merge = 'merge', allPartitions = True)
205
206 # 17) Configure histograms with number of masked Tile cells on the fly vs lumi block per partition
207 titleMaskCellLB = 'Number of masked cells on the fly'
208 titleMaskCellLB += ';LumiBlock;Number of masked cells'
209 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileMaskCellLB', path = 'Tile/Cell',
210 xvalue = 'lumiBlock', value = 'nMaskedCells', title = titleMaskCellLB,
211 xbins = 1000, xmin = -0.5, xmax = 999.5, type='TProfile', run = run, triggers = [],
212 subDirectory = False, perPartition = True, perSample = False, perGain = False,
213 opt = 'kAddBinsDynamically', merge = 'merge', allPartitions = True)
214
215 # 18) Configure histograms with number of masked Tile channels on the fly due to bad DQ status vs lumi block per partition
216 titleMaskDueDQ = 'Number of masked channels on the fly due to bad DQ status'
217 titleMaskDueDQ += ';LumiBlock;Number of masked channels'
218 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileMaskChannelDueDQvsLB', path = 'Tile/Cell',
219 xvalue = 'lumiBlock', value = 'nMaskedChannelsDueDQ', title = titleMaskDueDQ,
220 xbins = 1000, xmin = -0.5, xmax = 999.5, type='TProfile', run = run, triggers = [],
221 subDirectory = False, perPartition = True, perSample = False, perGain = False,
222 opt = 'kAddBinsDynamically', merge = 'merge', allPartitions = True)
223
224 # 19) Configure histograms with number of masked Tile cells on the fly due to bad DQ status vs lumi block per partition
225 titleMaskCellDueDQ = 'Number of masked cells on the fly due to bad DQ status'
226 titleMaskCellDueDQ += ';LumiBlock;Number of masked cells'
227 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileMaskedCellDueDQvsLB', path = 'Tile/Cell',
228 xvalue = 'lumiBlock', value = 'nMaskedCellsDueDQ', title = titleMaskCellDueDQ,
229 xbins = 1000, xmin = -0.5, xmax = 999.5, type='TProfile', run = run, triggers = [],
230 subDirectory = False, perPartition = True, perSample = False, perGain = False,
231 opt = 'kAddBinsDynamically', merge = 'merge', allPartitions = True)
232
233 if kwargs['fillChannelTimeHistograms']:
234 # Configure histograms with Tile channel time per partition and sample
235 titleChanTimeSamp = 'Channel Time, E_{ch} > %s MeV;time [ns]' % (kwargs['EnergyThresholdForTime'] / MeV)
236 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileChannelTime',
237 xvalue = 'time', title = titleChanTimeSamp, path = 'Tile/Cell',
238 xbins = 121, xmin = -60.5, xmax = 60.5, type='TH1D',
239 run = run, triggers = l1Triggers, subDirectory = True,
240 perPartition = True, perSample = True, perGain = False)
241
242 # 20) Configure histograms with energy difference between Tile cells' PMTs per partition and sample
243 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileCellEneDiff', xvalue = 'energyDiff',
244 title = 'Energy difference [MeV] between PMTs;Energy difference [MeV]',
245 path = 'Tile/Cell', xbins = 50, xmin = -1000., xmax = 1000., type='TH1D',
246 run = run, triggers = l1Triggers, subDirectory = True,
247 perPartition = True, perSample = True, perGain = False)
248
249 # 21) Configure histograms with time difference between Tile cells' PMTs per partition and sample
250 titleTimeDiffSamp = 'Time difference [ns] between PMTs with '
251 titleTimeDiffSamp += 'E_{ch} > %s MeV' % (kwargs['EnergyThresholdForTime'] / MeV)
252 titleTimeDiffSamp += ';time [ns]'
253 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileCellTimeDiff',
254 xvalue = 'timeDiff', title = titleTimeDiffSamp, path = 'Tile/Cell',
255 xbins = 50, xmin = -10., xmax = 10., type='TH1D', run = run, triggers = l1Triggers,
256 subDirectory = True, perPartition = True, perSample = True, perGain = False)
257
258 # Configure histograms with number of Tile cells vs lumiBlock per partition
259 titleCellsNumber = 'Tile Cells number per luminosity block;LumiBlock;Number of reconstructed cells'
260 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileCellsNumberLB', opt = 'kAddBinsDynamically', merge = 'merge',
261 xvalue = 'lumiBlock', value = 'nCells', title = titleCellsNumber, path = 'Tile/Cell',
262 xbins = 1000, xmin = -0.5, xmax = 999.5, type='TProfile', run = run, triggers = l1Triggers,
263 subDirectory = True, perPartition = True, perSample = False, perGain = False, allPartitions = True)
264
265 # 22) Configure histograms with number of Tile cells over threshold vs BCID per partition
266 titleCellsOvThrBCID = 'Tile Cell Occupancy over Threshold %s MeV' % (kwargs['EnergyThresholdForTime'] / MeV)
267 titleCellsOvThrBCID += ';BCID;Average number of cells over threshold'
268 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileCellOccOvThrBCID',
269 xvalue = 'BCID', value = 'nCells', title = titleCellsOvThrBCID, path = 'Tile/Cell',
270 xbins = 3565, xmin = 0., xmax = 3565., type='TProfile', run = run, triggers = l1Triggers,
271 subDirectory = True, perPartition = True, perSample = False, perGain = False, allPartitions = True)
272
273 # 23) Configure histograms with number of Tile E cell's energy per partition
274 titleEvEnergy = 'Tile Event SampE Energy;Event Energy [MeV]'
275 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileCellEventEnergy_SampE',
276 xvalue = 'energy', title = titleEvEnergy, path = 'Tile/Cell',
277 xbins = 120, xmin = -2000., xmax = 10000., type='TH1D', run = run, triggers = l1Triggers,
278 subDirectory = True, perPartition = True, perSample = False, perGain = False, allPartitions = True)
279
280 # 24) Configure histograms with number of Tile E cell's energy
281 titleSynch = 'Tile Time of Flight - Tile measured;Time of Flight - Tile measured [ns]'
282 addTile1DHistogramsArray(helper, tileCellMonAlg, name = 'TileCellSynch',
283 xvalue = 'timeDifference', title = titleSynch, path = 'Tile/Cell',
284 xbins = 50, xmin = -100., xmax = 100., type='TH1D', run = run, triggers = l1Triggers,
285 subDirectory = False, perPartition = False, perSample = False, perGain = False)
286
287
288 from TileMonitoring.TileMonitoringCfgHelper import addTileEtaPhiMapsArray
289
290 # 25) Configure histograms with everage Tile cell energy vs eta and phy per sample
291 addTileEtaPhiMapsArray(helper, tileCellMonAlg, name = 'TileCellEneEtaPhi', type='TProfile2D',
292 value = 'energy', title = 'Energy Average depostion [MeV]',
293 path = 'Tile/Cell', run = run, triggers = l1Triggers, perSample = True)
294
295 # 26) Configure histograms with number of Tile cells over threshold vs eta and phy per sample
296 titleEtaPhiOvThr = ('Position of cells over threshold %s MeV' % kwargs['EnergyThreshold'])
297 addTileEtaPhiMapsArray(helper, tileCellMonAlg, name = 'TileCellEtaPhiOvThr',
298 type='TH2D', title = titleEtaPhiOvThr, path = 'Tile/Cell',
299 run = run, triggers = l1Triggers, perSample = True)
300
301
302 from TileMonitoring.TileMonitoringCfgHelper import addTileModuleArray
303
304 # 27) Configure histograms with energy difference between Tile cell's PMTs vs module per partition
305 titleEnergyBal = "Cell's PMTs Energy Balance"
306 titleEnergyBal += ";;Energy balance between cell's PMTs (u-d)/(u+d)"
307 addTileModuleArray(helper, tileCellMonAlg, name = 'TileCellEnergyBalance', type='TProfile',
308 title = titleEnergyBal, path = 'Tile/Cell', value = 'energyBalance', run = run)
309
310
311 # 28) Configure histograms with time difference between Tile cell's PMTs vs module per partition
312 titleTimeBal = "Cell's PMTs Time Difference with "
313 titleTimeBal += 'E_{ch} > %s MeV' % (kwargs['EnergyThresholdForTime'] / MeV)
314 titleTimeBal += ";;Time balance between cell's PMTs [ns]"
315 addTileModuleArray(helper, tileCellMonAlg, name = 'TileCellTimeBalance', type='TProfile',
316 title = titleTimeBal, path = 'Tile/Cell', value = 'timeBalance', run = run)
317
318
319 if kwargs['fillGapScintilatorHistograms']:
320 rangeScalFactor = [1, 1, 1, 1]
321 gapScintDimentions = [2, Tile.MAX_DRAWER, 4] # EBA and EBC, 64 modules, E1-E4
322 gapScintArray = helper.addArray(gapScintDimentions, tileCellMonAlg, 'TileGapScintilatorEnergy', topPath = 'Tile/Cell')
323 for postfix, tool in gapScintArray.Tools.items():
324 ros, module, cellIdx = [int(x) for x in postfix.split('_')[1:]]
325 ros += 3
326
327 partionName = getPartitionName(ros)
328 moduleName = Tile.getDrawerString(ros, module)
329
330 fullName = f'energy;TileGapScintilatorEnergy{moduleName}_E{cellIdx + 1}'
331 fullTitle = f'Run {run} {moduleName} E{cellIdx + 1}: Energy;Energy [MeV]'
332
333 tool.defineHistogram(fullName, title = fullTitle, path = f'GapScint/{partionName}', type = 'TH1F',
334 xbins = 1000, xmin = -1.5, xmax = 25000.5 * rangeScalFactor[cellIdx])
335
336
337 accumalator = helper.result()
338 result.merge(accumalator)
339 return result
340
Generated on for ATLAS Offline Software by doxygen 1.14.0