From 87f87ffb95d76720e5daaa11f047cdc680874dbf Mon Sep 17 00:00:00 2001
From: sebastianp <sebastian.priebe@rub.de>
Date: Thu, 4 Aug 2016 16:07:01 +0200
Subject: [PATCH] [task] rewriting some functions to make them more flexible

---
 pylot/core/pick/compare.py | 49 +++++++++++++++++++++++++++++++-------
 1 file changed, 40 insertions(+), 9 deletions(-)

diff --git a/pylot/core/pick/compare.py b/pylot/core/pick/compare.py
index 1ef358cb..a3d85d69 100644
--- a/pylot/core/pick/compare.py
+++ b/pylot/core/pick/compare.py
@@ -345,13 +345,13 @@ class PDFstatistics(object):
     '''
     def __init__(self, directory):
         self.directory = directory
-        self.arraylen = 0
-        self.stations = {}
-        self.p_std = {}
-        self.s_std = {}
-        self.theta0 = []
-        self.theta1 = []
-        self.theta2 = []
+        #self.arraylen = 0
+       # self.stations = {}
+       # self.p_std = {}
+        #self.s_std = {}
+       # self.theta0 = []
+       # self.theta1 = []
+       # self.theta2 = []
 
 
     def readTheta(self, arname, dir, fnpattern):
@@ -375,6 +375,31 @@ class PDFstatistics(object):
     def makeFileList(self, fn_pattern='*'):
         self.evtlist = glob.glob1((os.path.join(self.directory)), '*.xml')
 
+    def nextPDF(self):
+
+        for evt in self.evtlist:
+            self.getPDFDict(self.directory, evt)
+            for station, pdfs in self.pdfdict.pdf_data.items():
+                try:
+                    yield pdfs['P']
+                except KeyError:
+                    yield np.nan
+                try:
+                    yield pdfs['S']
+                except KeyError:
+                    yield np.nan
+
+    def getQD(self,value):
+        pdfgen = self.nextPDF()
+        QDlist = []
+        for pdf in pdfgen:
+            if pdf == np.nan:
+                continue
+            QD = pdf.quantile_distance(value)
+            QDlist.append(QD)
+        return QDlist
+
+
 
     def getData(self):
         for evt in self.evtlist:
@@ -470,6 +495,12 @@ class PDFstatistics(object):
             fid.write(str(val)+'\n')
         fid.close()
 
+def main():
+    root_dir ='/home/sebastianp/Codetesting/xmls/'
+    Insheim = PDFstatistics(root_dir)
+    Insheim.makeFileList()
+    qdlist = Insheim.getQD(0.3)
+    print qdlist
 
-
-#if __name__ == "__main__":
\ No newline at end of file
+if __name__ == "__main__":
+    main()
\ No newline at end of file