aboutsummaryrefslogtreecommitdiff
AgeCommit message (Collapse)Author
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
2020-11-16DM_ENERGIES -> DM_MASSEStlatorre
2020-11-16update chi2 coverage bins to 10tlatorre
2020-11-16don't apply retrigger cut to MCtlatorre
2020-11-16loop over MC filenamestlatorre
2020-11-16update dm-search scripttlatorre
2020-11-03forgot to commit utils.pytlatorre
2020-11-03fix bug in get_events()tlatorre
This commit fixes an unnecessary call to set the index on the ev dataframe. This was causing issues when trying to process run 11903 since it didn't have any fits to merge into the ev dataframe.
2020-11-03initial commit of a script to do the dark matter searchtlatorre
2020-11-01update cat-grid-jobs to reprocess all zdabstlatorre
2020-11-01update qhl ratiotlatorre
2020-11-01update upper bound for external muon scaletlatorre
2020-11-01accidentally committed a print statementtlatorre
generated by cgit (git 2.34.1) at 2025-05-06 17:22:52 +0000