74 gRunManager.Initialize()
75 gRunManagerKernel.RunInitialization()
79 for n
in range(-3, 4):
80 for i
in range(10,99):
81 elist.append(i/10.*10.**n *MeV)
86 sys.stdout = mycout = StringIO()
87 xsection_list= g4py.emcalculator.CalculatePhotonCrossSection(material,
93 for x
in xsection_list:
94 xlist_tot.append((x[0]/MeV, x[1][
"tot"]/(cm2/g)))
95 xlist_comp.append((x[0]/MeV, x[1][
"compt"]/(cm2/g)))
96 xlist_pe.append((x[0]/MeV, x[1][
"phot"]/(cm2/g)))
97 xlist_conv.append((x[0]/MeV, x[1][
"conv"]/(cm2/g)))
100 global myCanvas, aplot, bplot, cplot, dplot
103 "Cross Section (cm^{2}/g)")
105 "Cross Section (cm^{2}/g)", 1)
107 "Cross Section (cm^{2}/g)", 7)
109 "Cross Section (cm^{2}/g)", 3)
110 myCanvas.SaveAs(
"/tmp/cs.png")