python.TileCalibTools.TileASCIIParser2 Class Reference

Inheritance diagram for python.TileCalibTools.TileASCIIParser2:

Collaboration diagram for python.TileCalibTools.TileASCIIParser2:

Public Member Functions
def	__init__ (self, fileName, calibId="", readGain=True)
def	getData (self, ros, drawer, channel, adc, iov=(MAXRUN, MAXLBK))
def	getDict (self)
def	channel2PMT (self, ros, drawer, chan)
def	PMT2channel (self, partition, pmt)

Private Attributes
	__dataDict
	__manyIOVs
	__readGain

Detailed Description

This is a class capable of parsing TileCal conditions data stored in
ASCII files. This version of parser can be used when mutiple IOVs are
given in the file. First column is (run,lumi) pair in this case

Definition at line 1218 of file TileCalibTools.py.

Constructor & Destructor Documentation

__init__()

def python.TileCalibTools.TileASCIIParser2.__init__	(		self,
			fileName,
			calibId = "",
			readGain = True
	)

Input:
- fileName          : input file name
- calibId           : like Ped, Las, ... or (r,l) or (run,lumi) but can be empty string as well
- readGain          : if False, no gain field in input file

Definition at line 1226 of file TileCalibTools.py.

    def __init__(self, fileName, calibId="", readGain=True):
        """
        Input:
        - fileName          : input file name
        - calibId           : like Ped, Las, ... or (r,l) or (run,lumi) but can be empty string as well
        - readGain          : if False, no gain field in input file
        """
 
        TileCalibLogger.__init__(self,"TileASCIIParser2")
        self.__dataDict = {}
        self.__manyIOVs = (calibId=="(run,lumi)" or calibId=="(r,l)" )
        self.__readGain = readGain
        iov=(0,0)
        gain=-1
 
        try:
            lines = open(fileName,"r").readlines()
        except Exception as e:
            self.log().error( "TileCalibASCIIParser2::ERROR: Problem opening input file:" )
            self.log().error( e )
            return
 
        self.log().info("Parsing file %s",fileName)
        if len(calibId)>0:
            self.log().info("Looking for prefix %s",calibId)
 
        for line in lines:
            fields = line.strip().split()
            #=== ignore empty and comment lines
            if not len(fields)          :
                continue
            if fields[0].startswith("#"):
                continue
 
            #=== read in fields
            if len(calibId)>0:
                pref = fields[0]
                frag = fields[1]
                chan = fields[2]
                if str(chan)[0:2].lower() == "pm":
                    chan = self.PMT2channel(frag,fields.pop(2))
                if readGain:
                    gain = fields[3]
                    data = fields[4:]
                else:
                    data = fields[3:]
 
                #=== check for correct calibId
                if self.__manyIOVs:
                    iov=tuple(int(i) for i in pref[1:-1].split(","))
                    if len(iov)!=2 or pref[0]!="(" or pref[-1]!=")":
                        raise Exception("%s is not %s IOV" % (pref,calibId))
                elif pref!=calibId:
                    raise Exception("%s is not calibId=%s" % (pref,calibId))
            else:
                frag = fields[0]
                chan = fields[1]
                if str(chan)[0:2].lower() == "pm":
                    chan = self.PMT2channel(frag,fields.pop(2))
                if readGain:
                    gain  = fields[2]
                    data = fields[3:]
                else:
                    data = fields[2:]
 
            #=== decode fragment
            if frag.startswith('0x') or frag.startswith('-0x'):
                frg = int(frag,16)
                ros = frg>>8
                if frg<0:
                    mod = (-frg)&255
                else:
                    mod = frg&255
            elif (frag.startswith("AUX") or
                  frag.startswith("LBA") or
                  frag.startswith("LBC") or
                  frag.startswith("EBA") or
                  frag.startswith("EBC") or
                  frag.startswith("ALL") or
                  frag.startswith("XXX") ):
                part_dict = {'AUX':0,'LBA':1,'LBC':2,'EBA':3,'EBC':4,'ALL':5,'XXX':-1}
                partname = str(frag[0:3])
                ros=part_dict[partname]
                mod = int(frag[3:])-1
            else:
                raise Exception("Unknown fragment %s" % frag)
 
            chn = int(chan)
            adc = int(gain)
 
            #=== fill dictionary
            if ros<0:
                rosmin=0
                rosmax=5
            elif ros>=5:
                rosmin=1
                rosmax=5
            else:
                rosmin=ros
                rosmax=ros+1
 
            if mod<0 or mod>=64:
                modmin=0
                modmax=64
            else:
                modmin=mod
                modmax=mod+1
 
            allchannels=True
            if chn<-2:
                chnmin=0
                chnmax=-chn
            elif chn<0:
                chnmin=0
                chnmax=48
                allchannels=(chn==-1) # if chn=-2 only connected channels will be updated
            else:
                chnmin=chn
                chnmax=chn+1
 
            if adc<-1:
                adcmin=0
                adcmax=-adc
            elif adc<0:
                adcmin=0
                adcmax=2
            else:
                adcmin=adc
                adcmax=adc+1
 
            for ros in range(rosmin,rosmax):
                for mod in range(modmin,modmax):
                    for chn in range(chnmin,chnmax):
                        if allchannels or self.channel2PMT(ros,mod,chn)>0:
                            for adc in range (adcmin,adcmax):
                                dictKey = (ros,mod,chn,adc)
                                if self.__manyIOVs:
                                    if dictKey in self.__dataDict:
                                        self.__dataDict[dictKey] += [(iov,data)]
                                    else:
                                        self.__dataDict[dictKey] = [(iov,data)]
                                else:
                                    self.__dataDict[dictKey] = data
 

Member Function Documentation

channel2PMT()

def python.TileCalibTools.TileASCIIParser2.channel2PMT	(		self,
			ros,
			drawer,
			chan
	)

Definition at line 1388 of file TileCalibTools.py.

    def channel2PMT(self,ros,drawer,chan):
        "Convert channel numbet to PMT number, negative for disconnected channels"
        "This takes ros [1-4], drawer [0-63], chan [0-47]"
 
        chan2PMT_LB=[  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12,
                      13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24,
                      27, 26, 25, 30, 29, 28,-33,-32, 31, 36, 35, 34,
                      39, 38, 37, 42, 41, 40, 45,-44, 43, 48, 47, 46 ]
 
        chan2PMT_EB=[  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12,
                      13, 14, 15, 16, 17, 18,-19,-20, 21, 22, 23, 24,
                     -27,-26,-25,-31,-32,-28, 33, 29, 30,-36,-35, 34,
                      44, 38, 37, 43, 42, 41,-45,-39,-40,-48,-47,-46 ]
 
        chan2PMT_Sp=[ -1, -2, -3, -4,  5,  6,  7,  8,  9, 10, 11, 12, # first 4 do not exist
                      13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, # PMT 19 and 20 exist
                     -27,-26,-25,-31,-32,-28, 33, 29, 30,-36,-35, 34,
                      44, 38, 37, 43, 42, 41,-45,-39,-40,-48,-47,-46 ]
 
        if ros <= 2:
            pmt = chan2PMT_LB[chan]
        elif (ros == 3 and drawer == 14) or (ros == 4 and drawer == 17):
            pmt = chan2PMT_Sp[chan]
        else:
            pmt = chan2PMT_EB[chan]
 
        return pmt
 

getData()

def python.TileCalibTools.TileASCIIParser2.getData	(		self,
			ros,
			drawer,
			channel,
			adc,
			iov = (MAXRUN,MAXLBK)
	)

Definition at line 1371 of file TileCalibTools.py.

    def getData(self, ros, drawer, channel, adc, iov=(MAXRUN,MAXLBK)):
        dictKey = (int(ros), int(drawer), int(channel), int(adc))
        data = self.__dataDict.get(dictKey,[])
        if self.__manyIOVs and len(data)>0:
            before= [i for i in sorted(data) if i[0] <= iov ]
            if len(before)>0:
                data = before[-1][1]
            else:
                data = []
        return data
 

getDict()

def python.TileCalibTools.TileASCIIParser2.getDict

(

self

)

Definition at line 1383 of file TileCalibTools.py.

    def getDict(self):
        import copy
        return copy.deepcopy(self.__dataDict)
 

PMT2channel()

def python.TileCalibTools.TileASCIIParser2.PMT2channel	(		self,
			partition,
			pmt
	)

Definition at line 1417 of file TileCalibTools.py.

    def PMT2channel(self,partition,pmt):
        "Convert PMT number to channel numbet"
        "This takes partition (LBA,LBC,EBA,EBC) and pmt [1-48]"
 
        chan2PMT_LB=[  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12,
                      13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24,
                      27, 26, 25, 30, 29, 28, 33, 32, 31, 36, 35, 34,
                      39, 38, 37, 42, 41, 40, 45, 44, 43, 48, 47, 46 ]
 
        chan2PMT_EB=[  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12,
                      13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24,
                      27, 26, 25, 31, 32, 28, 33, 29, 30, 36, 35, 34,
                      44, 38, 37, 43, 42, 41, 45, 39, 40, 48, 47, 46 ]
 
        chan = -1
        pm=abs(int(pmt))
 
        if pm>0 and pm<=48:
            if str(partition)[0].upper() == "E":
                chan = chan2PMT_EB.index(pm)
            else:
                chan = chan2PMT_LB.index(pm)
 
        return chan
 
#======================================================================
#===
#=== TileASCIIParser3
#===
#======================================================================
 
#______________________________________________________________________

Member Data Documentation

__dataDict

python.TileCalibTools.TileASCIIParser2.__dataDict

private

Definition at line 1235 of file TileCalibTools.py.

__manyIOVs

python.TileCalibTools.TileASCIIParser2.__manyIOVs

private

Definition at line 1236 of file TileCalibTools.py.

__readGain

python.TileCalibTools.TileASCIIParser2.__readGain

private

Definition at line 1237 of file TileCalibTools.py.

---

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

• TileCalibTools.py