I've refactored the code for rqpy.limit.drde_gauss_smear2d so that the code does a better job at calculating the smeared dR/dE using a 2D normal distribution. I also fixed a bug that was affecting the smearing. The explicit changes are detailed below.
remove rqpy.limit._limit._largest_et function
change calculation of the integral bounds in the trigger energy estimator
make sure that the true energy is integrated in 1 eV bins
change the heaviside theta function so that it is evaluated as 1 at threshold
fix bug in the lambda function defining the double integral so that the true energy is inputted correctly
I've done some testing for different levels of smearing and compared to data, and the code seems to be working well.
I've refactored the code for
rqpy.limit.drde_gauss_smear2dso that the code does a better job at calculating the smeared dR/dE using a 2D normal distribution. I also fixed a bug that was affecting the smearing. The explicit changes are detailed below.rqpy.limit._limit._largest_etfunctionI've done some testing for different levels of smearing and compared to data, and the code seems to be working well.