[Gate-users] Positron range measurement

nini latiff katoshka_7880 at yahoo.com
Thu Sep 1 02:35:28 CEST 2016


Hi everyone, I am Rukiah. I am new with this GATE. Currently I working to measure positron range. I already put my tracklength actor so that I can obtain positron range output. When I run the macrofile below, visualization is good. However, I can not open the MyOutputFile.root as the message says 'failed to read the file type data'. Please help regarding this matter as I have no background in programming. Thank you.

#=====================================================# VISUALISATION#=====================================================

#/vis/disable
/vis/open                           OGLS/vis/viewer/set/viewpointThetaPhi   60 60/vis/viewer/zoom                    1.0/vis/drawVolume/vis/scene/add/trajectories 

/control/execute                  Verbose.mac

#=====================================================# GEOMETRY#=====================================================

/gate/geometry/setMaterialDatabase    GateMaterials.db

##     W O R L D#/gate/world/geometry/setXLength       40. cm/gate/world/geometry/setYLength       40. cm/gate/world/geometry/setZLength       40. cm#-------------------oooooOOOOO00000OOOOOooooo---------------------##                                                                 ##     D E F I N I T I O N   A N D   D E S C R I T I O N           ##             O F   Y O U R   P H A N T O M                       ##                                                                 ##-------------------oooooOOOOO00000OOOOOooooo---------------------#

/gate/world/daughters/name phantom/gate/world/daughters/insert cylinder/gate/phantom/setMaterial Water/gate/phantom/vis/forceWireframe/gate/phantom/vis/setColor red/gate/phantom/geometry/setRmax   5 cm/gate/phantom/geometry/setHeight 10 cm/gate/phantom/attachPhantomSD
#=====================================================# PHYSICS#=====================================================#/gate/physics/addProcess PhotoElectric#/gate/physics/processes/PhotoElectric/setModel StandardModells
#/gate/physics/addProcess Compton#/gate/physics/processes/Compton/setModel StandardModel
#/gate/physics/addProcess RayleighScattering#/gate/physics/processes/RayleighScattering/setModel PenelopeModel
#/gate/physics/addProcess ElectronIonisation#/gate/physics/processes/ElectronIonisation/setModel StandardModel e-#/gate/physics/processes/ElectronIonisation/setModel StandardModel e+
#/gate/physics/addProcess Bremsstrahlung#/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-#/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+
/gate/physics/addProcess PositronAnnihilation
#/gate/physics/addProcess MultipleScattering e+#/gate/physics/addProcess MultipleScattering e-
/gate/physics/processList Enabled/gate/physics/processList Initialized
/gate/run/initialize#=====================================================# CUTS#=====================================================
#/gate/physics/Gamma/SetCutInRegion      crystal 1.0 cm#/gate/physics/Electron/SetCutInRegion   crystal 1.0 cm#/gate/physics/Positron/SetCutInRegion   crystal 1.0 cm

#/gate/physics/Gamma/SetCutInRegion      phantom 0.1 mm#/gate/physics/Electron/SetCutInRegion   phantom 0.1 mm/gate/physics/Positron/SetCutInRegion   phantom 10 cm
/gate/physics/SetMaxStepSizeInRegion    phantom 1 cm
#TrackLengthActor
/gate/actor/addActor  TrackLengthActor  MyActor/gate/actor/MyActor/save  MyOutputFile.root/gate/actor/MyActor/setLmin 0 mm/gate/actor/MyActor/setLmax 10 cm/gate/actor/MyActor/setNumberOfBins 200
#-------------------oooooOOOOO00000OOOOOooooo---------------------##                                                                 ##                   D E F I N I T I O N   O F                     ##            Y O U R   V E R B O S I T Y   L E V E L              ##                                                                 ##-------------------oooooOOOOO00000OOOOOooooo---------------------#
/control/execute                    Verbose.mac
#=====================================================# INITIALISATION#=====================================================
/gate/run/initialize
#-------------------oooooOOOOO00000OOOOOooooo---------------------##                                                                 ##                   D E F I N I T I O N   O F                     ##                Y O U R   O U T P U T   F I L E                  ##                                                                 ##-------------------oooooOOOOO00000OOOOOooooo---------------------#

#  R A N D O M
# JamesRandom Ranlux64 MersenneTwister/gate/random/setEngineName JamesRandom#/gate/random/setEngineSeed default#/gate/random/setEngineSeed auto#/gate/random/setEngineSeed 123456789/gate/random/setEngineSeed default#/gate/random/resetEngineFrom fileName/gate/random/verbose 1

#       S E T U P  - R O O T  F I L E/gate/output/root/enable/gate/output/root/setFileName MySimu/gate/output/root/setRootSinglesAdderFlag 1/gate/output/root/setRootSinglesReadoutFlag 1/gate/output/root/setRootHitFlag 1/gate/output/root/setRootSinglesFlag 1/gate/output/root/setRootCoincidencesFlag 1
#-------------------oooooOOOOO00000OOOOOooooo---------------------##                                                                 ##         P O S I T R O N   E M I S S I O N   S O U R C E         ##                                                                 ##                Ga-68              ##                                                                 #         #  NB :  In this example, each isotope source is generated        ##    in a sphere placed at the center of the FOV.           ##                The sphere radius is set at 0.5 cm.              ##                                                                 #   #-------------------oooooOOOOO00000OOOOOooooo---------------------#
## For Ga68#
/gate/source/addSource sourceGa68/gate/source/sourceGa68/setActivity 10. becquerel/gate/source/sourceGa68/gps/particle e+/gate/source/sourceGa68/setForcedUnstableFlag true/gate/source/sourceGa68/gps/energytype Ga68/gate/source/sourceGa68/setForcedHalfLife 6586 s/gate/source/sourceGa68/gps/type Volume/gate/source/sourceGa68/gps/shape Sphere/gate/source/sourceGa68/gps/radius 0.01 mm/gate/source/sourceGa68/gps/centre 0. 0. 0. cm/gate/source/sourceGa68/vis/setColor red/gate/source/sourceGa68/gps/confine NULL /gate/source/sourceGa68/gps/angtype iso


/gate/source/list#-------------------oooooOOOOO00000OOOOOooooo---------------------##                                                                 ##               S T A R T   A C Q U I S I T I O N                 ##                 T I M E   P A R A M E T E R S                   ##                                                                 ##-------------------oooooOOOOO00000OOOOOooooo---------------------#
## EXPERIMENT#/gate/application/setTimeSlice     120.  s/gate/application/setTimeStart     0.   s/gate/application/setTimeStop      240.  s
## LET'S RUN THE SIMULATION!#/gate/application/startDAQ#/gate/output/allowNoOutput


RukiahPhD CandidateThe University of SydneyAustralia 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20160901/d9ad2176/attachment-0001.html>


More information about the Gate-users mailing list