Index: trunk/MagicSoft/Mars/Changelog
===================================================================
--- trunk/MagicSoft/Mars/Changelog	(revision 7382)
+++ trunk/MagicSoft/Mars/Changelog	(revision 7384)
@@ -19,4 +19,23 @@
                                                  -*-*- END OF LINE -*-*-
 
+ 2005/10/10 Thomas Bretz
+
+   * mbase/MMath.[h,cc]:
+     - added new function SignificanceLiMaExc
+
+   * mbase/MTime.[h,cc]:
+     - added new Setter SetAxisTime
+
+   * mcalib/MCalibCalcFromPast.cc:
+     - fixed some format in output
+
+   * mhcalib/MHCalibrationCam.cc:
+     - fixed a typo
+
+   * mhflux/MAlphaFitter.[h,cc]:
+     - added new data member fSiginificanceExc
+
+
+
  2005/10/01 Thomas Bretz
 
@@ -30,5 +49,5 @@
 
    * ganymed.cc:
-     - implemented twonew options --ind and --ins
+     - implemented two new options --ind and --ins
 
    * datacenter/scripts/runganymed:
Index: trunk/MagicSoft/Mars/mbase/MMath.cc
===================================================================
--- trunk/MagicSoft/Mars/mbase/MMath.cc	(revision 7382)
+++ trunk/MagicSoft/Mars/mbase/MMath.cc	(revision 7384)
@@ -133,4 +133,17 @@
 // --------------------------------------------------------------------------
 //
+// Return Li/Ma (5) for the error of the excess, under the assumption that
+// the existance of a signal is already known.
+//
+Double_t MMath::SignificanceLiMaExc(Double_t s, Double_t b, Double_t alpha)
+{
+    Double_t Ns = s - alpha*b;
+    Double_t sN = s - alpha*alpha*b;
+
+    return Ns<0 || sN<0 ? 0 : Ns/TMath::Sqrt(sN);
+}
+
+// --------------------------------------------------------------------------
+//
 // Returns: 2/(sigma*sqrt(2))*integral[0,x](exp(-(x-mu)^2/(2*sigma^2)))
 //
Index: trunk/MagicSoft/Mars/mbase/MMath.h
===================================================================
--- trunk/MagicSoft/Mars/mbase/MMath.h	(revision 7382)
+++ trunk/MagicSoft/Mars/mbase/MMath.h	(revision 7384)
@@ -22,4 +22,5 @@
     Double_t SignificanceLiMa(Double_t s, Double_t b, Double_t alpha=1);
     Double_t SignificanceLiMaSigned(Double_t s, Double_t b, Double_t alpha=1);
+    Double_t SignificanceLiMaExc(Double_t s, Double_t b, Double_t alpha=1);
 
     void ReducePrecision(Float_t &val);
Index: trunk/MagicSoft/Mars/mbase/MTime.cc
===================================================================
--- trunk/MagicSoft/Mars/mbase/MTime.cc	(revision 7382)
+++ trunk/MagicSoft/Mars/mbase/MTime.cc	(revision 7384)
@@ -230,4 +230,13 @@
 {
     return (GetMjd()-49718)*kDay/1000;
+}
+
+// --------------------------------------------------------------------------
+//
+// Counterpart of GetAxisTime
+//
+void MTime::SetAxisTime(Double_t time)
+{
+    SetMjd(time*1000/kDay+49718);
 }
 
Index: trunk/MagicSoft/Mars/mbase/MTime.h
===================================================================
--- trunk/MagicSoft/Mars/mbase/MTime.h	(revision 7382)
+++ trunk/MagicSoft/Mars/mbase/MTime.h	(revision 7384)
@@ -87,4 +87,5 @@
     Bool_t   SetMjd(UInt_t mjd, ULong_t ms, UInt_t ns=0);
     void     SetMjd(Double_t m);
+    void     SetAxisTime(Double_t time);
 
     // Getter functions
Index: trunk/MagicSoft/Mars/mcalib/MCalibCalcFromPast.cc
===================================================================
--- trunk/MagicSoft/Mars/mcalib/MCalibCalcFromPast.cc	(revision 7382)
+++ trunk/MagicSoft/Mars/mcalib/MCalibCalcFromPast.cc	(revision 7384)
@@ -438,5 +438,5 @@
     *fLog << inf << endl;
     *fLog << GetDescriptor() << " execution statistics:" << endl;
-    *fLog << " " << setw(7) << fNumFails << " (" << Form("%5.1f", 100.*fNumFails/fNumCam) << "%) updates failed." << endl;
+    *fLog << " " << setfill(' ') << setw(7) << fNumFails << " (" << Form("%5.1f", 100.*fNumFails/fNumCam) << "%) updates failed." << endl;
     *fLog << endl;
 
Index: trunk/MagicSoft/Mars/mhcalib/MHCalibrationCam.cc
===================================================================
--- trunk/MagicSoft/Mars/mhcalib/MHCalibrationCam.cc	(revision 7382)
+++ trunk/MagicSoft/Mars/mhcalib/MHCalibrationCam.cc	(revision 7384)
@@ -967,5 +967,5 @@
     return kTRUE;
 
-  *fLog << inf << GetDescriptor() << ": Number of event used to fill histograms == " << fCurrentNumEvts << endl;
+  *fLog << inf << GetDescriptor() << ": Number of events used to fill histograms == " << fCurrentNumEvts << endl;
 
   if (fHiGainArray->GetSize() == 0 && fLoGainArray->GetSize() == 0)
Index: trunk/MagicSoft/Mars/mhflux/MAlphaFitter.cc
===================================================================
--- trunk/MagicSoft/Mars/mhflux/MAlphaFitter.cc	(revision 7382)
+++ trunk/MagicSoft/Mars/mhflux/MAlphaFitter.cc	(revision 7384)
@@ -60,4 +60,5 @@
 {
     fSignificance=0;
+    fSignificanceExc=0;
     fEventsExcess=0;
     fEventsSignal=0;
@@ -216,4 +217,5 @@
     fEventsExcess     = fEventsSignal-fEventsBackground;
     fSignificance     = MMath::SignificanceLiMaSigned(fEventsSignal, fEventsBackground);
+    fSignificanceExc  = MMath::SignificanceLiMaExc(fEventsSignal, fEventsBackground);
 
     if (TMath::IsNaN(fSignificance))
@@ -250,4 +252,5 @@
     fScaleFactor      = alpha;
     fSignificance     = MMath::SignificanceLiMaSigned(fEventsSignal, fEventsBackground/alpha, alpha);
+    fSignificanceExc  = MMath::SignificanceLiMaExc(fEventsSignal, fEventsBackground/alpha, alpha);
 
     if (TMath::IsNaN(fSignificance))
@@ -304,4 +307,5 @@
     // Result
     f.fSignificance     = fSignificance;
+    f.fSignificanceExc  = fSignificanceExc;
     f.fEventsExcess     = fEventsExcess;
     f.fEventsSignal     = fEventsSignal;
Index: trunk/MagicSoft/Mars/mhflux/MAlphaFitter.h
===================================================================
--- trunk/MagicSoft/Mars/mhflux/MAlphaFitter.h	(revision 7382)
+++ trunk/MagicSoft/Mars/mhflux/MAlphaFitter.h	(revision 7384)
@@ -54,4 +54,5 @@
     // Result
     Double_t fSignificance;     // significance of signal
+    Double_t fSignificanceExc;  // significance of excess
     Double_t fEventsExcess;     // calculated number of excess events (signal-bg)
     Double_t fEventsSignal;     // calculated number of signal events
@@ -154,4 +155,5 @@
 
     Double_t GetSignificance() const       { return fSignificance; }
+    Double_t GetSignificanceExc() const    { return fSignificanceExc; }
     Double_t GetChiSqSignal() const        { return fChiSqSignal; }
     Double_t GetChiSqBg() const            { return fChiSqBg; }