Skip to content

Commit

Permalink
relion_reconstruct: repair --ewald --newbox.
Browse files Browse the repository at this point in the history
Old codes for --ewald --newbox cropped particles in the real space after
applying CTFP or Q but have been broken (probably) since
the introduction of optics groups in RELION 3.1.

Now, box size differences between the input particles
(read from an optics table) and the output box (requested by --newbox)
are handled in BackProjector by applying a scaling matrix to
the rotation matrix with obsModel.applyScaleDifference().
This works with/without --ewald.

Thanks to Martin Obr for the bug report and testing.
  • Loading branch information
biochem-fan committed May 10, 2023
1 parent c3480bd commit e16f796
Showing 1 changed file with 0 additions and 8 deletions.
8 changes: 0 additions & 8 deletions src/reconstructor.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -877,14 +877,6 @@ void Reconstructor::applyCTFPandCTFQ(MultidimArray<Complex> &Fin, CTF &ctf, Four

softMaskOutsideMap(Iapp, ROUND(mask_diameter/(angpix*2.)), (RFLOAT)width_mask_edge);

// Re-box to a smaller size if necessary....
if (newbox > 0 && newbox < YSIZE(Fin))
{
Iapp.setXmippOrigin();
Iapp.window(FIRST_XMIPP_INDEX(newbox), FIRST_XMIPP_INDEX(newbox),
LAST_XMIPP_INDEX(newbox), LAST_XMIPP_INDEX(newbox));

}
// Back into Fourier-space
transformer.FourierTransform(Iapp, Fapp, false); // false means: leave Fapp in the transformer
CenterFFTbySign(Fapp);
Expand Down

0 comments on commit e16f796

Please sign in to comment.