Changeset 1652 for trunk/MagicSoft/Mars/mhist/MHFadcPix.cc

Timestamp:

11/20/02 12:18:01 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHFadcPix.cc (modified) (8 diffs)

trunk/MagicSoft/Mars/mhist/MHFadcPix.cc

@@ -16 +16 @@
 !
 !
-!   Author(s): Thomas Bretz  12/2000 <mailto:tbretz@uni-sw.gwdg.de>
+!   Author(s): Thomas Bretz  12/2000 <mailto:tbretz@astro.uni-wuerzburg.de>
 !   Author(s): Harald Kornmayer 1/2001
 !
-!   Copyright: MAGIC Software Development, 2000-2001
+!   Copyright: MAGIC Software Development, 2000-2002
 !
 !
@@ -28 +28 @@
 // MHFadcPix
 //
-// This container stores a hostogram to display an Fadc Spekrtum.
+// This container stores a histogram to display an Fadc Spektrum.
 // The spektrum of all the values measured by the Fadcs.
 //
 ///////////////////////////////////////////////////////////////////////
-
 #include "MHFadcPix.h"
 
-#include <TH1.h>
 #include <TPad.h>
 
 #include "MH.h"
+#include "MBinning.h"
+
+#include "MRawEvtData.h"
+#include "MRawEvtPixelIter.h"
 
 ClassImp(MHFadcPix);
@@ -46 +48 @@
 // Creates the histograms for lo and hi gain of one pixel
 //
-MHFadcPix::MHFadcPix(UInt_t pixid) : fPixId(pixid)
+MHFadcPix::MHFadcPix(UInt_t pixid, Type_t t)
+    : fPixId(pixid), fType(t)
 {
-    fHistHi =  new TH1F(pixid ? Form("HiGain%03d", pixid) : "HiGain",
-                        pixid ? Form("Hi Gain Pixel #%d", pixid) : "Hi Gain Samples",
-                        256, -.5, 255.5);
+    fHistHi.SetName(pixid ? Form("HiGain%03d", pixid) : "HiGain");
+    fHistHi.SetTitle(pixid ? Form("Hi Gain Pixel #%d", pixid) : "Hi Gain Samples");
+    fHistHi.SetDirectory(NULL);
 
-    fHistHi->SetDirectory(NULL);
-    fHistHi->SetXTitle("Signal/FADC Units");
-    fHistHi->SetYTitle("Count");
+    fHistLo.SetName(pixid ? Form("LoGain%03d", pixid) : "LoGain");
+    fHistLo.SetTitle(pixid ? Form("Lo Gain Pixel #%d", pixid) : "Lo Gain Samples");
+    fHistLo.SetDirectory(NULL);
 
-    fHistLo = new TH1F(pixid ? Form("LoGain%03d", pixid) : "LoGain",
-                       pixid ? Form("Lo Gain Pixel #%d", pixid) : "Lo Gain Samples",
-                       256, -.5, 255.5);
+    if (fType==kValue)
+    {
+        fHistHi.SetXTitle("Signal/FADC Units");
+        fHistLo.SetXTitle("Signal/FADC Units");
+        fHistHi.SetYTitle("Count");
+        fHistLo.SetYTitle("Count");
 
-    fHistLo->SetDirectory(NULL);
-    fHistLo->SetXTitle("Signal/FADC Units");
-    fHistLo->SetYTitle("Count");
+        MBinning bins;
+        bins.SetEdges(255, -.5, 255.5);
+
+        bins.Apply(fHistHi);
+        bins.Apply(fHistLo);
+    }
+    else
+    {
+        fHistHi.SetXTitle("Time/FADC Slices");
+        fHistLo.SetXTitle("Time/FADC Slices");
+        fHistLo.SetYTitle("Sum Signal/FADC Units");
+        fHistHi.SetYTitle("Sum Signal/FADC Units");
+    }
 }
 
-// --------------------------------------------------------------------------
-MHFadcPix::~MHFadcPix()
+void MHFadcPix::Init(Byte_t nhi, Byte_t nlo)
 {
-    delete fHistHi;
-    delete fHistLo;
+    MBinning bins;
+
+    bins.SetEdges(nhi, -.5, -.5+nhi);
+    bins.Apply(fHistHi);
+
+    bins.SetEdges(nlo, -.5, -.5+nlo);
+    bins.Apply(fHistLo);
 }
 
-// --------------------------------------------------------------------------
-void MHFadcPix::FillHi(Byte_t i)
+Bool_t MHFadcPix::Fill(const MRawEvtData &evt)
 {
-    fHistHi->Fill(i);
-}
+    MRawEvtPixelIter pixel((MRawEvtData*)&evt);
 
-// --------------------------------------------------------------------------
-void MHFadcPix::FillLo(Byte_t i)
-{
-    fHistLo->Fill(i);
+    if (!pixel.Jump(fPixId))
+        return kTRUE;
+
+    const Int_t nhisamples = evt.GetNumHiGainSamples();
+
+    if (fType==kValue)
+        for (Int_t i=0; i<nhisamples; i++)
+            fHistHi.Fill(pixel.GetHiGainSamples()[i]);
+    else
+        for (Int_t i=0; i<nhisamples; i++)
+            fHistHi.Fill(i, pixel.GetHiGainSamples()[i]);
+
+    if (!pixel.HasLoGain())
+        return kTRUE;
+
+    const Int_t nlosamples = evt.GetNumLoGainSamples();
+
+    if (fType==kValue)
+        for (Int_t i=0; i<nlosamples; i++)
+            fHistLo.Fill(pixel.GetLoGainSamples()[i]);
+    else
+        for (Int_t i=0; i<nlosamples; i++)
+            fHistLo.Fill(i, pixel.GetLoGainSamples()[i]);
+
+    return kTRUE;
 }
 
@@ -87 +126 @@
 void MHFadcPix::DrawHi()
 {
-    fHistHi->Draw();
+    fHistHi.Draw();
 }
 
@@ -93 +132 @@
 void MHFadcPix::DrawLo()
 {
-    fHistLo->Draw();
+    fHistLo.Draw();
 }
 
@@ -105 +144 @@
     MHFadcPix &pix = *(MHFadcPix*)TObject::Clone();
 
-    pix.fHistHi->SetDirectory(NULL);
-    pix.fHistLo->SetDirectory(NULL);
+    pix.fHistHi.SetDirectory(NULL);
+    pix.fHistLo.SetDirectory(NULL);
 
     return &pix;
@@ -116 +155 @@
     if (!gPad)
     {
-        const char *name  = StrDup(fPixId ? Form("Pixel #%d", fPixId) : "Pixel");
+        const char *name = StrDup(fPixId ? Form("Pixel #%d", fPixId) : "Pixel");
         MH::MakeDefCanvas(name, fPixId ? Form("%s FADC Samples", name) : "FADC Samples");
         delete [] name;
@@ -124 +163 @@
 
     gPad->cd(1);
-    fHistHi->Draw();
+    fHistHi.Draw();
 
     gPad->cd(2);
-    fHistLo->Draw();
+    fHistLo.Draw();
 }