PythonAnalysis/test/trackerHits.py

# first load cmstools and ROOT classes
from cmstools import *
from ROOT import *

# opening file and accessing branches
print "Opening SimHit file"
file = TFile("simevent.root")
events = file.Get("Events")
branch = events.GetBranch("PSimHit_r_TrackerHitsTIBLowTof.obj")

simHit = std.vector(PSimHit)()
branch.SetAddress(simHit)

histo = TH1F("tofhits", "Tof of hits", 100, -0.5, 50)

# loop over all events 
for ev in all(events):
    branch.GetEntry(ev)
    for hit in all(simHit):
        histo.Fill(hit.timeOfFlight())

hFile = TFile("histo.root", "RECREATE")
histo.Write()

gROOT.SetBatch()
gROOT.SetStyle("Plain")

c = TCanvas()
histo.Draw()
c.SaveAs("tofhits.jpg")
Revision:	1.2
Committed:	Wed Mar 22 13:49:03 2006 UTC (19 years, 1 month ago) by hegner
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	CMSSW_1_3_6, CMSSW_1_3_1_HLT6, CMSSW_1_3_5, CMSSW_1_3_1_HLT5, CMSSW_1_3_1_HLT4, CMSSW_1_3_1_HLT3, CMSSW_1_3_4, CMSSW_1_3_3, CMSSW_1_3_2, CMSSW_1_2_6, CMSSW_1_2_5, CMSSW_1_3_1, CMSSW_1_2_4, CMSSW_1_3_0, CMSSW_1_3_0_pre7, CMSSW_1_3_0_pre6, CMSSW_1_2_0_4821, CMSSW_1_2_0_g4_82p01, CMSSW_1_3_0_pre4, CMSSW_1_2_3, CMSSW_1_3_0_pre3, CMSSW_1_2_0_g4_82, CMSSW_1_3_0_SLC4_pre2, CMSSW_1_3_0_pre2, CMSSW_1_2_2, CMSSW_1_2_1, CMSSW_1_3_0_SLC4_pre1, CMSSW_1_3_0_pre1, CMSSW_1_2_0_g4_81, CMSSW_1_2_0_SL4, CMSSW_1_2_0, CMSSW_1_2_0_pre9, CMSSW_1_2_0_pre8_g4_81, CMSSW_1_2_0_pre8, CMSSW_1_1_2, CMSSW_1_2_0_pre7, CMSSW_1_2_0_pre6, CMSSW_1_2_0_pre5, CMSSW_1_2_0_pre4, CMSSW_1_1_1, CMSSW_1_2_0_pre3, CMSSW_1_0_6, CMSSW_1_2_0_pre2, CMSSW_1_0_5, CMSSW_1_1_0, CMSSW_1_2_0_pre1, CMSSW_1_0_4, CMSSW_1_1_0_pre4, CMSSW_1_0_3, CMSSW_1_0_0_g4_81, CMSSW_1_1_0_pre3, CMSSW_1_0_2, CMSSW_1_0_1, CMSSW_1_1_0_pre2, CMSSW_1_0_0, CMSSW_1_0_0_pre5, CMSSW_1_1_0_pre1, CMSSW_0_8_4, CMSSW_1_0_0_pre4, CMSSW_1_0_0_pre3, CMSSW_0_9_2, CMSSW_0_8_3, CMSSW_1_0_0_pre2, CMSSW_0_9_0_pre2_g4_81, CMSSW_1_0_0_pre1, CMSSW_0_9_1, CMSSW_0_8_2, CMSSW_0_9_0, CMSSW_0_9_0_pre3, CMSSW_0_9_0_pre2, CMSSW_0_9_0_pre1, CMSSW_0_8_1, CMSSW_0_8_0, CMSSW_0_8_0_pre5, CMSSW_0_8_0_pre4, CMSSW_0_8_0_pre3, CMSSW_0_7_2, CMSSW_0_8_0_pre2, CMSSW_0_7_1, CMSSW_0_8_0_pre1, CMSSW_0_7_0, CMSSW_0_7_0_pre6, CMSSW_0_7_0_pre5, CMSSW_0_7_0_pre4, CMSSW_0_7_0_pre3, CMSSW_0_6_1, CMSSW_0_7_0_pre2, V00-00-04, CMSSW_0_7_0_pre1, CMSSW_0_6_0, CMSSW_0_6_0_pre7, CMSSW_0_6_0_pre6, CMSSW_0_6_0_pre5, CMSSW_0_6_0_pre4, CMSSW_0_6_0_pre3, CMSSW_0_6_0_pre2, V00-00-03, V00-00-02
Changes since 1.1:	+9 -14 lines
Log Message:	change to new startup and tab completion for not loaded libs use of new iterator generator
#	User	Rev	Content
1	hegner	1.2	# first load cmstools and ROOT classes
2			from cmstools import *
3	hegner	1.1	from ROOT import *
4
5			# opening file and accessing branches
6			print "Opening SimHit file"
7			file = TFile("simevent.root")
8			events = file.Get("Events")
9			branch = events.GetBranch("PSimHit_r_TrackerHitsTIBLowTof.obj")
10	hegner	1.2
11	hegner	1.1	simHit = std.vector(PSimHit)()
12			branch.SetAddress(simHit)
13
14	hegner	1.2	histo = TH1F("tofhits", "Tof of hits", 100, -0.5, 50)
15	hegner	1.1
16	hegner	1.2	# loop over all events
17			for ev in all(events):
18	hegner	1.1	branch.GetEntry(ev)
19	hegner	1.2	for hit in all(simHit):
20			histo.Fill(hit.timeOfFlight())
21	hegner	1.1
22			hFile = TFile("histo.root", "RECREATE")
23	hegner	1.2	histo.Write()
24	hegner	1.1
25			gROOT.SetBatch()
26			gROOT.SetStyle("Plain")
27
28			c = TCanvas()
29			histo.Draw()
30			c.SaveAs("tofhits.jpg")