Changeset 5544 for trunk/MagicSoft/Mars

Timestamp:

12/01/04 20:49:05 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mhcalib/MHGausEvents.cc (modified) (2 diffs)
• msignal/MExtractTimeAndChargeSpline.cc (modified) (13 diffs)
• msignal/MExtractTimeAndChargeSpline.h (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -22 +22 @@
 
  2004/12/01: Markus Gaug
+
+   * MExtractTimeAndChargeSpline.[h,cc]:
+     - introduced a random iterator for the "pedestal" calculation of 
+       the amplitude spline.
+
 
    * MExtractTimeAndChargeDigitalFilter.cc:

trunk/MagicSoft/Mars/mhcalib/MHGausEvents.cc

@@ -487 +487 @@
   
   Int_t win = 1;
+  Int_t nofit = 0;
 
   if (option.Contains("events"))
-    win += 1;
+    {
+      option.ReplaceAll("events","");      
+      win += 1;
+    }
   if (option.Contains("fourier"))
-    win += 2;
+    {
+      option.ReplaceAll("fourier","");      
+      win += 2;
+    }
+  
   if (IsEmpty())
     win--;
 
+  if (option.Contains("nofit"))
+    {
+      option.ReplaceAll("nofit","");
+      nofit++;
+    }
+  
   pad->SetBorderMode(0);
   pad->Divide(1,win);
@@ -510 +524 @@
 
       fHGausHist.Draw(option);
-      
-      if (fFGausFit)
-        {
-          fFGausFit->SetLineColor(IsGausFitOK() ? kGreen : kRed);
-          fFGausFit->Draw("same");
-        }
+
+      if (!nofit)
+        if (fFGausFit)
+          {
+            fFGausFit->SetLineColor(IsGausFitOK() ? kGreen : kRed);
+            fFGausFit->Draw("same");
+          }
     }
   

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSpline.cc

@@ -177 +177 @@
 //
 MExtractTimeAndChargeSpline::MExtractTimeAndChargeSpline(const char *name, const char *title) 
-    : fAbMax(0.), fAbMaxPos(0.), fHalfMax(0.)
+    : fAbMax(0.), fAbMaxPos(0.), fHalfMax(0.), fRandomIter(0)
 {
 
@@ -251 +251 @@
 // --------------------------------------------------------------------------
 //
-// ReInit
-//
-// Calls:
-// - MExtractTimeAndCharge::ReInit(pList);
-// - Deletes all arrays, if not NULL
-// - Creates new arrays according to the extraction range
-//
-Bool_t MExtractTimeAndChargeSpline::ReInit(MParList *pList)
+// InitArrays
+//
+// Gets called in the ReInit() and initialized the arrays
+//
+Bool_t MExtractTimeAndChargeSpline::InitArrays()
 {
 
+  Int_t range = fHiGainLast - fHiGainFirst + 1 + fHiLoLast;
+
+  fHiGainSignal     .Set(range);
+  fHiGainFirstDeriv .Set(range);
+  fHiGainSecondDeriv.Set(range);
+
+  range = fLoGainLast - fLoGainFirst + 1;
+
+  fLoGainSignal     .Set(range);
+  fLoGainFirstDeriv .Set(range);
+  fLoGainSecondDeriv.Set(range);
+
+  fHiGainSignal     .Reset();
+  fHiGainFirstDeriv .Reset();
+  fHiGainSecondDeriv.Reset();
+
+  fLoGainSignal     .Reset();
+  fLoGainFirstDeriv .Reset();
+  fLoGainSecondDeriv.Reset();
+  
   if (IsExtractionType(kAmplitude))
     {
@@ -280 +297 @@
       fWindowSizeLoGain  = (Int_t)(fRiseTime + fFallTime+1);
     }
-
-  if (!MExtractTimeAndCharge::ReInit(pList))
-    return kFALSE;
-
-  return kTRUE;
-}
-
-// --------------------------------------------------------------------------
-//
-// InitArrays
-//
-// Gets called in the ReInit() and initialized the arrays
-//
-Bool_t MExtractTimeAndChargeSpline::InitArrays()
-{
-
-  Int_t range = fHiGainLast - fHiGainFirst + 1 + fHiLoLast;
-
-  fHiGainSignal     .Set(range);
-  fHiGainFirstDeriv .Set(range);
-  fHiGainSecondDeriv.Set(range);
-
-  range = fLoGainLast - fLoGainFirst + 1;
-
-  fLoGainSignal     .Set(range);
-  fLoGainFirstDeriv .Set(range);
-  fLoGainSecondDeriv.Set(range);
-
-  fHiGainSignal     .Reset();
-  fHiGainFirstDeriv .Reset();
-  fHiGainSecondDeriv.Reset();
-
-  fLoGainSignal     .Reset();
-  fLoGainFirstDeriv .Reset();
-  fLoGainSecondDeriv.Reset();
 
   return kTRUE;
@@ -415 +397 @@
   if (IsNoiseCalculation() && IsExtractionType(kAmplitude))
     {
-      //
-      // Take the spline value at the middle of the third slice (to avoid egde effects)
-      // 
-      sum = 0.5*fHiGainSignal[2]
-          + 0.5*fHiGainSignal[3]
-          + (-0.375)*fHiGainSecondDeriv[2] 
-          + (-0.375)*fHiGainSecondDeriv[3];
+      if (fRandomIter == (TMath::Floor(1./fResolution)))
+        fRandomIter = 0;
+
+      const Float_t b = fRandomIter * fResolution;
+      const Float_t a = 1. - b;
+
+      fRandomIter++;
+
+      sum = a*fHiGainSignal[1]
+          + b*fHiGainSignal[2]
+          + (a*a*a-a)*fHiGainSecondDeriv[1] 
+          + (b*b*b-b)*fHiGainSecondDeriv[2];
       return;
     }
@@ -531 +518 @@
   // Try a better precision. 
   //
-  const Float_t up = fAbMaxPos+step-0.035;
-  const Float_t lo = fAbMaxPos-step+0.035;
+  const Float_t up = fAbMaxPos+step - 1.5*fResolution;
+  const Float_t lo = fAbMaxPos-step + 1.5*fResolution;
   const Float_t maxpossave = fAbMaxPos;
   
@@ -539 +526 @@
   b     = x - lower;
  
-  step  = 0.025; // step size of 83 ps 
+  step  = fResolution; // step size of 83 ps 
  
   while (x<up)
@@ -562 +549 @@
 
   //
-  // Second, try from time down to time-0.2 in steps of 0.025.
+  // Second, try from time down to time-0.2 in steps of fResolution.
   //
   x     = maxpossave;
@@ -571 +558 @@
   // which requires new setting of klocont and khicont
   //
-  if (x < klo + 0.02)
+  if (x < klo + fResolution/2.)
     {
       klo--;
@@ -605 +592 @@
     {
       time = (Float_t)fHiGainFirst + fAbMaxPos;
-      dtime = 0.025;
+      dtime = fResolution;
     }
   else
@@ -892 +879 @@
   // Try a better precision. 
   //
-  const Float_t up = fAbMaxPos+step-0.035;
-  const Float_t lo = fAbMaxPos-step+0.035;
+  const Float_t up = fAbMaxPos+step - 1.5*fResolution;
+  const Float_t lo = fAbMaxPos-step + 1.5*fResolution;
   const Float_t maxpossave = fAbMaxPos;
   
@@ -900 +887 @@
   b     = x - lower;
  
-  step  = 0.025; // step size of 83 ps 
+  step  = fResolution; // step size of fResolution (33 ps )
  
   while (x<up)
@@ -932 +919 @@
   // which requires new setting of klocont and khicont
   //
-  if (x < klo + 0.02)
+  if (x < klo + fResolution/2.)
     {
       klo--;
@@ -966 +953 @@
     {
       time = (Float_t)fLoGainFirst + fAbMaxPos;
-      dtime = 0.02;
+      dtime = fResolution;
     }
   else

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSpline.h

@@ -39 +39 @@
   Float_t fFallTime;                     // The usual fall time of the pulse
 
+  UInt_t  fRandomIter;                   // Counter used to randomize weights for noise calculation
+  
   Byte_t  fFlags;                        // Bit-field to hold the time extraction types
-  
-  Bool_t ReInit( MParList *pList );
   
   Int_t  ReadEnv(const TEnv &env, TString prefix, Bool_t print);