aboutsummaryrefslogtreecommitdiff
path: root/utils/dm-search
AgeCommit message (Collapse)Author
2021-01-15fix x axistlatorre
2021-01-15fix y axis againtlatorre
2021-01-15fix y axistlatorre
2021-01-14calculate total number of atmospheric events by summing flux_weightstlatorre
2021-01-05add --fast argument to dm-searchtlatorre
2021-01-05add --fasttlatorre
2021-01-04get rid of nhit_threshtlatorre
2021-01-03add --universe argument to dm-searchtlatorre
2021-01-03reduce memory usage by creating weights dict earlytlatorre
2021-01-03cache results from get_events()tlatorre
2020-12-25don't exit if some runs don't have MCtlatorre
2020-12-22move check earliertlatorre
2020-12-22fix typo from last committlatorre
2020-12-22go back to checking rhdr instead of evtlatorre
2020-12-21fix another rhdr checktlatorre
2020-12-21fix rhdr checktlatorre
2020-12-21only concat files with at least 1 eventtlatorre
2020-12-21use correct bins when calculation thresholdtlatorre
2020-12-21update dm-search to only look for mediator masses up to 1 GeVtlatorre
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-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-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-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-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-17small fixestlatorre
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-16fix latex errortlatorre
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-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-16don't apply retrigger cut to MCtlatorre
2020-11-16loop over MC filenamestlatorre
generated by cgit (git 2.34.1) at 2025-05-07 10:47:40 +0000