lumiPlot.py

Go to the documentation of this file.

#!/usr/bin/env python
#########################################################################
# Command to produce plots with selected conditions and styles          #
# output                                                                # 
#                                                                       #
# Author:      Zhen Xie                                                 #
#########################################################################
#
import os,os.path,sys,datetime,time,csv,re
from RecoLuminosity.LumiDB import argparse,lumiTime,matplotRender,sessionManager,lumiParameters,RegexValidator,revisionDML,lumiCalcAPI,normDML
import matplotlib
from matplotlib.figure import Figure

def parseInputFiles(inputfilename):
    '''
    output ({run:[cmsls,cmsls,...]},[[resultlines]])
    '''
    selectedrunlsInDB={}
    resultlines=[]
    p=inputFilesetParser.inputFilesetParser(inputfilename)
    runlsbyfile=p.runsandls()
    selectedProcessedRuns=p.selectedRunsWithresult()
    selectedNonProcessedRuns=p.selectedRunsWithoutresult()
    resultlines=p.resultlines()
    for runinfile in selectedNonProcessedRuns:
        selectedrunlsInDB[runinfile]=runlsbyfile[runinfile]
    return (selectedrunlsInDB,resultlines)
##############################
## ######################## ##
## ## ################## ## ##
## ## ## Main Program ## ## ##
## ## ################## ## ##
## ######################## ##
##############################

if __name__=='__main__':
    referenceLabel='Delivered'#from which label,lumi unit should be calculated
    labels=['Delivered','Recorded']#default label order
    allowedActions = ['run','time','fill','perday','instpeakperday']
    beamChoices=['PROTPHYS','IONPHYS','PAPHYS']
    allowedscales=['linear','log','both']
    beamModeChoices = [ "stable", "quiet", "either"]
    amodetagChoices = [ "PROTPHYS","IONPHYS",'PAPHYS' ]
    actiontofilebasemap={'time':'lumivstime','run':'lumivsrun','fill':'lumivsfill','perday':'lumiperday','instpeakperday':'lumipeak'}
    #############################################################
    # parse figure configuration if found
    #############################################################
    currentdir=os.getcwd()
    rcparamfile='.lumiplotrc'
    mplrcdir=matplotlib.get_configdir()
    mpllumiconfig=os.path.join(mplrcdir,rcparamfile)
    locallumiconfig=os.path.join(currentdir,rcparamfile)
    figureparams={'sizex':7.5,'sizey':5.7,'dpi':135}
    if os.path.exists(locallumiconfig):
        import ConfigParser
        try:
            config=ConfigParser.RawConfigParser()
            config.read(locallumiconfig)
            figureparams['sizex']=config.getfloat('sizex')
            figureparams['sizey']=config.getfloat('sizey')
            figureparams['dpi']=config.getint('dpi')
        except ConfigParser.NoOptionError:
            pass
    elif os.path.exists(mpllumiconfig):
        import ConfigParser
        try:
            config=ConfigParser.RawConfigParser()
            config.read(mpllumiconfig)
            figureparams['sizex']=config.getfloat('sizex')
            figureparams['sizey']=config.getfloat('sizey')
            figureparams['dpi']=config.getint('dpi')
        except ConfigParser.NoOptionError:
            pass
    ###############################
    # parse arguments
    ###############################  
    #
    # basic arguments
    #
    parser = argparse.ArgumentParser(prog=os.path.basename(sys.argv[0]),
                                     description="Plot luminosity as function of the time variable of choice",
                                     formatter_class=argparse.ArgumentDefaultsHelpFormatter)
    parser.add_argument('-c',
                        dest='connect',
                        action='store',
                        help='connect string to lumiDB',
                        default='frontier://LumiCalc/CMS_LUMI_PROD')
    parser.add_argument('-P',
                        dest='authpath',
                        action='store',
                        help='path to authentication file')
    parser.add_argument('--datatag',
                        dest='datatag',
                        action='store',
                        default=None,
                        help='data tag')
    parser.add_argument('--normtag',
                        dest='normtag',
                        action='store',
                        default=None,
                        help='norm tag')
    #
    #optional arg to select exact run and ls
    #
    parser.add_argument('-i',
                        dest='inputfile',
                        action='store',
                        help='run/ls selection file')
    #
    #optional arg to select exact hltpath or pattern
    #
    parser.add_argument('--hltpath',
                        dest='hltpath',
                        action='store',
                        default=None,
                        required=False,
                        help='specific hltpath or hltpath pattern to calculate the effectived luminosity (optional)')
    #
    #optional args to filter run/ls
    #
    parser.add_argument('-b',
                        dest='beamstatus',
                        action='store',
                        help='selection criteria beam status')
    #
    #optional args to filter *runs*, they do not select on LS level.
    #
    parser.add_argument('-f','--fill',
                        dest='fillnum',
                        action='store',
                        type=int,
                        help='specific fill',
                        default=None)
    parser.add_argument('--amodetag',
                        dest='amodetag',
                        action='store',
                        help='accelerator mode')
    parser.add_argument('--beamenergy',
                        dest='beamenergy',
                        action='store',
                        type=float,
                        default=3500,
                        help='beamenergy (in GeV) selection criteria')
    parser.add_argument('--beamfluctuation',
                        dest='beamfluctuation',
                        action='store',
                        type=float,
                        default=0.2,
                        help='allowed fraction of beamenergy to fluctuate')
    parser.add_argument('--begin',
                        dest='begintime',
                        action='store',
                        default='03/01/10 00:00:00',
                        type=RegexValidator.RegexValidator("^\d\d/\d\d/\d\d \d\d:\d\d:\d\d$","must be form mm/dd/yy hh:mm:ss"),
                        help='min run start time,mm/dd/yy hh:mm:ss')
    parser.add_argument('--end',
                        dest='endtime',
                        action='store',
                        type=RegexValidator.RegexValidator("^\d\d/\d\d/\d\d \d\d:\d\d:\d\d$","must be form mm/dd/yy hh:mm:ss"),
                        help='max run start time,mm/dd/yy hh:mm:ss')
    #
    #frontier config options
    #
    parser.add_argument('--siteconfpath',
                        dest='siteconfpath',
                        action='store',
                        help='specific path to site-local-config.xml file, default to $CMS_PATH/SITECONF/local/JobConfig, if path undefined, fallback to cern proxy&server')
    #
    #plot options
    #
    ####input/output file names
    parser.add_argument('--inplotdata',
                        dest='inplot',
                        action='store',
                        help='existing base plot(s) in text format')
    parser.add_argument('--outplotdata',
                        dest='outplot',
                        action='store',
                        help='output plot. By default, a text dump of the plot is produced.')
    ####graphic switches
    parser.add_argument('--with-annotation',
                        dest='withannotation',
                        action='store_true',
                        help='annotate boundary run numbers')
    parser.add_argument('--interactive',
                        dest='interactive',
                        action='store_true',
                        help='graphical mode to draw plot in a QT pannel.')
    parser.add_argument('--without-textoutput',
                        dest='withoutTextoutput',
                        action='store_true',
                        help='not to write out text output file')
    parser.add_argument('--without-png',
                        dest='withoutpng',
                        action='store_true',
                        help='do not produce')
    parser.add_argument('--yscale',
                        dest='yscale',
                        action='store',
                        default='linear',
                        help='y_scale[linear,log,both]')
    parser.add_argument('--lastpointfromdb',
                        dest='lastpointfromdb',
                        action='store_true',
                        help='the last point in the plot is always recaculated from current values in LumiDB')
    
    ####correction switches
    parser.add_argument('--without-correction',
                        dest='withoutNorm',
                        action='store_true',
                        help='without any correction/calibration')    

    ###general switches
    parser.add_argument('--verbose',
                        dest='verbose',
                        action='store_true',
                        help='verbose mode, print result also to screen')
    parser.add_argument('--debug',
                        dest='debug',
                        action='store_true',
                        help='debug')
    parser.add_argument('action',
                        choices=allowedActions,
                        help='type of plots')
    options=parser.parse_args()
    if options.yscale=='both' and options.interactive:
        print '--interactive mode can draw either yscale log or linear at a time'
        exit(0)
    outplotfilename = options.outplot
    if not outplotfilename:
        outplotfilename=actiontofilebasemap[options.action]
    outtextfilename = outplotfilename+'.csv'
    if options.withoutTextoutput:
        outtextfilename=None
    lumip=lumiParameters.ParametersObject()
    svc=sessionManager.sessionManager(options.connect,
                                      authpath=options.authpath,
                                      siteconfpath=options.siteconfpath,
                                      debugON=options.debug)
    session=svc.openSession(isReadOnly=True,cpp2sqltype=[('unsigned int','NUMBER(10)'),('unsigned long long','NUMBER(20)')])
    lslength=lumip.lslengthsec()
    begtime=options.begintime
    endtime=options.endtime
    lut=lumiTime.lumiTime()
    if not endtime:
        endtime=lut.DatetimeToStr(datetime.datetime.utcnow(),customfm='%m/%d/%y %H:%M:%S')
    pbeammode=None
    if options.beamstatus=='stable':
        pbeammode='STABLE BEAMS'
    resultlines=[]
    inplots=[]
    #######################################################################
    ##process old plot csv files,if any, skipping #commentlines
    #######################################################################
    inplot=options.inplot
    if inplot:
        inplots=inplot.split('+')
        for ip in inplots:
            f=open(ip,'r')
            reader=csv.reader(f,delimiter=',')
            for row in reader:
                if '#' in row[0]:continue
                resultlines.append(row)
    #####################################
    ##find runs need to read from DB
    #####################################
    irunlsdict={}#either from -i or from other selection options
    rruns=[]
    reqTrg=False
    reqHlt=False
    if options.hltpath:
        reqHlt=True
    fillrunMap={}
    lastDrawnRun=None
    maxDrawnDay=None
    lastDrawnFill=None
    newFirstDay=None
    newFirstRun=None
    newFirstFill=None
    session.transaction().start(True)
    if options.action=='perday' or options.action=='instpeakperday':
        maxDrawnDay=int(lut.StrToDatetime(begtime,customfm='%m/%d/%y %H:%M:%S').date().toordinal()) #default maxDrawnDay is begtime
        if resultlines: #if there's old result lines
            for drawnDay in [ int(t[0]) for t in resultlines]:
                if drawnDay>maxDrawnDay: #all required days are already drawn
                    maxDrawnDay=drawnDay
        newFirstDay=maxDrawnDay
        midnight=datetime.time()
        begT=datetime.datetime.combine(datetime.date.fromordinal(newFirstDay),midnight)
        begTStr=lut.DatetimeToStr(begT,customfm='%m/%d/%y %H:%M:%S')
        #find runs not in old plot
        runlist=lumiCalcAPI.runList(session.nominalSchema(),options.fillnum,runmin=None,runmax=None,startT=begTStr,stopT=endtime,l1keyPattern=None,hltkeyPattern=None,amodetag=options.amodetag,nominalEnergy=options.beamenergy,energyFlut=options.beamfluctuation,requiretrg=reqTrg,requirehlt=reqHlt)
    if options.action=='run' or options.action=='time':
        lastDrawnRun=132000
        if resultlines:#if there's old plot, start to count runs only after that
            lastDrawnRun=max([int(t[0]) for t in resultlines])
        newFirstRun=lastDrawnRun+1
        if options.lastpointfromdb:
            newFirstRun=lastDrawnRun
        runlist=lumiCalcAPI.runList(session.nominalSchema(),options.fillnum,runmin=newFirstRun,runmax=None,startT=begtime,stopT=endtime,l1keyPattern=None,hltkeyPattern=None,amodetag=options.amodetag,nominalEnergy=options.beamenergy,energyFlut=options.beamfluctuation,requiretrg=reqTrg,requirehlt=reqHlt)
    if options.action=='fill':
        lastDrawnFill=1000
        lastDrawnRun=132000
        if resultlines:
            lastDrawnFill=max([int(t[0]) for t in resultlines])        
            lastDrawnRun=min([int(t[1]) for t in resultlines if int(t[0])==lastDrawnFill])
        newFirstRun=lastDrawnRun+1
        if options.lastpointfromdb:
            newFirstRun=lastDrawnRun
        runlist=lumiCalcAPI.runList(session.nominalSchema(),runmin=newFirstRun,runmax=None,startT=begtime,stopT=endtime,l1keyPattern=None,hltkeyPattern=None,amodetag=options.amodetag,nominalEnergy=options.beamenergy,energyFlut=options.beamfluctuation,requiretrg=reqTrg,requirehlt=reqHlt)
        fillrunMap=lumiCalcAPI.fillrunMap(session.nominalSchema(),runmin=newFirstRun,runmax=None,startT=begtime,stopT=endtime)

    if options.inputfile:
        (irunlsdict,iresults)=parseInputFiles(options.inputfile)
        #apply further filter only if specified
        if options.fillnum or options.begin or options.end or options.amodetag or options.beamenergy:
            rruns=[val for val in runlist if val in irunlsdict.keys()]
            for selectedrun in irunlsdict.keys():#if there's further filter on the runlist,clean input dict
                if selectedrun not in rruns:
                    del irunlsdict[selectedrun]
        else:
            rruns=irunlsdict.keys()
    else:
        for run in runlist:
            irunlsdict[run]=None
        rruns=irunlsdict.keys()
            
    if options.verbose:
        print '[INFO] runs from db: ',irunlsdict
        if lastDrawnRun:
            print '[INFO] last run in old plot: ',lastDrawnRun
            print '[INFO] first run from DB in fresh plot: ',newFirstRun
        if maxDrawnDay:
            print '[INFO] last day in old plot: ',maxDrawnDay
            print '[INFO] first day from DB in fresh plot: ',newFirstDay

    if len(rruns)==0:
        if len(resultlines)!=0:
            print '[INFO] drawing all from old plot data'
        else:
            print '[INFO] found no old nor new data, do nothing'
            exit(0)

    GrunsummaryData=lumiCalcAPI.runsummaryMap(session.nominalSchema(),irunlsdict)
    #####################################
    #resolve data ids
    #####################################
    datatagname=options.datatag
    if not datatagname:
        (datatagid,datatagname)=revisionDML.currentDataTag(session.nominalSchema())
        dataidmap=revisionDML.dataIdsByTagId(session.nominalSchema(),datatagid,runlist=rruns,withcomment=False)
        #{run:(lumidataid,trgdataid,hltdataid,())}
    else:
        dataidmap=revisionDML.dataIdsByTagName(session.nominalSchema(),datatagname,runlist=rruns,withcomment=False)
    #
    # check normtag and get norm values if required
    #
    normname='NONE'
    normid=0
    normvalueDict={}
    if not options.withoutNorm:
        normname=options.normtag
        if not normname:
            normmap=normDML.normIdByType(session.nominalSchema(),lumitype='HF',defaultonly=True)
            if len(normmap):
                normname=normmap.keys()[0]
                normid=normmap[normname]
        else:
            normid=normDML.normIdByName(session.nominalSchema(),normname)
        if not normid:
            raise RuntimeError('[ERROR] cannot resolve norm/correction')
            sys.exit(-1)
        normvalueDict=normDML.normValueById(session.nominalSchema(),normid) #{since:[corrector(0),{paramname:paramvalue}(1),amodetag(2),egev(3),comment(4)]}

    fig=Figure(figsize=(figureparams['sizex'],figureparams['sizey']),dpi=figureparams['dpi'])
    m=matplotRender.matplotRender(fig)
    logfig=Figure(figsize=(figureparams['sizex'],figureparams['sizey']),dpi=figureparams['dpi'])
    mlog=matplotRender.matplotRender(logfig)

    if not options.hltpath:
        lumibyls=lumiCalcAPI.lumiForIds(session.nominalSchema(),irunlsdict,dataidmap,runsummaryMap=GrunsummaryData,beamstatusfilter=pbeammode,normmap=normvalueDict,lumitype='HF')
    else:
        referenceLabel='Recorded'
        hltname=options.hltpath
        hltpat=None
        if hltname=='*' or hltname=='all':
            hltname=None
        elif 1 in [c in hltname for c in '*?[]']: #is a fnmatch pattern
            hltpat=hltname
            hltname=None
        lumibyls=lumiCalcAPI.effectiveLumiForIds(session.nominalSchema(),irunlsdict,dataidmap,runsummaryMap=GrunsummaryData,beamstatusfilter=pbeammode,normmap=normvalueDict,hltpathname=hltname,hltpathpattern=hltpat)
    session.transaction().commit()
    rawdata={}
    #
    # start to plot
    #
    if options.action=='run':
        for run in sorted(lumibyls):
            rundata=lumibyls[run]
            if not options.hltpath:
                if len(rundata)!=0:
                    rawdata.setdefault('Delivered',[]).append((run,sum([t[5] for t in rundata if t[5]])))
                    rawdata.setdefault('Recorded',[]).append((run,sum([t[6] for t in rundata if t[6]])))
            else:
                labels=['Recorded']
                if len(rundata)!=0:
                    pathdict={}#{pathname:[eff,]}
                    rawdata.setdefault('Recorded',[]).append((run,sum([t[6] for t in rundata if t[6]])))
                    for perlsdata in rundata:
                        effdict=perlsdata[8]
                        pathnames=effdict.keys()
                        for thispath in pathnames:
                            pathdict.setdefault(thispath,[]).append(effdict[thispath][3])
                    for thispath in pathdict.keys():
                        labels.append(thispath)
                        rawdata.setdefault(thispath,[]).append((run,sum([t for t in pathdict[thispath] if t])))
        if options.yscale=='linear':
            m.plotSumX_Run(rawdata,resultlines,textoutput=outtextfilename,yscale='linear',referenceLabel=referenceLabel,labels=labels)
        elif options.yscale=='log':
            m.plotSumX_Run(rawdata,resultlines,textoutput=outtextfilename,yscale='log',referenceLabel=referenceLabel,labels=labels)
        else:
            m.plotSumX_Run(rawdata,resultlines,textoutput=outtextfilename,yscale='linear',referenceLabel=referenceLabel,labels=labels)
            m.plotSumX_Run(rawdata,resultlines,textoutput=outtextfilename,yscale='log',referenceLabel=referenceLabel,labels=labels)
    if options.action=='fill':
        for fill in sorted(fillrunMap):
            for run in fillrunMap[fill]:
                if not lumibyls.has_key(run): continue
                rundata=lumibyls[run]
                if not options.hltpath:
                    if len(rundata)!=0:
                        rawdata.setdefault('Delivered',[]).append((fill,run,sum([t[5] for t in rundata if t[5]])))
                        rawdata.setdefault('Recorded',[]).append((fill,run,sum([t[6] for t in rundata if t[6]])))
                else:
                    labels=['Recorded']
                    if len(rundata)!=0:
                        pathdict={}#{pathname:[eff,]}
                        rawdata.setdefault('Recorded',[]).append((fill,run,sum([t[6] for t in rundata if t[6]])))
                        for perlsdata in rundata:
                            effdict=perlsdata[8]
                            pathnames=effdict.keys()
                            for thispath in pathnames:
                                pathdict.setdefault(thispath,[]).append(effdict[thispath][3])
                        for thispath in pathdict.keys():
                            labels.append(thispath)
                            rawdata.setdefault(thispath,[]).append((fill,run,sum([t for t in pathdict[thispath] if t])))
        if options.yscale=='linear':
            m.plotSumX_Fill(rawdata,resultlines,textoutput=outtextfilename,yscale='linear',referenceLabel=referenceLabel)
        elif options.yscale=='log':
            m.plotSumX_Fill(rawdata,resultlines,textoutput=outtextfilename,yscale='log',referenceLabel=referenceLabel)
        else:
            m.plotSumX_Fill(rawdata,resultlines,textoutput=outtextfilename,yscale='linear',referenceLabel=referenceLabel)
            m.plotSumX_Fill(rawdata,resultlines,textoutput=outtextfilename,yscale='log',referenceLabel=referenceLabel)
    if options.action=='time':
        for run in sorted(lumibyls):
            rundata=lumibyls[run]
            if not options.hltpath:
                if len(rundata)!=0:
                    rawdata.setdefault('Delivered',[]).append((run,rundata[0][2],rundata[-1][2],sum([t[5] for t in rundata if t[5]])))
                    rawdata.setdefault('Recorded',[]).append((run,rundata[0][2],rundata[-1][2],sum([t[6] for t in rundata if t[6]])))
            else:
                labels=['Recorded']
                if len(rundata)!=0:
                    pathdict={}#{pathname:[eff,]}
                    rawdata.setdefault('Recorded',[]).append((run,rundata[0][2],rundata[-1][2],sum([t[6] for t in rundata if t[6]])))
                    for perlsdata in rundata:
                        effdict=perlsdata[8]
                        pathnames=effdict.keys()
                        for thispath in pathnames:
                            pathdict.setdefault(thispath,[]).append(effdict[thispath][3])
                    for thispath in pathdict.keys():
                        labels.append(thispath)
                        rawdata.setdefault(thispath,[]).append((run,rundata[0][2],rundata[-1][2],sum([t for t in pathdict[thispath] if t])))
        if options.yscale=='linear':
            m.plotSumX_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='linear',referenceLabel=referenceLabel)
        elif options.yscale=='log':
            mlog.plotSumX_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='log',referenceLabel=referenceLabel)
        else:
            m.plotSumX_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='linear',referenceLabel=referenceLabel)
            mlog.plotSumX_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='log',referenceLabel=referenceLabel)
    if options.action=='perday':
        daydict={}
        for run in sorted(lumibyls):
            rundata=lumibyls[run]
            for lsdata in rundata:
                lumilsnum=lsdata[0]
                lsTS=lsdata[2]
                dellum=lsdata[5]
                reclum=lsdata[6]
                daynumber=lsTS.date().toordinal()
                daydict.setdefault(daynumber,[]).append((run,lumilsnum,dellum,reclum))
        for day in sorted(daydict):
            daydata=daydict[day]
            daybeg=str(daydata[0][0])+':'+str(daydata[0][1])
            dayend=str(daydata[-1][0])+':'+str(daydata[-1][1])
            daydel=sum([t[2] for t in daydata if t[2]])
            dayrec=sum([t[3] for t in daydata if t[3]])
            rawdata.setdefault('Delivered',[]).append((day,daybeg,dayend,daydel))
            rawdata.setdefault('Recorded',[]).append((day,daybeg,dayend,dayrec))
        #print 'rawdata ',rawdata
        if options.yscale=='linear':
            m.plotPerdayX_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='linear',referenceLabel=referenceLabel)
        elif options.yscale=='log':
            mlog.plotPerdayX_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='log',referenceLabel=referenceLabel)
        else:
            m.plotPerdayX_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='linear',referenceLabel=referenceLabel)
            mlog.plotPerdayX_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='log',referenceLabel=referenceLabel)
    if options.action=='instpeakperday':
        daydict={}#{daynumber:[(runnumber,lumilsnum,inst),..]}
        for run in sorted(lumibyls):
            rundata=lumibyls[run]
            for lsdata in rundata:
                lumilsnum=lsdata[0]
                lsTS=lsdata[2]
                instlum=lsdata[5]/lslength
                daynumber=lsTS.date().toordinal()
                daydict.setdefault(daynumber,[]).append((run,lumilsnum,instlum))
        for day in sorted(daydict):
            daydata=daydict[day]
            daymax_val=0.0
            daymax_run=0
            daymax_ls=0
            for datatp in daydata:
                if datatp[2]>daymax_val:
                    daymax_val=datatp[2]
                    daymax_run=datatp[0]
                    daymax_ls=datatp[1]
            rawdata.setdefault('Delivered',[]).append((day,daymax_run,daymax_ls,daymax_val))
        if options.yscale=='linear':
            m.plotPeakPerday_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='linear',referenceLabel=referenceLabel)
        elif options.yscale=='log':
            mlog.plotPeakPerday_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='log',referenceLabel=referenceLabel)
        else:
            m.plotPeakPerday_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='linear',referenceLabel=referenceLabel)
            mlog.plotPeakPerday_Time(rawdata,resultlines,minTime=begtime,maxTime=endtime,textoutput=outtextfilename,yscale='log',referenceLabel=referenceLabel)
    if options.action=='inst':
        thisfillnumber=fillrunMap.keys()[0]
        starttime=0
        stoptime=0
        rawxdata=[]
        rawydata={}
        for run,rundata in lumibyls.items():
            rundata.sort()
            totlumils=0
            totcmsls=0
            starttime=rundata[0][2]
            stoptime=rundata[-1][2]
            for lsdata in rundata:
                lumilsnum=lsdata[0]
                totlumils+=1
                cmslsnum=lsdata[1]
                if cmslsnum!=0:
                    totcmsls+=1
                lsTS=lsdata[2]
                dellumi=lsdata[5]
                reclumi=lsdata[6]
                rawydata.setdefault('Delivered',[]).append(dellumi)
                rawydata.setdefault('Recorded',[]).append(reclumi)
            rawxdata=[run,thisfillnumber,starttime,stoptime,totlumils,totcmsls]
        m.plotInst_RunLS(rawxdata,rawydata,textoutput=None)
        
    if options.yscale=='linear':
        if options.interactive:
            m.drawInteractive()
            exit(0)
        else:
            if not options.withoutpng:
                m.drawPNG(outplotfilename+'.png')
                exit(0)
    elif options.yscale=='log':
        if options.interactive:
            mlog.drawInteractive()
            exit(0)
        else:
            if not options.withoutpng:
                mlog.drawPNG(outplotfilename+'_log.png')
                exit(0)
    else:
        if options.interactive:
            print 'cannot draw both log and linear from interactive'
            exit(0)
        if not options.withoutpng:
            m.drawPNG(outplotfilename+'.png')
            mlog.drawPNG(outplotfilename+'_log.png')
            exit(0)