[gate-users] patch for STIR for recon Gate output
Qianmei Zhang
qzhang at triumf.ca
Fri Jan 7 23:16:27 CET 2005
Hello Claude and Gate Users
I have a question of how to get the Ecat7 sinogram output from certain
block rings in central.
Since i want to explore whether the random scattering events can be
supressed by adding more block rings in the Ecat 953 system (two block
rings). These added block detector rings act like this: if two gamma
rays were detected ( no matter whether it is true or random) by the two
block rings and other gamma rays were detected by the added block rings
within the coincidence time window, the event would not be recorded.
The images reconstruction only use the two block rings(same as Ecat953), not the added
block rings.
In this case, I have to get the Ecat7 sinogram just from the central two
block rings after coincidence sorter in Gate. I tried to open the source code
of GatetoSinogram.cc in line 136
m_ringNb = blockComponent->GetLinearRepeatNumber() * crystalComponent->GetRepeatNumber(2);
and GatetoSinogram.cc in line 284
G4int RingNb = blockComponent->GetLinearRepeatNumber() * crystalComponent->GetRepeatNumber(2);
Can i just change GetLinearRepeatNumber() to 2 and temporary make Gate again ?
or should i change
all setMaker->GetRingNb() to 15 in GatetoSinogram.cc ?
Is there any other parameter i should change? Or is there any other way
to handle this question. Since i am new in Gate and not familiar to C,
any comments are very welcomed. Thanks.
Qianmei
More information about the Gate-users
mailing list