Changeset 5610 for trunk

Timestamp:

12/16/04 20:14:29 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSpline.cc (modified) (11 diffs)

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSpline.cc

@@ -292 +292 @@
       fNumHiGainSamples  = fRiseTime + fFallTime;
       fNumLoGainSamples  = fLoGainLast ? fNumHiGainSamples + 1. : 0.;
+      fNumLoGainSamples  *= 0.3;      
+
       fSqrtHiGainSamples = TMath::Sqrt(fNumHiGainSamples);
       fSqrtLoGainSamples = TMath::Sqrt(fNumLoGainSamples);
@@ -712 +714 @@
   const Byte_t *end = first + range;
   Byte_t       *p  = first;
-  Int_t  count     = 0;
 
   const Float_t pedes  = ped.GetPedestal();
@@ -724 +725 @@
   fAbMaxPos     = 0.;
   Int_t  maxpos = 0;
+  Int_t  count  = 0;
 
   //
@@ -818 +820 @@
         }
       
+      if (maxpos > range-2)
+        CalcIntegralLoGain(sum, (Float_t)range - fRiseTime - fFallTime-1., (Float_t)range - 0.001);
+      else
+        CalcIntegralLoGain(sum, 0.001, fRiseTime + fFallTime + 1.);
+
+      return;
+    }
+
+  if (maxpos < (Int_t)(fRiseTime+2.))
+    {
+      time =  IsExtractionType(kMaximum) 
+        ? (Float_t)(fLoGainFirst + maxpos) 
+        : (Float_t)(fLoGainFirst + maxpos - 1);
+
       if (maxpos > range-2)
         CalcIntegralLoGain(sum, (Float_t)range - fRiseTime - fFallTime-1., (Float_t)range - 0.001);
@@ -1041 +1057 @@
         }
 
-      time  = x + (Float_t)fLoGainFirst;
+      time  = x + (Int_t)fLoGainFirst;
       dtime = fResolution;
     }
@@ -1067 +1083 @@
         }
       CalcIntegralLoGain(sum, start, last);
-    }
+      //      *fLog << inf << time << "  " << sum << "  " << start << "  " << last << endl;
+    }
+
 }
 
@@ -1084 +1102 @@
   Int_t khi = klo+1;
 
-  Float_t up = TMath::Ceil(start);
   Float_t lo = TMath::Floor(start);
+  Float_t up = lo + 1.;
 
   const Int_t m = int((start-klo)/step);
@@ -1140 +1158 @@
   Int_t khi = klo+1;
 
-  Float_t up = TMath::Ceil(start);
   Float_t lo = TMath::Floor(start);
+  Float_t up = lo + 1.;
 
   const Int_t m = int((start-klo)/step);
@@ -1149 +1167 @@
   Float_t a = up-start;
   Float_t b = start-lo;
-
-  Int_t cnt = 0;
 
   while (1)
@@ -1173 +1189 @@
             + (b*b*b-b)*fHiGainSecondDeriv[khi];
           
-          cnt++;
         }
 
@@ -1184 +1199 @@
       b = 0.;
     }
+  
 }