Changeset 18038 for trunk/Mars/msimcamera/MSimCamera.cc

Timestamp:

12/15/14 09:45:18 (10 years ago)

Author:

Jens Buss

Message:

Reintegrating branch for feature ResidualTimeSpread into Branch

Location:

trunk/Mars

Files:

• . (modified) (1 prop)
• msimcamera/MSimCamera.cc (modified) (4 diffs)

trunk/Mars/msimcamera/MSimCamera.cc

@@ -180 +180 @@
     }
  */
+    fResidualTimeSpread = (MParameterD*)pList->FindObject("ResidualTimeSpread");
+    if (!fResidualTimeSpread)
+    {
+        *fLog << err << "ResidualTimeSpread [MParameterD] not found... aborting." << endl;
+        return kFALSE;
+    }
 
     // Create it here to make sure that MGeomApply will set the correct size
@@ -305 +311 @@
     fCamera->SetValidRange(TMath::FloorNint(pr), TMath::CeilNint(nlen+pl));
 
+    Double_t timeoffset[npix];
+
+
     // Add electronic noise to empty channels
     for (UInt_t i=0; i<npix; i++)
     {
+
+        // Get the ResidualTimeSpread Parameter
+        const Double_t residualTimeSpread = fResidualTimeSpread->GetVal();
+
+        // Jens Buss on residual time spread:
+        // randomly draw an additional time offset to be added to the arrivaltime 
+        // from a gaussian normal distribution with a given standard deviation 
+        timeoffset[i] = gRandom->Gaus(0.0, residualTimeSpread);
         const MPedestalPix &pix = (*fElectronicNoise)[i];
 
@@ -405 +422 @@
     //--------------------------------------------------------------------------
 
+    
     // Simulate pulses
     for (Int_t i=0; i<num; i++)
@@ -419 +437 @@
         // column to hold the offsets in ns.
         t = t + freq*fFixTimeOffsetsBetweenPixelsInNs->fM[idx][0];
+
+        // Jens Buss on residual time spread:
+        // add random time offset to the arrivaltimes
+        t = t + timeoffset[idx];
 
         // FIXME: Time jitter?