Changeset 18100 for branches/Mars_McMismatchStudy

Timestamp:

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

Author:

ghughes

Message:

Added code to reconstruct the energy using MC lookup tables

Location:

branches/Mars_McMismatchStudy

Files:

• fact/analysis/mc/ganymed.C (modified) (3 diffs)
• manalysis/AnalysisLinkDef.h (modified) (1 diff)
• manalysis/MEnergyTable.cc (added)
• manalysis/MEnergyTable.h (added)

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

branches/Mars_McMismatchStudy/manalysis/AnalysisLinkDef.h

@@ -25 +25 @@
 #pragma link C++ class MMcCalibrationUpdate+;
 
+#pragma link C++ class MEnergyTable+;
+
 #endif