Changeset 17682 for trunk/Mars/msimcamera/MSimTrigger.cc

Timestamp:

04/23/14 13:53:35 (11 years ago)

Author:

ftemme

Message:

Added the clipping of the trigger signal

File:

• trunk/Mars/msimcamera/MSimTrigger.cc (modified) (10 diffs)

trunk/Mars/msimcamera/MSimTrigger.cc

@@ -97 +97 @@
 //
 MSimTrigger::MSimTrigger(const char *name, const char *title)
-    : fCamera(0), fPulsePos(0), fTrigger(0), fRunHeader(0),
-    fEvtHeader(0), fElectronicNoise(0), fGain(0),
-    fDiscriminatorThreshold(-1), fDigitalSignalLength(8), fCoincidenceTime(0.5),
-    fShiftBaseline(kTRUE), fUngainSignal(kTRUE), fSimulateElectronics(kTRUE),
-    fMinMultiplicity(-1)
+    : fCamera(0),
+      fPulsePos(0),
+      fTrigger(0),
+      fRunHeader(0),
+      fEvtHeader(0),
+      fElectronicNoise(0),
+      fGain(0),
+      fDiscriminatorThreshold(-1),
+      fDigitalSignalLength(8),
+      fCoincidenceTime(0.5),
+      fShiftBaseline(kTRUE),
+      fUngainSignal(kTRUE),
+      fSimulateElectronics(kTRUE),
+      fMinMultiplicity(-1),
+      fCableDelay(21),
+      fCableDamping(0.)     // default Damping Set to zero, so users, who do not set
+                            // the CableDamoing in the ceres.rc do not see a difference.
+
 {
     fName  = name  ? name  : "MSimTrigger";
@@ -385 +398 @@
     *fLog << "Using discriminator threshold of " << fDiscriminatorThreshold << endl;
 
+    *fLog << "Using fCableDelay " << fCableDelay << endl;
+    *fLog << "Using fCableDamping " << fCableDamping << endl;
+
+    ////////////////////////////////
+    // open some output files for debugging
+//    patch_file.open("/home/fact_opr/patch_file.csv", ios_base::out);
+//    clipped_file.open("/home/fact_opr/clipped_file.csv", ios_base::out);
+//    digital_file.open("/home/fact_opr/digital_file.csv", ios_base::out);
+//    ratescan_file.open("/home/fact_opr/ratescan_file.csv", ios_base::out);
+
+
+
     return kTRUE;
 }
@@ -497 +522 @@
     // summed overwrite with a newly allocated set of analog channels
     MAnalogChannels *patches = fCamera;
+    MAnalogChannels *raw_patches = fCamera;
     if (!empty)
     {
@@ -502 +528 @@
 
         patches = new MAnalogChannels(npatch, fCamera->GetNumSamples());
-        for (UInt_t i=0; i<npatch; i++)
-        {
-            const MArrayI &row = fRouteAC.GetRow(i);
-            for (UInt_t j=0; j<row.GetSize(); j++)
+        raw_patches = new MAnalogChannels(npatch, fCamera->GetNumSamples());
+        for (UInt_t patch_id=0; patch_id<npatch; patch_id++)
+        {
+            const MArrayI &row = fRouteAC.GetRow(patch_id);
+            for (UInt_t pxl_in_patch=0; pxl_in_patch<row.GetSize(); pxl_in_patch++)
             {
-                const UInt_t idx = row[j];
+                const UInt_t pixel_id = row[pxl_in_patch];
 
                 // FIXME: Shrinking the mapping table earlier (e.g.
                 //        ReInit) would avoid a lot of if's
-                if (idx<fCamera->GetNumChannels())
-                    (*patches)[i].AddSignal((*fCamera)[idx]);
+                if (pixel_id<fCamera->GetNumChannels())
+                {
+                    (*raw_patches)[patch_id].AddSignal((*fCamera)[pixel_id]);
+                    (*patches)[patch_id].AddSignal((*fCamera)[pixel_id]);
+                    (*patches)[patch_id].AddSignal((*fCamera)[pixel_id], fCableDelay, fCableDamping);
+                }
             }
         }
     }
+
+    // DN: 20140219 Ratescan:
+    //
+    //
+//    for (UInt_t patch_id=0; patch_id<npatch; patch_id++)
+//    {
+//        MAnalogSignal current_patch = (*raw_patches)[patch_id];
+//        float max = current_patch[0];
+//        for (UInt_t i=1; i<current_patch.GetSize(); i++)
+//        {
+//            if (current_patch[i] > max)
+//            {
+//                max = current_patch[i];
+//            }
+//        }
+//        ratescan_file << max << " ";
+//    }
+//    ratescan_file << endl;
+
+//    // DN 20131108: DEBUGGING:
+//    for (UInt_t patch_id=0; patch_id<npatch; patch_id++)
+//    {
+//        for (UInt_t slice=0; slice<fCamera->GetNumSamples(); slice++)
+//        {
+//            patch_file << (*raw_patches)[patch_id][slice] << "\t";
+//            clipped_file << (*patches)[patch_id][slice] << "\t";
+//        }
+//        patch_file << endl;
+//        clipped_file << endl;
+//    }
+
+
 
     // FIXME: Write patches
@@ -529 +592 @@
         const Double_t offset = fElectronicNoise ? (*fElectronicNoise)[i].GetPedestal() : 0;
         const Double_t gain   = fGain            ? (*fGain)[i].GetPedestal()            : 1;
-        ttls.AddAt((*patches)[i].Discriminate(fDiscriminatorThreshold*gain+offset, fDigitalSignalLength), i);
+        ttls.AddAt(
+                    (*patches)[i].Discriminate(
+                        fDiscriminatorThreshold*gain+offset,                // treshold
+                        Double_t(fCableDelay),                              // start
+                        Double_t(fCamera->GetNumSamples() - fCableDelay),   // end
+                        //fDigitalSignalLength                              // time-over-threshold, or fixed-length?
+                        -1                                                  // -1 = time-over-threshold
+                        ),
+                    i
+                    );
     }
 
@@ -537 +609 @@
     if (patches!=fCamera)
         delete patches;
+    if (raw_patches!=fCamera)
+        delete raw_patches;
 
     // =================== Simulate coincidences ======================
@@ -557 +631 @@
     Int_t rmlo = 0;
     Int_t rmhi = 0;
+
+//    cout << "MSimTrigger::fMinMultiplicity = " << fMinMultiplicity << endl;
+//    cout << "MSimTrigger::fCoincidenceTime = " << fCoincidenceTime << endl;
+//    cout << "fCoincidenceMap.GetEntries() = " << fCoincidenceMap.GetEntries() << endl;
+//    cout << "MSimTrigger::fCableDelay = " << fCableDelay << endl;
+//    cout << "MSimTrigger::fCableDamping = " << fCableDamping << endl;
+//    cout << "min:" << min << endl;
+//    cout << "max:" << max << endl;
 
     for (int j=0; j<fCoincidenceMap.GetEntries(); j++)
@@ -624 +706 @@
         arr->Compress();
 
+//        cout << "ttls(j=" << j << "):";
+//        TObjArray *arr2 = static_cast<TObjArray*>(ttls[j]);
+//        TIter Nexty(arr);
+//        MDigitalSignal *ttly = 0;
+//        while ((ttly=static_cast<MDigitalSignal*>(Nexty())))
+//        {
+//            cout << "|"<< ttly->GetStart() << ", " << ttly->GetLength();
+//        }
+//        cout << endl;
+
+
         // If we have at least one trigger keep the earliest one.
         // FIXME: The triggers might be ordered in time automatically:
@@ -674 +767 @@
 }
 
+Int_t MSimTrigger::PostProcess()
+{
+//    patch_file.close();
+//    clipped_file.close();
+//    digital_file.close();
+//    ratescan_file.close();
+    return kTRUE;
+}
+
+
 // --------------------------------------------------------------------------
 //
@@ -728 +831 @@
     }
 
+    if (IsEnvDefined(env, prefix, "CableDelay", print))
+    {
+        rc = kTRUE;
+        fCableDelay = GetEnvValue(env, prefix, "CableDelay", fCableDelay);
+    }
+
+    if (IsEnvDefined(env, prefix, "CableDamping", print))
+    {
+        rc = kTRUE;
+        fCableDamping = GetEnvValue(env, prefix, "CableDamping", fCableDamping);
+    }
+
+
+
     return rc;
 }