[Gate-users] Multi-System Digitizer

Wed Aug 5 06:22:57 CEST 2020

Hello,

Having done further investigation into why I am getting this error in Gate
v9.0, I tried running my code in Gate v8.2 and found in that version, the
segmentation violation does not occur. Unfortunately, that doesn't
completely solve the problem as the ring repeater for the second system
does not work the way I expect it to. Both systems have two heads but while
SpectHead1 has two sets of projection data, SpectHead2 only has one set of
full projection data and the second set is blank.
If I change the ring repeater for the first system (SpectHead1) to have 3
heads, then there are three sets of projection data produced for that
system. It also tried to give SpectHead2 three sets of projection data
(even though its ring repeater was set to two repeats) but now two of them
are blank. If I change the ring repeater of SpectHead2, nothing changes in
regard to the projection data.
In all cases, the visualisation shows the number of heads and positions
that I expect. I don't get any errors when I run the simulation in v8.2.

I've included a zip file with the macros I'm using to test the system. Is
this a problem with the way the ring repeater is implemented in v8.2, or a
problem in how I am using it?

Cheers,
Miriam

On Tue, Jul 28, 2020 at 3:49 PM Miriam Dixon <miriam.gatemc at gmail.com>
wrote:

> Hello GATE-users,
>
> I am trying to set up a four-head SPECT system by creating two 2-head
> systems, SpectHead1 and SpectHead2, both using the built in SPECThead
> system. I have found the online documentation outlining how to do this, but
> I am getting an error I am not sure how to fix. The error does not occur if
> I run the simulation only using SpectHead1 or only using SpectHead2, but
> does occur if they are both included. This indicates there is something
> wrong in how they interact within the simulation - does anyone have a
> suggestion of what I could do to fix this?
>
> The error message is:
> *** Break *** segmentation violation
>
> The lines below might hint at the cause of the crash.
>
> #5  0x0000564fdfdac8ed in GateDeadTime::ProcessOnePulse(GatePulse const*,
> GatePulseList&) ()
>
> #6  0x0000564fe01d7594 in
> GateVPulseProcessor::ProcessPulseList(GatePulseList const*) ()
>
> #7  0x0000564fe0138f9d in GatePulseProcessorChain::ProcessPulseList() ()
>
> #8  0x0000564fdfdc64ac in GateDigitizer::DigitizePulses() ()
>
> #9  0x0000564fdfdc6875 in GateDigitizer::Digitize() ()
>
> #10 0x0000564fdff4c545 in GateOutputMgr::RecordEndOfEvent(G4Event const*)
> ()
>
> #11 0x0000564fdfd39a5b in GateEventAction::EndOfEventAction(G4Event
> const*) ()
>
> #12 0x00007fdabbec3be2 in G4EventManager::DoProcessing(G4Event*) () at
> /usr/share/geant4/geant4-install/lib/libG4event.so
>
> #13 0x00007fdabc163585 in G4RunManager::ProcessOneEvent(int) () at
> /usr/share/geant4/geant4-install/lib/libG4run.so
>
> #14 0x00007fdabc16167b in G4RunManager::BeamOn(int, char const*, int) ()
> at /usr/share/geant4/geant4-install/lib/libG4run.so
>
> #15 0x0000564fe01f8bf1 in GateApplicationMgr::StartDAQ() ()
>
> #16 0x0000564fe01ffb50 in
> GateApplicationMgrMessenger::SetNewValue(G4UIcommand*, G4String) ()
>
> #17 0x00007fdab5e3e3bc in G4UIcommand::DoIt(G4String) () at
> /usr/share/geant4/geant4-install/lib/libG4intercoms.so
>
> #18 0x00007fdab5e5fc4c in G4UImanager::ApplyCommand(char const*) () at
> /usr/share/geant4/geant4-install/lib/libG4intercoms.so
>
> #19 0x00007fdab5e29c77 in G4UIbatch::ExecCommand(G4String const&) () at
> /usr/share/geant4/geant4-install/lib/libG4intercoms.so
>
> #20 0x00007fdab5e2bc9e in G4UIbatch::SessionStart() () at
> /usr/share/geant4/geant4-install/lib/libG4intercoms.so
>
> #21 0x00007fdab5e6114c in G4UImanager::ExecuteMacroFile(char const*) () at
> /usr/share/geant4/geant4-install/lib/libG4intercoms.so
>
> #22 0x00007fdab5e4b9cb in G4UIcontrolMessenger::SetNewValue(G4UIcommand*,
> G4String) () at /usr/share/geant4/geant4-install/lib/libG4intercoms.so
>
> #23 0x00007fdab5e3e3bc in G4UIcommand::DoIt(G4String) () at
> /usr/share/geant4/geant4-install/lib/libG4intercoms.so
>
> #24 0x00007fdab5e5fc4c in G4UImanager::ApplyCommand(char const*) () at
> /usr/share/geant4/geant4-install/lib/libG4intercoms.so
>
> #25 0x0000564fdfb48c16 in main ()
>
>
> I have included an outline of how I created the two systems below, if that
> helps identify the problem.
>
>
> Cheers,
>
> Miriam
>
>
> Below is a brief outline of what I have done to create the two systems:
> *## Detector System 1*
> /gate/world/daughters/name SpectHead1
> /gate/world/daughters/systemType    SPECThead
> /gate/world/daughters/insert box
> /gate/SpectHead1/geometry/setXLength 57.6 cm
> /gate/SpectHead1/geometry/setYLength 44.6 cm
> /gate/SpectHead1/geometry/setZLength 18 cm
> /gate/SpectHead1/vis/setColor white
> /gate/SpectHead1/setMaterial Air
> /gate/SpectHead1/vis/forceWireframe
> /gate/SpectHead1/repeaters/insert ring
> /gate/SpectHead1/ring/setPoint1 0 0 0
> /gate/SpectHead1/ring/setPoint2 0 1 0
> /gate/SpectHead1/ring/enableAutoRotation
> /gate/SpectHead1/ring/setFirstAngle 0 deg
> /gate/SpectHead1/ring/setAngularSpan 76 deg
> /gate/SpectHead1/ring/setRepeatNumber 2
>
> *## Detector System 2*
> /gate/world/daughters/name SpectHead2
> /gate/world/daughters/systemType SPECThead
> /gate/world/daughters/insert box
> /gate/SpectHead2/geometry/setXLength 57.6 cm
> /gate/SpectHead2/geometry/setYLength 44.6 cm
> /gate/SpectHead2/geometry/setZLength 18 cm
> /gate/SpectHead2/vis/setColor white
> /gate/SpectHead2/setMaterial Air
> /gate/SpectHead2/vis/forceWireframe
> /gate/SpectHead2/repeaters/insert ring
> /gate/SpectHead2/ring/setPoint1 0 0 0
> /gate/SpectHead2/ring/setPoint2 0 1 0
> /gate/SpectHead2/ring/enableAutoRotation
> /gate/SpectHead2/ring/setFirstAngle 180 deg
> /gate/SpectHead2/ring/setAngularSpan 76 deg
> /gate/SpectHead2/ring/setRepeatNumber 2
>
> *....further details about the head setup (can be provided if needed)....*
>
> *## Digitizer System 1 *
> /gate/digitizer/Singles/insert systemFilter
> /gate/digitizer/Singles/systemFilter/selectSystem SpectHead1
> /gate/digitizer/Singles/insert adder
> /gate/digitizer/Singles/insert readout
> /gate/digitizer/Singles/readout/setDepth 1
> /gate/digitizer/Singles/insert spblurring
> /gate/digitizer/Singles/spblurring/setSpresolution 3.97 mm
> /gate/digitizer/Singles/spblurring/verbose 0
> *##....further details about the digitizer (can be provided if needed)....*
> /gate/digitizer/Singles/insert deadtime
> /gate/digitizer/Singles/deadtime/setDeadTime 1650 ns
> /gate/digitizer/Singles/deadtime/setMode paralysable
> /gate/digitizer/Singles/deadtime/chooseDTVolume crystal_volume
>
> *## Digitizer System 2*
> /gate/digitizer/name    Singles2
> /gate/digitizer/insert  singleChain
> /gate/digitizer/Singles2/insert systemFilter
> /gate/digitizer/Singles2/systemFilter/selectSystem SpectHead2
> /gate/digitizer/Singles2/name adder2
> /gate/digitizer/Singles2/insert adder
> /gate/digitizer/Singles2/name readout2
> /gate/digitizer/Singles2/insert readout
> /gate/digitizer/Singles2/readout2/setDepth 1
> /gate/digitizer/Singles2/name spblurring2
> /gate/digitizer/Singles2/insert spblurring
> /gate/digitizer/Singles2/spblurring2/setSpresolution 3.97 mm
> /gate/digitizer/Singles2/spblurring2/verbose 0
> *##....further details about the digitizer (can be provided if needed)....
> *
> /gate/digitizer/Singles2/name deadtime2
> /gate/digitizer/Singles2/insert deadtime
> /gate/digitizer/Singles2/deadtime2/setDeadTime 1650 ns
> /gate/digitizer/Singles2/deadtime2/setMode paralysable
> /gate/digitizer/Singles2/deadtime2/chooseDTVolume crystal_volume2
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200805/8f70a067/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SPECT4Head.zip
Type: application/x-zip-compressed
Size: 14500 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200805/8f70a067/attachment-0001.bin>