Changeset 2124 for trunk/MagicSoft/Mars/macros/estfit.C

Timestamp:

05/20/03 11:08:23 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/estfit.C (modified) (9 diffs)

trunk/MagicSoft/Mars/macros/estfit.C

@@ -33 +33 @@
 
 // --------------------------------------------------------------------------
-
-void estfit(Bool_t evalenergy=kFALSE)
+//
+//  0: fit impact parameter only
+//  1: fit energy only
+//  2: fit all parameters with respect to the energy resolution
+//
+void estfit(Int_t evalenergy=0)
 {
     //
@@ -48 +52 @@
     eest.InitMapping(&matrix);
 
-    MReadTree read("Events", "~/ct1/MC_ON2.root");
+    MReadTree read("Events", "MC_ON2_short.root");
     read.DisableAutoScheme();
 
@@ -92 +96 @@
     //
     TArrayD fA(5);
-    fA[0] = -3907.74; //4916.4;    //-2414.75;
-    fA[1] = 1162.3;   //149.549;   // 1134.28;
-    fA[2] = 199.351;  //-558.209;  // 132.932;
-    fA[3] = 0.403192; //0.270725;  //0.292845;
-    fA[4] = 121.921;  //107.001;   // 107.001;
-
     TArrayD fB(7);
-    fB[0] = -100;
-    fB[1] = 10;
-    fB[2] = 0.1;
-    fB[3] = 10;
-    fB[4] = -1;
-    fB[5] = -0.1;
-    fB[6] = -0.1;
-
+    fA[0] =  392957;
+    fA[1] =  135;
+    fA[2] =  -37449;
+    fA[3] =  0.3464;
+    fA[4] =  1;
+    fB[0] =  1;
+    fB[1] =  1;
+    fB[2] =  1;
+    fB[3] =  1;
+    fB[4] =  1;
+    fB[5] =  1;
+    fB[6] =  1;
+     /*
+     fA[0] = 6122.97;
+     fA[1] = 144.889;
+     fA[2] = -601.256;
+     fA[3] = 0.00171985;
+     fA[4] = 116.451;
+     fB[0] =  -2368.21;
+     fB[1] =  1186.26;
+     fB[2] =  0.147235;
+     fB[3] =  144.49;
+     fB[4] =  42.7681;
+     fB[5] = -0.757817;
+     fB[6] =  0.182482;
+     */
+    /*
+     TArrayD fA(5);
+    fA[0] = 6122.97;
+    fA[1] = 144.889;
+    fA[2] = -601.256;
+    fA[3] = 0.00171985;
+    fA[4] = 116.451;
+
+    TArrayD fB(7);
+    fB[0] =  -10445.5;
+    fB[1] =  2172.05;
+    fB[2] =  0.69;
+    fB[3] =  491.2;
+    fB[4] =  4.71444;
+    fB[5] = -0.0331926;
+    fB[6] = -0.014833;
+    */
     // Set starting values and step sizes for parameters
     for (Int_t i=0; i<fA.GetSize(); i++)
@@ -120 +153 @@
 
         Bool_t rc = minuit.DefineParameter(i, name, vinit, step, limlo, limup);
-        if (evalenergy)
+        if (evalenergy==1)
+        {
             minuit.FixParameter(i);
+            cout << "Fixed Parameter #" << i << endl;
+        }
 
         if (!rc)
@@ -142 +178 @@
 
         Bool_t rc = minuit.DefineParameter(i+fA.GetSize(), name, vinit, step, limlo, limup);
-        if (!evalenergy)
+        if (evalenergy==0)
+        {
             minuit.FixParameter(i+fA.GetSize());
+            cout << "Fixed Parameter #" << i+fA.GetSize() << endl;
+        }
 
         if (!rc)
@@ -156 +195 @@
     //
     minuit.SetObjectFit(&evtloop);
+
+    cout << endl << "Fitting procedure running..." << endl;
 
     TStopwatch clock;
@@ -164 +205 @@
     Bool_t rc = minuit.Migrad();
     gLog.SetNullOutput(kFALSE);
-
     if (rc)
     {
@@ -176 +216 @@
     cout << endl;
 
-    for (Int_t i=(evalenergy?fA.GetSize():0); i<(evalenergy?fA.GetSize()+fB.GetSize():fA.GetSize()); i++)
+    for (Int_t i=(evalenergy==1?fA.GetSize():0); i<(evalenergy>0?fA.GetSize()+fB.GetSize():fA.GetSize()); i++)
+    //for (Int_t i=0; i<fA.GetSize()+fB.GetSize(); i++)
     {
         Double_t val;
@@ -182 +223 @@
 
         if (!minuit.GetParameter(i, val, er))
-        {
             cout << "Error getting parameter #" << i << endl;
-            return;
-        }
 
         cout << i << ":  " << val << "  +-  " << er << endl;