Changeset 14482 for fact/tools/pyscripts/pyfact

Timestamp:

10/17/12 14:01:07 (12 years ago)

Author:

neise

Message:

implemented spike destruction algo from TB into calfactfits class

File:

• fact/tools/pyscripts/pyfact/calfactfits.h (modified) (1 diff)

fact/tools/pyscripts/pyfact/calfactfits.h

@@ -239 +239 @@
                                         //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];
+                    }
+                }
+                
                         }
                 }