TileDCSDataPlotter.TileDCSDataPlotter Class Reference

Inheritance diagram for TileDCSDataPlotter.TileDCSDataPlotter:

Collaboration diagram for TileDCSDataPlotter.TileDCSDataPlotter:

Public Types
typedef HLT::TypeInformation::for_each_type_c< typenameEDMLIST::map, my_functor, my_result<>, my_arg< HLT::TypeInformation::get_cont, CONTAINER > >::type	result

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

Public Attributes
	useCool = useCool or hvonly
	useTestBeam = useTestBeam
	cmd = argv[1]
str	drawer = argv[2]
	varExp = argv[3]
	iovBeg = int(time.mktime(time.strptime(beg,"%Y-%m-%d %H:%M:%S")))
	iovEnd = int(time.mktime(time.strptime(end,"%Y-%m-%d %H:%M:%S")))
	dataGrabber = TileDCSDataGrabber(dbSource, logLvl, dbstring, putDuplicates, self.iovBeg)
	info = self.dataGrabber.info
str	cutExp = ""
str	outName = self.drawer+"_"+self.varExp[:100]+"__"+beg+"__"+end
str	rangeStr = "FROM: "+time.asctime(timeBegInfo)+" UNTIL: "+time.asctime(timeEndInfo)

Detailed Description

Definition at line 15 of file TileDCSDataPlotter.py.

Member Typedef Documentation

result

typedef HLT::TypeInformation::for_each_type_c<typenameEDMLIST::map,my_functor,my_result<>,my_arg<HLT::TypeInformation::get_cont,CONTAINER>>::type master_search< typenameEDMLIST::map, HLT::TypeInformation::get_cont, CONTAINER >::result

inherited

Definition at line 90 of file EDM_MasterSearch.h.

Constructor & Destructor Documentation

__init__()

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

Definition at line 18 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()

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

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

Definition at line 243 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()

TileDCSDataPlotter.TileDCSDataPlotter.getDiffPlot	(		self,
			opt )

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

Definition at line 439 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()

TileDCSDataPlotter.TileDCSDataPlotter.getDistributionPlot

(

self

)

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

Definition at line 551 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()

TileDCSDataPlotter.TileDCSDataPlotter.getTimelinePlot

(

self

)

Returns a TCanvas with the variable plotted vs time.

Definition at line 336 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()

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

Definition at line 148 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()

TileDCSDataPlotter.TileDCSDataPlotter.parseVarExpression	(		self,
			varExpression )

Definition at line 78 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()

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

Print mean values from histograms

Definition at line 699 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()

TileDCSDataPlotter.TileDCSDataPlotter.printDistributionStatAll

(

self

)

Print mean values from histograms

Definition at line 673 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()

TileDCSDataPlotter.TileDCSDataPlotter.scan

(

self

)

Scan ROOT TTree and show all values on the screen

Definition at line 267 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()

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

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

Definition at line 288 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 = argv[1]

Definition at line 34 of file TileDCSDataPlotter.py.

cutExp

TileDCSDataPlotter.TileDCSDataPlotter.cutExp = ""

Definition at line 45 of file TileDCSDataPlotter.py.

dataGrabber

TileDCSDataPlotter.TileDCSDataPlotter.dataGrabber = TileDCSDataGrabber(dbSource, logLvl, dbstring, putDuplicates, self.iovBeg)

Definition at line 42 of file TileDCSDataPlotter.py.

drawer

TileDCSDataPlotter.TileDCSDataPlotter.drawer = argv[2]

Definition at line 35 of file TileDCSDataPlotter.py.

info

TileDCSDataPlotter.TileDCSDataPlotter.info = self.dataGrabber.info

Definition at line 43 of file TileDCSDataPlotter.py.

iovBeg

TileDCSDataPlotter.TileDCSDataPlotter.iovBeg = int(time.mktime(time.strptime(beg,"%Y-%m-%d %H:%M:%S")))

Definition at line 39 of file TileDCSDataPlotter.py.

iovEnd

TileDCSDataPlotter.TileDCSDataPlotter.iovEnd = int(time.mktime(time.strptime(end,"%Y-%m-%d %H:%M:%S")))

Definition at line 40 of file TileDCSDataPlotter.py.

outName

TileDCSDataPlotter.TileDCSDataPlotter.outName = self.drawer+"_"+self.varExp[:100]+"__"+beg+"__"+end

Definition at line 50 of file TileDCSDataPlotter.py.

rangeStr

TileDCSDataPlotter.TileDCSDataPlotter.rangeStr = "FROM: "+time.asctime(timeBegInfo)+" UNTIL: "+time.asctime(timeEndInfo)

Definition at line 60 of file TileDCSDataPlotter.py.

useCool

TileDCSDataPlotter.TileDCSDataPlotter.useCool = useCool or hvonly

Definition at line 21 of file TileDCSDataPlotter.py.

useTestBeam

TileDCSDataPlotter.TileDCSDataPlotter.useTestBeam = useTestBeam

Definition at line 24 of file TileDCSDataPlotter.py.

varExp

TileDCSDataPlotter.TileDCSDataPlotter.varExp = argv[3]

Definition at line 36 of file TileDCSDataPlotter.py.

---

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

• TileDCSDataPlotter.py