Changeset 9318 for trunk/MagicSoft/Mars/msimcamera

Timestamp:

02/11/09 17:35:54 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msimcamera

Files:

• MSimBundlePhotons.cc (modified) (2 diffs)
• MSimReadout.cc (modified) (6 diffs)
• MSimReadout.h (modified) (1 diff)
• MSimTrigger.cc (modified) (7 diffs)
• MSimTrigger.h (modified) (1 diff)

trunk/MagicSoft/Mars/msimcamera/MSimBundlePhotons.cc

@@ -89 +89 @@
 
     // Read the look-up table
-    if (fLut.ReadFile(fFileName)<0)
+    fLut.Delete();
+    if (!fFileName.IsNull() && fLut.ReadFile(fFileName)<0)
         return kFALSE;
 
     // If the table is empty remove this task from the tasklist
     if (fLut.IsEmpty())
+    {
+        *fLog << inf << "Look-up table to bundle photons empty... skipping." << endl;
         return kSKIP;
+    }
 
     // Now invert the tablee. Otherwise we have to do a lot of
@@ -103 +107 @@
     if (!fLut.HasConstantLength() && fLut.GetMaxEntries()!=1)
         return kFALSE;
+
+    *fLog << inf << "Using look-up table from " << fFileName << endl;
 
     return kTRUE;

trunk/MagicSoft/Mars/msimcamera/MSimReadout.cc

@@ -31 +31 @@
 // system.
 //
+// You can give a conversion factor from the unitx of your analog signal
+// to the units of your adc. This is a fixed factor because it is just
+// a matter of what the meaning of an adc count is, nothing which could
+// jitter or is a real part of the electronics. Such effects should
+// be simulated somewhere else.
+//
 //
 //  Input Containers:
@@ -73 +79 @@
 //
 MSimReadout::MSimReadout(const char* name, const char *title)
-: fRunHeader(0), fEvtHeader(0), fCamera(0), fPulsePos(0), fTrigger(0), fData(0)
+    : fRunHeader(0), fEvtHeader(0), fCamera(0), fPulsePos(0), fTrigger(0), fData(0),
+    fConversionFactor(1)
 {
     fName  = name  ? name  : "MSimReadout";
@@ -135 +142 @@
 
     fRunHeader->InitPixels(cam->GetNumPixels());
-    fRunHeader->SetValidMagicNumber();
-    fRunHeader->SetSourceInfo("MonteCarlo");  // "Muon" from first event??
-    fRunHeader->SetReadyToSave();
 
     fData->InitRead(fRunHeader);
@@ -209 +213 @@
     }
 
-    const Float_t gain      = 64./1;
-    const Float_t offset    = 128;
+    const Float_t offset    = 0;//128;
     const UInt_t  max       = fData->GetMax();
 
@@ -230 +233 @@
         {
             Float_t slice = j+trig>=(Int_t)sig.GetSize() ? offset :
-                sig[j+trig] * gain + offset;
+                sig[j+trig] * fConversionFactor + offset;
 
             // FIXME: Handle/Implement saturation!
@@ -257 +260 @@
     return kTRUE;
 }
+
+// --------------------------------------------------------------------------
+//
+// Read the parameters from the resource file.
+//
+//  ConversionFactor: 1
+//
+Int_t MSimReadout::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
+{
+    Bool_t rc = kFALSE;
+    if (IsEnvDefined(env, prefix, "ConversionFactor", print))
+    {
+        rc = kTRUE;
+        fConversionFactor = GetEnvValue(env, prefix, "ConversionFactor", fConversionFactor);
+    }
+
+    return rc;
+}

trunk/MagicSoft/Mars/msimcamera/MSimReadout.h

@@ -24 +24 @@
     MRawEvtData      *fData;       //! Digitized FADC signal
 
+    Double_t fConversionFactor;    // Conversion factor (arbitrary) from analog signal to FADC counts
+
     // MTask
     Int_t  PreProcess(MParList *pList);
     Int_t  Process();
     Bool_t ReInit(MParList *pList);
+
+    // MParContainer
+    Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print=kFALSE);
 
 public:

trunk/MagicSoft/Mars/msimcamera/MSimTrigger.cc

@@ -82 +82 @@
 #include "MTriggerPattern.h"
 
+#include "MPedestalCam.h"
+#include "MPedestalPix.h"
+
 ClassImp(MSimTrigger);
 
@@ -91 +94 @@
 //
 MSimTrigger::MSimTrigger(const char *name, const char *title)
-    : fCamera(0), fPulsePos(0), fTrigger(0), fRunHeader(0), fEvtHeader(0),
-    fDiscriminatorThreshold(-1), fDigitalSignalLength(8), fCoincidenceTime(0.5)
+    : 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)
 {
     fName  = name  ? name  : "MSimTrigger";
@@ -166 +171 @@
     }
 
+    fElectronicNoise = 0;
+    if (fShiftBaseline)
+    {
+        fElectronicNoise = (MPedestalCam*)pList->FindObject("ElectronicNoise", "MPedestalCam");
+        if (!fElectronicNoise)
+        {
+            *fLog << err << "ElectronicNoise [MPedestalCam] not found... aborting." << endl;
+            return kFALSE;
+        }
+        *fLog << inf << "Baseline will be shifted back to 0 for discriminator." << endl;
+    }
+
+    fGain = 0;
+    if (fUngainSignal)
+    {
+        fGain = (MPedestalCam*)pList->FindObject("Gain", "MPedestalCam");
+        if (!fGain)
+        {
+            *fLog << err << "Gain [MPedestalCam] not found... aborting." << endl;
+            return kFALSE;
+        }
+        *fLog << inf << "Discriminator will be multiplied by applied gain." << endl;
+    }
+
     fRunHeader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
     if (!fRunHeader)
@@ -177 +206 @@
         return kFALSE;
 
-
-    if (fRouteAC.ReadFile(fNameRouteAC)<0)
-        return kFALSE;
-
-    if (fCoincidenceMap.ReadFile(fNameCoincidenceMap)<0)
-        return kFALSE;
+    fRouteAC.Delete();
+    if (!fNameRouteAC.IsNull() && fRouteAC.ReadFile(fNameRouteAC)<0)
+        return kFALSE;
+
+    fCoincidenceMap.Delete();
+    if (!fNameCoincidenceMap.IsNull() && fCoincidenceMap.ReadFile(fNameCoincidenceMap)<0)
+        return kFALSE;
+
+    // ---------------- Consistency checks ----------------------
+
+    if (!fRouteAC.IsEmpty() && !fCoincidenceMap.IsEmpty() &&
+        fCoincidenceMap.GetMaxIndex()>fRouteAC.GetNumRows()-1)
+    {
+        *fLog << err;
+        *fLog << "ERROR - AC routing produces " << fRouteAC.GetNumRows() << " analog channels," << endl;
+        *fLog << "        but the coincidence map expects at least " << fCoincidenceMap.GetMaxIndex()+1 << " channels." << endl;
+        return kERROR;
+    }
 
     if (fDiscriminatorThreshold<=0)
@@ -190 +231 @@
     }
 
-    *fLog << inf << "Using discriminator threshold of " << fDiscriminatorThreshold << endl;
+    if (fElectronicNoise && !fRouteAC.IsEmpty() && !fRouteAC.IsDefault())
+    {
+        // FIXME: Apply to analog channels when summing
+        *fLog << warn << "WARNING - A baseline shift doesn't make sense for sum-channels... reset." << endl;
+        fElectronicNoise = 0;
+    }
+
+    if (fGain && !fRouteAC.IsEmpty() && !fRouteAC.IsDefault())
+    {
+        // FIXME: Apply to analog channels when summing
+        *fLog << warn << "WARNING - Ungain doesn't make sense for sum-channels... reset." << endl;
+        fGain = 0;
+    }
+
+
+    // ---------------- Information output ----------------------
+
+    *fLog << inf;
+
+    if (fRouteAC.IsEmpty())
+        *fLog << "Re-routing/summing of analog channels before discriminator switched off." << endl;
+    else
+        *fLog << "Using " << fNameRouteAC << " for re-routing/summing of analog channels before discriminator." << endl;
+
+    if (fCoincidenceMap.IsEmpty())
+        *fLog << "No coincidences of digital channels will be checked. Signal-above-threshold trigger applied." << endl;
+    else
+        *fLog << "Using " << fNameCoincidenceMap << " to check for coincidences of the digital channels." << endl;
+
+    *fLog << "Using discriminator threshold of " << fDiscriminatorThreshold << endl;
 
     return kTRUE;
@@ -204 +274 @@
     // ================== Simulate channel bundling ====================
 
-    const UInt_t npatch = fRouteAC.GetEntriesFast();
-
-    MAnalogChannels patches(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++)
-        {
-            // FIXME: Simulate clipping
-            const UInt_t idx = row[j];
-            patches[i].AddSignal((*fCamera)[idx]);
+    // FIXME: Before we can bundle the channels we have to make a copy
+    //        and simulate clipping
+
+    // Check if routing should be done
+    const Bool_t empty = fRouteAC.IsEmpty();
+
+    // If no channels are summed the number of patches stays the same
+    const UInt_t npatch = empty ? fCamera->GetNumChannels() : fRouteAC.GetEntriesFast();
+
+    // Use the given analog channels as default out. If channels are
+    // summed overwrite with a newly allocated set of analog channels
+    MAnalogChannels *patches = fCamera;
+    if (!empty)
+    {
+        // FIXME: Can we add gain and offset here into a new container?
+
+        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++)
+            {
+                const UInt_t idx = row[j];
+                (*patches)[i].AddSignal((*fCamera)[idx]);
+            }
         }
     }
@@ -227 +310 @@
 
     for (UInt_t i=0; i<npatch; i++)
-        ttls.AddAt(patches[i].Discriminate(fDiscriminatorThreshold, fDigitalSignalLength), i);
+    {
+        // FIXME: What if the gain was also allpied to the baseline?
+        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);
+    }
 
     // FIXME: Write TTLs!
 
+    // If analog channels had been newly allocated free memmory
+    if (patches!=fCamera)
+        delete patches;
+
     // =================== Simulate coincidences ======================
+
+    // If the map is empty we create a one-pixel-coincidence map
+    // FIMXE: This could maybe be accelerated if the Clone can be
+    //        omitted in the loop
+    if (fCoincidenceMap.IsEmpty())
+        fCoincidenceMap.SetDefault(npatch);
 
     // Calculate the minimum and maximum time for a valid trigger

trunk/MagicSoft/Mars/msimcamera/MSimTrigger.h

@@ -15 +15 @@
 class MRawEvtHeader;
 class MRawRunHeader;
+class MPedestalCam;
 
 class MSimTrigger : public MTask
 {
 private:
-    MAnalogChannels *fCamera;         //! The analog input channels
-    MParameterD     *fPulsePos;       //! The intended pulse positon
-    MParameterD     *fTrigger;        //! The trigger position w.r.t. the analog channels
-    MRawRunHeader   *fRunHeader;      //! The run header storing infos about the digitization
-    MRawEvtHeader   *fEvtHeader;      //! The event header storing the trigger information
+    MAnalogChannels *fCamera;           //! The analog input channels
+    MParameterD     *fPulsePos;         //! The intended pulse positon
+    MParameterD     *fTrigger;          //! The trigger position w.r.t. the analog channels
+    MRawRunHeader   *fRunHeader;        //! The run header storing infos about the digitization
+    MRawEvtHeader   *fEvtHeader;        //! The event header storing the trigger information
+    MPedestalCam    *fElectronicNoise;  //! Electronic noise (for baseline correction)
+    MPedestalCam    *fGain;             //! Gain of the pulses
 
-    MLut fRouteAC;                    // Combinination map for the AC channels
-    MLut fCoincidenceMap;             // channels for which digital coincidence is checked
+    MLut fRouteAC;                      // Combinination map for the AC channels
+    MLut fCoincidenceMap;               // channels for which digital coincidence is checked
 
-    TString fNameRouteAC;             // Name for the AC routing
-    TString fNameCoincidenceMap;      // Name for the coincidence mape
+    TString fNameRouteAC;               // Name for the AC routing
+    TString fNameCoincidenceMap;        // Name for the coincidence mape
 
-    Float_t fDiscriminatorThreshold;  // Discriminator threshold
-    Float_t fDigitalSignalLength;     // Length of the output of the discriminator
-    Float_t fCoincidenceTime;         // Minimum coincidence time (gate)
+    Float_t fDiscriminatorThreshold;    // Discriminator threshold
+    Float_t fDigitalSignalLength;       // Length of the output of the discriminator
+    Float_t fCoincidenceTime;           // Minimum coincidence time (gate)
+
+    Bool_t  fShiftBaseline;             // Shift the baseline back to 0 for the threshold (needs ElectronicNoise [MPedestalCam])
+    Bool_t  fUngainSignal;              // "Remove" the gain from the signal (needs Gain [MPedestalCam])
 
     // MSimTrigger