<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"MS Gothic";
panose-1:2 11 6 9 7 2 5 8 2 4;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
@font-face
{font-family:"\@MS Gothic";
panose-1:2 11 6 9 7 2 5 8 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0064C8;
text-decoration:none none;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#0064C8;
text-decoration:none none;}
p
{mso-style-priority:99;
margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:Consolas;
mso-fareast-language:EN-GB;}
span.EmailStyle22
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:black;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="#0064C8" vlink="#0064C8">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US">Dear mas,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US">Yes, I agree with your statements. I’m glad it helped!<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US">Best regards,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US">David<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US">PS: I think it is possible to define bins of arbitrary size for the np.histogram2d function as follows:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US">xbins=[-250, -0.001, 0.001, 250]<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US">ybins=[-250, -0.001, 0.001, 250]<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US">That would give you a small bin around the center and a few larger ones around it. Of course, it depends on what you want
to do.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif"> mas <fight1_fight2@yahoo.co.jp>
<br>
<b>Sent:</b> Tuesday, 6 September 2022 9:58 AM<br>
<b>To:</b> David Leibold <D.Leibold@tudelft.nl><br>
<b>Cc:</b> gate-users@lists.opengatecollaboration.org<br>
<b>Subject:</b> Re: RE: RE: Re: [Gate-users] How to measure the X-rays scattered by the phantom when they enter the detector?<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p>Dear David,<o:p></o:p></p>
<p><o:p> </o:p></p>
<p>Thank you for your very detailed explanation.<o:p></o:p></p>
<p>I followed your instructions and was able to get the results I wanted.<o:p></o:p></p>
<p>By defining the bin width very small, you mean that there is only a primary x-ray there, right?<o:p></o:p></p>
<p>And when the bin width is quite coarse , the coordinates of (0,0) include primary and some scatter rays (Compton scattering).<o:p></o:p></p>
<p>Then, by subtracting the direct ray count obtained from the (0, 0) coordinates, it is possible to calculate the scattered ray count at (0, 0) when the bin width is coarse.<o:p></o:p></p>
<p><o:p> </o:p></p>
<p>Thank you very much!<o:p></o:p></p>
<p>mas<o:p></o:p></p>
<p><o:p> </o:p></p>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">----- Original Message -----<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><o:p> </o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>From: </strong>"David Leibold" <<a href="mailto:D.Leibold@tudelft.nl">D.Leibold@tudelft.nl</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>To: </strong>"mas" <<a href="mailto:fight1_fight2@yahoo.co.jp">fight1_fight2@yahoo.co.jp</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Cc: </strong>"<a href="mailto:gate-users@lists.opengatecollaboration.org">gate-users@lists.opengatecollaboration.org</a>" <<a href="mailto:gate-users@lists.opengatecollaboration.org">gate-users@lists.opengatecollaboration.org</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Date: </strong>2022/09/05 <span style="font-family:"MS Gothic"">
月</span> 18:00<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Subject: </strong>RE: RE: Re: [Gate-users] How to measure the X-rays scattered by the phantom when they enter the detector?<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><o:p> </o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><o:p> </o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">Dear mas,</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">If we talk about discretized pixels/bins, then I agree: the bin size of the histogram that I chose in the example is quite coarse (depending, of course, on the source-detector-distance), so I agree that it might
happen that scattered rays still end up in the centre bin.</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">Now let’s consider the actual intersection point of an X-ray with the detector plane, so before binning:</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">- In case of Compton scattering, the conservation of impulse dictates that the scattered photon ends up outside (0,0).
</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">- In case of Rayleigh scattering, it is (I think) possible that the scattered photon ends up in (0,0) as well, but since Rayleigh scattering conserves the energy of the photon, a “Rayleigh scattered photon
without change in direction” is indistinguishable from a “primary photon”. (By the way, Rayleigh scatter can usually be neglected in the energy range of diagnostic X-ray energies.)</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">Coming back to the problem that a coarse pixel grid might include scatter in the centre pixel: You can then either choose a very fine pixel grid, or define your bins such that there is one very small bin around
the centre (0,0) and then a large bin for the rest.</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">To get a feeling for what is a good pixel/bin size, you can for example plot a histogram of the distance between the intersection point of a photon with the detector and the origin. Choose a very fine grid for
the histogram. If you histogram bin size is small enough, you will then see that there is a small area around the origin without any counts in your histogram.</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">I hope this answers your questions, if not, let me know!</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">Best regards,</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">David</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>From:</strong> mas <<a href="mailto:fight1_fight2@yahoo.co.jp">fight1_fight2@yahoo.co.jp</a>>
<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Sent:</strong> Monday, 5 September 2022 9:45 AM<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>To:</strong> David Leibold <<a href="mailto:D.Leibold@tudelft.nl">D.Leibold@tudelft.nl</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Cc:</strong> <a href="mailto:gate-users@lists.opengatecollaboration.org">
gate-users@lists.opengatecollaboration.org</a><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Subject:</strong> Re: RE: Re: [Gate-users] How to measure the X-rays scattered by the phantom when they enter the detector?<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">Dear David,<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">Thank you for your very detailed explanation.<span style="font-family:"MS Gothic""> </span>And thank you for teaching me how to open root in python.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">This was very helpful!<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">I tried as you told me. Then I extracted the primary and scatter x-ray data from ”phaseSpace.root".<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">However, one question remained.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">It was assumed that the point hit by the pencil beam (eg P(0,0)) is only primary x-ray.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">However, in actually, I think that scattered rays with small scattering angles also enter (0,0).<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">As a result, the point in (x,y)=(0,0) will sum of the primary and scatter. (eg P(0,0)+S(0,0))<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">In such a case, does it mean that primary and scatter cannot be separated in the phase space actor?<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">I also want to get the scatter component at coordinates (0,0).<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">I'm sorry if my thinking is wrong.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">Thank you<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">mas<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">----- Original Message -----<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>From: </strong>"David Leibold" <<a href="mailto:D.Leibold@tudelft.nl" target="_blank">D.Leibold@tudelft.nl</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>To: </strong>"mas" <<a href="mailto:fight1_fight2@yahoo.co.jp" target="_blank">fight1_fight2@yahoo.co.jp</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Cc: </strong>"<a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank">gate-users@lists.opengatecollaboration.org</a>" <<a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank">gate-users@lists.opengatecollaboration.org</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Date: </strong>2022/09/02 <span style="font-family:"MS Gothic"">
金</span> 16:07<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Subject: </strong>RE: Re: [Gate-users] How to measure the X-rays scattered by the phantom when they enter the detector?<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">Dear Mas,</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">You mention that you got a one-dimensional profile, and you show a screenshot of the root output. However, please note that you select only the x-coordinates of the root output, that means, you select the x-coordinate
of all particles intersecting your detector plane. Since you select only one variable, you will obtain a one-dimensional output.
</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">The phase space actor does not have a built-in function that separates scatter and primaries; you have to do this yourself.</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">What you have to do is the following (and I hope that also answers your questions 2-3 as well): You have to extract the (x,y) coordinates for all events, and then filter the events based on their distance to the
centre of the detector. (I assume that the centre of the detector is the point at which the pencil beam is directed.) To be honest, I have never really used the root object browser, so others might share more light on this. For my data analysis I use Python.
There is a Python library called “</span><a href="https://urldefense.com/v3/__https:/uproot.readthedocs.io/en/latest/__;!!PAKc-5URQlI!4TACOtbbEWaSmINTkHUyQlT5ut0PBTLOoLSOiu362wQGtm8P0OPqzNqcB5JU8UtZY2ovvrMSeogoJOQDJq8PAEzFh4tM0w$" target="_blank"><span style="color:black">uproot</span></a><span style="color:black">”,
which enables you to read the root files.</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">Here's some example code to read the root file that is output by the phase space actor; in this example the file is called “phaseSpace.root”.</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">import numpy as np<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">import uproot<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">import matplotlib.pyplot as plt<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"># Load the root file:<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">with uproot.open("phaseSpace.root") as rootfile:<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> data = rootfile['PhaseSpace'].arrays(filter_name="*", library='np’)<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # root files have a tree structure, like a trunk with many branches.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # The “trunk” of a root file that is output by a phase space actor<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # is called ‘PhaseSpace’. With the function ‘arrays’ we can<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # access all the branches that we want. You can specify the<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # desired branches with the argument ‘filter_name’. If you use the<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # wildcard ‘*’, then all branches are loaded. If you just want the<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # x and y coordinate, you can do<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # filter_name = [‘X’, ‘Y’]<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # . The argument ‘library’ specifies whether you want to use numpy<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # or pandas. The former is better in terms of performance, the<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # latter is nicer for visualisation.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"># Create a 2D histogram:<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">xbins = np.arange(-240, 241, 1)<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">ybins = np.arange(-240, 241, 1)<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">histData, histXbins, histYbins = np.histogram2d(data['X'], data['Y'], bins=[xbins, ybins])<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"># Plot the result:<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">fig, ax = plt.subplots()<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">ax.imshow(histData, vmin=0, vmax=2)<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # Please note that the origin of the plot does not coincide with the<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> # data origin, you would have to shift that manually<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><img border="0" width="675" height="472" style="width:7.0277in;height:4.9166in" id="_x0000_i1025" src="cid:image001.png@01D8C1DF.C4BC6B60"><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">(That little dot in the centre of the figure is not dust on your screen, it’s actually a very faint signal…)</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">Based the 2D histogram, you can now say that all events that are outside, for example, a 0.1 mm region around the origin, are scatter.</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">I hope that helps! If you need more help, let me know.</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">Best regards,</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">David</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black">PS: Sorry for writing your name incorrectly in my last mail!</span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><span style="color:black"> </span><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>From:</strong> mas <<a href="mailto:fight1_fight2@yahoo.co.jp" target="_blank">fight1_fight2@yahoo.co.jp</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Sent:</strong> Friday, 2 September 2022 8:27 AM<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>To:</strong> David Leibold <<a href="mailto:D.Leibold@tudelft.nl" target="_blank">D.Leibold@tudelft.nl</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Cc:</strong> <a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank">
gate-users@lists.opengatecollaboration.org</a><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Subject:</strong> Re: Re: [Gate-users] How to measure the X-rays scattered by the phantom when they enter the detector?<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">Dear David,<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">Thanks to your reply, I understood that I should use the phase space actor.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">After reading the documentation, I tried using phase space actor.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">I got a one-dimensional profile.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">It was certainly like a point spread function.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">However, some questions remained.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">1)The output is a one-dimensional profile, but how do I know which location on the detector is profiled?<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">2)Which function of the phase space actor is used to separate scattered and direct lines?<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">[/gate/actor/MyActor/storeOutgoingParticles true] is used?<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">3)How do I output the actual values from the root file, e.g. a txt file? I am sorry that this is not the original question.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">Thank you.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">mas<o:p></o:p></p>
</blockquote>
<pre> <o:p></o:p></pre>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">----- Original Message -----<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>From: </strong>"David Leibold" <<a href="mailto:D.Leibold@tudelft.nl" target="_blank">D.Leibold@tudelft.nl</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>To: </strong>"mas" <<a href="mailto:fight1_fight2@yahoo.co.jp" target="_blank">fight1_fight2@yahoo.co.jp</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Cc: </strong>"<a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank">gate-users@lists.opengatecollaboration.org</a>" <<a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank">gate-users@lists.opengatecollaboration.org</a>><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Date: </strong>2022/09/02 <span style="font-family:"MS Gothic"">
金</span> 00:35<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><strong>Subject: </strong>Re: [Gate-users] How to measure the X-rays scattered by the phantom when they enter the detector?<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">Dear Masa,<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">May I summarise your requirements:<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">1) You mention that you want to simulate a CsI detector, but with 100 % detection efficiency and no interaction within the detector. <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">2) You would like to distinguish between scattered and primary rays.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">3) A monochromatic pencil beam enters a phantom.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">My suggestion:<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">I would suggest using a phase space actor (see the <a href="https://urldefense.com/v3/__https:/opengate.readthedocs.io/en/latest/tools_to_interact_with_the_simulation_actors.html*phase-space-actor__;Iw!!PAKc-5URQlI!_P9nlwtwj9dIeV2hHQtSbBZiOOlf2u0klOOSbhmmNrMhNxMxAJN4SgVVz_zaiEj0SQYvfHgXq2pNEPKPSyXWmZGQZIKttg$" target="_blank">documentation
here</a>), which registers each particle's energy, position, velocity etc. crossing its surface. Since you do not want to simulate interaction in your CsI detector material, there is no need to use it in the first place.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">Since you are using a pencil beam, all primary photons will go straight through your phantom, intersecting your detector at, for example, point (0,0). All scattered radiation will intersect your detector away from point (0,0). Hence, you
can filter out the scatter and primary rays purely based on where they end up on your detector.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">To set up your detector, use e.g. a box volume with the desired size, fill it with vacuum, and attach a phase space actor to it. <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">If you need more help, please let us know.<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">Best regards,<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">David <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">On 30 Aug2022, at 00:46, mas <<a href="mailto:fight1_fight2@yahoo.co.jp" target="_blank">fight1_fight2@yahoo.co.jp</a>> wrote:<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<pre>Hi GATE users.<o:p></o:p></pre>
<pre>I am masa@japan and new to Gate.<o:p></o:p></pre>
<pre> <o:p></o:p></pre>
<pre>I am currently working on a simple simulation.<o:p></o:p></pre>
<pre>I am trying to find the point spread function of the scattered rays detected by a CsI detector by injecting a pencil beam of monochromatic X-rays into an acrylic phantom.<o:p></o:p></pre>
<pre> <o:p></o:p></pre>
<pre>Assuming a detector detection efficiency of 100%, I want to output the spread of scattered radiation when it enters the detector.<o:p></o:p></pre>
<pre>No interaction occurs with the detector.<span style="font-family:"MS Gothic""> </span>I also want to output scattered rays and primary rays separately.<o:p></o:p></pre>
<pre>How should I set up the detector in such a case?<o:p></o:p></pre>
<pre> <o:p></o:p></pre>
<pre>Should I use Fluence Actor for detection?<o:p></o:p></pre>
<pre>Or should I use a sensitive detector?<o:p></o:p></pre>
<pre> <o:p></o:p></pre>
<pre>I would appreciate it if you could help me.<o:p></o:p></pre>
<pre> <o:p></o:p></pre>
<pre>Thank you.<o:p></o:p></pre>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">_______________________________________________<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal">Gate-users mailing list<o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><a href="mailto:Gate-users@lists.opengatecollaboration.org" target="_blank">Gate-users@lists.opengatecollaboration.org</a><o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><a href="https://urldefense.com/v3/__http:/lists.opengatecollaboration.org/mailman/listinfo/gate-users__;!!PAKc-5URQlI!5R6DdcwuD-nXX0bmvN3RckxZ2HEAgen9ExhbGuDcz_52hS9V1Rohx8vBP7HBR0msh9Bx0lNmj5dJ866ftlukHHjC8jk0_A$" target="_blank">https://urldefense.com/v3/__http://lists.opengatecollaboration.org/mailman/listinfo/gate-users__;!!PAKc-5URQlI!5R6DdcwuD-nXX0bmvN3RckxZ2HEAgen9ExhbGuDcz_52hS9V1Rohx8vBP7HBR0msh9Bx0lNmj5dJ866ftlukHHjC8jk0_A$</a> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"> <o:p></o:p></p>
</blockquote>
<blockquote style="border:none;border-left:solid #CCCCCC 3.0pt;padding:0cm 0cm 0cm 5.0pt;margin-left:7.5pt;margin-right:0cm">
<p class="MsoNormal"><o:p> </o:p></p>
</blockquote>
</div>
</body>
</html>