Index: trunk/MagicSoft/Mars/mbadpixels/BadPixelsLinkDef.h =================================================================== --- trunk/MagicSoft/Mars/mbadpixels/BadPixelsLinkDef.h (revision 4793) +++ trunk/MagicSoft/Mars/mbadpixels/BadPixelsLinkDef.h (revision 4841) @@ -13,4 +13,5 @@ #pragma link C++ class MMcBadPixelsSet+; +#pragma link C++ class MHBadPixels+; #endif Index: trunk/MagicSoft/Mars/mbadpixels/MHBadPixels.cc =================================================================== --- trunk/MagicSoft/Mars/mbadpixels/MHBadPixels.cc (revision 4841) +++ trunk/MagicSoft/Mars/mbadpixels/MHBadPixels.cc (revision 4841) @@ -0,0 +1,201 @@ +/* ======================================================================== *\ +! +! * +! * This file is part of MARS, the MAGIC Analysis and Reconstruction +! * Software. It is distributed to you in the hope that it can be a useful +! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. +! * It is distributed WITHOUT ANY WARRANTY. +! * +! * Permission to use, copy, modify and distribute this software and its +! * documentation for any purpose is hereby granted without fee, +! * provided that the above copyright notice appear in all copies and +! * that both that copyright notice and this permission notice appear +! * in supporting documentation. It is provided "as is" without express +! * or implied warranty. +! * +! +! +! Author(s): Thomas Bretz 04/2003 +! +! Copyright: MAGIC Software Development, 2000-2003 +! +! +\* ======================================================================== */ + +///////////////////////////////////////////////////////////////////////////// +// +// MHBadPixels +// +//////////////////////////////////////////////////////////////////////////// +#include "MHBadPixels.h" + +#include + +#include "MPointingPos.h" +#include "MBadPixelsCam.h" +#include "MGeomCam.h" +#include "MPedPhotCam.h" +#include "MParList.h" +#include "MBinning.h" + +#include "MLog.h" +#include "MLogManip.h" + +ClassImp(MHBadPixels); + +using namespace std; + +// ------------------------------------------------------------------------- +// +// Default Constructor. +// +MHBadPixels::MHBadPixels(const char *name, const char *title) +{ + fName = name ? name : "MHBadPixels"; + fTitle = title ? title : "Histogram for Bad Pixels vs. Theta"; + + fBadId.SetName("2D-IdBadPixels"); + fBadId.SetTitle("2D-IdBadPixels"); + fBadId.SetDirectory(NULL); + fBadId.SetXTitle("\\Theta [\\circ]"); + fBadId.SetYTitle("pixel Id"); + + fBadN.SetName("2D-NBadPixels"); + fBadN.SetTitle("2D-NBadPixels"); + fBadN.SetDirectory(NULL); + fBadN.SetXTitle("\\Theta [\\circ]"); + fBadN.SetYTitle("number of bad pixels"); + + fNamePedPhotCam = "MPedPhotCamFromData"; +} + +// -------------------------------------------------------------------------- +// +// Set the binnings and prepare the filling of the histogram +// +Bool_t MHBadPixels::SetupFill(const MParList *plist) +{ + MGeomCam *fCam = (MGeomCam*)plist->FindObject(AddSerialNumber("MGeomCam")); + if (!fCam) + { + *fLog << err << "MHBadPixels::SetupFill; MGeomCam not found... aborting." << endl; + return kFALSE; + } + + fPointPos = (MPointingPos*)plist->FindObject("MPointingPos"); + if (!fPointPos) + { + *fLog << err << "MPointingPos not found... aborting." << endl; + return kFALSE; + } + + + fPedPhot = (MPedPhotCam*)plist->FindObject(AddSerialNumber(fNamePedPhotCam), "MPedPhotCam"); + if (!fPedPhot) + { + *fLog << err << AddSerialNumber(fNamePedPhotCam) + << "[MPedPhotCam] not found... aborting." << endl; + return kFALSE; + } + //fPedPhot->InitSize(fCam->GetNumPixels()); + + + // 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 = fPedPhot->GetSize()+1; + + MBinning binspix("BinningPixel"); + binspix.SetEdges(npix1, -0.5, npix1-0.5); + + // Set binnings in histograms + SetBinning(&fBadId, binstheta, &binspix); + SetBinning(&fBadN, binstheta, &binspix); + + return kTRUE; +} + +// ------------------------------------------------------------------------ +// +// Drawing function. It creates its own canvas. +// +void MHBadPixels::Draw(Option_t *option) +{ + TVirtualPad *pad = gPad ? gPad : MakeDefCanvas(this); + pad->SetBorderMode(0); + + pad->Divide(2,2); + + TH1D *h; + + pad->cd(1); + fBadId.Draw(option); + + pad->cd(2); + fBadN.Draw(option); + + pad->cd(3); + gPad->SetBorderMode(0); + h = ((TH2*)&fBadId)->ProjectionY("ProjY-pixId", -1, 9999, ""); + h->SetDirectory(NULL); + h->SetTitle("Distribution of bad pixel Id"); + h->SetXTitle("Id of bad pixel"); + h->SetYTitle("No. of events"); + h->Draw(option); + h->SetBit(kCanDelete); + + pad->cd(4); + gPad->SetBorderMode(0); + h = ((TH2*)&fBadN)->ProjectionY("ProjY-pixN", -1, 9999, ""); + h->SetDirectory(NULL); + h->SetTitle("Distribution of no.of bad pixels"); + h->SetXTitle("No. of bad pixels"); + h->SetYTitle("No. of events"); + h->Draw(option); + h->SetBit(kCanDelete); + + pad->Modified(); + pad->Update(); +} + +Bool_t MHBadPixels::Fill(const MParContainer *par, const Stat_t w) +{ + if (!par) + return kFALSE; + + Double_t theta = fPointPos->GetZd(); + + const MBadPixelsCam *fBadPixels = (MBadPixelsCam*)par; + + UShort_t entries = fCam->GetNumPixels(); + UInt_t nb = 0; + for (UInt_t i=0; i +#endif + +class MPedPhotCam; +class MPointingPos; +class MParList; +class MGeomCam; + +class MHBadPixels : public MH +{ +private: + MGeomCam *fCam; //! + MPedPhotCam *fPedPhot; //! + MPointingPos *fPointPos; //! + + TString fNamePedPhotCam; //! name of the 'MPedPhotCam' container + + TH2D fBadId; // 2D-histogram : pixel Id vs. Theta + TH2D fBadN; // 2D-histogram : no.of bad pixels vs. Theta + +public: + MHBadPixels(const char *name=NULL, const char *title=NULL); + + void SetNamePedPhotCam(const char *name) { fNamePedPhotCam = name; } + + const TH2D *GetBadId() { return &fBadId; } + const TH2D *GetBadId() const { return &fBadId; } + + const TH2D *GetBadN() { return &fBadN; } + const TH2D *GetBadN() const { return &fBadN; } + + TH2 *GetBadIdByName(const TString name) { return &fBadId; } + TH2 *GetBadNByName(const TString name) { return &fBadN; } + + void Draw(Option_t* option = ""); + Bool_t SetupFill(const MParList *plist); + Bool_t Fill(const MParContainer *par, const Stat_t w=1); + + ClassDef(MHBadPixels, 1) // Histogram of bad pixel Id vs. Theta +}; + +#endif + + Index: trunk/MagicSoft/Mars/mbadpixels/Makefile =================================================================== --- trunk/MagicSoft/Mars/mbadpixels/Makefile (revision 4793) +++ trunk/MagicSoft/Mars/mbadpixels/Makefile (revision 4841) @@ -20,5 +20,6 @@ # INCLUDES = -I. -I../mbase -I../manalysis -I../mgeom -I../mgui \ - -I../mmc -I../mpedestal -I../msignal + -I../mmc -I../mpedestal -I../msignal -I../mpointing \ + -I../mhbase # MBadPixelsCalc (manalysis): MPedPhotCam, MSigmabar # MBadPixelsCalc (mgeom): MGeomCam @@ -35,5 +36,6 @@ MBadPixelsCalc.cc \ MBadPixelsTreat.cc \ - MMcBadPixelsSet.cc + MMcBadPixelsSet.cc \ + MHBadPixels.cc ############################################################