zhetadelta/greenburstAux
1
0
Fork 0
Code Issues Pull requests Actions Packages Projects Releases Wiki Activity

edited thresholds and pulsar period

Browse source
This commit is contained in:
Abyss Halley 2025-08-22 13:38:15 -04:00
parent 31be71c5ae
commit 0f86d58b20
1 changed files with 5 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

10
singlepulse.py
View file

@ -29,12 +29,12 @@ main_dir = os.path.join("/ldata","trunk")
pulsar_dir_list = [os.path.join(main_dir, name) for name in 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"]] ["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 #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. #AH: Increasing FT and TT
flux_threshold = 35 #minimum flux in mJy 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. 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 #Set options
detrend_ft = True 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) t_sys = 18 #Fig 3 of the GBT document: https://greenbankobservatory.org/wp-content/uploads/2017/07/GBTpg.pdf)
freq_mhz = 1400 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() gsm = GlobalSkyModel()
t_sky = gsm.get_sky_temperature(sky_coord, freq_mhz) t_sky = gsm.get_sky_temperature(sky_coord, freq_mhz)
#t_sky = 4.9805965 #t_sky = 4.9805965