Changeset 8709 for trunk/MagicSoft/Mars/macros

Timestamp:

08/25/07 16:30:26 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/train/traindisp.C (modified) (4 diffs)

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

@@ -88 +88 @@
 
     /*
-     -------------------- Magic-Cuts ----------------------
+     // -------------------- Magic-Cuts ----------------------
      MFMagicCuts cuts;
      cuts.SetHadronnessCut(MFMagicCuts::kArea);
@@ -111 +111 @@
      opt.AddPreCut(&cuts);
 
-     -------------------- Energy Slope --------------------
-     MFEnergySlope slope(-2.8);
-     opt.AddPreCut(&slope);
-
-     -------------------- Other cuts ----------------------
+     // -------------------- Other cuts ----------------------
      opt.AddPreCut("MNewImagePar.fLeakage1<0.0001");
      opt.AddPreCut("MHillas.fSize>200");
@@ -121 +117 @@
      opt.AddPreCut("MPointingPos.fZd<25");
 
-     ------------------ Zd distribution -------------------
-     TFile file("ganymed00001111.root");
-
-     MStatusArray arr;
-     if (arr.Read()<=0)
-        return;
-     TH1D *vstime = (TH1D*)arr.FindObjectInCanvas("Theta",  "TH1D", "OnTime");
-     if (!vstime)
-         return;
-
-     MMcSpectrumWeight weight;
-     weight.SetWeightsZd(vstime);
-     opt.AddPreTask(&weight);
-
-     ---------------------- Histogram  --------------------
-
+     // ---------------------- Histogram  --------------------
      MHn hist("MyHist", "Energy Residual (lg E_{est} - lg E_{mc})");
 
@@ -157 +138 @@
      opt.AddTestTask(&fill);
 
-     ------------------------------------------------------
+     // -------------------- 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);
+
+     // ------------------------------------------------------
     */