; Prepare files for input to refinement
;
; Adapted from copyin.pam
; NOTE: can combine data sets by using input register, by appending to defocus-group list
;
; I/O PARAMETERS AND FILES ARE SET HERE:
;
; ------------------------- Parameter -------------------------
[firstgrp] = 1 ; first defocus group (can be used to append to data set)
; ------------------------ Input files ------------------------
FR G
[defgrps]sel_group_cclim ; Group selection file
; DF_GROUP NUMPARTS AVGDEFOCUS
; used used
FR G
[initial_sel_particles]sel_particles_{***[oldgrp]} ; Particle selection files
FR G
[initial_unaligned_images]../Alignment/data*** ; Unaligned stacked image files
FR G
[initial_aligned_images]../Alignment/dala01_*** ; Aligned stacked image files
; ----------------------- Output files -----------------------
FR G
[good_dir]../Reconstruction ; Output directory
FR G
[order_select][good_dir]/good-sel_group_cclim ; Group selection file
; ; (may contain appended groups)
FR G
[sel_particles][good_dir]/good-sel_particles_*** ; Particle selection files
FR G
[unaligned_images]../Alignment/good-data*** ; Unaligned stacked image files
FR G
[start_aligned_images]../Alignment/good-dala01_*** ; Aligned stacked image files
; ---------------------- END BATCH HEADER ----------------------
MD
TR OFF ; Decreases results file output
MD
VB OFF ; Decreases results file output
MD
SET MP
0
vm
echo "Preparing refinement files"; date
VM ; Crete refinement input directory
mkdir -p [good_dir]
UD N [num-grp] ; Get total number of defocus groups
[defgrps]
DO LB1 [i]=1,[num-grp] ; Loop over defocus groups ---------------------
UD IC [i], [oldgrp],x99,[defocus-value]
[defgrps] ; Doc file listing groups (input)
[newgrp] = [oldgrp] + [firstgrp] - 1
; clean up pre-existing files
de
[unaligned_images][newgrp]@
de
[start_aligned_images][newgrp]@
de
[sel_particles][newgrp]
ud n [num-particles]
[initial_sel_particles]
sd [i],[newgrp],[num-particles],[defocus-value]
[order_select]
; loop though particles
do lb2 [slice-counter]=1,[num-particles]
; get slice#
ud ic [slice-counter], [old-slice] ; , [particle-num]
[initial_sel_particles] ; WAS [df_dir][oldgrp]/[initial_sel_particles]
; write to selection doc
sd [slice-counter], [slice-counter], [old-slice] ; , [particle-num]
[sel_particles][newgrp]
; copy from total-particle stacks to good-particle stacks
cp
[initial_unaligned_images][oldgrp]@{******[old-slice]}
[unaligned_images][newgrp]@{******[slice-counter]}
cp
[initial_aligned_images][oldgrp]@{******[old-slice]}
[start_aligned_images][newgrp]@{******[slice-counter]}
lb2
; end particle-loop
ud ice
[initial_sel_particles]
sd e
[sel_particles][newgrp]
VM
echo " Accumulated group: {***[newgrp]}"
LB1
UD ICE ; Close doc file access
[defgrps]
sd e
[order_select]
;doc ren
;[order_select]
;[order_select]_renum
; renumber
at it
[order_select]
[order_select]
SD / GROUP PARTICLES DEFOCUS
[order_select]
sd e
[order_select]
VM
echo 'Done'; date
EN D
;
; Modified 2011-01-15
; 2011-01-07 (trs) -- refinement batch files using fewer and different inputs
; 2011-01-05 (trs) -- using more general SEL_PARTICLES as input
; 2009-01-12 (trs) -- can combine data sets by starting from defocus-group >1
; 2007-03-21 (trs) -- bug fix
; 2007-01-30 (trs) -- adapted from copyin.pam