<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none"><!--P{margin-top:0;margin-bottom:0;} @font-face
{font-family:"Cambria Math"}
@font-face
{font-family:Calibri}
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif}
a:link, span.MsoHyperlink
{color:#0563C1;
text-decoration:underline}
a:visited, span.MsoHyperlinkFollowed
{color:#954F72;
text-decoration:underline}
span.EmailStyle17
{font-family:"Calibri",sans-serif;
color:windowtext}
.MsoChpDefault
{font-family:"Calibri",sans-serif}
@page WordSection1
{margin:1.0in 1.0in 1.0in 1.0in}
div.WordSection1
{}--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Hi Sadek,</p>
<p><br>
</p>
<p>If you look at an MLEM equation, it does not make sense to include normalization factor in the system matrix in the right part of the equation but not in the left part (where the sensitivity is taken into account).</p>
<p><br>
</p>
<p>When you perform a list-mode reconstruction, a sensitivity image has to be computed from all possible lines of response that may contribute to the list-mode file.</p>
<p><br>
</p>
<p>So there is a sensitivity image; either you provided it, or it is computed by looping on all crystal pairs or on events from a normalization datafile if provided.</p>
<p><br>
</p>
<p>If you include normalization correction factors in the list-mode file, it means that you are able to know the normalization correction factor for any possible crystal pair admissible in the list-mode file.</p>
<p><br>
</p>
<p>From these normalization factors, you have to build a normalization datafile that contains the list of any admissible line of response with the associated normalization factor (this type of datafile is described in the documentation).</p>
<p><br>
</p>
<p>Then you provide this normalization datafile using '-norm' option to castor, so that it will loop on it to compute the sensitivity image.<br>
</p>
<p><br>
</p>
<p>And the error message will disappear.</p>
<p><br>
</p>
<p>Simon<br>
</p>
<p><br>
</p>
<div style="color: rgb(33, 33, 33);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>De :</b> Castor-users <castor-users-bounces@lists.castor-project.org> de la part de Sadek Nehmeh <san2028@med.cornell.edu><br>
<b>Envoyé :</b> jeudi 22 octobre 2020 05:09<br>
<b>À :</b> castor-users@lists.castor-project.org<br>
<b>Objet :</b> [Castor-users] Normalization</font>
<div> </div>
</div>
<div>
<div class="WordSection1">
<p class="MsoNormal">Hi All,</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I am trying to reconstruct a dataset with the normalization correction factors included in the listmode data. I am getting the following error when I submit for reconstruction:</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">***** The following error message was issued by MPI instance 0</p>
<p class="MsoNormal">***** oSensitivityGenerator::InitializeNormalizationFiles() -> Normalization correction is included in the data file while it is not in the sensitivity computation !</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Here is the reconstruction command:</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">castor-recon -df LM_df.Cdh -opti MLEM -it 2:4 -proj joseph -dim 256,256,309 -vox 2.73,2.73,1.6 -dout my_images</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Please notice that I don’t have a sensitivity image.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Any help would be very appreciated.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Regards,</p>
<p class="MsoNormal">Sadek</p>
</div>
</div>
</div>
</body>
</html>