From 82efc26628859461e72b4fa314e1f7d50719daef Mon Sep 17 00:00:00 2001
From: voulot_d <voulot_d@sf-lc6a.psi.ch>
Date: Thu, 24 Aug 2017 08:27:54 +0200
Subject: [PATCH] Startup

---
 devices/CurrentCamera.properties       | 22 +++++-----
 script/Diagnostics/ICT_timing.py       | 59 ++++++++++++++++++++++++++
 script/Diagnostics/ICT_timing_rough.py | 39 +++++++++++++++++
 3 files changed, 109 insertions(+), 11 deletions(-)
 create mode 100644 script/Diagnostics/ICT_timing.py
 create mode 100644 script/Diagnostics/ICT_timing_rough.py

diff --git a/devices/CurrentCamera.properties b/devices/CurrentCamera.properties
index 966ef9b..5e8adb8 100644
--- a/devices/CurrentCamera.properties
+++ b/devices/CurrentCamera.properties
@@ -1,16 +1,16 @@
-#Wed Aug 23 13:29:11 CEST 2017
+#Wed Aug 23 19:15:45 CEST 2017
 colormap=Flame
 colormapAutomatic=false
-colormapMax=4718.0
-colormapMin=62.0
+colormapMax=1624.0
+colormapMin=67.0
 flipHorizontally=false
 flipVertically=false
 grayscale=false
-imageHeight=2148
-imageWidth=2272
+imageHeight=1680
+imageWidth=1744
 invert=false
-regionStartX=129
-regionStartY=7
+regionStartX=433
+regionStartY=241
 rescaleFactor=1.0
 rescaleOffset=0.0
 roiHeight=-1
@@ -21,9 +21,9 @@ rotation=0.0
 rotationCrop=false
 scale=1.0
 serverURL=localhost\:10000
-spatialCalOffsetX=-1057.4895329398094
-spatialCalOffsetY=-1420.5549062527236
-spatialCalScaleX=-8.510638153514359
-spatialCalScaleY=-8.235817137431614
+spatialCalOffsetX=-213.4660832767771
+spatialCalOffsetY=-193.62057589115275
+spatialCalScaleX=-8.512333108246718
+spatialCalScaleY=-8.592910815665192
 spatialCalUnits=mm
 transpose=false
diff --git a/script/Diagnostics/ICT_timing.py b/script/Diagnostics/ICT_timing.py
new file mode 100644
index 0000000..e323c77
--- /dev/null
+++ b/script/Diagnostics/ICT_timing.py
@@ -0,0 +1,59 @@
+import  ch.psi.pshell.epics.ChannelDouble as ChannelDouble
+
+run("CPython/wrapper")
+
+#t = ChannelDouble("delay (us)", "SINEG01-DICT009-EVR0:Pul0-Delay-SP")
+#t = ChannelDouble("fine delay (ns)", "SINEG01-DICT009-EVR0:UnivDlyModule0-Delay0-SP")
+#q = ChannelDouble("charge (pC)", "SINEG01-DICT215:B1_CHARGE")
+#V = ChannelDouble("DC output (V)", "SINEG01-DICT215:B1_VOLTAGE")
+#t = ChannelDouble("delay (us)", "S10DI01-DICT113-EVR0:Pul0-Delay-SP")
+#t = ChannelDouble("fine delay (ns)", "S10DI01-DICT113-EVR0:UnivDlyModule0-Delay0-SP")
+#q = ChannelDouble("charge (pC)", "S10DI01-DICT025:B1_CHARGE")
+#V = ChannelDouble("DC output (V)", "S10DI01-DICT025:B1_VOLTAGE")
+#t = ChannelDouble("delay (us)", "SARMA01-DICT482-EVR0:Pul0-Delay-SP")
+#t = ChannelDouble("fine delay (ns)", "SARMA01-DICT482-EVR0:UnivDlyModule0-Delay0-SP")
+#q = ChannelDouble("charge (pC)", "SARMA01-DICT090:B1_CHARGE")
+#t = ChannelDouble("delay (us)", "SARBD01-DICT599-EVR0:Pul0-Delay-SP")
+t = ChannelDouble("fine delay (ns)", "SARBD01-DICT599-EVR0:UnivDlyModule0-Delay0-SP")
+q = ChannelDouble("charge (pC)", "SARBD01-DICT030:B1_CHARGE")
+V = ChannelDouble("DC output (V)", "SARBD01-DICT030:B1_VOLTAGE")
+t.initialize()
+q.initialize()
+V.initialize()
+
+t0 = t.read()
+
+#start   = 1003.000
+#stop    = 1003.080
+#step    = 0.001
+#nb      = 1
+#lat     = 0.300
+start   = 2.2
+stop    = 12.7
+step    = 0.1
+nb      = 5
+lat     = 0.150
+
+
+try:
+    q_averager = create_averager(q, nb, lat)
+    V_averager = create_averager(V, nb, lat)
+    V_averager.monitored=True
+    r = lscan(t, (q_averager, V_averager), start, stop, step, latency=lat)
+    delay     = r.getPositions(0)
+    charge    = [val.mean for val in r.getReadable(0)]
+    chargerms = [val.stdev for val in r.getReadable(0)]
+    volt      = [val.mean for val in r.getReadable(1)]
+    voltrms   = [val.stdev for val in r.getReadable(1)]
+finally:
+    t.write(t0)
+    t.close()
+    q.close()
+    V.close()
+try:
+    i_max = charge.index(max(charge))
+    min_i, max_i = max(i_max-10, 0), min(i_max+11, len(charge))
+    (delay_charge_max, charge_max, delay_charge_fit, charge_fit)  = extremum(delay[min_i:max_i], charge[min_i:max_i])
+    print delay_charge_max, charge_max
+except:
+    raise Exception("Fit failure")
diff --git a/script/Diagnostics/ICT_timing_rough.py b/script/Diagnostics/ICT_timing_rough.py
new file mode 100644
index 0000000..9a83159
--- /dev/null
+++ b/script/Diagnostics/ICT_timing_rough.py
@@ -0,0 +1,39 @@
+import  ch.psi.pshell.epics.ChannelDouble as ChannelDouble
+
+dev_name = "SINEG01-DICT009"
+base_name = "SINEG01-DICT215"
+
+#t = ChannelDouble("delay (us)", "SINEG01-DICT009-EVR0:Pul0-Delay-SP")
+#q = ChannelDouble("charge (pC)", "SINEG01-DICT215:B1_CHARGE")
+#V = ChannelDouble("DC output (V)", "SINEG01-DICT215:B1_VOLTAGE")
+t = ChannelDouble("delay (us)", "S10DI01-DICT113-EVR0:Pul0-Delay-SP")
+q = ChannelDouble("charge (pC)", "S10DI01-DICT025:B1_CHARGE")
+V = ChannelDouble("DC output (V)", "S10DI01-DICT025:B1_VOLTAGE")
+t.initialize()
+q.initialize()
+V.initialize()
+
+t0 = t.read()
+
+start   = 1001.340
+stop    = 1001.400
+step    = 0.001
+nb      = 1
+lat     = 0.300
+
+
+try:
+    q_averager = create_averager(q, nb, lat)
+    V_averager = create_averager(V, nb, lat)
+    V_averager.monitored=True
+    r = lscan(t, (q_averager, V_averager), start, stop, step, latency=lat)
+    delay     = r.getPositions(0)
+    charge    = [val.mean for val in r.getReadable(0)]
+    chargerms = [val.stdev for val in r.getReadable(0)]
+    volt      = [val.mean for val in r.getReadable(1)]
+    voltrms   = [val.stdev for val in r.getReadable(1)]
+finally:
+    t.write(t0)
+    t.close()
+    q.close()
+    V.close()