1 | # $Id: readBremLPM.py,v 1.2 2008/08/22 08:16:17 schaelic Exp $ |
---|
2 | # GEANT4 tag $Name: geant4-09-04-ref-00 $ |
---|
3 | # |
---|
4 | # works together with BremLPMTest.cc |
---|
5 | # |
---|
6 | from ROOT import gROOT, gApplication, TFile, TCanvas |
---|
7 | |
---|
8 | # set nice style |
---|
9 | gROOT.SetStyle("Plain") |
---|
10 | |
---|
11 | # create 6 pads |
---|
12 | c1 = TCanvas('c1','CrossSection') |
---|
13 | c1.Divide(3,2) |
---|
14 | |
---|
15 | c2 = TCanvas('c2','DEDX') |
---|
16 | c2.Divide(3,2) |
---|
17 | |
---|
18 | #open root file |
---|
19 | f=TFile('eBremRel01.root') |
---|
20 | t=f.Get('info') |
---|
21 | zlink=t.GetLeaf('Z') |
---|
22 | t.Scan() |
---|
23 | |
---|
24 | #plot cross section histograms |
---|
25 | for i in range(t.GetEntries()): |
---|
26 | t.GetEntry(i) |
---|
27 | Z=int(zlink.GetValue()) |
---|
28 | c1.cd(i+1) |
---|
29 | c1.GetPad(i+1).SetLogx() |
---|
30 | |
---|
31 | gR=f.Get('modelR;'+str(i+1)) |
---|
32 | gR.GetHistogram().SetTitle('Z='+str(Z)) |
---|
33 | gR.Draw('Alp'); |
---|
34 | minR=gR.GetHistogram().GetMinimum() |
---|
35 | maxR=gR.GetHistogram().GetMaximum() |
---|
36 | |
---|
37 | g1=f.Get('model1;'+str(i+1)) |
---|
38 | g1.Draw('lp'); |
---|
39 | min1=g1.GetHistogram().GetMinimum() |
---|
40 | max1=g1.GetHistogram().GetMaximum() |
---|
41 | |
---|
42 | gR.SetMaximum(1.1*max(max1,maxR)) |
---|
43 | gR.SetMinimum(0.9*min(min1,minR)) |
---|
44 | print Z |
---|
45 | |
---|
46 | #plot dEdx histograms |
---|
47 | for i in range(t.GetEntries()): |
---|
48 | t.GetEntry(i) |
---|
49 | Z=int(zlink.GetValue()) |
---|
50 | c2.cd(i+1) |
---|
51 | c2.GetPad(i+1).SetLogx() |
---|
52 | c2.GetPad(i+1).SetLogy() |
---|
53 | |
---|
54 | gR=f.Get('xDEDXR;'+str(i+1)) |
---|
55 | gR.GetHistogram().SetTitle('Z='+str(Z)) |
---|
56 | gR.Draw('Alp'); |
---|
57 | minR=gR.GetHistogram().GetMinimum() |
---|
58 | maxR=gR.GetHistogram().GetMaximum() |
---|
59 | |
---|
60 | g1=f.Get('xDEDX1;'+str(i+1)) |
---|
61 | g1.Draw('lp'); |
---|
62 | min1=g1.GetHistogram().GetMinimum() |
---|
63 | max1=g1.GetHistogram().GetMaximum() |
---|
64 | |
---|
65 | # gR.SetMaximum(1.e-3) |
---|
66 | gR.SetMinimum(1.e-6*max(max1,maxR)) |
---|
67 | print Z |
---|
68 | |
---|
69 | #done |
---|
70 | gApplication.Run() |
---|