Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

 The following plot shows the correlation between HEX8 and LEX1 energies between 2 and 1000 MeV. Energy depositions from all 3072 crystal end are plotted on the same histogram. Red points represent the raw adc values multiplyed by linear calibration coefficients (used normally by onboard software) i.e without any correction of nonlinearity. Blue points represents reconstructed energies at each crystal end, i.e. using nonlinearuty correction by the function measured by charge injection calibration.

 
One can see that non-linearity correction is up to 10 % and definitely inproves the situation. Without correction LEX1 energy would be ~10% bigger than HEX8 energy at the end of the range.
Unfortunately, due to significant noise in HEX8 range (~6 MeV sigma) we cannot use this method to study nonlinearity below 3*sig ~ 20 MeV, as profile will be biased by LAC threshold. But we can make profile plot of HEX8/LEX1 ratio for energies between 20 30 MeV and 1000 MeV in log scale. Please note that this plot shows inversed quantity with respect to Philippe's plot of residual nonlinearity, becasue I use HEX8/LEX1 ratio, so if LEX1 is overestimated, then the ratio on my plot would be less than 1.
 
  Image Added
 
Red points shows the HEX8/LEX1 ratio without nonlinearity correction, blue points - after standard nonlinearity correction with charge injection spline used actualy in reconstruction.
 
One can see that after nonlinearity correction the nonlinearity (difference in ratio between 30 MeV and 1000 MeV decreases from 12% to 2% and changes sign, so we've overcorrected the real nonlinearity.

The residual nonlinearity, representad by blue points have the shape similar to Philippe's curve for CU (change of slope at ~200 MeV), but smaller amplitude and opposite direction 

So, I confirm that the effect similar to one seen by Philippe in CU, exist in LAT, but it has different amplitude and direction, so it should be measured and corrected independently from beam test data.

But for calibration unit the correction proposed by Philippe could be applied and should be implemented in reconstruction.