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

switched to full generation of sigproc data

Browse source
This commit is contained in:
Sakimori 2025-07-21 16:46:53 -04:00
parent bb36dffab5
commit a071724015
No known key found for this signature in database
1 changed files with 21 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

28
signalgen.py
View file

@ -126,13 +126,27 @@ def addBurst(values):
show_dynamic(injectedSpectra, f"{basename} Dynamic Spectra and {values.dm} DM Pulse", save=True) show_dynamic(injectedSpectra, f"{basename} Dynamic Spectra and {values.dm} DM Pulse", save=True)
#now generate new filterbank file #now generate new filterbank file
newName = f"{basename}_injected.fil" newName = f"{basename}_injected.fil"
#build header from original file sigprocObj = make_sigproc_object(
headerDic = filterbankObj.your_header.__dict__ rawdatafile=newName, #d
#remove extraneous keys source_name="TEMP", #d
for key in ['filelist', 'filename']: nchans=filterbankObj.nchans,
headerDic.pop(key, None) foff=filterbankObj.your_header.foff, # MHz
headerDic['rawdatafile'] = newName fch1=filterbankObj.your_header.fch1, # MHz
sigprocObj = make_sigproc_object(**headerDic) tsamp=filterbankObj.your_header.tsamp, # seconds
tstart=filterbankObj.your_header.tstart, # MJD
src_raj=filterbankObj.your_header.ra_deg,
src_dej=filterbankObj.your_header.dec_deg,
machine_id=0,
nbeams=0,
ibeam=0,
nbits=filterbankObj.your_header.nbits,
nifs=1,
barycentric=0,
pulsarcentric=0,
data_type=0,
az_start=-1,
za_start=-1,
)
sigprocObj.write_header(newName) sigprocObj.write_header(newName)
sigprocObj.append_spectra(injectedSpectra, newName) sigprocObj.append_spectra(injectedSpectra, newName)
logger.info(f"{newName} successfully written.") logger.info(f"{newName} successfully written.")