Changeset 1415 for trunk/MagicSoft/Mars/mhist/MHAlphaEnergyTime.cc

Timestamp:

07/16/02 15:35:15 (22 years ago)

Author:

wittek

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHAlphaEnergyTime.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/mhist/MHAlphaEnergyTime.cc

@@ -68 +68 @@
     fHist.SetTitle("3D-plot of alpha, E-est, time");
     fHist.SetXTitle("\\alpha [\\circ]");
-    fHist.SetYTitle("E_{est} [GeV]");
+    fHist.SetYTitle("E-est [GeV]            ");
     fHist.SetZTitle("time [s]");
 }
@@ -81 +81 @@
    if (!fEnergy)
    {
-       *fLog << err << dbginf << "MEnergyEst not found... aborting." << endl;
+       *fLog << err << dbginf << "MHAlphaEnergyTime : MEnergyEst not found... aborting." << endl;
        return kFALSE;
    }
@@ -88 +88 @@
    if (!fTime)
    {
-       *fLog << err << dbginf << "MTime not found... aborting." << endl;
+       *fLog << err << dbginf << "MHAlphaEnergyTime : MTime not found... aborting." << endl;
        return kFALSE;
    }
 
    MBinning* binsenergy = (MBinning*)plist->FindObject("BinningE");
-   MBinning* binsalpha  = (MBinning*)plist->FindObject("BinningAlpha");
+   MBinning* binsalphaflux  = (MBinning*)plist->FindObject("BinningAlphaFlux");
    MBinning* binstime   = (MBinning*)plist->FindObject("BinningTime");
-   if (!binsenergy || !binsalpha || !binstime)
+   if (!binsenergy || !binsalphaflux || !binstime)
    {
-       *fLog << err << dbginf << "At least one MBinning not found... aborting." << endl;
+       *fLog << err << dbginf << "MHAlphaEnergyTime : At least one MBinning not found... aborting." << endl;
        return kFALSE;      
    }
 
-   SetBinning(&fHist, binsalpha, binsenergy, binstime);
+   SetBinning(&fHist, binsalphaflux, binsenergy, binstime);
 
    fHist.Sumw2();
@@ -116 +116 @@
     MHillasSrc &hil = *(MHillasSrc*)par;
 
-    fHist.Fill(hil.GetAlpha(), fEnergy->GetEnergy(), 0.0001*fTime->GetTimeLo());
+    fHist.Fill(fabs(hil.GetAlpha()), fEnergy->GetEnergy(), 0.0001*fTime->GetTimeLo());
     return kTRUE;
 }
@@ -147 +147 @@
 
     h->SetTitle("Distribution of E-est [GeV]");
-    h->SetXTitle("E_{est} [GeV]");
+    h->SetXTitle("E-est [GeV]            ");
     h->SetYTitle("Counts");
 
@@ -226 +226 @@
 }
 
-// --------------------------------------------------------------------------
-//
-// Calculate the histogram as the difference of two histograms :
-//          fHist(gamma) = h1(source) - h2(antisource)
-// 
-void MHAlphaEnergyTime::Subtract(const TH3D *h1, const TH3D *h2)
-{
-  //    MH::SetBinning(&fHist, (TH1*)h1);
-
-  //    fHist.Sumw2();
-    fHist.Add((TH1*)h1, (TH1*)h2, 1, -1); // ROOT: FIXME!
-}
-
-// --------------------------------------------------------------------------
-//
-// Integrate fHist(gamma) in the alpha range (lo, up)
-// 
-TH2D *MHAlphaEnergyTime::GetAlphaProjection(Axis_t lo, Axis_t up)
-{
-    if (up < lo)
+
+// --------------------------------------------------------------------------
+//
+// Integrate fHist     (Alpha,E-est,Time) over the Time to get
+//           fAlphaEest(Alpha,E-est)
+// 
+TH2D *MHAlphaEnergyTime::IntegrateTime(const char *title, Bool_t draw)
+{
+    Int_t nzbins = fHist.GetNbinsZ();
+    TAxis &axez  = *fHist.GetZaxis();
+    axez.SetRange(1,nzbins);
+
+    TH2D &fAlphaEest = *(TH2D *)fHist.Project3D("exy");
+
+    fAlphaEest.SetTitle(title);
+    fAlphaEest.SetXTitle("E-est [GeV]            ");
+    fAlphaEest.SetYTitle("\\alpha  [  \\circ]");
+
+    if (draw == kTRUE)
     {
-        *fLog << err << fName << ": Alpha projection not possible: lo=" << lo << " up=" << up << endl;
-        return NULL;
+      TCanvas &c = *MakeDefCanvas(title, title);
+
+      gROOT->SetSelectedPad(NULL);
+
+      fAlphaEest.DrawCopy();
+      gPad->SetLogx();
+
+      c.Modified();
+      c.Update();
     }
 
-    TAxis &axe = *fHist.GetXaxis();
-
-    Int_t ilo = axe.FindFixBin(lo);
-    Int_t iup = axe.FindFixBin(up);
-
-    const Double_t epslo1 = lo-axe.GetBinLowEdge(ilo);
-    const Double_t epslo2 = axe.GetBinUpEdge(ilo)-lo;
-
-    const Double_t epsup1 = up-axe.GetBinLowEdge(iup);
-    const Double_t epsup2 = axe.GetBinUpEdge(iup)-up;
-
-    const Double_t epslo = epslo1<epslo2 ? epslo1 : epslo2;
-    const Double_t epsup = epsup1<epsup2 ? epsup1 : epsup2;
-
-    if (epslo1>epslo2)
-        ilo++;
-
-    if (epsup1<epsup2)
-        iup--;
-
-    if (epslo>0.01*axe.GetBinWidth(ilo) || epsup>0.01*axe.GetBinWidth(iup))
+    return &fAlphaEest;
+}
+
+// --------------------------------------------------------------------------
+//
+// Integrate fHist     (Alpha,E-est,Time) over E-est to get
+//           fAlphaTime(Alpha,Time)
+// 
+TH2D *MHAlphaEnergyTime::IntegrateEest(const char *title, Bool_t draw)
+{
+    Int_t nybins = fHist.GetNbinsY();
+    TAxis &axey  = *fHist.GetYaxis();
+    axey.SetRange(1,nybins);
+
+    TH2D &fAlphaTime = *(TH2D *)fHist.Project3D("exz");
+
+    fAlphaTime.SetTitle(title);
+    fAlphaTime.SetXTitle("Time [s]");
+    fAlphaTime.SetYTitle("\\alpha  [  \\circ]");
+
+    if (draw == kTRUE)
     {
-        *fLog << err << fName << ": binning is not adequate for the requested projection:" << endl;
-        *fLog << "Please specify a lower or upper limit which is not more than 1% away from a bin edge" << endl;
-        *fLog << " epslo = " << epslo << endl;
-        *fLog << " epsup = " << epsup << endl;
-        *fLog << " dwl   = " << axe.GetBinWidth(ilo) << endl;
-        *fLog << " dwu   = " << axe.GetBinWidth(iup) << endl;
-        return NULL;
+      TCanvas &c = *MakeDefCanvas(title, title);
+
+      gROOT->SetSelectedPad(NULL);
+
+      fAlphaTime.DrawCopy();
+
+      c.Modified();
+      c.Update();
     }
 
-    axe.SetRange(ilo, iup);
-
-    TH2D &h2D = *(TH2D *)fHist.Project3D("ezy");
-
-    h2D.SetTitle("2D-plot  of time vs. E-est");
-    h2D.SetXTitle("E-est [GeV]            ");
-    h2D.SetYTitle("time [s]");
-
-    return &h2D;
-}
-
-//---------------------------------------------------------
-//
-// Draw the projected histogram
-//
-TH2D *MHAlphaEnergyTime::DrawAlphaProjection(Axis_t lo, Axis_t up, Option_t *opt)
-{
-    TH2D *h2D = GetAlphaProjection(lo, up);
-
-    if (!h2D)
-        return NULL;
-
-    char txt[100];
-    sprintf(txt, "No.of Gammas vs. E-est and Time (%.1f < alpha < %.1f deg)", lo, up);
-
-    //    TCanvas *c = MakeDefCanvas("AlphaEnergyTime", "2D histogram of gamma signal in energy and time");
-    TCanvas &c = *MakeDefCanvas("AlphaEnergyTime", txt);
-
-    c.Divide(2, 2);
-
-    gROOT->SetSelectedPad(NULL);
-
-    TH1 *h;
-
-    c.cd(1);
-    h = h2D->ProjectionX("Eest", -1, 9999, "E");
-    h->SetTitle("Distribution of E-est [GeV]");
-    h->SetXTitle("E-est [GeV]            ");
-    h->SetYTitle("Counts");
-
-    h->Draw(opt);
-    h->SetBit(kCanDelete);
-    gPad->SetLogx();
-
-    c.cd(2);
-    h = h2D->ProjectionY("time", -1, 9999, "E");
-    h->SetTitle("Distribution of time [s]");
-    h->SetXTitle("time [s]");
-    h->SetYTitle("Counts");
-
-    h->Draw(opt);
-    h->SetBit(kCanDelete);
-
-    c.cd(3);
-
-    h2D->DrawCopy(opt);
-    gPad->SetLogx();
-
-    c.Modified();
-    c.Update();
-
-    return h2D;
-}
-
+    return &fAlphaTime;
+}
+
+// --------------------------------------------------------------------------
+//
+// Integrate fHist (Alpha,E-est,Time) over Eest and Time to get
+//           fAlpha(Alpha)
+// 
+TH1D *MHAlphaEnergyTime::IntegrateEestTime(const char *title, Bool_t draw)
+{
+    Int_t nybins = fHist.GetNbinsY();
+    TAxis &axey  = *fHist.GetYaxis();
+    axey.SetRange(1,nybins);
+
+    Int_t nzbins = fHist.GetNbinsZ();
+    TAxis &axez  = *fHist.GetZaxis();
+    axez.SetRange(1,nzbins);
+
+    TH1D &fAlpha = *(TH1D *)fHist.Project3D("ex");
+
+    fAlpha.SetTitle(title);
+    fAlpha.SetXTitle("\\alpha  [  \\circ]");
+    fAlpha.SetYTitle("Counts");
+
+    if (draw == kTRUE)
+    {
+      TCanvas &c = *MakeDefCanvas(title, title);
+
+      gROOT->SetSelectedPad(NULL);
+
+      fAlpha.DrawCopy();
+
+      c.Modified();
+      c.Update();
+    }
+
+    return &fAlpha;
+}
+
+
+
+
+
+
+