Source
2863
2863
2864
2864
os << "[" << itsImageName << "]" ;
2865
2865
os << " Peak residual (max,min) " ;
2866
2866
if( minresmask!=0.0 || maxresmask!=0.0 )
2867
2867
{ os << "within mask : (" << maxresmask << "," << minresmask << ") "; }
2868
2868
os << "over full image : (" << maxres << "," << minres << ")" << LogIO::POST;
2869
2869
2870
2870
os << "[" << itsImageName << "] Total Model Flux : " << getModelFlux() << LogIO::POST;
2871
2871
2872
2872
2873
+
Record* regionPtr=0;
2874
+
String LELmask("");
2875
+
Record thestats = SDMaskHandler::calcImageStatistics(*residual(), LELmask, regionPtr, True);
2876
+
Array<Double> maxs, mins;
2877
+
thestats.get(RecordFieldId("max"), maxs);
2878
+
thestats.get(RecordFieldId("min"), mins);
2879
+
os << LogIO::DEBUG1 << "Max : " << maxs << LogIO::POST;
2880
+
os << LogIO::DEBUG1 << "Min : " << mins << LogIO::POST;
2881
+
2873
2882
}
2874
2883
2875
2884
// Calculate the total model flux
2876
2885
Float SIImageStore::getMaskSum()
2877
2886
{
2878
2887
LogIO os( LogOrigin("SIImageStore","getMaskSum",WHERE) );
2879
2888
2880
2889
LatticeExprNode msum( sum( *mask() ) );
2881
2890
Float masksum = msum.getFloat();
2882
2891