mergeEnergyRamps Namespace Reference

Classes
class	GainsFromOracle
class	GainsFromSqlite
class	L1CaloGeometryConvertor

Functions
def	WriteSqlite (name, input_dict)
def	merge_gains (gains1, gains2, gains3, reference_gains, forced_list, geometry_map, writeAllChannels)

Variables
	parser
	action
	type
	dest
	help
	options
	args
	geometry_convertor
	gains_1
	gains_2
	gains_3
	gains_reference
	gains_to_load
	output_file_name

Function Documentation

merge_gains()

def mergeEnergyRamps.merge_gains	(		gains1,
			gains2,
			gains3,
			reference_gains,
			forced_list,
			geometry_map,
			writeAllChannels
	)

Definition at line 421 of file mergeEnergyRamps.py.

def merge_gains(gains1,gains2,gains3,reference_gains,forced_list,geometry_map,writeAllChannels):
 
     output_gains={}
 
     n_files_Tile        = 0
     n_files_LowEta_EMEC = 0
     n_files_HighEta_EMB = 0
 
#loop over gains3, fill in
 
     if gains3 is not None:
       good_gains=gains3.getGoodGains()
       print (" Using run ", gains3.run_nr, " run strategy= ", gains3.strategy       )
 
       if gains3.strategy == "GainOne":
         n_files_Tile = n_files_Tile + 1
       if gains3.strategy == "GainOneOvEmecFcalLowEta":
         n_files_LowEta_EMEC =  n_files_LowEta_EMEC + 1
       if gains3.strategy == "GainOneOvEmbFcalHighEta":
         n_files_HighEta_EMB =  n_files_HighEta_EMB + 1
 
 
       for ppm_channel in good_gains.keys():
         gain       = good_gains[ppm_channel][0]
         error_code = good_gains[ppm_channel][1]
         rec_chan   = geometry_map.getReceiverfromPPM(ppm_channel,gains3.strategy)
         output_gains[rec_chan]=[gain,error_code]
     else:
       print ("Ignoring File 3, probably not specified"      )
       
#loop over gains2, fill in
 
     if gains2 is not None:
       good_gains=gains2.getGoodGains()
       print (" Using run ", gains2.run_nr, " run strategy= ", gains2.strategy       )
 
       if gains2.strategy == "GainOne":
         n_files_Tile = n_files_Tile + 1
       if gains2.strategy == "GainOneOvEmecFcalLowEta":
         n_files_LowEta_EMEC =  n_files_LowEta_EMEC + 1
       if gains2.strategy == "GainOneOvEmbFcalHighEta":
         n_files_HighEta_EMB =  n_files_HighEta_EMB + 1
 
       for ppm_channel in good_gains.keys():
         gain       = good_gains[ppm_channel][0]
         error_code = good_gains[ppm_channel][1]
         rec_chan   = geometry_map.getReceiverfromPPM(ppm_channel,gains2.strategy)
         output_gains[rec_chan]=[gain,error_code]
     else:
       print ("Ignoring File 2, probably not specified"      )
     
       
#loop over gains 1, fill in
 
     if gains1 is not None:
       good_gains=gains1.getGoodGains()
       print (" Using run ", gains1.run_nr, " run strategy= ", gains1.strategy       )
 
       if gains1.strategy == "GainOne":
         n_files_Tile = n_files_Tile + 1
       if gains1.strategy == "GainOneOvEmecFcalLowEta":
         n_files_LowEta_EMEC =  n_files_LowEta_EMEC + 1
       if gains1.strategy == "GainOneOvEmbFcalHighEta":
         n_files_HighEta_EMB =  n_files_HighEta_EMB + 1
 
       for ppm_channel in good_gains.keys():
         gain       = good_gains[ppm_channel][0]
         error_code = good_gains[ppm_channel][1]
         rec_chan   = geometry_map.getReceiverfromPPM(ppm_channel,gains1.strategy)
         output_gains[rec_chan]=[gain,error_code]
     else:
       print ("Ignoring File 1, probably not specified"      )
     
 
#read in forced list, overwrite
 
     if forced_list is not None:
       myfile = open(forced_list,'r')
       for line in myfile.readlines():
         line.rstrip()
         line.lstrip()
         line_cont  = line.split(' ')
         line_cont  = [iii for iii in line_cont if not iii == ''] 
         rec_chan   = line_cont[0]
         gain       = float(line_cont[1])
         error_code = 10                    # code for forced channels
         output_gains[rec_chan]=[gain,error_code]
         print ("forcing channel ", rec_chan, " to value ", gain )
     else:
       print ("Ignoring forced channel list")
 
#check channels that haven't been found yet, take them from reference
 
     if writeAllChannels:
 
       print ("Adding gains for missing channels from Oracle"     )
       missing_channels = geometry_convertor.getMissingReceiverChannels(output_gains.keys())    
       default_gains = reference_gains.getGoodGains()
       for channel in missing_channels:
         gain = default_gains[channel]
         error_code = 100
         output_gains[channel]=[gain,error_code]
 
#       print ("missing receiver=",channel, "PPM=", geometry_map.getPPMfromReceiver(channel) )
 
     if (not n_files_Tile == 1) or (not  n_files_LowEta_EMEC == 1) or (not n_files_HighEta_EMB == 1):
       print ("WARNING! input files do not allow to calibrate all partitions, using defaults where needed" )
 
     return output_gains
 
 
 

WriteSqlite()

def mergeEnergyRamps.WriteSqlite	(		name,
			input_dict
	)

Definition at line 241 of file mergeEnergyRamps.py.

def WriteSqlite(name,input_dict):
  
  UNIX2COOL = 1000000000
  
  dbSvc = cool.DatabaseSvcFactory.databaseService()
  connectString = 'sqlite://;schema='+name+';dbname=L1CALO'
# folder_name="/TRIGGER/Receivers/Factors/CalibGains"
 
 
  print ('\nrecreating database file:',name)
  dbSvc.dropDatabase( connectString )
  db = dbSvc.createDatabase( connectString )
 
  spec = cool.RecordSpecification()
  spec.extend("factor", cool.StorageType.Float)
  spec.extend("status", cool.StorageType.UInt32 )
  folderSpec = cool.FolderSpecification(cool.FolderVersioning.SINGLE_VERSION, spec)
 
  now = int(time.time())
 
  since = now*UNIX2COOL
  until = cool.ValidityKeyMax
  db.createFolderSet('/TRIGGER')
  db.createFolderSet('/TRIGGER/Receivers')
  db.createFolderSet('/TRIGGER/Receivers/Factors')
 
  folder_description = '<timeStamp>time</timeStamp><addrHeader><address_header service_type="71" clid="1238547719"/></addrHeader><typeName>CondAttrListCollection</typeName>'
  f = db.createFolder( "/TRIGGER/Receivers/Factors/CalibGains", folderSpec, folder_description)
 
  print (" Now creating sqlite file for ", len(input_dict.keys()), " channels")
  for i in input_dict.keys():
    data = cool.Record( spec )
    data['factor'] = input_dict[i][0]
    data['status'] = input_dict[i][1]
    f.storeObject(since,until, data, int(i,16) )
 
  db.closeDatabase()
 

Variable Documentation

action

mergeEnergyRamps.action

Definition at line 539 of file mergeEnergyRamps.py.

args

mergeEnergyRamps.args

Definition at line 546 of file mergeEnergyRamps.py.

dest

mergeEnergyRamps.dest

Definition at line 539 of file mergeEnergyRamps.py.

gains_1

mergeEnergyRamps.gains_1

Definition at line 554 of file mergeEnergyRamps.py.

gains_2

mergeEnergyRamps.gains_2

Definition at line 559 of file mergeEnergyRamps.py.

gains_3

mergeEnergyRamps.gains_3

Definition at line 564 of file mergeEnergyRamps.py.

gains_reference

mergeEnergyRamps.gains_reference

Definition at line 569 of file mergeEnergyRamps.py.

gains_to_load

mergeEnergyRamps.gains_to_load

Definition at line 577 of file mergeEnergyRamps.py.

geometry_convertor

mergeEnergyRamps.geometry_convertor

Definition at line 548 of file mergeEnergyRamps.py.

help

mergeEnergyRamps.help

Definition at line 539 of file mergeEnergyRamps.py.

options

mergeEnergyRamps.options

Definition at line 546 of file mergeEnergyRamps.py.

output_file_name

mergeEnergyRamps.output_file_name

Definition at line 580 of file mergeEnergyRamps.py.

parser

mergeEnergyRamps.parser

Definition at line 537 of file mergeEnergyRamps.py.

type

mergeEnergyRamps.type

Definition at line 539 of file mergeEnergyRamps.py.