[Gate-users] Gate-users Digest, Vol 116, Issue 24
gaojuan at ihep.ac.cn
gaojuan at ihep.ac.cn
Mon Jun 13 04:34:06 CEST 2016
Hello,
I used VoxelizedPhantom and VoxelizedSource to simulation, and then I cannot find Compton interactions and Rayleigh interactions in phantoms. Who can tell me Why? The next is the Physics:
#
# P H Y S I C S
#
/gate/physics/addProcess PhotoElectric
/gate/physics/processes/PhotoElectric/setModel StandardModel
/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
# Cuts for particle in NEMACylinder
/gate/physics/Electron/SetCutInRegion LYSO 1.0 cm
/gate/physics/Positron/SetCutInRegion LYSO 1.0 cm
/gate/physics/Gamma/SetCutInRegion hof_brain 0.1 mm
/gate/physics/Electron/SetCutInRegion hof_brain 0.1 mm
/gate/physics/Positron/SetCutInRegion hof_brain 0.1 mm
/gate/physics/SetMaxStepSizeInRegion hof_brain 0.01 mm
=================================================
高娟
核技术应用研究中心
地址：主楼C226室
电话：88235869
=================================================
Juan Gao
Institute of High Energy Physics
Chinese of Academy of Sciences
Beijing, China
tel:86-010-88236364
From: gate-users-request
Date: 2016-01-26 21:22
To: gate-users
Subject: Gate-users Digest, Vol 116, Issue 24
Send Gate-users mailing list submissions to
gate-users at lists.opengatecollaboration.org
To subscribe or unsubscribe via the World Wide Web, visit
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
or, via email, send a message with subject or body 'help' to
gate-users-request at lists.opengatecollaboration.org
You can reach the person managing the list at
gate-users-owner at lists.opengatecollaboration.org
When replying, please edit your Subject line so it is more specific
than "Re: Contents of Gate-users digest..."
Today's Topics:
1. Fwd: Problem with Gate simulation (David Sarrut)
----------------------------------------------------------------------
Message: 1
Date: Tue, 26 Jan 2016 15:42:21 +0100
From: David Sarrut <David.Sarrut at creatis.insa-lyon.fr>
To: gate-users <gate-users at lists.opengatecollaboration.org>
Cc: mart.vallieres at gmail.com, ielnaqa at med.umich.edu
Subject: [Gate-users] Fwd: Problem with Gate simulation
Message-ID:
<CAKTJC4AaygxCSRBki7vd3SeV3Wx5f5NsKKXj1w7k8v0HMn92fw at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
Hello,
hard to say ... It seems fine for me. Except: that we recommend to use the
physic list builder
http://wiki.opengatecollaboration.org/index.php/Users_Guide_V7.1:Setting_up_the_physics#New_physics_list_mechanism
instead of the manual phys list.
To hunt the bug, try to start with a very very simple macro that works,
than add progressively new elements until it fails.
David
On Tue, Jan 26, 2016 at 1:58 AM, Martin Valli?res <mart.vallieres at gmail.com>
wrote:
> Hi Dr. Sarrut,
>
> For some unknown reason, the email I sent to the GATE users list gets
> blocked and not sent, so this is why I am sending this email to you only
> until I solve the mail issue.
>
> I am now having a huge segmentation fault using GATE, and I have not much
> ideas where to start to debug that. Attached are some files that may help
> in describing the bug.
>
> I would so much appreciate if you could help me, or point me in the good
> direction.
>
> Best regards,
> -------------------------------------------------------------
> Martin Valli?res, PhD candidate
> Medical Physics Unit, McGill University
> martin.carrier-vallieres at mail.mcgill.ca
> 1-514-839-0248
>
> Cedars Cancer Centre, DS1.7116
> McGill University Health Centre - Glen Campus
> 1001 boulevard D?carie
> Montreal (QC), H4A 3J1
> -------------------------------------------------------------
>
--
David Sarrut, Phd
Directeur de recherche CNRS
CREATIS, UMR CNRS 5220, Inserm U1206
Centre de lutte contre le cancer L?on B?rard
28 rue La?nnec, 69373 Lyon cedex 08
Tel : 04 78 78 51 51 / 06 74 72 05 42
http://www.creatis.insa-lyon.fr/~dsarrut
_________________________________
"2 + 2 = 5, for extremely large values of 2"
_________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20160126/0aad19fa/attachment.html>
-------------- next part --------------
*** Break *** segmentation violation
===========================================================
There was a crash.
This is the entire stack trace of all threads:
===========================================================
#0 0x000000307e2ac61e in waitpid () from /lib64/libc.so.6
#1 0x000000307e23e609 in do_system () from /lib64/libc.so.6
#2 0x00002b5e83070448 in TUnixSystem::StackTrace() () from /gs/project/bnq-630-aa/martin/PET_SIM/root_v5.34/lib/libCore.so
#3 0x00002b5e8306f2c3 in TUnixSystem::DispatchSignals(ESignals) () from /gs/project/bnq-630-aa/martin/PET_SIM/root_v5.34/lib/libCore.so
#4 <signal handler called>
#5 0x000000308869d4eb in std::basic_string<char, std::char_traits<char>, std::allocator<char> >::assign(std::basic_string<char, std::char_traits<char>, std::allocator<char> > const&) () from /usr/lib64/libstdc++.so.6
#6 0x00000000007318e3 in GateAnalysis::RecordEndOfEvent(G4Event const*) ()
#7 0x0000000000725fba in GateOutputMgr::RecordEndOfEvent(G4Event const*) ()
#8 0x0000000000774ae9 in GateEventAction::EndOfEventAction(G4Event const*) ()
#9 0x00002b5e7f2e6689 in G4EventManager::DoProcessing(G4Event*) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4event.so
#10 0x00002b5e7f08ef17 in G4RunManager::ProcessOneEvent(int) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4run.so
#11 0x00002b5e7f08e5d5 in G4RunManager::DoEventLoop(int, char const*, int) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4run.so
#12 0x00002b5e7f08e56c in G4RunManager::BeamOn(int, char const*, int) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4run.so
#13 0x000000000084a1d2 in GateApplicationMgr::StartDAQ() ()
#14 0x000000000089e9c8 in GateApplicationMgrMessenger::SetNewValue(G4UIcommand*, G4String) ()
#15 0x00002b5e8224c2cb in G4UIcommand::DoIt(G4String) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#16 0x00002b5e822596b7 in G4UImanager::ApplyCommand(char const*) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#17 0x00002b5e8223eaa7 in G4UIbatch::ExecCommand(G4String const&) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#18 0x00002b5e8223fa3b in G4UIbatch::SessionStart() () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#19 0x00002b5e82257313 in G4UImanager::ExecuteMacroFile(char const*) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#20 0x00002b5e82252658 in G4UIcontrolMessenger::SetNewValue(G4UIcommand*, G4String) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#21 0x00002b5e8224c2cb in G4UIcommand::DoIt(G4String) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#22 0x00002b5e822596b7 in G4UImanager::ApplyCommand(char const*) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#23 0x00000000004ef83d in main ()
===========================================================
The lines below might hint at the cause of the crash.
If they do not help you then please submit a bug report at
http://root.cern.ch/bugs. Please post the ENTIRE stack trace
from above as an attachment in addition to anything else
that might help us fixing this issue.
===========================================================
#5 0x000000308869d4eb in std::basic_string<char, std::char_traits<char>, std::allocator<char> >::assign(std::basic_string<char, std::char_traits<char>, std::allocator<char> > const&) () from /usr/lib64/libstdc++.so.6
#6 0x00000000007318e3 in GateAnalysis::RecordEndOfEvent(G4Event const*) ()
#7 0x0000000000725fba in GateOutputMgr::RecordEndOfEvent(G4Event const*) ()
#8 0x0000000000774ae9 in GateEventAction::EndOfEventAction(G4Event const*) ()
#9 0x00002b5e7f2e6689 in G4EventManager::DoProcessing(G4Event*) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4event.so
#10 0x00002b5e7f08ef17 in G4RunManager::ProcessOneEvent(int) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4run.so
#11 0x00002b5e7f08e5d5 in G4RunManager::DoEventLoop(int, char const*, int) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4run.so
#12 0x00002b5e7f08e56c in G4RunManager::BeamOn(int, char const*, int) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4run.so
#13 0x000000000084a1d2 in GateApplicationMgr::StartDAQ() ()
#14 0x000000000089e9c8 in GateApplicationMgrMessenger::SetNewValue(G4UIcommand*, G4String) ()
#15 0x00002b5e8224c2cb in G4UIcommand::DoIt(G4String) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#16 0x00002b5e822596b7 in G4UImanager::ApplyCommand(char const*) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#17 0x00002b5e8223eaa7 in G4UIbatch::ExecCommand(G4String const&) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#18 0x00002b5e8223fa3b in G4UIbatch::SessionStart() () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#19 0x00002b5e82257313 in G4UImanager::ExecuteMacroFile(char const*) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#20 0x00002b5e82252658 in G4UIcontrolMessenger::SetNewValue(G4UIcommand*, G4String) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#21 0x00002b5e8224c2cb in G4UIcommand::DoIt(G4String) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#22 0x00002b5e822596b7 in G4UImanager::ApplyCommand(char const*) () from /gs/project/bnq-630-aa/martin/PET_SIM/geant4.10.01-install/lib64/libG4intercoms.so
#23 0x00000000004ef83d in main ()
===========================================================
-------------- next part --------------
[G4]
[G4] *************************************************************
[G4] Geant4 version Name: geant4-10-01 (5-December-2014)
[G4] Copyright : Geant4 Collaboration
[G4] Reference : NIM A 506 (2003), 250-303
[G4] WWW : http://cern.ch/geant4
[G4] *************************************************************
[G4]
[Core-0] Initialization of geometry
[Core-0] Initialization of physics
[Core-0] Initialization of actors
[Core-0]
[Core-0] **********************************************************************
[Core-0] GATE version name: gate_v7.1
[Core-0] Copyright : OpenGATE Collaboration
[Core-0] Reference : Phys. Med. Biol. 49 (2004) 4543-4561
[Core-0] Reference : Phys. Med. Biol. 56 (2011) 881-901
[Core-0] WWW : http://www.opengatecollaboration.org
[Core-0] **********************************************************************
[Core-0]
[Core-0] You are using Geant4 version 10.1.0
/control/alias simNumber 1
/control/alias simTimeSlice 0.083333
/control/alias simTimeStart 0.000000
/control/alias simTimeStop 0.083333
/control/alias seed 733251318
[Core-0] Starting macro simuMacro.mac
CREATED THE COMMAND /gate/systems/cylindricalPET/addAnewRsector
### WARNING ### setSkipEqualMaterials at false !! The Geant4 method is not safe since the release 9.5 - Need to be fixed
### WARNING ### setSkipEqualMaterials at false !! The Geant4 method is not safe since the release 9.5 - Need to be fixed
List of Enabled processes:
=== Process: RayleighScattering ===
Particle: gamma
* Model(s):
- PenelopeModel
=== Process: eMultipleScattering ===
Particle: e+
* Model(s):
- <!> *** Warning *** <!> No model selected!
Particle: e-
* Model(s):
- <!> *** Warning *** <!> No model selected!
=== Process: Compton ===
Particle: gamma
* Model(s):
- StandardModel
=== Process: ElectronIonisation ===
Particle: e+
* Model(s):
- StandardModel
Particle: e-
* Model(s):
- StandardModel
=== Process: Bremsstrahlung ===
Particle: e+
* Model(s):
- StandardModel
Particle: e-
* Model(s):
- StandardModel
=== Process: PositronAnnihilation ===
Particle: e+
=== Process: PhotoElectric ===
Particle: gamma
* Model(s):
- StandardModel
List of initialized processes:
[Core-0] Initialization of geometry
Header read from '/gs/project/bnq-630-aa/martin/SIMU_PET/test/simuPhantom.h33'
Data file name '/gs/project/bnq-630-aa/martin/SIMU_PET/test/simuPhantom.i33'
Nb of planes: 47
Nb of pixels per plane: 91 91
Pixel size: 5.46875 5.46875
Slice thickness: 3.27
Matrix size: 497.656 497.656
Data type: FLOAT
Bytes per pixel: 4
Data byte order: LITTLEENDIAN
Machine byte order: LITTLEENDIAN
[Core-0] Initialization of physics
[Physic-0] WARNING: manual physic lists are being deprecated.
Please, use physic list builder mechanism instead. Related documentation can be found at:
http://wiki.opengatecollaboration.org/index.php/Users_Guide_V7.0:Setting_up_the_physics
[Core-0] Initialization of actors
------------------------------------------------------------------------------------------------
WARNING: Macro commands related to voxelized source description have been modified in GATE V7.1.
Older ones are being deprecated and will be removed from the next release.
Please, have a look at the related documentation at:
http://wiki.opengatecollaboration.org/index.php/Users_Guide_V7.1:Voxelized_Source_and_Phantom
------------------------------------------------------------------------------------------------
GateSourceVoxelInterfileReader::ReadFile : fileName: /gs/project/bnq-630-aa/martin/SIMU_PET/test/simuSource.h33
Header read from '/gs/project/bnq-630-aa/martin/SIMU_PET/test/simuSource.h33'
Data file name '/gs/project/bnq-630-aa/martin/SIMU_PET/test/simuSource.i33'
Nb of planes: 47
Nb of pixels per plane: 46 46
Pixel size: 5.46875 5.46875
Slice thickness: 3.27
Matrix size: 251.562 251.562
Data type: UNSIGNED INTEGER
Bytes per pixel: 2
Data byte order: LITTLEENDIAN
Machine byte order: LITTLEENDIAN
[Acquisition-0]
[Acquisition-0] ============= Source initialization =============
[Acquisition-0]
[Acquisition-0] ============= Acquisition starts! =============
[Acquisition-0] Simulation start time = 0 sec
[Acquisition-0] Simulation end time = 0.083333 sec
[Acquisition-0] Simulation will have = 1 run(s)
[Acquisition-0] Slice 0 from 0 to 0.083333 s [slice=0.083333 s], final stop at 0.083333 s.
RayleighScattering: for gamma, applyCuts: 1 SubType= 11 BuildTable= 1
Lambda table from 100 eV to 100 keV, 7 bins per decade, spline: 0
LambdaPrime table from 100 keV to 10 TeV in 56 bins
===== EM models for the G4Region DefaultRegionForTheWorld ======
PenRayleigh : Emin= 0 eV Emax= 100 GeV
Compton: for gamma, applyCuts: 1 SubType= 13 BuildTable= 1
Lambda table from 100 eV to 1 MeV, 7 bins per decade, spline: 1
LambdaPrime table from 1 MeV to 10 TeV in 49 bins
===== EM models for the G4Region DefaultRegionForTheWorld ======
Klein-Nishina : Emin= 0 eV Emax= 10 TeV
PhotoElectric: for gamma, applyCuts: 1 SubType= 12 BuildTable= 0
LambdaPrime table from 200 keV to 10 TeV in 54 bins
===== EM models for the G4Region DefaultRegionForTheWorld ======
PhotoElectric : Emin= 0 eV Emax= 10 TeV AngularGenSauterGavrila
msc: for e- SubType= 10
RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1
===== EM models for the G4Region DefaultRegionForTheWorld ======
UrbanMsc : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV
ElectronIonisation: for e- SubType= 2
dE/dx and range tables from 100 eV to 10 TeV in 77 bins
Lambda tables from threshold to 10 TeV, 7 bins per decade, spline: 1
finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
===== EM models for the G4Region DefaultRegionForTheWorld ======
MollerBhabha : Emin= 0 eV Emax= 10 TeV
Bremsstrahlung: for e- SubType= 3
dE/dx and range tables from 100 eV to 10 TeV in 77 bins
Lambda tables from threshold to 10 TeV, 7 bins per decade, spline: 1
LPM flag: 1 for E > 1 GeV
===== EM models for the G4Region DefaultRegionForTheWorld ======
eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
msc: for e+ SubType= 10
RangeFactor= 0.04, stepLimitType: 1, latDisplacement: 1
===== EM models for the G4Region DefaultRegionForTheWorld ======
UrbanMsc : Emin= 0 eV Emax= 10 TeV Table with 77 bins Emin= 100 eV Emax= 10 TeV
ElectronIonisation: for e+ SubType= 2
dE/dx and range tables from 100 eV to 10 TeV in 77 bins
Lambda tables from threshold to 10 TeV, 7 bins per decade, spline: 1
finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
===== EM models for the G4Region DefaultRegionForTheWorld ======
MollerBhabha : Emin= 0 eV Emax= 10 TeV
Bremsstrahlung: for e+ SubType= 3
dE/dx and range tables from 100 eV to 10 TeV in 77 bins
Lambda tables from threshold to 10 TeV, 7 bins per decade, spline: 1
LPM flag: 1 for E > 1 GeV
===== EM models for the G4Region DefaultRegionForTheWorld ======
eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
annihil: for e+, integral: 1 , applyCuts: 1 SubType= 5 BuildTable= 0
===== EM models for the G4Region DefaultRegionForTheWorld ======
eplus2gg : Emin= 0 eV Emax= 10 TeV
-------------- next part --------------
/vis/disable
/gate/output/verbose 2
# MATERIALS
/gate/geometry/setMaterialDatabase /gs/project/bnq-630-aa/martin/SIMU_PET/test/GateMaterials.db
## SCANNER GEOMETRY
# WORLD
/gate/world/geometry/setXLength 3.000000 mm
/gate/world/geometry/setYLength 3.000000 mm
/gate/world/geometry/setZLength 3.000000 mm
/gate/world/setMaterial Air
# SYSTEM
/gate/world/daughters/name cylindricalPET
/gate/world/daughters/insert cylinder
/gate/cylindricalPET/geometry/setRmin 443.000000 mm
/gate/cylindricalPET/geometry/setRmax 503.000000 mm
/gate/cylindricalPET/geometry/setHeight 157.000800 mm
/gate/cylindricalPET/setMaterial Air
/gate/cylindricalPET/vis/forceWireframe
# RSECTOR
/gate/cylindricalPET/daughters/name RSECTOR
/gate/cylindricalPET/daughters/insert box
/gate/RSECTOR/placement/setTranslation 458.000000 0 0 mm
/gate/RSECTOR/geometry/setXLength 30.000000 mm
/gate/RSECTOR/geometry/setYLength 78.258700 mm
/gate/RSECTOR/geometry/setZLength 157.000800 mm
/gate/RSECTOR/setMaterial Glass
/gate/RSECTOR/vis/forceWireframe
/gate/systems/cylindricalPET/rsector/attach RSECTOR
# MODULE
/gate/RSECTOR/daughters/name MODULE
/gate/RSECTOR/daughters/insert box
/gate/MODULE/geometry/setXLength 30.000000 mm
/gate/MODULE/geometry/setYLength 39.008500 mm
/gate/MODULE/geometry/setZLength 39.008500 mm
/gate/MODULE/setMaterial PTFE
/gate/MODULE/vis/forceWireframe
/gate/systems/cylindricalPET/module/attach MODULE
# CRYSTAL
/gate/MODULE/daughters/name CRYSTAL
/gate/MODULE/daughters/insert box
/gate/CRYSTAL/geometry/setXLength 30.000000 mm
/gate/CRYSTAL/geometry/setYLength 6.300000 mm
/gate/CRYSTAL/geometry/setZLength 6.300000 mm
/gate/CRYSTAL/setMaterial BGO
/gate/CRYSTAL/vis/forceWireframe
/gate/systems/cylindricalPET/crystal/attach CRYSTAL
# REPEAT CRYSTAL
/gate/CRYSTAL/repeaters/insert cubicArray
/gate/CRYSTAL/cubicArray/setRepeatNumberX 1
/gate/CRYSTAL/cubicArray/setRepeatNumberY 6
/gate/CRYSTAL/cubicArray/setRepeatNumberZ 6
/gate/CRYSTAL/cubicArray/setRepeatVector 0. 6.541700 6.541700 mm
# REPEAT MODULE
/gate/MODULE/repeaters/insert cubicArray
/gate/MODULE/cubicArray/setRepeatNumberX 1
/gate/MODULE/cubicArray/setRepeatNumberY 2
/gate/MODULE/cubicArray/setRepeatNumberZ 4
/gate/MODULE/cubicArray/setRepeatVector 0. 39.250200 39.250200 mm
# REPEAT RSECTOR
/gate/RSECTOR/repeaters/insert ring
/gate/RSECTOR/ring/setRepeatNumber 35
# ATTACH CRYSTAL SD
/gate/CRYSTAL/attachCrystalSD
## PHANTOM GEOMETRY: CT VOLUME
/gate/geometry/setMaterialDatabase /gs/project/bnq-630-aa/martin/SIMU_PET/test/patient-HUmaterials.db
/gate/world/daughters/name voxelizedPhantom
/gate/world/daughters/insert ImageRegularParametrisedVolume
/gate/voxelizedPhantom/geometry/setImage /gs/project/bnq-630-aa/martin/SIMU_PET/test/simuPhantom.h33
/gate/voxelizedPhantom/geometry/setHUToMaterialFile /gs/project/bnq-630-aa/martin/SIMU_PET/test/patient-HU2mat.txt
/gate/voxelizedPhantom/setSkipEqualMaterials 0
/gate/voxelizedPhantom/placement/setTranslation 0.000000 -0.000000 0.000000 mm
/gate/voxelizedPhantom/attachPhantomSD
## PHYSICS
/gate/physics/addProcess PhotoElectric
/gate/physics/processes/PhotoElectric/setModel StandardModel
/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 eMultipleScattering e+
/gate/physics/addProcess eMultipleScattering e-
/gate/physics/processList Enabled
/gate/physics/processList Initialized
## INITIALIZATION
/gate/run/initialize
## DIGITIZERS
/gate/digitizer/Singles/insert adder
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 2
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.150000
/gate/digitizer/Singles/blurring/setEnergyOfReference 511.000000 keV
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 350.000000 keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 650.000000 keV
/gate/digitizer/Coincidences/setWindow 10.000000 ns
## SOURCE
/gate/source/addSource voxelizedSource voxel
/gate/source/voxelizedSource/reader/insert interfile
/gate/source/voxelizedSource/interfileReader/translator/insert linear
/gate/source/voxelizedSource/interfileReader/linearTranslator/setScale 1 Bq
/gate/source/voxelizedSource/interfileReader/readFile /gs/project/bnq-630-aa/martin/SIMU_PET/test/simuSource.h33
/gate/source/voxelizedSource/setPosition -43.750000 -185.937500 -76.845000 mm
/gate/source/voxelizedSource/gps/particle e+
/gate/source/voxelizedSource/gps/energytype Fluor18
/gate/source/voxelizedSource/setForcedUnstableFlag true
/gate/source/voxelizedSource/setForcedHalfLife 6586.2 s
/gate/source/voxelizedSource/gps/angtype iso
## OUTPUTS
# BINARY OUTPUT FORMAT
/gate/output/root/enable
/gate/output/root/setFileName /gs/project/bnq-630-aa/martin/SIMU_PET/test/simu{simNumber}
/gate/output/root/setRootHitFlag 0
/gate/output/root/setRootSinglesFlag 0
/gate/output/root/setRootCoincidencesFlag 1
## EXPERIMENT
# SIMULATION TIME
/gate/application/setTimeSlice {simTimeSlice} s
/gate/application/setTimeStart {simTimeStart} s
/gate/application/setTimeStop {simTimeStop} s
# RANDOM GENERATOR
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed {seed}
# START ACQUISITION
/gate/application/startDAQ
exit
------------------------------
_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
End of Gate-users Digest, Vol 116, Issue 24
*******************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20160613/1baebe99/attachment-0001.html>
More information about the Gate-users
mailing list