casa5
Public
Commits
| 132 132 | modelflux = itsImages->getModelFlux(); |
| 133 133 | |
| 134 134 | startpeakresidual = peakresidual; |
| 135 135 | startmodelflux = modelflux; |
| 136 136 | |
| 137 137 | // returns as an Array but itsImages is already single plane so |
| 138 138 | // the return rms contains only a single element |
| 139 139 | robustrms = itsImages->calcRobustRMS(); |
| 140 140 | //Float nsigma = 150.0; // will set by user, fixed for 3sigma for now. |
| 141 141 | Float nsigma = loopcontrols.getNsigma(); |
| 142 + | os<<"CURRENT NSIGMA="<<nsigma<<LogIO::POST; |
| 142 143 | Float nsigmathresh = nsigma * (Float)robustrms(IPosition(1,0)); |
| 143 144 | |
| 144 145 | Float thresholdtouse; |
| 145 146 | if (nsigma>0.0) { |
| 146 147 | thresholdtouse = max( nsigmathresh, loopcontrols.getCycleThreshold()); |
| 147 148 | } |
| 148 149 | else { |
| 149 150 | thresholdtouse = loopcontrols.getCycleThreshold(); |
| 150 151 | } |
| 151 152 | os << LogIO::DEBUG1<<"loopcontrols.getCycleThreshold()="<<loopcontrols.getCycleThreshold()<<LogIO::POST; |