Index: trunk/MagicSoft/Mars/manalysis/MPadSchweizer.cc =================================================================== --- trunk/MagicSoft/Mars/manalysis/MPadSchweizer.cc (revision 2019) +++ (revision ) @@ -1,794 +1,0 @@ -/* ======================================================================== *\ -! -! * -! * 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): Robert Wagner 10/2002 -! Author(s): Wolfgang Wittek 02/2003 -! -! Copyright: MAGIC Software Development, 2000-2003 -! -! -\* ======================================================================== */ - -///////////////////////////////////////////////////////////////////////////// -// -// MPadSchweizer -// -// This task applies padding such that for a given pixel and for a given -// Theta bin the resulting distribution of the pedestal sigma is identical -// to the distributions given by fHSigmaPixTheta and fHDiffPixTheta. -// -// The number of photons, its error and the pedestal sigmas are altered. -// On average, the number of photons added is zero. -// -// The formulas used can be found in Thomas Schweizer's Thesis, -// Section 2.2.1 -// -// There are 2 options for the padding : -// -// 1) fPadFlag = 1 : -// Generate first a Sigmabar using the 2D-histogram Sigmabar vs. Theta -// (fHSigmaTheta). Then generate a pedestal sigma for each pixel using -// the 3D-histogram Theta, pixel no., Sigma^2-Sigmabar^2 -// (fHDiffPixTheta). -// -// This is the preferred option as it takes into account the -// correlations between the Sigma of a pixel and Sigmabar. -// -// 2) fPadFlag = 2 : -// Generate a pedestal sigma for each pixel using the 3D-histogram -// Theta, pixel no., Sigma (fHSigmaPixTheta). -// -// -// The padding has to be done before the image cleaning because the -// image cleaning depends on the pedestal sigmas. -// -// For random numbers gRandom is used. -// -// This implementation has been tested for CT1 data. For MAGIC some -// modifications are necessary. -// -///////////////////////////////////////////////////////////////////////////// -#include "MPadSchweizer.h" - -#include - -#include -#include -#include -#include -#include - -#include "MBinning.h" -#include "MSigmabar.h" -#include "MMcEvt.hxx" -#include "MLog.h" -#include "MLogManip.h" -#include "MParList.h" -#include "MGeomCam.h" - -#include "MCerPhotPix.h" -#include "MCerPhotEvt.h" - -#include "MPedestalCam.h" -#include "MPedestalPix.h" - -ClassImp(MPadSchweizer); - -// -------------------------------------------------------------------------- -// -// Default constructor. -// -MPadSchweizer::MPadSchweizer(const char *name, const char *title) -{ - fName = name ? name : "MPadSchweizer"; - fTitle = title ? title : "Task for the padding (Schweizer)"; -} - -// -------------------------------------------------------------------------- -// -// Destructor. STLL DOESN'T DESTRUCT EVERYTHING -// -MPadSchweizer::~MPadSchweizer() -{ - *fLog << all << "WARNING: ~MPadSchweizer called: Potential Memory Leak" << endl; -} - -// -------------------------------------------------------------------------- -// -// Set the references to the histograms to be used in the padding -// DOCUMENTATION of the arguments MISSING -// -void MPadSchweizer::SetHistograms(TH2D *hist2, TH3D *hist3, TH3D *hist3Diff) -{ - fHSigmaTheta = hist2; - fHSigmaPixTheta = hist3; - fHDiffPixTheta = hist3Diff; - - fHSigmaTheta->SetDirectory(NULL); - fHSigmaPixTheta->SetDirectory(NULL); - fHDiffPixTheta->SetDirectory(NULL); - - Print(); -} - -// -------------------------------------------------------------------------- -// -// Set the option for the padding -// -// There are 2 options for the padding : -// -// 1) fPadFlag = 1 : -// Generate first a Sigmabar using the 2D-histogram Sigmabar vs. Theta -// (fHSigmaTheta). Then generate a pedestal sigma for each pixel using -// the 3D-histogram Theta, pixel no., Sigma^2-Sigmabar^2 -// (fHDiffPixTheta). -// -// This is the preferred option as it takes into account the -// correlations between the Sigma of a pixel and Sigmabar. -// -// 2) fPadFlag = 2 : -// Generate a pedestal sigma for each pixel using the 3D-histogram -// Theta, pixel no., Sigma (fHSigmaPixTheta). -// -void MPadSchweizer::SetPadFlag(Int_t padflag) -{ - fPadFlag = padflag; - *fLog << "MPadSchweizer::SetPadFlag(); choose option " << fPadFlag << endl; -} - -// -------------------------------------------------------------------------- -// -// Set the pointers and prepare the histograms -// -Bool_t MPadSchweizer::PreProcess(MParList *pList) -{ - fMcEvt = (MMcEvt*)pList->FindObject("MMcEvt"); - if (!fMcEvt) - { - *fLog << err << dbginf << "MMcEvt not found... aborting." << endl; - return kFALSE; - } - - fPed = (MPedestalCam*)pList->FindObject("MPedestalCam"); - if (!fPed) - { - *fLog << err << "MPedestalCam not found... aborting." << endl; - return kFALSE; - } - - fCam = (MGeomCam*)pList->FindObject("MGeomCam"); - if (!fCam) - { - *fLog << err << "MGeomCam not found (no geometry information available)... aborting." << endl; - return kFALSE; - } - - fEvt = (MCerPhotEvt*)pList->FindObject("MCerPhotEvt"); - if (!fEvt) - { - *fLog << err << "MCerPhotEvt not found... aborting." << endl; - return kFALSE; - } - - fSigmabar = (MSigmabar*)pList->FindCreateObj("MSigmabar"); - if (!fSigmabar) - { - *fLog << err << "MSigmabar not found... aborting." << endl; - return kFALSE; - } - - // Get Theta Binning - MBinning* binstheta = (MBinning*)pList->FindObject("BinningTheta"); - if (!binstheta) - { - *fLog << err << "BinningTheta not found... aborting." << endl; - return kFALSE; - } - - // Get Sigma Binning - MBinning* binssigma = (MBinning*)pList->FindObject("BinningSigmabar"); - if (!binssigma) - { - *fLog << err << "BinningSigmabar not found... aborting." << endl; - return kFALSE; - } - - // Get binning for (sigma^2-sigmabar^2) - MBinning* binsdiff = (MBinning*)pList->FindObject("BinningDiffsigma2"); - if (!binsdiff) - { - *fLog << err << "BinningDiffsigma2 not found... aborting." << endl; - return kFALSE; - } - - // Get binning for pixel number - UInt_t npix = fPed->GetSize(); - MBinning binspix("BinningPixel"); - MBinning* binspixel = &binspix; - binspixel->SetEdges(npix, -0.5, ((float)npix)-0.5 ); - - - //-------------------------------------------------------------------- - // histograms for checking the padding - // - // plot pedestal sigmas - fHSigmaPedestal = new TH2D("SigPed","Sigma: after vs. before padding", - 100, 0.0, 5.0, 100, 0.0, 5.0); - fHSigmaPedestal->SetXTitle("orig. Pedestal sigma"); - fHSigmaPedestal->SetYTitle("padded Pedestal sigma"); - - // plot no.of photons (before vs. after padding) - fHPhotons = new TH2D("Photons","Photons: after vs.before padding", - 100, -10.0, 90.0, 100, -10, 90); - fHPhotons->SetXTitle("no.of photons before padding"); - fHPhotons->SetYTitle("no.of photons after padding"); - - // plot of added NSB - fHNSB = new TH1D("NSB","Distribution of added NSB", 100, -10.0, 10.0); - fHNSB->SetXTitle("no.of added NSB photons"); - fHNSB->SetYTitle("no.of pixels"); - - fHSigbarTheta = new TH2D(); - fHSigbarTheta->SetNameTitle("fHSigmabar", - "2D : Sigmabar, \\Theta (after padding)"); - MH::SetBinning(fHSigbarTheta, binstheta, binssigma); - fHSigbarTheta->SetXTitle("Theta"); - fHSigbarTheta->SetYTitle("Sigmabar"); - - fHSigPixTh = new TH3D(); - fHSigPixTh->SetNameTitle("fHSigPixTh","3D: \\Theta, pixel no., Sigma (after padding)"); - MH::SetBinning(fHSigPixTh, binstheta, binspixel, binssigma); - fHSigPixTh->SetXTitle("\\Theta [\\circ]"); - fHSigPixTh->SetYTitle("pixel number"); - fHSigPixTh->SetZTitle("Sigma"); - - fHDiffPixTh = new TH3D(); - fHDiffPixTh->SetNameTitle("fHDiffPixTh","3D: \\Theta, pixel no., Sigma^2-Sigmabar^2 (after padding)"); - MH::SetBinning(fHDiffPixTh, binstheta, binspixel, binsdiff); - fHDiffPixTh->SetXTitle("\\Theta [\\circ]"); - fHDiffPixTh->SetYTitle("pixel number"); - fHDiffPixTh->SetZTitle("Sigma^2-Sigmabar^2"); - - //-------------------------------------------------------------------- - - memset(fErrors, 0, sizeof(fErrors)); - - return kTRUE; -} - -// -------------------------------------------------------------------------- -// -// Do the Padding -// idealy the events to be padded should have been generated without NSB -// -Bool_t MPadSchweizer::Process() -{ - //*fLog << "Entry MPadSchweizer::Process();" << endl; - - Int_t rc=0; - - const UInt_t npix = fEvt->GetNumPixels(); - - fSigmabar->Calc(*fCam, *fPed, *fEvt); - //*fLog << "before padding : " << endl; - //fSigmabar->Print(""); - - - //$$$$$$$$$$$$$$$$$$$$$$$$$$ - // to simulate the situation that before the padding the NSB and - // electronic noise are zero : set Sigma = 0 for all pixels - //for (UInt_t i=0; ioperator[](i); - // Int_t j = pix.GetPixId(); - - // MPedestalPix &ppix = fPed->operator[](j); - // ppix.SetMeanRms(0.0); - //} - //$$$$$$$$$$$$$$$$$$$$$$$$$$ - - //------------------------------------------- - // Calculate average sigma of the event - // - Double_t sigbarold = fSigmabar->Calc(*fCam, *fPed, *fEvt); - //fSigmabar->Print(""); - - if (sigbarold > 0) - { - //*fLog << "MPadSchweizer::Process(); Sigmabar of event to be padded is > 0 : " - // << SigbarOld << ". Stop event loop " << endl; - // input data should have Sigmabar = 0; stop event loop - - rc = 1; - fErrors[rc]++; - return kCONTINUE; - } - - const Double_t theta = kRad2Deg*fMcEvt->GetTelescopeTheta(); - // *fLog << "Theta = " << Theta << endl; - - - //------------------------------------------- - // for the current Theta, - // generate a Sigmabar according to the histogram fHSigmaTheta - // - Double_t sigmabar=0; - Int_t binNumber = fHSigmaTheta->GetXaxis()->FindBin(theta); - - if ( binNumber < 1 || binNumber > fHSigmaTheta->GetNbinsX() ) - { - *fLog << "MPadSchweizer::Process(); binNumber out of range : Theta, binNumber = " - << theta << ", " << binNumber << "; Skip event " << endl; - // event cannot be padded; skip event - - rc = 2; - fErrors[rc]++; - return kCONTINUE; - } - else - { - TH1D* hsigma = - fHSigmaTheta->ProjectionY("", binNumber, binNumber, ""); - if ( hsigma->GetEntries() == 0 ) - { - *fLog << "MPadSchweizer::Process(); projection for Theta bin " - << binNumber << " has no entries; Skip event " << endl; - // event cannot be padded; skip event - delete hsigma; - - rc = 3; - fErrors[rc]++; - return kCONTINUE; - } - else - { - sigmabar = hsigma->GetRandom(); - - //*fLog << "Theta, bin number = " << Theta << ", " << binNumber - // << ", Sigmabar = " << Sigmabar << endl; - } - delete hsigma; - } - - - //------------------------------------------- - - //*fLog << "MPadSchweizer::Process(); SigbarOld, Sigmabar = " - // << SigbarOld << ", "<< Sigmabar << endl; - - // Skip event if target Sigmabar is <= SigbarOld - if (sigmabar <= sigbarold) - { - *fLog << "MPadSchweizer::Process(); target Sigmabar is less than SigbarOld : " - << sigmabar << ", " << sigbarold << ", Skip event" << endl; - - rc = 4; - fErrors[rc]++; - return kCONTINUE; - } - - - //------------------------------------------- - // - // Calculate average number of NSB photons to be added (lambdabar) - // from the value of Sigmabar, - // - making assumptions about the average electronic noise (elNoise2) and - // - using a fixed value (F2excess) for the excess noise factor - - Double_t F2excess = 1.3; - - const Double_t sigmabar2 = sigmabar*sigmabar; - - Int_t binTheta = fHDiffPixTheta->GetXaxis()->FindBin(theta); - if (binTheta != binNumber) - { - cout << "The binnings of the 2 histograms are not identical; aborting" - << endl; - return kERROR; - } - - // Get RMS of (Sigma^2-Sigmabar^2) in this Theta bin. - // The average electronic noise (to be added) has to be well above this RMS, - // otherwise the electronic noise of an individual pixel (elNoise2Pix) - // may become negative - TH1D* hnoise = fHDiffPixTheta->ProjectionZ("", binTheta, binTheta, - 0, 9999, ""); - - const Double_t elNoise2 = TMath::Min(hnoise->GetRMS(1), - sigmabar2 - sigbarold*sigbarold); // [photons] - - delete hnoise; - - //*fLog << "elNoise2 = " << elNoise2 << endl; - - const Double_t lambdabar = (sigmabar2 - sigbarold*sigbarold - elNoise2) - / F2excess; // [photons] - - // This value of lambdabar is the same for all pixels; - // note that lambdabar is normalized to the area of pixel 0 - - //---------- start loop over pixels --------------------------------- - // do the padding for each pixel - // - // pad only pixels - which are used (before image cleaning) - // - Double_t sig = 0; - Double_t sigma2 = 0; - Double_t diff = 0; - Double_t addSig2 = 0; - Double_t elNoise2Pix = 0; - - - for (UInt_t i=0; iGetPixRatio(j); - - MPedestalPix &ppix = (*fPed)[j]; - Double_t oldsigma = ppix.GetMeanRms(); - - //--------------------------------- - // throw the Sigma for this pixel - // - Int_t binPixel = fHDiffPixTheta->GetYaxis()->FindBin( (Double_t)j ); - - Int_t count; - Bool_t ok=kFALSE; - TH1D* hDiff; - TH1D* hSig; - - switch (fPadFlag) - { - case 1 : - // throw the Sigma for this pixel from the distribution fHDiffPixTheta - hDiff = fHDiffPixTheta->ProjectionZ("", binTheta, binTheta, - binPixel, binPixel, ""); - if ( hDiff->GetEntries() == 0 ) - { - *fLog << "MPadSchweizer::Process(); projection for Theta bin " - << binTheta << " and pixel bin " << binPixel - << " has no entries; aborting " << endl; - - rc = 5; - fErrors[rc]++; - return kCONTINUE; - } - - count = 0; - for (Int_t m=0; m<20; m++) - { - count += 1; - diff = hDiff->GetRandom(); - // the following condition ensures that elNoise2Pix > 0.0 - if ( (diff+sigmabar2-oldsigma*oldsigma/area- - lambdabar*F2excess) > 0 ) - { - ok = kTRUE; - break; - } - } - if (!ok) - { - *fLog << "Theta, j, count, Sigmabar, Diff = " << theta << ", " - << j << ", " << count << ", " << sigmabar << ", " - << diff << endl; - diff = lambdabar*F2excess + oldsigma*oldsigma/area - - sigmabar2; - } - delete hDiff; - sigma2 = diff + sigmabar2; - break; - - case 2 : - // throw the Sigma for this pixel from the distribution fHSigmaPixTheta - hSig = fHSigmaPixTheta->ProjectionZ("", binTheta, binTheta, - binPixel, binPixel, ""); - if ( hSig->GetEntries() == 0 ) - { - *fLog << "MPadSchweizer::Process(); projection for Theta bin " - << binTheta << " and pixel bin " << binPixel - << " has no entries; aborting " << endl; - - rc = 6; - fErrors[rc]++; - return kCONTINUE; - } - - count = 0; - for (Int_t m=0; m<20; m++) - { - count += 1; - sig = hSig->GetRandom(); - sigma2 = sig*sig/area; - // the following condition ensures that elNoise2Pix > 0.0 - if ( (sigma2-oldsigma*oldsigma/area-lambdabar*F2excess) > 0.0 ) - { - ok = kTRUE; - break; - } - } - if (!ok) - { - *fLog << "Theta, j, count, Sigmabar, Sig = " << theta << ", " - << j << ", " << count << ", " << sigmabar << ", " - << sig << endl; - sigma2 = lambdabar*F2excess + oldsigma*oldsigma/area; - } - delete hSig; - break; - } - - //--------------------------------- - // get the additional sigma^2 for this pixel (due to the padding) - addSig2 = sigma2*area - oldsigma*oldsigma; - - - //--------------------------------- - // get the additional electronic noise for this pixel - elNoise2Pix = addSig2 - lambdabar*F2excess*area; - - - //--------------------------------- - // throw actual number of additional NSB photons (NSB) - // and its RMS (sigmaNSB) - Double_t NSB0 = gRandom->Poisson(lambdabar*area); - Double_t arg = NSB0*(F2excess-1) + elNoise2Pix; - Double_t sigmaNSB0; - - if (arg >= 0) - { - sigmaNSB0 = sqrt( arg ); - } - else - { - *fLog << "MPadSchweizer::Process(); argument of sqrt < 0 : " - << arg << endl; - sigmaNSB0 = 0.0000001; - } - - - //--------------------------------- - // smear NSB0 according to sigmaNSB0 - // and subtract lambdabar because of AC coupling - Double_t NSB = gRandom->Gaus(NSB0, sigmaNSB0) - lambdabar*area; - - //--------------------------------- - - // add additional NSB to the number of photons - Double_t oldphotons = pix.GetNumPhotons(); - Double_t newphotons = oldphotons + NSB; - pix.SetNumPhotons( newphotons ); - - const Double_t areasqrt = sqrt(area); - - fHNSB->Fill( NSB/areasqrt ); - fHPhotons->Fill( oldphotons/areasqrt, newphotons/areasqrt ); - - - // error: add sigma of padded noise quadratically - Double_t olderror = pix.GetErrorPhot(); - Double_t newerror = sqrt( olderror*olderror + addSig2 ); - pix.SetErrorPhot( newerror ); - - - Double_t newsigma = sqrt( oldsigma*oldsigma + addSig2 ); - ppix.SetMeanRms( newsigma ); - - fHSigmaPedestal->Fill( oldsigma, newsigma ); - fHSigPixTh-> Fill( theta, (Double_t) j, newsigma ); - } - //---------- end of loop over pixels --------------------------------- - - // Calculate Sigmabar again and crosscheck - - Double_t sigbarnew = fSigmabar->Calc(*fCam, *fPed, *fEvt); - //*fLog << "after padding : " << endl; - //fSigmabar->Print(""); - - fHSigbarTheta->Fill( theta, sigbarnew ); - - // this loop is only for filling the histogram fHDiffPixTh - for (UInt_t i=0; iGetPixRatio(j); - - MPedestalPix &ppix = (*fPed)[j]; - Double_t newsigma = ppix.GetMeanRms(); - - fHDiffPixTh->Fill( theta, (Double_t) j, - newsigma*newsigma/area-sigbarnew*sigbarnew); - } - - - //*fLog << "Exit MPadSchweizer::Process();" << endl; - - rc = 0; - fErrors[rc]++; - - return kTRUE; - -} - -// -------------------------------------------------------------------------- -// -// -Bool_t MPadSchweizer::PostProcess() -{ - if (GetNumExecutions() != 0) - { - - *fLog << inf << endl; - *fLog << GetDescriptor() << " execution statistics:" << endl; - *fLog << dec << setfill(' '); - *fLog << " " << setw(7) << fErrors[1] << " (" << setw(3) - << (int)(fErrors[1]*100/GetNumExecutions()) - << "%) Evts skipped due to: Sigmabar_old > 0" << endl; - - *fLog << " " << setw(7) << fErrors[2] << " (" << setw(3) - << (int)(fErrors[2]*100/GetNumExecutions()) - << "%) Evts skipped due to: Zenith angle out of range" << endl; - - *fLog << " " << setw(7) << fErrors[3] << " (" << setw(3) - << (int)(fErrors[3]*100/GetNumExecutions()) - << "%) Evts skipped due to: No data for generating Sigmabar" << endl; - - *fLog << " " << setw(7) << fErrors[4] << " (" << setw(3) - << (int)(fErrors[4]*100/GetNumExecutions()) - << "%) Evts skipped due to: Target sigma <= Sigmabar_old" << endl; - - *fLog << " " << setw(7) << fErrors[5] << " (" << setw(3) - << (int)(fErrors[5]*100/GetNumExecutions()) - << "%) Evts skipped due to: No data for generating Sigma^2-Sigmabar^2" << endl; - - *fLog << " " << setw(7) << fErrors[6] << " (" << setw(3) - << (int)(fErrors[6]*100/GetNumExecutions()) - << "%) Evts skipped due to: No data for generating Sigma" << endl; - - *fLog << " " << fErrors[0] << " (" - << (int)(fErrors[0]*100/GetNumExecutions()) - << "%) Evts survived the padding!" << endl; - *fLog << endl; - - } - - //--------------------------------------------------------------- - TCanvas &c = *(MH::MakeDefCanvas("PadSchweizer", "", 900, 1200)); - c.Divide(3, 4); - gROOT->SetSelectedPad(NULL); - - - c.cd(1); - fHSigmaTheta->SetDirectory(NULL); - fHSigmaTheta->SetTitle("2D : Sigmabar, \\Theta (reference sample)"); - fHSigmaTheta->DrawCopy(); - fHSigmaTheta->SetBit(kCanDelete); - - //c.cd(1); - //fHSigmaPixTheta->DrawCopy(); - //fHSigmaPixTheta->SetBit(kCanDelete); - - c.cd(4); - fHSigbarTheta->SetDirectory(NULL); - fHSigbarTheta->DrawCopy(); - fHSigbarTheta->SetBit(kCanDelete); - - - c.cd(7); - fHNSB->SetDirectory(NULL); - fHNSB->DrawCopy(); - fHNSB->SetBit(kCanDelete); - - - //-------------------------------------------------------------------- - // draw the 3D histogram : Theta, pixel, Sigma^2-Sigmabar^2 - - TH2D *l; - - c.cd(2); - l = (TH2D*) ((TH3*)fHDiffPixTh)->Project3D("zx"); - l->SetDirectory(NULL); - l->SetTitle("Sigma^2-Sigmabar^2 vs. \\Theta (all pixels)"); - l->SetXTitle("\\Theta [\\circ]"); - l->SetYTitle("Sigma^2-Sigmabar^2"); - - l->DrawCopy("box"); - l->SetBit(kCanDelete);; - - c.cd(5); - l = (TH2D*) ((TH3*)fHDiffPixTh)->Project3D("zy"); - l->SetDirectory(NULL); - l->SetTitle("Sigma^2-Sigmabar^2 vs. pixel number (all \\Theta)"); - l->SetXTitle("pixel"); - l->SetYTitle("Sigma^2-Sigmabar^2"); - - l->DrawCopy("box"); - l->SetBit(kCanDelete);; - - c.cd(8); - ((TH2*)fHDiffPixTh)->DrawCopy(); - - - //-------------------------------------------------------------------- - // draw the 3D histogram : Theta, pixel, Sigma - - TH2D *k; - - c.cd(3); - k = (TH2D*) ((TH3*)fHSigPixTh)->Project3D("zx"); - k->SetDirectory(NULL); - k->SetTitle("Sigma vs. \\Theta (all pixels)"); - k->SetXTitle("\\Theta [\\circ]"); - k->SetYTitle("Sigma"); - - k->DrawCopy("box"); - k->SetBit(kCanDelete); - - c.cd(6); - k = (TH2D*) ((TH3*)fHSigPixTh)->Project3D("zy"); - k->SetDirectory(NULL); - k->SetTitle("Sigma vs. pixel number (all \\Theta)"); - k->SetXTitle("pixel"); - k->SetYTitle("Sigma"); - - k->DrawCopy("box"); - k->SetBit(kCanDelete);; - - c.cd(9); - ((TH2*)fHSigPixTh)->DrawCopy(); - - //-------------------------------------------------------------------- - - c.cd(10); - fHSigmaPedestal->SetDirectory(NULL); - fHSigmaPedestal->DrawCopy(); - fHSigmaPedestal->SetBit(kCanDelete); - - c.cd(11); - fHPhotons->SetDirectory(NULL); - fHPhotons->DrawCopy(); - fHPhotons->SetBit(kCanDelete); - - //-------------------------------------------------------------------- - - return kTRUE; -} - -// -------------------------------------------------------------------------- - Index: trunk/MagicSoft/Mars/manalysis/MPadSchweizer.h =================================================================== --- trunk/MagicSoft/Mars/manalysis/MPadSchweizer.h (revision 2019) +++ (revision ) @@ -1,75 +1,0 @@ -#ifndef MARS_MPadSchweizer -#define MARS_MPadSchweizer - -#ifndef MARS_MTask -#include "MTask.h" -#endif - -#ifndef MARS_MH -#include "MH.h" -#endif - -class TH1D; -class TH2D; -class TH3D; - -class MGeomCam; -class MCerPhotEvt; -class MPedestalCam; -class MMcEvt; -class MSigmabar; -class MParList; - -class MPadSchweizer : public MTask -{ -private: - MGeomCam *fCam; - MCerPhotEvt *fEvt; - MSigmabar *fSigmabar; - MMcEvt *fMcEvt; - MPedestalCam *fPed; - - Int_t fPadFlag; - Int_t fRunType; - Int_t fGroup; - - Int_t fErrors[7]; - - // plots used for the padding - TH2D *fHSigmaTheta; // 2D-histogram (sigmabar vs. Theta) - TH3D *fHSigmaPixTheta; // 3D-histogram (Theta, pixel, sigma) - TH3D *fHDiffPixTheta; // 3D-histogram (Theta, pixel, sigma^2-sigmabar^2) - - // plots for checking the padding - TH2D *fHSigmaPedestal; // 2D-histogram : pedestal sigma after - // versus before padding - TH2D *fHPhotons; // 2D-histogram : no.of photons after - // versus before padding - TH1D *fHNSB; // 1D-histogram : additional NSB - - -public: - MPadSchweizer(const char *name=NULL, const char *title=NULL); - ~MPadSchweizer(); - - void SetHistograms(TH2D *hist2, TH3D *hist3, TH3D *hist3Diff); - - Bool_t PreProcess(MParList *pList); - Bool_t Process(); - Bool_t PostProcess(); - - void SetPadFlag(Int_t padflag); - - ClassDef(MPadSchweizer, 0) // task for the padding (Schweizer) -}; - -#endif - - - - - - - - -