<div dir="ltr"><div>Hi,</div><div>Please stick to the mailing list.</div><div>Sorry, this code is quite a mess... You can find how this is computed here:</div><div><a href="https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/include/GateFixedForcedDetectionFunctors.hh#L345">https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/include/GateFixedForcedDetectionFunctors.hh#L345</a></div><div>So bin[0] is from 0 to 0.5*energyResolvedBinSize keV, bin[1] from 0.5*energyResolvedBinSize to 1.5*energyResolvedBinSize keV etc.</div><div>Simon<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 29, 2017 at 1:57 PM, Clemens S. <span dir="ltr"><<a href="mailto:clem.schmid@gmail.com" target="_blank">clem.schmid@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Simon,<div><br></div><div>concerning the resolved bin size: I don't understand the output of the energy resolved images exactly. They seem to have one more entry than expected when looking at the maximum spectrum energy divided by energyResolvedBinSize.</div><div><br></div><div>For example: </div><div><font face="monospace, monospace">    max(spectrum) = 141 keV</font></div><div><font face="monospace, monospace">    energyResolvedBinSize = 1 keV</font></div><div><font face="monospace, monospace">        => first image dimension = 142</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="arial, helvetica, sans-serif">Can you tell me which energy corresponds to   </font><font face="monospace, monospace">image[i] </font><font face="arial, helvetica, sans-serif">exactly? I had a look at the FFD source code but I could not find an answer so far.</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Thank you very much for your help!</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Kind regards</font></div><span class="HOEnZb"><font color="#888888"><div><font face="arial, helvetica, sans-serif">Clemens</font></div><div><font face="monospace, monospace"><br></font></div><div><br></div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-10-04 10:34 GMT+02:00 Clemens S. <span dir="ltr"><<a href="mailto:clem.schmid@gmail.com" target="_blank">clem.schmid@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Simon,<div><br></div><div>this sound perfect! Thank you very much.</div><div><br></div><div>Kind regards</div><span class="m_-2350030350452461067HOEnZb"><font color="#888888"><div>Clemens</div></font></span></div><div class="m_-2350030350452461067HOEnZb"><div class="m_-2350030350452461067h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-10-04 10:31 GMT+02:00 Simon Rit <span dir="ltr"><<a href="mailto:simon.rit@creatis.insa-lyon.fr" target="_blank">simon.rit@creatis.insa-lyon.f<wbr>r</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hi,<br></div>What I do is record one spectral projection by setting the energyResolvedBinSize to 1 keV and a detector response d(e)=1 with a text file that looks like this:<br>0 1<br>1 1<br></div>I obtain then a 3 image for each run where the 3rd dimension is the energy. I use this result offline to create the projection of each threshold separately.<br></div>I don't know when this option was added wrt to Gate 7.2, you might have to update your Gate.<br></div>Simon<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_-2350030350452461067m_-410477377591758832h5">On Wed, Oct 4, 2017 at 10:18 AM, Clemens S. <span dir="ltr"><<a href="mailto:clem.schmid@gmail.com" target="_blank">clem.schmid@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_-2350030350452461067m_-410477377591758832h5"><div dir="ltr">Dear fellow GATE users,<div><br></div><div>I am simulating spectral CT with <b>Fixed Forced detection</b> in GATE 7.2. I am simulating one physical detector which has several energy threshold, i.e. <b>several detector responses</b>. Each detector response produces its own output (primary and scatter images).</div><div><br></div><div>Up to now, I'm using one FFD actor for each threshold/response. This gives correct results, but is not optimally efficient since the interaction points are projected onto the detectors indepentently instead of projecting once (conserving energy information) and <b>then</b> applying the detector responses. (or I think that's what's happening)</div><div><br></div><div><b>My question is</b>: Is there a (more) efficient way to simulate several detector responses with FFD?</div><div><br></div><div>Thank you for your help!</div><div><br></div><div>Kind regards</div><div>Clemens Schmid</div></div>
<br></div></div>______________________________<wbr>_________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org" target="_blank">Gate-users@lists.opengatecolla<wbr>boration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank">http://lists.opengatecollabora<wbr>tion.org/mailman/listinfo/gate<wbr>-users</a><br></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>