Changeset 6988 for trunk/MagicSoft/Mars/mhflux/MHAlpha.h

Timestamp:

05/02/05 10:17:41 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhflux/MHAlpha.h (modified) (4 diffs)

trunk/MagicSoft/Mars/mhflux/MHAlpha.h

@@ -24 +24 @@
 class MPointingPos;
 
-/*
-class MAlphaFitter : public MParContainer
-{
-private:
-    Float_t fSigInt;
-    Float_t fSigMax;
-    Float_t fBgMin;
-    Float_t fBgMax;
-    Int_t   fPolynom;
-
-    Double_t fSignificance;
-    Double_t fExcessEvents;
-    Double_t fChiSqSignal;
-    Double_t fChiSqBg;
-    Double_t fSigmaGaus;
-    Double_t fIntegralMax;
-
-public:
-    MAlphaFitter() : fSigInt(10), fSigMax(75), fBgMin(45), fBgMax(85), fPolynom(2)
-    {
-    }
-
-    MAlphaFitter(const MAlphaFitter &f)
-    {
-        f.Copy(*this);
-    }
-
-    void Copy(TObject &o) const
-    {
-        MAlphaFitter &f = static_cast<MAlphaFitter&>(o);
-
-        f.fSigInt  = fSigInt;
-        f.fSigMax  = fSigMax;
-        f.fBgMin   = fBgMin;
-        f.fBgMax   = fBgMax;
-        f.fPolynom = fPolynom;
-    }
-
-    void SetSignalIntegralMax(Float_t s) { fSigInt = s; }
-    void SetSignalFitMax(Float_t s)      { fSigMax = s; }
-    void SetBackgroundFitMin(Float_t s)  { fBgMin = s; }
-    void SetBackgroundFitMax(Float_t s)  { fBgMax = s; }
-    void SetNumPolynom(Int_t s)          { fPolynom = s; }
-
-    Double_t GetExcessEvents() const { return fExcessEvents; }
-    Double_t GetSignificance() const { return fSignificance; }
-    Double_t GetChiSqSignal() const  { return fChiSqSignal; }
-    Double_t GetChiSqBg() const      { return fChiSqBg; }
-    Double_t GetSigmaGaus() const    { return fSigmaGaus; }
-
-    void PaintResult(Float_t x=0.04, Float_t y=0.94, Float_t size=0.035) const;
-
-    Bool_t Fit(TH1D &h, Bool_t paint=kFALSE);
-    ClassDef(MAlphaFitter, 1)
-};
-*/
-
 class MHAlpha : public MH
 {
+protected:
+    TH3D    fHist;              // Alpha vs. theta and energy
+    TH1D    fHistTime;          //! temporary histogram to get alpha vs. time
+
+    TString fNameParameter;
+
+    MParContainer *fParameter;  //!
+
 private:
     const TH3D *fOffData;
@@ -88 +39 @@
     MAlphaFitter fFit;          // SEEMS THAT STREAMER HAS SOME PROBLEMS... MAYBE IF FUNC IS USED AT THE SAME TIME FOR FITS (PAINT)
 
-    TH3D fHAlpha;               // Alpha vs. theta and energy
     TH1D fHEnergy;              // excess events vs energy
     TH1D fHTheta;               // excess events vs theta
     TH1D fHTime;                // excess events vs time;
-    TH1D fHAlphaTime;           //! temporary histogram to get alpha vs. time
 
-    MParameterD  *fResult;      //!
-    MParameterD  *fEnergy;      //!
-    MHillas      *fHillas;      //!
-    MPointingPos *fPointPos;    //!
+    MParameterD   *fResult;     //!
+    MParameterD   *fEnergy;     //!
+    MHillas       *fHillas;     //!
+    MPointingPos  *fPointPos;   //!
 
     MTime   *fTimeEffOn;        //! Time to steer filling of fHTime
@@ -116 +65 @@
     Int_t fMap[5];              //!
 
-    void FitEnergySpec(Bool_t paint=kFALSE);
     Float_t FitEnergyBins(Bool_t paint=kFALSE);
     void FitThetaBins(Bool_t paint=kFALSE);
@@ -128 +76 @@
     Int_t DistancetoPrimitive(Int_t px, Int_t py) { return 0; }
 
+    virtual Bool_t      GetParameter(const MParList &pl);
+    virtual Double_t    GetVal() const;
+    virtual const char *GetParameterRule() const
+    {
+        return "MHillasSrc.fAlpha";
+    }
+
 public:
     MHAlpha(const char *name=NULL, const char *title=NULL);
 
+    // MH
     Bool_t SetupFill(const MParList *pl);
     Bool_t Fill(const MParContainer *par, const Stat_t w=1);
     Bool_t Finalize();
 
-    const TH3D &GetHist() const { return fHAlpha; }
+    // MHAlpha
+    const TH3D &GetHist() const { return fHist; }
     const MAlphaFitter &GetAlphaFitter() const { return fFit; }
 
     const TH1D &GetHEnergy() const { return fHEnergy; }
 
+    // Setter
+    void SetNameParameter(const char *name) { fNameParameter=name; }
     void SetOffData(const MHAlpha &h)
     {
-        fOffData        = &h.fHAlpha;
+        fOffData        = &h.fHist;
         fForceUsingSize =  h.fForceUsingSize;
         fNumTimeBins    =  h.fNumTimeBins;