; Window particles from tilt pairs.  Pawel A. Penczek 04/2002
;x50 is the size of the windowed image.
; noise is the noise image.
x50=105
x51=INT(X50/2)+1
CP
../mics/noise
_8
;
bl
_6
x50,x50
n
0.0
ma
_6
_3
50.0
d
e
1.0
x51,x51 
;
; particle counter
X91=0
; loop over micrographs
DO LB21 X90=1,3
; X87  - untilted micrograph
RR X87
8583,8585,8589
; X88  - tilted micrograph
RR X88
8582,8584,8588
;
;x55 is the number of particles selected from the micrographs.
UD N,X55
../mics/dct{***X90}
;
do lb1 x19=1,x55
ud ic,x19,x20,x21,x22
../mics/dcu{***X90}
x30=x21-x51
x31=x22-x51
ud ic,x19,x24,x25,x26
../mics/dct{***X90}
;untilted particles
x34=x25-x51
x35=x26-x51
wi
../mics/mic{****X87}
_1
x50,x50
x30,x31
ra
_1
_2
;
X91=X91+1
; ce fit
ce fit
_8
_2
_3
../win/unt{****x91}
;tilted particles
wi
../mics/mic{****X88}
_1
x50,x50
x34,x35
ra
_1
_2
; ce fit
ce fit
_8
_2
_3
../win/tlt{****x91}
lb1
ud ice
../mics/dcu{***X90}
ud ice
../mics/dct{***X90}
LB21
en