-
Notifications
You must be signed in to change notification settings - Fork 77
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
Wrong base nsb level used in tuning at waveform level #1271
Comments
This reliable way exists in the most recent version of sim_telarray and I already opened a corresponding issue in the lst1-sim-config repo two years ago: cta-observatory/lst-sim-config#29 With the metadata, you always get the final configuration value of the configuration item for each telescope. |
Thanks @maxnoe . Sadly it seems the LSTProd2 was not done with the latest version. We have no |
Looking a bit more, I have found a way to differentiate between the wrong and correct NSB levels in the config history. The correct telescope wise entries of 'NIGHTSKY_BACKGROUND' are preceded by an entry 'STORE_PHOTOELECTRONS' while the wrong ones are preceded by 'MIN_PHOTOELECTRONS' in both prod 5 and LSTprod2. I could thus adapt the implementation of |
Hi @gabemery, the current implementation does work for LSTProd2, right? To make it work for Prod5 too, since the goal is to achieve the same level of fluctuations in MC and data, isn't it possible to obtain empirically the level of fluctuations from the MC waveforms? |
@moralejo Current implementation is wrong for LSTProd2. It was implemented in prod5 (and was correct then). The issue is that the level of fluctuation in MC vs data informs on the multiplicative factor to the simulated NSB required to reach the level of data. So knowing the original pulse rate in MC is mandatory to know the number of pulse that needs to be added. |
btw. there is no need for you to do this by hand, the
|
By the way, @gabemery, does the method assume that all noise is from NSB? It is certainly dominant, but there is also electronic noise. |
@moralejo Yes it does assume it is all NSB. If there is a non-negligible contribution of the electronic noise to the waveform variance it would thus be biased. |
The tuning at waveform level is done by extracting the relative NSB increase required and getting the original simulated NSB from a Header in simtel files.
Original implementation extracted all occurrences of 'NIGHTSKY_BACKGROUND' in yield_all_subobjects(EventIOFile(filename), [History, HistoryConfig]).
When tested on prod5_tran80 files the output was :
[0.25307, 0.24586, 0.25307, 0.25307, 0.25307, 0.1, 0.12]
where index 1 was the value for LST-1, index 2 for LST-2, ... index 5 for MAGIC 1 and index 6 for MAGIC 2.
Index 0 was ignored and assumed to be some unused default value.
It was then correct to use tel_id to access the correct original NSB level.
This is not the case with the current MC production which has an additional value : [0.24586, 0.217, 0.24586, 0.1, 0.12].
Reported by @jstvdk
The preferred way to solve this issue would be to have a reliable way to find the NSB level effectively used for each telescope.
The text was updated successfully, but these errors were encountered: