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

added gnuplot script and modified some things

Browse Source
This commit is contained in:
caes 2016-06-21 05:07:31 -04:00
parent 9055ee9b6e
commit 7e53d636e4
3 changed files with 36 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
inputfiles/calibration_final
View File

@ -1,10 +1,10 @@
2 2
data/calibration/lc_drive.dat 0 data/calibration/lc_drive.dat 0
data/calibration/lc_1.dat 0 data/calibration/lc_1.dat 0
8 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 8 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603
0 0
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
0:1 1 0:1 1
0 0
0 0
100 50 50 mcmc.dat 100 50 50 mcmc.dat

19
scripts/gnuplot/psd_freq Normal file
View File

@ -0,0 +1,19 @@
# plot 'file' with errorbars
# plot 'file' using 1:2:(sqrt($1)) with xerrorbars
# plot 'file' using 1:2:($1-$3):($1+$3):4:5 with xyerrorbars
#plot 'tmp.sourcePSD' using 1:2:3:4 with xyerrorbars
plot 'tmp.reprocPSD' using 1:2:3:4 with xyerrorbars
#plot 'tmp.ccPSD' using 1:2:3:4 with xyerrorbars
set logscale xy
set xrange [0.005:0.603]
# set yrange [:1]
# (x, y, ydelta),
# (x, y, ylow, yhigh),
# (x, y, xdelta),
# (x, y, xlow, xhigh),
# (x, y, xdelta, ydelta), or
# (x, y, xlow, xhigh, ylow, yhigh).
pause -1

24
scripts/power_lag_plot.pl
View File

@ -112,10 +112,10 @@ foreach ( sort { $a <=> $b } keys %function_bin ) {
my $μ = my $σ = <$outputfile>; my $μ = my $σ = <$outputfile>;
$μ =~ s/^([\-\+e0-9\.]+)+\s+[\-\+e0-9\.]+\s*$/$1/; $μ =~ s/^([\-\+e0-9\.]+)+\s+[\-\+e0-9\.]+\s*$/$1/;
$σ =~ s/^[\-\+e0-9\.]+\s+([\-\+e0-9\.]+)\s*$/$1/; $σ =~ s/^[\-\+e0-9\.]+\s+([\-\+e0-9\.]+)\s*$/$1/;
if ($μ<0) { # if ($μ<0) {
$μ = abs($μ); # $μ = abs($μ);
$function_bin{$_}{"φdiff_μ"} = PI; # $function_bin{$_}{"φdiff_μ"} = PI;
} # }
$function_bin{$_}{"source_PSD_μ"} = $μ; $function_bin{$_}{"source_PSD_μ"} = $μ;
$function_bin{$_}{"source_PSD_σ"} = $σ; $function_bin{$_}{"source_PSD_σ"} = $σ;
if ($debug) { if ($debug) {
@ -135,10 +135,10 @@ foreach ( sort { $a <=> $b } keys %function_bin ) {
my $μ = my $σ = <$outputfile>; my $μ = my $σ = <$outputfile>;
$μ =~ s/^([\-\+e0-9\.]+)+\s+[\-\+e0-9\.]+\s*$/$1/; $μ =~ s/^([\-\+e0-9\.]+)+\s+[\-\+e0-9\.]+\s*$/$1/;
$σ =~ s/^[\-\+e0-9\.]+\s+([\-\+e0-9\.]+)\s*$/$1/; $σ =~ s/^[\-\+e0-9\.]+\s+([\-\+e0-9\.]+)\s*$/$1/;
if ($μ<0) { # if ($μ<0) {
$μ = abs($μ); # $μ = abs($μ);
$function_bin{$_}{"φdiff_μ"} = PI; # $function_bin{$_}{"φdiff_μ"} = PI;
} # }
$function_bin{$_}{"reproc_PSD_μ"} = $μ; $function_bin{$_}{"reproc_PSD_μ"} = $μ;
$function_bin{$_}{"reproc_PSD_σ"} = $σ; $function_bin{$_}{"reproc_PSD_σ"} = $σ;
if ($debug) { if ($debug) {
@ -159,10 +159,10 @@ foreach ( sort { $a <=> $b } keys %function_bin ) {
my $μ = my $σ = <$outputfile>; my $μ = my $σ = <$outputfile>;
$μ =~ s/^([\-\+e0-9\.]+)+\s+[\-\+e0-9\.]+\s*$/$1/; $μ =~ s/^([\-\+e0-9\.]+)+\s+[\-\+e0-9\.]+\s*$/$1/;
$σ =~ s/^[\-\+e0-9\.]+\s+([\-\+e0-9\.]+)\s*$/$1/; $σ =~ s/^[\-\+e0-9\.]+\s+([\-\+e0-9\.]+)\s*$/$1/;
if ($μ<0) { # if ($μ<0) {
$μ = abs($μ); # $μ = abs($μ);
$function_bin{$_}{"φdiff_μ"} = PI; # $function_bin{$_}{"φdiff_μ"} = PI;
} # }
$function_bin{$_}{"cc_PSD_μ"} = $μ; $function_bin{$_}{"cc_PSD_μ"} = $μ;
$function_bin{$_}{"cc_PSD_σ"} = $σ; $function_bin{$_}{"cc_PSD_σ"} = $σ;
if ($debug) { if ($debug) {