Changeset 6954 for trunk/MagicSoft/Mars/mhflux

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhflux

Files:

• MHCollectionArea.h (modified) (3 diffs)
• MHEnergyEst.cc (modified) (1 diff)
• MHEnergyEst.h (modified) (1 diff)

trunk/MagicSoft/Mars/mhflux/MHCollectionArea.h

@@ -36 +36 @@
 
     void Calc(TH2D &hsel, TH2D &hall);
+    void CalcEfficiency();
 
 public:
@@ -43 +44 @@
     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 +53 @@
     TH2D &GetHistAll()             { return fHistAll; }
     const TH2D &GetHistAll() const { return fHistAll; }
+    const TH1D &GetHEnergy() const { return fHEnergy; }
 
 /*

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)
 {

trunk/MagicSoft/Mars/mhflux/MHEnergyEst.h

@@ -49 +49 @@
 
     const TH3D &GetHEnergy() const { return fHEnergy; }
+    void GetWeights(TH1D &hist) const;
 
     void InitMapping(MHMatrix *mat);