RELION実装

2.0.3をベースに読み解いたソースコードになります。

RELION: Polisher

particle_polisher.cpp

void ParticlePolisher::run()
{

▸---// Fit straight lines through all beam-induced translations
       // ビーム照射により誘起された移動を通して、直線にフィットする
▸---if (fitting_mode != NO_FIT)
▸---▸---fitMovementsAllMicrographs();

▸---// Perform single-frame reconstructions and estimate dose-dependent B-factors
       // フレーム毎の三次元再構成とＢファクターの計算
▸---if (do_weighting)
▸---▸---calculateAllSingleFrameReconstructionsAndBfactors();
                                                                                                                                                         
▸---// Make a logfile in pdf format

　　// ファイルに書き出す

▸---generateLogFilePDF(); 

▸---// Write out the intermediately polished particles

　　// 選択された粒子を書き出す。

▸---polishParticlesAllMicrographs();

▸---// Now reconstruct with all polished particles: two independent halves, FSC-weighting of the sum of the two...

　　// 2つの独立した選ばれた粒子から再構成して、FSCを計算する。

▸---reconstructShinyParticlesAndFscWeight(1);

▸---if (verb > 0) 
▸---▸---std::cout << " done!" << std::endl;

}