ATLAS Offline Software
Functions
TileTBMonitorAlgorithm Namespace Reference

Functions

def TileTBMonitoringConfig (flags, fragIDs=[0x100, 0x101, 0x200, 0x201, 0x402], **kwargs)
 

Function Documentation

◆ TileTBMonitoringConfig()

def TileTBMonitorAlgorithm.TileTBMonitoringConfig (   flags,
  fragIDs = [0x100, 0x101, 0x200, 0x201, 0x402],
**  kwargs 
) 
Function to configure TileTBMonitorAlgorithm algorithm in the monitoring system.

Definition at line 12 of file TileTBMonitorAlgorithm.py.

12 def TileTBMonitoringConfig(flags, fragIDs=[0x100, 0x101, 0x200, 0x201, 0x402], **kwargs):
13 
14  ''' Function to configure TileTBMonitorAlgorithm algorithm in the monitoring system.'''
15 
16  from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
17  result = ComponentAccumulator()
18 
19  from TileGeoModel.TileGMConfig import TileGMCfg
20  result.merge(TileGMCfg(flags))
21 
22  from LArGeoAlgsNV.LArGMConfig import LArGMCfg
23  result.merge(LArGMCfg(flags))
24 
25  from TileConditions.TileCablingSvcConfig import TileCablingSvcCfg
26  result.merge(TileCablingSvcCfg(flags))
27 
28  from AthenaMonitoring import AthMonitorCfgHelper
29  helper = AthMonitorCfgHelper(flags, 'TileTBMonitoring')
30 
31  from AthenaConfiguration.ComponentFactory import CompFactory
32  tileTBMonAlg = helper.addAlgorithm(CompFactory.TileTBMonitorAlgorithm, 'TileTBMonAlg')
33 
34  tileTBMonAlg.TriggerChain = ''
35 
36  kwargs.setdefault('BeamEnergy', flags.Beam.Energy)
37  kwargs.setdefault('CellEnergyThreshold', 0.1 * GeV)
38  energyThreshold = kwargs['CellEnergyThreshold']
39 
40  masked = ['LBA01 0 ' + ",".join([str(channel) for channel in range(0, 48)]),
41  'LBA01 1 ' + ",".join([str(channel) for channel in range(0, 48)]),
42  'LBC01 0 36,37,38,39,40,41',
43  'LBC01 1 36,37,38,39,40,41']
44  kwargs.setdefault('Masked', masked)
45 
46  from TileCalibBlobObjs.Classes import TileCalibUtils as Tile
47 
48  modules = []
49  if fragIDs:
50  for fragID in fragIDs:
51  ros = fragID >> 8
52  drawer = fragID & 0x3F
53  modules += [Tile.getDrawerString(ros, drawer)]
54  else:
55  for ros in range(1, Tile.MAX_ROS):
56  for drawer in range(0, Tile.MAX_DRAWER):
57  fragIDs += [(ros << 8) | drawer]
58  modules += [Tile.getDrawerString(ros, drawer)]
59 
60  tileTBMonAlg.TileFragIDs = fragIDs
61 
62  for k, v in kwargs.items():
63  setattr(tileTBMonAlg, k, v)
64 
65  run = str(flags.Input.RunNumbers[0])
66 
67  # Configure histogram with TileTBMonAlg algorithm execution time
68  executeTimeGroup = helper.addGroup(tileTBMonAlg, 'TileTBMonExecuteTime', 'TestBeam')
69  executeTimeGroup.defineHistogram('TIME_execute', path='', type='TH1F',
70  title='Time for execute TileTBMonAlg algorithm;time [#mus]',
71  xbins=100, xmin=0, xmax=10000)
72 
73  totalEnergyGroup = helper.addGroup(tileTBMonAlg, 'TileTBTotalEventEnergy', 'TestBeam')
74  totalEnergyGroup.defineHistogram('energy;TileTBTotalEventEnergy', path='', type='TH1F',
75  title=f'Run {run}: Total TileCal Event Energy;Event Energy [pC]',
76  xbins=400, xmin=-2, xmax=200)
77 
78  hotCellAGroup = helper.addGroup(tileTBMonAlg, 'TileTBHotCellA_LBC02', 'TestBeam')
79  hotCellAGroup.defineHistogram('tower;TileTBHotCellA_LBC02', path='', type='TH1F',
80  title=f'Run {run} LBC02: Tile TB Hot Cell A;Tower',
81  xbins=10, xmin=-0.5, xmax=9.5)
82 
83  CtotGroup = helper.addGroup(tileTBMonAlg, 'TileTBCtot', 'TestBeam')
84  CtotGroup.defineHistogram('Ctot;TileTBCtot', path='', type='TH1F',
85  title=f'Run {run} LBC02: Tile TB Ctot;Ctot',
86  xbins=100, xmin=0.03, xmax=0.22)
87 
88  ClongGroup = helper.addGroup(tileTBMonAlg, 'TileTBClong', 'TestBeam')
89  ClongGroup.defineHistogram('Clong;TileTBClong', path='', type='TH1F',
90  title=f'Run {run} LBC02: Tile TB Clong;Clong',
91  xbins=100, xmin=0., xmax=1.8)
92 
93  CtotVsClongGroup = helper.addGroup(tileTBMonAlg, 'TileTBCtotVsClong', 'TestBeam')
94  CtotVsClongGroup.defineHistogram('Clong,Ctot;TileTBCtotVsClong', path='', type='TH2F',
95  title=f'Run {run} LBC02: Tile TB Ctot Vs Clong;Clong;Ctot',
96  xbins=100, xmin=0., xmax=1.8, ybins=100, ymin=0.03, ymax=0.22)
97 
98 
99  maxTotalEnergy = 150
100  nCellsVsEnergyGroup = helper.addGroup(tileTBMonAlg, 'TileTBCellsNumberVsTotalEnergy', 'TestBeam')
101  nCellsVsEnergyGroup.defineHistogram('nCells,energy;TileTBCellsNumberVsTotalEnergy', path='', type='TH2F',
102  title=f'Run {run}: Tile Event energy [C side] vs # cells with energy > {energyThreshold} pC;# Cells;Energy [pC]',
103  xbins=25, xmin=-0.5, xmax=24.5, ybins=maxTotalEnergy, ymin=0.0, ymax=maxTotalEnergy)
104 
105  hitMapGroup = helper.addGroup(tileTBMonAlg, 'TileTBHitMap', 'TestBeam')
106  hitMapGroup.defineHistogram('side,module,energy;TileTBHitMap', path='', type='TProfile2D',
107  title=f'Run {run}: Tile TB setup map with average energy',
108  xlabels=['A Side', 'C side'], ylabels=['LB01', 'LB02', 'EB03'],
109  xbins=2, xmin=-0.5, xmax=1.5, ybins=3, ymin=-0.5, ymax=2.5)
110 
111  timeArray = helper.addArray([modules], tileTBMonAlg, 'TileTBChannelTime', topPath='TestBeam')
112  for postfix, tool in timeArray.Tools.items():
113  moduleName = postfix[1:]
114  name = f'channel,time;TileTBChannelTime_{moduleName}'
115  fullTitle = f'Run {run} {moduleName}: Tile TB channel average time from cells;Channel;Time [ns]'
116  tool.defineHistogram(name, title=fullTitle, path='', type='TProfile',
117  xbins=Tile.MAX_CHAN, xmin=-0.5, xmax=Tile.MAX_CHAN-0.5)
118 
119  maxSample = 3
120  maxTowerLB = 10
121  xCellLB = [[[] for tower in range(0, maxTowerLB)] for sample in range(0, maxSample)]
122  yCellLB = [[[] for tower in range(0, maxTowerLB)] for sample in range(0, maxSample)]
123 
124  periodWidthLB = 18.22
125  # Approximate number of periods in Tile Cells in LB per sampling and tower
126  nPeriodsLB = [[14, 13, 14, 14, 15, 16, 16, 17, 19, 16], # A1-A10
127  [16, 15, 16, 16, 17, 18, 18, 20, 18], # B1-B9
128  [18, 18, 18, 19, 19, 20, 22, 20], # C1-C8
129  [20, 0, 41, 0, 43, 0, 50]] # D0-D3
130 
131  yLB = [0, 300, 690, 1140, 1520]
132  for sampleIndex in range(0, len(nPeriodsLB)):
133  sample = sampleIndex if sampleIndex < 2 else sampleIndex - 1
134  cellOffsetX = 0.0
135  for tower in range(0, len(nPeriodsLB[sampleIndex])):
136  cellWidth = nPeriodsLB[sampleIndex][tower] * periodWidthLB
137  x1 = cellOffsetX + 9 # Approximate center of the period
138  x2 = cellOffsetX + cellWidth
139  while x1 < x2:
140  xCellLB[sample][tower] += [x1]
141  yCellLB[sample][tower] += [yLB[sampleIndex]]
142  x1 += periodWidthLB
143 
144  cellOffsetX += cellWidth
145 
146  tileTBMonAlg.xCellLongBarrelSampleA = xCellLB[0]
147  tileTBMonAlg.xCellLongBarrelSampleBC = xCellLB[1]
148  tileTBMonAlg.xCellLongBarrelSampleD = xCellLB[2]
149 
150  tileTBMonAlg.yCellLongBarrelSampleA = yCellLB[0]
151  tileTBMonAlg.yCellLongBarrelSampleBC = yCellLB[1]
152  tileTBMonAlg.yCellLongBarrelSampleD = yCellLB[2]
153 
154  maxTowerEB = 16
155  periodWidthEB = 18.28
156  xCellEB = [[[] for tower in range(0, maxTowerEB)] for sample in range(0, maxSample)]
157  yCellEB = [[[] for tower in range(0, maxTowerEB)] for sample in range(0, maxSample)]
158 
159  yEB = [0, 300, 690, 840, 840, 1140, 1520]
160  # Approximate number of periods in Tile Cells in EB per sampling and tower
161  nPeriodsEB = [[17, 0, 0, 9, 25, 28, 30, 48], # A12-A16
162  [17, 0, 16, 27, 30, 32, 35], # B11-B15
163  [12, 5, 16, 27, 30, 32, 35], # C10, B11-B15
164  [12, 5], # C10
165  [0, 17, 65, 0, 75], # D5,D6
166  [17, 0, 65, 0, 75]] # D4,D5,D6
167 
168  samples = [0, 1, 1, 1, 2, 2]
169  for sampleIndex in range(0, len(nPeriodsEB)):
170  sample = samples[sampleIndex]
171  cellOffsetX = 0.0
172  for tower in range(0, len(nPeriodsEB[sampleIndex])):
173  cellWidth = nPeriodsEB[sampleIndex][tower] * periodWidthEB
174  x1 = cellOffsetX + 9 # Approximate center of the period
175  x2 = cellOffsetX + cellWidth
176  while x1 < x2:
177  xCellEB[sample][tower + 8] += [x1]
178  yCellEB[sample][tower + 8] += [yEB[sampleIndex]]
179  x1 += periodWidthEB
180 
181  cellOffsetX += cellWidth
182 
183  tileTBMonAlg.xCellExtendedBarrelSampleA = xCellEB[0]
184  tileTBMonAlg.xCellExtendedBarrelSampleBC = xCellEB[1]
185  tileTBMonAlg.xCellExtendedBarrelSampleD = xCellEB[2]
186 
187  tileTBMonAlg.yCellExtendedBarrelSampleA = yCellEB[0]
188  tileTBMonAlg.yCellExtendedBarrelSampleBC = yCellEB[1]
189  tileTBMonAlg.yCellExtendedBarrelSampleD = yCellEB[2]
190 
191  xBinsLB = [x * periodWidthLB for x in range(0, sum(nPeriodsLB[0]) + 1)]
192  yBinsLB = [0, 300, 690, 1140, 1520]
193 
194  xBinsEB = [x * periodWidthEB for x in range(0, sum(nPeriodsEB[0]) + 1)]
195  yBinsEB = [0, 300, 690, 840, 1140, 1520]
196  cellMapArray = helper.addArray([modules], tileTBMonAlg, 'TileTBCellMap', topPath='TestBeam')
197  for postfix, tool in cellMapArray.Tools.items():
198  moduleName = postfix[1:]
199  name = f'x,y,energy;TileTBCellMap_{moduleName}'
200  fullTitle = f'Run {run} {moduleName}: Tile TB cell map with average energy;z [mm];x [mm];Energy [pc]'
201  tool.defineHistogram(name, title=fullTitle, path='', type='TProfile2D',
202  xbins=xBinsLB if moduleName.startswith('LB') else xBinsEB,
203  ybins=yBinsLB if moduleName.startswith('LB') else yBinsEB)
204 
205  accumalator = helper.result()
206  result.merge(accumalator)
207  return result
208 
209 
python.JetAnalysisCommon.ComponentAccumulator
ComponentAccumulator
Definition: JetAnalysisCommon.py:302
convertTimingResiduals.sum
sum
Definition: convertTimingResiduals.py:55
TileTBMonitorAlgorithm.TileTBMonitoringConfig
def TileTBMonitoringConfig(flags, fragIDs=[0x100, 0x101, 0x200, 0x201, 0x402], **kwargs)
Definition: TileTBMonitorAlgorithm.py:12
plotBeamSpotVxVal.range
range
Definition: plotBeamSpotVxVal.py:195
LArGMConfig.LArGMCfg
def LArGMCfg(flags)
Definition: LArGMConfig.py:8
TCS::join
std::string join(const std::vector< std::string > &v, const char c=',')
Definition: Trigger/TrigT1/L1Topo/L1TopoCommon/Root/StringUtils.cxx:10
str
Definition: BTagTrackIpAccessor.cxx:11
python.TileCablingSvcConfig.TileCablingSvcCfg
def TileCablingSvcCfg(flags)
Definition: TileCablingSvcConfig.py:11
TileGMConfig.TileGMCfg
def TileGMCfg(flags)
Definition: TileGMConfig.py:7
Generated on Sun Dec 22 2024 22:09:03 for ATLAS Offline Software by   doxygen 1.8.18