Changeset 8708 for trunk/MagicSoft/Mars/macros/train/trainenergy.C

Timestamp:

08/24/07 19:00:47 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/train/trainenergy.C (modified) (2 diffs)

trunk/MagicSoft/Mars/macros/train/trainenergy.C

@@ -118 +118 @@
 
     /*
-     -------------------- Energy Slope --------------------
-     // This is a way to throw away events to a different slope
-     MFEnergySlope slope(-4.0); // New slope for mc spectrum
-     opt.AddPreCut(&slope);     // throw away events to change slope
-
-     // This is a way to weight the events to a different spectrum
-     MMcSpectrumWeight weight;
-     weight.SetFormula("pow(X/300, -2.31-0.26*log10(X/300))");
-     opt.SetWeights(&weight);
-
-     -------------------- Other cuts ----------------------
-     opt.AddPreCut("MHillasSrc.fDist*MGeomCam.fConvMm2Deg<1.0");
-     opt.AddPreCut("MHillas.fSize>200");
-
-     ------------------ Zd distribution -------------------
-     TFile file("ganymed00001111.root");
-
-     MStatusArray arr;
-     if (arr.Read()<=0)
-        return;
-     TH1D *vstime = (TH1D*)arr.FindObjectInCanvas("Theta",  "TH1D", "OnTime");
-     if (!vstime)
-         return -1;
-
-     MMcSpectrumWeight weight;
-     weight.SetWeightsZd(vstime);
-     opt.AddPreTask(&weight);
-
-     ---------------------- Histogram  --------------------
+     // ---------------------- Histogram  --------------------
      MHn hist("MyHist", "Energy Residual (lg E_{est} - lg E_{mc})");
 
@@ -161 +133 @@
      opt.AddTestTask(&fill);
 
-     ------------------------------------------------------
+     // -------------------- Other cuts ----------------------
+     opt.AddPreCut("MHillasSrc.fDist*MGeomCam.fConvMm2Deg<1.0");
+     opt.AddPreCut("MHillas.fSize>200");
+
+     // -------------------- Energy Slope --------------------
+     // Note, that weight normally doesn't improve anything here.
+     // This is a way to throw away events to a different slope
+     MFEnergySlope slope(-4.0); // New slope for mc spectrum
+     opt.AddPreCut(&slope);     // throw away events to change slope
+
+     // This is a way to weight the events to a different spectrum
+     MMcSpectrumWeight weight;
+     weight.SetFormula("pow(X/300, -2.31-0.26*log10(X/300))");
+     opt.SetWeights(&weight);
+
+     // ------------------ Zd distribution -------------------
+     TFile file("ganymed00001111.root");
+
+     MStatusArray arr;
+     if (arr.Read()<=0)
+        return;
+     TH1D *vstime = (TH1D*)arr.FindObjectInCanvas("Theta",  "TH1D", "OnTime");
+     if (!vstime)
+         return -1;
+
+     MMcSpectrumWeight weight;
+     weight.SetWeightsZd(vstime);
+     opt.AddPreTask(&weight);
+
+     // ------------------------------------------------------
     */