Commit b5f02160 authored by Pilar Cossio's avatar Pilar Cossio
Browse files

Fixing Norm in B prior

parent 37f46999
......@@ -1313,9 +1313,13 @@ int bioem_param::CalculateRefCTF()
// ********** Calculating normalized volumen element *********
if(!printModel){
param_device.volu = voluang * (myfloat_t) param_device.GridSpaceCenter * pixelSize * (myfloat_t) param_device.GridSpaceCenter * pixelSize
* gridCTF_phase * gridCTF_amp * gridEnvelop / ((2.f * (myfloat_t) param_device.maxDisplaceCenter+1.)) / (2.f * (myfloat_t) (param_device.maxDisplaceCenter + 1.)) / ((myfloat_t) numberGridPointsCTF_amp * gridCTF_amp )
/ ((myfloat_t) numberGridPointsCTF_phase * gridCTF_phase) / ((myfloat_t) numberGridPointsEnvelop * gridEnvelop ) / sqrt(2 * M_PI ) / param_device.sigmaPriorbctf ;
// All priors (uniform or not) normalized to 1
// The volume is the grid-spacing of the parameter / normalization
// the order is angles, displacement, ctf amplitud (all uniform) then env b & phase (non uniform) the sqrt(2) cancel out (see SI)
param_device.volu = voluang *
(myfloat_t) param_device.GridSpaceCenter * pixelSize * (myfloat_t) param_device.GridSpaceCenter * pixelSize / ((2.f * (myfloat_t) param_device.maxDisplaceCenter+1.)) / (2.f * (myfloat_t) (param_device.maxDisplaceCenter + 1.))
/ (myfloat_t) numberGridPointsCTF_amp
* gridEnvelop * gridCTF_phase / M_PI / param_device.sigmaPriorbctf / param_device.sigmaPriordefo ;
// cout << "VOLU " << param_device.volu << " " << gridCTF_amp << "\n";
// *** Number of total pixels***
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment