Changeset 18100 for branches/Mars_McMismatchStudy/fact/analysis/mc

Timestamp:

02/03/15 08:27:19 (10 years ago)

Author:

ghughes

Message:

Added code to reconstruct the energy using MC lookup tables

File:

• branches/Mars_McMismatchStudy/fact/analysis/mc/ganymed.C (modified) (3 diffs)

branches/Mars_McMismatchStudy/fact/analysis/mc/ganymed.C

@@ -23 +23 @@
     write->AddContainer("MTime",          "Events", kFALSE);
     write->AddContainer("MMcEvt",         "Events", kFALSE);
-    write->AddContainer("MEnergyEst",     "Events", kFALSE);
+    //write->AddContainer("MEnergyEst",     "Events", kFALSE);
+    //write->AddContainer("MEnergyTable",   "Events", kFALSE);
     write->AddContainer("DataType",       "Events");
     write->AddContainer("Weight",         "Events");
@@ -170 +171 @@
     energyest.SetRule("MHillas.fSize^0.815*pow(10, 0.0199*(MHillasSrc.fDist*1.17e-2 - (-6.95))^2*(MHillasSrc.fDist*1.17e-2>(-6.95)) + 0.5997 *(MHillasSrc.fDist*1.17e-2 - 0.7568)^2 * (MHillasSrc.fDist*1.17e-2>0.7568) )");
     //    energyest.SetRule("1737 * (MHillas.fSize / 300.)^0.810");
+
+    //MEnergyTable energytable;
+    // Tables available from GH
+    //energytable.SetTableFile("Table_standard_test.root");
 
     // Cuts before writing
@@ -283 +288 @@
     tlist2.AddToList(&cuts);
     tlist2.AddToList(&energyest);
+    //tlist2.AddToList(&energytable);
     tlist2.AddToList(write0);
     tlist2.AddToList(&cont0);  // Post/Pre shown cuts