Changeset 7832

Timestamp:

08/02/06 08:03:32 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• NEWS (modified) (1 diff)
• mhist/MHCamEvent.cc (modified) (1 diff)
• msignal/MExtractTimeAndCharge.cc (modified) (1 diff)
• msignal/MExtractTimeAndChargeDigitalFilter.cc (modified) (3 diffs)
• msignal/MSignalCam.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -18 +18 @@
 
                                                  -*-*- END OF LINE -*-*-
+
+ 2006/08/02 Thomas Bretz
+
+   * mhist/MHCamEvent.cc:
+     - removed some obsolete includes
+
+   * msignal/MExtractTimeAndCharge.cc:
+     - reformatted position of a parenthesis
+
+   * msignal/MExtractTimeAndChargeDigitalFilter.cc:
+     - changed default lo-gain offset to new value 1.05 which was
+       determined plotting arrival time vs. signal height.
+     - added fHiGainFirst to time (as it is in the spline extractor)
+
+   * msignal/MSignalCam.cc:
+     - after extesive tests using the arrival time vs signal height
+       histograms, new threshold levels (independent of the pixel size)
+       have been determined for the pulse positions
+       (50phe for pulse, hi-gain and 320phe for lo-gain)
+
+
 
  2006/08/01 Thomas Bretz

trunk/MagicSoft/Mars/NEWS

@@ -37 +37 @@
 
    - callisto: Implemented two new tabs, the pulse position of pulses
-     not saturating the hi-gain (and above 15phe) and the pulse-position
+     not saturating the hi-gain (and above 50phe) and the pulse-position
      saturating the hi-gain (means the pulse position of signal extracted
-     from the lo-gain) This can be used to check (and/or correct) the
-     hi-/lo-gain offset.
+     from the lo-gain and above 320phe) This can be used to check (and/or
+     correct) the hi-/lo-gain offset.
+
+   - callisto: After extensive tests a new threshold level for determination
+     of the pulse position has been found at 50phe. This is due to the
+     fact that for lower pulses the arrival time changes its average
+     behaviour dramatically. It tends to the average of the extraction
+     range and starts sticking to non-floating point numbers (5, 6, 7, ...)
 
    - star: changed the fit for the ffective on time such that initial

trunk/MagicSoft/Mars/mhist/MHCamEvent.cc

@@ -54 +54 @@
 #include "MHCamEvent.h"
 
-#include <TCanvas.h>
-#include <TPaveStats.h>
-
 #include "MLog.h"
 #include "MLogManip.h"

trunk/MagicSoft/Mars/msignal/MExtractTimeAndCharge.cc

@@ -226 +226 @@
 
           if ( fLoGainFirst <= fLoGainLast-fWindowSizeLoGain)
-            {
+          {
               const Bool_t logainabflag = (higainabflag + pixel.GetNumHiGainSamples()) & 0x1;
               FindTimeAndChargeLoGain(pixel.GetLoGainSamples()+fLoGainFirst,

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilter.cc

@@ -88 +88 @@
 const Int_t  MExtractTimeAndChargeDigitalFilter::fgSignalStartBinHiGain    =  4;
 const Int_t  MExtractTimeAndChargeDigitalFilter::fgSignalStartBinLoGain    =  4;
-const Float_t MExtractTimeAndChargeDigitalFilter::fgOffsetLoGain           =  1.7; // 5 ns
+const Float_t MExtractTimeAndChargeDigitalFilter::fgOffsetLoGain           =  1.05;
 const Float_t MExtractTimeAndChargeDigitalFilter::fgLoGainStartShift       = -1.8;
 
@@ -113 +113 @@
     SetSignalStartBin();
 
-    SetOffsetLoGain(fgOffsetLoGain);
-    SetLoGainStartShift(fgLoGainStartShift);
+    SetOffsetLoGain(fgOffsetLoGain);          // Order between both
+    SetLoGainStartShift(fgLoGainStartShift);  // is important
 }
 
@@ -480 +480 @@
 
   // here, the first high-gain slice is already included in the fTimeShiftHiGain
-  time = fTimeShiftHiGain + max_p - Float_t(t_iter)/fBinningResolutionHiGain;
+//  time = fTimeShiftHiGain + max_p - Float_t(t_iter)/fBinningResolutionHiGain;
+  time = max_p + fTimeShiftHiGain + (Float_t)fHiGainFirst /* this shifts the time to the start of the rising edge */
+      - ((Float_t)t_iter)/fBinningResolutionHiGain;
 
   const Float_t timefineadjust = time_sum/sum;

trunk/MagicSoft/Mars/msignal/MSignalCam.cc

@@ -643 +643 @@
         break;
 
-    case 7: // pulse position above 15phe
+    case 7: // pulse position above 50phe
         // The number of photons is not scaled with the ratio because
         // otherwise to many large pixels survive (maybe because the
         // fluctuations scale different than expected)
-        if (pix->IsPixelUnmapped() || pix->GetNumPhotons()*ratio<15)
+        if (pix->IsPixelUnmapped() || pix->GetNumPhotons()<50)
             return kFALSE;
         val = pix->GetArrivalTime();
@@ -657 +657 @@
 
     case 10: // lo gain time
-        if (pix->IsPixelUnmapped() || !pix->IsLoGainUsed())
+        if (pix->IsPixelUnmapped() || !pix->IsLoGainUsed() ||
+            pix->GetNumPhotons()<320)
             return kFALSE;
         val = pix->GetArrivalTime();
@@ -667 +668 @@
         // fluctuations scale different than expected)
         if (pix->IsPixelUnmapped() || pix->IsLoGainUsed() ||
-            pix->GetNumPhotons()*ratio<15)
+            pix->GetNumPhotons()<50)
             return kFALSE;
         val = pix->GetArrivalTime();