Changeset 2852 for trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.h

Timestamp:

01/19/04 23:02:57 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.h (modified) (10 diffs)

trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.h

@@ -38 +38 @@
   TH1I* fHBlindPixelChargevsN;     //-> Summed Charge vs. Event Nr. 
   
-  TF1 *fSinglePheFit;
-  TF1 *fTimeGausFit;  
-  TF1 *fSinglePhePedFit;
+  TF1 fSinglePheFit;
+  TF1 fTimeGausFit;  
+  TF1 fSinglePhePedFit;
 
   Axis_t  fBlindPixelChargefirst;
@@ -49 +49 @@
   void DrawLegend();
 
-  TPaveText *fFitLegend;
+  TPaveText *fFitLegend;                  //!
   Bool_t fFitOK;  
   
@@ -81 +81 @@
   ~MHCalibrationBlindPixel();
 
-  typedef Double_t (*BlindPixelFitFunc)(Double_t *, Double_t *);
-
   Bool_t FillBlindPixelCharge(Float_t q)             { return fHBlindPixelCharge->Fill(q) > -1;  }  
   Bool_t FillBlindPixelTime(Int_t t)                 { return fHBlindPixelTime->Fill(t) > -1;  }
   Bool_t FillBlindPixelChargevsN(Stat_t rq, Int_t t) { return fHBlindPixelChargevsN->Fill(t,rq) > -1;  }  
   
+
+  //Getters
   const Double_t GetLambda()         const { return fLambda; }
-  const Double_t GetLambdaCheck()    const { return fLambdaCheck; }
+  const Double_t GetLambdaCheck()     const { return fLambdaCheck; }
   const Double_t GetMu0()            const { return fMu0; }
   const Double_t GetMu1()            const { return fMu1; }
@@ -95 +95 @@
 
   const Double_t GetLambdaErr()      const { return fLambdaErr; }
-  const Double_t GetLambdaCheckErr() const { return fLambdaCheckErr; }
+  const Double_t GetLambdaCheckErr()  const { return fLambdaCheckErr; }
   const Double_t GetMu0Err()         const { return fMu0Err; }
   const Double_t GetMu1Err()         const { return fMu1Err; }
@@ -102 +102 @@
 
   const Double_t GetChiSquare()      const { return fChisquare; }
-  const Double_t GetProb()           const { return fProb;      }  
-  const Int_t    GetNdf()            const { return fNdf;       }   
+  const Double_t GetProb()         const { return fProb;      }  
+  const Int_t    GetNdf()          const { return fNdf;       }   
 
   const Double_t GetMeanTime()      const { return fMeanTime; }
@@ -110 +110 @@
   const Double_t GetSigmaTimeErr()    const { return fSigmaTimeErr; }
 
-  Bool_t SimulateSinglePhe(Double_t lambda,
-                           Double_t mu0,
-                           Double_t mu1,
-                           Double_t sigma0,
-                           Double_t sigma1);
-  
+  const Bool_t IsFitOK()                { return fFitOK; }
+
+  // Draws
+  TObject *DrawClone(Option_t *option="") const;
+  void Draw(Option_t *option="");
+
+  // Fits
+  enum FitFunc_t  { kEPoisson4, kEPoisson5, kEPoisson6, kEPoisson7, kEPolya, kEMichele };
+
+private:
+  FitFunc_t fFitFunc;
+
+public:
   Bool_t FitSinglePhe(Axis_t rmin=0, Axis_t rmax=0, Option_t *opt="RL0+Q");
   Bool_t FitTime(Axis_t rmin=0., Axis_t rmax=0.,Option_t *opt="R0+Q");
-
-  void ChangeFitFunc(BlindPixelFitFunc fitfunc, Int_t par=6);
-
+  void ChangeFitFunc(FitFunc_t func)      { fFitFunc = func;  }
+  
+  // Simulation
+  Bool_t SimulateSinglePhe(Double_t lambda,
+                           Double_t mu0,Double_t mu1,
+                           Double_t sigma0,Double_t sigma1);
+  
+  // Others
   void CutAllEdges();
 
-  TObject *DrawClone(Option_t *option="") const;
-  void Draw(Option_t *option="");
-
-  Bool_t IsFitOK() { return fFitOK; }
-
 private:
-  
-  BlindPixelFitFunc fgSinglePheFitFunc;     //! In the beginning, 
-  Int_t             fgSinglePheFitNPar;     //! we want to be flexible using different functions
-
-  inline static Double_t fAna(Double_t *x, Double_t *par)
+
+  void InitFit(TF1& f, Axis_t min, Axis_t max);
+  void ExitFit(TF1& f);  
+  
+  inline static Double_t fFitFuncMichele(Double_t *x, Double_t *par)
     {
 
@@ -189 +196 @@
     }
     
-  inline static Double_t fKto4(Double_t *x, Double_t *par)
+  inline static Double_t fPoissonKto4(Double_t *x, Double_t *par)
     {
 
@@ -246 +253 @@
 
   
-  inline static Double_t fKto5(Double_t *x, Double_t *par)
+  inline static Double_t fPoissonKto5(Double_t *x, Double_t *par)
     {
       
@@ -311 +318 @@
   
   
-  inline static Double_t fKto6(Double_t *x, Double_t *par)
+  inline static Double_t fPoissonKto6(Double_t *x, Double_t *par)
     {
       
@@ -382 +389 @@
     }
   
-  ClassDef(MHCalibrationBlindPixel, 1) 
+  ClassDef(MHCalibrationBlindPixel, 1)  // Histograms from the Calibration Blind Pixel
 };