Changeset 3651

Timestamp:

04/05/04 16:48:56 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• macros/calibration.C (modified) (1 diff)
• manalysis/MGeomApply.cc (modified) (2 diffs)
• manalysis/MHPedestalCam.cc (modified) (6 diffs)
• manalysis/MHPedestalPix.h (modified) (1 diff)
• mcalib/CalibLinkDef.h (modified) (1 diff)
• mcalib/MCalibrationChargeCam.h (modified) (1 diff)
• mcalib/MCalibrationChargePix.h (modified) (3 diffs)
• mcalib/MCalibrationPedCam.h (added)
• mcalib/Makefile (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -24 +24 @@
    * mcalib/Makefile
    * mcalib/CalibLinkDef.h
+   * manalysis/MGeomApply.cc
      - new class to store information obtained from MHPedestalCam
 

trunk/MagicSoft/Mars/macros/calibration.C

@@ -59 +59 @@
   if (!pedloop.Process())
     return;
+
+  //
+  // Create a empty Parameter List and an empty Task List
+  // The tasklist is identified in the eventloop by its name
+  //
+  MParList  plist0;
+  MTaskList tlist0;
+  plist0.AddToList(&tlist0);
+  
+  //
+  // Now setup the tasks and tasklist for the pedestals:
+  // ---------------------------------------------------
+  //
+  MReadMarsFile read("Events");
+  read.DisableAutoScheme();
+  static_cast<MRead&>(read).AddFiles(pruns);
+
+  MGeomApply      geomapl;
+  MExtractSignal  sigcalc;
+  
+  //
+  // Additionally to calculating the pedestals, 
+  // you can fill histograms and look at them
+  //
+  MFillH fill("MHPedestalCam", "MExtractedSignalCam");
+  
+  tlist0.AddToList(&read);
+  tlist0.AddToList(&geomapl);
+  tlist0.AddToList(&sigcalc);
+  tlist0.AddToList(&fill);
+  
+  MGeomCamMagic  geomcam;
+  MHPedestalCam  hpedcam;
+  plist0.AddToList(&geomcam);
+  plist0.AddToList(&pedloop.GetPedestalCam());
+  plist0.AddToList(&badcam);
+  plist0.AddToList(&hpedcam);
+  
+  //
+  // Create and setup the eventloop
+  //
+  MEvtLoop evtloop0;
+  
+  evtloop0.SetParList(&plist0);
+  evtloop0.SetDisplay(display);
+  
+  //
+  // Execute first analysis
+  //
+  if (!evtloop0.Eventloop())
+    return;
+
+  tlist0.PrintStatistics();
+
+  for (Int_t aidx=0;aidx<2;aidx++)
+    {
+      hpedcam.GetAverageHiGainArea(aidx).DrawClone("fourierevents");
+      hpedcam.GetAverageLoGainArea(aidx).DrawClone("fourierevents");
+    }
+  
+  for (Int_t sector=1;sector<7;sector++)
+    {
+      hpedcam.GetAverageHiGainSector(sector).DrawClone("fourierevents");
+      hpedcam.GetAverageLoGainSector(sector).DrawClone("fourierevents");
+    }
+  
 
   // 

trunk/MagicSoft/Mars/manalysis/MGeomApply.cc

@@ -74 +74 @@
 #include "MCalibrationRelTimeCam.h"
 #include "MCalibrationQECam.h"
+#include "MCalibrationPedCam.h"
 #include "MPedPhotCam.h"
 #include "MExtractedSignalCam.h"
@@ -156 +157 @@
       (MCalibrationQECam*)pList->FindObject(AddSerialNumber("MCalibrationQECam"));
     if (qe)
-        qe->InitSize(cam->GetNumPixels());
-
+      {
+        qe->InitSize           ( cam->GetNumPixels()  );
+        qe->InitAverageAreas   ( cam->GetNumAreas()   );
+        qe->InitAverageSectors ( cam->GetNumSectors() );
+      }
+    
+    MCalibrationPedCam *pcam = 
+      (MCalibrationPedCam*)pList->FindObject(AddSerialNumber("MCalibrationPedCam"));
+    if (pcam)
+      {
+        pcam->InitSize           ( cam->GetNumPixels()  );
+        pcam->InitAverageAreas   ( cam->GetNumAreas()   );
+        pcam->InitAverageSectors ( cam->GetNumSectors() );
+      }
+    
     MPedPhotCam *pedphot = (MPedPhotCam*)pList->FindObject(AddSerialNumber("MPedPhotCam"));
     if (pedphot)

trunk/MagicSoft/Mars/manalysis/MHPedestalCam.cc

@@ -90 +90 @@
 /////////////////////////////////////////////////////////////////////////////
 #include "MHPedestalCam.h"
+#include "MHPedestalPix.h"
 
 #include "MLog.h"
@@ -96 +97 @@
 #include "MParList.h"
 
-#include "MHPedestalPix.h"
-
 #include "MExtractedSignalCam.h"
 #include "MExtractedSignalPix.h"
@@ -106 +105 @@
 #include "MGeomCam.h"
 #include "MGeomPix.h"
+
+#include "MCalibrationPedCam.h"
 
 ClassImp(MHPedestalCam);
@@ -133 +134 @@
 // - MPedestalCam
 // - MExtractedSignalCam
+//
+// Searches or creates:
+// - MCalibrationPedCam
 //
 // Retrieves from MExtractedSignalCam:
@@ -179 +183 @@
     }
 
-  Float_t sliceshi = signal->GetNumUsedHiGainFADCSlices();
-  Float_t sliceslo = signal->GetNumUsedLoGainFADCSlices();
-
-  if (sliceshi == 0.)
-    {
-      gLog << err << "Number of used signal slices in MExtractedSignalCam is zero  ... abort." 
-           << endl;
-      return kFALSE;
-    }
-
-  if (fExtractHiGainSlices != 0. && sliceshi != fExtractHiGainSlices )
-    {
-      gLog << err << "Number of used High Gain signal slices changed in MExtractedSignalCam  ... abort." 
-           << endl;
-      return kFALSE;
-    }
-
-  if (fExtractLoGainSlices != 0. && sliceslo != fExtractLoGainSlices )
-    {
-      gLog << err << "Number of used Low Gain signal slices changed in MExtractedSignalCam  ... abort." 
-           << endl;
-      return kFALSE;
-    }
-
-  fExtractHiGainSlices = sliceshi;
-  fExtractLoGainSlices = sliceslo;
-
   const Int_t npixels  = fGeom->GetNumPixels();
   const Int_t nsectors = fGeom->GetNumSectors();
   const Int_t nareas   = fGeom->GetNumAreas();
+
+  fCam = (MCalibrationCam*)pList->FindObject("MCalibrationPedCam");
+  if (!fCam)
+    {
+      fCam = (MCalibrationCam*)pList->FindCreateObj(AddSerialNumber("MCalibrationPedCam"));
+      if (!fCam)
+        {
+          gLog << err << "Cannot find nor create MCalibrationPedCam ... abort." << endl;
+          return kFALSE;
+        }
+      else 
+        {
+          fCam->InitSize(npixels);
+          fCam->InitAverageAreas(nareas);
+          fCam->InitAverageSectors(nsectors);
+        }
+    }
+  
+  Float_t sliceshi = signal->GetNumUsedHiGainFADCSlices();
+  Float_t sliceslo = signal->GetNumUsedLoGainFADCSlices();
+
+  if (sliceshi == 0.)
+    {
+      gLog << err << "Number of used signal slices in MExtractedSignalCam is zero  ... abort." 
+           << endl;
+      return kFALSE;
+    }
+
+  if (fExtractHiGainSlices != 0. && sliceshi != fExtractHiGainSlices )
+    {
+      gLog << err << "Number of used High Gain signal slices changed in MExtractedSignalCam  ... abort." 
+           << endl;
+      return kFALSE;
+    }
+
+  if (fExtractLoGainSlices != 0. && sliceslo != fExtractLoGainSlices )
+    {
+      gLog << err << "Number of used Low Gain signal slices changed in MExtractedSignalCam  ... abort." 
+           << endl;
+      return kFALSE;
+    }
+
+  fExtractHiGainSlices = sliceshi;
+  fExtractLoGainSlices = sliceslo;
 
   if (fHiGainArray->GetEntries()==0)
@@ -431 +452 @@
 {
 
-  FitHiGainArrays((MCalibrationCam&)(*fCam),*fBadPixels, 
+  FitHiGainArrays((*fCam),*fBadPixels, 
                   MBadPixelsPix::kHiGainNotFitted,
                   MBadPixelsPix::kHiGainOscillating);
-  FitLoGainArrays((MCalibrationCam&)(*fCam),*fBadPixels,
+  FitLoGainArrays((*fCam),*fBadPixels,
                   MBadPixelsPix::kLoGainNotFitted,
                   MBadPixelsPix::kLoGainOscillating);

trunk/MagicSoft/Mars/manalysis/MHPedestalPix.h

@@ -11 +11 @@
 private:
 
-  static const Int_t   fgChargeNbins;        // Default for fNbins          (now set to: 450  )
-  static const Axis_t  fgChargeFirst;        // Default for fFirst          (now set to: -0.5  )
-  static const Axis_t  fgChargeLast;         // Default for fLast           (now set to: 449.5)
+  static const Int_t   fgChargeNbins;        // Default for fNbins   (now set to: 450  )
+  static const Axis_t  fgChargeFirst;        // Default for fFirst   (now set to: -0.5 )
+  static const Axis_t  fgChargeLast;         // Default for fLast    (now set to: 449.5)
 
   Float_t fNSlices;                         // Number of FADC slices summed in extraction

trunk/MagicSoft/Mars/mcalib/CalibLinkDef.h

@@ -12 +12 @@
 #pragma link C++ class MCalibrationRelTimeCam+;
 #pragma link C++ class MCalibrationRelTimePix+;
+#pragma link C++ class MCalibrationPedCam+;
 #pragma link C++ class MCalibrationQECam+;
 #pragma link C++ class MCalibrationQEPix+;

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCam.h

@@ -12 +12 @@
 class MCalibrationChargePINDiode;
 class MCalibrationChargePix;
-class MGeomCam;
-class MBadPixelsCam;
 class MCalibrationChargeCam : public MCalibrationCam
 {

trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.h

@@ -28 +28 @@
   Float_t fConversionHiLo;                  // Conversion factor betw. Hi Gain and Lo Gain  
   Float_t fConversionHiLoVar;               // Variance Conversion factor betw. Hi and Lo Gain
-  Float_t fLoGainPedRms;                    // Pedestal RMS of Low Gain
-  Float_t fLoGainPedRmsVar;                 // Pedestal RMS Variance of Low Gain
+  Float_t fLoGainPedRmsSquare;              // Pedestal RMS square of Low Gain
+  Float_t fLoGainPedRmsSquareVar;           // Pedestal RMS square Variance of Low Gain
   Float_t fMeanConversionFFactorMethod;     // Conversion factor (F-factor method)
   Float_t fMeanConversionBlindPixelMethod;  // Conversion factor (Blind Pixel method)
@@ -40 +40 @@
   Float_t fPheFFactorMethodVar;             // Variance number of Phe's (F-factor method)
   Float_t fPheFFactorMethodLimit;           // Min. number Photo-electrons for pix to be accepted.
-  Float_t fRSigma;                          // Reduced sigma
-  Float_t fRSigmaVar;                       // Variance Reduced sigma
+  Float_t fRSigmaSquare;                    // Square of Reduced sigma
+  Float_t fRSigmaSquareVar;                 // Variance Reduced sigma
   Float_t fSigmaConversionFFactorMethod;    // Sigma Conversion factor (F-factor method)
   Float_t fSigmaConversionBlindPixelMethod; // Sigma Conversion factor (Blind Pixel method)
@@ -57 +57 @@
   enum  { kBlindPixelMethodValid, kFFactorMethodValid,
           kPINDiodeMethodValid, kCombinedMethodValid }; // Possible bits to be set
+
+  const Float_t GetConversionHiLoRelVar()  const;
+  const Float_t GetChargeRelVar()  const;
 
 public:

trunk/MagicSoft/Mars/mcalib/Makefile

@@ -45 +45 @@
            MCalibrationQECam.cc \
            MCalibrationQEPix.cc  \
+           MCalibrationPedCam.cc \
            MCalibrationChargeCalc.cc \
            MCalibrationChargeCam.cc \