Changeset 2749

Timestamp:

12/23/03 16:26:30 (21 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft

Files:

• Mars/Changelog (modified) (1 diff)
• Mars/manalysis/MExtractedSignalPix.cc (modified) (1 diff)
• Mars/manalysis/MGeomApply.cc (modified) (2 diffs)
• Mars/manalysis/MMcCalibrationCalc.cc (modified) (7 diffs)
• Mars/manalysis/MMcCalibrationCalc.h (modified) (2 diffs)
• include-Classes/MMcFormat/MMcFadcHeader.hxx (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -4 +4 @@
 
                                                  -*-*- END OF LINE -*-*-
+
+ 2003/12/23: Abelardo Moralejo
+
+   * manalysis/MExtractedSignalPix.cc
+     - in Clear(), initialized fExtractedSignalHiGain to 0 instead 
+       of -1. In MC files run without noise, the ADC info of pixels 
+       with no C-photons is missing. Later noise will be added to all 
+       pixels in the padding, and "empty" ones must have a 0 before 
+       the addition of this noise.
+
+   * manalysis/MGeomApply.cc
+     - initialized size of MPedPhotCam.
+
+   * manalysis/MMcCalibrationCalc.[cc,h]
+     - Removed setter function SetConversionHiLo: the ratio from high 
+       to low gain is now read from MMcFadcHeader, and set in ReInit.
+     - Fixed bug: fADC2PhOuter was initialized before the MMcFadcHeader
+       was available. Now it is done in ReInit.
+     - fPedPhotCam is now initialized in PreProcess and not in ReInit.
+     - fGeom is now set in ReInit (there should be one MGeomCam per
+       run header in the camera file).
+
+   * mmc/MMcFadcHeader.hxx
+     - added getter function GetLow2HighGain()
+
 
  2003/12/23: Thomas Bretz

trunk/MagicSoft/Mars/manalysis/MExtractedSignalPix.cc

@@ -64 +64 @@
 {
 
-  fExtractedSignalHiGain = -1.;
+  fExtractedSignalHiGain = 0.;
+  //  fExtractedSignalHiGain = -1.;
   fExtractedSignalHiGainError = -1.;
+
   fIsLoGainUsed = kFALSE;
   fNumHiGainSaturated = 0;

trunk/MagicSoft/Mars/manalysis/MGeomApply.cc

@@ -58 +58 @@
 #include "MPedestalCam.h"
 #include "MCalibrationCam.h"
+#include "MPedPhotCam.h"
 #include "MBlindPixels.h"
 
@@ -122 +123 @@
 
 
+    MPedPhotCam *pedphot = (MPedPhotCam*)pList->FindObject(AddSerialNumber("MPedPhotCam"));
+    if (pedphot)
+        pedphot->InitSize(cam->GetNumPixels());
+
+
     MBlindPixels *bnd = (MBlindPixels*)pList->FindObject(AddSerialNumber("MBlindPixels"));
     if (bnd)

trunk/MagicSoft/Mars/manalysis/MMcCalibrationCalc.cc

@@ -101 +101 @@
 {
   fADC2PhInner = x;
-  fADC2PhOuter = x * 
-    fHeaderFadc->GetAmplitud() / fHeaderFadc->GetAmplitudOuter();
 
   return;
@@ -125 +123 @@
     }
 
-  if ( ! pList->FindCreateObj(AddSerialNumber("MPedPhotCam")))
+  fPedPhotCam = (MPedPhotCam*) pList->FindCreateObj(AddSerialNumber("MPedPhotCam"));
+  if ( ! fPedPhotCam)
     {
       *fLog << err << dbginf << "Cannot create MPedPhotCam... aborting." << endl;
@@ -131 +130 @@
     }
 
-  fGeom = (MGeomCam*) pList->FindObject(AddSerialNumber("MGeomCam"));
-  if ( ! fGeom )
-    {
-      *fLog << err << dbginf << "Cannot find MGeomCam... aborting." << endl;
-      return kFALSE;
-    }
-
   fSignalCam = (MExtractedSignalCam*) pList->FindObject(AddSerialNumber("MExtractedSignalCam"));
   if ( ! fSignalCam)
@@ -145 +137 @@
     }
 
+    fCalCam = (MCalibrationCam*)pList->FindObject(AddSerialNumber("MCalibrationCam"));
+
+    if (!fCalCam)
+    {
+        *fLog << err << dbginf << "Could not create MCalibrationCam... aborting. " << endl;
+        return kFALSE;
+    }
+
   return kTRUE;
 
@@ -163 +163 @@
 
     //
-    // Now check the existence of all necessary containers. This has
-    // to be done only if this is a MC file.
-    //
+    // Now check the existence of all necessary containers.
+    //
+
+    fGeom = (MGeomCam*) pList->FindObject(AddSerialNumber("MGeomCam"));
+    if ( ! fGeom )
+      {
+        *fLog << err << dbginf << "Cannot find MGeomCam... aborting." << endl;
+        return kFALSE;
+      }
 
     fHeaderFadc = (MMcFadcHeader*)pList->FindObject(AddSerialNumber("MMcFadcHeader"));
@@ -174 +180 @@
     }
 
-    fCalCam = (MCalibrationCam*)pList->FindObject(AddSerialNumber("MCalibrationCam"));
-
-    if (!fCalCam)
-    {
-        *fLog << err << dbginf << "Could not create MCalibrationCam... aborting. " << endl;
-        return kFALSE;
-    }
-
-    fPedPhotCam = (MPedPhotCam*)pList->FindObject(AddSerialNumber("MPedPhotCam"));
-
-    if (!fPedPhotCam)
-    {
-        *fLog << err << dbginf << "Could not create MPedPhotCam... aborting. " << endl;
-        return kFALSE;
-    }
-
-    fPedPhotCam->InitSize(fGeom->GetNumPixels());
+    //
+    // Set the ADC to photons conversion factor for outer pixels:
+    //
+    fADC2PhOuter = fADC2PhInner *
+      (fHeaderFadc->GetAmplitud() / fHeaderFadc->GetAmplitudOuter());
+
+    //
+    // Set the conversion factor between high and low gain:
+    //
+    fConversionHiLo = fHeaderFadc->GetLow2HighGain();
 
     return kTRUE;
@@ -231 +230 @@
         //
 
-        Float_t adc2phot = (fGeom->GetPixRatio(i) > fGeom->GetPixRatio(0))?
+        Float_t adc2phot = (fGeom->GetPixRatio(i) < fGeom->GetPixRatio(0))?
           fADC2PhOuter : fADC2PhInner;
 

trunk/MagicSoft/Mars/manalysis/MMcCalibrationCalc.h

@@ -28 +28 @@
     Float_t fADC2PhInner; // Conversion factor from ADC counts to photons
     Float_t fADC2PhOuter; // for inner and outer pixels.
-
     Float_t fConversionHiLo;
 
@@ -35 +34 @@
 
     void SetADC2PhInner(Float_t x);
-    void SetConversionHiLo(Float_t x) { fConversionHiLo = x; }
 
     ClassDef(MMcCalibrationCalc, 0)   // Task which obtains, for MC files, the pedestal mean and rms, and the calibration factor from ADC counts to photons. 

trunk/MagicSoft/include-Classes/MMcFormat/MMcFadcHeader.hxx

@@ -83 +83 @@
   Float_t GetAmplitud() const { return fAmplFadc; }
   Float_t GetAmplitudOuter() const { return fAmplFadcOuter; }
+  Float_t GetLow2HighGain() const { return fLow2HighGain; }
 
   UInt_t GetNumPixel() const { return MFADC_CHANNELS; }