Changeset 7193 for trunk/MagicSoft/Mars/mhflux/MHDisp.h

Timestamp:

07/14/05 18:45:35 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhflux/MHDisp.h (modified) (3 diffs)

trunk/MagicSoft/Mars/mhflux/MHDisp.h

@@ -5 +5 @@
 #include "MHFalseSource.h"
 #endif
+#ifndef ROOT_TH2
+#include <TH2.h>
+#endif
+#ifndef ROOT_TVector2
+#include <TVector2.h>
+#endif
 
 class MParList;
 class MHillasExt;
+class MHillasSrc;
 class MParameterD;
 
@@ -13 +20 @@
 {
 private:
-    MHillasExt  *fHilExt;  //!
-    MParameterD *fDisp;    //!
+    MHillasExt   *fHilExt;  //!
+    MHillasSrc   *fHilSrc;  //!
+    MParameterD  *fDisp;    //!
 
-    Double_t     fM3lCut;  //!
-    Double_t     fXi;      //!
-    Double_t     fXiTheta; //!
+    TH2D         fHistBg;
 
-    //Bool_t       fIsWobble;
+    TH2D         fHistBg1;
+    TH2D         fHistBg2;
+
+    TVector2     fFormerSrc;
+    Bool_t       fHalf;
+
+    Double_t     fSmearing;
+    Bool_t       fWobble;
 
     // MHDisp
     Double_t GetOffSignal(TH1 &h) const;
+    Double_t DeltaPhiSrc(const TVector2 &v) const;
+
+    void Profile(TH1 &h2, const TH2 &h1, Axis_t x0=0, Axis_t y0=0) const;
+    void MakeSignificance(TH2 &s, const TH2 &h1, const TH2 &h2, const Double_t scale=1) const;
+    void MakeDot(TH2 &h2) const;
+
+    Double_t Sq(Double_t x, Double_t y) const { return x*x+y*y; }
 
 public:
     MHDisp(const char *name=NULL, const char *title=NULL);
 
+    // MHDisp
+    void SetSmearing(Double_t s=-1) { fSmearing=s; }
+    void SetWobble(Bool_t w=kTRUE)  { fWobble=w;   }
+
+    // MHFalseSource (to be moved!)
+    void SetOffData(const MHFalseSource &fs)
+    {
+        MHFalseSource::SetOffData(fs);
+        if (dynamic_cast<const MHDisp*>(&fs))
+        {
+            const MHDisp &h = dynamic_cast<const MHDisp&>(fs);
+            fWobble   = h.fWobble;
+            fSmearing = h.fSmearing;
+
+            h.fHistBg1.Copy(fHistBg1);
+            h.fHistBg2.Copy(fHistBg2);
+
+            fHistBg1.SetDirectory(0);
+            fHistBg2.SetDirectory(0);
+        }
+    }
+
     // MH
     Bool_t SetupFill(const MParList *pList);
     Bool_t Fill(const MParContainer *par, const Stat_t w=1);
+    Bool_t Finalize();
 
     // TObject
@@ -37 +80 @@
 
     // MParContainer
-    //Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print=kFALSE);
+    Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print=kFALSE);
 
-    ClassDef(MHDisp, 1) //3D-histogram in alpha, x and y
+    ClassDef(MHDisp, 2) //Class to provide a Disp map
 };