Changeset 1227 for trunk/MagicSoft/Mars/mhist/MHMcCollectionArea.cc

Timestamp:

02/28/02 16:36:54 (23 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHMcCollectionArea.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mhist/MHMcCollectionArea.cc

@@ -30 +30 @@
 
 #include "MH.h"
+#include "MHMcEfficiency.h"
 
 ClassImp(MHMcCollectionArea);
@@ -192 +193 @@
 }
 
+void MHMcCollectionArea::Calc(MHMcEfficiency *heff)
+{
+    //
+    //  now calculate the Collection area for different
+    //  energies
+    //
+    TH2D &h = (TH2D&)*heff->GetHist();
+
+    const Int_t nbinx = h.GetXaxis()->GetNbins();
+    const Int_t nbiny = h.GetYaxis()->GetNbins();
+
+    for (Int_t ix=1; ix<=nbinx; ix++)
+    {
+        Double_t errA = 0;
+        Double_t colA = 0;
+
+        for (Int_t iy=1; iy<=nbiny; iy++)
+        {
+            TAxis *yaxis = h.GetYaxis();
+
+            const Double_t r1  = yaxis->GetBinLowEdge(iy);
+            const Double_t r2  = yaxis->GetBinLowEdge(iy+1);
+
+            const Double_t A   = TMath::Pi() * (r2*r2 - r1*r1);
+
+            const Double_t eff = h.GetCellContent(ix, iy);
+            const Double_t err = h.GetCellError(ix, iy);
+
+            colA += eff*A;
+            errA += A*A * err*err;
+        }
+
+        fHistCol->SetBinContent(ix, colA);
+        fHistCol->SetBinError(ix, sqrt(errA));
+    }
+}
+
 // --------------------------------------------------------------------------
 //
@@ -201 +239 @@
     // Description!
 
+    /*
     //
     //  first of all calculate the efficency
@@ -258 +297 @@
         fHistCol->SetBinContent(ix, colA);
         fHistCol->SetBinError(ix, sqrt(errA));
-    }
+        }
+        */
+    MHMcEfficiency heff;
+    heff.Calc(fHistSel, fHistAll);
+
+    Calc(&heff);
 
     SetReadyToSave();