Changeset 9181 for trunk/MagicSoft/Mars/mjobs

Timestamp:

11/28/08 12:07:19 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mjobs/MJSpectrum.cc (modified) (7 diffs)

trunk/MagicSoft/Mars/mjobs/MJSpectrum.cc

@@ -879 +879 @@
 TString MJSpectrum::FormFlux(const TF1 &f, const char *unit)
 {
-    Double_t p0 = f.GetParameter(0);
-    Double_t p1 = f.GetParameter(1);
-
-    Double_t e0 = f.GetParError(0);
-    Double_t e1 = f.GetParError(1);
+    Double_t p0 = -f.GetParameter(0);
+    Double_t p1 =  f.GetParameter(1);
+
+    Double_t e0 =  f.GetParError(0);
+    Double_t e1 =  f.GetParError(1);
 
     MMath::Format(p0, e0);
@@ -894 +894 @@
                        FormFloat(p1/exp).Data(), FormFloat(e1/exp).Data(), i);
 
-    str += Form("#left(#frac{E}{%s}#right)^{%s #pm %s}", unit,
+    str += Form("#left(#frac{E}{%s}#right)^{-%s #pm %s}", unit,
                 FormFloat(p0).Data(), FormFloat(e0).Data());
 
@@ -1013 +1013 @@
     // collection area and cancel out.
     //spectrum.Multiply(&weights);
-    spectrum.SetNameTitle("Flux", "Spectrum (after unfolding)");
+    spectrum.SetNameTitle("Flux", "Spectrum");
     spectrum.SetBit(TH1::kNoStats);
 
@@ -1728 +1728 @@
     // ---------------------------------------------------------
 
+    MFillH fillsp("MHSrcPosCam", "MSrcPosCam", "FillSrcPosCam");
     MFillH fill0a(&henergy, "", "FillEventDist");
     MFillH fill1a("MHHillasMCPre  [MHHillas]",      "MHillas",      "FillHillasPre");
@@ -1737 +1738 @@
     MFillH fill7a("MHNewParMCPost [MHNewImagePar]", "MNewImagePar", "FillNewParPost");
     MFillH fill8a("ExcessMC       [MH3]",           "",             "FillExcessMC");
+    fillsp.SetNameTab("SrcPos");
     fill0a.SetNameTab("EvtDist");
     fill1a.SetNameTab("PreCut");
@@ -1754 +1756 @@
     fill7a.SetWeight();
     fill8a.SetWeight();
+    fillsp.SetWeight();
 
     // FIXME: To be done: A task checking the lower 1% after the lower
@@ -1785 +1788 @@
     //tlist2.AddToList(&hcalc2);
     tlist2.AddToList(&weight);
+    tlist2.AddToList(&fillsp);
     tlist2.AddToList(&fill1a);
     tlist2.AddToList(&fill30);