Source
133
133
134
134
if( validMask ) peakresidual = itsImages->getPeakResidualWithinMask();
135
135
else peakresidual = itsImages->getPeakResidual();
136
136
modelflux = itsImages->getModelFlux();
137
137
138
138
startpeakresidual = peakresidual;
139
139
startmodelflux = modelflux;
140
140
141
141
//Float nsigma = 150.0; // will set by user, fixed for 3sigma for now.
142
142
Float nsigma = loopcontrols.getNsigma();
143
-
os<<"robustrms nelements="<<robustrms.nelements()<<LogIO::POST;
143
+
//os<<"robustrms nelements="<<robustrms.nelements()<<LogIO::POST;
144
144
Float nsigmathresh;
145
145
if (robustrms.nelements()==0) {
146
146
nsigmathresh = 0.0;
147
147
} else{
148
148
nsigmathresh = nsigma * (Float)robustrms(IPosition(1,0));
149
149
}
150
150
151
151
Float thresholdtouse;
152
152
if (nsigma>0.0) {
153
153
// returns as an Array but itsImages is already single plane so