Index: /fact/tools/pyscripts/sandbox/dneise/introductory_course/outline.txt
===================================================================
--- /fact/tools/pyscripts/sandbox/dneise/introductory_course/outline.txt	(revision 14098)
+++ /fact/tools/pyscripts/sandbox/dneise/introductory_course/outline.txt	(revision 14098)
@@ -0,0 +1,29 @@
+FACT, Python & pyscripts
+introductory course
+
+This might in the end be the outline of some meetings at the TU Dortmund.
+The aim of these rather informal meetings is to give an introduction on the use
+of Python for looking at, and even analyse the FACT Data.
+
+The first session will deal with formal aspects such as:
+    * How to start
+        * programs to install
+        * source code to find
+        * data to find
+    
+    * How to plot
+        * Using Plotter classes, in plotters.py
+        * Using matplotlib.pyplot
+
+    * RawFile access
+        * Using RawData class, in pyfact.py
+        * Removing the DRS spike artifacts
+        * Waveform treatment Using Filters in fir_filter.py
+
+
+    * Our first Signal Extractors
+        * e.g. GlobalMaxFinder
+
+
+    * Our first Image Parameter Extractors
+        * e.g. SimpleArea, SimpleSize
Index: /fact/tools/pyscripts/sandbox/dneise/introductory_course/produce_npz.py
===================================================================
--- /fact/tools/pyscripts/sandbox/dneise/introductory_course/produce_npz.py	(revision 14098)
+++ /fact/tools/pyscripts/sandbox/dneise/introductory_course/produce_npz.py	(revision 14098)
@@ -0,0 +1,77 @@
+#!/usr/bin/python -tt
+#
+# Dominik Neise
+# TU Dortmund
+#
+# example for storing numpy arrays containing FACT analysis results
+# to file
+
+#import this at first, please
+from pyfact import RawData
+
+import os.path
+import numpy as np
+
+from drs_spikes import DRSSpikes
+from fir_filter import SlidingAverage
+from extractor  import GlobalMaxFinder
+from cleaners   import AmplitudeCleaner
+
+##############################################################################
+confirm_next_step = False# this is for user interaction
+
+data_file_name = '/media/daten_platte/FACT/data/20120229_144.fits.gz'
+if not os.path.isfile(data_file_name):
+    print 'not able to find file:', data_file_name
+    sys.exit(-1)
+calib_file_name = '/media/daten_platte/FACT/data/20120229_132.drs.fits.gz'
+if not os.path.isfile(calib_file_name ):
+    print 'not able to find file:', calib_file_name 
+    sys.exit(-1)
+
+run = RawData(data_file_name, calib_file_name)
+despike = DRSSpikes()
+smooth = SlidingAverage(8)
+extract = GlobalMaxFinder(40,200)
+cleaner = AmplitudeCleaner(45, 18)
+
+areas = []
+sizes = []
+for data,startcell,tt in run:
+    # trigger type 4 means 'physics event'
+    if tt==4:
+        data = despike(data)
+        data = smooth(data)
+        amplitude, time_of_max = extract(data)
+        survivors = cleaner(amplitude, return_bool_mask=False)
+
+        # this is not python like, should be done in a single line
+        # adding up the amplitude of all survivors
+        size = 0
+        for pixel in survivors:
+            size += amplitude[pixel]
+        
+        area = len(survivors)
+        
+        if area > 0:
+            areas.append( area )
+            sizes.append( size )
+        
+        # This is for ---------- USER INTERACTION -------------------------
+        if confirm_next_step:
+            user_input = raw_input("'q'-quit, 'r'-run, anything else goes one step")
+            number=None
+            try:
+                number=int(user_input)
+            except:
+                number=None
+            if user_input.find('q') != -1:
+                sys.exit(0)
+            elif user_input.find('r') != -1:
+                confirm_next_step = False
+            elif number!=None:
+                run += number
+        # ---------------END OF USER INTERACTION -------------------------
+
+outfname = raw_input("please enter name for outputfile")
+np.savez(outfname, Area=areas, Size=sizes)