Changeset 2972


Ignore:
Timestamp:
01/30/04 10:19:25 (21 years ago)
Author:
gaug
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/Changelog

    r2971 r2972  
    2121     - fill Blind Pixel with information about its fitted pedestal if
    2222       available
     23     - initialize charges histogram with a bigger range
     24     - use pedestal information to constrain the pedestal fit ranges
    2325
    2426
  • trunk/MagicSoft/Mars/mcalib/MCalibrationCalc.cc

    r2956 r2972  
    114114using namespace std;
    115115
     116const Int_t MCalibrationCalc::fBlindPixelId = 559;
     117const Int_t MCalibrationCalc::fPINDiodeId   = 9999;
     118
    116119// --------------------------------------------------------------------------
    117120//
     
    121124    : fPedestals(NULL), fCalibrations(NULL), fSignals(NULL),
    122125      fRawEvt(NULL), fRunHeader(NULL), fArrivalTime(NULL), fEvtTime(NULL)
    123   //      fBlindPixelId(559), fPINDiodeId(9999)
    124126{
    125127
     
    461463        {
    462464         
    463         case kBlindPixelId:
     465        case fBlindPixelId:
    464466         
    465467          if (!blindpixel.FillCharge(sumhi))
     
    479481          break;
    480482         
    481         case kPINDiodeId:
     483        case fPINDiodeId:
    482484
    483485          if (!pindiode.FillCharge(sumhi))
     
    583585      *fLog << inf << GetDescriptor() << ": Fitting the Blind Pixel" << endl;
    584586
     587      //
     588      // retrieve mean and sigma of the blind pixel pedestal,
     589      // so that we can use it for the fit
     590      //
     591      if (fPedestals->IsUseHists())
     592        {
     593          //
     594          // retrieve the pedestal pix of the blind pixel
     595          //
     596          MPedestalPix &ped = (*fPedestals)[fBlindPixelId];
     597          //
     598          // retrieve the blind pixel histogram container
     599          //
     600          MHCalibrationBlindPixel *hist = blindpixel.GetHist();
     601          //
     602          // Set the corresponding values
     603          //
     604          const Float_t peddiff     = ped.GetMean()
     605                                    - ped.GetPedestal()*fSignals->GetNumUsedFADCSlices();
     606          const Float_t pederr      = ped.GetMeanErr();
     607          const Float_t pedsigma    = ped.GetSigma();
     608          const Float_t pedsigmaerr = ped.GetSigmaErr();
     609         
     610          hist->SetMeanPedestal(peddiff);
     611          hist->SetMeanPedestalErr(pederr);
     612          hist->SetSigmaPedestal(pedsigma);
     613          hist->SetSigmaPedestalErr(pedsigmaerr);
     614        }
     615     
    585616      if (!blindpixel.FitCharge())
    586617        {
  • trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.h

    r2931 r2972  
    6464  Double_t  fLambdaCheckErr;
    6565
     66  Double_t  fMeanPedestal;
     67  Double_t  fMeanPedestalErr;
     68  Double_t  fSigmaPedestal;
     69  Double_t  fSigmaPedestalErr;
     70 
    6671public:
    6772
     
    7681  Bool_t FillBlindPixelChargevsN(Stat_t rq, Int_t t);
    7782 
    78 
    79   //Getters
     83  // Setters
     84  void SetMeanPedestal(const Float_t f)      { fMeanPedestal = f;     }
     85  void SetMeanPedestalErr(const Float_t f)   { fMeanPedestalErr = f;  }
     86  void SetSigmaPedestal(const Float_t f)     { fSigmaPedestal = f;    }
     87  void SetSigmaPedestalErr(const Float_t f)  { fSigmaPedestalErr = f; }
     88 
     89  // Getters
    8090  const Double_t GetLambda()         const { return fLambda;         }
    8191  const Double_t GetLambdaCheck()    const { return fLambdaCheck;    }
     
    224234      Double_t arg = 0.;
    225235     
    226       //      Double_t mu0 = 0.;
    227236      Double_t mu0 = par[1];
    228237      Double_t mu1 = par[2];
Note: See TracChangeset for help on using the changeset viewer.