diff --git a/script/ImageStats.py b/script/ImageStats.py
index b4d2ccb..a357e61 100644
--- a/script/ImageStats.py
+++ b/script/ImageStats.py
@@ -10,7 +10,7 @@ def get_centroid(source):
         return None
     op = show_panel(bi, "Original")
     ip = load_image(bi)  
-    plot(get_histogram(ip))  
+    plot(get_histogram(ip), title = "Histogram")  
     grayscale(ip)    
     invert(ip)
     gaussian_blur(ip)    
@@ -18,11 +18,10 @@ def get_centroid(source):
     
     #binary_erode(ip)        
     show_panel(ip.getBufferedImage(), "Image")
-    (results,output_img)=analyse_particles(ip, 1000,20000, print_table=True)    
+    (results,output_img)=analyse_particles(ip, 2000,20000, exclude_edges=False, print_table=True)    
     op.clearOverlays()
     show_panel(output_img.getBufferedImage(), "Outlines")
-    if results.size()>0:
-        print results        
+    if results.size()>0:        
         centroid  = (results.getValue("XM",0), results.getValue("YM",0))
         ov = Crosshairs(Pen(java.awt.Color.ORANGE),  java.awt.Point(int(centroid[0]),int(centroid[1])), java.awt.Dimension(15,15))
         op.addOverlay(ov)
@@ -39,23 +38,22 @@ class SimulatedSource(Filter):
         if img is None: 
             return None 
         ip = load_image(img)         
-        pad_h = int((random.random()-0.5) * 300)
-        pad_v = int((random.random()-0.5) * 300)
-        print "Pad = " , (pad_h, pad_v)
+        pad_h = int((random.random()-0.5) * 500)
+        pad_v = int((random.random()-0.5) * 500)
+        #print "Pad = " , (pad_h, pad_v)
         ip = pad_image(ip, -pad_h, pad_h, pad_v, -pad_v, fill_color = Color.BLACK)
         return ip.getBufferedImage() 
         #return img
         
     def waitNext(self, timeout):
-        #self.pushImage(self.process(self.img,None))
         self.pushImage(self.process(self.img, None))
-        print "---> " , self.getImage()
         
         
 simulated_source = SimulatedSource()           
 simulated_source.img=None
 simulated_source.initialize()
 image.addListener(simulated_source)
+show_panel(simulated_source)
 
        
 class ImageStats(DeviceBase):
@@ -115,9 +113,9 @@ class ImageStats(DeviceBase):
     def doClose(self):        
         pass    
 
+"""
 if __name__ == "__builtin__":
-    print get_centroid(simulated_source)    
-    show_panel(simulated_source)
+    print get_centroid(simulated_source)        
     
     add_device(ImageStats("image_stats", simulated_source), True)
     image_stats.enableBackground(False)    
@@ -126,3 +124,4 @@ if __name__ == "__builtin__":
         print image_stats.take(), image_stats.com_x_mean.read(), image_stats.com_y_mean.read() 
         time.sleep(1)
 
+"""
\ No newline at end of file