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