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")