-
Notifications
You must be signed in to change notification settings - Fork 35
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
Improve debug output direction fitter #549
base: develop
Are you sure you want to change the base?
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
See comment regarding the standard deviation - I don't think this is correct right now. Other than that, this looks good to me. Documentation failing is due to the make_docs script, I'll look into that next week as well. Thanks Christian!
mean_sim_az = np.angle(np.mean(temp_complex_representation.real)+ 1j * np.mean(temp_complex_representation.imag)) | ||
std_sim_az = np.angle(np.std(temp_complex_representation.real)+ 1j * np.std(temp_complex_representation.imag)) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
First line could in principle be simplified (the two means are additive), and I'm not sure the second one is correct - Wikipedia backed up my some toy Monte Carlo suggests
std_sim_az = np.sqrt(np.log(1/np.mean(temp_complex_representation)**2))
try: | ||
import NuRadioReco.eventbrowser.dataprovider_root | ||
except: | ||
print("can't import RNO-G root reader") |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The dataprovider_root has not yet been updated for the new RNO-G reader, and we've already removed the old one, leading to this error. This is already fixed on rnog_eventbrowser
, I can make a PR next week. But it's probably okay to leave it like this.
improves output of correlation DirectionFitter to calcualate angular differences better.
Also fixed hardcoded RNO-G specific import