
<html>

  <head>


    <meta http-equiv="content-type" content="text/html; charset=UTF-8">

  </head>

  <body>

    Dear Castor Developers,<br>

    <br>

    first of all, thanks for making Castor open source and maintaining

    the code in the first place. It has already helped me<br>

    to find Bugs in my own reconstruction when I was comparing my

    results to the castor result.<br>

    <br>

    Well I guess I get right to the point then. <br>

    I use very simple GATE data to make PET reconstructions and compare

    them to my own reconstruction.<br>

    If I understand multiray correctly, then I think in

    iProjectorIncrementalSiddonMulti::ProjectWithoutTOF() the line 282:<br>

    <br>

    <i> if( alphaMax <= alphaMin ) return 0;</i><br>

    <br>

    will cause the function to return prematurely without having

    iterated over all lines (I mean multiray lines, not oProjectionLine

    objects).<br>

    This can cause for example low sensitivity values on the edges (I

    added a screenshot to the attachment). <br>

    When for example you use 16 lines, but because of the<br>

    random position of each line, line no. 2 might not hit a voxel, the

    function will return and stop iterating over <br>

    the remaining 14 lines.<br>

    <br>

    Or am I missing something? If this is the case, I apologize in

    advance.<br>

    <br>

    <br>

    Best regards,<br>

    Andreas<br>

  </body>

</html>