aboutsummaryrefslogtreecommitdiff
AgeCommit message (Collapse)Author
2018-10-03add \ntlatorre
2018-10-03add a README filetlatorre
2018-10-03update .gitignoretlatorre
2018-10-03move python scripts into utils/ directorytlatorre
2018-10-03add a script to plot fit resultstlatorre
2018-10-03add a python script to generate SNOMAN command filestlatorre
2018-10-03make sure wavelengths are uniformly spaced since we're using interp1d()tlatorre
2018-10-02update get_quantum_efficiency() to use interp1d()tlatorre
This speeds up the "boot up" time from ~30 seconds to ~12 seconds.
2018-10-02update MIN_THETA0 to 0.02tlatorre
2018-10-02update charge fraction to 0.5tlatorre
2018-10-02update optics to use LETA constants for the D2O phasetlatorre
This commit updates the PMT response bank and absorption lengths for H2O and D2O based on the values in mcprod which were used for the LETA analysis.
2018-10-01flag PMTs with the KPF_DIS bit settlatorre
Doh! I was previously using the tube number variable!
2018-10-01fall back to old scattering rms calculation when n = 0tlatorre
2018-10-01update negative log likelihood for path coefficientstlatorre
2018-10-01loop over all normal PMTs when calculating the expected number of photonstlatorre
Previously we ignored PMTs which were flagged when computing the expected number of PE for each PMT, but since we calculate the amount of reflected light here we need to include even PMTs which are offline (since they still reflect light).
2018-10-01use the PMT response table to calculate the amount of reflected lighttlatorre
To calculate the expected number of photons from reflected light we now integrate over the track and use the PMT response table to calculate what fraction of the light is reflected. Previously we were just using a constant fraction of the total detected light which was faster since we only had to integrate over the track once, but this should be more accurate.
2018-10-01add absorption length for acrylictlatorre
2018-10-01fix a bug when computing the absorption length in H2O and D2Otlatorre
Previously I was interpolating the absorption lengths using interp1d() but that only works when the x array is uniform. Since the wavelengths are not spaced uniformly, we have to use the GSL interpolation routines.
2018-09-26speed up fast likelihood calculationtlatorre
This commit updates the fast likelihood calculation to use the identity sin(a-b) = sin(a)*cos(b) - cos(a)*sin(b) to speed up the fast likelihood calculation.
2018-09-26speed up fast likelihood calculationtlatorre
This commit speeds up the fast likelihood calculation by avoiding calls to trigonometric functions where possible. Specifically we calculate sin(a) = sqrt(1-pow(cos(a),2)); instead of sin(a) = sin(acos(cos(a)));
2018-09-26increase number of points in cos(theta) interpolation to 1000tlatorre
2018-09-25update indirect scattering PDF start timetlatorre
Currently the PDF for scattered light is modelled as a flat distribution starting at some time t. Previously I was using the mean hit time for all PMTs, however this should really be a flat distribution in the time *residual* after the main peak. Therefore, the PDF now starts at the estimated time for direct photons.
2018-09-25update likelihood calculation to use PMT_TTS macrotlatorre
I accidentally hardcoded the single PE TTS to 1.5 ns in the likelihood calculation.
2018-09-25update integration bounds in likelihood calculationtlatorre
This commit updates the bounds of the track integration in the likelihood function to integrate up to 1 meter around the point at which the PMT is at the Cerenkov angle from the track. This fixes an issue I was seeing where a *very* small change in the fit paramters would cause the likelihood to jump by a large amount. I eventually tracked it down to the same issue I was seeing before which I solved by splitting up the integration into two intervals. However that fix did not seem to completely fix the issue. Based on initial tests with 500 MeV muons, this fix seems to do a much better job.
2018-09-25increase maxeval to 20 for the "quick" minimization phasetlatorre
When testing out the fitter on 500 MeV muons, there was at least one event which started out at a position very far from the true position. This event had a secondary electron like ring which is what I suspect caused the fit to start out in a position far from the true position. This fix correctly starts the minimization close to the true position. In the future I should look at updating get_direction() so that it finds the largest ring direction instead of just doing a weighted average of all the vectors from the position to the PMTs.
2018-09-21update likelihood function to include the probability of the path coefficientstlatorre
2018-09-21split up the track integral into two intervalstlatorre
This commit updates the likelihood calculation to split up the track integral into two intervals in some cases. I noticed when fitting some events that the likelihood value would change drastically for a very small change in the fit parameters. I eventually tracked it down to the fact that the track integral was occasionally returning a very small charge for a PMT which should have a very high charge. This was happening because the region of the track which was hitting the PMT was very small and the cquad integration routine was completely skipping it. The solution to this problem is a bit of a hack, but it seems to work. I first calculate where along the track (for a straight track) the PMT would be at the Cerenkov angle from the track. If this point is somewhere along the track then we split up the integral into two intervals: one going from the start of the track to this point and the other from the point to the end of the track. Since the cquad routine always samples points near the end of the intervals this should prevent it from completely skipping over the point in the track where the integrand is non-zero.
2018-09-21increase MAX_PE to 10000tlatorre
2018-09-20delete unused variable distancetlatorre
2018-09-20don't include the OWL PMTs in the likelihood calculationtlatorre
For some reason the OWL tubes have 9999.00 for the x, y, and z coordinates of the normal vector in the PMT file. For now, I'm just going to remove them from the likelihood calculation.
2018-09-20make sure direction vector is normalized in path_eval()tlatorre
2018-09-20add particle id code to output filetlatorre
2018-09-20add git SHA1 hash to output filetlatorre
2018-09-20add time elapsed to the output filetlatorre
2018-09-20add a command line option to only fit the first eventtlatorre
This commit adds the ability to run the fit program with the --skip-second-event flag to only fit the first event after a MAST bank. This way we avoid fitting secondaries like Michel electrons when fitting MC events.
2018-09-20add absorption lengths for D2O and H2O weighted by the Cerenkov spectrum and ↵tlatorre
the quantum efficiency
2018-09-19change output file format to YAMLtlatorre
2018-09-18stop fitting when the likelihood difference is less than 1e-2tlatorre
2018-09-18speed likelihood calculation up a bittlatorre
2018-09-18free muon_energy struct in test_muon_get_energy()tlatorre
2018-09-18fix heap overflow in interp1d()tlatorre
This commit fixes a potential heap overflow in interp1d() which could occur if x was very close to the last value in the xp array. The bounds check is now performed on the index rather than the x values.
2018-09-18fix typotlatorre
2018-09-18update PMT TTS and dark rate from pmt_response.dattlatorre
2018-09-18fix typotlatorre
2018-09-18fix PMTR bank offsetstlatorre
2018-09-18fix offsettlatorre
2018-09-18add functions to get the PMT reflectivity from the PMTR banktlatorre
2018-09-18add free_charge() to free memory used to interpolate the charge distributionstlatorre
2018-09-18add a ctrl-c handlertlatorre
2018-09-18update CHARGE_FRACTIONtlatorre
This commit updates the CHARGE_FRACTION value to now represent approximately the fraction of light reflected from each PMT. It also updates the value to be closer to the true value based on a couple of fits.
generated by cgit (git 2.34.1) at 2025-05-08 14:06:46 +0000