Changeset 7152 for trunk/MagicSoft/Mars/mtools

Timestamp:

06/13/05 10:33:55 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtools

Files:

• MChisqEval.cc (modified) (8 diffs)
• MChisqEval.h (modified) (1 diff)

trunk/MagicSoft/Mars/mtools/MChisqEval.cc

@@ -25 +25 @@
 /////////////////////////////////////////////////////////////////////////////
 //
+// MChisqEval
+//
+// Evaluates a chisq from one or two MParameterD calculated in an eventloop.
 //
 /////////////////////////////////////////////////////////////////////////////
@@ -30 +33 @@
 
 #include <fstream>
+
+#include "MLog.h"
+#include "MLogManip.h"
 
 #include "MDataChain.h"
@@ -40 +46 @@
 using namespace std;
 
-void MChisqEval::StreamPrimitive(ofstream &out) const
-{
-    out << "   MChisqEval " << GetUniqueName() << ";";
-    if (fData0)
-        out << "   " << GetUniqueName() << ".SetY1(\"" << fData0->GetRule() << "\");" << endl;
-    if (fData1)
-        out << "   " << GetUniqueName() << ".SetY1(\"" << fData1->GetRule() << "\");" << endl;
-}
+const TString MChisqEval::gsDefName  = "MChisqEval";
+const TString MChisqEval::gsDefTitle = "Evaluate a chisq";
 
-MChisqEval::MChisqEval(const char *name, const char *title) : fData0(NULL), fData1(NULL)// : fMatrix(mat), fColumn(col), fEvalE(evale)
+MChisqEval::MChisqEval(const char *name, const char *title) : fData0(NULL), fData1(NULL), fNameResult("MinimizationValue")
 {
     fName  = name  ? name  : gsDefName.Data();
@@ -55 +55 @@
 }
 
-MChisqEval::MChisqEval(MData *y1, const char *name, const char *title) : fData0(NULL), fData1(NULL)// : fMatrix(mat), fColumn(col), fEvalE(evale)
+MChisqEval::MChisqEval(MData *y1, const char *name, const char *title) : fData0(NULL), fData1(NULL), fNameResult("MinimizationValue")
 {
     fName  = name  ? name  : gsDefName.Data();
@@ -62 +62 @@
 }
 
-MChisqEval::MChisqEval(MData *y1, MData *y2, const char *name, const char *title) : fData0(NULL), fData1(NULL)// : fMatrix(mat), fColumn(col), fEvalE(evale)
+MChisqEval::MChisqEval(MData *y1, MData *y2, const char *name, const char *title) : fData0(NULL), fData1(NULL), fNameResult("MinimizationValue")
 {
     fName  = name  ? name  : gsDefName.Data();
@@ -123 +123 @@
             return kFALSE;
 
-    fResult = (MParameterD*)plist->FindCreateObj("MParameterD", "MFitResult");
+    fResult = (MParameterD*)plist->FindCreateObj("MParameterD", fNameResult);
     if (!fResult)
         return kFALSE;
@@ -136 +136 @@
 
     const Double_t dy  = y2-y1;
-    const Double_t err = fData1 ? y1*y1 : 1;
+    const Double_t err = fData1 ? y1*y2 : 1;
 
     fChisq += dy*dy/err;
@@ -144 +144 @@
 Int_t MChisqEval::PostProcess()
 {
-    fChisq /= GetNumExecutions();
+    if (GetNumExecutions()>0)
+        fChisq /= GetNumExecutions();
 
     fResult->SetVal(fChisq);
+
+    *fLog << inf << GetDescriptor() << ": Result=" << fChisq << endl;
 
     return kTRUE;
 }
 
-const TString MChisqEval::gsDefName  = "MChisqEval";
-const TString MChisqEval::gsDefTitle = "Evaluate a chisq";
+void MChisqEval::StreamPrimitive(ofstream &out) const
+{
+    out << "   MChisqEval " << GetUniqueName() << ";";
+    if (fData0)
+        out << "   " << GetUniqueName() << ".SetY1(\"" << fData0->GetRule() << "\");" << endl;
+    if (fData1)
+        out << "   " << GetUniqueName() << ".SetY1(\"" << fData1->GetRule() << "\");" << endl;
+}
 

trunk/MagicSoft/Mars/mtools/MChisqEval.h

@@ -21 +21 @@
     MData   *fData1; // Data Member two (measured data)
 
-    // --------------------------------------------------------------------------
-    //
-    // Implementation of SavePrimitive. Used to write the call to a constructor
-    // to a macro. In the original root implementation it is used to write
-    // gui elements to a macro-file.
-    //
+    TString fNameResult;
+
     void StreamPrimitive(ofstream &out) const;