#!/usr/bin/python
import os,sys,popen2
import numpy

if len(sys.argv)!=2:
    print 'get_mean_rms.py  genxxx'
    sys.exit(-1)

#dir = sys.argv[1]
#if os.path.isdir(dir):
#    pass
#else:
#    print 'ERROR:',dir,'is not a valid directory'
#    sys.exit(-1)

sz = []
#o,i = popen2.popen2('zcat '+dir+'/*.log.gz | grep Spin')
o,i = popen2.popen2('zcat '+sys.argv[1]+'run*.log.gz | grep Spin')
for l in o.readlines():
    ls = l.split()
    if len(ls)==4 and ls[0]=='Spin':
        sz.append(float(ls[3]))
print 'spins =',sz
print '\nn=',len(sz)
print '\nmean =',numpy.mean(sz)
print 'std (x1e-4) =',numpy.std(sz)*1E4