Index: /fact/tools/pyscripts/fsc3.py =================================================================== --- /fact/tools/pyscripts/fsc3.py (revision 12610) +++ /fact/tools/pyscripts/fsc3.py (revision 12610) @@ -0,0 +1,128 @@ +#!/usr/bin/python26 +# +# FACT SLOW CONTROL DATA +# FSC_CONTROL_TEMPERATURE DATA +# +# Plot temperatures versus time +# Save plots to disk +# +# Werner Lustermann +# ETH Zurich + +import pyfits # module for reading fits files + +import numpy as np +import numpy.ma as ma # masked array module + +import matplotlib.pyplot as plt # for plotting + +class FSC_TEMPERATURE( object ): + + def __init__(self, dpath = '../temperature/', dfile = '20111103.FSC_CONTROL_TEMPERATURE.fits', plot_groups = True, show_plots = True ): + + self.dfile = dfile + self.dpath = dpath + self.fname = dpath + dfile + + self.Tgroups = ['T_sens', 'T_ps', 'T_crate', 'T_eth', 'T_back', 'T_aux' ] + self.rData() + + def rData( self ): + + try: + hdu = pyfits.open( self.fname ) + except: + print 'file: ', fname, ' does not exist! ==> BREAK!!!' + exit( 'unknown data file' ) + + self.tbdata = hdu[1].data + hdu.close() + + t = self.tbdata.field( 'Time' ) + tt = np.modf( t )[0] * 24. + # self.time = np.roll( tt, len( t[ t<12.] ) ) + for i in range( len( tt[ tt >= 12.] ) + 1, len( tt ) ): tt[i] += 24. + self.time = tt + + def tplot_group( self, mT, gr ): + """ plot maxumum, minimum and mean of a group of sensors """ + ax = plt.subplot(111) + ax.plot( self.time, np.max(mT, 1), 'b.', label = 'max', markersize = 1 ) + ax.plot( self.time, np.min(mT, 1), 'b.', label = 'min', markersize = 1 ) + ax.plot( self.time, np.mean(mT, 1), 'r.', label = 'mean', markersize = 1 ) + + ax.grid( 'on' ) + plt.title( 'Temperatures: ' + gr + ' ' + self.dfile[0:12], fontsize = 18 ) + plt.xlabel( 'time', fontsize = 14 ) + plt.ylabel( 'T in degC', fontsize = 14 ) + #plt.legend( loc = 4 ) + + # Shink current axis by 20% + box = ax.get_position() + ax.set_position([box.x0, box.y0, box.width * 0.9, box.height]) + plt.legend( fancybox=True, loc='center left', bbox_to_anchor=(1, 0.5), numpoints = 7 ) +# + def plot_groups( self ): + for i,group in enumerate( self.Tgroups ): + T = self.tbdata.field( group ) + maska = self.tbdata.field( group ) == 0.0 + maskb = self.tbdata.field( group ) > 50. + mask = maska | maskb + mT = ma.masked_array( T, mask ) + + plt.figure() + + self.tplot_group( mT, group ) + plt.savefig( self.dpath + self.dfile[0:-5] + '_' + group + '.png') + + +# +# + def plot_means( self ): + + fig = plt.figure() + ax = plt.subplot(111) + box = ax.get_position() + ax.set_position([box.x0, box.y0, box.width * 0.97, box.height]) + + + for i,group in enumerate( self.Tgroups ): + + #print 'i = ', i, group + T = self.tbdata.field( group ) + maska = self.tbdata.field( group ) == 0.0 + maskb = self.tbdata.field( group ) > 50. + mask = maska | maskb + mT = ma.masked_array( T, mask ) + + ax.plot( self.time, np.mean(mT, 1) , '.', label = group, markersize = 2 ) + + ax.grid( 'on' ) + plt.title( 'Mean Temperatures: ' + self.dfile[0:12], fontsize = 18 ) + plt.xlabel( 'time', fontsize = 14 ) + plt.ylabel( 'T in degC', fontsize = 14 ) + + # Shink current axis by 20% + box = ax.get_position() + ax.set_position([box.x0, box.y0, box.width * 0.97, box.height]) + plt.legend( fancybox=True, loc='center left', bbox_to_anchor=(1, 0.5), numpoints = 7 ) + plt.savefig( self.dpath + self.dfile[0:-5] + '_means' + '.png' ) + + + def show_plot( self ): + plt.show() + +if __name__ == '__main__': + + print 'FACT FSC_CONTROL TEMPERATURE DATA' + + fsct = FSC_TEMPERATURE() + fsct.plot_means() + fsct.plot_groups() + fsct.show_plot() + + + + + + Index: /fact/tools/pyscripts/fscvolt1.py =================================================================== --- /fact/tools/pyscripts/fscvolt1.py (revision 12610) +++ /fact/tools/pyscripts/fscvolt1.py (revision 12610) @@ -0,0 +1,98 @@ +#!/usr/bin/python26 +# +# FACT SLOW CONTROL DATA +# FSC_CONTROL_TEMPERATURE DATA +# +# Plot temperatures versus time +# Save plots to disk +# +# Werner Lustermann +# ETH Zurich + +import pyfits # module for reading fits files + +import numpy as np +import numpy.ma as ma # masked array module + +import matplotlib.pyplot as plt # for plotting + +class FSC_VOLTAGE( object ): + + def __init__(self, dpath = '../voltage/', dfile = '20111102.FSC_CONTROL_VOLTAGE.fits', plot_groups = True, show_plots = True ): + + self.dfile = dfile + self.dpath = dpath + self.fname = dpath + dfile + + self.Dgroups = ['FAD_Ud', 'FAD_U+', 'FAD_U-', 'FPA_Ud', 'FPA_U+', 'FPA_U-', 'ETH_U', 'FTM_U', 'FFC_U', 'FLP_U' ] + # self.Dgroups = ['FAD_Ud', 'FAD_U+', 'FAD_U-', 'FPA_Ud', 'FPA_U+', 'FPA_U-', 'ETH_U', 'FTM_U' ] + + self.rData() + + def rData( self ): + + try: + hdu = pyfits.open( self.fname ) + except: + print 'file: ', fname, ' does not exist! ==> BREAK!!!' + exit( 'unknown data file' ) + + self.tbdata = hdu[1].data + hdu.close() + + self.time = self.tbdata.field( 'Time' ) +# +# +# + def vplot_all( self ): + + for i, group in enumerate( self.Dgroups ): + + #print 'i = ', i, group + T = self.tbdata.field( group ) + # maska = self.tbdata.field( group ) == 0.0 + # maskb = self.tbdata.field( group ) > 50. + # mask = maska | maskb + # mT = ma.masked_array( T, mask ) + mT = T + + plt.figure(111) + ax = plt.subplot(111) + if i < 8: + ax.plot( self.time, np.mean(mT,1) , '.', label = group, markersize = 2 ) + else: + ax.plot( self.time, mT , '.', label = group, markersize = 2 ) + + ax.grid( 'on' ) + plt.title( 'DC Voltages: ' + self.dfile[0:12], fontsize = 18 ) + plt.xlabel( 'time', fontsize = 14 ) + plt.ylabel( 'Voltage in V', fontsize = 14 ) + + # Shink current axis by 20% + box = ax.get_position() + ax.set_position([box.x0, box.y0, box.width * 0.97, box.height]) + plt.legend( fancybox=True, loc='center left', bbox_to_anchor=(1, 0.5), numpoints = 7 ) + plt.savefig( self.dfile[0:-5] + '.png' ) + + + + def tshow_plots( self ): + plt.show() + +if __name__ == '__main__': + + print 'FACT FAD_CONTROL TEMPERATURE DATA' + + fadt = FSC_VOLTAGE() + fadt.vplot_all() + #fadt.tshow_plots() + + #fadt = FSC_VOLTAGE() + #fadt.vplot_all() + #fadt.tshow_plots() + + + + + + Index: /fact/tools/pyscripts/paux.py =================================================================== --- /fact/tools/pyscripts/paux.py (revision 12610) +++ /fact/tools/pyscripts/paux.py (revision 12610) @@ -0,0 +1,60 @@ +#!/usr/bin/python26 + +import glob +from os import path +from os import system + +#import pyfits + +from fscvolt1 import FSC_VOLTAGE +from fsc3 import FSC_TEMPERATURE + +class PAUX( object ): + + def __init__( self, year = 2011, month = 11, day = 0, dpoff = '../' ): + + self.year = str( year ) + self.month = str( month ) + self.day = str( day ) + self.dpoffset = dpoff + + def anatemp( self ): + + # fsc temperatures + self.dpath = '../temperature/' + self.searchpath = self.dpath + self.year + self.month + '*.FSC_CONTROL_TEMPERATURE.fits' + print 'searchpath', self.searchpath + dflist = glob.glob( self.searchpath ) + print 'dflist', dflist, 'shape: ', len(dflist) + for i, file in enumerate( dflist ): + file = file[15:] + print 'file: ', i, ': ', file + tdata = FSC_TEMPERATURE( dpath = self.dpath, dfile = file ) + tdata.plot_means() + tdata.plot_groups() + + + def anavolt( self ): + + # fsc voltages + self.dpath = '../voltage/' + self.searchpath = self.dpath + self.year + self.month + '*.FSC_CONTROL_VOLTAGE.fits' + # print 'searchpath', self.searchpath + dflist = glob.glob( self.searchpath ) + #print 'dflist', dflist, 'shape: ', len(dflist) + for i, file in enumerate( dflist ): + print 'file: ', i, ': ', file + vdata = FSC_VOLTAGE( dpath = self.dpath, dfile = file ) + vdata.vplot_all() + del vdata + + + + +if __name__ == '__main__': + + print 'FACT SLOW CONTROL DATA: PAUX' + + paux = PAUX(2011, 11) + paux.anatemp() +