Version used for data taking,

Also added some basic py scripts for presenting the data

Author: Sam

camac_control.py

@@ -21,6 +21,9 @@ class CAMACError(Exception):
     pass
 
 
+class CAMACreadfail(Exception):
+    pass
+
 def status(sock): 
     """
     returns the status as received from the socket, sock. 
@@ -52,13 +55,16 @@ def status(sock):
 
     return res
 
-def status_check(data):
+def status_check(data, location=""):
     """
     Checks the status and qresp for the passed packed"""
-    if data["status"][1] != 1:
-        raise CAMACError("Status failure code: %i" % data["status"][1])
-    elif data["qresp"][1]  != 0x8000:
-        raise CAMACError("Qresp failure %i" % data["qresp"][1])
+    if (data["status"][1] == 90) and (data["qresp"][1] == 0):
+        # warning only (no xresp)
+        return
+    elif data["status"][1] != 1:
+        raise CAMACError("Status failure code: %i %s" % (data["status"][1], location))
+    elif (data["qresp"][1] != 0x8000):
+        raise CAMACError("Qresp failure %i %s" % (data["qresp"][1], location))
 
 def cccc(sock):
     """
@@ -159,9 +165,9 @@ def rsync(sock, maxTries = 10):
     "readClrGrp1":       (2,  "subaddr"),
     "readClrGrp2":       (3,  "subaddr"),
     "testLAM":           (8,  ),
-    "clearGrp1":         (9,  "subaddr"),
+    "clearGrp1":         (9,  ),
     "clearLAM":          (10, ),
-    "clearGrp2":         (11, "subaddr"),
+    "clearGrp2":         (11, ),
     "overWriteGrp1":     (16, "subaddr", "data"),
     "overWriteGrp2":     (17, "subaddr", "data"),
     "maskOverWriteGrp1": (18, "subaddr", "data"), 
@@ -181,13 +187,13 @@ def naf(n, f, a = None, data = None ):
     # get the function info required (number and whether an address etc are needed) 
     funcCode = (f, ) if str(f).isdigit() else camacFunction[f] 
     if (("subaddr" in funcCode) and (a == None)):
-        raise Exception("function requires a subaddress")
+        raise CAMACError("function requires a subaddress")
     elif (("data" in funcCode) and (data == None)):
-        raise Exception("function requires a data")
+        raise CAMACError("function requires data")
 
     a = a if a else 0
     data = pack('H', data) if (data >= 0) else b''
-    # print funcCode[0], n, a
+    # print 'n is', n, 'a is', a, 'f is', funcCode[0], '\n'
     res = funcCode[0] <<10 | n <<5 | a << 1 | 0 
     res = pack('H', res)
     return res + data
@@ -217,28 +223,29 @@ def main_test(sock):
     ccci(sock)
     print("\ncompleted i")
 
-    cssa(sock, n=22, f=26, a=0)
+    print cssa(sock, n=22, f='clearLAM', a=0, data=0)
     # cssa(sock, n=22, f=26, a=1)
-        # cssa(sock, n=22, f=9, a=0)
-        # cssa(sock, n=22, f=11, a=0)
-    count = 0
-    while(True):
-        t = cssa(sock, n=22, f=8, a=0)
-        cssa(sock, n= 22, f='clearLAM')
-        cssa(sock = sock, n = 23, f = 16, a = 0, data=1) 
-        cssa(sock = sock, n = 23, f = 16, a = 0, data=0)
-        sleep(1)
-        if t["qresp"][1] != 0:
-            # print "wow!"
-            # print count
-            # print t
-            dat = cssa(sock, 21, "readGrp1", 8)
-            cssa (sock, 21, "clearGrp1", 8)
-            print dat["data"]
-        else:
-            count += 1
-            
-        
+    # cssa(sock, n=22, f=9, a=0)
+    # cssa(sock, n=22, f=11, a=0)
+    # count = 0
+    # while(True):
+    #   t = cssa(sock, n=22, f=8, a=0)
+    #   cssa(sock, n= 22, f='clearLAM')
+    #   break
+    #   cssa(sock = sock, n = 23, f = 16, a = 0, data=1) 
+    #   cssa(sock = sock, n = 23, f = 16, a = 0, data=0)
+    #   sleep(1)
+    #   if t["qresp"][1] != 0:
+    #       # print "wow!"
+    #       # print count
+    #       # print t
+    #       dat = cssa(sock, 21, "readGrp1", 8)
+    #       cssa (sock, 21, "clearGrp1", 8)
+    #       print dat["data"]
+    #   else:
+    #       count += 1
+    #           
+    #       
 
     # msg = gettop(flags = 0x0300, ops = 5, cmd = "cmd_attach_lam", mod = 21)
     # msg += naf(22, 0, 0) + naf(21,0,8) + naf(21, 0, 9)

display.py

@@ -0,0 +1,39 @@
+#! /usr/bin/python
+
+"""
+reads in the data file and plots it in ROOT
+"""
+
+import ROOT
+
+if __name__ == '__main__':        
+    file_numbers = (262,) #range(209)
+    file_fmt = 'data/test_%03i.txt'
+    
+    # open the files and read them into a histogram
+    hist = []
+    for ch in range(4):
+        h = ROOT.TH1F("ch%i" % ch, "Beam on 1nA pos: (0, 20)", 2**12, 0, 2**12 + 1)
+        hist.append(h)
+    
+    for n in file_numbers:
+        file = open(file_fmt % n, "r")
+        for line in file:
+            if '-' in line: continue
+            
+            dat = line.split()
+            for ch in range(4):
+                hist[ch].Fill(float(dat[ch]))
+    
+    canvas = ROOT.TCanvas("dat", "dat")
+    canvas.Divide(2, 2)
+    
+    for pad in range (1, 5):
+        canvas.cd(pad)
+        hist[pad - 1].Draw()
+        
+    run = True
+    while run:
+        d = raw_input()
+        if d:
+            run = False

displayXY.py

@@ -0,0 +1,37 @@
+#! /usr/bin/python
+
+"""
+reads in the data file and plots it in ROOT
+"""
+
+import ROOT
+
+if __name__ == '__main__':        
+    file_numbers = (2,3,4)
+    file_fmt = 'test_dat/test_%03i.txt'
+    
+    # open the files and read them into a histogram
+    hist = []
+    for ch in range(1, 4):
+        h = ROOT.TH2F("ch1&ch%i" % ch, "pedestal_data", 2000, 0, 2000, 2000, 0, 2000)
+        hist.append(h)
+    
+    for n in file_numbers:
+        file = open(file_fmt % n, "r")
+        for line in file:
+            dat = line.split()
+            for ch in range(1, 4):
+                hist[ch - 1].Fill(float(dat[0]), float(dat[ch]))
+    
+    canvas = ROOT.TCanvas("dat", "dat")
+    canvas.Divide(2, 2)
+    
+    for pad in range (1, 4):
+        canvas.cd(pad)
+        hist[pad - 1].Draw()
+        
+    run = True
+    while run:
+        d = raw_input()
+        if d:
+            run = False

fileUpdater.py

@@ -33,24 +33,33 @@ class FileUpdater(file):
     is also plotted in n_data root histograms, the assumption being
     that each n_data corresponds to a channel"""
     file_number = 0 # class attribute
-    canvas = None,
+    canvas = None
 
     def __init__(self, file_path = 'r_%03i.txt', mode="w", n_dat=4, header = '', TH1args=def_TH1, update_on=100):
-        """if file_path contains '%' then that is formatted as the file_number
+        """
+        if file_path contains '%' then that is formatted as the file_number
         multiple uses of '%' cause an error. If included as an argument
         the header is added to the beginning of the file. Arguments for initialising 
         the histogram should be included in the dictionary TH1F_kargs"""
+        
         # set instance attributes
         self.n_ch = n_dat
-        self.n_run = FileUpdater.file_number
         self.entries = 0
         self.update_on = update_on
         FileUpdater.file_number +=1
-        if '%' in file_path: file_path = file_path % self.n_run
+        if '%' in file_path: 
+            while os.path.isfile(file_path % FileUpdater.file_number): 
+                FileUpdater.file_number += 1
+            file_path = file_path % FileUpdater.file_number
+        self.n_run = FileUpdater.file_number
 
         # open the file
+        print "creating file %s" % file_path
+        print repr(file_path), repr(mode)
         file.__init__(self, file_path, mode)
-        if header: file.write(self, header)
+        if header: 
+            if header[-1] != '\n': header += '\n'
+            file.write(self, header)
 
         # initilse the histograms
         self.hist = []
@@ -62,7 +71,7 @@ def __init__(self, file_path = 'r_%03i.txt', mode="w", n_dat=4, header = '', TH1
         FileUpdater.canvas.Divide(*frameFactorer(self.n_ch))
         self.update()
 
-    def write(self, dataStr, sep=None):
+    def write(self, dataStr, sep=None, alert = True):
         """write the dataStr to file then splits it using sep to plot"""
         dat = dataStr.split(sep)
         if len(dat) != self.n_ch: 
@@ -80,8 +89,9 @@ def write(self, dataStr, sep=None):
                 for h in self.hist:
                     del h
                     del FileUpdater.canvas
-        
+        self.entries += 1
         if (self.entries % self.update_on == 0):
+            # print "%i complete, updating" % (self.entries/self.update_on)
             self.update()
 
     def update(self):
@@ -95,6 +105,14 @@ def update(self):
 
         FileUpdater.canvas.Update()
 
+    def closeFooter(self, footer):
+        """
+        writes a footer then closes the file"""
+        if footer:
+            if footer[-1] != '\n': footer = footer + '\n'
+            file.write(self, footer)
+        self.close()
+    
     def close(self):
         """close the file and delete the histograms"""
         FileUpdater.canvas = 0
@@ -140,7 +158,7 @@ def string_stripper(s):
 
     n_files = 2
     for i in range (n_files):
-        t = FileUpdater("test%i.txt")
+        t = FileUpdater("test%i.txt",header='hello')
 
         for i in range(n_tests):
             print "0"
@@ -162,4 +180,4 @@ def string_stripper(s):
     except Exception, e:    
         pass
     finally:
-        t.close()
+        t.close("hello chicken!")