From 0f86d58b2033d69a5207378deab82063744fcdb9 Mon Sep 17 00:00:00 2001 From: Abyss Halley Date: Fri, 22 Aug 2025 13:38:15 -0400 Subject: [PATCH] edited thresholds and pulsar period --- singlepulse.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/singlepulse.py b/singlepulse.py index 92bb06e..26c25a4 100644 --- a/singlepulse.py +++ b/singlepulse.py @@ -29,12 +29,12 @@ main_dir = os.path.join("/ldata","trunk") pulsar_dir_list = [os.path.join(main_dir, name) for name in ["data_2025-04-24_07-36-04", "data_2025-04-29_07-50-16", "data_2025-04-30_07-53-07", "data_2025-04-30_08-18-17", "data_2025-05-01_07-47-34"]] #AH: old variable: pulsar_dir -pulsar_period = 1.34789947290 * 1000 #period of the pulsar. AH: in ms +pulsar_period = 0.5 * 1000 #period of the pulsar. AH: in ms -#AH: Leaving these as default for now. -flux_threshold = 35 #minimum flux in mJy +#AH: Increasing FT and TT +flux_threshold = 80 #minimum flux in mJy distance_peaks = 6 #Number of samples/bins each pulse needs to be away from another pulse to be considered a detection. -time_threshold = 0.1 #error interval for matching pulses only at pulsar spin period apart. Set to 1 if don't want to apply this filter +time_threshold = 0.2 #error interval for matching pulses only at pulsar spin period apart. Set to 1 if don't want to apply this filter #Set options detrend_ft = True @@ -47,7 +47,7 @@ save_path_list = [os.path.join(name, "calib-pngs", "") for name in pulsar_dir_li t_sys = 18 #Fig 3 of the GBT document: https://greenbankobservatory.org/wp-content/uploads/2017/07/GBTpg.pdf) freq_mhz = 1400 -sky_coord = SkyCoord(339.492 * u.deg, 0.406 * u.deg, frame="galactic") #AH: For J1643-4522 +sky_coord = SkyCoord(339.492 * u.deg, 0.406 * u.deg, frame="galactic") #AH: For B1641-45 gsm = GlobalSkyModel() t_sky = gsm.get_sky_temperature(sky_coord, freq_mhz) #t_sky = 4.9805965