Changeset 3957 for trunk/MagicSoft/Mars/mhist

Timestamp:

05/04/04 15:27:25 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

• MHCamera.cc (modified) (4 diffs)
• MHCamera.h (modified) (3 diffs)

trunk/MagicSoft/Mars/mhist/MHCamera.cc

@@ -34 +34 @@
 // To change the scale to a logarithmic scale SetLogy() of the Pad.
 //
+// You can correct for the abberation. Assuming that the distance
+// between the mean position of the light distribution and the position
+// of a perfect reflection on a perfect mirror in the distance r on
+// the camera plane is dr it is  d = a*dr  while a is the abberation
+// constant (for the MAGIC mirror it is roughly 0.0713). You can
+// set this constant by calling SetAbberation(a) which results in a
+// 'corrected' display (all outer pixels are shifted towards the center
+// of the camera to correct for this abberation)
+//
 // Be carefull: Entries in this context means Entries/bin or Events
 //
@@ -102 +111 @@
 //  Default Constructor. To be used by the root system ONLY.
 //
-MHCamera::MHCamera() : TH1D(), fGeomCam(NULL), fColors(kItemsLegend)
+MHCamera::MHCamera() : TH1D(), fGeomCam(NULL), fColors(kItemsLegend), fAbberation(0)
 {
     Init();
@@ -114 +123 @@
 //
 MHCamera::MHCamera(const MGeomCam &geom, const char *name, const char *title)
-: fGeomCam(NULL), fColors(kItemsLegend)
+: fGeomCam(NULL), fColors(kItemsLegend), fAbberation(0)
 {
     //fGeomCam = (MGeomCam*)geom.Clone();
@@ -759 +768 @@
         const MGeomPix &pix = (*fGeomCam)[i];
 
-        const Float_t x = pix.GetX()*conv;
-        const Float_t y = pix.GetY()*conv;
-        const Float_t d = pix.GetD()*conv;
+        Float_t x = pix.GetX()*conv/(fAbberation+1);
+        Float_t y = pix.GetY()*conv/(fAbberation+1);
+        Float_t d = pix.GetD()*conv;
 
         if (!isbox)

trunk/MagicSoft/Mars/mhist/MHCamera.h

@@ -47 +47 @@
 //#endif
 
+    Float_t fAbberation;
+
     void Init();
 
@@ -188 +190 @@
     void     Draw(Option_t *option="");
     TObject *DrawClone(Option_t *option="") const;
-    void     DrawProjection (Int_t fit=0) const;        
-    void     DrawRadialProfile()           const;       
+    void     DrawProjection (Int_t fit=0) const;
+    void     DrawRadialProfile()           const;
 
     void     SavePrimitive(ofstream &out, Option_t *);
@@ -209 +211 @@
     Bool_t   IsFreezed() const { return TestBit(kFreezed); }
     //void  SetOptStat(Int_t os=-1) { fOptStat = os; } // *MENU*
+
+    void     SetAbberation(Float_t f=0.0713) { fAbberation=f; } // *MENU*
 
     void     AddNotify(TObject *event);