Ignore:
Timestamp:
10/17/06 18:18:40 (18 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/mcalib
Files:
4 edited

Legend:

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

    r5701 r8106  
    11/* ======================================================================== *\
     2! $Name: not supported by cvs2svn $:$Id: MCalibrationBlindPix.cc,v 1.15 2006-10-17 17:15:59 tbretz Exp $
     3! --------------------------------------------------------------------------
    24!
    35! *
     
    375377{
    376378
    377   const Float_t err = fAttErr[fColor];
    378 
    379   if (err < 0.)
    380     return -1.;
    381  
    382   return err*err*2.3;
     379  const Float_t ferr = fAttErr[fColor];
     380
     381  if (ferr < 0.)
     382    return -1.;
     383 
     384  return ferr*ferr*2.3;
    383385}
    384386
  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

    r7881 r8106  
    11/* ======================================================================== *\
     2! $Name: not supported by cvs2svn $:$Id: MCalibrationChargeCalc.cc,v 1.166 2006-10-17 17:15:59 tbretz Exp $
     3! --------------------------------------------------------------------------
    24!
    35! *
     
    18071809     
    18081810      MCalibrationChargePix &pix   = (MCalibrationChargePix&)(*chargecam)[i];
    1809       MCalibrationQEPix     &qepix = (MCalibrationQEPix&)    (*qecam)    [i];
     1811      MCalibrationQEPix     &qpix = (MCalibrationQEPix&)    (*qecam)    [i];
    18101812      MBadPixelsPix         &bad   =                         (*badcam)   [i];
    18111813
     
    18181820      const Float_t qerelvar = avphotrelvar +  pix.GetPheFFactorMethodRelVar();
    18191821
    1820       qepix.SetQEFFactor    ( qe            , fPulserColor );
    1821       qepix.SetQEFFactorVar ( qerelvar*qe*qe, fPulserColor );     
    1822       qepix.SetFFactorMethodValid(  kTRUE   , fPulserColor );
    1823 
    1824       if (!qepix.UpdateFFactorMethod( qecam->GetPlexiglassQE() ))
     1822      qpix.SetQEFFactor    ( qe            , fPulserColor );
     1823      qpix.SetQEFFactorVar ( qerelvar*qe*qe, fPulserColor );
     1824      qpix.SetFFactorMethodValid(  kTRUE   , fPulserColor );
     1825
     1826      if (!qpix.UpdateFFactorMethod( qecam->GetPlexiglassQE() ))
    18251827        *fLog << warn << GetDescriptor()
    18261828              << ": Cannot update Quantum efficiencies with the F-Factor Method" << endl;
     
    19491951    {
    19501952     
    1951       MCalibrationChargePix &pix   = (MCalibrationChargePix&)(*chargecam)[i];
    1952       MCalibrationQEPix     &qepix = (MCalibrationQEPix&)    (*qecam)    [i];
    1953       MBadPixelsPix         &bad   =                         (*badcam)   [i];
     1953      MCalibrationChargePix &pix  = (MCalibrationChargePix&)(*chargecam)[i];
     1954      MCalibrationQEPix     &qpix = (MCalibrationQEPix&)    (*qecam)    [i];
     1955      MBadPixelsPix         &bad  =                         (*badcam)   [i];
    19541956
    19551957      if (bad.IsUnsuitable(MBadPixelsPix::kUnsuitableRun))
    19561958        {
    1957           qepix.SetFFactorMethodValid(kFALSE,fPulserColor);
     1959          qpix.SetFFactorMethodValid(kFALSE,fPulserColor);
    19581960          pix.SetFFactorMethodValid(kFALSE);
    19591961          continue;
  • trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCam.cc

    r7059 r8106  
    11/* ======================================================================== *\
     2! $Name: not supported by cvs2svn $:$Id: MCalibrationChargeCam.cc,v 1.69 2006-10-17 17:18:40 tbretz Exp $
     3! --------------------------------------------------------------------------
    24!
    35! *
     
    680682
    681683
    682 Bool_t MCalibrationChargeCam::GetConversionFactorFFactor(Int_t ipx, Float_t &mean, Float_t &err, Float_t &ffactor)
     684Bool_t MCalibrationChargeCam::GetConversionFactorFFactor(Int_t ipx, Float_t &mean, Float_t &ferr, Float_t &ffactor)
    683685{
    684686
     
    691693
    692694  mean    = conv;
    693   err     = pix.GetMeanConvFADC2PheErr();
     695  ferr    = pix.GetMeanConvFADC2PheErr();
    694696  ffactor = pix.GetMeanFFactorFADC2Phot();
    695697
  • 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.