Ignore:
Timestamp:
10/17/06 18:18:40 (18 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

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

    r7188 r8106  
    11/* ======================================================================== *\
     2! $Name: not supported by cvs2svn $:$Id: MCalibrationIntensityChargeCam.cc,v 1.23 2006-10-17 17:15:59 tbretz Exp $
     3! --------------------------------------------------------------------------
    24!
    35! *
     
    266268      sigerr[cnt]    = pheerr;
    267269
    268       Double_t sig  = 0.;
     270      Double_t sig1 = 0.;
    269271      Double_t sig2 = 0.;
    270272      Int_t    num  = 0;
    271273
    272       for (Int_t i=0; i<cam->GetSize(); i++)
    273         {
    274           const MCalibrationChargePix &pix = (MCalibrationChargePix&)(*cam)[i];
     274      for (Int_t j=0; j<cam->GetSize(); j++)
     275        {
     276          const MCalibrationChargePix &pix = (MCalibrationChargePix&)(*cam)[j];
    275277          //
    276278          // Don't use bad pixels
     
    280282          //
    281283          //
    282           if (aidx != geom[i].GetAidx())
     284          if (aidx != geom[j].GetAidx())
    283285            continue;
    284286         
    285           sig  += pix.GetConvertedMean();
     287          sig1 += pix.GetConvertedMean();
    286288          sig2 += pix.GetConvertedMean() * pix.GetConvertedMean();
    287289          num++;
     
    290292      if (num > 1)
    291293        {
    292           sig          /= num;
    293 
    294           Double_t var = (sig2 - sig*sig*num) / (num-1);
    295           var /= sig*sig;
     294          sig1 /= num;
     295
     296          Double_t var = (sig2 - sig1*sig1*num) / (num-1);
     297          var /= sig1*sig1;
    296298          var += pherelvar;
    297299
    298           phepersig[cnt] = phe/sig;
     300          phepersig[cnt] = phe/sig1;
    299301          if (var > 0.)
    300302            phepersigerr[cnt] = TMath::Sqrt(var) * phepersig[cnt];
     
    10571059      for (Int_t i=0; i<GetSize(); i++)
    10581060        {
    1059           MCalibrationChargeCam *cam = (MCalibrationChargeCam*)GetCam(i);
     1061          MCalibrationChargeCam *ccam = (MCalibrationChargeCam*)GetCam(i);
    10601062          //
    10611063          // Get the calibration pix from the calibration cam
    10621064          //
    1063           MCalibrationChargePix &pix = (MCalibrationChargePix&)(*cam)[npix];
     1065          MCalibrationChargePix &pix = (MCalibrationChargePix&)(*ccam)[npix];
    10641066          //
    10651067          // Don't use bad pixels
     
    10961098          if (option.Contains("conversionfactor"))
    10971099          {
    1098               const MCalibrationChargePix &apix = (MCalibrationChargePix&)cam->GetAverageArea(0);
     1100              const MCalibrationChargePix &apix = (MCalibrationChargePix&)ccam->GetAverageArea(0);
    10991101              pvar = apix.GetPheFFactorMethod()/pix.GetConvertedMean();
    11001102          }
Note: See TracChangeset for help on using the changeset viewer.