Commits
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 | |