Changeset 14487

Timestamp:

10/17/12 16:54:03 (12 years ago)

Author:

neise

Message:

debugged pointer f**k in implementation of despiker in calfactfits

File:

• fact/tools/pyscripts/pyfact/calfactfits.h (modified) (3 diffs)

fact/tools/pyscripts/pyfact/calfactfits.h

@@ -225 +225 @@
                         return false;
                 }
-                else {
+                else 
+        {
                         UInt_t drs_calib_offset;
             double raw, raw_bsl, shifted;
@@ -241 +242 @@
 
                 // stolen from TBs MARS ... DrsCalib.h or something like that
-                double *p = npcaldata;
+                 double *p = npcaldata+(pixel*data_nroi);
+
                 // second loop for despiking
                 for (size_t i=1; i<data_nroi-2; i++)
@@ -247 +249 @@
                     if (p[i]-p[i-1]>25 && p[i]-p[i+1]>25)
                     {
+                        //std::cout << "single cand:\t" << pixel << "\t" << i << std::endl;
                         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)
                     {
+                        //std::cout << "double cand:\t" << pixel << "\t" << i << std::endl;
                         p[i] = (p[i-1]+p[i+2])/2;
                         p[i+1] = p[i];