mergeEnergyRamps.L1CaloGeometryConvertor Class Reference

Collaboration diagram for mergeEnergyRamps.L1CaloGeometryConvertor:

Public Member Functions
	__init__ (self)
	LoadReceiverPPMMap (self)
	getPPMfromReceiver (self, ReceiverId)
	getReceiverfromPPM (self, PPMId, strategy_string=None)
	getCoolEm (self, i_eta, i_phi)
	getCoolHad (self, i_eta, i_phi)
	isCoolEm (self, CoolId)
	isCoolHad (self, CoolId)
	getEtaBin (self, CoolId)
	getPhiBin (self, CoolId)
	getMissingReceiverChannels (self, channel_list)
	getReceiverCMCP (self, ReceiverId)
	isPPMFCAL (self, CoolId)
	isPPMOverlap (self, CoolId)
	getOverlapLayer (self, RecCoolId)
	getFCAL23RecEta (self, RecCoolId)

Public Attributes
	coolIdPath = ROOT.PathResolver.find_calib_file("TrigT1Calo/COOLIdDump_v1.txt")
dict	list_of_channels_em = {}
dict	list_of_channels_had = {}
dict	receiver_to_ppm_map = {}
int	UNIX2COOL = 1000000000

Detailed Description

Definition at line 13 of file mergeEnergyRamps.py.

Constructor & Destructor Documentation

__init__()

mergeEnergyRamps.L1CaloGeometryConvertor.__init__

(

self

)

Definition at line 15 of file mergeEnergyRamps.py.

     def __init__(self):
          self.coolIdPath=ROOT.PathResolver.find_calib_file("TrigT1Calo/COOLIdDump_v1.txt")
          input = open(self.coolIdPath)
          self.list_of_channels_em={}
          self.list_of_channels_had={}
 
          for line in input.readlines():
               parts = line.split(' ')
               emCool = parts[4].rstrip()
               hadCool = parts[5].rstrip()
               self.list_of_channels_em[(parts[0],parts[1])]  = '0x'+emCool
               self.list_of_channels_had[(parts[0],parts[1])] = '0x'+hadCool
     
          input.close()
 

Member Function Documentation

getCoolEm()

mergeEnergyRamps.L1CaloGeometryConvertor.getCoolEm	(		self,
			i_eta,
			i_phi )

Definition at line 111 of file mergeEnergyRamps.py.

     def getCoolEm(self,i_eta,i_phi):
          if (str(i_eta),str(i_phi)) in self.list_of_channels_em:
              cool = self.list_of_channels_em[(str(i_eta),str(i_phi))]
              cool.rstrip()
              cool.lstrip()
              return (cool)
          else:
              return ('')         
     
     

getCoolHad()

mergeEnergyRamps.L1CaloGeometryConvertor.getCoolHad	(		self,
			i_eta,
			i_phi )

Definition at line 121 of file mergeEnergyRamps.py.

     def getCoolHad(self,i_eta,i_phi):
          if (str(i_eta),str(i_phi)) in self.list_of_channels_had:
              cool = self.list_of_channels_had[(str(i_eta),str(i_phi))]
              cool.rstrip()
              cool.lstrip()
              return (cool)
          else:
              return ('')         
     

getEtaBin()

mergeEnergyRamps.L1CaloGeometryConvertor.getEtaBin	(		self,
			CoolId )

Definition at line 136 of file mergeEnergyRamps.py.

     def getEtaBin(self,CoolId):
          if self.isCoolEm(CoolId):
            channel = [item[0] for item in self.list_of_channels_em.items() if item[1]==CoolId]
            return int(channel[0][0])
          elif self.isCoolHad(CoolId):
            channel = [item[0] for item in self.list_of_channels_had.items() if item[1]==CoolId]
            return int(channel[0][0])
          else:
            return -1
 

getFCAL23RecEta()

mergeEnergyRamps.L1CaloGeometryConvertor.getFCAL23RecEta	(		self,
			RecCoolId )

Definition at line 215 of file mergeEnergyRamps.py.

     def getFCAL23RecEta(self,RecCoolId):
 
       ppm_id = self.getPPMfromReceiver(RecCoolId)
 
       if (not self.isPPMFCAL(ppm_id)) or (not self.isCoolHad(ppm_id)):
         return None
       eta_bin = self.getEtaBin(ppm_id)  
 
       RecCoolInt = int(RecCoolId,16)
       if RecCoolInt%2 == 1:
         isRecOdd = True
       else:
         isRecOdd = False
 
       if eta_bin>0:
         if isRecOdd:
           return 'LowEta'
         else:
           return 'HighEta'
       else:
         if isRecOdd:
           return 'HighEta'
         else:
           return 'LowEta' 
         

getMissingReceiverChannels()

mergeEnergyRamps.L1CaloGeometryConvertor.getMissingReceiverChannels	(		self,
			channel_list )

Definition at line 156 of file mergeEnergyRamps.py.

     def getMissingReceiverChannels(self, channel_list):
 
       missing_channels= [channel for channel in self.receiver_to_ppm_map.keys() if channel not in channel_list]
       return missing_channels
 
 

getOverlapLayer()

mergeEnergyRamps.L1CaloGeometryConvertor.getOverlapLayer	(		self,
			RecCoolId )

Definition at line 197 of file mergeEnergyRamps.py.

     def getOverlapLayer(self,RecCoolId):
 
       ppm_id = self.getPPMfromReceiver(RecCoolId)
 
       if not self.isPPMOverlap(ppm_id):
         return None
 
       cabling = self.getReceiverCMCP(RecCoolId)
       if cabling[0] < 2:                 # unconnected channel has barrel crate nr.
         return 'Unconnected'
       elif cabling[2] == 0:
         return 'EMEC'
       elif cabling[2] == 2:
         return 'EMB'
       else:
         print ("Error in GetOverlapLayer, can't determine layer!")
         return None        
 

getPhiBin()

mergeEnergyRamps.L1CaloGeometryConvertor.getPhiBin	(		self,
			CoolId )

Definition at line 146 of file mergeEnergyRamps.py.

     def getPhiBin(self,CoolId):
          if self.isCoolEm(CoolId):
            channel = [item[0] for item in self.list_of_channels_em.items() if item[1]==CoolId]
            return int(channel[0][1])
          elif self.isCoolHad(CoolId):
            channel = [item[0] for item in self.list_of_channels_had.items() if item[1]==CoolId]
            return int(channel[0][1])
          else:
            return -1
 

getPPMfromReceiver()

mergeEnergyRamps.L1CaloGeometryConvertor.getPPMfromReceiver	(		self,
			ReceiverId )

Definition at line 71 of file mergeEnergyRamps.py.

     def getPPMfromReceiver(self,ReceiverId):
     
       if ReceiverId in self.receiver_to_ppm_map:
         return self.receiver_to_ppm_map[ReceiverId]
       else:
         return ''   
 

getReceiverCMCP()

mergeEnergyRamps.L1CaloGeometryConvertor.getReceiverCMCP	(		self,
			ReceiverId )

Definition at line 162 of file mergeEnergyRamps.py.

     def getReceiverCMCP(self,ReceiverId):
     
       recI=int(ReceiverId,16)
 
       crate = recI/1024
       recI = recI - crate*1024
 
       module = recI/64
       recI = recI - module*64
 
       conn = recI/16
       recI = recI - conn*16
 
       pair = recI
 
       return [crate,module,conn,pair]
     

getReceiverfromPPM()

mergeEnergyRamps.L1CaloGeometryConvertor.getReceiverfromPPM	(		self,
			PPMId,
			strategy_string = None )

Definition at line 78 of file mergeEnergyRamps.py.

     def getReceiverfromPPM(self,PPMId,strategy_string=None):
        
       ReceiverChannels = [item[0] for item in self.receiver_to_ppm_map.items() if item[1]==PPMId]       
 
       if strategy_string is None:
         print (" Warning! in getReceiverfromPPM no runtype give, using default!")
         return ReceiverChannels[0]
 
       if self.isPPMFCAL(PPMId) and self.isCoolHad(PPMId):      # pick correct FCAL23 channel
 
         if strategy_string == "GainOneOvEmbFcalHighEta":
           for channel in ReceiverChannels:
             if self.getFCAL23RecEta(channel) == 'HighEta':
                return channel
         if strategy_string == "GainOneOvEmecFcalLowEta":
           for channel in ReceiverChannels:
             if self.getFCAL23RecEta(channel) == 'LowEta':
                return channel
 
       elif self.isPPMOverlap(PPMId):
 
         if strategy_string == "GainOneOvEmbFcalHighEta":
           for channel in ReceiverChannels:
             if self.getOverlapLayer(channel) == 'EMB':
                return channel
         if strategy_string == "GainOneOvEmecFcalLowEta":
           for channel in ReceiverChannels:
             if self.getOverlapLayer(channel) == 'EMEC':
                return channel
 
       else:
         return ReceiverChannels[0]
 

isCoolEm()

mergeEnergyRamps.L1CaloGeometryConvertor.isCoolEm	(		self,
			CoolId )

Definition at line 130 of file mergeEnergyRamps.py.

     def isCoolEm(self,CoolId):
          return (CoolId in self.list_of_channels_em.values())
 

isCoolHad()

mergeEnergyRamps.L1CaloGeometryConvertor.isCoolHad	(		self,
			CoolId )

Definition at line 133 of file mergeEnergyRamps.py.

     def isCoolHad(self,CoolId):
          return (CoolId in self.list_of_channels_had.values())
 

isPPMFCAL()

mergeEnergyRamps.L1CaloGeometryConvertor.isPPMFCAL	(		self,
			CoolId )

Definition at line 179 of file mergeEnergyRamps.py.

     def isPPMFCAL(self,CoolId):
     
       eta_bin = self.getEtaBin(CoolId)  
 
       if eta_bin >= 32 or eta_bin <= -36:
         return True
       else:
         return False
 
 

isPPMOverlap()

mergeEnergyRamps.L1CaloGeometryConvertor.isPPMOverlap	(		self,
			CoolId )

Definition at line 189 of file mergeEnergyRamps.py.

     def isPPMOverlap(self,CoolId):
 
       eta_bin = self.getEtaBin(CoolId)  
       if self.isCoolEm(CoolId) is True and (eta_bin == 14 or eta_bin == -15):
         return True
       else:
         return False
 

LoadReceiverPPMMap()

mergeEnergyRamps.L1CaloGeometryConvertor.LoadReceiverPPMMap

(

self

)

Definition at line 30 of file mergeEnergyRamps.py.

     def LoadReceiverPPMMap(self):
     
         self.receiver_to_ppm_map={}
         self.UNIX2COOL = 1000000000
 
         # get database service and open database
         dbSvc = cool.DatabaseSvcFactory.databaseService()
 
         dbString = 'oracle://ATLAS_COOLPROD;schema=ATLAS_COOLONL_TRIGGER;dbname=CONDBR2'
         try:
           db = dbSvc.openDatabase(dbString, False)        
         except Exception as e:
           print ('Error: Problem opening database', e)
           sys.exit(1)
 
         folder_name = "/TRIGGER/Receivers/RxPpmIdMap"
         folder=db.getFolder(folder_name)
       
         startUtime = int(time.time())
         endUtime = int(time.time())
         startValKey = startUtime * self.UNIX2COOL
         endValKey = endUtime * self.UNIX2COOL
         chsel = cool.ChannelSelection(0,sys.maxsize)
 
         try:
           itr=folder.browseObjects(startValKey, endValKey, chsel)
         except Exception as e:
           print (e)
           sys.exit(1)
 
         for row in itr:
           ReceiverId = hex(int(row.channelId()))
           payload = row.payload()
           PPMId = hex(int(payload['ppmid']))
           self.receiver_to_ppm_map[ReceiverId]= PPMId
  
#         print (self.receiver_to_ppm_map)
         # close database
         db.closeDatabase()
 
 

Member Data Documentation

coolIdPath

mergeEnergyRamps.L1CaloGeometryConvertor.coolIdPath = ROOT.PathResolver.find_calib_file("TrigT1Calo/COOLIdDump_v1.txt")

Definition at line 16 of file mergeEnergyRamps.py.

list_of_channels_em

dict mergeEnergyRamps.L1CaloGeometryConvertor.list_of_channels_em = {}

Definition at line 18 of file mergeEnergyRamps.py.

list_of_channels_had

dict mergeEnergyRamps.L1CaloGeometryConvertor.list_of_channels_had = {}

Definition at line 19 of file mergeEnergyRamps.py.

receiver_to_ppm_map

dict mergeEnergyRamps.L1CaloGeometryConvertor.receiver_to_ppm_map = {}

Definition at line 32 of file mergeEnergyRamps.py.

UNIX2COOL

int mergeEnergyRamps.L1CaloGeometryConvertor.UNIX2COOL = 1000000000

Definition at line 33 of file mergeEnergyRamps.py.

---

The documentation for this class was generated from the following file:

• mergeEnergyRamps.py