-
Notifications
You must be signed in to change notification settings - Fork 272
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Change DL1 charge calibration in event to be a single, relative coefficient per channel / pixel #1397
Comments
So this came from the equation I believe is used in HESS, but with the naming changed: |
I believe they are kept in separate values as they are updated on different timescales |
After implementing #1463 and discussing with @watsonjj, we had the following idea for a simplification: The This would contain the total necesseray corrections for an individual event, which would make the code in The complexity would then be shifted where it actually belongs: filling the correct value for an event from different monitoring data. Application is simple, calculating the correct calibration value is hard. So the complexity should be where it is conceptually. |
Since in the CTAO data model, R1 data will already contain the absolute calibration and we will only need to apply a relative update computed from interleaved calibration events, we should really implement the above conclusion. |
Note that this is not 100% correct: in R1 calibration there will be a rough DC/PE factor applied, but it doesn't mean we will not correct that factor in the offline analysis with better calibration information. The main reason these are separated is that they are really different things: one is a conversion factor from DC to PE, so the units are important there. It uses no inter-pixel information, only the position of the single-PE peak in the charge distribution. The flatfielding factor is something else: it corrects for residual differences between pixels, just to make doing things like Hillas parameterization easier (at the expense of losing some information about the pixel's gain factor). It should be a numer very close to 1.0 if the pixels are well-matched in the camera. However, there are reconstruction techniques that use only the PE/DC factor and need to know its value (e.g. Model++), and need to be able to "un-correct" the flatfielding. So I would disagree with this issue partially (I agree the naming and one being multiplicative and the other a divisor is bad): |
Similar comment to my comment in the PR: this issue is about changing the field in the event data structure, not the monitoring data. I think it is good to move the complexity of combining the different monitored quantities like gain and other relative pixel differences to the monitoring code that computes the correction for the current event and only have the result of that computation in the event data structure. |
The documentation of these factors makes no mention that in fact, that the
absolute_factor
is used not as factor but as divisor.they are applied as :
But the documentation of these values never mentions that and I think the naming is very confusing.
Also, why do we have these two factors in the same step? Why are these different calibration constants?
The text was updated successfully, but these errors were encountered: