L1RpcTriggerAnalysis/macros/makePressureGraph.py

#!/usr/bin/env python

import os, commands, sys
from array import array
import ROOT

if __name__ == "__main__":
    if len(sys.argv[1:]) != 1:
        print "\nUsage "+sys.argv[0]+" <file>"
        print "where <file> is a text file with 2 columns:"
        print "RUN_NUMBER  PRESSURE[mbar]\n"
        exit()

    ###Read data from txt file
    #file = open('pressure_twiki_160406-172268.dat')
    file = open(sys.argv[1])
    counter = 0
    runs = []
    pressure = []
    for line in file:
        print line.split()
        runs.append(float(line.split()[0]))
        if(len(line.split())>1):               
               pressure.append(float(line.split()[1]))
        else:
               pressure.append(0.0)               
               counter+=1

    x = array("d",runs)
    y = array("d",pressure)

    ###Define efficiency map histogram
    myFile = ROOT.TFile("PressureGraph.root","RECREATE")
    myGraph = ROOT.TGraph(len(runs),x,y)

    myGraph.Write()
    myFile.Write()
    myFile.Close()
Revision:	1.2
Committed:	Sat Sep 10 15:14:22 2011 UTC (13 years, 7 months ago) by cwiok
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	Mikolaj_cmssw52x, Kasia_30_04_2012, Mikolaj_08_10_2011, Mikolaj_05_10_2011, Mikolaj_14_09_2011
Changes since 1.1:	+27 -20 lines
Log Message:	Updated job submission, job merging and plotting scripts
#	User	Rev	Content
1	akalinow	1.1	#!/usr/bin/env python
2
3	cwiok	1.2	import os, commands, sys
4	akalinow	1.1	from array import array
5			import ROOT
6
7	cwiok	1.2	if __name__ == "__main__":
8			if len(sys.argv[1:]) != 1:
9			print "\nUsage "+sys.argv[0]+" <file>"
10			print "where <file> is a text file with 2 columns:"
11			print "RUN_NUMBER PRESSURE[mbar]\n"
12			exit()
13	akalinow	1.1
14	cwiok	1.2	###Read data from txt file
15			#file = open('pressure_twiki_160406-172268.dat')
16			file = open(sys.argv[1])
17			counter = 0
18			runs = []
19			pressure = []
20			for line in file:
21			print line.split()
22			runs.append(float(line.split()[0]))
23			if(len(line.split())>1):
24	akalinow	1.1	pressure.append(float(line.split()[1]))
25	cwiok	1.2	else:
26	akalinow	1.1	pressure.append(0.0)
27	cwiok	1.2	counter+=1
28	akalinow	1.1
29	cwiok	1.2	x = array("d",runs)
30			y = array("d",pressure)
31	akalinow	1.1
32	cwiok	1.2	###Define efficiency map histogram
33			myFile = ROOT.TFile("PressureGraph.root","RECREATE")
34			myGraph = ROOT.TGraph(len(runs),x,y)
35	akalinow	1.1
36	cwiok	1.2	myGraph.Write()
37			myFile.Write()
38			myFile.Close()