diff --git a/plot_fsc.gpl b/plot_fsc.gpl index 3da9fd8f3b6bc93c8ede4257a3d702d96e18215e..486dcd7087c2728517c142e6dc0bab186b825835 100644 --- a/plot_fsc.gpl +++ b/plot_fsc.gpl @@ -2,7 +2,7 @@ set terminal pdfcairo size 8.5cm, 5.25cm font 'helvetica,12' enhanced unset key set output 'plots/fscs.pdf' -p for [f in system('ls -v [0-9]*/*threshold.fsc')] f t f, 0.143 +p for [f in system('ls -v [0-9]*/*threshold.fsc')] f t f, 0.5 set output 'plots/resolution.pdf' resolution=system('cat resolution.txt') set logscale x 10 diff --git a/postprocess.sh b/postprocess.sh index 8e9afb3fe0ee9c63484f8347681c7f4a0e76d392..7c4b001b7b2fe25af6c51bacaec70182c8c8a6a2 100755 --- a/postprocess.sh +++ b/postprocess.sh @@ -28,8 +28,8 @@ ${bindir}/gmconvert/gmconvert -igmm ${n}/${n}.gmm -imap threshold.map -omap ${n} e2proc3d.py ${map_name} ${n}/${n}.gmm.mrc.fsc --calcfsc ${n}/${n}.gmm.mrc e2proc3d.py threshold.map ${n}/${n}.gmm.mrc.threshold.fsc --calcfsc ${n}/${n}.gmm.mrc e2proc3d.py ${n}/${n}.gmm.mrc tmp.map --calcsf ${n}/${n}.gmm.mrc.sf -awk '$2<=0.143 && y>0.143{print 1.0/(x+(0.143-y)*($1-x)/($2-y))} {x=$1;y=$2}' < ${n}/${n}.gmm.mrc.fsc > ${n}/${n}.gmm.resolution -awk '$2<=0.143 && y>0.143{print 1.0/(x+(0.143-y)*($1-x)/($2-y))} {x=$1;y=$2}' < ${n}/${n}.gmm.mrc.threshold.fsc > ${n}/${n}.gmm.threshold.resolution +awk '$2<=0.5 && y>0.5{print 1.0/(x+(0.5-y)*($1-x)/($2-y))} {x=$1;y=$2}' < ${n}/${n}.gmm.mrc.fsc > ${n}/${n}.gmm.resolution +awk '$2<=0.5 && y>0.5{print 1.0/(x+(0.5-y)*($1-x)/($2-y))} {x=$1;y=$2}' < ${n}/${n}.gmm.mrc.threshold.fsc > ${n}/${n}.gmm.threshold.resolution rm -f $n/sub* rm -rf {.,$n}/tmp*