Changeset 3551 for trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeCam.cc

Timestamp:

03/19/04 16:56:19 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeCam.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeCam.cc

@@ -65 +65 @@
 #include <TCanvas.h>
 #include <TPad.h>
+#include <TText.h>
+#include <TPaveText.h>
 
 #include "MLog.h"
@@ -151 +153 @@
     fNumOuterPixels = 0; 
     fNumExcluded    = 0;
+
+    fAverageInnerSat            = 0;
+    fAverageOuterSat            = 0;  
+    fAverageInnerPixSigma       = 0.;
+    fAverageOuterPixSigma       = 0.;  
+    fAverageInnerPixSigmaErr    = 0.;
+    fAverageOuterPixSigmaErr    = 0.;  
+    fAverageInnerPixRelSigma    = 0.;
+    fAverageOuterPixRelSigma    = 0.;  
+    fAverageInnerPixRelSigmaErr = 0.;
+    fAverageOuterPixRelSigmaErr = 0.;  
+
 }
 
@@ -169 +183 @@
   delete fAverageHiGainOuterPix;
   delete fAverageLoGainOuterPix;
+
 }
 
@@ -542 +557 @@
     FinalizeLoGainHists(*fAverageLoGainOuterPix,*fCam->GetAverageOuterPix(),*fCam->GetAverageOuterBadPix());
 
-    fCam->GetAverageInnerPix()->SetSigmaCharge   (fCam->GetAverageInnerPix()->GetSigmaCharge()   
-                                                  *TMath::Sqrt((Float_t)fNumInnerPixels));
-    fCam->GetAverageOuterPix()->SetSigmaCharge   (fCam->GetAverageOuterPix()->GetSigmaCharge()   
-                                                  *TMath::Sqrt((Float_t)fNumOuterPixels));
-    fCam->GetAverageInnerPix()->SetSigmaChargeErr(fCam->GetAverageInnerPix()->GetSigmaChargeErr()
-                                                  *TMath::Sqrt((Float_t)fNumInnerPixels));
-    fCam->GetAverageOuterPix()->SetSigmaChargeErr(fCam->GetAverageOuterPix()->GetSigmaChargeErr()
-                                                  *TMath::Sqrt((Float_t)fNumOuterPixels));
-    
+    FinalizeAveragePix(*fCam->GetAverageInnerPix(),fNumInnerPixels, 
+                       fAverageInnerPixSigma, fAverageInnerPixSigmaErr,
+                       fAverageInnerPixRelSigma, fAverageInnerPixRelSigmaErr,
+                       fAverageInnerSat);
+    FinalizeAveragePix(*fCam->GetAverageOuterPix(),fNumOuterPixels, 
+                       fAverageOuterPixSigma, fAverageOuterPixSigmaErr,
+                       fAverageOuterPixRelSigma, fAverageOuterPixRelSigmaErr,
+                       fAverageOuterSat);
+
     return kTRUE;
 }
@@ -693 +708 @@
       }
 }    
+
+void MHCalibrationChargeCam::FinalizeAveragePix(MCalibrationChargePix &pix, Int_t npix, 
+                                                Float_t &sigma, Float_t &sigmaerr, 
+                                                Float_t &relsigma, Float_t &relsigmaerr, 
+                                                Bool_t &b)
+{
+
+  if (pix.IsHiGainSaturation())
+    b = kTRUE;
+
+  sigma    = pix.GetSigmaCharge    () * TMath::Sqrt((Float_t)npix);
+  sigmaerr = pix.GetSigmaChargeErr () * TMath::Sqrt((Float_t)npix);
+
+  relsigma     = sigma / pix.GetMeanCharge();
+
+  relsigmaerr  = sigmaerr*sigmaerr / sigma / sigma;
+  relsigmaerr += pix.GetMeanChargeErr()*pix.GetMeanChargeErr() / pix.GetMeanCharge() / pix.GetMeanCharge();
+
+  relsigmaerr *= relsigma;
+  relsigmaerr  = TMath::Sqrt(relsigmaerr);
+
+  pix.SetSigmaCharge   (sigma);
+  pix.SetSigmaChargeErr(sigmaerr);
+
+}
+
 
 Bool_t MHCalibrationChargeCam::GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type) const
@@ -729 +770 @@
   pad->cd(1);
 
-  if (!fAverageHiGainInnerPix->IsEmpty())
-    gPad->SetLogy();
   fAverageHiGainInnerPix->Draw(opt);
+  
+  if (!fAverageInnerSat)
+    DrawAverageSigma(fAverageInnerSat, 1,
+                     fAverageInnerPixSigma, fAverageInnerPixSigmaErr,
+                     fAverageInnerPixRelSigma, fAverageInnerPixRelSigmaErr);
 
   pad->cd(2);
 
-  if (!fAverageLoGainInnerPix->IsEmpty())
-    gPad->SetLogy();
   fAverageLoGainInnerPix->Draw(opt);
 
+  if (fAverageInnerSat)
+    DrawAverageSigma(fAverageInnerSat, 1,
+                     fAverageInnerPixSigma, fAverageInnerPixSigmaErr,
+                     fAverageInnerPixRelSigma, fAverageInnerPixRelSigmaErr);
+
   pad->cd(3);
 
-  if (!fAverageHiGainOuterPix->IsEmpty())
-    gPad->SetLogy();
   fAverageHiGainOuterPix->Draw(opt);
 
+  if (!fAverageOuterSat)
+    DrawAverageSigma(fAverageOuterSat, 0,
+                     fAverageOuterPixSigma, fAverageOuterPixSigmaErr,
+                     fAverageOuterPixRelSigma, fAverageOuterPixRelSigmaErr);
+
   pad->cd(4);
 
-  if (!fAverageLoGainOuterPix->IsEmpty())
-    gPad->SetLogy();
   fAverageLoGainOuterPix->Draw(opt);
 
-}
+  if (fAverageOuterSat)
+    DrawAverageSigma(fAverageOuterSat, 0,
+                     fAverageOuterPixSigma, fAverageOuterPixSigmaErr,
+                     fAverageOuterPixRelSigma, fAverageOuterPixRelSigmaErr);
+}
+
+void MHCalibrationChargeCam::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 ? "Inner" : "Outer",
+                                 " Pixels ", sat ? "Low Gain" : "High Gain");
+      TText *txt1 = text->AddText(line1.Data());
+      const TString line2 = Form("Sigma per Pixel: %2.2f #pm %2.2f",sigma,sigmaerr);
+      TText *txt2 = text->AddText(line2.Data());
+      const TString line3 = Form("Rel. Sigma per Pixel: %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);
+    }
+}
+