diff --git a/script/ScriptingExample.py b/script/ScriptingExample.py
index 2df47ffa..9be60e53 100644
--- a/script/ScriptingExample.py
+++ b/script/ScriptingExample.py
@@ -20,8 +20,6 @@ class PseudoDevice(Writable):
 pseudo_device = PseudoDevice()
 #add_device(PseudoDevice("pseudo_device"), True)
 
-        
-
 
 
 ascan((ManipulatorPhi, ManipulatorTheta, pseudo_device), (Counts, Scienta.image), (9.5,7.5,50.0), (10.5,8.5,70.0), (5,5,1), zigzag = True, before_read=before_readout, after_read = after_readout)
diff --git a/script/local.py b/script/local.py
index 04a3d952..a60b142c 100644
--- a/script/local.py
+++ b/script/local.py
@@ -3,6 +3,8 @@ import ch.psi.pshell.device.Readable.ReadableArray as ReadableArray
 import ch.psi.pshell.device.Readable.ReadableCalibratedArray as ReadableCalibratedArray
 import ch.psi.pshell.device.ArrayCalibration as ArrayCalibration
 import ch.psi.utils.Str
+from mathutils import estimate_peak_indexes, fit_gaussians, create_fit_point_list, Gaussian
+import java.awt.Color as Color
 
 
 beam_ok = True
@@ -259,4 +261,83 @@ controller.dataManager.setAlias(Scienta.spectrum, "ScientaSpectrum")
 
 
 #Additional device configuration
-ManipulatorPhi.trustedWrite = False
\ No newline at end of file
+ManipulatorPhi.trustedWrite = False
+
+
+def fit(ydata, xdata = None):
+    """
+    """
+    if xdata is None:
+        xdata = frange(0, len(ydata), 1)
+    max_y= max(ydata)   
+    index_max = ydata.index(max_y)   
+    max_x= xdata[index_max]
+    print "Max index:" + str(index_max),
+    print " x:" + str(max_x),
+    print " y:" + str(max_y)
+    gaussians = fit_gaussians(ydata, xdata, [index_max,])
+    (norm, mean, sigma) = gaussians[0]
+    p = plot([ydata],["data"],[xdata], title="Fit" )[0]
+    fitted_gaussian_function = Gaussian(norm, mean, sigma)   
+    scale_x = [float(min(xdata)), float(max(xdata)) ]
+    points = max((len(xdata)+1), 100)
+    resolution = (scale_x[1]-scale_x[0]) / points   
+    fit_y = []
+    fit_x = frange(scale_x[0],scale_x[1],resolution, True)
+    for x in fit_x:
+        fit_y.append(fitted_gaussian_function.value(x))   
+    p.addSeries(LinePlotSeries("fit"))
+    p.getSeries(1).setData(fit_x, fit_y)   
+    
+    if abs(mean - xdata[index_max]) < ((scale_x[0] + scale_x[1])/2):
+        print "Mean -> " +  str(mean)
+        p.addMarker(mean, None, "Mean="+str(round(norm,2)), Color.MAGENTA.darker())
+        return (norm, mean, sigma)
+    else:
+        p.addMarker(max_x, None, "Max="+str(round(max_x,2)), Color.GRAY)
+        print "Invalid gaussian fit: " +  str(mean)
+        return (None, None, None)
+
+
+
+def elog(title, message, attachments = [], author = None, category = "Info", domain = "", logbook = "SwissFEL commissioning data", encoding=1):
+    """
+    Add entry to ELOG.
+    """
+    if author is None:
+        author = "pshell" #controller.getUser().name
+    typ = "pshell"
+    entry = ""
+   
+    cmd = 'G_CS_ELOG_add -l "' + logbook+ '" '
+    cmd = cmd + '-a "Author=' + author + '" '
+    cmd = cmd + '-a "Type=' + typ + '" '
+    cmd = cmd + '-a "Entry=' + entry + '" '
+    cmd = cmd + '-a "Title=' + title + '" '
+    cmd = cmd + '-a "Category=' + category + '" '
+    cmd = cmd + '-a "Domain=' + domain + '" '
+    for attachment in attachments:
+        cmd = cmd + '-f "' + attachment + '" '
+    cmd = cmd + '-n ' + str(encoding)
+    cmd = cmd + ' "' + message + '"'
+    #print cmd
+    #os.system (cmd)
+    #print os.popen(cmd).read()  
+    import subprocess
+    proc = subprocess.Popen(cmd, stdout=subprocess.PIPE, shell=True)
+    (out, err) = proc.communicate()
+    if (err is not None) and err!="":
+        raise Exception(err)
+    print out
+
+def get_plot_snapshots(title = None, file_type = "jpg", temp_path = controller.setup.getContextPath()):
+    """
+    Returns list with file names of plots snapshots from a plotting context.
+    """
+    sleep(0.02) #Give some time to plot to be finished - it is not sync  with acquisition
+    ret = []
+    for p in get_plots(title):        
+        file_name = os.path.abspath(temp_path + "/" + p.getTitle() + "." + file_type)
+        p.saveSnapshot(file_name , file_type)
+        ret.append(file_name)
+    return ret