<div dir="ltr">Ilya,<div><br></div><div>We just learned that the combination of lead glass and lucite contributions to the light yield is a particularly simple one, with the coefficients of the lucite part being zero.</div><div><br><div>-Richard</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 13, 2019 at 1:41 PM Sean Dobbs <<a href="mailto:sdobbs@fsu.edu">sdobbs@fsu.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Richard,<br>
<br>
One correction: by default the energy deposited in the light guides<br>
is not added to the main block energy in mcsmear, because its effect<br>
has not yet been studied in detail.<br>
<br>
---Sean<br>
<br>
On Wed, Mar 13, 2019 at 1:38 PM Richard Jones <<a href="mailto:richard.t.jones@uconn.edu" target="_blank">richard.t.jones@uconn.edu</a>> wrote:<br>
><br>
> Hello Ilya,<br>
><br>
> I agree, although I want to add the light is generated separately in the lead glass blocks and the lucite light guides, in the MC simulation. These two are combined into a single light output in mcsmear, in a treatment that is the same for both hdgeant and hdgeant4.<br>
><br>
> -Richard J.<br>
><br>
> On Wed, Mar 13, 2019 at 12:58 PM Ilya Larin <<a href="mailto:ilarin@jlab.org" target="_blank">ilarin@jlab.org</a>> wrote:<br>
>><br>
>><br>
>> Hi Liping,<br>
>><br>
>><br>
>> I guess final linearity corrections might be involved in<br>
>><br>
>> mcsmear and reconstruction to explain the observed<br>
>><br>
>> positive deviation in reconstructed energy at high energy end.<br>
>><br>
>> I don't think it is an effect of light guides or limited<br>
>> transparency of the blocks.<br>
>><br>
>> Regards,<br>
>> Ilya<br>
>><br>
>> ________________________________<br>
>> От: Gan, Liping <<a href="mailto:ganl@uncw.edu" target="_blank">ganl@uncw.edu</a>><br>
>> Отправлено: 13 марта 2019 г. 12:50<br>
>> Кому: Richard Jones; Ilya Larin<br>
>> Копия: Liping Gan; Matthew Shepherd; <a href="mailto:halld-cal@jlab.org" target="_blank">halld-cal@jlab.org</a><br>
>> Тема: Re: updated study of Fcal response in MC<br>
>><br>
>><br>
>> Hi, Richard,<br>
>><br>
>><br>
>> In your result, the energy deviation for the reconstructed photon in FCAL reaches more than 100 MeV. I would like to ask Ilya to comment on this since he has done a lot of simulations on calorimeter.<br>
>><br>
>><br>
>> Liping<br>
>><br>
>><br>
>><br>
>> ________________________________<br>
>> From: Richard Jones <<a href="mailto:richard.t.jones@uconn.edu" target="_blank">richard.t.jones@uconn.edu</a>><br>
>> Sent: Tuesday, March 12, 2019 2:16 PM<br>
>> To: Gan, Liping<br>
>> Cc: Liping Gan; Matthew Shepherd; <a href="mailto:halld-cal@jlab.org" target="_blank">halld-cal@jlab.org</a><br>
>> Subject: Re: updated study of Fcal response in MC<br>
>><br>
>> Hello Liping,<br>
>><br>
>> The response of the calorimeter itself in the MC does not have these two inflection points that you see in these plots. Those are coming from the interaction between the non-linear response in the simulation coupled with the (non-equivalent) non-linear corrections being applied in mcsmear and the cluster energy reconstruction. This plot shows the net effect of all 3 (MC + mcsmear + recon).<br>
>><br>
>> -Richard<br>
>><br>
>> On Tue, Mar 12, 2019 at 1:13 PM Gan, Liping <<a href="mailto:ganl@uncw.edu" target="_blank">ganl@uncw.edu</a>> wrote:<br>
>><br>
>> Hi, Richard,<br>
>><br>
>><br>
>> Thanks for the excellent result. I just wonder if there is any explanation for the small bump at 1.5-2.5 GeV in you plots? Thanks!<br>
>><br>
>><br>
>> Liping<br>
>><br>
>> ________________________________<br>
>> From: Richard Jones <<a href="mailto:richard.t.jones@uconn.edu" target="_blank">richard.t.jones@uconn.edu</a>><br>
>> Sent: Tuesday, March 12, 2019 10:37:29 AM<br>
>> To: Liping Gan; Matthew Shepherd<br>
>> Cc: <a href="mailto:halld-cal@jlab.org" target="_blank">halld-cal@jlab.org</a><br>
>> Subject: updated study of Fcal response in MC<br>
>><br>
>> Hello Liping and all,<br>
>><br>
>> At your suggestion, I went back and redid my study of reconstructed showers in the FCAL from MC, comparing hdgeant and hdgeant4, this time covering a larger energy range up to 8 GeV. For the results, see the latest document uploaded to docdb.<br>
>><br>
>> <a href="https://gcc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fhalldweb.jlab.org%2Fdoc-private%2FDocDB%2FShowDocument%3Fdocid%3D3852&data=02%7C01%7Challd-cal%40jlab.org%7C8e06c29b6cc34350d8c208d6a7db8e2e%7Cb4d7ee1f4fb34f0690372b5b522042ab%7C1%7C0%7C636880958759491455&sdata=GJyGwJFUTqdMjooogDPFVqWlUtOHp6q2UWHY7d4Nvcs%3D&reserved=0" originalSrc="https://halldweb.jlab.org/doc-private/DocDB/ShowDocument?docid=3852" shash="GyCgaFlU8olQpVsbvFD44b27MBrR1kfNk25kzUeEiWJbsPK0AgO/BfbPZhRf07PHKdqJm4HmT2jRCtXHwaz9sbu1Ypbrw4lE26WKUXf7aq/wTm4EjIp2bk2x8BjKhFO/amevhttPktx9+IHmP87dRnn2GS8cOdJgX5wSq4Q1efo=" rel="noreferrer" target="_blank">https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fhalldweb.jlab.org%2Fdoc-private%2FDocDB%2FShowDocument%3Fdocid%3D3852&data=02%7C01%7Crichard.t.jones%40uconn.edu%7C6ae138e9de9e431492a908d6a7db2443%7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C636880956985861736&sdata=T3HhVQgaeIAC44ubTXVsgB3JXW7WYuUUj9a3NZfxN3M%3D&reserved=0</a><br>
>><br>
>> Here is a summary of what I found.<br>
>><br>
>> I found a small bug in my fitting code that shifted the Ereconstructed down by 16 MeV relative to the generated in my earlier studies. This made the negative y-intercept in the Ereconstructed vs Egenerated appear worse than it actually is.<br>
>> I adjusted the output from shower particles in hdgeant4 relative to hdgeant down by 0.8% as we discussed a the last calorimeter meeting. This change is visible in the new results. I now claim that the agreement in the shower energy response between hdgeant and hdgeant4 is very good.<br>
>> With the new larger simulated energy range, my statistics of 100K events total is limiting my statistical accuracy in studying the non-linearity. Since there is no visible variation with polar angle seen in these individual plots, I added at the end a plot with the full statistics over the range 3-11 degrees, which has less statistical scatter and reveals the nonlinear behavior better.<br>
>> Above 4 GeV there is a very large rise in the reconstructed energy. Overall the nonlinear correction needs work at both the low and high energy ends of the Fcal spectrum.<br>
>><br>
>> -Richard Jones<br>
><br>
> _______________________________________________<br>
> Halld-cal mailing list<br>
> <a href="mailto:Halld-cal@jlab.org" target="_blank">Halld-cal@jlab.org</a><br>
> <a href="https://gcc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.proofpoint.com%2Fv2%2Furl%3Fu%3Dhttps-3A__mailman.jlab.org_mailman_listinfo_halld-2Dcal%26d%3DDwICAg%26c%3DHPMtquzZjKY31rtkyGRFnQ%26r%3DbxTPW7N21WY8eJ2MkW85CQ%26m%3DbouDUuc6DyFXdu5NFXB5Z9lKF6BEkByPb1KeQQSg7OQ%26s%3D8SkCF38MWuPDef6OIjS0x2dCrljm_jf7sqkd73-ZY8s%26e&data=02%7C01%7Challd-cal%40jlab.org%7C8e06c29b6cc34350d8c208d6a7db8e2e%7Cb4d7ee1f4fb34f0690372b5b522042ab%7C1%7C0%7C636880958759491455&sdata=XJayioacmWjk3gOduLIRewZ7SspijK54qDxkRM7VvS4%3D&reserved=0" originalSrc="https://urldefense.proofpoint.com/v2/url?u=https-3A__mailman.jlab.org_mailman_listinfo_halld-2Dcal&d=DwICAg&c=HPMtquzZjKY31rtkyGRFnQ&r=bxTPW7N21WY8eJ2MkW85CQ&m=bouDUuc6DyFXdu5NFXB5Z9lKF6BEkByPb1KeQQSg7OQ&s=8SkCF38MWuPDef6OIjS0x2dCrljm_jf7sqkd73-ZY8s&e" shash="KzKKLvGA76Kh2Bn+o8LKhtR8qLK2csLFTQREfysLJ2vdpLuq0NwJaLgOLQsgYjSOYcsW1FMXrfX4uYJl6oeW1ZjU0AhTPi7YSfsaVPJagwxfKLZdCCYSz8DQDJjjLgfwemyXyp+PauNdZrfYhf5AFV7lgpVD96wN9e+CFzXAqGM=" rel="noreferrer" target="_blank">https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.proofpoint.com%2Fv2%2Furl%3Fu%3Dhttps-3A__mailman.jlab.org_mailman_listinfo_halld-2Dcal%26d%3DDwICAg%26c%3DHPMtquzZjKY31rtkyGRFnQ%26r%3DbxTPW7N21WY8eJ2MkW85CQ%26m%3DbouDUuc6DyFXdu5NFXB5Z9lKF6BEkByPb1KeQQSg7OQ%26s%3D8SkCF38MWuPDef6OIjS0x2dCrljm_jf7sqkd73-ZY8s%26e&data=02%7C01%7Crichard.t.jones%40uconn.edu%7C6ae138e9de9e431492a908d6a7db2443%7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C636880956985861736&sdata=HCDjV1mkHcR3qqI1wO%2F%2Fj%2FOCmJuh72uK%2FZNXuZXmhvM%3D&reserved=0=</a><br>
</blockquote></div>