[Halld-cal] updated study of Fcal response in MC

[Richard Jones]

Hello Liping,

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).

-Richard

On Tue, Mar 12, 2019 at 1:13 PM Gan, Liping <ganl at uncw.edu> wrote:

> Hi, Richard,
>
>
> 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!
>
>
> Liping
> ------------------------------
> *From:* Richard Jones <richard.t.jones at uconn.edu>
> *Sent:* Tuesday, March 12, 2019 10:37:29 AM
> *To:* Liping Gan; Matthew Shepherd
> *Cc:* halld-cal at jlab.org
> *Subject:* updated study of Fcal response in MC
>
> Hello Liping and all,
>
> 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.
>
> https://halldweb.jlab.org/doc-private/DocDB/ShowDocument?docid=3852
> <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%7C20629a2299a9413396c708d6a716d9ba%7Cb4d7ee1f4fb34f0690372b5b522042ab%7C1%7C0%7C636880113921184968&sdata=HIcji%2FTj7Ml35qbyScD%2B7kRj%2FlkwyQcsMC6cdq01bBU%3D&reserved=0>
>
> Here is a summary of what I found.
>
>    1. 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.
>    2. 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.
>    3. 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.
>    4. 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.
>
> -Richard Jones
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.jlab.org/pipermail/halld-cal/attachments/20190312/dc5c67d0/attachment-0001.html>