aboutsummaryrefslogtreecommitdiff
AgeCommit message (Collapse)Author
2020-12-21update dm-search to only look for mediator masses up to 1 GeVtlatorre
2020-12-21add high orphan run listtlatorre
2020-12-21add dark matter run liststlatorre
2020-12-21add muon neutrino run liststlatorre
2020-12-21update bins for muon histogramstlatorre
This commit updates the bins for the muon histograms to not go too far above 2 GeV. The reason is that at these energies most muons will exit the detetor and my energy reconstruction doesn't work too well. I also updated chi2 and dm-search to use the flux_weight when sampling from the MC.
2020-12-16use a hash to merge weights with MC datatlatorre
2020-12-16print out livetime from run_info.log in dm-searchtlatorre
2020-12-16python 3 fixtlatorre
2020-12-15add code to reweight the tau neutrino eventstlatorre
This commit updates the code to reweight the MC data from tau neutrinos since I stupidly simulated the muon neutrino flux instead of the tau neutrino flux.
2020-12-09update discovery threshold againtlatorre
2020-12-09ensure discovery threshold is at least 0tlatorre
2020-12-09speed up chi2 and dm-searchtlatorre
2020-12-09use nanpercentile in get_events()tlatorre
2020-12-09discard negative weightstlatorre
2020-12-09fix bug introduced in previous committlatorre
This commit fixes a bug I introduced earlier in chi2 and dm-search since we want to remove runs not in the MC for both the signal and atmospheric samples.
2020-12-09fix bug when running chi2 --pulltlatorre
2020-12-09update logger.py to be python3 compatibletlatorre
This commite updates the Logger class to use str instead of basestring.
2020-12-09don't delete mc data from runs not in data when running chi2 --coveragetlatorre
2020-12-08speed up dm-searchtlatorre
2020-12-08use weighted MC when running dm-search --testtlatorre
2020-12-08fix more typostlatorre
2020-12-08fix another typotlatorre
2020-12-08fix typo from last committlatorre
2020-12-08update dm-search and chi2 to refit with GENIE systematicstlatorre
2020-12-08don't sample DM events when running --testtlatorre
2020-12-08update dm-search to fix bugtlatorre
This commit updates dm-search to fix a bug where I was returning lists from get_limits() but then comparing the best fit and discovery lists as if they were numpy arrays. I also updated how I calculate the expected number of events based on the results from doing a toy MC.
2020-12-08update chi2 script to pick the best universetlatorre
2020-11-30update code to work with python3tlatorre
This commit updates the python code to work with python 3 and with a newer version of matplotlib. - zip_longest -> izip_longest - fix tick marks for log plots - scipy.misc -> scipy.special
2020-11-30call get_events() with mc=True in submit-grid-jobstlatorre
2020-11-30get rid of superfluous --data argumenttlatorre
2020-11-30update submit-grid-jobs-queuetlatorre
2020-11-29add a script to calculate the cc, nc, and total event ratetlatorre
This commit adds a script called print-event-rate which calculates the event rate per year from the GENIE MCPL files. The livetime comes from the autosno run_info.log file. The output of the script is a latex table that I included in my thesis.
2020-11-17update burst cuttlatorre
2020-11-17small fixestlatorre
2020-11-16add script to plot distribution of number of high nhit orphanstlatorre
2020-11-16require signal events to have nhit_cal/nhit > 0.8tlatorre
2020-11-16add a new burst cuttlatorre
2020-11-16add salt pmt response file to submit-grid-jobs-queuetlatorre
2020-11-16add --run-list argument to chi2 and dm-searchtlatorre
This commit updates the chi2 and dm-search scripts to add the ability to pass a run list on the comand line.
2020-11-16add script to plot neutron nhit and delta t distributionstlatorre
2020-11-16bump version to 3tlatorre
2020-11-16fix latex errortlatorre
2020-11-16update fit to fit events from the salt phasetlatorre
This commit updates both the PMT response and optics code to be able to load constants from the salt phase.
2020-11-16fix bug by converting list to arraytlatorre
2020-11-16calculate livetime based on pulse gt and 10 MHz clocktlatorre
This commit updates get_events() to calculate the livetime based on both the number of pulse gt events and the 10 MHz clock and to return it in a dictionary stored with the dataframe. I also update dm-search so that the results are now reported as a function of events/cm^3/s. Also updated radius cut to be the AV radius.
2020-11-16add jdy, ut1, ut2, dte, and hmsc to hdf5 output filetlatorre
2020-11-16fix integrate_mc_hist() by dividing by bin widthtlatorre
2020-11-16use nlopt to find xopttlatorre
This commit updates both the chi2 and dm-search scripts to run nlopt at the end of the MCMC to find the best fit point. The reason for this is that for the dm-search script, we rely on the best fit point being correct when doing the discovery threshold analysis. To really get the best fit point with the MCMC we would need to run a ridiculous number of steps, so it's better to run a fewer amount of steps to get close and then run the minimization from there.
2020-11-16return np.inf if parameters are out of boundstlatorre
2020-11-16update lower bound for external muons and dm to 0tlatorre
generated by cgit (git 2.34.1) at 2025-05-06 15:06:37 +0000