Changeset 6954 for trunk/MagicSoft/Mars/mhflux/MHEnergyEst.cc

Timestamp:

04/18/05 20:11:09 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhflux/MHEnergyEst.cc (modified) (1 diff)

trunk/MagicSoft/Mars/mhflux/MHEnergyEst.cc

@@ -219 +219 @@
 }
 
+void MHEnergyEst::GetWeights(TH1D &hist) const
+{
+    // Project into EnergyEst_ey
+    TH1D *h1 = (TH1D*)fHEnergy.Project3D("rtlprmft_ex");
+    TH1D *h2 = (TH1D*)fHEnergy.Project3D("rtlprmft_ey");
+
+    h1->Copy(hist);
+    hist.Divide(h2);
+
+    delete h1;
+    delete h2;
+
+    hist.SetNameTitle("EnergyRatio", "Ratio between estimated and monte carlo energy");
+    hist.SetXTitle("E [GeV]");
+    hist.SetYTitle("N_{est}/N_{mc} [1]");
+    hist.SetDirectory(0);
+}
+
 void MHEnergyEst::Paint(Option_t *opt)
 {