flux_i_air_vertical = 0.0004890 * 100e-6 flux_i_air_horizontal = 0.002120 * 100e-6 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor_air_vertical = sensitivity / (photodiode_solid_angle * flux_i_air_vertical) intensity_factor_air_horizontal = sensitivity / (photodiode_solid_angle * flux_i_air_horizontal) data_45_horizontal_air = make_data_by_run( "HP and VP for 45 deg in air.txt", 0, 80, intensity_factor_air_horizontal)[1] points_45_horizontal_air = make_points(data_45_horizontal_air[0], 0, 47, 1, 1, data_45_horizontal_air[1], 405, photodiode_solid_angle, photodiode_angular_width, "45 degrees horizontal water") data_45_vertical_air = make_data_by_run("HP and VP for 45 deg in air.txt", 0, 80, intensity_factor_air_vertical)[0] points_45_vertical_air = make_points(data_45_vertical_air[0], 0, 47, 1, 0, data_45_vertical_air[1], 405, photodiode_solid_angle, photodiode_angular_width,
sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor_horizontal_air = sensitivity / ( photodiode_solid_angle * flux_i_with_polarizer_horizontal_air) intensity_factor_horizontal_water = sensitivity / ( photodiode_solid_angle * flux_i_with_polarizer_horizontal_water) intensity_factor_horizontal_mineral_oil = sensitivity / ( photodiode_solid_angle * flux_i_with_polarizer_horizontal_mineral_oil) horizontal = True if horizontal: data_30_horizontal_air = make_data_by_run( "without_slit_7_18_horizontal/75, 60, 45, and 30 deg in air.txt", 0, 85, intensity_factor_horizontal_air)[0] data_45_horizontal_air = make_data_by_run( "without_slit_7_18_horizontal/75, 60, 45, and 30 deg in air.txt", 0, 85, intensity_factor_horizontal_air)[1] data_60_horizontal_air = make_data_by_run( "without_slit_7_18_horizontal/75, 60, 45, and 30 deg in air.txt", 0, 85, intensity_factor_horizontal_air)[2] data_75_horizontal_air = make_data_by_run( "without_slit_7_18_horizontal/75, 60, 45, and 30 deg in air.txt", 0, 85, intensity_factor_horizontal_air)[3] points_30_horizontal_air = make_points( data_30_horizontal_air[0], 0, 30, 1, 1, data_30_horizontal_air[1], 405, photodiode_solid_angle, photodiode_angular_width, "30 degrees horizontal air without slit")
make_all_points = True if make_all_points: make_points_100 = True if make_points_100: # volts * amps/volt flux_i = 0.002582 * 100e-6 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor = sensitivity / (photodiode_solid_angle * flux_i) data_30 = make_data_by_run( "glycerol_tests/glycerol 30,45,60,52.5,75.txt", -90, 90, intensity_factor)[4] data_45 = make_data_by_run( "glycerol_tests/glycerol 30,45,60,52.5,75.txt", -90, 90, intensity_factor)[3] data_60 = make_data_by_run( "glycerol_tests/glycerol 30,45,60,52.5,75.txt", -90, 90, intensity_factor)[2] data_52_5 = make_data_by_run( "glycerol_tests/glycerol 30,45,60,52.5,75.txt", -90, 90, intensity_factor)[1] data_75 = make_data_by_run( "glycerol_tests/glycerol 30,45,60,52.5,75.txt", -90, 90, intensity_factor)[0] points_100_30 = make_points(data_30[0], 0, 30, 1.47399, 0.5,
pol_arr = run_data_transpose[2] leg_arr = run_data_transpose[3] use_arr = run_data_transpose[4] # note that data is read out in reverse order of data taking sub_avg = False avg_angle = 10. points_arr = [] ii=0 for entry in file_data: filename = entry[0] flux_i = entry[1] # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor = sensitivity / (photodiode_solid_angle * flux_i) data_arr_i = make_data_by_run(filename, -35, 90, intensity_factor) for data in data_arr_i: if not get_entry_or_default(use_arr, ii, default=False): ii += 1 continue if sub_avg: y_avg = TSTR_fit.average_by_angle(data[0], data[1], avg_angle) y_sub = [data_i-avg_i for data_i, avg_i in zip(data[1],y_avg)] data_add = [data[0],y_avg] else: data_add = data leg_i = get_entry_or_default(leg_arr, ii, default="") theta_i = get_entry_or_default(theta_i_arr, ii, default=0.) n_i = get_entry_or_default(n_arr, ii, default=1.) pol_i = get_entry_or_default(pol_arr, ii, default=0.5) points_arr += make_points(data_add[0], 0, theta_i, n_i, pol_i, data_add[1], wavelength,
# these are from a different day, we didn't measure it for this test, we average them to compare flux_i_air_vertical_sample_1 = 0.0005240 * 100e-6 flux_i_air_horizontal_sample_1 = 0.002138 * 100e-6 flux_i = (flux_i_air_horizontal_sample_1 + flux_i_air_vertical_sample_1) / 2. # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor = sensitivity / (photodiode_solid_angle * flux_i) data_75 = make_data_by_run( "tube bumps analysis data/relative to normal angles.txt", -90, 90, intensity_factor)[0] data_45 = make_data_by_run( "tube bumps analysis data/relative to normal angles.txt", -90, 90, intensity_factor)[1] data_60 = make_data_by_run( "tube bumps analysis data/relative to normal angles.txt", -90, 90, intensity_factor)[2] data_30 = make_data_by_run( "tube bumps analysis data/relative to normal angles.txt", -90, 90, intensity_factor)[3] points_30 = make_points(data_30[0], 0, 30, 1, 0.5, data_30[1], 405, photodiode_solid_angle, photodiode_angular_width, "30 degrees") points_45 = make_points(data_45[0], 0, 45, 1, 0.5, data_45[1], 405,
photodiode_angular_width = photodiode_radius / distance_from_sample_to_photodiode make_all_points = True if make_all_points: # volts * amps/volt flux_i = 0.00618 * 100e-6 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor = sensitivity / (photodiode_solid_angle * flux_i) full_data = make_data_by_run("bubble tests all 3_12_18.txt", -90, 90, intensity_factor) full_data.reverse() test_1 = True if test_1: test_1_data = full_data[:] before_bubbles = test_1_data[:4] after_bubbles = test_1_data[4:6] after_sonic = test_1_data[6:] before_x = [] before_y = [] for data in before_bubbles: x_data = data[0] for x in x_data:
photodiode_angular_width = photodiode_radius / distance_from_sample_to_photodiode make_all_points = True if make_all_points: # volts * amps/volt flux_i = 0.005830 * 100e-6 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor = sensitivity / (photodiode_solid_angle * flux_i) data_fine_30 = make_data_by_run( "10_18_diffuse_reflectors/diffuse reflector 1500 grit in mineral oil 30, 45, 60, 75.txt", -90, 90, intensity_factor)[3] data_fine_45 = make_data_by_run( "10_18_diffuse_reflectors/diffuse reflector 1500 grit in mineral oil 30, 45, 60, 75.txt", -90, 90, intensity_factor)[2] data_fine_60 = make_data_by_run( "10_18_diffuse_reflectors/diffuse reflector 1500 grit in mineral oil 30, 45, 60, 75.txt", -90, 90, intensity_factor)[1] data_fine_75 = make_data_by_run( "10_18_diffuse_reflectors/diffuse reflector 1500 grit in mineral oil 75 full range.txt", -90, 90, intensity_factor)[0] data_coarse_30 = make_data_by_run( "10_18_diffuse_reflectors/diffuse reflector 120 grit in mineral oil 30 45 60.txt", -90, 90, intensity_factor)[2] data_coarse_45 = make_data_by_run(
# volts * amps/volt flux_i = 0.00238 * 100e-6 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor_mineral_oil_vertical = polarizer_transmission_vertical * sensitivity / ( photodiode_solid_angle * flux_i) intensity_factor_mineral_oil_horizontal = polarizer_transmission_vertical * sensitivity / ( photodiode_solid_angle * flux_i) horizontal = True if horizontal: data_30_horizontal_air = make_data_by_run( "horizontal_7_18/75, 60, 45, and 30 in air.txt", -90, 90, intensity_factor_air_horizontal)[0] data_45_horizontal_air = make_data_by_run( "horizontal_7_18/75, 60, 45, and 30 in air.txt", 0, 80, intensity_factor_air_horizontal)[1] data_60_horizontal_air = make_data_by_run( "horizontal_7_18/75, 60, 45, and 30 in air.txt", -90, 90, intensity_factor_air_horizontal)[2] data_75_horizontal_air = make_data_by_run( "horizontal_7_18/75, 60, 45, and 30 in air.txt", -90, 90, intensity_factor_air_horizontal)[3] points_30_horizontal_air = make_points(data_30_horizontal_air[0], 0, 30, 1, 1, data_30_horizontal_air[1], 405, photodiode_solid_angle,
photodiode_angular_width = photodiode_radius / distance_from_sample_to_photodiode make_all_points = True if make_all_points: # volts * amps/volt flux_i = 0.005830 * 100e-6 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor = sensitivity / (photodiode_solid_angle * flux_i) data_coarse_30 = \ make_data_by_run("10_18_diffuse_reflectors/diffuse reflector 120 grit in mineral oil 30 45 60.txt", -90, 90, intensity_factor)[2] data_coarse_45 = \ make_data_by_run("10_18_diffuse_reflectors/diffuse reflector 120 grit in mineral oil 30 45 60.txt", -90, 90, intensity_factor)[1] data_coarse_60 = \ make_data_by_run("10_18_diffuse_reflectors/diffuse reflector 120 grit in mineral oil 30 45 60.txt", -90, 90, intensity_factor)[0] points_coarse_30 = make_points(data_coarse_30[0], 0, 30, 1.47399, 0.5, data_coarse_30[1], 405, photodiode_solid_angle, photodiode_angular_width, "30 degrees, 120 grit sample, older data") points_coarse_45 = make_points(data_coarse_45[0], 0, 45, 1.47399, 0.5, data_coarse_45[1], 405, photodiode_solid_angle,
# volts * amps/volt flux_i_vertical = 0.002093 * 100e-6 flux_i_horizontal = 0.003907 * 100e-6 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor_vertical = sensitivity / (photodiode_solid_angle * flux_i_vertical) intensity_factor_horizontal = sensitivity / (photodiode_solid_angle * flux_i_horizontal) data_45_vertical = make_data_by_run( "10_12_polarization_data/45 degree vertical and horizontal mineral oil.txt", -90, 90, intensity_factor_vertical)[1] data_45_horizontal = make_data_by_run( "10_12_polarization_data/45 degree vertical and horizontal mineral oil.txt", -90, 90, intensity_factor_horizontal)[0] data_60_horizontal = make_data_by_run( "10_12_polarization_data/60 degree horizontal and vertical mineral oil.txt", -90, 90, intensity_factor_horizontal)[1] data_60_vertical = make_data_by_run( "10_12_polarization_data/60 degree horizontal and vertical mineral oil.txt", -90, 90, intensity_factor_vertical)[0] data_75_vertical = make_data_by_run( "10_12_polarization_data/75 degree vertical and horizontal mineral oil.txt", -90, 90, intensity_factor_vertical)[1] data_75_horizontal = make_data_by_run( "10_12_polarization_data/75 degree vertical and horizontal mineral oil.txt",
make_all_points = True if make_all_points: make_air_points = True if make_air_points: # volts * amps/volt flux_i = 0.00205 * 100e-6 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor = sensitivity / (photodiode_solid_angle * flux_i) data_30 = make_data_by_run("7_18_plots/30 deg in air with slit.txt", 10, 85, intensity_factor)[0] data_45 = make_data_by_run("7_18_plots/45 deg in air with slit.txt", 20, 80, intensity_factor)[0] data_60 = make_data_by_run("7_18_plots/60 in air with slit.txt", 0, 80, intensity_factor)[0] data_75 = make_data_by_run("7_18_plots/75 in air with slit.txt", 0, 80, intensity_factor)[0] points_30 = make_points(data_30[0], 0, 30, 1, 0.5, data_30[1], 405, photodiode_solid_angle, "30 degrees air with slit") points_45 = make_points(data_45[0], 0, 45, 1, 0.5, data_45[1], 405, photodiode_solid_angle, "45 degrees air with slit") points_60 = make_points(data_60[0], 0, 60, 1, 0.5, data_60[1], 405, photodiode_solid_angle,
# volts * amps/volt flux_i_no_tube = 0.005850 * 100e-6 flux_i_with_tube = 0.005940 * 100e-6 flux_i_older = 0.002603 * 100e-6 # from 9/12/2017 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor_no_tube = sensitivity / (photodiode_solid_angle * flux_i_no_tube) intensity_factor_with_tube = sensitivity / (photodiode_solid_angle * flux_i_with_tube) intensity_factor_older = sensitivity / (photodiode_solid_angle * flux_i_older) data_old_walls = make_data_by_run("background in mineral oil no sample old walls new walls lens tubes.txt", -90, 90, intensity_factor_no_tube)[2] data_new_walls = make_data_by_run("background in mineral oil no sample old walls new walls lens tubes.txt", -90, 90, intensity_factor_no_tube)[1] data_with_tube = make_data_by_run("background in mineral oil no sample old walls new walls lens tubes.txt", -90, 90, intensity_factor_with_tube)[0] data_older_without_cone = make_data_by_run("mineral oil background with and without cone.txt", -90, 90, intensity_factor_older)[0] data_older_with_cone = make_data_by_run("mineral oil background with and without cone.txt", -90, 90, intensity_factor_older)[1] points_old_walls = make_points(data_old_walls[0], 0, 90, 1.47399, 0.5, data_old_walls[1], 405, photodiode_solid_angle, photodiode_angular_width, "cardboard walls") points_new_walls = make_points(data_new_walls[0], 0, 90, 1.47399, 0.5, data_new_walls[1], 405, photodiode_solid_angle, photodiode_angular_width, "new walls") points_with_tube = make_points(data_with_tube[0], 0, 90, 1.47399, 0.5, data_with_tube[1], 405, photodiode_solid_angle, photodiode_angular_width, "new walls with cone of darkness") points_older_without_cone = make_points(data_older_without_cone[0], 0, 90, 1.47399, 0.5, data_older_without_cone[1], 405, photodiode_solid_angle, photodiode_angular_width, "older without cone") points_older_with_cone = make_points(data_older_with_cone[0], 0, 90, 1.47399, 0.5, data_older_with_cone[1],
photodiode_angular_width = photodiode_radius / distance_from_sample_to_photodiode make_all_points = True if make_all_points: # volts * amps/volt flux_i = 0.006083 * 100e-6 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor = sensitivity / (photodiode_solid_angle * flux_i) full_data = make_data_by_run("scraping_tests.txt", -90, 90, intensity_factor) full_data.reverse() test_1 = True if test_1: test_1_data = full_data[:7] before_data = test_1_data[:1] after_data = test_1_data[1:] before_x = [] before_y = [] for data in before_data: x_data = data[0] for x in x_data: before_x.append(x)
photodiode_solid_angle * flux_i) intensity_factor_water_horizontal = polarizer_transmission_vertical * sensitivity / ( photodiode_solid_angle * flux_i) # volts * amps/volt flux_i = 0.00238 * 100e-6 # amps per volt during measurements sensitivity = 100 * 1e-9 # product of this and measured voltage is (flux/str)/flux_i, flux in units of amps # intensity_factor * V = (V * sensitivity / photodiode_solid_angle) / flux_i intensity_factor_mineral_oil_vertical = polarizer_transmission_vertical * sensitivity / ( photodiode_solid_angle * flux_i) intensity_factor_mineral_oil_horizontal = polarizer_transmission_vertical * sensitivity / ( photodiode_solid_angle * flux_i) data_horizontal = make_data_by_run("45_deg_polarization_test.txt", 0, 80, intensity_factor_air_horizontal)[0] data_vertical = make_data_by_run("45_deg_polarization_test.txt", 0, 80, intensity_factor_air_vertical)[1] points_horizontal = make_points(data_horizontal[0], 0, 45, 1, 1, data_horizontal[1], 405, photodiode_solid_angle, "horizontal") points_vertical = make_points(data_vertical[0], 0, 45, 1, 0, data_vertical[1], 405, photodiode_solid_angle, "vertical") plot = True if plot: plot_with_TSTR_fit( points_vertical + points_horizontal, "Different Polarizations in Air, 45 degrees Without Slit")