Changeset 5507 for trunk/MagicSoft/Mars/mpedestal

Timestamp:

11/29/04 20:58:20 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mpedestal

Files:

• MExtractPedestal.cc (modified) (5 diffs)
• MPedCalcFromLoGain.cc (modified) (2 diffs)
• MPedCalcPedRun.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/mpedestal/MExtractPedestal.cc

@@ -151 +151 @@
 #include "MExtractPedestal.h"
 #include "MExtractTimeAndCharge.h"
+#include "MExtractTimeAndChargeDigitalFilter.h"
 
 #include "MParList.h"
@@ -420 +421 @@
   
   if (fExtractor)
-    fExtractor->InitArrays();
-  
+    {
+      fExtractor->InitArrays();
+      SetExtractWindow(fExtractor->GetHiGainFirst(),(Int_t)fExtractor->GetNumHiGainSamples());
+      if (fExtractor->InheritsFrom("MExtractTimeAndChargeDigitalFilter"))
+        {
+          MExtractTimeAndChargeDigitalFilter* df = (MExtractTimeAndChargeDigitalFilter*)fExtractor;
+          if (!df->IsWeightsSet())
+            if (!df->ReadWeightsFile(df->GetNameWeightsFile()))
+              return kFALSE;
+        }
+    }
+
   Print();
 
@@ -510 +521 @@
 
   // 4. Scale the mean, variance and AB-noise to the number of slices:
-  ped    /= fExtractWinSize;
-  var    /= fExtractWinSize;
-  abOffs /= fExtractWinSize;
+  ped    /= fExtractor ? fExtractor->GetNumHiGainSamples() : fExtractWinSize;
+  var    /= fExtractor ? fExtractor->GetNumHiGainSamples() : fExtractWinSize;
+  abOffs /= fExtractor ? fExtractor->GetNumHiGainSamples() : fExtractWinSize;
   
   // 5. Calculate the RMS from the Variance:
@@ -549 +560 @@
   // 4. Scale the mean, variance and AB-noise to the number of slices:
   //
-  ped    /= fExtractWinSize;
-  var    /= fExtractWinSize;
-  abOffs /= fExtractWinSize;
+  ped    /= fExtractor ? fExtractor->GetNumHiGainSamples() : fExtractWinSize;
+  var    /= fExtractor ? fExtractor->GetNumHiGainSamples() : fExtractWinSize;
+  abOffs /= fExtractor ? fExtractor->GetNumHiGainSamples() : fExtractWinSize;
   //
   // 5. Scale the mean, variance and AB-noise to the number of pixels:
@@ -596 +607 @@
   // 4. Scale the mean, variance and AB-noise to the number of slices:
   //
-  ped    /= fExtractWinSize;
-  var    /= fExtractWinSize;
-  abOffs /= fExtractWinSize;
+  ped    /= fExtractor ? fExtractor->GetNumHiGainSamples() : fExtractWinSize;
+  var    /= fExtractor ? fExtractor->GetNumHiGainSamples() : fExtractWinSize;
+  abOffs /= fExtractor ? fExtractor->GetNumHiGainSamples() : fExtractWinSize;
   //
   // 5. Scale the mean, variance and AB-noise to the number of pixels:

trunk/MagicSoft/Mars/mpedestal/MPedCalcFromLoGain.cc

@@ -137 +137 @@
 
 #include "MExtractTimeAndCharge.h"
-#include "MExtractTimeAndChargeDigitalFilter.h"
 
 #include "MParList.h"
@@ -243 +242 @@
 
   if (fExtractor)
-    {
-      SetExtractWindow(fExtractor->GetHiGainFirst(),(Int_t)fExtractor->GetNumHiGainSamples());
-      lastavailable = losamples-1;
-      if (fExtractor->InheritsFrom("MExtractTimeAndChargeDigitalFilter"))
-        {
-          MExtractTimeAndChargeDigitalFilter* df = (MExtractTimeAndChargeDigitalFilter*)fExtractor;
-          if (!df->IsWeightsSet())
-            if (!df->ReadWeightsFile(df->GetNameWeightsFile()))
-              return kFALSE;
-        }
-    }
+    lastavailable = losamples-1;
 
   // If the size is not yet set, set the size

trunk/MagicSoft/Mars/mpedestal/MPedCalcPedRun.cc

@@ -122 +122 @@
 
 #include "MExtractTimeAndCharge.h"
-#include "MExtractTimeAndChargeDigitalFilter.h"
 
 #include "MParList.h"
@@ -204 +203 @@
 {
 
-  if (fExtractor)
-    {
-      SetExtractWindow(fExtractor->GetHiGainFirst(),(Int_t)fExtractor->GetNumHiGainSamples());
-      if (fExtractor->InheritsFrom("MExtractTimeAndChargeDigitalFilter"))
-        {
-          MExtractTimeAndChargeDigitalFilter* df = (MExtractTimeAndChargeDigitalFilter*)fExtractor;
-          if (!df->IsWeightsSet())
-            if (!df->ReadWeightsFile(df->GetNameWeightsFile()))
-              return kFALSE;
-        }
-    }
-  
   MExtractPedestal::ReInit(pList);
 
@@ -470 +457 @@
 {
 
+  if (fRunHeader->GetRunNumber() == 38996)
+    return kTRUE;
+  
   if (fRunHeader->GetRunNumber() < gkFirstRunWithFinalBits)
     return kFALSE;