[Gate-users] coincidences sorter /depth 5
Ashley Devarim
ashley_devarim at hotmail.com
Sat Jan 12 16:33:51 CET 2008
Hi Simon,
Thank you very much for your answer ....
Here the Macro (I added boxa (module) and boxb (submodule) in purpose even if they are useless but the problem persists when I specify coincidences sorter /depth 5) I take for the readout the depth value 5 and it's OK !!!!! even with 3 system levels, but not for coincidence sorter....
######################################################################################################################
###### VISUALIZATION ########
/vis/open OGLSX
/vis/viewer/reset
/vis/viewer/set/viewpointThetaPhi 60 60
/vis/viewer/zoom 3
/vis/viewer/set/style surface
/vis/drawVolume
/tracking/storeTrajectory 1
/vis/scene/endOfEventAction accumulate
/vis/viewer/update
/gate/geometry/enableAutoUpdate
#### WORLD #####
/gate/world/geometry/setXLength 50. cm
/gate/world/geometry/setYLength 50. cm
/gate/world/geometry/setZLength 50. cm
#### System ####
/gate/world/daughters/name cylindricalPET
/gate/world/daughters/insert cylinder
/gate/cylindricalPET/setMaterial Air
/gate/cylindricalPET/geometry/setRmax 95 mm
/gate/cylindricalPET/geometry/setRmin 80 mm
/gate/cylindricalPET/geometry/setHeight 50 mm
/gate/cylindricalPET/vis/setColor white
/gate/cylindricalPET/vis/forceWireframe
/gate/cylindricalPET/daughters/name box1
/gate/cylindricalPET/daughters/insert box
/gate/box1/placement/setTranslation 85. 0 0 mm
/gate/box1/geometry/setXLength 10 mm
/gate/box1/geometry/setYLength 24 mm
/gate/box1/geometry/setZLength 4.4 mm
/gate/box1/setMaterial Air
/gate/box1/vis/setColor grey
/gate/box1/vis/forceWireframe
/gate/box1/daughters/name boxa
/gate/box1/daughters/insert box
/gate/boxa/geometry/setXLength 10 mm
/gate/boxa/geometry/setYLength 24 mm
/gate/boxa/geometry/setZLength 4.4 mm
/gate/boxa/setMaterial Air
/gate/boxa/vis/setColor grey
/gate/boxa/vis/forceWireframe
/gate/boxa/daughters/name boxb
/gate/boxa/daughters/insert box
/gate/boxb/geometry/setXLength 10 mm
/gate/boxb/geometry/setYLength 24 mm
/gate/boxb/geometry/setZLength 4.4 mm
/gate/boxb/setMaterial Air
/gate/boxb/vis/setColor grey
/gate/boxb/vis/forceWireframe
/gate/boxb/daughters/name box2
/gate/boxb/daughters/insert box
/gate/box2/geometry/setXLength 10 mm
/gate/box2/geometry/setYLength 8 mm
/gate/box2/geometry/setZLength 4.4 mm
/gate/box2/setMaterial Air
/gate/box2/vis/setColor grey
/gate/box2/vis/forceWireframe
/gate/box2/daughters/name box3
/gate/box2/daughters/insert box
/gate/box3/placement/setTranslation 0 0 -1.1 mm
/gate/box3/geometry/setXLength 10. mm
/gate/box3/geometry/setYLength 8 mm
/gate/box3/geometry/setZLength 2.2 mm
/gate/box3/setMaterial LSO
/gate/box3/vis/setColor white
/gate/box3/vis/forceSolid
/gate/box2/daughters/name box4
/gate/box2/daughters/insert box
/gate/box4/placement/setTranslation 0 0 1.1 mm
/gate/box4/geometry/setXLength 10 mm
/gate/box4/geometry/setYLength 8 mm
/gate/box4/geometry/setZLength 2.2 mm
/gate/box4/setMaterial BGO
/gate/box4/vis/setColor green
/gate/box4/vis/forceSolid
/gate/box1/repeaters/insert ring
/gate/box1/ring/setRepeatNumber 15
/gate/box1/repeaters/insert linear
/gate/box1/linear/setRepeatNumber 4
/gate/box1/linear/setRepeatVector 0. 0. 8. mm
/gate/box2/repeaters/insert cubicArray
/gate/box2/cubicArray/setRepeatNumberX 1
/gate/box2/cubicArray/setRepeatNumberY 3
/gate/box2/cubicArray/setRepeatNumberZ 1
/gate/box2/cubicArray/setRepeatVector 0. 8. 0. mm
/gate/systems/cylindricalPET/rsector/attach box1
/gate/systems/cylindricalPET/module/attach boxa
/gate/systems/cylindricalPET/submodule/attach boxb
/gate/systems/cylindricalPET/crystal/attach box2
/gate/systems/cylindricalPET/layer0/attach box3
/gate/systems/cylindricalPET/layer1/attach box4
/gate/box3/attachCrystalSD
/gate/box4/attachCrystalSD
###################################################
# physics processes
/gate/physics/gamma/selectPhotoelectric standard
/gate/physics/gamma/selectCompton standard
/gate/physics/gamma/selectGammaConversion standard
/gate/physics/setElectronCut 30 m
/gate/physics/setXRayCut 1. GeV
/gate/physics/setDeltaRayCut 1. GeV
###################################################
# INITIALIZE
/gate/geometry/enableAutoUpdate
/run/initialize
##############################################
# Setting-up the digitizer
/gate/digitizer/Singles/insert adder
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 5
# Energy blurring
/gate/digitizer/Singles/insert localBlurring
/gate/digitizer/Singles/localBlurring/chooseNewVolume box3
/gate/digitizer/Singles/localBlurring/box3/setResolution 0
/gate/digitizer/Singles/localBlurring/box3/setEnergyOfReference 511. keV
/gate/digitizer/Singles/insert localBlurring
/gate/digitizer/Singles/localBlurring/chooseNewVolume box4
/gate/digitizer/Singles/localBlurring/box4/setResolution 0
/gate/digitizer/Singles/localBlurring/box4/setEnergyOfReference 511. keV
# Energy Window
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 300. keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 600. keV
# Coincidence sorter
/gate/digitizer/Coincidences/setWindow 10. ns
/gate/digitizer/Coincidences/setDepth 5
###################################################
# Source
/gate/source/addSource twogamma
/gate/source/twogamma/setActivity 100000. 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/type Point
/gate/source/twogamma/gps/confine NULL
/gate/source/twogamma/gps/angtype iso
# set min and max emission angles
/gate/source/twogamma/gps/mintheta 0. deg
/gate/source/twogamma/gps/maxtheta 180. deg
/gate/source/twogamma/gps/minphi 0. deg
/gate/source/twogamma/gps/maxphi 360. deg
/gate/source/list
####################################################
# ASCII Output format
/gate/output/ascii/setOutFileHitsFlag 1
/gate/output/ascii/setOutFileSinglesFlag 1
/gate/output/ascii/setOutFileCoincidencesFlag 1
# ROOT Output format
/gate/output/root/setFileName gatescan
/gate/output/root/setRootSinglesFlag 1
/gate/output/root/setRootCoincidencesFlag 1
/gate/output/root/setRootHitFlag 1
/gate/output/root/setRootNtupleFlag 0
/gate/application/setTimeSlice 0.1 s
/gate/application/setTimeStart 0. s
/gate/application/setTimeStop 0.2 s
# Start the acquisition
/gate/application/startDAQ
################################################################################################
I'm looking forward for your response
Thank you very much for your help
Ashley
________________________________
> Date: Sat, 12 Jan 2008 13:27:19 +0100
> From: gate.stute at gmail.com
> To: ashley_devarim at hotmail.com
> Subject: Re: [Gate-users] coincidences sorter /depth 5
> CC: gate-users at lists.healthgrid.org
>
> Hi,
>
>
> 2008/1/12, Ashley Devarim <ashley_devarim at hotmail.com>:
>
> Hi Gaters,
>
> Please, I'm using a phoswich detector in a cylindrical PET scanner. and I need to ask you a question, is it forbidden to use depth value = 5 in the coincidence sorter :
> /gate/digitizer/Coincidences/setDepth 5
>
> No it is not forbidden but if you have defined only 3 depths in your cylindricalPET (for example : the cylindricalPET, rsectors and crystals), depth number 5 does not exist. I am not sure of that and I can't verify it as I am not at work ;-) But if my memory is good, I already had this problem.
>
> If the problem persists, send your macro.
> Cheers,
> Simon
>
> because when i choose this value, the simulation doesn't finish, it stops and gives the following message :
>
>
> *** Break *** floating point exception
> Generating stack trace...
> 0x082625f8 in GateCoincidenceSorter::ComputeSectorID(GatePulse const&) + 0x188 from /usr/local/gate/bin/Linux-g++/Gate
> 0x0826354f in GateCoincidenceSorter::IsForbiddenCoincidence(GatePulse const&, GatePulse const&) + 0x15f from /usr/local/gate/bin/Linux-g++/Gate
> 0x08264bcc in GateCoincidenceSorter::CoincidentPulseIsValid(GateCoincidencePulse*, bool) + 0x11c from /usr/local/gate/bin/Linux-g++/Gate
> 0x082652e2 in GateCoincidenceSorter::ProcessSinglePulseList(GatePulseList*) + 0x642 from /usr/local/gate/bin/Linux-g++/Gate
> 0x081d2e4d in GateDigitizer::Digitize() + 0x1ad from /usr/local/gate/bin/Linux-g++/Gate
> 0x081ec1e8 in GateOutputMgr::RecordEndOfEvent(G4Event const*) + 0x58 from /usr/local/gate/bin/Linux-g++/Gate
> 0x081df89a in GateEventAction::EndOfEventAction(G4Event const*) + 0x2a from /usr/local/gate/bin/Linux-g++/Gate
> 0x087cc6f3 in G4EventManager::DoProcessing(G4Event*) + 0x5c3 from /usr/local/gate/bin/Linux-g++/Gate
> 0x086d493a in G4RunManager::DoEventLoop(int, char const*, int) + 0xba from /usr/local/gate/bin/Linux-g++/Gate
> 0x086d1f8d in G4RunManager::BeamOn(int, char const*, int) + 0x5d from /usr/local/gate/bin/Linux-g++/Gate
> 0x081c93bf in GateApplicationMgr::StartDAQ() + 0xdf from /usr/local/gate/bin/Linux-g++/Gate
> 0x081c9ac2 in GateApplicationMgrMessenger::SetNewValue(G4UIcommand*, G4String) + 0x122 from /usr/local/gate/bin/Linux-g++/Gate
> 0x0895005f in G4UIcommand::DoIt(G4String) + 0xa3f from /usr/local/gate/bin/Linux-g++/Gate
> 0x08956e57 in G4UImanager::ApplyCommand(char const*) + 0x4f7 from /usr/local/gate/bin/Linux-g++/Gate
> 0x08963601 in G4UIbatch::SessionStart() + 0x1e1 from /usr/local/gate/bin/Linux-g++/Gate
> 0x089557f1 in G4UImanager::ExecuteMacroFile(char const*) + 0x41 from /usr/local/gate/bin/Linux-g++/Gate
> 0x08964916 in G4UIcontrolMessenger::SetNewValue(G4UIcommand*, G4String) + 0xe6 from /usr/local/gate/bin/Linux-g++/Gate
> 0x0895005f in G4UIcommand::DoIt(G4String) + 0xa3f from /usr/local/gate/bin/Linux-g++/Gate
> 0x08956e57 in G4UImanager::ApplyCommand(char const*) + 0x4f7 from /usr/local/gate/bin/Linux-g++/Gate
> 0x08963601 in G4UIbatch::SessionStart() + 0x1e1 from /usr/local/gate/bin/Linux-g++/Gate
> 0x089557f1 in G4UImanager::ExecuteMacroFile(char const*) + 0x41 from /usr/local/gate/bin/Linux-g++/Gate
> 0x0822fbb1 in GateUIcontrolMessenger::LaunchMacroFile(G4String) + 0x91 from /usr/local/gate/bin/Linux-g++/Gate
> 0x0822fd08 in GateUIcontrolMessenger::SetNewValue(G4UIcommand*, G4String) + 0x98 from /usr/local/gate/bin/Linux-g++/Gate
> 0x0895005f in G4UIcommand::DoIt(G4String) + 0xa3f from /usr/local/gate/bin/Linux-g++/Gate
> 0x08956e57 in G4UImanager::ApplyCommand(char const*) + 0x4f7 from /usr/local/gate/bin/Linux-g++/Gate
> 0x081c7e0b in main + 0x34b from /usr/local/gate/bin/Linux-g++/Gate
> 0x0663df70 in __libc_start_main + 0xe0 from /lib/libc.so.6
> 0x081c6891 in TCanvasImp::Streamer(TBuffer&) + 0x35 from /usr/local/gate/bin/Linux-g++/Gate
> /usr/bin/Gate: line 7: 4778 Abandon $GATEDIR/bin/Linux-g++/Gate $1
>
>
> Do you know the explanation please ?
>
> Thank you very much
> Ashley
> _________________________________________________________________
> Viens lire ce que Père Noël est en train de faire! Pour les plus récentes nouvelles du Pôle Nord, visite demandeauperenoel.spaces.live.com
> http://demandeauperenoel.spaces.live.com/_______________________________________________
> Gate-users mailing list
> Gate-users at lists.healthgrid.org
> http://lists.healthgrid.org/mailman/listinfo/gate-users
_________________________________________________________________
Découvrez de nouvelles façons de rester en contact grâce à Windows Live! Visitez la Cité @ Live dès aujourd’hui!
http://www.tonadresselive.ca/?icid=LIVEIDFRCA006
More information about the Gate-users
mailing list