diff options
| -rwxr-xr-x | utils/dm-search | 4 | ||||
| -rwxr-xr-x | utils/plot-muons | 2 | ||||
| -rwxr-xr-x | utils/plot-neutrons | 10 |
3 files changed, 8 insertions, 8 deletions
diff --git a/utils/dm-search b/utils/dm-search index 6c78d4e..872004b 100755 --- a/utils/dm-search +++ b/utils/dm-search @@ -721,11 +721,11 @@ if __name__ == '__main__': fig = plt.figure() for color, dm_particle_id in zip(('C0','C1'),(2020,2222)): - plt.plot(DM_MASSES[dm_particle_id],np.array(limits[dm_particle_id])/fiducial_volume/livetime,color=color,label='$' + ''.join([particle_id[int(''.join(x))] for x in grouper(str(dm_particle_id),2)]) + '$') + plt.plot(DM_MASSES[dm_particle_id],np.array(limits[dm_particle_id])*100**3*3600*24*365/fiducial_volume/livetime,color=color,label='$' + ''.join([particle_id[int(''.join(x))] for x in grouper(str(dm_particle_id),2)]) + '$') plt.gca().set_xscale("log") despine(fig,trim=True) plt.xlabel("Energy (MeV)") - plt.ylabel("Event Rate Limit (events/$\mathrm{cm}^3$/s)") + plt.ylabel("Event Rate Limit (events/$\mathrm{m}^3$/year)") plt.legend() plt.tight_layout() diff --git a/utils/plot-muons b/utils/plot-muons index 20ee7df..2616ad2 100755 --- a/utils/plot-muons +++ b/utils/plot-muons @@ -95,7 +95,7 @@ if __name__ == '__main__': evs.append(get_events(df.filename.values, merge_fits=True, nhit_thresh=args.nhit_thresh)) ev = pd.concat(evs) ev = correct_energy_bias(ev) - ev_mc = get_events(args.mc, merge_fits=True, apply_nhit_trigger=False) + ev_mc = get_events(args.mc, merge_fits=True, mc=True) ev_mc = correct_energy_bias(ev_mc) # Drop events without fits diff --git a/utils/plot-neutrons b/utils/plot-neutrons index 18108f8..291b890 100755 --- a/utils/plot-neutrons +++ b/utils/plot-neutrons @@ -99,30 +99,30 @@ if __name__ == '__main__': print("neutrons with nhit > 100") print(atm[atm.nhit_neutron >= 100][['run','gtid_neutron','nhit_neutron','nhit_cal_neutron']]) - fig = plt.figure() + fig = plt.figure(1) plt.hist(atm.nhit_cal_neutron.values,bins=np.linspace(0,100,101),histtype='step',color='C0',label="Data") weights = np.tile(len(atm)/len(atm_mc),len(atm_mc)) plt.hist(atm_mc.nhit_cal_neutron.values,bins=np.linspace(0,100,101),weights=weights,histtype='step',color='C1',label="Monte Carlo") plt.xlabel("Nhit") despine(fig,trim=True) plt.tight_layout() - plt.legend() + plot_legend(1) if args.save: plt.savefig("neutron_nhit_cal.pdf") plt.savefig("neutron_nhit_cal.eps") else: plt.title("Neutron Nhit Distribution") - fig = plt.figure() + fig = plt.figure(2) dt = (atm.gtr_neutron - atm.gtr)/1e6; bins = np.linspace(20e-3,250,101) plt.hist(dt,bins=bins,histtype='step',color='C0',label="Data") dt = (atm_mc.gtr_neutron - atm_mc.gtr)/1e6; plt.hist(dt,bins=bins,weights=weights,histtype='step',color='C1',label="Monte Carlo") - plt.xlabel(r"$\Delta t (ms)$") + plt.xlabel(r"$\Delta$ t (ms)") despine(fig,trim=True) plt.tight_layout() - plt.legend() + plot_legend(1) if args.save: plt.savefig("neutron_delta_t.pdf") plt.savefig("neutron_delta_t.eps") |