ATLAS Offline Software
LArNoisyROMonAlg.py
Go to the documentation of this file.
1 #
2 # Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
3 #
4 
5 def LArNoisyROMonConfig(flags, inKey="",
6  NoisyFEBDefStr="(>30 chan with Q>4000)",
7  NoisyHVlineDefStr="",
8  MNBTightFEBDefStr="",
9  MNBTight_PsVetoFEBDefStr="",
10  MNBLooseFEBDefStr=""):
11 
12  from AthenaMonitoring import AthMonitorCfgHelper
13  helper = AthMonitorCfgHelper(flags,'LArNoisyROMonAlgCfg')
14 
15  from AthenaConfiguration.ComponentFactory import CompFactory
16  NoisyFEBDefStr="(>"+str(flags.LAr.NoisyRO.BadChanPerFEB)+" chan with Q>"+str(flags.LAr.NoisyRO.CellQuality)+")"
17  NoisyHVlineDefStr="(>"+str(flags.LAr.NoisyRO.BadHVCut)+" HVlines with >"+str(flags.LAr.NoisyRO.BadHVlineFrac)+" channels with Q>"+str(flags.LAr.NoisyRO.CellQuality)+")"
18 
19  MNBTightFEBDefStr="(>"+str(flags.LAr.NoisyRO.MNBTightCut)+" chan with Q>"+str(flags.LAr.NoisyRO.CellQuality)+")"
20  MNBTight_PsVetoFEBDefStr="(>"+str(flags.LAr.NoisyRO.MNBTight_PsVetoCut[0])+" chan with Q>"+str(flags.LAr.NoisyRO.CellQuality)+") + PS veto (<"+str(flags.LAr.NoisyRO.MNBTight_PsVetoCut[1])+" channels)"
21  MNBLooseFEBDefStr="(>"+str(flags.LAr.NoisyRO.MNBLooseCut)+" chan with Q>"+str(flags.LAr.NoisyRO.CellQuality)+")"
22 
23  return LArNoisyROMonConfigCore(helper,CompFactory.LArNoisyROMonAlg, flags, inKey, NoisyFEBDefStr, NoisyHVlineDefStr, MNBTightFEBDefStr, MNBTight_PsVetoFEBDefStr, MNBLooseFEBDefStr)
24 
25 
26 def LArNoisyROMonConfigCore(helper,algoinstance,flags,
27  inKey="",
28  NoisyFEBDefStr="(>30 chan with Q>4000)",
29  NoisyHVlineDefStr="",
30  MNBTightFEBDefStr="",
31  MNBTight_PsVetoFEBDefStr="",
32  MNBLooseFEBDefStr=""):
33 
34  # first configure known bad FEBs
35  from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
36  cfg=ComponentAccumulator()
37  from LArBadChannelTool.LArBadFebsConfig import LArKnownBadFebCfg, LArKnownMNBFebCfg
38  cfg.merge(LArKnownBadFebCfg(flags))
39  cfg.merge(LArKnownMNBFebCfg(flags))
40 
41  larNoisyROMonAlg = helper.addAlgorithm(algoinstance,'larNoisyROMonAlg')
42 
43  NoisyROMonGroupName="LArNoisyROMonGroup"
44 
45  from LArMonitoring.GlobalVariables import lArDQGlobals
46 
47  larNoisyROMonAlg.NoisyROGroupName=NoisyROMonGroupName
48  larNoisyROMonAlg.SubDetNames=lArDQGlobals.SubDet[0:2]
49  larNoisyROMonAlg.PartitionNames=lArDQGlobals.Partitions[0:4]
50 
51  larNoisyROMonAlg.storeLooseMNBFEBs=True
52  if inKey != "":
53  larNoisyROMonAlg.inputKey=inKey
54 
55  if not flags.Common.isOnline:
56  larNoisyROMonAlg.HVMapKey="LArHVIdMap"
57 
58  # variable for testing on ESD
59  try:
60  LArNoisyROMonForceTrigger
61  except NameError:
62  LArNoisyROMonForceTrigger = False
63 
64  EFNoiseBurstTriggersList=[
65  "HLT_larnoiseburst_L1XE60",
66  "HLT_larnoiseburst_L1XE70",
67  "HLT_larnoiseburst_L1XE80",
68  "HLT_larnoiseburst_L1J100",
69  "HLT_larnoiseburst_L1J40_XE60",
70  "HLT_larnoiseburst_L1J40_XE50",
71  "HLT_larpsall_L1J12_FIRSTEMPTY",
72  "HLT_larpsall_L1J30_FIRSTEMPTY",
73  "HLT_larpsallem_L1EM7_FIRSTEMPTY"
74  ]
75  L1NoiseBurstTriggersList = [
76  "L1_J75",
77  "L1_J100",
78  "L1_J40_XE50",
79  "L1_J40_XE60",
80  "L1_J10_EMPTY",
81  "L1_J30_EMPTY",
82  "L1_EM7_FIRSTEMPTY",
83  "L1_J12_FIRSTEMPTY",
84  "L1_J30_FIRSTEMPTY",
85  "L1_XE60",
86  "L1_XE70",
87  "L1_XE80"
88  ]
89 
90  if flags.DQ.useTrigger or LArNoisyROMonForceTrigger:
91  larNoisyROMonAlg.doTrigger = True
92  larNoisyROMonAlg.EFNoiseBurstTriggers = EFNoiseBurstTriggersList
93  larNoisyROMonAlg.L1NoiseBurstTriggers = L1NoiseBurstTriggersList
94  else:
95  larNoisyROMonAlg.doTrigger = False
96  larNoisyROMonAlg.EFNoiseBurstTriggers = []
97  larNoisyROMonAlg.L1NoiseBurstTriggers = []
98 
99  noisyROGroup = helper.addGroup(
100  larNoisyROMonAlg,
101  NoisyROMonGroupName,
102  '/LAr/NoisyRO/'
103  )
104 
105 
106  larNoisyRO_hist_path='Summary/'
107 
108  # first trees
109  noisyROGroup.defineTree('time,time_ns,algo;LArNoise',path=larNoisyRO_hist_path,
110  treedef='time/i:time_ns/i:algo/b')
111 
112  noisyROGroup.defineTree('candidate_MNB_time,candidate_MNB_time_ns,n_candidate_MNBTight_FEB,n_candidate_MNBTight_PsVeto_FEB,n_candidate_MNBLoose_FEB,v_candidate_MNBTightFEB,v_candidate_MNBTight_PsVetoFEB,v_candidate_MNBLooseFEB;CandidateMNBFebs',path=larNoisyRO_hist_path,
113  treedef='candidate_MNB_time/i:candidate_MNB_time_ns/i:n_candidate_MNBTight_FEB/i:n_candidate_MNBTight_PsVeto_FEB/i:n_candidate_MNBLoose_FEB/i:v_candidate_MNBTightFEB/vector<int>:v_candidate_MNBTight_PsVetoFEB/vector<int>:v_candidate_MNBLooseFEB/vector<int>')
114 
115  #then global histo
116  noisyROGroup.defineHistogram('n_noisyFEBs;NoisyFEB',
117  title='Number of noisy FEB '+ NoisyFEBDefStr + ' per event;# of noisy FEB',
118  type='TH1I',
119  path=larNoisyRO_hist_path,
120  xbins=lArDQGlobals.noisyFEB_Bins,xmin=lArDQGlobals.noisyFEB_Min,xmax=lArDQGlobals.noisyFEB_Max)
121  noisyROGroup.defineHistogram('n_noisyHVlines;NoisyHVlines',
122  title='Number of noisy HVlines '+ NoisyHVlineDefStr + ' per event;# of noisy HVlines',
123  type='TH1I',
124  path=larNoisyRO_hist_path,
125  xbins=lArDQGlobals.noisyFEB_Bins,xmin=lArDQGlobals.noisyFEB_Min,xmax=lArDQGlobals.noisyFEB_Max)
126  noisyROGroup.defineHistogram('LBN;LBN',type='TH1I',
127  title='Event counter per LB',
128  path=larNoisyRO_hist_path,
129  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
130 
131  for subdet in range(0,2):
132  hist_path='/LAr/NoisyRO/'+lArDQGlobals.SubDet[subdet]+'/'
133  slot_low = lArDQGlobals.FEB_Slot[lArDQGlobals.Partitions[subdet*2]][0] - 0.5
134  slot_up = lArDQGlobals.FEB_Slot[lArDQGlobals.Partitions[subdet*2]][1] + 0.5
135  slot_n = int(slot_up - slot_low)
136  ft_low = lArDQGlobals.FEB_Feedthrough[lArDQGlobals.Partitions[subdet*2]][0] - 0.5
137  ft_up = lArDQGlobals.FEB_Feedthrough[lArDQGlobals.Partitions[subdet*2]][1] + 0.5
138  ft_n = int(ft_up - ft_low)
139 
140  darray = helper.addArray([lArDQGlobals.Partitions[2*subdet:2*subdet+2]],larNoisyROMonAlg,lArDQGlobals.SubDet[subdet],topPath='/')
141  # Known bad FEBS
142  darray.defineHistogram('slotBad,FTBad;KnownBadFEB', title='Known Bad FEBs {0} ; Slot ; FT',
143  type='TH2I', path=hist_path,
144  xbins=slot_n,xmin=slot_low,xmax=slot_up,
145  ybins=ft_n, ymin=ft_low, ymax=ft_up,
146  merge='identical')
147 
148  darray.defineHistogram('slotMNB,FTMNB;MNBKnownFEB', title='Known MNB FEBs {0} ; Slot ; FT',
149  type='TH2I', path=hist_path,
150  xbins=slot_n,xmin=slot_low,xmax=slot_up,
151  ybins=ft_n, ymin=ft_low, ymax=ft_up,
152  merge='identical')
153 
154  # 2D plots of noisy FEBs with various cuts
155  darray.defineHistogram('slotNoisy,FTNoisy;NoisyFEBPerEvt', title='Yield of events with {0} FEB noisy -'+NoisyFEBDefStr+' (only vetoed events) ; Slot ; FT',
156  type='TH2I', path=hist_path,
157  xbins=slot_n,xmin=slot_low,xmax=slot_up,
158  ybins=ft_n, ymin=ft_low, ymax=ft_up)
159 
160  darray.defineHistogram('slotTight,FTTight;MNBTightFEBPerEvt', title='Yield of events with {0} FEB noisy -'+MNBTightFEBDefStr+' (only vetoed events) ; Slot ; FT',
161  type='TH2I', path=hist_path,
162  xbins=slot_n,xmin=slot_low,xmax=slot_up,
163  ybins=ft_n, ymin=ft_low, ymax=ft_up)
164 
165  darray.defineHistogram('slotTightCan,FTTightCan;CandidateMNBTightFEBPerEvt', title='Yield of events with {0} FEB noisy -'+MNBTightFEBDefStr+' (only vetoed events) ; Slot ; FT',
166  type='TH2I', path=hist_path,
167  xbins=slot_n,xmin=slot_low,xmax=slot_up,
168  ybins=ft_n, ymin=ft_low, ymax=ft_up)
169 
170  darray.defineHistogram('slot_PsVetoTight,FT_PsVetoTight;MNBTight_PsVetoFEBPerEvt', title='Yield of events with {0} FEB noisy -'+MNBTight_PsVetoFEBDefStr+' (only vetoed events) ; Slot ; FT',
171  type='TH2I', path=hist_path,
172  xbins=slot_n,xmin=slot_low,xmax=slot_up,
173  ybins=ft_n, ymin=ft_low, ymax=ft_up)
174 
175  darray.defineHistogram('slot_PsVetoTightCan,FT_PsVetoTightCan;CandidateMNBTight_PsVetoFEBPerEvt', title='Yield of events with {0} FEB noisy -'+MNBTight_PsVetoFEBDefStr+' (only vetoed events) ; Slot ; FT',
176  type='TH2I', path=hist_path,
177  xbins=slot_n,xmin=slot_low,xmax=slot_up,
178  ybins=ft_n, ymin=ft_low, ymax=ft_up)
179 
180  darray.defineHistogram('slotLoose,FTLoose;MNBLooseFEBPerEvt', title='Yield of events with {0} FEB noisy -'+MNBLooseFEBDefStr+' (only vetoed events) ; Slot ; FT',
181  type='TH2I', path=hist_path,
182  xbins=slot_n,xmin=slot_low,xmax=slot_up,
183  ybins=ft_n, ymin=ft_low, ymax=ft_up)
184 
185  darray.defineHistogram('slotLooseCan,FTLooseCan;CandidateMNBLooseFEBPerEvt', title='Yield of events with {0} FEB noisy -'+MNBLooseFEBDefStr+' (only vetoed events) ; Slot ; FT',
186  type='TH2I', path=hist_path,
187  xbins=slot_n,xmin=slot_low,xmax=slot_up,
188  ybins=ft_n, ymin=ft_low, ymax=ft_up)
189 
190  # 1D plots of noisy events of various type
191  darray.defineHistogram('LBStd;NoisyEvent',type='TH1I', path=hist_path,
192  title='Yield of events flagged as RNB-Standard - {0} ; Luminosity Block; Number of events',
193  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
194 
195  darray.defineHistogram('LBSat;SaturatedNoisyEvent',type='TH1I', path=hist_path,
196  title='Yield of events flagged as RNB-Saturated - {0} ; Luminosity Block',
197  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
198 
199  darray.defineHistogram('LBMNBTight;MNBTightEvent',type='TH1I', path=hist_path,
200  title='Yield of events flagged as MNB-Tight - {0} ; Luminosity Block',
201  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
202 
203  darray.defineHistogram('LBMNBTight_PsVeto;MNBTight_PsVetoEvent',type='TH1I', path=hist_path,
204  title='Yield of events flagged as MNB-Tight_PsVeto - {0} ; Luminosity Block',
205  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
206 
207  darray.defineHistogram('LBMNBLoose;MNBLooseEvent',type='TH1I', path=hist_path,
208  title='Yield of events flagged as MNB-Loose - {0} ; Luminosity Block',
209  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
210 
211  # 1D plots of events still noisy after EventVeto
212  darray.defineHistogram('LBStd_Veto;NoisyEvent_TimeVeto',type='TH1I', path=hist_path,
213  title='Yield of events flagged as RNB-Standard not vetoed by time window - {0} ; Luminosity Block',
214  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
215 
216  darray.defineHistogram('LBSat_Veto;SaturatedNoisyEvent_TimeVeto',type='TH1I', path=hist_path,
217  title='Yield of events flagged as RNB-Saturated not vetoed by time window - {0} ; Luminosity Block',
218  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
219 
220  darray.defineHistogram('LBMNBTight_Veto;MNBTightEvent_TimeVeto',type='TH1I', path=hist_path,
221  title='Yield of events flagged as MNB-Tight not vetoed by time window - {0} ; Luminosity Block',
222  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
223 
224  darray.defineHistogram('LBMNBTight_PsVeto_Veto;MNBTight_PsVetoEvent_TimeVeto',type='TH1I', path=hist_path,
225  title='Yield of events flagged as MNB-Tight_PsVeto not vetoed by time window - {0} ; Luminosity Block',
226  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
227 
228  darray.defineHistogram('LBMNBLoose_Veto;MNBLooseEvent_TimeVeto',type='TH1I', path=hist_path,
229  title='Yield of events flagged as MNB-Loose not vetoed by time window - {0} ; Luminosity Block',
230  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
231 
232  # Trigger histos
233  if larNoisyROMonAlg.doTrigger:
234  siz=len(larNoisyROMonAlg.EFNoiseBurstTriggers)
235  darray.defineHistogram('Triggers;NoisyEventTrigger',type='TH1I', path=hist_path,
236  title='Trigger fired for RNB flagged events - {0} ; Special trigger fired',
237  xbins=siz+1,xmin=0.5,xmax=siz+1.5,
238  xlabels=larNoisyROMonAlg.EFNoiseBurstTriggers.append("NONE"))
239 
240  l1siz=len(larNoisyROMonAlg.L1NoiseBurstTriggers)
241  darray.defineHistogram('L1Triggers;NoisyEventL1Term',type='TH1I', path=hist_path,
242  title='L1 term fired for RNB flagged events - {0} ; Special trigger fired',
243  xbins=l1siz+1,xmin=0.5,xmax=l1siz+1.5,
244  xlabels=larNoisyROMonAlg.L1NoiseBurstTriggers.append("NONE"))
245  for subdet in range(0,4):
246  darray1 = helper.addArray([["HV"+ sd for sd in lArDQGlobals.Partitions[2*subdet:2*subdet+2]]],larNoisyROMonAlg,lArDQGlobals.SubDet[subdet],topPath='/')
247 
248  darray1.defineHistogram('LBHV;HVlineNoisyEvent',type='TH1I', path=hist_path,
249  title='Yield of events flagged by HVlines flag - {0} ; Luminosity Block; Number of events',
250  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
251 
252  darray1.defineHistogram('LBHV_Veto;HVlineNoisyEvent_TimeVeto',type='TH1I', path=hist_path,
253  title='Yield of events flagged by HVlines flag not vetoed by time window - {0} ; Luminosity Block',
254  xbins=lArDQGlobals.LB_Bins,xmin=lArDQGlobals.LB_Min,xmax=lArDQGlobals.LB_Max)
255 
256 
257  cfg.merge(helper.result())
258  return cfg
259 
260 
261 if __name__=='__main__':
262 
263  # Setup logs
264  from AthenaCommon.Constants import INFO
265 
266  # Set the Athena configuration flags
267  from AthenaConfiguration.AllConfigFlags import initConfigFlags
268  flags = initConfigFlags()
269  from LArMonitoring.LArMonConfigFlags import addLArMonFlags
270  flags.addFlagsCategory("LArMon", addLArMonFlags)
271 
272  nightly = '/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/Tier0ChainTests/q431/21.0/v1/'
273  file = 'myESD.pool.root'
274 
275  flags.Input.Files = [nightly+file]
276  flags.Input.isMC = False
277  flags.DQ.useTrigger = True
278 
279  flags.Output.HISTFileName = 'LArNoisyROMonitoringOutput.root'
280 
281  flags.Exec.OutputLevel=INFO
282  flags.lock()
283 
284  # Initialize configuration object, add accumulator, merge, and run.
285  from AthenaConfiguration.MainServicesConfig import MainServicesCfg
286  from AthenaPoolCnvSvc.PoolReadConfig import PoolReadCfg
287  cfg = MainServicesCfg(flags)
288  cfg.merge(PoolReadCfg(flags))
289 
290  # try NoisyRO algo
291  #first geometry
292  from LArGeoAlgsNV.LArGMConfig import LArGMCfg
293  from TileGeoModel.TileGMConfig import TileGMCfg
294  cfg.merge(LArGMCfg(flags))
295  cfg.merge(TileGMCfg(flags))
296  from LArCabling.LArCablingConfig import LArOnOffIdMappingCfg
297  cfg.merge(LArOnOffIdMappingCfg(flags))
298  from LArCalibUtils.LArHVScaleConfig import LArHVScaleCfg
299  cfg.merge(LArHVScaleCfg(flags))
300  # then NoisyROSummary creator
301  from LArCellRec.LArNoisyROSummaryConfig import LArNoisyROSummaryCfg
302  noisyROSumm = LArNoisyROSummaryCfg(flags)
303  noisyROSumm.getEventAlgo("LArNoisyROAlg").OutputKey="LArNoisyROSummary_recomputed"
304  cfg.merge(noisyROSumm)
305  # then monitoring algo
306  noisemon=LArNoisyROMonConfig(flags, inKey="LArNoisyROSummary_recomputed")
307  cfg.merge(noisemon)
308 
309  flags.dump()
310 
311  cfg.run(20) #use cfg.run() to run on all events
python.JetAnalysisCommon.ComponentAccumulator
ComponentAccumulator
Definition: JetAnalysisCommon.py:302
python.LArBadFebsConfig.LArKnownBadFebCfg
def LArKnownBadFebCfg(configFlags, tag=None, dbname=None)
Definition: LArBadFebsConfig.py:11
LArNoisyROMonAlg.LArNoisyROMonConfig
def LArNoisyROMonConfig(flags, inKey="", NoisyFEBDefStr="(>30 chan with Q>4000)", NoisyHVlineDefStr="", MNBTightFEBDefStr="", MNBTight_PsVetoFEBDefStr="", MNBLooseFEBDefStr="")
Definition: LArNoisyROMonAlg.py:5
LArNoisyROMonAlg.LArNoisyROMonConfigCore
def LArNoisyROMonConfigCore(helper, algoinstance, flags, inKey="", NoisyFEBDefStr="(>30 chan with Q>4000)", NoisyHVlineDefStr="", MNBTightFEBDefStr="", MNBTight_PsVetoFEBDefStr="", MNBLooseFEBDefStr="")
Definition: LArNoisyROMonAlg.py:26
python.MainServicesConfig.MainServicesCfg
def MainServicesCfg(flags, LoopMgr='AthenaEventLoopMgr')
Definition: MainServicesConfig.py:310
plotBeamSpotVxVal.range
range
Definition: plotBeamSpotVxVal.py:194
Constants
some useful constants -------------------------------------------------—
LArCablingConfig.LArOnOffIdMappingCfg
def LArOnOffIdMappingCfg(configFlags)
Definition: LArCablingConfig.py:62
LArGMConfig.LArGMCfg
def LArGMCfg(flags)
Definition: LArGMConfig.py:8
python.LArBadFebsConfig.LArKnownMNBFebCfg
def LArKnownMNBFebCfg(configFlags, tag=None, dbname=None)
Definition: LArBadFebsConfig.py:30
LArNoisyROSummaryConfig.LArNoisyROSummaryCfg
def LArNoisyROSummaryCfg(configFlags, **kwargs)
Definition: LArNoisyROSummaryConfig.py:9
python.CaloAddPedShiftConfig.int
int
Definition: CaloAddPedShiftConfig.py:45
python.AllConfigFlags.initConfigFlags
def initConfigFlags()
Definition: AllConfigFlags.py:19
str
Definition: BTagTrackIpAccessor.cxx:11
python.LArHVScaleConfig.LArHVScaleCfg
def LArHVScaleCfg(configFlags)
Definition: LArHVScaleConfig.py:8
python.PoolReadConfig.PoolReadCfg
def PoolReadCfg(flags)
Definition: PoolReadConfig.py:71
TileGMConfig.TileGMCfg
def TileGMCfg(flags)
Definition: TileGMConfig.py:7
Generated on Sun Oct 26 2025 21:13:25 for ATLAS Offline Software by   doxygen 1.8.18