Ignore:
Timestamp:
10/17/12 16:54:03 (12 years ago)
Author:
neise
Message:
debugged pointer f**k in implementation of despiker in calfactfits
File:
1 edited

Legend:

Unmodified
Added
Removed
  • fact/tools/pyscripts/pyfact/calfactfits.h

    r14482 r14487  
    225225                        return false;
    226226                }
    227                 else {
     227                else
     228        {
    228229                        UInt_t drs_calib_offset;
    229230            double raw, raw_bsl, shifted;
     
    241242
    242243                // stolen from TBs MARS ... DrsCalib.h or something like that
    243                 double *p = npcaldata;
     244                 double *p = npcaldata+(pixel*data_nroi);
     245
    244246                // second loop for despiking
    245247                for (size_t i=1; i<data_nroi-2; i++)
     
    247249                    if (p[i]-p[i-1]>25 && p[i]-p[i+1]>25)
    248250                    {
     251                        //std::cout << "single cand:\t" << pixel << "\t" << i << std::endl;
    249252                        p[i] = (p[i-1]+p[i+1])/2;
    250253                    }
    251254                    if (p[i]-p[i-1]>22 && fabs(p[i]-p[i+1])<4 && p[i+1]-p[i+2]>22)
    252255                    {
     256                        //std::cout << "double cand:\t" << pixel << "\t" << i << std::endl;
    253257                        p[i] = (p[i-1]+p[i+2])/2;
    254258                        p[i+1] = p[i];
Note: See TracChangeset for help on using the changeset viewer.