ATLAS Offline Software
Functions
TileTBBeamMonitorAlgorithm Namespace Reference

Functions

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

Function Documentation

◆ TileTBBeamMonitoringConfig()

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

Definition at line 12 of file TileTBBeamMonitorAlgorithm.py.

12 def TileTBBeamMonitoringConfig(flags, fragIDs=[0x100,0x101,0x200,0x201,0x402], **kwargs):
13 
14  ''' Function to configure TileTBBeamMonitorAlgorithm 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 TileConditions.TileCablingSvcConfig import TileCablingSvcCfg
23  result.merge(TileCablingSvcCfg(flags))
24 
25  from AthenaMonitoring import AthMonitorCfgHelper
26  helper = AthMonitorCfgHelper(flags, 'TileTBBeamMonitoring')
27 
28  from AthenaConfiguration.ComponentFactory import CompFactory
29  tileTBBeamMonAlg = helper.addAlgorithm(CompFactory.TileTBBeamMonitorAlgorithm, 'TileTBBeamMonAlg')
30 
31  tileTBBeamMonAlg.TriggerChain = ''
32 
33  kwargs.setdefault('CaloCellContainer', 'AllCalo')
34  cellContainer = kwargs['CaloCellContainer']
35 
36  kwargs.setdefault('TBperiod', 2023)
37  kwargs.setdefault('MaskMuonPMTs', [7])
38  kwargs.setdefault('TOFDifferencePairs', [[2, 1], [2, 3], [3, 1]])
39 
40  from TileCalibBlobObjs.Classes import TileCalibUtils as Tile
41 
42  modules = []
43  if fragIDs:
44  for fragID in fragIDs:
45  ros = fragID >> 8
46  drawer = fragID & 0x3F
47  modules += [Tile.getDrawerString(ros, drawer)]
48  else:
49  for ros in range(1, Tile.MAX_ROS):
50  for drawer in range(0, Tile.MAX_DRAWER):
51  fragIDs += [(ros << 8) | drawer]
52  modules += [Tile.getDrawerString(ros, drawer)]
53 
54  tileTBBeamMonAlg.TileFragIDs = fragIDs
55 
56  for k, v in kwargs.items():
57  setattr(tileTBBeamMonAlg, k, v)
58 
59  from TileMonitoring.TileTBBeamChambersCalibration import updateBeamChambersCalibrations
60  updateBeamChambersCalibrations(tileTBBeamMonAlg, flags.Input.RunNumbers[0])
61 
62  run = str(flags.Input.RunNumbers[0])
63 
64  # Configure histogram with TileTBBeamMonAlg algorithm execution time
65  executeTimeGroup = helper.addGroup(tileTBBeamMonAlg, 'TileTBBeamMonExecuteTime', 'TestBeam')
66  executeTimeGroup.defineHistogram('TIME_execute', path='BeamElements', type='TH1F',
67  title='Time for execute TileTBBeamMonAlg algorithm;time [#mus]',
68  xbins=100, xmin=0, xmax=10000)
69 
70  nMuonWallPMT = 12
71  muonWallPMTArray = helper.addArray([nMuonWallPMT], tileTBBeamMonAlg, 'MuonWallPMT', topPath='TestBeam')
72  for postfix, tool in muonWallPMTArray.Tools.items():
73  pmt = int(postfix[1:]) + 1
74  title = f'Run {run}: Muon Wall PMT{pmt} Amplitude;Amplitude [ADC];Counts'
75  name = f'amplitude;MuonWallPMT{pmt}'
76  tool.defineHistogram(name, title=title, path='BeamElements', type='TH1F',
77  xbins=410, xmin=0, xmax=4096)
78 
79  totalMuEnergyGroup = helper.addGroup(tileTBBeamMonAlg, 'TileTBTotalMuonEnergy', 'TestBeam')
80  totalMuEnergyGroup.defineHistogram('TotalMuonEnergy', path='BeamElements', type='TH1F',
81  title=f'Run {run}: Muon Wall Total Energy; [ADC]',
82  xbins=1500, xmin=0, xmax=10000)
83 
84  nScounters = 3
85  sCountersArray = helper.addArray([nScounters], tileTBBeamMonAlg, 'Scounter', topPath='TestBeam')
86  for postfix, tool in sCountersArray.Tools.items():
87  counter = int(postfix[1:]) + 1
88  title = f'Run {run}: S{counter} Counter Amplitude;Amplitude [ADC];Counts'
89  name = f'amplitude;S{counter}hist'
90  tool.defineHistogram(name, title=title, path='BeamElements', type='TH1F',
91  xbins=410, xmin=0, xmax=4096)
92 
93  nCherenkov = 3
94  cherenkovArray = helper.addArray([nCherenkov], tileTBBeamMonAlg, 'Cherenkov', topPath='TestBeam')
95  for postfix, tool in cherenkovArray.Tools.items():
96  cherenkov = int(postfix[1:]) + 1
97  title = f'Run {run}: Cherenkov {cherenkov} Amplitude;Amplitude [ADC];Counts'
98  name = f'amplitude;Cher{cherenkov}hist'
99  tool.defineHistogram(name, title=title, path='BeamElements', type='TH1F',
100  xbins=410, xmin=0, xmax=4096)
101 
102 
103  cherCompGroup = helper.addGroup(tileTBBeamMonAlg, 'CherCompare', 'TestBeam')
104  cherCompGroup.defineHistogram('amplitude1,amplitude2;CherCompare', path='BeamElements', type='TH2F',
105  title=f'Run {run}: Cherenkov2 vs Cherenkov1;Amplitude [ADC];Amplitude [ADC]',
106  xbins=410, xmin=0, xmax=4096, ybins=410, ymin=0, ymax=4096)
107 
108  nTOF = 3
109  tofArray = helper.addArray([nTOF], tileTBBeamMonAlg, 'TOF', topPath='TestBeam')
110  for postfix, tool in tofArray.Tools.items():
111  tof = int(postfix[1:]) + 1
112  title = f'Run {run}: TOF{tof};[ADC];Counts'
113  name = f'amplitude;TOF{tof}'
114  tool.defineHistogram(name, title=title, path='BeamElements', type='TH1F',
115  xbins=4096, xmin=-0.5, xmax=4095.5)
116 
117  tofDiffPairs = kwargs['TOFDifferencePairs']
118  tofDiffArray = helper.addArray([len(tofDiffPairs)], tileTBBeamMonAlg, 'TOFDiff', topPath='TestBeam')
119  for postfix, tool in tofDiffArray.Tools.items():
120  pairIdx = int(postfix.split('_').pop())
121  tof1,tof2 = tofDiffPairs[pairIdx]
122  title = f'Run {run}: TOF{tof1} - TOF{tof2} Amplitude difference;[ADC];Counts'
123  tool.defineHistogram(f'TOFDiff;TOFDiff{tof1}{tof2}', title=title, path='BeamElements',
124  type='TH1F', xbins=4096, xmin=-0.5, xmax=4095.5)
125 
126  cherenkovVsTOFArray = helper.addArray([nTOF, nCherenkov], tileTBBeamMonAlg, 'CherenkovVsTOF', topPath='TestBeam')
127  for postfix, tool in cherenkovVsTOFArray.Tools.items():
128  cherenkovTof = postfix.split('_')
129  cherenkov = int(cherenkovTof.pop()) + 1
130  tof = int(cherenkovTof.pop()) + 1
131  title = f'Run {run}: Cherenkov {cherenkov} Amplitude vs TOF{tof} Amplitude'
132  title += f';TOF{tof} Amplitude [ADC];Cherenkov{cherenkov} Amplitude [ADC]'
133  name = f'amplitudeTOF,amplitudeCherenkov;Cher{cherenkov}TOF{tof}'
134  tool.defineHistogram(name, title=title, path='BeamElements', type='TH2F',
135  xbins=410, xmin=0, xmax=4096, ybins=410, ymin=0, ymax=4096)
136 
137  pmtHitMapGroup = helper.addGroup(tileTBBeamMonAlg, 'PMTHitMap', 'TestBeam')
138  pmtHitMapGroup.defineHistogram('column,row,amplitude;PMTHitMap', path='BeamElements', type='TProfile2D',
139  title=f'Run {run}: Muon Wall PMT Hit Map',
140  xbins=4, xmin=0, xmax=4, ybins=2, ymin=0, ymax=2)
141 
142  nScaler = 3
143  scalerArray = helper.addArray([nScaler], tileTBBeamMonAlg, 'Scaler', topPath='TestBeam')
144  for postfix, tool in scalerArray.Tools.items():
145  scaler = int(postfix[1:]) + 1
146  title = f'Run {run}: Scaler S{scaler};Counts;# Events'
147  name = f'counts;Scaler{scaler}'
148  tool.defineHistogram(name, title=title, path='BeamElements', type='TH1F',
149  xbins=20000, xmin=-0.5, xmax=19999.5)
150 
151  scalerCoincedenceGroup = helper.addGroup(tileTBBeamMonAlg, 'Scaler12', 'TestBeam')
152  scalerCoincedenceGroup.defineHistogram('counts12;Scaler12', path='BeamElements', type='TH1F',
153  title=f'Run {run}: Scaler S1 and S2 coincedence;counts;# Events',
154  xbins=20000, xmin=-0.5, xmax=19999.5)
155 
156 
157  beamChambers = ['BC1', 'BC2']
158  beamChambersArray = helper.addArray([beamChambers], tileTBBeamMonAlg, 'BeamChamber', topPath='TestBeam')
159  for postfix, tool in beamChambersArray.Tools.items():
160  beamChamber = postfix[1:]
161 
162  for coordinate in ['X', 'Y']:
163  title = f'Run {run}: {beamChamber}{coordinate} Coordinate;{coordinate}[mm];Counts'
164  name = f'{beamChamber}{coordinate};{beamChamber}{coordinate}hist'
165  tool.defineHistogram(name, title=title, path='BeamElements', type='TH1F',
166  xbins=201, xmin=-100.5, xmax=100.5)
167 
168  tool.defineHistogram(f'{beamChamber}X,{beamChamber}Y;{beamChamber}Profile', path='BeamElements',
169  type='TH2D', title=f'Run {run}: {beamChamber} Beam Profile;X[mm];Y[mm]',
170  xbins=1000, xmin=-100, xmax=100, ybins=1000, ymin=-100, ymax=100)
171 
172  tool.defineHistogram(f'{beamChamber}Xsum,{beamChamber}Ysum;{beamChamber}ProfileSum', path='BeamElements',
173  type='TH2D', title=f'Run {run}: {beamChamber} Beam Profile Sum;X [mm];Y [mm]',
174  xbins=1000, xmin=-300, xmax=0, ybins=1000, ymin=-300, ymax=0)
175 
176 
177  impactProfileGroup = helper.addGroup(tileTBBeamMonAlg, 'ImpactProfile', 'TestBeam')
178  impactProfileGroup.defineHistogram('Ximp,Yimp;ImpactProfile', path='BeamElements', type='TH2F',
179  title=f'Run {run}: {beamChamber} Impact Profile;X [mm];Y [mm]',
180  xbins=200, xmin=-100, xmax=100, ybins=200, ymin=-100, ymax=100)
181 
182 
183  if cellContainer:
184  cherenkovEnergyArray = helper.addArray([nCherenkov], tileTBBeamMonAlg, 'CherenkovVsEnergy', topPath='TestBeam')
185  for postfix, tool in cherenkovEnergyArray.Tools.items():
186  cherenkov = int(postfix[1:]) + 1
187  title = f'Run {run}: Cherenkov {cherenkov} Amplitude vs Total Energy;Energy [pC]; Amplitude [ADC]'
188  name = f'totalEnergy,amplitude;Cher{cherenkov}Energy'
189  tool.defineHistogram(name, title=title, path='BeamElements', type='TH2F',
190  xbins=150, xmin=0, xmax=150, ybins=410, ymin=0, ymax=4096)
191 
192  cellEneXimpGroup = helper.addGroup(tileTBBeamMonAlg, 'CellEnergyImpactX', 'TestBeam')
193  cellEneXimpGroup.defineHistogram('Ximp,cellEnergy;CellEnergyImpactX', path='BeamElements', type='TH2F',
194  title=f'Run {run}: Maximum Cell Energy vs Impact X;X [mm];Cell Energy [pc]',
195  xbins=200, xmin=-100, xmax=100, ybins=150, ymin=0, ymax=150)
196 
197  cellEneYimpGroup = helper.addGroup(tileTBBeamMonAlg, 'CellEnergyImpactY', 'TestBeam')
198  cellEneYimpGroup.defineHistogram('Yimp,cellEnergy;CellEnergyImpactY', path='BeamElements', type='TH2F',
199  title=f'Run {run}: Maximum Cell Energy vs Impact Y;Y [mm];Cell Energy [pc]',
200  xbins=200, xmin=-100, xmax=100, ybins=150, ymin=0, ymax=150)
201 
202  totalEneXimpGroup = helper.addGroup(tileTBBeamMonAlg, 'TotalEnergyImpactX', 'TestBeam')
203  totalEneXimpGroup.defineHistogram('Ximp,totalEnergy;TotalEnergyImpactX', path='BeamElements', type='TH2F',
204  title=f'Run {run}: Total Energy vs Impact X;X [mm];Total Energy [pc]',
205  xbins=200, xmin=-100, xmax=100, ybins=150, ymin=0, ymax=150)
206 
207  totalEneYimpGroup = helper.addGroup(tileTBBeamMonAlg, 'TotalEnergyImpactY', 'TestBeam')
208  totalEneYimpGroup.defineHistogram('Yimp,totalEnergy;TotalEnergyImpactY', path='BeamElements', type='TH2F',
209  title=f'Run {run}: Total Energy vs Impact Y;Y [mm];Total Energy [pc]',
210  xbins=200, xmin=-100, xmax=100, ybins=150, ymin=0, ymax=150)
211 
212 
213  totalEneS1Group = helper.addGroup(tileTBBeamMonAlg, 'ScinCalEnergy', 'TestBeam')
214  totalEneS1Group.defineHistogram('amplitude,totalEnergy;ScinCalEnergy', path='BeamElements', type='TH2F',
215  title=f'Run {run}: Total Energy vs S1 amplitude; Amplitude [ADC];Total Energy [pc]',
216  xbins=410, xmin=0, xmax=4096, ybins=150, ymin=0, ymax=150)
217 
218 
219  accumalator = helper.result()
220  result.merge(accumalator)
221  return result
222 
223 
python.JetAnalysisCommon.ComponentAccumulator
ComponentAccumulator
Definition: JetAnalysisCommon.py:302
CaloCellPos2Ntuple.int
int
Definition: CaloCellPos2Ntuple.py:24
TileTBBeamMonitorAlgorithm.TileTBBeamMonitoringConfig
def TileTBBeamMonitoringConfig(flags, fragIDs=[0x100, 0x101, 0x200, 0x201, 0x402], **kwargs)
Definition: TileTBBeamMonitorAlgorithm.py:12
plotBeamSpotVxVal.range
range
Definition: plotBeamSpotVxVal.py:195
TileTBBeamChambersCalibration.updateBeamChambersCalibrations
def updateBeamChambersCalibrations(alg, run)
Definition: TileTBBeamChambersCalibration.py:92
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 Tue Jul 2 2024 22:28:31 for ATLAS Offline Software by   doxygen 1.8.18