Changeset 1300 for trunk/MagicSoft/Mars/mhist

Timestamp:

04/25/02 11:21:46 (23 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

• MH.cc (modified) (8 diffs)
• MH.h (modified) (2 diffs)
• MH3.cc (modified) (1 diff)
• MHMcCollectionArea.cc (modified) (2 diffs)
• MHMcCollectionArea.h (modified) (1 diff)

trunk/MagicSoft/Mars/mhist/MH.cc

@@ -50 +50 @@
 
 #include <TH1.h>
+#include <TH2.h>
+#include <TH3.h>
 #include <TCanvas.h>
 
@@ -136 +138 @@
 }
 
-void MH::SetBinning(TH1 *h, const MBinning *binsx, const MBinning *binsy)
+void MH::SetBinning(TH2 *h, const MBinning *binsx, const MBinning *binsy)
 {
     TAxis &x = *h->GetXaxis();
@@ -169 +171 @@
 }
 
-void MH::SetBinning(TH1 *h, const MBinning *binsx, const MBinning *binsy, const MBinning *binsz)
+void MH::SetBinning(TH3 *h, const MBinning *binsx, const MBinning *binsy, const MBinning *binsz)
 {
     //
@@ -214 +216 @@
 }
 
-void MH::SetBinning(TH1 *h, const TArrayD *binsx, const TArrayD *binsy)
+void MH::SetBinning(TH2 *h, const TArrayD *binsx, const TArrayD *binsy)
 {
     MBinning bx;
@@ -223 +225 @@
 }
 
-void MH::SetBinning(TH1 *h, const TArrayD *binsx, const TArrayD *binsy, const TArrayD *binsz)
+void MH::SetBinning(TH3 *h, const TArrayD *binsx, const TArrayD *binsy, const TArrayD *binsz)
 {
     MBinning bx;
@@ -245 +247 @@
 }
 
-void MH::SetBinning(TH1 *h, const TAxis *binsx, const TAxis *binsy)
+void MH::SetBinning(TH2 *h, const TAxis *binsx, const TAxis *binsy)
 {
     const Int_t nx = binsx->GetNbins();
@@ -260 +262 @@
 }
 
-void MH::SetBinning(TH1 *h, const TAxis *binsx, const TAxis *binsy, const TAxis *binsz)
+void MH::SetBinning(TH3 *h, const TAxis *binsx, const TAxis *binsy, const TAxis *binsz)
 {
     const Int_t nx = binsx->GetNbins();
@@ -281 +283 @@
 void MH::SetBinning(TH1 *h, TH1 *x)
 {
-    SetBinning(h, x->GetXaxis(), x->GetYaxis(), x->GetZaxis());
-}
+    if (h->InheritsFrom(TH3::Class()) && x->InheritsFrom(TH3::Class()))
+    {
+        SetBinning((TH3*)h, x->GetXaxis(), x->GetYaxis(), x->GetZaxis());
+        return;
+    }
+    if (h->InheritsFrom(TH2::Class()) && x->InheritsFrom(TH2::Class()))
+    {
+        SetBinning((TH2*)h, x->GetXaxis(), x->GetYaxis());
+        return;
+    }
+    if (h->InheritsFrom(TH1::Class()) && x->InheritsFrom(TH1::Class()))
+    {
+        SetBinning(h, x->GetXaxis());
+        return;
+    }
+}

trunk/MagicSoft/Mars/mhist/MH.h

@@ -7 +7 @@
 
 class TH1;
+class TH2;
+class TH3;
 class TAxis;
 class TArrayD;
@@ -28 +30 @@
 
     static void SetBinning(TH1 *h, const MBinning *binsx);
-    static void SetBinning(TH1 *h, const MBinning *binsx, const MBinning *binsy);
-    static void SetBinning(TH1 *h, const MBinning *binsx, const MBinning *binsy, const MBinning *binsz);
+    static void SetBinning(TH2 *h, const MBinning *binsx, const MBinning *binsy);
+    static void SetBinning(TH3 *h, const MBinning *binsx, const MBinning *binsy, const MBinning *binsz);
 
     static void SetBinning(TH1 *h, const TArrayD *binsx);
-    static void SetBinning(TH1 *h, const TArrayD *binsx, const TArrayD *binsy);
-    static void SetBinning(TH1 *h, const TArrayD *binsx, const TArrayD *binsy, const TArrayD *binsz);
+    static void SetBinning(TH2 *h, const TArrayD *binsx, const TArrayD *binsy);
+    static void SetBinning(TH3 *h, const TArrayD *binsx, const TArrayD *binsy, const TArrayD *binsz);
 
     static void SetBinning(TH1 *h, const TAxis *binsx);
-    static void SetBinning(TH1 *h, const TAxis *binsx, const TAxis *binsy);
-    static void SetBinning(TH1 *h, const TAxis *binsx, const TAxis *binsy, const TAxis *binsz);
+    static void SetBinning(TH2 *h, const TAxis *binsx, const TAxis *binsy);
+    static void SetBinning(TH3 *h, const TAxis *binsx, const TAxis *binsy, const TAxis *binsz);
 
     static void SetBinning(TH1 *h, TH1 *x);

trunk/MagicSoft/Mars/mhist/MH3.cc

@@ -229 +229 @@
     case 2:
         fHist->SetTitle(title+" (2D)");
-        SetBinning(fHist, binsx, binsy);
+        SetBinning((TH2*)fHist, binsx, binsy);
         return kTRUE;
     case 3:
         fHist->SetTitle(title+" (3D)");
-        SetBinning(fHist, binsx, binsy, binsz);
+        SetBinning((TH3*)fHist, binsx, binsy, binsz);
         return kTRUE;
     }

trunk/MagicSoft/Mars/mhist/MHMcCollectionArea.cc

@@ -196 +196 @@
 }
 
+// --------------------------------------------------------------------------
+//
+//  Calculate the Efficiency (collection area) and set the 'ReadyToSave'
+//  flag
+//
+void MHMcCollectionArea::CalcEfficiency()
+{
+    MHMcEfficiency heff;
+    heff.Calc(*fHistSel, *fHistAll);
+
+    Calc(heff);
+}
+
+void MHMcCollectionArea::CalcEfficiency(UInt_t numevts, Float_t angle)
+{
+    // Here we estimate the total number of showers in each energy bin
+    // known the energy range and spectral index of the generated showers
+    // (set now by hand since the info is not available in the header!)
+
+    TH1D &histsel = *fHistSel->ProjectionX();
+
+    TH1D histall;
+
+    TAxis &xaxis = *histsel.GetXaxis();
+
+    MH::SetBinning(&histall, &xaxis);
+    MH::SetBinning(fHistCol, &xaxis);
+
+    const Float_t emin  = 10.;
+    const Float_t emax  = 30000.;  // Energies in GeV.
+    const Float_t index = 2.6;     // Differential spectral Index.
+
+    const Float_t expo = 1.-index;
+
+    const Float_t k = (Float_t)numevts / (pow(emax,expo) - pow(emin,expo));
+
+    const Int_t nbinx = xaxis.GetNbins();
+
+    for (Int_t i=1; i<=nbinx; i++)
+    {
+        const Float_t e1 = histall.GetBinLowEdge(i);
+        const Float_t e2 = histall.GetBinLowEdge(i+1);
+
+        if (e1 < emin || e2 > emax)
+            continue;
+
+        const Float_t events = k * (pow(e2, expo) - pow(e1, expo));
+
+        histall.SetBinContent(i, events);
+        histall.SetBinError(i, sqrt(events));
+
+    }
+
+    // -----------------------------------------------------------
+
+    // Impact parameter range:
+    const Float_t r1 = 0;
+    const Float_t r2 = 400;
+
+    const Float_t dr = TMath::Pi() * (r2*r2 - r1*r1);
+
+    angle *= TMath::Pi()/180;
+
+    for (Int_t ix=1; ix<=nbinx; ix++)
+    {
+        const Float_t Na = histall.GetBinContent(ix);
+
+        if (Na <= 0)
+            continue;
+
+        const Float_t Ns = histsel.GetBinContent(ix);
+
+        // Since Na is an estimate of the total number of showers generated
+        // in the energy bin, it may happen that Ns (triggered showers) is
+        // larger than Na. In that case, the bin is skipped:
+
+        if (Na < Ns)
+            continue;
+
+        const Double_t eff = Ns/Na;
+
+        const Double_t err = sqrt((1.-eff)*Ns)/Na;
+
+        const Float_t area = dr * cos(angle);
+
+        fHistCol->SetBinContent(ix, eff*area);
+        fHistCol->SetBinError(ix, err*area);
+    }
+
+    delete &histsel;
+
+    SetReadyToSave();
+}
+
+void MHMcCollectionArea::Calc(const MHMcEnergyImpact &mcsel, const MHMcEnergyImpact &mcall)
+{
+    MHMcEfficiency heff;
+    heff.Calc(*mcsel.GetHist(), *mcall.GetHist());
+
+    Calc(heff);
+}
+
 void MHMcCollectionArea::Calc(const MHMcEfficiency &heff)
 {
@@ -242 +344 @@
 //  flag
 //
-void MHMcCollectionArea::CalcEfficiency()
-{
-    MHMcEfficiency heff;
-    heff.Calc(*fHistSel, *fHistAll);
-
-    Calc(heff);
-}
-
-void MHMcCollectionArea::Calc(const MHMcEnergyImpact &mcsel, const MHMcEnergyImpact &mcall)
-{
-    MHMcEfficiency heff;
-    heff.Calc(*mcsel.GetHist(), *mcall.GetHist());
-
-    Calc(heff);
-}
+/*
+void MHMcCollectionArea::Calc(MHMcEnergyImpact &mcsel, UInt_t numevents, Float_t angle)
+{
+}
+*/

trunk/MagicSoft/Mars/mhist/MHMcCollectionArea.h

@@ -38 +38 @@
 
     void CalcEfficiency();
+    void CalcEfficiency(UInt_t allevts, Float_t theta);
 
     void Calc(const MHMcEnergyImpact &mcsel, const MHMcEnergyImpact &mcall);