; Compute the FSC
;
; SOURCE: fsc.spi
;
; PURPOSE: Compute the FSC.
;
; MASTER COPY: spider/docs/techs/recon/newprogs/
;
; I/O PARAMETERS AND FILES ARE SET HERE:
;
; ----------------- Input files -----------------------------------

FR G
[order]sel_group_cclim            ; Selection file for defocus groups

FR G
[voleven]df{***[grp]}/vol01_even  ; Volume from even numbered particles

FR G
[volodd]df{***[grp]}/vol01_odd    ; Volume from odd numbered particles

; -------------------------- Output files ---------------------------

FR G
[out]df{***[grp]}/doccmp01        ; Output doc file (one / defocus group)

; -------------- END BATCH HEADER ---------------------------------

FR G
[tmp]tmpres89             ; Temp doc file (deleted)

UD N [ngrp]               ; Find number of defocus groups
[order]

DO [grp]= 1,[ngrp]        ; Loop over all defocus groups
   DE
   [out]
   DE
   [tmp]

   RF 3
   [voleven]              ; First input volume
   [volodd]               ; Second input volume
   (0.5)                  ; Ring width
   (0.2,2)                ; Lower, upper scale factors
   C                      ; C = missing cone
   (90)                   ; Max. tilt angle
   (3)                    ; Factor for noise comparison
   [tmp]                  ; Output doc file

   DOC REN                ; Renumber doc file keys
   [tmp]                  ; Input  doc file
   [tmp]                  ; Output doc file

   UD N [nkeys]           ; Get number of keys in tmp file
   [tmp]

   ; Title for output doc file
   SD /    NORM'D FREQ    DPR        FSC        FSCCRIT     VOXELS
   [out]

   DO [v61] = 1,[nkeys]
      UD [v61],[v62],[v63],[v64],[v65],[v66]
      [tmp]

      SD [v61,[v6]2,[v63],[v64],[v65],[v66]
      [out]
   ENDDO

   UD E
ENDDO

DE
[tmp]

EN D
;