Index: fact/tools/pyscripts/examples/test_dneise.py
===================================================================
--- fact/tools/pyscripts/examples/test_dneise.py	(revision 12817)
+++ fact/tools/pyscripts/examples/test_dneise.py	(revision 12817)
@@ -0,0 +1,52 @@
+#!/usr/bin/python
+#
+# Dominik Neise
+# TU Dortmund
+#
+# test for freshly implemented methods in pyfact.
+# subject to frequent change!
+
+from pyfact import *
+
+dfname = '/data03/fact-construction/raw/2011/11/24/20111124_121.fits'
+calfname = '/data03/fact-construction/raw/2011/11/24/20111124_111.drs.fits'
+# access the data
+rd = rawdata( dfname, calfname )
+
+print 'dfname:  ', dfname
+print 'calfname:', calfname
+print 'NEvents: ', rd.NEvents
+
+def loop_acal( Neve = 1000 ):
+    """
+    bla
+    """
+    if rd.NEvents < Neve:
+        print 'data file contains not enough events'
+        exit( 0 )
+    else:
+        for i in range( Neve ):
+        #if np.mod(i,1000):
+        #    print 'Event: ', i
+            rd.next()
+            
+            rd.filterSlidingAverage()
+            print rd.smoothData
+            rd.filterCFD()
+            print rd.cfdData
+            rd.findPeak()
+            print rd.maxPos
+            print rd.maxAmp
+            rd.sumAroundPeak()
+            print rd.integral
+            #print len(rd.maxPos)
+            #print len(rd.maxAmp)
+            print 'Trigger Type', rd.trigType
+
+print (rd.filterSlidingAverage.__doc__)
+print (rd.filterCFD.__doc__)
+print (rd.findPeak.__doc__)
+print (rd.sumAroundPeak.__doc__)
+
+loop_acal(2)
+
Index: fact/tools/pyscripts/examples/test_fits.py
===================================================================
--- fact/tools/pyscripts/examples/test_fits.py	(revision 12817)
+++ fact/tools/pyscripts/examples/test_fits.py	(revision 12817)
@@ -0,0 +1,27 @@
+#!/usr/bin/python
+#
+#
+from ROOT import *
+
+gSystem.Load('fits_h.so')
+
+from ROOT import *
+
+# now we have full access to all ROOT classes and to the class fits
+
+# access a raw Data file
+
+dfname = '/data03/fact-construction/raw/2011/11/24/20111124_121.fits'
+df = fits( dfname )
+
+NROI  = df.GetUInt('NROI')
+NEvents = df.GetNumRows()
+
+print 'accessing data file: ', dfname
+print 'NROI: ', NROI
+print 'NEvents: ', NEvents
+
+# up to here everything is fine
+# the question now is how to really access the data
+
+df.SetPtrAddress('EventNum', a)
Index: fact/tools/pyscripts/examples/test_save.py
===================================================================
--- fact/tools/pyscripts/examples/test_save.py	(revision 12817)
+++ fact/tools/pyscripts/examples/test_save.py	(revision 12817)
@@ -0,0 +1,27 @@
+#!/usr/bin/python
+#
+# Werner Lustermann
+# ETH Zurich
+#
+# test / example: save multiple arrays into a file
+#
+import numpy as np
+
+a = np.zeros( 10 )
+b = np.ones( (10, 3) )
+
+print 'a', a
+print 'b', b
+
+np.savez( 'save.npz', amplitude = a, time = b )
+
+f = np.load( 'save.npz' )
+
+print 'f.files'
+print f.files
+
+ampl = f['amplitude']
+time = f['time']
+
+print 'ampl', ampl
+print 'time', time
Index: fact/tools/pyscripts/examples/where.py
===================================================================
--- fact/tools/pyscripts/examples/where.py	(revision 12817)
+++ fact/tools/pyscripts/examples/where.py	(revision 12817)
@@ -0,0 +1,20 @@
+#!/usr/bin/python
+#
+# test of numpy where
+import numpy as np
+
+
+
+NPIX = 20
+a = np.random.rand( NPIX )
+
+
+
+a_gt_thr1 = np.where( a > 0.5 )
+a_gt_thr2 = np.where( a > 0.9 )
+
+for index in a_gt_thr1[0] :
+    print index
+
+print 'a_gt_thr1', a_gt_thr1
+print 'a_gt_thr2', a_gt_thr2