Science/psdlag-agn
1
0
Fork 0
mirror of https://asciireactor.com/otho/psdlag-agn.git synced 2024-11-25 00:45:06 +00:00
Code Issues Packages Projects Releases Wiki Activity

working on python stuff clag

Browse Source
This commit is contained in:
caes 2017-01-09 20:19:49 -05:00
parent 803df17307
commit 8829072a49
1 changed files with 8 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
scripts/analyze_lightcurve.py
View File

@ -20,8 +20,9 @@ lc1 = np.loadtxt(args[0])
# works if first two entries represent minimum spacing, from example # works if first two entries represent minimum spacing, from example
# dt = lc1[1,0] - lc1[0, 0] # dt = lc1[1,0] - lc1[0, 0]
# Time resolution determined from inspection and testing. # Time resolution determined from inspection and testing. This script
dt = 0.01 # does not expect evenly spaced data in time.
dt = 0.1
_ = plot(lc1[:,0], lc1[:,1]) _ = plot(lc1[:,0], lc1[:,1])
@ -36,10 +37,12 @@ lc1_strength = [lc1[i, 1] for i in index]
lc1_strength_err = [lc1[i, 2] for i in index] lc1_strength_err = [lc1[i, 2] for i in index]
# This would work if both curves are in same file # This would work if both curves are in same file
# Lc2 = [lc1[i, 3] for i in index] # lc2 = [lc1[i, 3] for i in index]
#Lc2e = [lc1[i, 4] for i in index] #Lc2e = [lc1[i, 4] for i in index]
# Load second light curve
lc2 = np.loadtxt(args[1])
#### Get the psd for the first light curve #### #### Get the psd for the first light curve ####
@ -86,7 +89,7 @@ errorbar(fqd[1:-1], psd1[1:-1], yerr=psd1e[1:-1], fmt='o', ms=10, label='fit')
## Now do the second light curve ## Now do the second light curve
P2 = clag.clag('psd', lc1_time, Lc2, Lc2e, dt, fqL) P2 = clag.clag('psd', lc1_time, lc2, Lc2e, dt, fqL)
psd2, psd2e = clag.optimize(P2, inpars) psd2, psd2e = clag.optimize(P2, inpars)
@ -99,7 +102,7 @@ psd2, psd2e = clag.optimize(P2, inpars)
### We also give it the calculated psd values as input ### ### We also give it the calculated psd values as input ###
Cx = clag.clag('cxd', Cx = clag.clag('cxd',
[list(i) for i in zip(lc1_time,lc1_time)], [list(i) for i in zip(lc1_time,lc1_time)],
[list(i) for i in zip(lc1_strength,Lc2)], [list(i) for i in zip(lc1_strength,lc2)],
[list(i) for i in zip(lc1_strength_err,Lc2e)], [list(i) for i in zip(lc1_strength_err,Lc2e)],
dt, fqL, psd1, psd2) dt, fqL, psd1, psd2)