Index: fact/tools/pyscripts/pyfact/calfactfits.h
===================================================================
--- fact/tools/pyscripts/pyfact/calfactfits.h	(revision 14481)
+++ fact/tools/pyscripts/pyfact/calfactfits.h	(revision 14482)
@@ -239,4 +239,21 @@
 					//Note: data_nroi=calib_nroi, calib_blm_size=calib_gm_size
 				}
+
+                // stolen from TBs MARS ... DrsCalib.h or something like that
+                double *p = npcaldata;
+                // second loop for despiking
+                for (size_t i=1; i<data_nroi-2; i++)
+                {
+                    if (p[i]-p[i-1]>25 && p[i]-p[i+1]>25)
+                    {
+                        p[i] = (p[i-1]+p[i+1])/2;
+                    }
+                    if (p[i]-p[i-1]>22 && fabs(p[i]-p[i+1])<4 && p[i+1]-p[i+2]>22)
+                    {
+                        p[i] = (p[i-1]+p[i+2])/2;
+                        p[i+1] = p[i];
+                    }
+                }
+                
 			}
 		}