TileDCSDataPlotter.TileDCSDataPlotter Class Reference

Inheritance diagram for TileDCSDataPlotter.TileDCSDataPlotter:

Collaboration diagram for TileDCSDataPlotter.TileDCSDataPlotter:

Public Member Functions
def	__init__ (self, argv, useCool, useTestBeam, verbose, dbstring=None, putDuplicates=False, hvonly=False)
def	parseVarExpression (self, varExpression)
def	getTree (self, drawer=None, var=None, lastOnly=False, firstAndLast=False)
def	cutReplace (self, incut, var)
def	scan (self)
def	scantree (self, t, treevar, cut)
def	getTimelinePlot (self)
def	getDiffPlot (self, opt)
def	getDistributionPlot (self)
def	printDistributionStatAll (self)
def	printDistributionStat (self, lastOnly=False)

Public Attributes
	useCool
	useTestBeam
	cmd
	drawer
	varExp
	iovBeg
	iovEnd
	dataGrabber
	info
	cutExp
	outName
	rangeStr

Detailed Description

Definition at line 16 of file TileDCSDataPlotter.py.

Constructor & Destructor Documentation

__init__()

def TileDCSDataPlotter.TileDCSDataPlotter.__init__	(		self,
			argv,
			useCool,
			useTestBeam,
			verbose,
			dbstring = None,
			putDuplicates = False,
			hvonly = False
	)

Definition at line 19 of file TileDCSDataPlotter.py.

    def __init__( self, argv, useCool, useTestBeam, verbose, dbstring=None, putDuplicates=False, hvonly=False ):
 
        dbSource = "ORACLE"
        self.useCool = useCool or hvonly
        if useCool:
            dbSource = "COOL"
        self.useTestBeam = useTestBeam
        if useTestBeam:
            dbSource = "TESTBEAM"
        logLvl = logging.WARNING
        if verbose:
            logLvl = logging.DEBUG
 
#        self.dataGrabber = TileDCSDataGrabber(dbSource, logLvl, dbstring, putDuplicates)
#        self.info = self.dataGrabber.info
 
        self.cmd     = argv[1]
        self.drawer  = argv[2]
        self.varExp  = argv[3]
        beg          = argv[4]
        end          = argv[5]
        self.iovBeg  = int(time.mktime(time.strptime(beg,"%Y-%m-%d %H:%M:%S")))
        self.iovEnd  = int(time.mktime(time.strptime(end,"%Y-%m-%d %H:%M:%S")))
 
        self.dataGrabber = TileDCSDataGrabber(dbSource, logLvl, dbstring, putDuplicates, self.iovBeg)
        self.info = self.dataGrabber.info
 
        self.cutExp  = ""
        if len(argv) >6:
            self.cutExp = argv[6]
        beg = beg.replace(' ','_')
        end = end.replace(' ','_')
        self.outName   = self.drawer+"_"+self.varExp[:100]+"__"+beg+"__"+end
        self.outName   = self.outName.replace('/',':')
        self.outName   = self.outName.replace('(','0')
        self.outName   = self.outName.replace(')','0')
        if len(argv) >7:
            self.outName = argv[7]
        print ("---> OUTNAME: " , self.outName)
 
        timeBegInfo = time.localtime(self.iovBeg)
        timeEndInfo = time.localtime(self.iovEnd)
        self.rangeStr = "FROM: "+time.asctime(timeBegInfo)+"  UNTIL: "+time.asctime(timeEndInfo)
 
        part=["LBA","LBC","EBA","EBC"]
        if self.drawer == "ALL":
            self.drawer=""
            for p in range(4):
                for m in range(64):
                    self.drawer+="%s%2.2d," % (part[p],m+1)
            self.drawer=self.drawer[:-1]
        elif self.drawer in part:
            p=self.drawer
            self.drawer=""
            for m in range(64):
                self.drawer+="%s%2.2d," % (p,m+1)
            self.drawer=self.drawer[:-1]
 
 

Member Function Documentation

cutReplace()

def TileDCSDataPlotter.TileDCSDataPlotter.cutReplace	(		self,
			incut,
			var
	)

Relace "ALL" keywords in cut variable by appropriate names from var variable

Definition at line 244 of file TileDCSDataPlotter.py.

    def cutReplace(self, incut, var):
        """
        Relace "ALL" keywords in cut variable by appropriate names from var variable
 
        """
 
        cut = incut
        for vv in var.split("-"):
            if "hvOut" in vv:
                if "Set" in vv:
                    cut = cut.replace("ALL_SETHV",vv)
                    cut = cut.replace("ALL_HVSET",vv)
                    vvv=vv.replace("Set.hv","hv")
                    cut = cut.replace("ALL_HV",vvv)
                else:
                    vvv=vv.replace("hv","Set.hv")
                    cut = cut.replace("ALL_SETHV",vvv)
                    cut = cut.replace("ALL_HVSET",vvv)
                    cut = cut.replace("ALL_HV",vv)
        cut = cut.replace("ALL",var)
        return cut
 
 

getDiffPlot()

def TileDCSDataPlotter.TileDCSDataPlotter.getDiffPlot	(		self,
			opt
	)

Returns a TCanvas with 2D plot showing difference between last and first entries.

Definition at line 440 of file TileDCSDataPlotter.py.

    def getDiffPlot(self,opt):
        """
        Returns a TCanvas with 2D plot showing difference between last and first entries.
 
        """
 
        ROOT.gStyle.SetOptStat(0)
        tree = self.getTree(firstAndLast=(opt==1))
 
        if tree.GetEntries()==1:
            print ("ERROR: Only one data point, need at least 2")
            print ("---> Try increasing the time span")
            sys.exit(1)
 
        if len(self.cutExp) and "ALL" not in self.cutExp:
            print ("Cut expression is ignored in diff plot")
 
        #=== extract the values
        leaves=tree.GetListOfLeaves()
        xaxis=[]
        yaxis=[]
        for leaf in leaves:
            z=leaf.GetTitle().split('.')
            if len(z)==2 or (len(z)==3 and z[1]=="Set"):
                x=z[0]
                y=z[1]
                if len(z)==3:
                    y+= "."+z[2]
                if not y[-2].isdigit() and y[-1].isdigit():
                    y=y[:-2]+'0'+y[-1]
                if x not in xaxis:
                    xaxis+=[x]
                if y not in yaxis:
                    yaxis+=[y]
        xaxis.sort()
        yaxis.sort()
        nx=len(xaxis)
        ny=len(yaxis)
        if opt==1:
            pref="Diff "
        elif opt==2:
            pref="Diff_Prof "
        elif opt==3:
            pref="Prof "
        else:
            pref=""
        if opt==1:
            hist=ROOT.TH2F("hist2D",pref+self.outName,nx,-0.5,nx-0.5,ny,-0.5,ny-0.5)
        else:
            hist=ROOT.TProfile2D("hist2D",pref+self.outName,nx,-0.5,nx-0.5,ny,-0.5,ny-0.5)
 
        XA=hist.GetXaxis()
        n=0
        for i in range(nx):
            if n:
                XA.SetBinLabel(i + 1, '')
            else:
                XA.SetBinLabel(i + 1, xaxis[i])
            n=1-n
 
        YA=hist.GetYaxis()
        for i in range(ny):
            YA.SetBinLabel(i + 1, yaxis[i])
 
        tree.GetEntry(0)
        val = {}
        for leaf in leaves:
            if opt>2:
                val[leaf.GetTitle()] = 0
            else:
                val[leaf.GetTitle()] = leaf.GetValue()
 
        ne=tree.GetEntries()
        if opt==1:
            n1=ne-1
        elif opt==2:
            n1=1
        else:
            n1=0
        bar = progressBar(n1,ne, 78)
        for n in range(n1,ne):
            tree.GetEntry(n)
            bar.update(n)
            for leaf in leaves:
                v1=val[leaf.GetTitle()]
                z=leaf.GetTitle().split('.')
                if len(z)==2 or (len(z)==3 and z[1]=="Set"):
                    x=z[0]
                    y=z[1]
                    if len(z)==3:
                        y+= "."+z[2]
                    if not y[-2].isdigit() and y[-1].isdigit():
                        y=y[:-2]+'0'+y[-1]
                    v2=leaf.GetValue()
                    if "ALL" in self.cutExp:
                        var=str(v2)
                        vardiff=str(v2-v1)
                        cut = self.cutExp.replace("ALL_DIFF",vardiff)
                        cut = self.cutReplace(cut,var)
                        ok = eval(cut)
                        #print (cut)
                        #print (ok)
                    else:
                        ok = True
                    if ok:
                        hist.Fill(xaxis.index(x),yaxis.index(y),v2-v1)
        bar.done()
 
        return hist
 
 

getDistributionPlot()

def TileDCSDataPlotter.TileDCSDataPlotter.getDistributionPlot

(

self

)

Returns a TCanvas with the variable distribution plotted (using weights)

Definition at line 552 of file TileDCSDataPlotter.py.

    def getDistributionPlot(self):
        """
        Returns a TCanvas with the variable distribution plotted (using weights)
 
        """
 
        t = self.getTree()
 
        hh=[]
        ih=0
        color = [2,4,3,5,6,7]
        xmin=1000000000
        xmax=-1000000000
        ymin=1000000000
        ymax=-1000000000
        hmin=1000000000
        hmax=-1000000000
        dim=0
        xtit=""
        ytit=""
        for var in sorted(self.varExp.split(",")):
            cut=self.cutExp
            dim = var.count(':')
            if dim>0:
                varX = var.split(":")[1]
                varY = var.split(":")[0]
                if len(cut):
                    cutY = self.cutReplace(cut,varY)
                    cut = self.cutReplace(cut,varX)
                    cut = "("+cut+") && ("+cutY+")"
            else:
                varX=var
                varY=""
                if len(cut):
                    cut = self.cutReplace(cut,varX)
            if cut!="":
                cut = "weight*("+cut+")"
            else:
                cut = "weight"
            #print (var)
            #print (cut)
 
            t.Draw(var,cut,"goff")
            h = ROOT.gDirectory.Get("htemp")
            h.SetXTitle(var.split(":")[0])
            h.SetTitle(self.rangeStr)
            h.SetName("TileDCSDataPlotter"+str(ih))
            #print ("n=",h.GetEntries())
 
            if h.GetEntries()>0:
                xmi = h.GetXaxis().GetBinLowEdge(1)
                xma = h.GetXaxis().GetBinLowEdge(h.GetNbinsX())+h.GetXaxis().GetBinWidth(h.GetNbinsX())
                if xmi<xmin:
                    xmin=xmi
                if xma>xmax:
                    xmax=xma
                if dim>0:
                    ymi = h.GetYaxis().GetBinLowEdge(1)
                    yma = h.GetYaxis().GetBinLowEdge(h.GetNbinsY())+h.GetYaxis().GetBinWidth(h.GetNbinsY())
                    if ymi<ymin:
                        ymin=ymi
                    if yma>ymax:
                        ymax=yma
                hmi = h.GetMinimum()
                hma = h.GetMaximum()
                if hmi<hmin:
                    hmin=hmi
                if hma>hmax:
                    hmax=hma
                xtit+=varX+","
                if dim>0:
                    ytit+=varY+","
 
            h.SetLineColor(color[ih%len(color)])
            h.SetMarkerColor(color[ih%len(color)])
 
            ih+=1
            hh+=[h]
 
        if len(hh)>1:
            if xmin>xmax:
                xmin=0
                xmax=0
            delta=(xmax-xmin)/20.
            if delta<=0:
                delta=0.5
            xmin-=delta
            xmax+=delta
            if ymin>ymax:
                ymin=0
                ymax=0
            delta=(ymax-ymin)/20.
            if delta<=0:
                delta=0.5
            ymin-=delta
            ymax+=delta
            if hmin>hmax:
                hmin=0
                hmax=1
            delta=(hmax-hmin)/20.
            if delta<=0:
                delta=0.5
            if hmin!=0:
                hmin-=delta
            hmax+=delta
            if dim>0:
                h0 = TH2D("hDCS","dummy",2,xmin,xmax,2,ymin,ymax)
                h0.SetYTitle(ytit[:-1])
            else:
                h0 = TH1D("hDCS","dummy",2,xmin,xmax)
            h0.SetXTitle(xtit[:-1])
            h0.SetTitle(self.rangeStr)
            h0.SetMaximum(hmax)
            if hmin!=0:
                h0.SetMinimum(hmin)
            hh=[h0]+hh
            ROOT.gStyle.SetOptStat(0)
 
        return hh
 
 

getTimelinePlot()

def TileDCSDataPlotter.TileDCSDataPlotter.getTimelinePlot

(

self

)

Returns a TCanvas with the variable plotted vs time.

Definition at line 337 of file TileDCSDataPlotter.py.

    def getTimelinePlot(self):
        """
        Returns a TCanvas with the variable plotted vs time.
 
        """
 
        ROOT.gStyle.SetOptStat(0)
        t = self.getTree()
 
        if t.GetEntries()==1:
            print ("ERROR: Only one data point, need at least 2")
            print ("---> Try increasing the time span")
            sys.exit(1)
 
        #=== extract the values
        hh=[]
        ih=0
        hmin=1000000000
        hmax=-1000000000
        color = [2,4,3,5,6,7]
        for var in self.varExp.split(","):
            cut=self.cutReplace(self.cutExp,var)
            #print (var)
            #print (cut)
            t.Draw(var+":EvTime",cut,"goff")
            n = t.GetSelectedRows()
            if len(cut)>0:
                print ("var=",var,"   \tcut=",cut,"   \tnpoints=",n)
            else:
                print ("var=",var,"   \tnpoints=",n)
            if n>0: # at least one point passed the cut
                x1 = t.GetV2()
                y1 = t.GetV1()
                x = TArrayD(n)
                y = TArrayD(n)
                for i in range(n):
                    x[i] = x1[i]
                    y[i] = y1[i]
 
                t.Draw(var+":EvTime","","goff") # the same but without cut
                n0 = t.GetSelectedRows()
                x0 = t.GetV2()
 
                #=== fix for **** root time convention and add end time
                offset = 788918400 # = 1.1.1995(UTC), the root offset
                xarr = TArrayD(n0+1)
                for i in range(n0):
                    xarr[i] = x0[i]-offset
                xarr[n0] = self.iovEnd-offset
                #=== create and fill histogram with values
                title = var
                h = TH1D("hDCS"+str(ih),title,n0,xarr.GetArray())
                for i in range(n):
                    h.Fill(x[i]-offset,y[i])
                    if y[i]>hmax:
                        hmax=y[i]
                    if y[i]<hmin:
                        hmin=y[i]
 
                #=== set time display
                sec_min = 60
                sec_hrs = sec_min*60
                sec_day = sec_hrs*24
                sec_mon = sec_day*30
                timeDiff = self.iovEnd-self.iovBeg
                h.GetXaxis().SetTimeFormat("%S")
                h.SetXTitle("time [seconds]")
                if timeDiff > 12*sec_mon:
                    h.GetXaxis().SetTimeFormat("%m/%y")
                    h.SetXTitle("time [month/year]")
                #elif timeDiff >3.33448*sec_mon:
                elif timeDiff > 6*sec_day:
                    h.GetXaxis().SetTimeFormat("%d/%m")
                    h.SetXTitle("time [day/month]")
                elif timeDiff > 2*sec_day:
                    h.GetXaxis().SetTimeFormat("%d(%H)")
                    h.SetXTitle("time [day(hour)]")
                elif timeDiff > 2*sec_hrs:
                    h.GetXaxis().SetTimeFormat("%H:%M")
                    h.SetXTitle("time [hour:min]")
                elif   timeDiff > 2*sec_min:
                    h.GetXaxis().SetTimeFormat("%M")
                    h.SetXTitle("time [min]")
                h.GetXaxis().SetTimeDisplay(1)
                h.SetLineColor(color[ih%len(color)])
 
                ih+=1
                hh+=[h]
 
        if len(hh)>1:
            delta=(hmax-hmin)/20.
            if delta<=0:
                delta=0.5
            hmin-=delta
            hmax+=delta
            for h in hh:
                h.SetMaximum(hmax)
                h.SetMinimum(hmin)
 
        return hh
 
 

getTree()

def TileDCSDataPlotter.TileDCSDataPlotter.getTree	(		self,
			drawer = None,
			var = None,
			lastOnly = False,
			firstAndLast = False
	)

Definition at line 149 of file TileDCSDataPlotter.py.

    def getTree(self, drawer=None, var=None,lastOnly=False,firstAndLast=False):
 
        #=== parse for variables needed
        varList = self.parseVarExpression(self.varExp+" "+self.cutExp)
 
        #=== get the variable
        if drawer is None:
            t = self.dataGrabber.getDCSTree(self.iovBeg, self.iovEnd, self.drawer.split(","), varList, lastOnly, firstAndLast)
        elif var is None:
            drl=[drawer]
            t = self.dataGrabber.getDCSTree(self.iovBeg, self.iovEnd, drl, varList, lastOnly, firstAndLast)
        else:
            drl=[drawer]
            varl=[var]
            t = self.dataGrabber.getDCSTree(self.iovBeg, self.iovEnd, drl, varl, lastOnly, firstAndLast)
 
        if t.GetEntries()==0:
            print ("ERROR: No data in time interval!")
            sys.exit(1)
        else:
            print ("Found number of entries: ", t.GetEntries())
 
        #=== append drawer to all variables
        if self.cmd!="tree" and drawer is None and var is None:
            for var in varList:
                if "," in self.drawer:
                    newvar=""
                    newcut=""
                    if var in self.cutExp:
                        cut=self.cutExp
                    else:
                        cut=None
                    for dr in self.drawer.split(","):
                        newvar += dr+"."+var+","
                        if cut:
                            newcut += dr+"."+cut+","
                    self.varExp = self.varExp.replace(var,newvar[:-1])
                    if cut:
                        self.cutExp = newcut[:-1]
                    self.varExp = self.varExp.replace("Set."+newvar[:-1],"Set."+var)
                    if cut and "Set." in cut and "Set." not in var:
                        self.cutExp=cut
                else:
                    self.varExp = self.varExp.replace(var,self.drawer+"."+var)
                    self.cutExp = self.cutExp.replace(var,self.drawer+"."+var)
                    self.varExp = self.varExp.replace(self.drawer+"."+self.drawer,self.drawer)
                    self.cutExp = self.cutExp.replace(self.drawer+"."+self.drawer,self.drawer)
                    self.varExp = self.varExp.replace("Set."+self.drawer,"Set")
                    self.cutExp = self.cutExp.replace("Set."+self.drawer,"Set")
            if "ALL_LVPS_AI" in self.varExp:
                newvar=""
                for dr in self.drawer.split(","):
                    for var in list(self.info.vars_LVPS_AI.keys()):
                        newvar += dr+"."+var+","
                self.varExp = self.varExp.replace("ALL_LVPS_AI",newvar[:-1])
 
            if "ALL_LVPS_STATES" in self.varExp:
                newvar=""
                for dr in self.drawer.split(","):
                    for var in list(self.info.vars_LVPS_STATES.keys()):
                        newvar += dr+"."+var+","
                self.varExp = self.varExp.replace("ALL_LVPS_STATES",newvar[:-1])
 
            if "ALL_HV-ALL_SETHV" in self.varExp or "ALL_HV-ALL_HVSET" in self.varExp:
                newvar=""
                for dr in self.drawer.split(","):
                    for var in list(self.info.vars_HVSET.keys()):
                        if not ("Set.vFix" in var or "Set.hvIn" in var or "Set.volt" in var or "Set.temp" in var):
                            newvar += dr+"."+var.replace("Set.","")+"-"+dr+"."+var+","
                self.varExp = self.varExp.replace("ALL_HV-ALL_SETHV",newvar[:-1])
                self.varExp = self.varExp.replace("ALL_HV-ALL_HVSET",newvar[:-1])
            if "ALL_HVSET" in self.varExp or "ALL_SETHV" in self.varExp:
                newvar=""
                for dr in self.drawer.split(","):
                    for var in list(self.info.vars_HVSET.keys()):
                        if not self.useCool or not ("Set.vFix" in var or "Set.hvIn" in var or "Set.volt" in var or "Set.temp" in var):
                            newvar += dr+"."+var+","
                self.varExp = self.varExp.replace("ALL_HVSET",newvar[:-1])
                self.varExp = self.varExp.replace("ALL_SETHV",newvar[:-1])
            if "ALL_HV" in self.varExp and "ALL_HVSET" not in self.varExp :
                newvar=""
                for dr in self.drawer.split(","):
                    for var in list(self.info.vars_HV.keys()):
                        if not self.useCool or not ("hvIn" in var or "volt" in var or "temp" in var):
                            newvar += dr+"."+var+","
                self.varExp = self.varExp.replace("ALL_HV",newvar[:-1])
 
            print ("self.drawer: ", self.drawer)
            print ("self.varExp: ", self.varExp)
            print ("self.cutExp: ", self.cutExp)
 
        return t
 
 

parseVarExpression()

def TileDCSDataPlotter.TileDCSDataPlotter.parseVarExpression	(		self,
			varExpression
	)

Definition at line 79 of file TileDCSDataPlotter.py.

    def parseVarExpression(self, varExpression):
 
        knownVars = self.info.get_all_variables()
        varDict = {}
        for var in knownVars:
            pos = varExpression.find(var)
            if pos>-1:
                oldVar=""
                if pos in varDict:
                    #=== other variable found at same place?
                    #=== -> choose the longer one
                    oldVar = varDict[pos]
                    if len(oldVar) > len(var):
                        vtmp = var
                        var = oldVar
                        oldVar = vtmp
                varDict[pos] = var
                if oldVar!="":
                    #=== check if shorter variable is also present at another place
                    varExpr = varExpression
                    for v in list(varDict.values()):
                        varExpr=varExpr.replace(v," "*len(v))
                    oldPos = varExpr.find(oldVar)
                    if oldPos>-1:
                        varDict[oldPos] = oldVar
 
        #=== check for overlap
        positions = sorted(varDict.keys())
        posLength = len(positions)
        if posLength > 1:
            iPos = 1
            while iPos < posLength:
                if positions[iPos] < positions[iPos-1]+len(varDict[positions[iPos-1]]):
                    varDict.pop(positions[iPos])
                    positions.pop(iPos)
                    posLength=len(positions)
                else:
                    iPos+=1
 
        varList = list(varDict.values())
        if "ALL_LVPS_AI" in varExpression:
            varList.extend(list(self.info.vars_LVPS_AI.keys()))
        if "ALL_LVPS_STATES" in varExpression:
            varList.extend(list(self.info.vars_LVPS_STATES.keys()))
        if "ALL_HVSET" in varExpression or "ALL_SETHV" in varExpression:
            varList.extend(list(self.info.vars_HVSET.keys()))
            if self.useCool or ("ALL_HV-ALL_SETHV" in varExpression or "ALL_HV-ALL_HVSET" in varExpression):
                for i in range(2):
                    if "Set.vFix1%d" % (i+1) in varList:
                        varList.remove("Set.vFix1%d" % (i+1))
                for i in range(4):
                    varList.remove("Set.hvIn%d" % (i+1))
                for i in range(7):
                    varList.remove("Set.volt%d" % (i+1))
                for i in range(7):
                    varList.remove("Set.temp%d" % (i+1))
        if ("ALL_HV" in varExpression and "ALL_HVSET" not in varExpression) or "ALL_HV-ALL_HVSET" in varExpression:
            varList.extend(list(self.info.vars_HV.keys()))
            if self.useCool or ("ALL_HV-ALL_SETHV" in varExpression or "ALL_HV-ALL_HVSET" in varExpression):
                for i in range(4):
                    varList.remove("hvIn%d" % (i+1))
                for i in range(7):
                    varList.remove("volt%d" % (i+1))
                for i in range(7):
                    varList.remove("temp%d" % (i+1))
 
        return list(set(varList))
 
 

printDistributionStat()

def TileDCSDataPlotter.TileDCSDataPlotter.printDistributionStat	(		self,
			lastOnly = False
	)

Print mean values from histograms

Definition at line 700 of file TileDCSDataPlotter.py.

    def printDistributionStat(self, lastOnly=False):
        """
        Print mean values from histograms
 
        """
 
        if "ALL_LVPS_AI" in self.varExp:
            self.varExp = self.varExp.replace("ALL_LVPS_AI",",".join(sorted(self.info.vars_LVPS_AI.keys())))
        if "ALL_LVPS_STATES" in self.varExp:
            self.varExp = self.varExp.replace("ALL_LVPS_STATES",",".join(sorted(self.info.vars_LVPS_STATES.keys())))
        if "ALL_HV-ALL_SETHV" in self.varExp or "ALL_HV-ALL_HVSET" in self.varExp:
            vlist = []
            for var in list(self.info.vars_HVSET.keys()):
                if not ("Set.vFix" in var or "Set.hvIn" in var or "Set.volt" in var or "Set.temp" in var):
                    vlist += [var.replace("Set.","")+"-"+var]
            self.varExp = self.varExp.replace("ALL_HV-ALL_SETHV",",".join(sorted(vlist)))
            self.varExp = self.varExp.replace("ALL_HV-ALL_HVSET",",".join(sorted(vlist)))
        if "ALL_HVSET" in self.varExp or "ALL_SETHV" in self.varExp:
            vlist = list(self.info.vars_HVSET.keys())
            if self.useCool:
                for i in range(2):
                    if "Set.vFix1%d" % (i+1) in vlist:
                        vlist.remove("Set.vFix1%d" % (i+1))
                for i in range(4):
                    vlist.remove("Set.hvIn%d" % (i+1))
                for i in range(7):
                    vlist.remove("Set.volt%d" % (i+1))
                for i in range(7):
                    vlist.remove("Set.temp%d" % (i+1))
            self.varExp = self.varExp.replace("ALL_HVSET",",".join(sorted(vlist)))
            self.varExp = self.varExp.replace("ALL_SETHV",",".join(sorted(vlist)))
        if "ALL_HV" in self.varExp and "ALL_HVSET" not in self.varExp:
            vlist = list(self.info.vars_HV.keys())
            if self.useCool:
                for i in range(4):
                    vlist.remove("hvIn%d" % (i+1))
                for i in range(7):
                    vlist.remove("volt%d" % (i+1))
                for i in range(7):
                    vlist.remove("temp%d" % (i+1))
            self.varExp = self.varExp.replace("ALL_HV",",".join(sorted(vlist)))
        varlist=self.parseVarExpression(self.varExp)
        print ("var list is",varlist)
 
        reslist=[]
        cutlist = self.parseVarExpression(self.cutExp)
        print ("cut list is",cutlist)
        fulllist=sorted(list(set(varlist+cutlist)))
        print ("full list is",fulllist)
        t = None
        drlistprev=[]
        varlistprev=[]
        for drawer in self.drawer.split(","):
            drlist=[drawer]
            for v in self.varExp.split(","):
                varlist=self.parseVarExpression(v)
                varlist += cutlist
                varlist=sorted(list(set(varlist)))
                #for drawer in self.drawer.split(","):
                #    drlist=[drawer]
                if not (t and drlist==drlistprev and varlist==varlistprev):
                    t = self.dataGrabber.getDCSTree(self.iovBeg, self.iovEnd, drlist, varlist, lastOnly)
                    drlistprev=drlist
                    varlistprev=varlist
                if t and t.GetEntries()>0:
                    var = v
                    cut = self.cutReplace(self.cutExp,var)
                    for va in varlist:
                        var = var.replace(va,drawer+"."+va)
                        cut = cut.replace(va,drawer+"."+va)
                    var = var.replace(drawer+"."+drawer,drawer)
                    cut = cut.replace(drawer+"."+drawer,drawer)
                    var = var.replace("Set."+drawer,"Set")
                    cut = cut.replace("Set."+drawer,"Set")
                    print ("var is",var)
                    print ("cut is",cut)
                    t.Draw(var,cut,"goff")
                    h = ROOT.gDirectory.Get("htemp")
                    if h:
                        if cut!="":
                            reslist.append((var+" CUT "+cut,self.rangeStr,h.GetEntries(),h.GetMean(),h.GetRMS()))
                        else:
                            reslist.append((var,self.rangeStr,h.GetEntries(),h.GetMean(),h.GetRMS()))
                    else:
                        reslist.append((var,self.rangeStr,-1,0.0,0.0))
 
        if lastOnly:
            for res in reslist:
                print ("%-25s  %s        Last %-9.3f" % (res[0],res[1],res[3]))
        else:
            for res in reslist:
                print ("%-25s  %s        Nentries %-7d  Mean %-9.3f  RMS %-9.4f" % res)
        return 0
 
 
#================================================================================================
 
 
#==== interactive use

printDistributionStatAll()

def TileDCSDataPlotter.TileDCSDataPlotter.printDistributionStatAll

(

self

)

Print mean values from histograms

Definition at line 674 of file TileDCSDataPlotter.py.

    def printDistributionStatAll(self):
        """
        Print mean values from histograms
 
        """
 
        t = self.getTree()
        if self.cutExp!="":
            for v, c in zip(self.varExp.split(","), self.cutExp.split(",")):
                t.Draw(v,c,"goff")
                h = ROOT.gDirectory.Get("htemp")
 
                print ("%s  %s  %s  Nentries %d  Mean %7.3f  RMS %8.4f" % \
                    (v,c,self.rangeStr,h.GetEntries(),h.GetMean(),h.GetRMS()))
        else:
            cut = ""
            for v in self.varExp.split(","):
                t.Draw(v,cut,"goff")
                h = ROOT.gDirectory.Get("htemp")
                print ("%s  %s  Nentries %d  Mean %7.3f  RMS %8.4f" % \
                    (v,self.rangeStr,h.GetEntries(),h.GetMean(),h.GetRMS()))
 
        return 0
 
 

scan()

def TileDCSDataPlotter.TileDCSDataPlotter.scan

(

self

)

Scan ROOT TTree and show all values on the screen

Definition at line 268 of file TileDCSDataPlotter.py.

    def scan(self):
        """
        Scan ROOT TTree and show all values on the screen
 
        """
 
        t = self.getTree()
        if "ALL" in self.cutExp or self.varExp.count(',')>2 :
            #=== extract the values
            for var in self.varExp.split(","):
                treevar="EvTime:"+var
                cut=self.cutReplace(self.cutExp,var)
                self.scantree(t,treevar,cut)
        else:
            treevar="EvTime:"+self.varExp.replace(",",":")
            self.scantree(t,treevar,self.cutExp)
 
        return
 
 

scantree()

def TileDCSDataPlotter.TileDCSDataPlotter.scantree	(		self,
			t,
			treevar,
			cut
	)

Scan ROOT TTree and show up to 4 variables on the screen

Definition at line 289 of file TileDCSDataPlotter.py.

    def scantree(self,t,treevar,cut):
        """
        Scan ROOT TTree and show up to 4 variables on the screen
 
        """
 
        #t.Scan("-1000000000+"+treevar,cut)
        t.Draw(treevar,cut,"goff")
        n = t.GetSelectedRows()
        if n>0: # at least one point passed the cut
            vars = treevar.split(":")
            nv = len(vars)
            fmt = [ "%d" ]
            for i in range(1,nv):
                if "DAQ" in vars[i]:
                    fmt += [ "\t%6d\t\t" ]
                else:
                    fmt += [ "\t%6.2f\t" ]
            x = t.GetV1()
            v1 = None
            v2 = None
            v3 = None
            v=vars[0]
            if v=="EvTime": v += 13*" "
            if nv>1:
                v1 = t.GetV2()
                v += "\t"+vars[1]
                if nv>2:
                    v2 = t.GetV3()
                    v += "\t"+vars[2]
                    if nv>3:
                        v3 = t.GetV4()
                        v += "\t"+vars[3]
            print(v)
            for i in range(n):
                v=""
                if v1 is not None:
                    v += fmt[1] % v1[i]
                    if v2 is not None:
                        v += fmt[2] % v2[i]
                        if v3 is not None:
                            v += fmt[3] % v3[i]
                dt = datetime.datetime.fromtimestamp(int(x[i])).strftime('%Y-%m-%d %H:%M:%S')
                print(dt,v)
        return
 
 

Member Data Documentation

cmd

TileDCSDataPlotter.TileDCSDataPlotter.cmd

Definition at line 35 of file TileDCSDataPlotter.py.

cutExp

TileDCSDataPlotter.TileDCSDataPlotter.cutExp

Definition at line 46 of file TileDCSDataPlotter.py.

dataGrabber

TileDCSDataPlotter.TileDCSDataPlotter.dataGrabber

Definition at line 43 of file TileDCSDataPlotter.py.

drawer

TileDCSDataPlotter.TileDCSDataPlotter.drawer

Definition at line 36 of file TileDCSDataPlotter.py.

info

TileDCSDataPlotter.TileDCSDataPlotter.info

Definition at line 44 of file TileDCSDataPlotter.py.

iovBeg

TileDCSDataPlotter.TileDCSDataPlotter.iovBeg

Definition at line 40 of file TileDCSDataPlotter.py.

iovEnd

TileDCSDataPlotter.TileDCSDataPlotter.iovEnd

Definition at line 41 of file TileDCSDataPlotter.py.

outName

TileDCSDataPlotter.TileDCSDataPlotter.outName

Definition at line 51 of file TileDCSDataPlotter.py.

rangeStr

TileDCSDataPlotter.TileDCSDataPlotter.rangeStr

Definition at line 61 of file TileDCSDataPlotter.py.

useCool

TileDCSDataPlotter.TileDCSDataPlotter.useCool

Definition at line 22 of file TileDCSDataPlotter.py.

useTestBeam

TileDCSDataPlotter.TileDCSDataPlotter.useTestBeam

Definition at line 25 of file TileDCSDataPlotter.py.

varExp

TileDCSDataPlotter.TileDCSDataPlotter.varExp

Definition at line 37 of file TileDCSDataPlotter.py.

---

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

• TileDCSDataPlotter.py