Esempi in Python per EmPlot

Esempio n. 1

File: emcalc_gui.py Progetto: C-D-D/Geant4.9.6

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

Esempio n. 2

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

Esempio n. 3

File: emcalc_gui.py Progetto: C-D-D/Geant4.9.6

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

Esempio n. 4

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

Esempio n. 5

File: pplot.py Progetto: cesarotti/geant4-1

# 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.0 * 10.0 ** n * MeV)


# calculate cross sections
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 a plot
myCanvas = EmPlot.init_root()
aplot = EmPlot.make_plot(xlist_tot, "Photon Cross Section (" + material + ")", "Cross Section (cm^{2}/g)")
bplot = 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)

Esempio n. 6

File: eplot.py Progetto: Geant4/geant4

material = "G4_Cu"
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.0 * 10.0 ** n * MeV)

# calculate stopping power
pname = "e-"
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
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)

Esempio n. 7

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

Esempio n. 8

File: eplot.py Progetto: zoulianmp/mt-g4py

#!/usr/bin/python
# ==================================================================
# An example of ploting by EmCalculator
#
# Plotting photon cross sections and stopping power
# ==================================================================
from Geant4 import *
import g4py.ExN03pl
import g4py.emcalculator
import EmPlot

# initialize
EmPlot.Configure()

# user physics list
g4py.ExN03pl.Construct()

# target material
material = "G4_Cu"
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)

Esempio n. 9

File: pplot.py Progetto: zoulianmp/mt-g4py

#!/usr/bin/python
# ==================================================================
# An example of ploting by EmCalculator
#
# Plotting photon cross sections and stopping power
# ==================================================================
from Geant4 import *
import g4py.ExN03pl
import g4py.emcalculator
import EmPlot

# initialize
EmPlot.Configure()

# user physics list
g4py.ExN03pl.Construct()

# target material
material= "G4_Pb"
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)

Esempio n. 10

File: emcalc_gui.py Progetto: zoulianmp/mt-g4py

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