key_y, eventids=eventids_airplane, cmap=cmap, include_roi=False, lognorm=lognorm) #fig.savefig(key_x + key_y + '.svg', bbox_inches=0, transparent=False) fig.savefig('run%i_' % run + key_x + key_y + '.svg', bbox_inches=0, transparent=False) #import pdb; pdb.set_trace() else: ds = dataSlicer(runs, impulsivity_dset_key, time_delays_dset_key, map_direction_dset_key,\ curve_choice=0, trigger_types=trigger_types,included_antennas=[0,1,2,3,4,5,6,7],include_test_roi=False,\ cr_template_n_bins_h=200,cr_template_n_bins_v=200,\ impulsivity_n_bins_h=200,impulsivity_n_bins_v=200,\ time_delays_n_bins_h=150,time_delays_n_bins_v=150,min_time_delays_val=-200,max_time_delays_val=200,\ std_n_bins_h=200,std_n_bins_v=200,max_std_val=9,\ p2p_n_bins_h=128,p2p_n_bins_v=128,max_p2p_val=128,\ snr_n_bins_h=200,snr_n_bins_v=200,max_snr_val=35,\ n_phi=n_phi, range_phi_deg=range_phi_deg, n_theta=n_theta, range_theta_deg=range_theta_deg) plot_param_pairs = [['phi_best_h', 'elevation_best_h'], ['impulsivity_h', 'impulsivity_v']] for key_x, key_y in plot_param_pairs: print('Generating %s plot' % (key_x + ' vs ' + key_y)) fig, ax = ds.plotROI2dHist(key_x, key_y, cmap=cmap, include_roi=False, lognorm=lognorm)
impulsivity_dset_key = 'LPf_100.0-LPo_8-HPf_None-HPo_None-Phase_1-Hilb_0-corlen_65536-align_0-shortensignals-0-shortenthresh-0.70-shortendelay-10.00-shortenlength-90.00-sinesubtract_1' time_delays_dset_key = 'LPf_100.0-LPo_8-HPf_None-HPo_None-Phase_1-Hilb_0-corlen_65536-align_0-shortensignals-0-shortenthresh-0.70-shortendelay-10.00-shortenlength-90.00-sinesubtract_1' map_direction_dset_key = 'LPf_100.0-LPo_8-HPf_None-HPo_None-Phase_1-Hilb_0-upsample_32768-maxmethod_0-sinesubtract_1-deploy_calibration_15' #'LPf_100.0-LPo_8-HPf_None-HPo_None-Phase_1-Hilb_0-upsample_32768-maxmethod_0-sinesubtract_1-deploy_calibration_22-scope_belowhorizon'#'LPf_100.0-LPo_8-HPf_None-HPo_None-Phase_1-Hilb_0-upsample_32768-maxmethod_0-sinesubtract_1-deploy_calibration_15' trigger_types = [2] #[2] plot_maps = True sources_ENU, data_slicer_cut_dict = info.loadValleySourcesENU( ) #Plot all potential sources try: if False: ds = dataSlicer(runs, impulsivity_dset_key, time_delays_dset_key, map_direction_dset_key,\ curve_choice=0, trigger_types=trigger_types,included_antennas=[0,1,2,3,4,5,6,7],include_test_roi=False,\ cr_template_n_bins_h=200,cr_template_n_bins_v=200,\ impulsivity_n_bins_h=200,impulsivity_n_bins_v=200,\ time_delays_n_bins_h=500,time_delays_n_bins_v=500,min_time_delays_val=-140,max_time_delays_val=-90,\ std_n_bins_h=200,std_n_bins_v=200,max_std_val=9,\ p2p_n_bins_h=128,p2p_n_bins_v=128,max_p2p_val=128,\ snr_n_bins_h=200,snr_n_bins_v=200,max_snr_val=35) for source_key, cut_dict in data_slicer_cut_dict.items(): if source_key in include_sources: ds.addROI(source_key, cut_dict) ds.trackROICounts(roi_keys=None, time_bin_width_s=3600, plot_run_start_times=True) ds.trackROICounts(roi_keys=None, time_bin_width_s=3600, plot_run_start_times=False)