ATLAS Offline Software
Public Member Functions | Public Attributes | List of all members
UpdateHVReference.HVRefUpdator Class Reference
Collaboration diagram for UpdateHVReference.HVRefUpdator:

Public Member Functions

def __init__ (self)
 
def read_forced_list (self, forced_file_name)
 
def read_HV_reference (self, ref_file_name)
 
def read_HV_geometry (self, ref_file_name)
 
def update_HV (self)
 
def write_geometry (self, new_file_name)
 
def write_new_file (self, new_file_name)
 

Public Attributes

 reference_HV
 
 reference_AffectedCells
 
 forced_HV
 
 NLayers
 
 LayerName
 
 NCells
 
 UNIX2COOL
 

Detailed Description

Definition at line 12 of file UpdateHVReference.py.

Constructor & Destructor Documentation

◆ __init__()

def UpdateHVReference.HVRefUpdator.__init__ (   self)

Definition at line 14 of file UpdateHVReference.py.

14  def __init__(self):
15 
16  self.reference_HV = {}
17  self.reference_AffectedCells = {}
18  self.forced_HV = {}
19  self.NLayers = {}
20  self.LayerName = {}
21  self.NCells = {}
22  self.UNIX2COOL = 1000000000
23 

Member Function Documentation

◆ read_forced_list()

def UpdateHVReference.HVRefUpdator.read_forced_list (   self,
  forced_file_name 
)

Definition at line 24 of file UpdateHVReference.py.

24  def read_forced_list(self,forced_file_name):
25  print (" Loading forced file", forced_file_name)
26 
27  new_hv = []
28  myfile = open(forced_file_name,'r')
29  for line in myfile.readlines():
30  line.rstrip()
31  line.lstrip()
32  line_cont = line.split(' ')
33  line_cont = [iii for iii in line_cont if not iii == '']
34  rec_chan = line_cont[0]
35 
36  new_hv = [float(line_cont[iii]) for iii in range(2,6)]
37 
38  self.forced_HV[rec_chan]=new_hv
39  print ("updating channel ", rec_chan, " to HV ", new_hv )
40 
41  return
42 

◆ read_HV_geometry()

def UpdateHVReference.HVRefUpdator.read_HV_geometry (   self,
  ref_file_name 
)

Definition at line 102 of file UpdateHVReference.py.

102  def read_HV_geometry(self,ref_file_name):
103  print (" Reading HV geometry from file", ref_file_name)
104  # get database service and open database
105  dbSvc = cool.DatabaseSvcFactory.databaseService()
106 
107  dbString='sqlite://;schema='+ref_file_name+';dbname=L1CALO'
108  try:
109  db = dbSvc.openDatabase(dbString, False)
110  except Exception as e:
111  print ('Error: Problem opening database', e)
112  sys.exit(1)
113 
114  folder_name = '/TRIGGER/L1Calo/V1/Results/RxLayers'
115  folder=db.getFolder(folder_name)
116 
117  startUtime = int(time.time())
118  endUtime = int(time.time())
119  startValKey = startUtime * self.UNIX2COOL
120  endValKey = endUtime * self.UNIX2COOL
121  chsel = cool.ChannelSelection(0,sys.maxint)
122 
123  try:
124  itr=folder.browseObjects(startValKey, endValKey, chsel)
125  except Exception as e:
126  print (e)
127  sys.exit(1)
128 
129  for row in itr:
130  HV_layer_names=4*[0]
131  HV_Ncells = 4*[0]
132 
133  CoolId = hex(int(row.channelId()))
134  payload = row.payload()
135 
136  self.NLayers[CoolId] = struct.unpack('B',payload['NLayers'])[0]
137 
138  HV_layer_names[0] = struct.unpack('B',payload['Name1'])[0]
139  HV_layer_names[1] = struct.unpack('B',payload['Name2'])[0]
140  HV_layer_names[2] = struct.unpack('B',payload['Name3'])[0]
141  HV_layer_names[3] = struct.unpack('B',payload['Name4'])[0]
142 
143  HV_Ncells[0] = struct.unpack('B',payload['NCells1'])[0]
144  HV_Ncells[1] = struct.unpack('B',payload['NCells2'])[0]
145  HV_Ncells[2] = struct.unpack('B',payload['NCells3'])[0]
146  HV_Ncells[3] = struct.unpack('B',payload['NCells4'])[0]
147 
148  self.LayerName[CoolId] = HV_layer_names
149  self.NCells[CoolId] = HV_Ncells
150 
151 
152 # for iii in self.NLayers.keys():
153 # print (iii," NLayers: ",self.NLayers[iii], " Names: ",self.LayerName[iii], " NCells: ",self.NCells[iii])
154 
155 
156  # close database
157  db.closeDatabase()
158 
159  return 0
160 
161 
162 

◆ read_HV_reference()

def UpdateHVReference.HVRefUpdator.read_HV_reference (   self,
  ref_file_name 
)

Definition at line 43 of file UpdateHVReference.py.

43  def read_HV_reference(self,ref_file_name):
44  print (" Loading reference HV file", ref_file_name)
45  # get database service and open database
46  dbSvc = cool.DatabaseSvcFactory.databaseService()
47 
48  dbString='sqlite://;schema='+ref_file_name+';dbname=L1CALO'
49  try:
50  db = dbSvc.openDatabase(dbString, False)
51  except Exception as e:
52  print ('Error: Problem opening database', e)
53  sys.exit(1)
54 
55  folder_name = '/TRIGGER/L1Calo/V1/Results/HVCorrections'
56  folder=db.getFolder(folder_name)
57 
58  startUtime = int(time.time())
59  endUtime = int(time.time())
60  startValKey = startUtime * self.UNIX2COOL
61  endValKey = endUtime * self.UNIX2COOL
62  chsel = cool.ChannelSelection(0,sys.maxint)
63 
64  try:
65  itr=folder.browseObjects(startValKey, endValKey, chsel)
66  except Exception as e:
67  print (e)
68  sys.exit(1)
69 
70  for row in itr:
71  HV_corrections=4*[0]
72  Affected_cells = 4*[0]
73 
74  CoolId = hex(int(row.channelId()))
75  payload = row.payload()
76 
77  HV_corrections[0] = payload['LayerMean1']
78  HV_corrections[1] = payload['LayerMean2']
79  HV_corrections[2] = payload['LayerMean3']
80  HV_corrections[3] = payload['LayerMean4']
81 
82  Affected_cells[0] = struct.unpack('B',payload['AffectedCells1'])[0]
83  Affected_cells[1] = struct.unpack('B',payload['AffectedCells2'])[0]
84  Affected_cells[2] = struct.unpack('B',payload['AffectedCells3'])[0]
85  Affected_cells[3] = struct.unpack('B',payload['AffectedCells4'])[0]
86 
87  self.reference_AffectedCells[CoolId] = Affected_cells
88  self.reference_HV[CoolId] = HV_corrections
89 
90 # for iii in self.reference_AffectedCells.keys():
91 # print (iii," ",self.reference_AffectedCells[iii] )
92 
93 # for iii in self.reference_HV.keys():
94 # print (iii," ",self.reference_HV[iii] )
95 
96  # close database
97  db.closeDatabase()
98 
99  return 0
100 
101 

◆ update_HV()

def UpdateHVReference.HVRefUpdator.update_HV (   self)

Definition at line 163 of file UpdateHVReference.py.

163  def update_HV(self):
164 
165  print (" Updating internal HV list")
166 
167  for iii in self.forced_HV.keys():
168 
169  isNominal = False
170  if (self.forced_HV[iii][0] == 0 or self.forced_HV[iii][0] == 1.) and \
171  (self.forced_HV[iii][1] == 0 or self.forced_HV[iii][1] == 1.) and \
172  (self.forced_HV[iii][2] == 0 or self.forced_HV[iii][2] == 1.) and \
173  (self.forced_HV[iii][3] == 0 or self.forced_HV[iii][3] == 1.) :
174 
175  isNominal = True
176 
177  if isNominal: # if the channel is on nominal HV, doesn't go to the file
178  del self.reference_HV[iii]
179  del self.reference_AffectedCells[iii]
180  else: #
181  self.reference_HV[iii] = self.forced_HV[iii]
182  self.reference_AffectedCells[iii] = [99,99,99,99] # don't know what to put there, hope this is OK
183 
184 
185 # print ("***************************************************************")
186 # print (" Now changed HV dictionary")
187 # print ("***************************************************************")
188 # for iii in self.reference_HV.keys():
189 # print (iii," ",self.reference_HV[iii] )
190 
191  return
192 

◆ write_geometry()

def UpdateHVReference.HVRefUpdator.write_geometry (   self,
  new_file_name 
)

Definition at line 193 of file UpdateHVReference.py.

193  def write_geometry(self,new_file_name):
194  print (" Writing geometry to file " , new_file_name)
195 
196  dbSvc = cool.DatabaseSvcFactory.databaseService()
197  connectString = 'sqlite://;schema='+new_file_name+';dbname=L1CALO'
198 
199 
200  print ('Writing into database file',new_file_name)
201  db = dbSvc.openDatabase(connectString, False)
202 
203  spec = cool.RecordSpecification()
204  spec.extend('NLayers',cool.StorageType.UChar)
205 
206  spec.extend( 'Name1', cool.StorageType.UChar)
207  spec.extend( 'Name2', cool.StorageType.UChar)
208  spec.extend( 'Name3', cool.StorageType.UChar)
209  spec.extend( 'Name4', cool.StorageType.UChar)
210 
211  spec.extend('NCells1',cool.StorageType.UChar)
212  spec.extend('NCells2',cool.StorageType.UChar)
213  spec.extend('NCells3',cool.StorageType.UChar)
214  spec.extend('NCells4',cool.StorageType.UChar)
215 
216  now = int(time.time())
217 
218  since = now*self.UNIX2COOL
219  until = cool.ValidityKeyMax
220 
221  folder_description = "<timeStamp>time</timeStamp><addrHeader><address_header service_type=\"71\" clid=\"1238547719\"/></addrHeader><typeName>CondAttrListCollection</typeName>"
222  f = db.createFolder( "/TRIGGER/L1Calo/V1/Results/RxLayers" , spec,folder_description)
223 
224  for i in self.NLayers.keys():
225 
226  data = cool.Record( spec )
227 
228  data['NLayers'] = struct.pack('B',self.NLayers[i])
229  data['Name1'] = struct.pack('B',self.LayerName[i][0])
230  data['Name2'] = struct.pack('B',self.LayerName[i][1])
231  data['Name3'] = struct.pack('B',self.LayerName[i][2])
232  data['Name4'] = struct.pack('B',self.LayerName[i][3])
233 
234  data['NCells1'] = struct.pack('B',self.NCells[i][0])
235  data['NCells2'] = struct.pack('B',self.NCells[i][1])
236  data['NCells3'] = struct.pack('B',self.NCells[i][2])
237  data['NCells4'] = struct.pack('B',self.NCells[i][3])
238 
239  f.storeObject(since,until, data, int(i,16) )
240 
241  db.closeDatabase()
242 
243  return
244 
245 

◆ write_new_file()

def UpdateHVReference.HVRefUpdator.write_new_file (   self,
  new_file_name 
)

Definition at line 246 of file UpdateHVReference.py.

246  def write_new_file(self,new_file_name):
247 
248  print (" Writing output to file " , new_file_name)
249 
250  dbSvc = cool.DatabaseSvcFactory.databaseService()
251  connectString = 'sqlite://;schema='+new_file_name+';dbname=L1CALO'
252 
253  print ('recreating database file',new_file_name)
254  dbSvc.dropDatabase( connectString )
255  db = dbSvc.createDatabase( connectString )
256 
257  spec = cool.RecordSpecification()
258  spec.extend('RxMean',cool.StorageType.Float)
259 
260  spec.extend( 'AffectedCells1', cool.StorageType.UChar)
261  spec.extend( 'AffectedCells2', cool.StorageType.UChar)
262  spec.extend( 'AffectedCells3', cool.StorageType.UChar)
263  spec.extend( 'AffectedCells4', cool.StorageType.UChar)
264 
265  spec.extend('LayerMean1',cool.StorageType.Float)
266  spec.extend('LayerMean2',cool.StorageType.Float)
267  spec.extend('LayerMean3',cool.StorageType.Float)
268  spec.extend('LayerMean4',cool.StorageType.Float)
269 
270  now = int(time.time())
271 
272  since = now*self.UNIX2COOL
273  until = cool.ValidityKeyMax
274  db.createFolderSet('/TRIGGER')
275  db.createFolderSet('/TRIGGER/L1Calo')
276  db.createFolderSet('/TRIGGER/L1Calo/V1')
277  db.createFolderSet('/TRIGGER/L1Calo/V1/Results')
278 # db.createFolderSet('/TRIGGER/L1Calo/V1/Results/HVCorrections')
279 
280  folder_description = "<timeStamp>time</timeStamp><addrHeader><address_header service_type=\"71\" clid=\"1238547719\"/></addrHeader><typeName>CondAttrListCollection</typeName>"
281  f = db.createFolder( "/TRIGGER/L1Calo/V1/Results/HVCorrections" , spec,folder_description)
282 # f = db.createFolder( folder_name , spec, True )
283 
284  for i in self.reference_HV.keys():
285 
286  data = cool.Record( spec )
287 
288  MeanHV=0.
289  TotalNrCells=0
290 
291  for ilayer in range(4):
292  if self.NCells[i][ilayer]>0:
293  MeanHV += self.reference_HV[i][ilayer]*self.NCells[i][ilayer]
294  TotalNrCells += self.NCells[i][ilayer]
295 
296  MeanHV /= TotalNrCells
297 
298 # data['RxMean'] =(self.reference_HV[i][0]+self.reference_HV[i][1]+self.reference_HV[i][2]+self.reference_HV[i][3])/4.
299  data['RxMean'] = MeanHV
300 
301  data['LayerMean1'] = self.reference_HV[i][0]
302  data['LayerMean2'] = self.reference_HV[i][1]
303  data['LayerMean3'] = self.reference_HV[i][2]
304  data['LayerMean4'] = self.reference_HV[i][3]
305 
306  data['AffectedCells1'] = struct.pack('B',self.reference_AffectedCells[i][0])
307  data['AffectedCells2'] = struct.pack('B',self.reference_AffectedCells[i][1])
308  data['AffectedCells3'] = struct.pack('B',self.reference_AffectedCells[i][2])
309  data['AffectedCells4'] = struct.pack('B',self.reference_AffectedCells[i][3])
310 
311  f.storeObject(since,until, data, int(i,16) )
312 
313  db.closeDatabase()
314 
315  return
316 

Member Data Documentation

◆ forced_HV

UpdateHVReference.HVRefUpdator.forced_HV

Definition at line 18 of file UpdateHVReference.py.

◆ LayerName

UpdateHVReference.HVRefUpdator.LayerName

Definition at line 20 of file UpdateHVReference.py.

◆ NCells

UpdateHVReference.HVRefUpdator.NCells

Definition at line 21 of file UpdateHVReference.py.

◆ NLayers

UpdateHVReference.HVRefUpdator.NLayers

Definition at line 19 of file UpdateHVReference.py.

◆ reference_AffectedCells

UpdateHVReference.HVRefUpdator.reference_AffectedCells

Definition at line 17 of file UpdateHVReference.py.

◆ reference_HV

UpdateHVReference.HVRefUpdator.reference_HV

Definition at line 16 of file UpdateHVReference.py.

◆ UNIX2COOL

UpdateHVReference.HVRefUpdator.UNIX2COOL

Definition at line 22 of file UpdateHVReference.py.

The documentation for this class was generated from the following file:
CaloCellPos2Ntuple.int
int
Definition: CaloCellPos2Ntuple.py:24
plotBeamSpotVxVal.range
range
Definition: plotBeamSpotVxVal.py:195
python.processes.powheg.ZZ.ZZ.__init__
def __init__(self, base_directory, **kwargs)
Constructor: all process options are set here.
Definition: ZZ.py:18
Trk::open
@ open
Definition: BinningType.h:40
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:790
readCCLHist.float
float
Definition: readCCLHist.py:83
Generated on Thu Jul 4 2024 22:31:40 for ATLAS Offline Software by   doxygen 1.8.18