Index: /trunk/MagicSoft/Mars/Changelog =================================================================== --- /trunk/MagicSoft/Mars/Changelog (revision 3628) +++ /trunk/MagicSoft/Mars/Changelog (revision 3629) @@ -38,7 +38,14 @@ * mcalib/MHCalibrationChargeCam.[h,cc] - updated and enlarged documentation + - derives now from mcalib/MHCalibrationCam * mcalib/MCalibrationChargePix.[h,cc] - removed flag kLoGainSaturation, because it is already in MBadPixelsPix + + * mcalib/MCalibrationCam.[h,cc] + * mcalib/Makefile + * mcalib/CalibLinkDef.h + - new base class for all MHCalibration*Cam object, contains the average + pixel functionality Index: /trunk/MagicSoft/Mars/mcalib/CalibLinkDef.h =================================================================== --- /trunk/MagicSoft/Mars/mcalib/CalibLinkDef.h (revision 3628) +++ /trunk/MagicSoft/Mars/mcalib/CalibLinkDef.h (revision 3629) @@ -17,4 +17,5 @@ #pragma link C++ class MCalibrationChargePINDiode+; +#pragma link C++ class MHCalibrationCam+; #pragma link C++ class MHCalibrationChargeCam+; #pragma link C++ class MHCalibrationChargePix+; Index: /trunk/MagicSoft/Mars/mcalib/MHCalibrationCam.cc =================================================================== --- /trunk/MagicSoft/Mars/mcalib/MHCalibrationCam.cc (revision 3629) +++ /trunk/MagicSoft/Mars/mcalib/MHCalibrationCam.cc (revision 3629) @@ -0,0 +1,386 @@ +/* ======================================================================== *\ +! +! * +! * 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): Markus Gaug 02/2004 +! +! Copyright: MAGIC Software Development, 2000-2004 +! +! +\* ======================================================================== */ +///////////////////////////////////////////////////////////////////////////// +// +// MHCalibrationCam +// +// +///////////////////////////////////////////////////////////////////////////// +#include "MHCalibrationCam.h" + +#include +#include +#include +#include +#include + +#include "MLog.h" +#include "MLogManip.h" + +#include "MHGausEvents.h" + +#include "MGeomCam.h" +#include "MGeomPix.h" + +ClassImp(MHCalibrationCam); + +using namespace std; + +const Int_t MHCalibrationCam::fgPulserFrequency = 500; +// -------------------------------------------------------------------------- +// +// Default Constructor. +// +// Initializes and sets owner of: +// - fHiGainArray, fLoGainArray +// - fAverageHiGainAreas, fAverageLoGainAreas +// - fAverageHiGainSectors, fAverageLoGainSectors +// +// Initializes: +// - fPulserFrequency to fgPulserFrequency +// +MHCalibrationCam::MHCalibrationCam(const char *name, const char *title) +{ + fName = name ? name : "MHCalibrationCam"; + fTitle = title ? title : "Class to fill the calibration histograms "; + + fHiGainArray = new TObjArray; + fHiGainArray->SetOwner(); + + fLoGainArray = new TObjArray; + fLoGainArray->SetOwner(); + + fAverageHiGainAreas = new TObjArray; + fAverageHiGainAreas->SetOwner(); + + fAverageLoGainAreas = new TObjArray; + fAverageLoGainAreas->SetOwner(); + + fAverageHiGainSectors = new TObjArray; + fAverageHiGainSectors->SetOwner(); + + fAverageLoGainSectors = new TObjArray; + fAverageLoGainSectors->SetOwner(); + + SetPulserFrequency(); +} + +// -------------------------------------------------------------------------- +// +// Deletes the TClonesArray of: +// - fHiGainArray, fLoGainArray +// - fAverageHiGainAreas, fAverageLoGainAreas +// - fAverageHiGainSectors, fAverageLoGainSectors +// +MHCalibrationCam::~MHCalibrationCam() +{ + delete fHiGainArray; + delete fLoGainArray; + + delete fAverageHiGainAreas; + delete fAverageLoGainAreas; + + delete fAverageHiGainSectors; + delete fAverageLoGainSectors; +} + +// -------------------------------------------------------------------------- +// +// Get i-th High Gain pixel (pixel number) +// +MHGausEvents &MHCalibrationCam::operator[](UInt_t i) +{ + return *static_cast(fHiGainArray->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th High Gain pixel (pixel number) +// +const MHGausEvents &MHCalibrationCam::operator[](UInt_t i) const +{ + return *static_cast(fHiGainArray->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th Low Gain pixel (pixel number) +// +MHGausEvents &MHCalibrationCam::operator()(UInt_t i) +{ + return *static_cast(fLoGainArray->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th Low Gain pixel (pixel number) +// +const MHGausEvents &MHCalibrationCam::operator()(UInt_t i) const +{ + return *static_cast(fLoGainArray->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th High Gain pixel Area (area number) +// +MHGausEvents &MHCalibrationCam::GetAverageHiGainArea(UInt_t i) +{ + return *static_cast(fAverageHiGainAreas->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th High Gain pixel Area (area number) +// +const MHGausEvents &MHCalibrationCam::GetAverageHiGainArea(UInt_t i) const +{ + return *static_cast(fAverageHiGainAreas->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th Low Gain pixel Area (area number) +// +MHGausEvents &MHCalibrationCam::GetAverageLoGainArea(UInt_t i) +{ + return *static_cast(fAverageLoGainAreas->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th Low Gain pixel Area (area number) +// +const MHGausEvents &MHCalibrationCam::GetAverageLoGainArea(UInt_t i) const +{ + return *static_cast(fAverageLoGainAreas->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th High Gain Sector (sector number) +// +MHGausEvents &MHCalibrationCam::GetAverageHiGainSector(UInt_t i) +{ + return *static_cast(fAverageHiGainSectors->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th High Gain Sector (sector number) +// +const MHGausEvents &MHCalibrationCam::GetAverageHiGainSector(UInt_t i) const +{ + return *static_cast(fAverageHiGainSectors->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th Low Gain Sector (sector number) +// +MHGausEvents &MHCalibrationCam::GetAverageLoGainSector(UInt_t i) +{ + return *static_cast(fAverageLoGainSectors->UncheckedAt(i)); +} + +// -------------------------------------------------------------------------- +// +// Get i-th Low Gain Sector (sector number) +// +const MHGausEvents &MHCalibrationCam::GetAverageLoGainSector(UInt_t i) const +{ + return *static_cast(fAverageLoGainSectors->UncheckedAt(i)); +} + + +// -------------------------------------------------------------------------- +// +// Our own clone function is necessary since root 3.01/06 or Mars 0.4 +// I don't know the reason. +// +// Creates new MHCalibrationCam +// Deletes the TObjArray's and Clones them individually +// Copies the TArray's +// +TObject *MHCalibrationCam::Clone(const char *) const +{ + + const Int_t nhi = fHiGainArray->GetEntries(); + const Int_t nlo = fLoGainArray->GetEntries(); + const Int_t navhi = fAverageHiGainAreas->GetEntries(); + const Int_t navlo = fAverageLoGainAreas->GetEntries(); + const Int_t nsehi = fAverageHiGainSectors->GetEntries(); + const Int_t nselo = fAverageLoGainSectors->GetEntries(); + + // + // FIXME, this might be done faster and more elegant, by direct copy. + // + MHCalibrationCam *cam = new MHCalibrationCam(); + + cam->fHiGainArray->Expand(nhi); + cam->fLoGainArray->Expand(nlo); + cam->fAverageHiGainAreas->Expand(navhi); + cam->fAverageLoGainAreas->Expand(navlo); + cam->fAverageHiGainSectors->Expand(nsehi); + cam->fAverageLoGainSectors->Expand(nselo); + + for (int i=0; ifHiGainArray)[i]; + (*cam->fHiGainArray)[i] = (*fHiGainArray)[i]->Clone(); + } + for (int i=0; ifLoGainArray)[i]; + (*cam->fLoGainArray)[i] = (*fLoGainArray)[i]->Clone(); + } + for (int i=0; ifAverageHiGainAreas)[i]; + (*cam->fAverageHiGainAreas)[i] = (*fAverageHiGainAreas)[i]->Clone(); + } + for (int i=0; ifAverageLoGainAreas)[i]; + (*cam->fAverageLoGainAreas)[i] = (*fAverageLoGainAreas)[i]->Clone(); + } + for (int i=0; ifAverageHiGainSectors)[i]; + (*cam->fAverageHiGainSectors)[i] = (*fAverageHiGainSectors)[i]->Clone(); + } + for (int i=0; ifAverageLoGainSectors)[i]; + (*cam->fAverageLoGainSectors)[i] = (*fAverageLoGainSectors)[i]->Clone(); + } + + cam->fAverageAreaNum = fAverageAreaNum; + cam->fAverageAreaSat = fAverageAreaSat; + cam->fAverageAreaSigma = fAverageAreaSigma; + cam->fAverageAreaSigmaErr = fAverageAreaSigmaErr; + cam->fAverageAreaRelSigma = fAverageAreaRelSigma; + cam->fAverageAreaRelSigmaErr = fAverageAreaRelSigmaErr; + cam->fAverageSectorNum = fAverageSectorNum; + + return cam; +} + +// -------------------------------------------------------------------------- +// +// Dummy, needed by MCamEvent +// +Bool_t MHCalibrationCam::GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type) const +{ + return kTRUE; +} + +// -------------------------------------------------------------------------- +// +// What MHCamera needs in order to draw an individual pixel in the camera +// +void MHCalibrationCam::DrawPixelContent(Int_t idx) const +{ +} + +// ----------------------------------------------------------------------------- +// +// Default draw: +// +// Displays the averaged areas, both High Gain and Low Gain +// +// The following options can be chosen: +// +// Calls the Draw of the fAverageHiGainAreas and fAverageLoGainAreas objects with options +// +void MHCalibrationCam::Draw(const Option_t *opt) +{ + + const Int_t nareas = fAverageHiGainAreas->GetEntries(); + if (nareas == 0) + return; + + TVirtualPad *pad = gPad ? gPad : MH::MakeDefCanvas(this); + pad->SetBorderMode(0); + + pad->Divide(2,nareas); + + for (Int_t i=0; icd(2*(i+1)-1); + GetAverageHiGainArea(i).Draw(opt); + + if (!fAverageAreaSat[i]) + DrawAverageSigma(fAverageAreaSat[i], i, + fAverageAreaSigma[i], fAverageAreaSigmaErr[i], + fAverageAreaRelSigma[i], fAverageAreaRelSigmaErr[i]); + + pad->cd(2*(i+1)); + GetAverageLoGainArea(i).Draw(opt); + + if (fAverageAreaSat[i]) + DrawAverageSigma(fAverageAreaSat[i], i, + fAverageAreaSigma[i], fAverageAreaSigmaErr[i], + fAverageAreaRelSigma[i], fAverageAreaRelSigmaErr[i]); + } +} + +// ----------------------------------------------------------------------------- +// +// Default draw: +// +// Displays a TPaveText with the re-normalized sigmas of the average area +// +void MHCalibrationCam::DrawAverageSigma(Bool_t sat, Bool_t inner, + Float_t sigma, Float_t sigmaerr, + Float_t relsigma, Float_t relsigmaerr) const +{ + + if (sigma != 0) + { + + TPad *newpad = new TPad("newpad","transparent",0,0,1,1); + newpad->SetFillStyle(4000); + newpad->Draw(); + newpad->cd(); + + TPaveText *text = new TPaveText(sat? 0.1 : 0.35,0.7,sat ? 0.4 : 0.7,1.0); + text->SetTextSize(0.07); + const TString line1 = Form("%s%s%s",inner ? "Outer" : "Inner", + " Pixels ", sat ? "Low Gain" : "High Gain"); + TText *txt1 = text->AddText(line1.Data()); + const TString line2 = Form("#sigma per pix: %2.2f #pm %2.2f",sigma,sigmaerr); + TText *txt2 = text->AddText(line2.Data()); + const TString line3 = Form("Rel. #sigma per pix: %2.2f #pm %2.2f",relsigma,relsigmaerr); + TText *txt3 = text->AddText(line3.Data()); + text->Draw(""); + + text->SetBit(kCanDelete); + txt1->SetBit(kCanDelete); + txt2->SetBit(kCanDelete); + txt3->SetBit(kCanDelete); + newpad->SetBit(kCanDelete); + } +} + Index: /trunk/MagicSoft/Mars/mcalib/MHCalibrationCam.h =================================================================== --- /trunk/MagicSoft/Mars/mcalib/MHCalibrationCam.h (revision 3629) +++ /trunk/MagicSoft/Mars/mcalib/MHCalibrationCam.h (revision 3629) @@ -0,0 +1,100 @@ +#ifndef MARS_MHCalibrationCam +#define MARS_MHCalibrationCam + +#ifndef ROOT_TObjArray +#include +#endif + +#ifndef ROOT_TArrayI +#include +#endif + +#ifndef ROOT_TArrayF +#include +#endif + +#ifndef MARS_MH +#include "MH.h" +#endif +#ifndef MARS_MCamEvent +#include "MCamEvent.h" +#endif + +class TText; +class TArrayI; +class TArrayF; +class MHGausEvents; +class MHCalibrationCam : public MH, public MCamEvent +{ + +protected: + + static const Int_t fgPulserFrequency; // The default for fPulserFrequency + Int_t fPulserFrequency; // Light pulser frequency + + TObjArray *fHiGainArray; //-> Array of calibration pixels, one per pixel + TObjArray *fLoGainArray; //-> Array of calibration pixels, one per pixel + TObjArray *fAverageHiGainAreas; //-> Array of calibration pixels, one per pixel area + TObjArray *fAverageLoGainAreas; //-> Array of calibration pixels, one per pixel area + TObjArray *fAverageHiGainSectors; //-> Array of calibration pixels, one per camera sector + TObjArray *fAverageLoGainSectors; //-> Array of calibration pixels, one per camera sector + + TArrayI fAverageAreaNum; // Number of pixels in average pixels per area + TArrayI fAverageAreaSat; // Number of saturated slices in average pixels per area + TArrayF fAverageAreaSigma; // Re-normalized sigmas in average pixels per area + TArrayF fAverageAreaSigmaErr; // Errors of Re-normalized sigmas in average pixels per area + TArrayF fAverageAreaRelSigma; // Re-normalized relative sigmas in average pixels per area + TArrayF fAverageAreaRelSigmaErr; // Errors of Re-normalized relative sigmas in average pixels per area + TArrayI fAverageSectorNum; // Number of pixels in average pixels per sector + + void DrawAverageSigma(Bool_t sat, Bool_t inner, + Float_t sigma, Float_t sigmaerr, + Float_t relsigma, Float_t relsigmaerr) const; + +public: + + MHCalibrationCam(const char *name=NULL, const char *title=NULL); + ~MHCalibrationCam(); + + void SetPulserFrequency(const Int_t f=fgPulserFrequency) { fPulserFrequency = f; } + + MHGausEvents &operator[](UInt_t i); + const MHGausEvents &operator[](UInt_t i) const; + + MHGausEvents &operator()(UInt_t i); + const MHGausEvents &operator()(UInt_t i) const; + + MHGausEvents &GetAverageHiGainArea(UInt_t i); + const MHGausEvents &GetAverageHiGainArea(UInt_t i) const; + + MHGausEvents &GetAverageLoGainArea(UInt_t i); + const MHGausEvents &GetAverageLoGainArea(UInt_t i) const; + + MHGausEvents &GetAverageHiGainSector(UInt_t i); + const MHGausEvents &GetAverageHiGainSector(UInt_t i) const; + + MHGausEvents &GetAverageLoGainSector(UInt_t i); + const MHGausEvents &GetAverageLoGainSector(UInt_t i) const; + + // Clone + TObject *Clone(const char *) const; + + // Draw + virtual void Draw(const Option_t *opt); + + virtual Bool_t GetPixelContent ( Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type=0) const; + virtual void DrawPixelContent( Int_t num ) const; + + ClassDef(MHCalibrationCam, 1) // Base Histogram class for camera calibration +}; + +#endif + + + + + + + + + Index: /trunk/MagicSoft/Mars/mcalib/Makefile =================================================================== --- /trunk/MagicSoft/Mars/mcalib/Makefile (revision 3628) +++ /trunk/MagicSoft/Mars/mcalib/Makefile (revision 3629) @@ -49,4 +49,5 @@ MHCalibrationChargeBlindPix.cc \ MHCalibrationChargePix.cc \ + MHCalibrationCam.cc \ MHCalibrationChargeCam.cc \ MHCalibrationChargeHiGainPix.cc \