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

fixed cross spectrum errors

Browse Source
This commit is contained in:
caes 2017-06-22 15:19:38 -04:00
parent fdaea7f534
commit 630275a5a0
6 changed files with 6 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
scripts/compute_psdlags.sh
View File

@ -3,7 +3,7 @@
mkdir -p data/tables mkdir -p data/tables
mkdir -p logs mkdir -p logs
analysis_script="psdlag_5bin.py" analysis_script="psdlag_4bin.py"
ref_band="1367A" ref_band="1367A"
refpsd_tabfile=data/tables/psd_${ref_band}.tab refpsd_tabfile=data/tables/psd_${ref_band}.tab

1
scripts/psdlag_4bin.py
View File

@ -77,6 +77,7 @@ Cx = clag.clag('cxd10r',
# a good starting point generally # a good starting point generally
p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) ) p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) )
p, pe = clag.optimize(Cx, p) p, pe = clag.optimize(Cx, p)
p, pe = clag.errors(Cx, p, pe)
phi, phie = p[nfq:], pe[nfq:] phi, phie = p[nfq:], pe[nfq:]
lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd) lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd)

1
scripts/psdlag_5bin.py
View File

@ -80,6 +80,7 @@ echo_psd_err = p2e
Cx = clag.clag('cxd10r', [[t1,t2]], [[l1,l2]], [[l1e,l2e]], dt, fqL, p1, p2) Cx = clag.clag('cxd10r', [[t1,t2]], [[l1,l2]], [[l1e,l2e]], dt, fqL, p1, p2)
p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) ) # a good starting point generally p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) ) # a good starting point generally
p, pe = clag.optimize(Cx, p) p, pe = clag.optimize(Cx, p)
p, pe = clag.errors(Cx, p, pe)
phi, phie = p[nfq:], pe[nfq:] phi, phie = p[nfq:], pe[nfq:]
lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd) lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd)

1
scripts/psdlag_6bin.py
View File

@ -80,6 +80,7 @@ echo_psd_err = p2e
Cx = clag.clag('cxd10r', [[t1,t2]], [[l1,l2]], [[l1e,l2e]], dt, fqL, p1, p2) Cx = clag.clag('cxd10r', [[t1,t2]], [[l1,l2]], [[l1e,l2e]], dt, fqL, p1, p2)
p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) ) # a good starting point generally p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) ) # a good starting point generally
p, pe = clag.optimize(Cx, p) p, pe = clag.optimize(Cx, p)
p, pe = clag.errors(Cx, p, pe)
phi, phie = p[nfq:], pe[nfq:] phi, phie = p[nfq:], pe[nfq:]
lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd) lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd)

1
scripts/psdlag_7bin.py
View File

@ -79,6 +79,7 @@ echo_psd_err = p2e
Cx = clag.clag('cxd10r', [[t1,t2]], [[l1,l2]], [[l1e,l2e]], dt, fqL, p1, p2) Cx = clag.clag('cxd10r', [[t1,t2]], [[l1,l2]], [[l1e,l2e]], dt, fqL, p1, p2)
p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) ) # a good starting point generally p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) ) # a good starting point generally
p, pe = clag.optimize(Cx, p) p, pe = clag.optimize(Cx, p)
p, pe = clag.errors(Cx, p, pe)
phi, phie = p[nfq:], pe[nfq:] phi, phie = p[nfq:], pe[nfq:]
lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd) lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd)

1
scripts/psdlag_8bin.py
View File

@ -70,6 +70,7 @@ echo_psd_err = p2e
Cx = clag.clag('cxd10r', [[t1,t2]], [[l1,l2]], [[l1e,l2e]], dt, fqL, p1, p2) Cx = clag.clag('cxd10r', [[t1,t2]], [[l1,l2]], [[l1e,l2e]], dt, fqL, p1, p2)
p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) ) # a good starting point generally p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) ) # a good starting point generally
p, pe = clag.optimize(Cx, p) p, pe = clag.optimize(Cx, p)
p, pe = clag.errors(Cx, p, pe)
phi, phie = p[nfq:], pe[nfq:] phi, phie = p[nfq:], pe[nfq:]
lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd) lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd)