[Gate-users] statistical error

Ismet Isnaini ismeth74 at yahoo.com
Sat Jul 28 01:15:12 CEST 2012 

Dear GATE users,

I am trying to run my simulation on sensitivity profile of PET machine, but the result has great statistical error and distorted shape of triangular .

has any body got any idea what went wrong? 

below is the physics process I have, I using a line source without phantom

regards





#-----------------------------------------------------------------#
#                                                                 #
#     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   P H Y S I C S                            #
#                                                                 #
#-----------------------------------------------------------------#



########################################
#
#       NEW          P H Y S I C S
#
/gate/physics/addProcess PhotoElectric
/gate/physics/processes/PhotoElectric/setModel StandardModel

/gate/physics/addProcess Compton gamma

/gate/physics/processes/Compton/setModel StandardModel

/gate/physics/addProcess GammaConversion

/gate/physics/processes/GammaConversion/setModel StandardModel 


/gate/physics/addProcess RayleighScattering gamma
/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 e+

/gate/physics/processes/PositronAnnihilation/setModel StandardModel

/gate/physics/addProcess eMultipleScattering e+
/gate/physics/addProcess eMultipleScattering e-

/gate/physics/processList Enabled
/gate/physics/processList Initialized


/gate/physics/Gamma/SetCutInRegion      LAYER1 1.0 cm
/gate/physics/Electron/SetCutInRegion   LAYER1 1.0 cm
/gate/physics/Positron/SetCutInRegion   LAYER1 1.0 cm
/gate/physics/Gamma/SetCutInRegion      LAYER2 1.0 cm
/gate/physics/Electron/SetCutInRegion   LAYER2 1.0 cm
/gate/physics/Positron/SetCutInRegion   LAYER2 1.0 cm
/gate/physics/Gamma/SetCutInRegion      LAYER3 1.0 cm
/gate/physics/Electron/SetCutInRegion   LAYER3 1.0 cm
/gate/physics/Positron/SetCutInRegion   LAYER3 1.0 cm
/gate/physics/Gamma/SetCutInRegion      LAYER4 1.0 cm
/gate/physics/Electron/SetCutInRegion   LAYER4 1.0 cm
/gate/physics/Positron/SetCutInRegion   LAYER4 1.0 cm


#-----------------------------------------------------------------#
#                                                                 #
#     I N I T I A L I Z A T I O N   O F   S I M U L A T I O N     #
#                                                                 #
#-----------------------------------------------------------------#

#/gate/systems/cylindricalPET/verbose 0
#/gate/geometry/enableAutoUpdate
/gate/run/initialize

#/geometry/test/recursive_test


#-----------------------------------------------------------------#
#                                                                 #
#        D E F I N I T I O N   O F   A C Q U I S I T I O N        #
#  D I G I T I Z E R   &   C O I N C I D E N C E   S O R T E R    #
#                                                                 #
#-----------------------------------------------------------------#

#A D D E R
/gate/digitizer/Singles/insert adder


#R E A D O U T *******************************
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 3
######for depth 3:submodule ,2:module, 1:rsector, 0:cylinder

#E N E R G Y    B L U R R I N G
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.15
/gate/digitizer/Singles/blurring/setEnergyOfReference 511. keV

/gate/digitizer/Singles/insert timeResolution
/gate/digitizer/Singles/timeResolution/setTimeResolution 900 ps

# THRESHOLDER
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 400. keV
/gate/digitizer/Singles/thresholder/verbose 0

/gate/digitizer/Singles/insert deadtime
/gate/digitizer/Singles/deadtime/setDeadTime 250. ns
/gate/digitizer/Singles/deadtime/setMode paralysable
/gate/digitizer/Singles/deadtime/chooseDTVolume submodule

# UPHOLDER
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 600. keV

#       D E A D    T I M E 

/gate/digitizer/Singles/name deadtime_cassette
/gate/digitizer/Singles/insert deadtime
/gate/digitizer/Singles/deadtime_cassette/setDeadTime 256. ns
/gate/digitizer/Singles/deadtime_cassette/setMode nonparalysable
/gate/digitizer/Singles/deadtime_cassette/chooseDTVolume rsector

/gate/digitizer/Coincidences/setWindow 3. ns
/gate/digitizer/Coincidences/setOffset 0. ns
/gate/digitizer/Coincidences/describe
#/gate/digitizer/Coincidences/minSectorDifference 1
/gate/digitizer/Coincidences/MultiplesPolicy takeWinnerOfGoods

#-----------------------------------------------------------------#
#                                                                 #
#       D E F I N I T I O N   O F   O U T P U T   F I L E         #
#                                                                 #
#-----------------------------------------------------------------#

#C H A N G E    T H E     S E E D (1)  O R   N O T (0)
#/gate/output/root/setSaveRndmFlag 0
/gate/random/setEngineSeed auto

#S E T U P  - R O O T  F I L E
/gate/output/root/enable
/gate/output/root/setFileName trial
/gate/output/root/setRootHitFlag 0
/gate/output/root/setRootSinglesFlag 1
/gate/output/root/setRootCoincidencesFlag 1
#/gate/output/lmf disable

#S E T U P  - A S C I I  F I L E
/gate/output/ascii/disable



#-----------------------------------------------------------------#
#                                                                 #
#    D E F I N I T I O N   O F   V E R B O S I T Y   L E V E L    #
#                                                                 #
#-----------------------------------------------------------------#

/control/verbose 0
#/grdm/verbose 0
/run/verbose 0
/event/verbose 0
/tracking/verbose 0
/gate/application/verbose 0
/gate/generator/verbose 0
#/gate/stacking/verbose 0
#/gate/event/verbose 0
/gate/source/verbose 0


#-----------------------------------------------------------------#
#                                                                 #
#          D E F I N I T I O N   O F   S O U R C E S              #
#                                                                 #
#-----------------------------------------------------------------#

/gate/source/addSource twogamma
/gate/source/twogamma/setActivity 50000. becquerel
/gate/source/twogamma/setType backtoback
/gate/source/twogamma/gps/particle gamma
/gate/source/twogamma/gps/energytype Mono
/gate/source/twogamma/gps/monoenergy 511. keV
/gate/source/twogamma/gps/centre 0. 0. 0. cm
/gate/source/twogamma/gps/angtype iso
/gate/source/twogamma/gps/type Volume
/gate/source/twogamma/gps/shape Cylinder
/gate/source/twogamma/gps/radius 3.2 mm
/gate/source/twogamma/gps/halfz 87.5 cm

#/control/execute LineSource.mac


/gate/source/list

#-----------------------------------------------------------------#
#                                                                 #
#               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                   #
#                                                                 #
#-----------------------------------------------------------------#

/random/setSavingFlag 1 

#
#S E T    T I M E    P A R A M E T E R
#
#/gate/application/setTotalNumberOfPrimaries 10000000
/gate/application/setTimeSlice  50 ms
/gate/application/setTimeStart      0.   s
/gate/application/setTimeStop   37333  ms


#
#LET'S RUN THE SIMULATION!
#
/gate/application/startDAQ


exit

More information about the Gate-users mailing list