Ignore:
Timestamp:
07/01/03 10:46:07 (21 years ago)
Author:
wittek
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/manalysis/MSigmabar.cc

    r2237 r2255  
    136136
    137137        const Int_t idx = cerpix.GetPixId();
    138         const Double_t area = geom.GetPixRatio(idx);
    139 
    140138        if (idx == 0)
    141139        {
     
    145143        }
    146144
     145        // ratio is the area of pixel 0
     146        //          divided by the area of the current pixel
     147        const Double_t ratio = geom.GetPixRatio(idx);
     148
    147149        const MGeomPix &gpix = geom[idx];
    148150
     
    157159            continue;
    158160
    159         if (area < 1.5)
     161        //if (area < 1.5)
     162        if (ratio > 0.5)
    160163        {
    161164            innerPixels[sector]++;
    162             innerSquaredSum[sector]+= sigma*sigma / area;
     165            innerSquaredSum[sector]+= sigma*sigma * ratio;
    163166        }
    164167        else
    165168        {
    166169            outerPixels[sector]++;
    167             outerSquaredSum[sector]+= sigma*sigma / area;
     170            outerSquaredSum[sector]+= sigma*sigma * ratio;
    168171        }
    169172    }
     
    243246
    244247}
     248
     249
     250
     251
     252
     253
     254
     255
     256
Note: See TracChangeset for help on using the changeset viewer.