Changeset 2128 for trunk/MagicSoft/Mars/mhist/MHBlindPixels.cc

Timestamp:

05/21/03 12:21:58 (21 years ago)

Author:

wittek

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHBlindPixels.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/mhist/MHBlindPixels.cc

@@ +1 @@
+
 /* ======================================================================== *\
 !
@@ -32 +33 @@
 #include <TCanvas.h>
 
+#include "MMcEvt.hxx"
 #include "MBlindPixels.h"
+#include "MPedestalCam.h"
+#include "MParList.h"
+#include "MBinning.h"
+
+#include "MLog.h"
+#include "MLogManip.h"
 
 ClassImp(MHBlindPixels);
@@ -41 +49 @@
 //
 MHBlindPixels::MHBlindPixels(const char *name, const char *title)
-    : fHist(fName, fTitle, 577, -.5, 577-.5)
 {
     fName  = name  ? name  : "MHBlindPixels";
-    fTitle = title ? title : "Histogram for Blind Pixels";
+    fTitle = title ? title : "Histogram for Blind Pixels vs. Theta";
 
     //  - we initialize the histogram
-    //  - we have to give diferent names for the diferent histograms because
+    //  - we have to give different names for the different histograms because
     //    root don't allow us to have diferent histograms with the same name
 
-    fHist.SetName("1D-BlindPixels");
-    fHist.SetTitle(fTitle);
+    fBlindId.SetName("2D-IdBlindPixels");
+    fBlindId.SetTitle("2D-IdBlindPixels");
+    fBlindId.SetDirectory(NULL);
+    fBlindId.SetXTitle("\\Theta [\\circ]");
+    fBlindId.SetYTitle("pixel Id");
 
-    fHist.SetDirectory(NULL);
+    fBlindN.SetName("2D-NBlindPixels");
+    fBlindN.SetTitle("2D-NBlindPixels");
+    fBlindN.SetDirectory(NULL);
+    fBlindN.SetXTitle("\\Theta [\\circ]");
+    fBlindN.SetYTitle("number of blind pixels");
+}
 
-    fHist.SetXTitle("Id");
-    fHist.SetYTitle("Counts");
+// --------------------------------------------------------------------------
+//
+// Set the binnings and prepare the filling of the histogram
+//
+Bool_t MHBlindPixels::SetupFill(const MParList *plist)
+{
+    fMcEvt = (MMcEvt*)plist->FindObject("MMcEvt");
+    if (!fMcEvt)
+    {
+        *fLog << err << "MMcEvt not found... aborting." << endl;
+        return kFALSE;
+    }
+
+    fPed = (MPedestalCam*)plist->FindObject("MPedestalCam");
+    if (!fPed)
+    {
+        *fLog << err << "MPedestalCam not found... aborting." << endl;
+        return kFALSE;
+    }
+
+
+    // Get Theta Binning
+    MBinning* binstheta  = (MBinning*)plist->FindObject("BinningTheta", "MBinning");
+    if (!binstheta)
+    {
+        *fLog << err << "Object 'BinningTheta' [MBinning] not found... aborting" << endl;
+        return kFALSE;
+    }
+
+    // Get binning for pixel number
+    const UInt_t npix1 = fPed->GetSize()+1;
+
+    MBinning binspix("BinningPixel");
+    binspix.SetEdges(npix1, -0.5, npix1-0.5);
+
+    // Set binnings in histograms
+    SetBinning(&fBlindId, binstheta, &binspix);
+    SetBinning(&fBlindN,  binstheta, &binspix);
+
+    return kTRUE;
 }
 
@@ -68 +121 @@
     pad->SetBorderMode(0);
 
-    fHist.Draw(option);
+    pad->Divide(2,2);
+
+    TH1D *h;
+
+    pad->cd(1);
+    fBlindId.Draw(option);
+
+    pad->cd(2);
+    fBlindN.Draw(option);
+
+    pad->cd(3);
+    gPad->SetBorderMode(0);
+    h = ((TH2*)&fBlindId)->ProjectionY("ProjY-pixId", -1, 9999, "");
+    h->SetDirectory(NULL);
+    h->SetTitle("Distribution of blind pixel Id");
+    h->SetXTitle("Id of blind pixel");
+    h->SetYTitle("No. of events");
+    h->Draw(option);
+    h->SetBit(kCanDelete);
+
+    pad->cd(4);
+    gPad->SetBorderMode(0);
+    h = ((TH2*)&fBlindN)->ProjectionY("ProjY-pixN", -1, 9999, "");
+    h->SetDirectory(NULL);
+    h->SetTitle("Distribution of no.of blind pixels");
+    h->SetXTitle("No. of blind pixels");
+    h->SetYTitle("No. of events");
+    h->Draw(option);
+    h->SetBit(kCanDelete);
 
     pad->Modified();
@@ -79 +160 @@
         return kFALSE;
 
+    Double_t theta = fMcEvt->GetTelescopeTheta()*kRad2Deg;
     const MBlindPixels &bp = *(MBlindPixels*)par;
 
     // FIXME: Slow.
-    for (int i=0; i<577; i++)
+    const UInt_t npix = fPed->GetSize();
+
+    UInt_t nb = 0;
+    for (UInt_t i=0; i<npix; i++)
+    {
         if (bp.IsBlind(i))
-            fHist.Fill(i, w);
+        {
+          fBlindId.Fill(theta, i, w);
+          nb++;
+        }   
+    }
+    fBlindN.Fill(theta, nb, w);
 
     return kTRUE;
 }
+
+
+
+
+
+
+
+
+
+
+
+