emcalc_gui Namespace Reference

Classes
class	MainWindow
class	TextView

Functions
def	g4_configure ()
def	plot_charged (material, pname)
def	plot_gamma (material)
def	main ()

Function Documentation

g4_configure()

def emcalc_gui.g4_configure

(

)

Definition at line 20 of file emcalc_gui.py.

def g4_configure() :
  EmPlot.Configure()
  g4py.ExN03pl.Construct()

# -------------------------------------------------------------------
# plot for chaged particles
# -------------------------------------------------------------------

Here is the call graph for this function:

Here is the caller graph for this function:

main()

def emcalc_gui.main

(

)

Definition at line 342 of file emcalc_gui.py.

def main() :
  SetG4PyCoutDestination()

  default_stdout = sys.stdout
  global mycout
  sys.stdout = mycout = StringIO()

  # G4 setup
  g4_configure()

  # start GUI
  application = MainWindow()
  gtk.main()

  gTerminate()
  sys.stdout = default_stdout

Here is the call graph for this function:

plot_charged()

def emcalc_gui.plot_charged	(		material,
			pname
	)

Definition at line 27 of file emcalc_gui.py.

def plot_charged(material, pname) :
  EmPlot.SetMaterial(material)

  # initialize G4 kernel
  gRunManager.Initialize()
  gRunManagerKernel.RunInitialization()

  # energy
  elist= []
  for n in range(-3, 3):
    for i in range(10,99):
      elist.append(i/10.*10.**n *MeV)

  # calculate stopping power
  global mycout
  mycout.close()
  sys.stdout = mycout = StringIO()
  dedx_list= g4py.emcalculator.CalculateDEDX(pname, material, elist, 1)
  xlist_tot=[]
  xlist_ioni=[]
  xlist_brems=[]

  for x in dedx_list:
    xlist_tot.append((x[0], x[1]["tot"]/(MeV*cm2/g)))
    xlist_ioni.append((x[0], x[1]["ioni"]/(MeV*cm2/g)))
    xlist_brems.append((x[0], x[1]["brems"]/(MeV*cm2/g)))

  # make plot
  global myCanvas, aplot, bplot, cplot
  myCanvas = EmPlot.init_root()
  aplot = EmPlot.make_plot(xlist_tot, pname+" Stopping Power ("+material+")",
                           "dE/dX (MeV cm^{2}/g)")
  bplot = EmPlot.make_plot(xlist_ioni, "Stopping Power ("+material+")",
                           "dE/dX (MeV cm^{2}/g)", 1)
  cplot = EmPlot.make_plot(xlist_brems, "Stopping Power ("+material+")",
                           "dE/dX (MeV cm^{2}/g)", 3)
  myCanvas.SaveAs("/tmp/sp.png")


# -------------------------------------------------------------------
# plot for gamma
# -------------------------------------------------------------------

Here is the call graph for this function:

Here is the caller graph for this function:

plot_gamma()

def emcalc_gui.plot_gamma

(

material

)

Definition at line 69 of file emcalc_gui.py.

def plot_gamma(material) :
  EmPlot.SetMaterial(material)

  # initialize G4 kernel
  gRunManager.Initialize()
  gRunManagerKernel.RunInitialization()

  # energy
  elist= []
  for n in range(-3, 4):
    for i in range(10,99):
      elist.append(i/10.*10.**n *MeV)

  # calculate cross sections
  global mycout
  mycout.close()
  sys.stdout = mycout = StringIO()
  xsection_list= g4py.emcalculator.CalculatePhotonCrossSection(material,
                                                               elist, 1)
  xlist_tot=[]
  xlist_comp=[]
  xlist_pe=[]
  xlist_conv=[]
  for x in xsection_list:
    xlist_tot.append((x[0]/MeV, x[1]["tot"]/(cm2/g)))
    xlist_comp.append((x[0]/MeV, x[1]["compt"]/(cm2/g)))
    xlist_pe.append((x[0]/MeV, x[1]["phot"]/(cm2/g)))
    xlist_conv.append((x[0]/MeV, x[1]["conv"]/(cm2/g)))

  # make plots
  global myCanvas, aplot, bplot, cplot, dplot
  myCanvas = EmPlot.init_root()
  aplot = EmPlot.make_plot(xlist_tot,  "Photon Cross Section ("+material+")",
                             "Cross Section (cm^{2}/g)")
  bplo = EmPlot.make_plot(xlist_comp, "Photon Cross Section ("+material+")",
                             "Cross Section (cm^{2}/g)", 1)
  cplot = EmPlot.make_plot(xlist_pe,   "Photon Cross Section ("+material+")",
                             "Cross Section (cm^{2}/g)", 7)
  dplot = EmPlot.make_plot(xlist_conv, "Photon Cross Section ("+material+")",
                             "Cross Section (cm^{2}/g)", 3)
  myCanvas.SaveAs("/tmp/cs.png")


# ==================================================================
# GUI
# ==================================================================

Here is the call graph for this function:

Here is the caller graph for this function: