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

undid scale modification, run with original

Browse source
This commit is contained in:
Sakimori 2025-08-08 11:23:25 -04:00
parent 66886fac71
commit 8f8396cf95
No known key found for this signature in database
1 changed files with 1 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
signalgen.py
View file

@ -31,8 +31,6 @@ LOG_MIN_WIDTH = 0
LOG_MAX_WIDTH = 3 #in ms LOG_MAX_WIDTH = 3 #in ms
MAX_BURSTS = 4 MAX_BURSTS = 4
LOG_FIDUCIAL_WIDTH = -2.6 LOG_FIDUCIAL_WIDTH = -2.6
SCALE_SLOPE = 0.2
SCALE_OFFSET = 1
rng = np.random.default_rng() rng = np.random.default_rng()
#This function from https://josephwkania.github.io/will/examples/inject_pulse.html#Show-how-we-can-inject-a-pulse-into-a-GREENBURST-filterbank. #This function from https://josephwkania.github.io/will/examples/inject_pulse.html#Show-how-we-can-inject-a-pulse-into-a-GREENBURST-filterbank.
@ -176,7 +174,7 @@ def addBurst(values):
) )
#We need to scale the number of samples with pulse width to (hopefully) maintain constant peak flux #We need to scale the number of samples with pulse width to (hopefully) maintain constant peak flux
scaleFac = ((pWidth/(10**(LOG_FIDUCIAL_WIDTH)))*SCALE_SLOPE) + SCALE_OFFSET scaleFac = pWidth/(10**(LOG_FIDUCIAL_WIDTH))
scaledSamps = int(values.nsamp * scaleFac) scaledSamps = int(values.nsamp * scaleFac)
logger.info(f"Sampling pulse # {run+1} with width {round(pWidth,3)} s {scaledSamps} times.") logger.info(f"Sampling pulse # {run+1} with width {round(pWidth,3)} s {scaledSamps} times.")
pulse = pulseObj.sample_pulse(nsamp=scaledSamps) #30000 by default pulse = pulseObj.sample_pulse(nsamp=scaledSamps) #30000 by default