Ignore:
Timestamp:
05/10/05 12:13:12 (20 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mcalib/MMcCalibrationCalc.cc

    r4710 r7005  
    242242{
    243243    //
    244     // Exclude events with some saturated pixel
     244    // Exclude events with some high-gain saturated pixel
    245245    //
    246246    if (fPar->GetNumSatPixelsHG()>0)
     
    266266    //
    267267 
    268     const Float_t inner_photel = (Float_t) fMcEvt->GetPhotElfromShower() * innersize / size;
    269     const Float_t outer_photel = (Float_t) fMcEvt->GetPhotElfromShower() - inner_photel;
    270 
    271     const Float_t corrected_photel = inner_photel + outer_photel / fOuterPixelsLightCollection;
    272 
    273 
    274     //    fHistADC2PhotEl->Fill(TMath::Log10(fMcEvt->GetPhotElfromShower()/size));
     268    const Float_t corrected_photel = (Float_t) fMcEvt->GetPhotElfromShower() /
     269      (fOuterPixelsLightCollection + innersize / size * (1. - fOuterPixelsLightCollection));
    275270
    276271    fHistADC2PhotEl->Fill(TMath::Log10(corrected_photel/size));
Note: See TracChangeset for help on using the changeset viewer.