L1RpcTriggerAnalysis/macros/makePressureGraph.py

#!/usr/bin/env python

import os, commands
from array import array
import ROOT


###Read data from txt file
file = open('pressure_twiki_160406-171484.dat')
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.1
Committed:	Tue Jul 26 20:00:50 2011 UTC (13 years, 9 months ago) by akalinow
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	Artur_27_07_2011
Log Message:	* empty log message *
#	User	Rev	Content
1	akalinow	1.1	#!/usr/bin/env python
2
3			import os, commands
4			from array import array
5			import ROOT
6
7
8			###Read data from txt file
9			file = open('pressure_twiki_160406-171484.dat')
10			counter = 0
11			runs = []
12			pressure = []
13			for line in file:
14			print line.split()
15			runs.append(float(line.split()[0]))
16			if(len(line.split())>1):
17			pressure.append(float(line.split()[1]))
18			else:
19			pressure.append(0.0)
20			counter+=1
21
22			x = array("d",runs)
23			y = array("d",pressure)
24
25			###Define efficiency map histogram
26			myFile = ROOT.TFile("PressureGraph.root","RECREATE")
27			myGraph = ROOT.TGraph(len(runs),x,y)
28
29			myGraph.Write()
30			myFile.Write()
31			myFile.Close()