Changeset 1263

Timestamp:

03/20/02 10:36:39 (23 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• macros/flux.C (modified) (1 diff)
• mhist/MHHillas.cc (modified) (2 diffs)
• mhist/MHHillas.h (modified) (1 diff)
• mhist/MHStarMap.cc (modified) (1 diff)
• mhist/MHTimeDiffTheta.cc (modified) (2 diffs)
• mhist/MHTimeDiffTime.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/macros/flux.C

@@ -116 +116 @@
     fillsp.SetFilter(&fasrc);
     fillasp.SetFilter(&fsrc);
-    fillsptheta.SetFilter(&fsrc);
-    fillasptheta.SetFilter(&fasrc);
+    fillsptheta.SetFilter(&fasrc);
+    fillasptheta.SetFilter(&fsrc);
 
     MFillH fillontime ("EffOnTime  [MHTimeDiffTime]",  "MMcEvt");

trunk/MagicSoft/Mars/mhist/MHHillas.cc

@@ -53 +53 @@
 //
 MHHillas::MHHillas(const char *name, const char *title)
-    : fMm2Deg(1), fUseMmScale(kFALSE)
+    : fMm2Deg(-1), fUseMmScale(kFALSE)
 {
     //
@@ -152 +152 @@
 // --------------------------------------------------------------------------
 //
+// Use this function to setup your own conversion factor between degrees
+// and millimeters. The conversion factor should be the one calculated in
+// MGeomCam. Use this function with Caution: You could create wrong values
+// by setting up your own scale factor.
+//
+void MHHillas::SetMm2Deg(Float_t mmdeg)
+{
+    if (mmdeg<0)
+    {
+        *fLog << warn << dbginf << "Warning - Conversion factor < 0 - nonsense. Ignored." << endl;
+        return;
+    }
+
+    if (fMm2Deg>=0)
+        *fLog << warn << dbginf << "Warning - Conversion factor already set. Overwriting" << endl;
+
+    fMm2Deg = mmdeg;
+}
+
+// --------------------------------------------------------------------------
+//
+// With this function you can convert the histogram ('on the fly') between
+// degrees and millimeters.
+//
+void MHHillas::SetMmScale(Bool_t mmscale=kTRUE)
+{
+    if (fUseMmScale == mmscale)
+        return;
+
+    if (fMm2Deg<0)
+    {
+        *fLog << warn << dbginf << "Warning - Sorry, no conversion factor for conversion available." << endl;
+        return;
+    }
+
+    if (fUseMmScale)
+    {
+        fLength->GetXaxis()->SetTitle("Length [mm]");
+        fWidth->GetXaxis()->SetTitle("Width [mm]");
+
+        fLength->Scale(1./fMm2Deg);
+        fWidth->Scale(1./fMm2Deg);
+    }
+    else
+    {
+        fLength->GetXaxis()->SetTitle("Length [\\circ]");
+        fWidth->GetXaxis()->SetTitle("Width [\\circ]");
+
+        fLength->Scale(fMm2Deg);
+        fWidth->Scale(fMm2Deg);
+    }
+
+    fUseMmScale = mmscale;
+}
+
+// --------------------------------------------------------------------------
+//
 // Draw clones of all four histograms. So that the object can be deleted
 // and the histograms are still visible in the canvas.

trunk/MagicSoft/Mars/mhist/MHHillas.h

@@ -23 +23 @@
     ~MHHillas();
 
-    void SetMmScale(Bool_t mmscale=kTRUE) { fUseMmScale = mmscale; }
+    void SetMmScale(Bool_t mmscale=kTRUE);
+    void SetMm2Deg(Float_t mmdeg);
 
     Bool_t SetupFill(const MParList *pList);

trunk/MagicSoft/Mars/mhist/MHStarMap.cc

@@ -98 +98 @@
 
     const float m = tan(delta);
-    const float t = m*h.GetMeanX()-h.GetMeanY();
+    const float t = h.GetMeanY() - m*h.GetMeanX();
 
     if (m>-1 && m<1)

trunk/MagicSoft/Mars/mhist/MHTimeDiffTheta.cc

@@ -106 +106 @@
     // FIXME: ProjectionX,Y is not const within root
     //
-    TH1 *h;
+    TH1D *h;
 
     c->cd(1);
@@ -133 +133 @@
         MakeDefCanvas("DiffTimeTheta", "Distrib of Delta t, Theta");
 
-    TH1 *h;
+    TH1D *h;
 
     gPad->Divide(2,2);

trunk/MagicSoft/Mars/mhist/MHTimeDiffTime.cc

@@ -98 +98 @@
     // FIXME: ProjectionX,Y is not const within root
     //
-    TH1 *h;
+    TH1D *h;
 
     c->cd(1);
@@ -127 +127 @@
     gPad->Divide(2,2);
 
-    TH1 *h;
+    TH1D *h;
 
     gPad->cd(1);