Changeset 3636 for trunk/MagicSoft/Mars/mcalib/MHGausEvents.h

Timestamp:

04/03/04 17:27:50 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MHGausEvents.h (modified) (9 diffs)

trunk/MagicSoft/Mars/mcalib/MHGausEvents.h

@@ -20 +20 @@
 private:
 
-  const static Float_t  fgProbLimit;            // Default for fProbLimit (now 0.001)
-  const static Int_t    fgNDFLimit;             // Default for fNDFLimit  (now 2)
-  const static Int_t    fgPowerProbabilityBins; // Default for fPowerProbabilityBins (now 20)
-  const static Int_t    fgBinsAfterStripping;   // Default for fBinsAfterStripping (now 40)
-
+  const static Float_t  fgProbLimit;            //! Default for fProbLimit (now set to: 0.001)
+  const static Int_t    fgNDFLimit;             //! Default for fNDFLimit  (now set to: 2)
+  const static Float_t  fgPickupLimit;          //! Default for fPickupLimit (now set to: 5. )
+  const static Int_t    fgPowerProbabilityBins; //! Default for fPowerProbabilityBins (now set to: 20)
+  const static Int_t    fgBinsAfterStripping;   //! Default for fBinsAfterStripping (now set to: 40)
+  
   Int_t    fPowerProbabilityBins;      // Bins for the projected power spectrum
   Int_t    fBinsAfterStripping;        // Bins for the Gauss Histogram after stripping off the zeros at both ends
@@ -41 +42 @@
   Double_t fProb;                      //  Probability of the Gauss fit (derived from Chi-Square and NDF
 
-  enum { kGausFitOK, kExpFitOK, kFourierSpectrumOK }; // Bits to hold information about fit results 
+  enum { kGausFitOK, kExpFitOK, kFourierSpectrumOK,
+         kExcluded };                  //  Bits to hold information about fit results 
   
-  Byte_t fFlags;                       //   Byte to hold the bits fit result bits
+  Byte_t fFlags;                       //  Byte to hold the bits fit result bits
   
-  Int_t fCurrentSize;                  //   Current size of the array fEvents
-  
+  Int_t fCurrentSize;                  //  Current size of the array fEvents
+
 protected:
 
+  Int_t   fNbins;                      // Number histogram bins for fHGausHist (used by Init())
+  Axis_t  fFirst;                      // Lower histogram edge  for fHGausHist (used by Init()) 
+  Axis_t  fLast;                       // Upper histogram edge  for fHGausHist (used by Init()) 
+  Int_t   fPixId;                      //  Pixel ID 
+  
   TH1F    fHGausHist;                  // Histogram which should hold the Gaussian distribution
   TArrayF fEvents;                     // Array which holds the entries of GausHist
@@ -57 +64 @@
   Float_t  fProbLimit;                 // Probability limit for judgement if fit is OK 
   Int_t    fNDFLimit;                  // NDF limit for judgement if fit is OK
+  Float_t  fPickupLimit;               // Upper number of sigmas from the mean until events are considered as pickup
 
   Float_t *CreateEventXaxis(Int_t n);  //   Create an x-axis for the Event TGraphs
@@ -75 +83 @@
   virtual void Clear(Option_t *o="");
   virtual void Reset();  
-
+  virtual void InitBins();
+  
   // Setters
   void  SetEventFrequency(const Float_t f)   { fEventFrequency = f; }
@@ -84 +93 @@
   void  SetSigmaErr( const Double_t d )   { fSigmaErr = d;   }
   void  SetProb    ( const Double_t d )   { fProb     = d;   }
+  void  SetPixId    ( const Int_t i    )   { fPixId    = i;   }
+
+  void  SetNbins    ( const Int_t i    )   { fNbins    = i;   }  
+  void  SetFirst   ( const Double_t d )   { fFirst    = d;   }
+  void  SetLast    ( const Double_t d )   { fLast     = d;   }
   
+  void  SetNDFLimit(  const Int_t   lim=fgNDFLimit  ) {  fNDFLimit = lim;  }  
+  void  SetPickupLimit( const Float_t  lim =fgPickupLimit)  { fPickupLimit  = lim;   }
   void  SetProbLimit( const Float_t lim=fgProbLimit ) {  fProbLimit = lim; }
-  void  SetNDFLimit(  const Int_t   lim=fgNDFLimit  ) {  fNDFLimit = lim;  }  
 
   // Setters ONLY for MC:
-  void  SetGausFitOK(         const Bool_t b );
-  void  SetExpFitOK(          const Bool_t b );
-  void  SetFourierSpectrumOK( const Bool_t b );
+  void  SetGausFitOK(         const Bool_t b=kTRUE );
+  void  SetExpFitOK(          const Bool_t b=kTRUE );
+  void  SetFourierSpectrumOK( const Bool_t b=kTRUE );
+  void  SetExcluded         ( const Bool_t b=kTRUE );  
 
   // Getters
@@ -101 +117 @@
   const Double_t GetProb()          const { return fProb;      }
   const Int_t    GetNdf()           const;
+  const Double_t GetPickup()        const;
+  const Int_t    GetPixId()         const { return fPixId;     }
 
   const Double_t GetSlope()        const;
@@ -136 +154 @@
   const Bool_t IsEmpty()              const;
   const Bool_t IsFourierSpectrumOK()  const;
+  const Bool_t IsExcluded         ()  const;
 
   // Fill
@@ -145 +164 @@
   Bool_t FitGaus(Option_t *option="RQ0",
                  const Double_t xmin=0., const Double_t xmax=0.); // Fit the histogram HGausHist with a Gaussian
-  void BypassFit();   // Take mean and RMS from the histogram
-  
+  Bool_t RepeatFit(const Option_t *option="RQ0");      // Repeat fit within limits defined by fPickupLimit
+  void BypassFit();                                    // Take mean and RMS from the histogram
   
   // Draws
@@ -155 +174 @@
   
   // Miscelleaneous
+  virtual void ChangeHistId(Int_t id);   // Changes names and titles of the histogram
   void CreateFourierSpectrum();             // Create the fourier spectrum out of fEvents
   void CreateGraphEvents();                 // Create the TGraph fGraphEvents of fEvents