Index: trunk/MagicSoft/Mars/mhflux/MHCollectionArea.h
===================================================================
--- trunk/MagicSoft/Mars/mhflux/MHCollectionArea.h	(revision 6949)
+++ trunk/MagicSoft/Mars/mhflux/MHCollectionArea.h	(revision 6954)
@@ -36,4 +36,5 @@
 
     void Calc(TH2D &hsel, TH2D &hall);
+    void CalcEfficiency();
 
 public:
@@ -43,9 +44,8 @@
     Bool_t ReInit(MParList *pList);
     Bool_t Fill(const MParContainer *par, const Stat_t weight=1);
+    Bool_t Finalize() { CalcEfficiency(); return kTRUE; }
 
     void Draw(Option_t *option="");
     void Paint(Option_t *option="");
-
-    void CalcEfficiency();
 
     void SetHistAll(const TH2D &h) { h.Copy(fHistAll); fIsExtern=kTRUE; }
@@ -53,4 +53,5 @@
     TH2D &GetHistAll()             { return fHistAll; }
     const TH2D &GetHistAll() const { return fHistAll; }
+    const TH1D &GetHEnergy() const { return fHEnergy; }
 
 /*
Index: trunk/MagicSoft/Mars/mhflux/MHEnergyEst.cc
===================================================================
--- trunk/MagicSoft/Mars/mhflux/MHEnergyEst.cc	(revision 6949)
+++ trunk/MagicSoft/Mars/mhflux/MHEnergyEst.cc	(revision 6954)
@@ -219,4 +219,22 @@
 }
 
+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)
 {
Index: trunk/MagicSoft/Mars/mhflux/MHEnergyEst.h
===================================================================
--- trunk/MagicSoft/Mars/mhflux/MHEnergyEst.h	(revision 6949)
+++ trunk/MagicSoft/Mars/mhflux/MHEnergyEst.h	(revision 6954)
@@ -49,4 +49,5 @@
 
     const TH3D &GetHEnergy() const { return fHEnergy; }
+    void GetWeights(TH1D &hist) const;
 
     void InitMapping(MHMatrix *mat);