Ignore:
Timestamp:
03/12/04 16:07:58 (21 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

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

    r3471 r3487  
    153153        ULong_t sqsb  = 0;   // sum sqares signal+background
    154154
    155         Int_t sat = 0; // saturates?
    156         Int_t nb  = 0;
    157         Int_t nsb = 0;
     155        Int_t sat  = 0; // saturates?
     156        Int_t ishi = 0; // has a high content?
     157        Int_t nb   = 0;
     158        Int_t nsb  = 0;
    158159
    159160        if (*max==255)
    160161            sat++;
     162
     163        if (*max>80)
     164            ishi++;
    161165
    162166        while (ptr<end)
     
    177181        }
    178182
    179         if (sat>1)
     183        if (sat==0 && ishi)
     184        {
     185            // Area: x9
     186            ptr = pixel.GetLoGainSamples();
     187            max = ptr+pixel.GetIdxMaxLoGainSample();
     188            end = ptr+fRawEvt->GetNumLoGainSamples();
     189
     190            sumb = 0;   // sum background
     191            sqb  = 0;   // sum sqares background
     192            nb   = 0;
     193
     194            while (ptr<end)
     195            {
     196                // Background already calced from hi-gains!
     197                sumb += *ptr;
     198                sqb  += *ptr* *ptr;
     199                nb++;
     200            }
     201        }
     202
     203        if (sat>1 && !ishi)
    180204        {
    181205            // Area: x9
Note: See TracChangeset for help on using the changeset viewer.