[daip] IRING q?

[Eric Greisen]

Lawrence Rudnick wrote:
> Eric - what is the formal definition of sigma in IRING output?  thanks.  lr
> 

                  MEAN(K) = MEAN(K) + VAL
                  MEANSQ(K) = MEANSQ(K) + VAL * VAL

where K is the ring counter

             MEAN(K) = MEAN(K) / NP(K)
             MEANSQ(K) = MEANSQ(K) / NP(K)
             SIGMA = SQRT (MAX (0.0D0, MEANSQ(K) - MEAN(K) * MEAN(K)))
             IF (NNZ(K).GT.1) SIGMA = SIGMA / SQRT (NNZ(K)-1.0)
             ERR(K) = SIGMA
             ERRS(K) = SIGMA * SIGMA
             IF (K.GT.1) ERRS(K) = ERRS(K) + ERRS(K-1)**2
             ERRS(K) = SQRT (ERRS(K))

where ERR is printed if one is plotting the individual rings and ERRS is 
printed if one is doing the cumulative sum.

RMFIT got some minor improvements recently and 31DEC13 had a serious bug 
fix.

Eric