[Gate-users] Source histogram for CT

Tony Evershed a.n.z.evershed at qmul.ac.uk
Fri Jun 25 16:48:14 CEST 2010


Hi,
I’m trying to do cone-beam X-ray CT from a point source, and fell over at
the first stage (setting up a source). The source macro (called by a main
macro) uses the following code (taken pretty much wholesale from the classic
example_CT in the latest vGate and slimmed down, with the source type
changed to a point):

/gate/source/addSource xraygun
/gate/source/verbose 0
/gate/source/xraygun/setActivity 350000. becquerel
/gate/source/xraygun/gps/verbose 0
/gate/source/xraygun/gps/particle gamma
/gate/source/xraygun/gps/energytype Arb
/gate/source/xraygun/gps/histname arb
/gate/source/xraygun/gps/emin 10.00 keV
/gate/source/xraygun/gps/emax 40.00 keV

/gate/source/xraygun/gps/histpoint  0.0100  5060586
/gate/source/xraygun/gps/histpoint  0.0105  4547466
/gate/source/xraygun/gps/histpoint  0.0110  4274620
/gate/source/xraygun/gps/histpoint  0.0115  4489648
/gate/source/xraygun/gps/histpoint  0.0120  4733971
/gate/source/xraygun/gps/histpoint  0.0125  5089220
/gate/source/xraygun/gps/histpoint  0.0130  5486806
/gate/source/xraygun/gps/histpoint  0.0135  6030102
/gate/source/xraygun/gps/histpoint  0.0140  6665430
/gate/source/xraygun/gps/histpoint  0.0145  7140114
/gate/source/xraygun/gps/histpoint  0.0150  7674678
/gate/source/xraygun/gps/histpoint  0.0155  8086993
/gate/source/xraygun/gps/histpoint  0.0160  8551672
/gate/source/xraygun/gps/histpoint  0.0165  8868724
/gate/source/xraygun/gps/histpoint  0.0170  9215647
/gate/source/xraygun/gps/histpoint  0.0175  9438009
/gate/source/xraygun/gps/histpoint  0.0180  9682990
/gate/source/xraygun/gps/histpoint  0.0185  9834746
/gate/source/xraygun/gps/histpoint  0.0190  9998208
/gate/source/xraygun/gps/histpoint  0.0195  10080460
/gate/source/xraygun/gps/histpoint  0.0200  10180510
/gate/source/xraygun/gps/histpoint  0.0205  10205130
/gate/source/xraygun/gps/histpoint  0.0210  10241880
/gate/source/xraygun/gps/histpoint  0.0215  10223500
/gate/source/xraygun/gps/histpoint  0.0220   10216490
/gate/source/xraygun/gps/histpoint  0.0225  10156620
/gate/source/xraygun/gps/histpoint  0.0230  10105940
/gate/source/xraygun/gps/histpoint  0.0235  10022750
/gate/source/xraygun/gps/histpoint  0.0240  9940856
/gate/source/xraygun/gps/histpoint  0.0245  9828219
/gate/source/xraygun/gps/histpoint  0.0250  9723552
/gate/source/xraygun/gps/histpoint  0.0255  9603177
/gate/source/xraygun/gps/histpoint  0.0260  9485084
/gate/source/xraygun/gps/histpoint  0.0265  9349130
/gate/source/xraygun/gps/histpoint  0.0270  9218213
/gate/source/xraygun/gps/histpoint  0.0275  9078410
/gate/source/xraygun/gps/histpoint  0.0280  8940065
/gate/source/xraygun/gps/histpoint  0.0285  8793341
/gate/source/xraygun/gps/histpoint  0.0290  8654856
/gate/source/xraygun/gps/histpoint  0.0295  8503122
/gate/source/xraygun/gps/histpoint  0.0300  8357768
/gate/source/xraygun/gps/histpoint  0.0305  8207681
/gate/source/xraygun/gps/histpoint  0.0310  8067777
/gate/source/xraygun/gps/histpoint  0.0315  7917228
/gate/source/xraygun/gps/histpoint  0.0320  7772499
/gate/source/xraygun/gps/histpoint  0.0325  7629667
/gate/source/xraygun/gps/histpoint  0.0330  7487478
/gate/source/xraygun/gps/histpoint  0.0335  7345984
/gate/source/xraygun/gps/histpoint  0.0340  7209631
/gate/source/xraygun/gps/histpoint  0.0345  7071618
/gate/source/xraygun/gps/histpoint  0.0350  6935182
/gate/source/xraygun/gps/histpoint  0.0355  6797408
/gate/source/xraygun/gps/histpoint  0.0360  6662980
/gate/source/xraygun/gps/histpoint  0.0365  6534767
/gate/source/xraygun/gps/histpoint  0.0370  6407531
/gate/source/xraygun/gps/histpoint  0.0375  6280954
/gate/source/xraygun/gps/histpoint  0.0380  6159509
/gate/source/xraygun/gps/histpoint  0.0385  6036038
/gate/source/xraygun/gps/histpoint  0.0390  5916583
/gate/source/xraygun/gps/histpoint  0.0395  5798614
/gate/source/xraygun/gps/histpoint  0.0400  5684363
/gate/source/xraygun/gps/arbint Lin
/gate/source/xraygun/gps/type Point
/gate/source/xraygun/gps/mintheta 0  deg
/gate/source/xraygun/gps/maxtheta  6.8 deg
/gate/source/xraygun/gps/centre 0. 0. -15.0 cm
/gate/source/xraygun/gps/angtype iso
/gate/source/list

This raises the following error:

: double value expected.6
> ***** [G4-cerr] ***** Illegal parameter (1)
</gate/source/xraygun/gps/histpoint 0.0100 506586
[Gate] Sorry, error in a macro command : abort.

Am I missing something really obvious here that’s leading to this error?

Many thanks,
Tony Evershed



More information about the Gate-users mailing list