Tomo CTF refinement¶
In this step we will show how to apply the refinement of the different parameters related to the CTF estimation. That is, tilt series projections defocus and astigmatism, scale and the asymmetrical and symmetrical aberrations. For a full description of the arguments, check relion_tomo_refine_ctf program.
Reference map at bin 1 and FSC data¶
Before running any of the tomo refinement jobs (this Tomo CTF refinement or Tomo frame alignment jobs), the user should run a Tomo reconstruct particle to estimate the reference map to feed these type of jobs. The main reason is that projection and reconstruction algorithms fully match and we will prevent suboptimal results because artifacts related to pseudo-subtomogram construction in the average map from the 3D auto-refine job.
On the other hand, a post-processed FSC data for the reference half-maps is optionally required to estimate the SNR. In case it is not provided, these programs internally calculate it without phase randomization so SNR would be slightly optimistic. Since the FSC data estimation is also integrated in the Tomo reconstruct particle job, we recommend to estimate it before running tomo refinement jobs.
Note that, independently of the binning factor level you were processing in the previous 3D auto-refine job, Tomo CTF refinement and Tomo frame alignment protocols both process data in the original pixel size (binning 1). Therefore, the reference map should always be reconstructed at this binning level as well as proper reference and FSC masks should be used.
Running the job¶
As described in previous jobs, we’re using the Tomo reconstruct particle output optimisation set file as input.
At this point, that optimisation set file should include the initial tomogram set we imported, the
run_data.star file from the 3D auto-refine job as particle set,
ReconstructParticleTomo/jobXXX/half<12>.mrc files as half maps and
ReconstructParticleTomo/jobXXX/PostProcess/postprocess.star file as FSC data.
Note the reference mask file used in previous 3D auto-refine job should also be included in the optimisation set file. Otherwise, or in case that reference mask is designed for a binning factor other than 1, a reference mask should be specifically provided.
On the I/O tab of the Tomo CTF refinement job-type set:
- Input optimisation set::
- Reference mask (optional)::
(If optimisation set does not include it)
On the Defocus tab, set:
- Box size for estimation (pix):
- Refine defocus?:
- Defocus search range (Å):
- Do defocus regularisation?:
- Defocus regularsation lambda:
- Refine constrast scale?:
- Refine scale per frame?:
- Refine scale per tomogram?:
On the Aberrations tab set:
- Refine odd aberrations?:
- Order of odd aberrations:
- Refine even aberrations?:
- Order of even aberrations:
On the Running tab, set:
- Number of MPI procs::
- Number of threads::
With these running parameters, the process should take around 10 minutes to finish.
Analysing the results¶
If you check the output folder
CtfRefineTomo/job013 you will find new
particles.star files with the refined CTF, scale and Zernike aberrations. To assess the result, it is recommended to run a new Tomo reconstruct particle job, with FSC estimation, using the new parameters. Note this reference map will also be used as input for the next Tomo frame alignment run. Compared to the previous FSC estimation, we should observe a slight improvement in the middle and high frequency ranges.