[Gate-users] User defined energy histogram
venialgo at cae.cnea.gov.ar
venialgo at cae.cnea.gov.ar
Fri Aug 28 17:00:46 CEST 2009
Dear Nicholas Forwood,
Some time ago I've had the same problem but fortunately I've found some
information about GPS in internet.
Here is the link to GPS website, and there you can find all the
information that you need to create defined energy histogram.
http://reat.space.qinetiq.com/gps/new_gps_sum_files/gps_sum.htm
The problem with the GPS's online tutorial is that the commands are
slightly different. For instance, the command /gps/ang/type is explained
in the online tutorial, but in Gate you have to type
/gate/source/USER_SOURCE/gps/angtype iso
And here you have an example of Na-22 source, based on
http://ie.lbl.gov/toi/index.asp.
In that webpage there is a java applet which describes the pdf of the
possitron, and you can take points of that pdf to describe it in the macro
file.
I'm using one source to simulate the gamma emission, and another one to
simulate the positron emission.
# Na-22 punctual 7 uCi 99.94% for1274 keV gamma
/gate/source/addSource NaGamma
/gate/source/NaGamma/setActivity 259000 Bq
/gate/source/NaGamma/gps/centre 0. 0. 0. cm
/gate/source/NaGamma/gps/particle gamma
/gate/source/NaGamma/gps/energytype Mono
/gate/source/NaGamma/gps/monoenergy 1.274 MeV
/gate/source/NaGamma/gps/type Volume
/gate/source/NaGamma/gps/shape Sphere
/gate/source/NaGamma/gps/radius 1 mm
/gate/source/NaGamma/gps/angtype iso
/gate/source/NaGamma/gps/mintheta 0. deg
/gate/source/NaGamma/gps/maxtheta 180. deg
/gate/source/NaGamma/gps/minphi 0. deg
/gate/source/NaGamma/gps/maxphi 360. deg
/gate/source/list
# Na-22 punctual 7 uCi for possitron source.
# the pdf was token from: http://ie.lbl.gov/toi/
/gate/source/addSource NaPositron
/gate/source/NaPositron/setActivity 259000 Bq
/gate/source/NaPositron/gps/centre 0. 0. 0. cm
/gate/source/NaPositron/gps/particle e+
/gate/source/NaPositron/gps/energytype Arb
/gate/source/NaPositron/gps/histname arb
/gate/source/NaPositron/gps/emin 0.003 MeV
/gate/source/NaPositron/gps/emax 1.825 MeV
# I pick those pdf's points from the java applet.
/gate/source/NaPositron/gps/histpoint 0.003 73
/gate/source/NaPositron/gps/histpoint 0.024 942
/gate/source/NaPositron/gps/histpoint 0.055 1760
/gate/source/NaPositron/gps/histpoint 0.108 2560
/gate/source/NaPositron/gps/histpoint 0.134 2750
/gate/source/NaPositron/gps/histpoint 0.155 2830
/gate/source/NaPositron/gps/histpoint 0.176 2860
/gate/source/NaPositron/gps/histpoint 0.196 2840
/gate/source/NaPositron/gps/histpoint 0.228 2720
/gate/source/NaPositron/gps/histpoint 0.275 2390
/gate/source/NaPositron/gps/histpoint 0.315 2000
/gate/source/NaPositron/gps/histpoint 0.338 1750
/gate/source/NaPositron/gps/histpoint 0.400 1040
/gate/source/NaPositron/gps/histpoint 0.426 715
/gate/source/NaPositron/gps/histpoint 0.463 394
/gate/source/NaPositron/gps/histpoint 0.488 203
/gate/source/NaPositron/gps/histpoint 0.515 60
/gate/source/NaPositron/gps/histpoint 0.533 10
/gate/source/NaPositron/gps/histpoint 0.541 1.5
/gate/source/NaPositron/gps/histpoint 0.544 0.38
/gate/source/NaPositron/gps/histpoint 1.127 0.239
/gate/source/NaPositron/gps/histpoint 1.419 0.183
/gate/source/NaPositron/gps/histpoint 1.764 0.0008
/gate/source/NaPositron/gps/histpoint 1.817 0
/gate/source/NaPositron/gps/arbint Lin
/gate/source/NaPositron/gps/type Volume
/gate/source/NaPositron/gps/shape Sphere
/gate/source/NaPositron/gps/radius 1 mm
/gate/source/NaPositron/gps/angtype iso
/gate/source/NaPositron/gps/mintheta 0. deg
/gate/source/NaPositron/gps/maxtheta 180. deg
/gate/source/NaPositron/gps/minphi 0. deg
/gate/source/NaPositron/gps/maxphi 360. deg
/gate/source/list
> Hi everyone, I'm trying to compare the output of a Y90 ion source (which I
> mentioned in previous emails to the list) with the output of a beta source
> which
> has the same energy spectrum as Y90 (I looked it up on the internet). My
> problem
> is that I can't find a clear complete description of how to use the
> gps/histpoint
> command. I've entered what I believe to be the correct histogram macro but
> the
> program produces primary electrons with several hundred MeV energy. Where
> have I
> gone wrong?
>
> /gate/source/addSource Y90Source
> /gate/source/Y90Source/setActivity 10. becquerel
> /gate/source/Y90Source/gps/particle e-
> /gate/source/Y90source/gps/energytype User
> /gate/source/Y90Source/gps/histname energy
> /gate/source/Y90Source/gps/emin 0. keV
> /gate/source/Y90Source/gps/emax 2284. keV
> /gate/source/Y90Source/gps/histpoint 10. 0.51
> /gate/source/Y90Source/gps/histpoint 20. 0.53
> /gate/source/Y90Source/gps/histpoint 30. 0.55
> /gate/source/Y90Source/gps/histpoint 40. 0.56
> /gate/source/Y90Source/gps/histpoint 50. 0.58
> /gate/source/Y90Source/gps/histpoint 60. 0.59
> /gate/source/Y90Source/gps/histpoint 70. 0.61
> /gate/source/Y90Source/gps/histpoint 80. 0.62
> /gate/source/Y90Source/gps/histpoint 90. 0.63
> /gate/source/Y90Source/gps/histpoint 100. 0.65
> /gate/source/Y90Source/gps/histpoint 110. 0.66
> /gate/source/Y90Source/gps/histpoint 120. 0.67
> /gate/source/Y90Source/gps/histpoint 130. 0.68
> /gate/source/Y90Source/gps/histpoint 140. 0.69
> /gate/source/Y90Source/gps/histpoint 150. 0.71
> /gate/source/Y90Source/gps/histpoint 160. 0.72
> /gate/source/Y90Source/gps/histpoint 170. 0.73
> /gate/source/Y90Source/gps/histpoint 180. 0.74
> /gate/source/Y90Source/gps/histpoint 190. 0.75
> /gate/source/Y90Source/gps/histpoint 200. 0.76
> /gate/source/Y90Source/gps/histpoint 210. 0.77
> /gate/source/Y90Source/gps/histpoint 220. 0.78
> /gate/source/Y90Source/gps/histpoint 230. 0.78
> /gate/source/Y90Source/gps/histpoint 240. 0.79
> /gate/source/Y90Source/gps/histpoint 250. 0.8
> /gate/source/Y90Source/gps/histpoint 260. 0.81
> /gate/source/Y90Source/gps/histpoint 270. 0.82
> /gate/source/Y90Source/gps/histpoint 280. 0.82
> /gate/source/Y90Source/gps/histpoint 290. 0.83
> /gate/source/Y90Source/gps/histpoint 300. 0.84
> /gate/source/Y90Source/gps/histpoint 310. 0.85
> /gate/source/Y90Source/gps/histpoint 320. 0.85
> /gate/source/Y90Source/gps/histpoint 330. 0.86
> /gate/source/Y90Source/gps/histpoint 340. 0.87
> /gate/source/Y90Source/gps/histpoint 350. 0.87
> /gate/source/Y90Source/gps/histpoint 360. 0.88
> /gate/source/Y90Source/gps/histpoint 370. 0.88
> /gate/source/Y90Source/gps/histpoint 380. 0.89
> /gate/source/Y90Source/gps/histpoint 390. 0.89
> /gate/source/Y90Source/gps/histpoint 400. 0.9
> /gate/source/Y90Source/gps/histpoint 410. 0.9
> /gate/source/Y90Source/gps/histpoint 420. 0.91
> /gate/source/Y90Source/gps/histpoint 430. 0.91
> /gate/source/Y90Source/gps/histpoint 440. 0.92
> /gate/source/Y90Source/gps/histpoint 450. 0.92
> /gate/source/Y90Source/gps/histpoint 460. 0.93
> /gate/source/Y90Source/gps/histpoint 470. 0.93
> /gate/source/Y90Source/gps/histpoint 480. 0.94
> /gate/source/Y90Source/gps/histpoint 490. 0.94
> /gate/source/Y90Source/gps/histpoint 500. 0.94
> /gate/source/Y90Source/gps/histpoint 510. 0.95
> /gate/source/Y90Source/gps/histpoint 520. 0.95
> /gate/source/Y90Source/gps/histpoint 530. 0.95
> /gate/source/Y90Source/gps/histpoint 540. 0.96
> /gate/source/Y90Source/gps/histpoint 550. 0.96
> /gate/source/Y90Source/gps/histpoint 560. 0.96
> /gate/source/Y90Source/gps/histpoint 570. 0.97
> /gate/source/Y90Source/gps/histpoint 580. 0.97
> /gate/source/Y90Source/gps/histpoint 590. 0.97
> /gate/source/Y90Source/gps/histpoint 600. 0.97
> /gate/source/Y90Source/gps/histpoint 610. 0.98
> /gate/source/Y90Source/gps/histpoint 620. 0.98
> /gate/source/Y90Source/gps/histpoint 630. 0.98
> /gate/source/Y90Source/gps/histpoint 640. 0.98
> /gate/source/Y90Source/gps/histpoint 650. 0.98
> /gate/source/Y90Source/gps/histpoint 660. 0.99
> /gate/source/Y90Source/gps/histpoint 670. 0.99
> /gate/source/Y90Source/gps/histpoint 680. 0.99
> /gate/source/Y90Source/gps/histpoint 690. 0.99
> /gate/source/Y90Source/gps/histpoint 700. 0.99
> /gate/source/Y90Source/gps/histpoint 710. 0.99
> /gate/source/Y90Source/gps/histpoint 720. 0.99
> /gate/source/Y90Source/gps/histpoint 730. 0.99
> /gate/source/Y90Source/gps/histpoint 740. 1
> /gate/source/Y90Source/gps/histpoint 750. 1
> /gate/source/Y90Source/gps/histpoint 760. 1
> /gate/source/Y90Source/gps/histpoint 770. 1
> /gate/source/Y90Source/gps/histpoint 780. 1
> /gate/source/Y90Source/gps/histpoint 790. 1
> /gate/source/Y90Source/gps/histpoint 800. 1
> /gate/source/Y90Source/gps/histpoint 810. 1
> /gate/source/Y90Source/gps/histpoint 820. 1
> /gate/source/Y90Source/gps/histpoint 830. 1
> /gate/source/Y90Source/gps/histpoint 840. 1
> /gate/source/Y90Source/gps/histpoint 850. 1
> /gate/source/Y90Source/gps/histpoint 860. 1
> /gate/source/Y90Source/gps/histpoint 870. 1
> /gate/source/Y90Source/gps/histpoint 880. 1
> /gate/source/Y90Source/gps/histpoint 890. 1
> /gate/source/Y90Source/gps/histpoint 900. 1
> /gate/source/Y90Source/gps/histpoint 910. 1
> /gate/source/Y90Source/gps/histpoint 920. 1
> /gate/source/Y90Source/gps/histpoint 930. 1
> /gate/source/Y90Source/gps/histpoint 940. 0.99
> /gate/source/Y90Source/gps/histpoint 950. 0.99
> /gate/source/Y90Source/gps/histpoint 960. 0.99
> /gate/source/Y90Source/gps/histpoint 970. 0.99
> /gate/source/Y90Source/gps/histpoint 980. 0.99
> /gate/source/Y90Source/gps/histpoint 990. 0.99
> /gate/source/Y90Source/gps/histpoint 1000. 0.99
> /gate/source/Y90Source/gps/histpoint 1010. 0.98
> /gate/source/Y90Source/gps/histpoint 1020. 0.98
> /gate/source/Y90Source/gps/histpoint 1030. 0.98
> /gate/source/Y90Source/gps/histpoint 1040. 0.98
> /gate/source/Y90Source/gps/histpoint 1050. 0.98
> /gate/source/Y90Source/gps/histpoint 1060. 0.97
> /gate/source/Y90Source/gps/histpoint 1070. 0.97
> /gate/source/Y90Source/gps/histpoint 1080. 0.97
> /gate/source/Y90Source/gps/histpoint 1090. 0.96
> /gate/source/Y90Source/gps/histpoint 1100. 0.96
> /gate/source/Y90Source/gps/histpoint 1110. 0.96
> /gate/source/Y90Source/gps/histpoint 1120. 0.96
> /gate/source/Y90Source/gps/histpoint 1130. 0.95
> /gate/source/Y90Source/gps/histpoint 1140. 0.95
> /gate/source/Y90Source/gps/histpoint 1150. 0.95
> /gate/source/Y90Source/gps/histpoint 1160. 0.94
> /gate/source/Y90Source/gps/histpoint 1170. 0.94
> /gate/source/Y90Source/gps/histpoint 1180. 0.93
> /gate/source/Y90Source/gps/histpoint 1190. 0.93
> /gate/source/Y90Source/gps/histpoint 1200. 0.93
> /gate/source/Y90Source/gps/histpoint 1210. 0.92
> /gate/source/Y90Source/gps/histpoint 1220. 0.92
> /gate/source/Y90Source/gps/histpoint 1230. 0.91
> /gate/source/Y90Source/gps/histpoint 1240. 0.91
> /gate/source/Y90Source/gps/histpoint 1250. 0.9
> /gate/source/Y90Source/gps/histpoint 1260. 0.9
> /gate/source/Y90Source/gps/histpoint 1270. 0.89
> /gate/source/Y90Source/gps/histpoint 1280. 0.89
> /gate/source/Y90Source/gps/histpoint 1290. 0.88
> /gate/source/Y90Source/gps/histpoint 1300. 0.88
> /gate/source/Y90Source/gps/histpoint 1310. 0.87
> /gate/source/Y90Source/gps/histpoint 1320. 0.87
> /gate/source/Y90Source/gps/histpoint 1330. 0.86
> /gate/source/Y90Source/gps/histpoint 1340. 0.85
> /gate/source/Y90Source/gps/histpoint 1350. 0.85
> /gate/source/Y90Source/gps/histpoint 1360. 0.84
> /gate/source/Y90Source/gps/histpoint 1370. 0.84
> /gate/source/Y90Source/gps/histpoint 1380. 0.83
> /gate/source/Y90Source/gps/histpoint 1390. 0.82
> /gate/source/Y90Source/gps/histpoint 1400. 0.82
> /gate/source/Y90Source/gps/histpoint 1410. 0.81
> /gate/source/Y90Source/gps/histpoint 1420. 0.8
> /gate/source/Y90Source/gps/histpoint 1430. 0.79
> /gate/source/Y90Source/gps/histpoint 1440. 0.79
> /gate/source/Y90Source/gps/histpoint 1450. 0.78
> /gate/source/Y90Source/gps/histpoint 1460. 0.77
> /gate/source/Y90Source/gps/histpoint 1470. 0.76
> /gate/source/Y90Source/gps/histpoint 1480. 0.76
> /gate/source/Y90Source/gps/histpoint 1490. 0.75
> /gate/source/Y90Source/gps/histpoint 1500. 0.74
> /gate/source/Y90Source/gps/histpoint 1510. 0.73
> /gate/source/Y90Source/gps/histpoint 1520. 0.72
> /gate/source/Y90Source/gps/histpoint 1530. 0.71
> /gate/source/Y90Source/gps/histpoint 1540. 0.71
> /gate/source/Y90Source/gps/histpoint 1550. 0.7
> /gate/source/Y90Source/gps/histpoint 1560. 0.69
> /gate/source/Y90Source/gps/histpoint 1570. 0.68
> /gate/source/Y90Source/gps/histpoint 1580. 0.67
> /gate/source/Y90Source/gps/histpoint 1590. 0.66
> /gate/source/Y90Source/gps/histpoint 1600. 0.65
> /gate/source/Y90Source/gps/histpoint 1610. 0.64
> /gate/source/Y90Source/gps/histpoint 1620. 0.63
> /gate/source/Y90Source/gps/histpoint 1630. 0.62
> /gate/source/Y90Source/gps/histpoint 1640. 0.61
> /gate/source/Y90Source/gps/histpoint 1650. 0.6
> /gate/source/Y90Source/gps/histpoint 1660. 0.59
> /gate/source/Y90Source/gps/histpoint 1670. 0.58
> /gate/source/Y90Source/gps/histpoint 1680. 0.57
> /gate/source/Y90Source/gps/histpoint 1690. 0.56
> /gate/source/Y90Source/gps/histpoint 1700. 0.55
> /gate/source/Y90Source/gps/histpoint 1710. 0.54
> /gate/source/Y90Source/gps/histpoint 1720. 0.53
> /gate/source/Y90Source/gps/histpoint 1730. 0.52
> /gate/source/Y90Source/gps/histpoint 1740. 0.51
> /gate/source/Y90Source/gps/histpoint 1750. 0.5
> /gate/source/Y90Source/gps/histpoint 1760. 0.48
> /gate/source/Y90Source/gps/histpoint 1770. 0.47
> /gate/source/Y90Source/gps/histpoint 1780. 0.46
> /gate/source/Y90Source/gps/histpoint 1790. 0.45
> /gate/source/Y90Source/gps/histpoint 1800. 0.44
> /gate/source/Y90Source/gps/histpoint 1810. 0.43
> /gate/source/Y90Source/gps/histpoint 1820. 0.42
> /gate/source/Y90Source/gps/histpoint 1830. 0.41
> /gate/source/Y90Source/gps/histpoint 1840. 0.39
> /gate/source/Y90Source/gps/histpoint 1850. 0.38
> /gate/source/Y90Source/gps/histpoint 1860. 0.37
> /gate/source/Y90Source/gps/histpoint 1870. 0.36
> /gate/source/Y90Source/gps/histpoint 1880. 0.35
> /gate/source/Y90Source/gps/histpoint 1890. 0.34
> /gate/source/Y90Source/gps/histpoint 1900. 0.33
> /gate/source/Y90Source/gps/histpoint 1910. 0.31
> /gate/source/Y90Source/gps/histpoint 1920. 0.3
> /gate/source/Y90Source/gps/histpoint 1930. 0.29
> /gate/source/Y90Source/gps/histpoint 1940. 0.28
> /gate/source/Y90Source/gps/histpoint 1950. 0.27
> /gate/source/Y90Source/gps/histpoint 1960. 0.26
> /gate/source/Y90Source/gps/histpoint 1970. 0.25
> /gate/source/Y90Source/gps/histpoint 1980. 0.24
> /gate/source/Y90Source/gps/histpoint 1990. 0.22
> /gate/source/Y90Source/gps/histpoint 2000. 0.21
> /gate/source/Y90Source/gps/histpoint 2010. 0.2
> /gate/source/Y90Source/gps/histpoint 2020. 0.19
> /gate/source/Y90Source/gps/histpoint 2030. 0.18
> /gate/source/Y90Source/gps/histpoint 2040. 0.17
> /gate/source/Y90Source/gps/histpoint 2050. 0.16
> /gate/source/Y90Source/gps/histpoint 2060. 0.15
> /gate/source/Y90Source/gps/histpoint 2070. 0.14
> /gate/source/Y90Source/gps/histpoint 2080. 0.13
> /gate/source/Y90Source/gps/histpoint 2090. 0.12
> /gate/source/Y90Source/gps/histpoint 2100. 0.11
> /gate/source/Y90Source/gps/histpoint 2110. 0.1
> /gate/source/Y90Source/gps/histpoint 2120. 0.09
> /gate/source/Y90Source/gps/histpoint 2130. 0.08
> /gate/source/Y90Source/gps/histpoint 2140. 0.07
> /gate/source/Y90Source/gps/histpoint 2150. 0.06
> /gate/source/Y90Source/gps/histpoint 2160. 0.06
> /gate/source/Y90Source/gps/histpoint 2170. 0.05
> /gate/source/Y90Source/gps/histpoint 2180. 0.04
> /gate/source/Y90Source/gps/histpoint 2190. 0.03
> /gate/source/Y90Source/gps/histpoint 2200. 0.02
> /gate/source/Y90Source/gps/histpoint 2210. 0.02
> /gate/source/Y90Source/gps/histpoint 2220. 0.01
> /gate/source/Y90Source/gps/histpoint 2230. 0
>
> /gate/source/Y90source/gps/number 1
> /gate/source/Y90Source/setForcedUnstableFlag true
> /gate/source/Y90soutce/setForcedHalfLife 100. s
> /gate/source/Y90soutce/gps/type Point
>
> Nick Forwood
> Scientific Officer
> Dept of Nuclear Medicine
> St Vincent's Hospital
>
> PH: (02) 8382 2659
> FAX: (02) 8382 2619
>
> **********************************************************************
> This email and any files transmitted with it are confidential and
> intended solely for the use of the individual or entity to whom they
> are addressed. If you have received this email in error please notify
> the system manager.
>
> This footnote also confirms that this email message has been virus
> scanned and although no viruses were detected by the system, St Vincents &
> Mater Health Sydney accepts no liability for any consequential damage
> resulting from email containing any computer viruses.
>
> **********************************************************************
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.healthgrid.org
> http://lists.healthgrid.org/mailman/listinfo/gate-users
>
More information about the Gate-users
mailing list