Changeset 19160 for trunk/Mars

Timestamp:

08/14/18 10:20:02 (6 years ago)

Author:

Daniela Dorner

Message:

new disp parametrization for analysis + cut that had not been committed yet

File:

• trunk/Mars/fact/analysis/ganymed.C (modified) (6 diffs)

trunk/Mars/fact/analysis/ganymed.C

@@ -43 +43 @@
     // Crab (05h34'32" 22d00'52")
     MPointingPos source("MSourcePos");
-    //source.SetSkyPosition(5.5755, 22.0144); // ra[h], dec[deg]
+    //source.SetSkyPosition(5.5755, 22.0144); // ra[h], dec[deg] // Crab
     //source.SetSkyPosition(12.356089, 30.176897); // ra[h], dec[deg]//1218
     //source.SetSkyPosition(16.897867,39.760201);// mrk501
@@ -54 +54 @@
 
     // FIXME: If fPathIn read cuts and energy estimator from file!
-    MContinue contq("MImagePar.fNumIslands>3.5"
-                    "|| MNewImagePar.fNumUsedPixels<5.5"
+    MContinue contq("MNewImagePar.fNumUsedPixels<5.5"
+                 // "|| MImagePar.fNumIslands>3.5"
                     "|| MNewImagePar.fLeakage1>0.1",
                     "CutQ");
 
     MContinue cont0("MHillas.fSize<0"
-                    "|| MHillasExt.fSlopeSpreadWeighted>0.68"
-                    "|| MHillasExt.fSlopeSpreadWeighted<0.18"
-                    "|| log10(MHillas.GetArea*MGeomCam.fConvMm2Deg^2)>(log10(MHillas.fSize)-2)*1.1-1.55"
-                    "|| MNewImagePar.fConcCore < 0.13"
-                    "|| MNewImagePar.fConcCOG  < 0.15"
+                    "|| MHillas.GetArea > log10(MHillas.fSize)*898-1535 "
                     , "Cut0");
 
@@ -82 +78 @@
     // Parametrization of Disp
     // for new cleaning
-    param[0]  =  1.39252;  // constant
-    param[8]  =  0.154247; // slope
+    //param[0]  =  1.39252;  // constant
+    param[0]  =  1.34;  // constant
+    //param[8]  =  0.154247; // slope
+    param[8]  =  0.0755; // slope
     param[9]  =  1.67972;  // leak
     param[10] =  4.86232;  // leak
@@ -96 +94 @@
     param[7]  =  0.5;
     // ThetaSq-Cut
-    param[1]  =  0.11;  // 0.215
+    param[1]  =  0.156;  // 0.215
     // Area-Cut
     param[2]  =  0.215468;
@@ -114 +112 @@
     // --------------------------------------------------------------------------------
 
+    //MBinning bins3( 67, -0.005, 0.665, "BinningTheta",  "asin");
     MBinning bins3(100, -0.005, 0.995, "BinningTheta",  "asin");
     MBinning binsD( 25,   0,      2.5, "BinningDist");
@@ -466 +465 @@
         path += fname1;
 
-        d->SaveAs(path);
+        //d->SaveAs(path);
+        d->Write();
     }