8import g4py.NISTmaterials
19 g4py.NISTmaterials.Construct()
20 g4py.ezgeom.Construct()
26 material= gNistManager.FindOrBuildMaterial(material_name)
27 g4py.ezgeom.SetWorldMaterial(material)
34from math
import log, log10, sqrt, ceil, floor
35from array
import array
51 ixminlog= xminlog-0.5-xmargin
53 return [10**ixminlog, 10**ixmaxlog]
62 ROOT.gStyle.SetTextFont(82)
64 ROOT.gStyle.SetTitleFont(82,
"X")
65 ROOT.gStyle.SetTitleFontSize(0.04)
66 ROOT.gStyle.SetLabelFont(82,
"X")
67 ROOT.gStyle.SetTitleFont(82,
"Y")
68 ROOT.gStyle.SetLabelFont(82,
"Y")
71 ROOT.gStyle.SetErrorX(0)
73 canvas= ROOT.TCanvas(
"g4plot",
"g4plot", 620, 30, 600, 600)
84def make_plot(xlist, user_title, axis_titile, q_super_impose=0):
86 ekin_array, y_array = array(
'd'), array(
'd')
89 ekin_array.append(x[0])
101 if(q_super_impose==0):
104 frame= ROOT.TH1F(
"dumy", htit, 1, xrange[0], xrange[1]);
105 frame.SetMinimum(yrange[0]);
106 frame.SetMaximum(yrange[1]);
107 frame.SetXTitle(
"Kinetic Energy (MeV)")
108 frame.GetXaxis().SetLabelSize(0.025)
109 frame.GetXaxis().SetTitleSize(0.03)
110 frame.SetYTitle(axis_titile)
111 frame.GetYaxis().SetLabelSize(0.025)
112 frame.GetYaxis().SetTitleSize(0.03)
116 plot= ROOT.TGraph(len(ekin_array), ekin_array, y_array)
118 plot.SetLineColor(q_super_impose+1)
def make_plot(xlist, user_title, axis_titile, q_super_impose=0)
def plot_range(xmin, xmax, xmargin=0.)
def SetMaterial(material_name)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
T min(const T t1, const T t2)
brief Return the smallest of the two arguments