[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