Changeset 18479 for trunk/Mars

Timestamp:

04/30/16 18:38:36 (9 years ago)

Author:

tbretz

Message:

After some study this seems right now to be the most consistent implementation.

File:

• trunk/Mars/manalysis/MSoftwareTrigger.cc (modified) (5 diffs)

trunk/Mars/manalysis/MSoftwareTrigger.cc

@@ -106 +106 @@
     MArrayF buf(160*num);
 
+    double avg = 0;
+
     for (int hw=0; hw<1440; hw++)
     {
@@ -121 +123 @@
             continue;
 
-        const UInt_t sw = idx[hw];
-
-        const Float_t *raw = fSignal->GetSamples(sw);
+        const Float_t *raw = fSignal->GetSamples(idx[hw]);
 
         Float_t *sum = buf.GetArray()+(hw/9)*num;
@@ -154 +154 @@
 */
 
-    Float_t max = 0;
-    double avg = 0;
+    Float_t max = -50000;
 
     const UInt_t nsum = buf.GetSize();
@@ -163 +162 @@
         Float_t v[4] = { 0, 0, 0, 0 };
 
-        for (Float_t *ptr=sum+20; ptr<sum+num; ptr++)
+        for (Float_t *ptr=sum+15; ptr<sum+num; ptr++)
         {
-            *ptr -= 0.6 * ptr[-20];
-            avg += *ptr;
-            v[idx++%4] = *ptr;
+            ptr[0] -= 0.5*ptr[-15];
+
+            avg += ptr[0];
+            v[idx++%4] = ptr[0];
 
             const Float_t min = *std::min_element(v, v+4);
@@ -175 +175 @@
     }
 
-    avg /= num*(1440-9);  // excluding the 0's from the broken DRS board
+    // Question: Should we also keep position and patch id?
+
+    avg /= (num-15)*(1440-24);
 
     fTriggerSignal->SetVal(max);