Changeset 1463 for trunk/MagicSoft/Mars/mhist/MHHillasSrc.cc

Timestamp:

07/31/02 18:03:51 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHHillasSrc.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/mhist/MHHillasSrc.cc

@@ -68 +68 @@
     // connect all the histogram with the container fHist
     //
-    fAlpha    = new TH1F("Alpha",    "Alpha of Ellipse",             90,  0,  90);
-    fDist     = new TH1F("Dist",     "Dist of Ellipse",             100,  0, 600);
-    fHeadTail = new TH1F("HeadTail", "HeadTail of Ellipse",          45,  0,  90);
-    fCosDA    = new TH1F("CosDA",    "cos(Delta,Alpha) of Ellipse", 100, -1,   1);
+    fAlpha    = new TH1F("Alpha",    "Alpha of Ellipse",             90,    0,  90);
+    fDist     = new TH1F("Dist",     "Dist of Ellipse",             100,    0, 445);
+    fHeadTail = new TH1F("HeadTail", "HeadTail of Ellipse",         101, -445, 445);
+    fCosDA    = new TH1F("CosDA",    "cos(Delta,Alpha) of Ellipse", 101,   -1,   1);
 
     fAlpha->SetDirectory(NULL);
@@ -80 +80 @@
     fAlpha->SetXTitle("\\alpha [\\circ]");
     fDist->SetXTitle("Dist [mm]");
-    fHeadTail->SetXTitle("Head-Tail [\\circ]");
+    fHeadTail->SetXTitle("Head-Tail [mm]");
     fCosDA->SetXTitle("cos(\\delta,\\alpha) [mm]");
 
@@ -113 +113 @@
 Bool_t MHHillasSrc::SetupFill(const MParList *plist)
 {
-    const MBinning* binsa = (MBinning*)plist->FindObject("BinningAlpha");
-    const MBinning* binsd = (MBinning*)plist->FindObject("BinningDist");
-    if (!binsa || !binsd)
-    {
-        *fLog << err << dbginf << "At least one MBinning not found... aborting." << endl;
-        return kFALSE;
-    }
-
-    SetBinning(fAlpha, binsa);
-    SetBinning(fDist,  binsd);
-
     const MGeomCam *geom = (MGeomCam*)plist->FindObject("MGeomCam");
     if (!geom)
-    {
         *fLog << warn << dbginf << "No Camera Geometry available. Using mm-scale for histograms." << endl;
-        return kTRUE;
-    }
-
-    fDist->SetXTitle("Dist [\\circ]");
-
-    fMm2Deg = geom->GetConvMm2Deg();
-    fUseMmScale = kFALSE;
+    else
+    {
+        fMm2Deg = geom->GetConvMm2Deg();
+        SetMmScale(kFALSE);
+    }
+
+    ApplyBinning(*plist, "Alpha",    fAlpha);
+    ApplyBinning(*plist, "Dist",     fDist);
+    ApplyBinning(*plist, "HeadTail", fHeadTail);
 
     return kTRUE;
@@ -150 +140 @@
     fAlpha   ->Fill(fabs(h.GetAlpha()));
     fDist    ->Fill(fUseMmScale ? h.GetDist() : fMm2Deg*h.GetDist());
-    fHeadTail->Fill(h.GetHeadTail());
+    fHeadTail->Fill(fUseMmScale ? h.GetHeadTail() : fMm2Deg*h.GetHeadTail());
     fCosDA   ->Fill(h.GetCosDeltaAlpha());
 
     return kTRUE;
+}
+
+// --------------------------------------------------------------------------
+//
+// 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 MHHillasSrc::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 MHHillasSrc::SetMmScale(Bool_t mmscale)
+{
+    if (fUseMmScale == mmscale)
+        return;
+
+    if (fMm2Deg<0)
+    {
+        *fLog << warn << dbginf << "Warning - Sorry, no conversion factor for conversion available." << endl;
+        return;
+    }
+
+    const Double_t scale = mmscale ? 1./fMm2Deg : fMm2Deg;
+    MH::ScaleAxis(fDist,     scale);
+    MH::ScaleAxis(fHeadTail, scale);
+
+    if (mmscale)
+    {
+        fDist->SetXTitle("Dist [mm]");
+        fHeadTail->SetXTitle("Head-Tail [mm]");
+    }
+    else
+    {
+        fDist->SetXTitle("Dist [\\circ]");
+        fHeadTail->SetXTitle("Head-Tail [\\circ]");
+    }
+
+    fUseMmScale = mmscale;
 }
 
@@ -187 +232 @@
 
     c->cd(4);
+    gPad->SetLogy();
     fCosDA->DrawCopy();
 
@@ -221 +267 @@
 
     gPad->cd(2);
+    gPad->SetLogy();
     fCosDA->Draw();