Changeset 9308 for trunk/MagicSoft/Mars/macros/optim

Timestamp:

02/09/09 11:30:54 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/optim/optimdisp.C (modified) (3 diffs)

trunk/MagicSoft/Mars/macros/optim/optimdisp.C

@@ -18 +18 @@
     char *r = "M[0]*([0] + [1]*M[1] + [2]*M[2] + (M[3]>[3])*[4]*(M[3]-[3])^2)";
 
-    opt.FixParameter(0, 1.266195);
-    opt.FixParameter(1, 0.100577);
-    opt.FixParameter(2, 1.80309);
-    opt.FixParameter(3, 2.87177);
-    opt.FixParameter(4, 0.616823); 
+    opt.FixParameter(0, 1.15136);
+    opt.FixParameter(1, 0.0681437);
+    opt.FixParameter(2, 2.62932);
+    opt.FixParameter(3, 1.51279);
+    opt.FixParameter(4, 0.0507821);
 
     opt.AddPreCut("MNewImagePar.fLeakage1>0");
-    opt.AddPreCut("log10(MHillas.fSize)<3.2");
+    //opt.AddPreCut("log10(MHillas.fSize)<2.5"); //3.2
+
+    opt.AddPreCut("MHillasExt.fM3Long*sign(MHillasSrc.fCosDeltaAlpha)*MGeomCam.fConvMm2Deg>-0.07");
+    opt.AddPreCut("DataType.fVal>0.5");
 
     //opt.AddPreCut("abs(MHillasSrc.fDCA*MGeomCam.fConvMm2Deg)<0.2");
     //opt.AddPreCut("(MHillasSrc.fDist*MGeomCam.fConvMm2Deg-0.5)*7.2>MHillasExt.fSlopeLong*sign(MHillasSrc.fCosDeltaAlpha)/MGeomCam.fConvMm2Deg");
-    opt.AddPreCut("MHillasExt.fM3Long*sign(MHillasSrc.fCosDeltaAlpha)*MGeomCam.fConvMm2Deg>-0.07");
-    opt.AddPreCut("DataType.fVal>0.5");
 
     // -------------------- Run ----------------------------
@@ -66 +67 @@
      */
 
+    // opt.SetPathOut("optimdisp.root");
     opt.RunDisp("ganymed00000001-summary.root", r);
 }
@@ -73 +75 @@
    Par  |   0    1     2      3     4     |  Cut
  -------+---------------------------------+-----------------------
-  Fit 1 |  1.3  0.8  fix=0  fix=0  fix=0  | Leak1==0 lgSize<2.5
-  Fit 2 |  fix  fix  fix=0   2.4    0.3   | Leak1==0
+  Fit 1 |  1.3  0.1  fix=0  fix=0  fix=0  | Leak1==0 lgSize<2.5
+  Fit 2 |  fix  fix  fix=0   2.5    0.2   | Leak1==0
   Fit 2 |  fix  fix   1.8    fix    fix   | Leak1>0
   Fit 3 |  free free  fix    fix    fix   | -/-
+  Fit 4 |  fix  fix   fix    free   free  | -/-
+  Fit 5 |  fix  fix   free   fix    fix   | -/-
  -------+---------------------------------+-----------------------
 */