runs_old_sample = [ Run(path + "2018_05_31__13_25_38.txt"), Run(path + "2018_05_31__14_14_42.txt"), Run(path + "2018_05_31__14_18_09.txt"), Run(path + "2018_05_31__14_21_59.txt"), Run(path + "2018_05_31__14_25_21.txt"), Run(path + "2018_05_31__14_28_42.txt"), Run(path + "2018_05_31__14_32_14.txt") ] plot_runs(runs_nxt_85, title="NXT85 in air, adjusting sample height", labels=[ "height 1 \n(beam hitting center of sample)", "height 2", "height 3", "height 4", "height 5", "height 6", "height 7", "height 8", "height 9 \n(beam hitting near top of sample)" ], smooth=True, show=False) plot_runs(runs_m18_turn, title="M18 Turn in air, adjusting sample height", labels=[ "height 1 \n(beam hitting center of sample)", "height 2", "height 3", "height 4", "height 5", "height 6", "height 7 \n(beam hitting near top of sample)" ], smooth=True, show=False) plot_runs(runs_old_sample, title="Old Sample in air, adjusting sample height",
Run(path + "2018_06_13__16_47_03.txt") ] runs_marked = [ Run(path + "2018_06_14__12_32_29.txt"), Run(path + "2018_06_14__12_36_26.txt"), Run(path + "2018_06_14__12_40_10.txt"), Run(path + "2018_06_14__12_44_20.txt"), Run(path + "2018_06_14__12_48_09.txt") ] plot_runs( runs_unmarked, title= "LZ skived, unmarked side, not pushed all the way in, adjusting sample height", labels=[ "height 1 \n(beam hitting near center of sample)", "height 2", "height 3", "height 4\n(beam hitting near top of sample)" ], smooth=True, show=False) plot_runs(runs_unmarked_adjusted, title="LZ skived, unmarked side, adjusting sample height", labels=[ "height 1 \n(beam hitting near center of sample)", "height 2", "height 3", "height 4", "height 5\n(beam hitting near top of sample)" ], smooth=True, show=False) plot_runs(runs_marked, title="LZ skived, marked side, adjusting sample height",
Run(path + "2018_06_01__15_19_56.txt"), Run(path + "2018_06_01__15_31_34.txt") ] runs_pole = [ Run(path + "2018_06_04__14_48_21.txt"), Run(path + "2018_06_04__14_53_05.txt"), Run(path + "2018_06_04__14_57_33.txt"), Run(path + "2018_06_04__15_03_01.txt"), Run(path + "2018_06_04__15_07_27.txt") ] plot_runs( runs_no_pole, title= "NXT85 in air, 45 degrees, trying to repeat same sample height with no pole", labels=["test 1", "test 2", "test 3", "test 4"], smooth=True, show=False, ylabel="relative intensity") plot_runs( runs_pole, title= "NXT85 in air, 45 degrees, trying to repeat same sample height with pole", labels=["test 1", "test 2", "test 3", "test 4", "test 5"], smooth=True, show=False, ylabel="relative intensity") plt.show()
from file_reader import plot_runs, Run path = "new setup background tests/" no_walls = Run(path + "2018_05_25__13_44_05.txt") same_side = Run(path + "2018_05_25__13_57_52.txt") opposite_side = Run(path + "2018_05_25__14_06_33.txt") min_oil = Run(path + "2018_05_29__16_22_24.txt") runs = [no_walls, same_side, opposite_side, min_oil] plot_runs(runs, labels=[ "no walls in air", "walls opposite laser and detector in air", "walls opposite laser and on same side as detector in air", "new setup in mineral oil" ], show=False, rot=True) # copied from 10_16_background_tests.py, using data from background_measurements/background in mineral oil no sample old walls new walls lens tubes.txt distance_from_sample_to_photodiode = 5.435 photodiode_radius = (9 / 2.0) / 25.4 photodiode_solid_angle = np.pi * np.power(photodiode_radius, 2) / np.power( distance_from_sample_to_photodiode, 2) flux_i = 0.005940 * 100e-6 sensitivity = 100 * 1e-9 intensity_factor = sensitivity / (photodiode_solid_angle * flux_i) total_data = np.loadtxt( "old_setup_data/background in mineral oil no sample old walls new walls lens tubes.txt",
from file_reader import plot_runs, Run runs = [Run("3608219133.txt")] plot_runs(runs)
import matplotlib.pyplot as plt from file_reader import plot_runs, Run path = "laser_alignment_tests/" runs = [ Run(path + "2018_06_07__15_02_01.txt"), Run(path + "2018_06_07__15_15_04.txt"), Run(path + "2018_06_07__15_20_18.txt"), ] plot_runs(runs, title="Laser alignment tests", labels=["alignment 1", "alignment 2", "alignment 3"], smooth=True, show=False, ylabel="relative intensity") plt.show()
runs_unpolished = [ Run(path + "2018_06_07__15_53_10.txt"), Run(path + "2018_06_07__15_57_42.txt"), Run(path + "2018_06_07__16_01_32.txt"), Run(path + "2018_06_07__16_05_22.txt"), Run(path + "2018_06_07__16_09_12.txt"), Run(path + "2018_06_07__16_13_01.txt"), Run(path + "2018_06_07__16_17_38.txt"), Run(path + "2018_06_07__16_17_38.txt") ] runs_polished = [ Run(path + "2018_06_08__14_30_48.txt"), Run(path + "2018_06_08__14_34_56.txt"), Run(path + "2018_06_08__14_38_38.txt"), Run(path + "2018_06_08__14_42_37.txt"), Run(path + "2018_06_08__14_47_53.txt"), Run(path + "2018_06_08__14_52_12.txt") ] plot_runs(runs_unpolished, title="M18 turned, adjusting sample height", labels=["height 1 \n(beam hitting center of sample)", "height 2", "height 3", "height 4", "height 5", "height 6", "height 7", "height 8\n(beam hitting near top of sample)"], smooth=True, show=False) plot_runs(runs_polished, title="M18 polished, adjusting sample height", labels=["height 1 \n(beam hitting center of sample)", "height 2", "height 3", "height 4", "height 5", "height 6\n(beam hitting near top of sample)"], smooth=True, show=False) plt.show()