Changeset 18009 for trunk/Mars

Timestamp:

11/07/14 15:41:06 (10 years ago)

Author:

smueller

Message:

New feature: Fix temporal offsets in between the pixels can be simulated and defined in a text file. This is solving ticket #9

Location:

trunk/Mars

Files:

• . (modified) (1 prop)
• ceres.rc (modified) (1 diff)
• mfileio/FileIOLinkDef.h (modified) (1 diff)
• mfileio/MMatrix.cc (copied) (copied from branches/MarsFixTimeOffsetBranch/mfileio/MMatrix.cc )
• mfileio/MMatrix.h (copied) (copied from branches/MarsFixTimeOffsetBranch/mfileio/MMatrix.h )
• mfileio/Makefile (modified) (1 diff)
• mjobs/MJSimulation.cc (modified) (2 diffs)
• msimcamera/MSimCamera.cc (modified) (3 diffs)
• msimcamera/MSimCamera.h (modified) (2 diffs)
• resmc/fact/AllPhidoFiles_delays.csv (copied) (copied from branches/MarsFixTimeOffsetBranch/resmc/fact/AllPhidoFiles_delays.csv )
• resmc/fact/pixel_delays_ALL_ZERO.csv (copied) (copied from branches/MarsFixTimeOffsetBranch/resmc/fact/pixel_delays_ALL_ZERO.csv )

trunk/Mars/ceres.rc

@@ -233 +233 @@
 # Does not trigger anyway
 ContEmpty3.Condition: MPhotonEvent.GetNumPhotons<10
+
+MFixTimeOffset.FileName: resmc/fact/pixel_delays_ALL_ZERO.csv
+# MFixTimeOffset.FileName: resmc/fact/AllPhidoFiles_delays.csv
+
+# last line comment

trunk/Mars/mfileio/FileIOLinkDef.h

@@ -19 +19 @@
 #pragma link C++ class MWriteRootFile+;
 #pragma link C++ class MWriteFitsFile+;
+#pragma link C++ class MMatrix+;
 
 #endif

trunk/Mars/mfileio/Makefile

@@ -33 +33 @@
            MWriteFitsFile.cc \
            MTopFitsGroup.cc \
-           MFitsArray.cc
+           MFitsArray.cc \
+           MMatrix.cc
 
 ############################################################

trunk/Mars/mjobs/MJSimulation.cc

@@ -113 +113 @@
 #include "MParSpline.h"
 #include "MGeomCam.h"
+#include "MMatrix.h"
 
 #include "MPedestalCam.h"
@@ -718 +719 @@
     trigger.SetVal(0);
     plist.AddToList(&trigger);
+    
+    // -------------------------------------------------------------------
+    // Dominik and Sebastian on: fix time offsets
+    MMatrix fix_time_offsets_between_pixels_in_ns(
+        "MFixTimeOffset","titel"
+    );
+    plist.AddToList(&fix_time_offsets_between_pixels_in_ns);
 
     // -------------------------------------------------------------------

trunk/Mars/msimcamera/MSimCamera.cc

@@ -122 +122 @@
         return kFALSE;
     }
+    // -------------------------------------------------------------------
+    // Dominik Neise and Sebastian Mueller on fix time offsets:
+    // We obtain the fix temporal offsets for the FACT camera pixels out of 
+    // a text file. The textfile must be mentioned in the ceres.rc file.
+    // There are no default offsets on purporse. The filename must be specified
+    // in ceres.rc and the file must be parsed without errors and it must 
+    // provide exactly 1440 floating point numbers.   
+    fFixTimeOffsetsBetweenPixelsInNs = 
+    (MMatrix*)pList->FindObject("MFixTimeOffset");
+    if (!fFixTimeOffsetsBetweenPixelsInNs)
+    {
+        // the key value pair providing the text file is not present in the
+        // environment env.
+        *fLog << err << "In Source: "<< __FILE__ <<" in line: "<< __LINE__;
+        *fLog << " in function: "<< __func__ <<"\n";
+        *fLog << "MFixTimeOffset not found... aborting." << endl;
+        return kFALSE;
+
+    }
+    else if ( fFixTimeOffsetsBetweenPixelsInNs->fM.size() != 1440 )
+    {
+        // The number of time offsets must match the number of pixels in the
+        // FACT camera.
+        *fLog << err << "In Source: "<< __FILE__ <<" in line: "<< __LINE__;
+        *fLog << " in function: "<< __func__ <<"\n";
+        *fLog << "MFixTimeOffset has the wrong dimension! ";
+        *fLog << "There should be "<< 1440 <<" time offsets ";
+        *fLog << "(one for each pixel in FACT) but there are: ";
+        *fLog << fFixTimeOffsetsBetweenPixelsInNs->fM.size() << "! ";
+        *fLog << "... aborting." << endl;
+        return kFALSE;
+    }
+    // Check all entries for inf and nan. Those are not accepted here.
+    for( std::vector< double > row : fFixTimeOffsetsBetweenPixelsInNs->fM ){
+        for( double number : row){
+
+            if( std::isnan(number) || std::isinf(number) ){
+
+                *fLog << err << "In Source: "<< __FILE__ <<" in line: ";
+                *fLog << __LINE__;
+                *fLog << " in function: "<< __func__ <<"\n";
+                *fLog << "There is a non normal number in the fix temporal ";
+                *fLog << "pixel offsets. This is at least one number is ";
+                *fLog << "NaN or Inf. This here is >"<< number;
+                *fLog << "<... aborting." << endl;               
+                return kFALSE;
+            }
+        }
+    }
+    // -------------------------------------------------------------------
 /*
     fPulsePos = (MParameterD*)pList->FindObject("IntendedPulsePos", "MParameterD");
@@ -353 +403 @@
     }
 
+    //--------------------------------------------------------------------------
+
     // Simulate pulses
     for (Int_t i=0; i<num; i++)
@@ -359 +411 @@
 
         const UInt_t   idx = ph.GetTag();
-        const Double_t t   = (ph.GetTime()-fStat->GetTimeFirst())*freq+rndm;// - fSpline->GetXmin();
+        Double_t t   = (ph.GetTime()-fStat->GetTimeFirst())*freq+rndm;// - fSpline->GetXmin();
+
+        // Sebastian Mueller and Dominik Neise on fix time offsets:
+        // We add a fix temporal offset to the relative arrival time of the 
+        // individual pixel. The offsets are stored in the
+        // fFixTimeOffsetsBetweenPixelsInNs -> fM matrix. We identify the first
+        // column to hold the offsets in ns.
+        t = t + freq*fFixTimeOffsetsBetweenPixelsInNs->fM[idx][0];
 
         // FIXME: Time jitter?

trunk/Mars/msimcamera/MSimCamera.h

@@ -7 +7 @@
 
 #include "MArrayF.h"
+#include "MMatrix.h"
 
 class MMcEvt;
@@ -35 +36 @@
 
     MParameterD       *fCrosstalkCoeffParam;
+
+    MMatrix *fFixTimeOffsetsBetweenPixelsInNs; //! Container to store the fix temporal offsets for each pixel in ns
+
     MTruePhotonsPerPixelCont    *fTruePhotons;  //! Container to store the number of photons per pixel