Changeset 4152 for trunk

Timestamp:

05/24/04 18:13:47 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mcalib/MHCalibrationTestCam.cc (modified) (5 diffs)
• mcalib/MHCalibrationTestCam.h (modified) (3 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -26 +26 @@
         
  2004/05/24: Markus Gaug
+
+   * mcalib/MHCalibrationTestCam.[h,cc]
+     - number of photons/area are now stored in TArray's in order to 
+       retrieve them easier in later processing from stored files.
 
    * mcalib/MCalibrationCam.cc

trunk/MagicSoft/Mars/mcalib/MHCalibrationTestCam.cc

@@ -263 +263 @@
       }
   }
-    
-
+
+
+  fMeanMeanPhotPerArea.Set(nareas);   
+  fRmsMeanPhotPerArea .Set(nareas); 
+  fMeanSigmaPhotPerArea.Set(nareas);  
+  fRmsSigmaPhotPerArea.Set(nareas);
+ 
   return kTRUE;
 }
@@ -359 +364 @@
 {
 
+  TArrayI numaidx;
+  numaidx.Set(fGeom->GetNumAreas());
+
   for (Int_t i=0; i<fHiGainArray->GetSize(); i++)
     {
@@ -378 +386 @@
       hist.CreateFourierSpectrum();
       
-    }
-  
+      const Float_t area = (*fGeom)[i].GetA();
+      const Int_t   aidx = (*fGeom)[i].GetAidx();
+
+      fMeanMeanPhotPerArea[aidx]  += hist.GetMean() / area;
+      fRmsMeanPhotPerArea [aidx]  += hist.GetMean() / area * hist.GetMean()  / area;
+      fMeanSigmaPhotPerArea[aidx] += hist.GetSigma()/ area;
+      fRmsSigmaPhotPerArea [aidx] += hist.GetSigma()/ area * hist.GetSigma() / area;
+      numaidx[aidx]++;
+    }
+
+
   for (Int_t j=0; j<fAverageHiGainAreas->GetSize(); j++)
     {
@@ -395 +412 @@
       hist.CreateFourierSpectrum();
       
-
+      fRmsMeanPhotPerArea [j]  -= fMeanMeanPhotPerArea [j]*fMeanMeanPhotPerArea [j]/numaidx[j];
+      fRmsSigmaPhotPerArea[j]  -= fMeanSigmaPhotPerArea[j]*fMeanSigmaPhotPerArea[j]/numaidx[j];
+
+      fMeanMeanPhotPerArea [j]  /=  numaidx[j];
+      fMeanSigmaPhotPerArea[j]  /=  numaidx[j];
+      fRmsMeanPhotPerArea  [j]  /=  numaidx[j]-1.;
+      fRmsSigmaPhotPerArea [j]  /=  numaidx[j]-1.; 
   }
 
@@ -455 +478 @@
 // =================
 //
-// 7:  Fitted Mean Test Calibration (MHGausEvents::GetMean()) by MGeomPix::GetA()
-// 8:  Fitted Mean Error Calibration (MHGausEvents::GetMeanErr()) by MGeomPix::GetA()
-// 9:  Fitted Sigma Test Calibration (MHGausEvents::GetSigma()) by MGeomPix::GetA()
+// 7:  Fitted Mean Test Calibration   (MHGausEvents::GetMean())     by MGeomPix::GetA()
+// 8:  Fitted Mean Error Calibration  (MHGausEvents::GetMeanErr())  by MGeomPix::GetA()
+// 9:  Fitted Sigma Test Calibration  (MHGausEvents::GetSigma())    by MGeomPix::GetA()
 // 10: Fitted Sigma Error Calibration (MHGausEvents::GetSigmaErr()) by MGeomPix::GetA()
 //

trunk/MagicSoft/Mars/mcalib/MHCalibrationTestCam.h

@@ -4 +4 @@
 #ifndef MARS_MHCalibrationCam
 #include "MHCalibrationCam.h"
+#endif
+#ifndef ROOT_TArrayF
+#include "TArrayF.h"
 #endif
 
@@ -12 +15 @@
 private:
 
+  TArrayF fMeanMeanPhotPerArea;
+  TArrayF fRmsMeanPhotPerArea   ;
+  TArrayF fMeanSigmaPhotPerArea;
+  TArrayF fRmsSigmaPhotPerArea   ;
+  
   Bool_t ReInitHists(MParList *pList);
   Bool_t FillHists(const MParContainer *par, const Stat_t w=1);
@@ -25 +33 @@
   void DrawPixelContent(Int_t idx) const;
 
+  const Float_t GetMeanMeanPhotPerArea       ( const Int_t aidx ) const   { return fMeanMeanPhotPerArea [aidx]; }
+  const Float_t GetRmsMeanPhotPerArea        ( const Int_t aidx ) const   { return fRmsMeanPhotPerArea  [aidx]; }
+  const Float_t GetMeanSigmaPhotPerArea      ( const Int_t aidx ) const   { return fMeanSigmaPhotPerArea [aidx]; }
+  const Float_t GetRmsSigmaPhotPerArea       ( const Int_t aidx ) const   { return fRmsSigmaPhotPerArea  [aidx]; }
+
   void CalcAverageSigma();