PDRates.py

Go to the documentation of this file.

#! /usr/bin/env python
import sys, imp
import os
import commands

from optparse import OptionParser


# --
# -- Usage :
# -- Rate within a given PD :
# --       ./PDRates.py  --runNumber 135525 --PD Mu
# -- Plot the rate of a PD versus the LS number:
# --       ./PDRates.py  --runNumber 135525 --PD Mu --perLS
# -- Rates in all PDs:
# --       ./PDRates.py  --runNumber 135525
# --


parser = OptionParser(usage="Example: ./PDRates.py --runNumber 146644 --PD Mu --perLS")
parser.add_option("--runNumber",  dest="RunNumber",  help="run number", type="int", metavar="RunNumber")
parser.add_option("--PD",dest="PrimaryDataset",help="PD name",metavar="PrimaryDataset")
parser.add_option("--perLS",action="store_true",default=False,dest="perLS",help="plot the rate vs the LS number",metavar="perLS")
parser.add_option("--logy",action="store_true",default=False,dest="logy",help="log scale for the y axis",metavar="logy")
parser.add_option("--saveplot",action="store_true",default=False,dest="saveplot",help="save the plot as a tmp.eps",metavar="saveplot")
(options, args) = parser.parse_args()


def PrimaryDatasets(Run):
    # -- List of Primary Datasets:

    dbs_cmd = """ dbs search --query='find primds.name 
            where run=%d and dataset like */RAW' """ % (Run)
    rows = commands.getoutput(dbs_cmd)
    lines = rows.split("\n")
    j=0
    print "\nThe primary datasets for this run are: \n"
    for Line in lines:
        j=j+1
        if j <=4:
            continue
        print Line
        line=Line.split()
        Datasets.append(line[0])
    print " "



def RateInPD(Run,PrimaryDataset,lsMin,lsMax,printLS=False):
    dbs_cmd = """ dbs search --query='find file,lumi,file.numevents, file.size
        where run=%d and dataset like /%s/*/RAW 
        and lumi >= %d and lumi <= %d
        and file.status=VALID '""" % (Run, PrimaryDataset,lsMin, lsMax)
    rows = commands.getoutput(dbs_cmd)
    lines = rows.split("\n")
    j=0
    LumiSections = []
    Files = []
    Evts = 0
    Size = 0
    LSinFile = {}
    NumberOfLSInFile = {}
    for Line in lines:
        j=j+1
        if j <=4:
            continue
        line=Line.split()
        LS = line[1]
        file = line[0]
        Nevts = int(line[2])
        size = int(line[3])
        LSinFile[LS] = file
        if LumiSections.count(LS) == 0:
            LumiSections.append(LS)
        if Files.count(file) == 0:
            Files.append(file)
            Evts += Nevts
            Size += size
            NumberOfLSInFile[file] =1
        else:
            NumberOfLSInFile[file] = NumberOfLSInFile[file]+1
        RatePerLS[LS] = Nevts
            
    Number_of_LS = len(LumiSections)
    LS_Length = 23.3
    if Run < 125100:
        LS_Length = 93.3
    rate = Evts / (Number_of_LS * LS_Length)
    if Evts > 0:
        size_per_event = (Size / Evts) / 1000.
    else:
        size_per_event=-1
    print "Rate in \t",PrimaryDataset,"\t is : \t",rate," Hz", " \t size is : \t",size_per_event, "kB / event "

    lsmin=9999999
    lsmax=-1
    for (LS,file) in LSinFile.iteritems():
        nls = NumberOfLSInFile[file]
        RatePerLS[LS] = RatePerLS[LS] / nls
        RatePerLS[LS] = RatePerLS[LS] / LS_Length
        if int(LS) > lsmax:
            lsmax=int(LS)
        if int(LS) < lsmin:
            lsmin=int(LS)
    if printLS:
        print "lsmin lsmax",lsmin,lsmax
        for ls in range(lsmin,lsmax):
            if not repr(ls) in RatePerLS.keys():
                RatePerLS[LS] = 0
                print "Missing LS ",ls


if __name__ == "__main__":

  if not options.RunNumber:
    print "wrong usage"
    exit(2)

  lsMin = -1
  lsMax = 999999

# -- does not work yet.
#  if options.lsMin:
#   lsMin = options.lsMin
#  if options.lsMax:
#   lsMax = options.lsMax


  print "\nRun Number: ",options.RunNumber
  if options.PrimaryDataset:
    Run = options.RunNumber
    PrimaryDataset = options.PrimaryDataset
    RatePerLS = {}
    if not options.perLS:
        RateInPD(Run,PrimaryDataset,lsMin, lsMax, False)
    if options.perLS:
        average = 0
        nLS_within_range = 0
        RateInPD(Run,PrimaryDataset,lsMin, lsMax, True)
        RatesTmp = open("rates_tmp.txt","w")
        #RatesTmpSort = open("rates_tmp_sort.txt","w")
        for (LS, rate) in RatePerLS.iteritems():
            RatesTmp.write(LS+"\t"+repr(rate)+"\n")
            #if int(LS) >=  lsMin and int(LS) <= lsMax:
                #nLS_within_range =nLS_within_range +1
                #average = average + rate
        #print "Average rate within ",options.lsMin," and ",options.lsMax," is: ",average/nLS_within_range
        #if os.path.exists("./rates_tmp_sort.txt"):
            #os.system("rm rates_tmp_sort.txt")
        #os.system("sort -n rates_tmp.txt > rates_tmp_sort.txt")
        RatesTmp.close()
        
        TempFile = open("tmp.gnuplot.txt","w")
        if options.logy:
            TempFile.write("set logscale y \n")
        if options.saveplot:
            TempFile.write("set term postscript eps enhanced \n")
            TempFile.write("set output \"tmp.eps\" \n")
        st_title = " \"Rates in PrimaryDataset " + PrimaryDataset+ " in Run " + repr(Run)+ "\" "
        TempFile.write("set title " + st_title + "\n")
                TempFile.write(" set pointsize 2. \n")
                TempFile.write(" set nokey \n")
        TempFile.write(" set xlabel \"LS number\" \n")
                TempFile.write(" set ylabel \"Rate (Hz)\" \n")
        TempFile.write(" plot \"rates_tmp.txt\" using 1:2 title 'Rate per LS' \n")
        if not options.saveplot:
                    TempFile.write(" pause -1")
        else:
            print "The plot is saved under tmp.eps"
                TempFile.close()

        os.system("gnuplot tmp.gnuplot.txt")


  else:
    Run = options.RunNumber
    Datasets = []
    PrimaryDatasets(Run)
    for PrimaryDataset in Datasets:
        RatePerLS = {}
        RateInPD(Run,PrimaryDataset,lsMin,lsMax)