[Gate-users] Energy spectrum
Konstantinos P. Chatzipapas, MSc
konhat88 at yahoo.gr
Tue Sep 26 15:45:04 CEST 2017
When I had tried the physicslist way of physics I had some faults and couldn't work with it. So finally I used this way for physics:
/gate/physics/addProcess PhotoElectric/gate/physics/processes/PhotoElectric/setModel StandardModel/gate/physics/addProcess Compton
/gate/physics/processes/Compton/setModel PenelopeModel/gate/physics/addProcess RayleighScattering
/gate/physics/processes/RayleighScattering/setModel PenelopeModel/gate/physics/addProcess ElectronIonisation
/gate/physics/processes/ElectronIonisation/setModel StandardModel e-/gate/physics/addProcess Bremsstrahlung
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-/gate/physics/addProcess eMultipleScattering e-
/gate/physics/processList Enabled
/gate/physics/processList Initialized
Try this, with also the basic commands for the actor. You can always change the models later, when you make it work. ----------------------------------------------------- Konstantinos Chatzipapas
Medical Physicist, MSc
University of Patras
web page: Physics Magazine
Στις 4:15 μ.μ. Τρίτη, 26 Σεπτεμβρίου 2017, ο/η tony youness <tony_youness at live.com> έγραψε:
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# VERBOSE and VISUALISATION
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
/control/execute mac/verbose.mac
#/control/execute mac/visu.mac
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# A specific material database is created
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
/gate/geometry/setMaterialDatabase data/Materials.db
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# GEOMETRY
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# WORLD
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
/gate/world/setMaterial Air
/gate/world/geometry/setXLength 5.0 m
/gate/world/geometry/setYLength 5.0 m
/gate/world/geometry/setZLength 5.0 m
/gate/world/vis/setVisible 0
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# PATIENT/PHANTOM
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
/control/execute mac/Waterphantom_DSP95.mac
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# Accelerator GEOMETRY
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
/control/execute mac/Linac.mac
#=====================================================
# OUTPUT ACTORS
#=====================================================
/gate/actor/addActor EnergySpectrumActor spectrum
/gate/actor/spectrum/attachTo layer3
/gate/actor/spectrum/addFilter particleFilter
/gate/actor/spectrum/particleFilter/addParticle e-
/gate/actor/spectrum/save output/WaterSpec.root
/gate/actor/spectrum/energySpectrum/setEmin 0 eV
/gate/actor/spectrum/energySpectrum/setEmax 6 MeV
/gate/actor/spectrum/energySpectrum/setNumberOfBins 300
#=====================================================
# PHYSICS
#=====================================================
/gate/physics/addPhysicsList emlivermore
/gate/physics/Gamma/SetCutInRegion world 1 mm
/gate/physics/Electron/SetCutInRegion world 1 mm
/gate/physics/Positron/SetCutInRegion world 1 mm
/gate/physics/Gamma/SetCutInRegion patient 0.1 mm
/gate/physics/Electron/SetCutInRegion patient 0.1 mm
/gate/physics/Positron/SetCutInRegion patient 0.1 mm
/gate/physics/Gamma/SetCutInRegion mylinac 1 mm
/gate/physics/Electron/SetCutInRegion mylinac 1 mm
/gate/physics/Positron/SetCutInRegion mylinac 1 mm
/gate/physics/displayCuts
/gate/physics/print output/physics.txt
#=====================================================
# INITIALISATION
#=====================================================
/gate/run/initialize
#=====================================================
# INPUT BEAM
#=====================================================
/gate/source/addSource MyBeam phaseSpace
/gate/source/MyBeam/addPhaseSpaceFile data/IAEA.IAEAphsp
/gate/source/MyBeam/attachTo outPhSp
#=====================================================
# START BEAMS
#=====================================================
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto
/control/execute mac/timeslices.mac
/gate/application/setTotalNumberOfPrimaries 10000000
/gate/application/start
On 26/09/17 3:03 pm, Konstantinos P. Chatzipapas, MSc wrote:
If you can, please attach the mac file that you use. ----------------------------------------------------- Konstantinos Chatzipapas
Medical Physicist, MSc
University of Patras
web page: Physics Magazine
Στις 3:57 μ.μ. Τρίτη, 26 Σεπτεμβρίου 2017, ο/η tony youness<tony_youness at live.com> έγραψε:
Hi,I just tried to stick with only basic commands and still getting the same error.
On 26/09/17 2:35 pm, Konstantinos P. Chatzipapas, MSc wrote:
Hi,
Have you try to only use the basic commands that are given by the guide?? Do they work??
Energy spectrum
This actor builds four histograms: the initial kinetic energy of each track (energySpectrum), the energy deposition per event (edepHisto), the energy deposition per track (edepTrackHisto), the energy loss per track (eLossHisto). These histograms are stored in a root file. They takes into account the weight of particles. They are three commands to define the boundaries and the binning of the energy spectrum and three commands to define the boundaries and the binning of the energy loss histograms (edepHisto, edepTrackHisto, eLossHisto)./gate/actor/addActor EnergySpectrumActor MyActor
/gate/actor/MyActor/save MyOutputFile.root
/gate/actor/MyActor/energySpectrum/setEmin 0 eV
/gate/actor/MyActor/energySpectrum/setEmax 10 GeV
/gate/actor/MyActor/energySpectrum/setNumberOfBins 200
/gate/actor/MyActor/energyLossHisto/setEmin 0 eV
/gate/actor/MyActor/energyLossHisto/setEmax 15 MeV
/gate/actor/MyActor/energyLossHisto/setNumberOfBins 120
If you haven't try them, try. Because I see that you are trying with extra commands, such as AttachTo, Filter etc. Regards,----------------------------------------------------- Konstantinos Chatzipapas
Medical Physicist, MSc
University of Patras
web page: Physics Magazine
Στις 3:23 μ.μ. Τρίτη, 26 Σεπτεμβρίου 2017, ο/η tony youness<tony_youness at live.com> έγραψε:
Dear David,
-Gate v7.2
-Root v5.34/30
Tony
On 26/09/17 2:17 pm, David Boersma wrote:
> Dear Tony,
>
> For us to reproduce the error:
> - Which version of GATE?
> - Which version of ROOT?
>
> DavidB
>
> Den 26/09/2017 kl. 14:03, skrev tony youness:
>> Dear Maxime,
>>
>> Thank you for your response.
>>
>> I just tried to increase the Emin to 1 MeV and reduced the number of
>> bins to 100 (just for testing), and I still get the same error...
>>
>> Tony
>>
>>
>> On 26/09/17 1:51 pm, Maxime Chauvin wrote:
>>> Dear Tony,
>>>
>>> what I can tell is that the error seems to come from Root (TH1D).
>>> I don’t know what cause the error but maybe the number of bins is
>>> too high or the Emin is too small…
>>>
>>> Maxime
>>>
>>>> On 26 Sep 2017, at 13:17, tony youness <tony_youness at live.com> wrote:
>>>>
>>>> Dear fellow Gate users,
>>>>
>>>> I am trying to generate the electrons energy spectrum at 10 cm in
>>>> water
>>>> phantom. I used the Energy spectrum actor as follow:
>>>>
>>>> /gate/actor/addActor EnergySpectrumActor spectrum
>>>> /gate/actor/spectrum/attachTo layer3
>>>> /gate/actor/spectrum/addFilter particleFilter
>>>> /gate/actor/spectrum/particleFilter/addParticle e-
>>>> /gate/actor/spectrum/save output/WaterSpec.root
>>>> /gate/actor/spectrum/energySpectrum/setEmin 0 eV
>>>> /gate/actor/spectrum/energySpectrum/setEmax 6 MeV
>>>> /gate/actor/spectrum/energySpectrum/setNumberOfBins 300
>>>>
>>>> layer3 is a small water volume created in the water phantom at 10 cm.
>>>> When i launch the simulation i get the following error:
>>>>
>>>> Error in <TH1D::Rebin>: Illegal value of ngroup=10815379
>>>> Error in <TH1D::Rebin>: Illegal value of ngroup=10815379
>>>> Error in <TH2D::Rebin>: Illegal value of nxgroup=10815379
>>>> Error in <TH1D::Rebin>: Illegal value of ngroup=10815379
>>>> Error in <TH1D::Rebin>: Illegal value of ngroup=10815379
>>>> WARNING: Scene "none" not found.
>>>> /vis/scene/list to see scenes.
>>>>
>>>>
>>>>
>>>> ### === G4UAtomicDeexcitation::InitialiseForNewRun()
>>>> ### === PIXE model for hadrons: Empirical
>>>> ### === PIXE model for e+-: Livermore
>>>> Error in <TH2D::Fill>: Invalid signature - do nothing
>>>> Error in <TH2D::Fill>: Invalid signature - do nothing
>>>> Error in <TH2D::Fill>: Invalid signature - do nothing
>>>> Error in <TH2D::Fill>: Invalid signature - do nothing
>>>> Error in <TH2D::Fill>: Invalid signature - do nothing
>>>> Error in <TH2D::Fill>: Invalid signature - do nothing
>>>> Error in <TH2D::Fill>: Invalid signature - do nothing
>>>> ...
>>>>
>>>> Any ideas where this problem could come from?
>>>>
>>>> Regards,
>>>>
>>>> Tony Younes
>>>>
>>>> _______________________________________________
>>>> Gate-users mailing list
>>>> Gate-users at lists.opengatecollaboration.org
>>>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>>
>> _______________________________________________
>> Gate-users mailing list
>> Gate-users at lists.opengatecollaboration.org
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>>
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20170926/266788d0/attachment-0001.html>
More information about the Gate-users
mailing list