set time resolution for the sts pixel digitizer

sim/detectors/sts/CbmStsDigitizePixel.cxx

@@ -37,6 +37,8 @@
 #include "FairRuntimeDb.h"
 #include <Logger.h>
 
+#include <TRandom.h>
+
 // Includes from CbmRoot
 #include "CbmStsDigi.h"
 #include "CbmStsParSetModule.h"
@@ -72,11 +74,11 @@ using namespace CbmSts;
 ClassImp(CbmStsDigitizePixel);
 
 // -----   Standard constructor   ------------------------------------------
-CbmStsDigitizePixel::CbmStsDigitizePixel(Double_t pitchXcm, Double_t pitchYcm, Double_t pitchTns)
+CbmStsDigitizePixel::CbmStsDigitizePixel(Double_t pitchXcm, Double_t pitchYcm, Double_t resolutionTns)
   : CbmDigitize<CbmStsDigi>("StsDigitizePixel")
   , fPitchXcm(pitchXcm)
   , fPitchYcm(pitchYcm)
-  , fPitchTns(pitchTns)
+  , fResolutionTns(resolutionTns)
 {
 }
 // -------------------------------------------------------------------------
@@ -173,9 +175,10 @@ void CbmStsDigitizePixel::Exec(Option_t* /*opt*/)
     UInt_t address   = static_cast<UInt_t>(point->GetDetectorID());
     UShort_t channel = 0;
     double timef     = fCurrentEventTime + point->GetTime();
-    timef            = (floor(timef / fPitchTns) + 0.5) * fPitchTns;
-    Long64_t time    = Long64_t(round(timef));
-    UShort_t adc     = 30;
+    timef += gRandom->Gaus(0, fResolutionTns);
+    Long64_t time = Long64_t(round(timef));
+    if (time < 0) { time = 0; }
+    UShort_t adc = 30;
     assert(time >= 0);
 
     // Create digi and (if required) match and send them to DAQ
@@ -228,13 +231,13 @@ void CbmStsDigitizePixel::InitParams()
     UInt_t nAsicChannels = 128;   // Number of readout channels per ASIC
 
     // --- ASIC parameters
-    UShort_t nAdc      = 32;                     // Number of ADC channels (5 bit)
-    Double_t dynRange  = 75000.;                 // Dynamic range [e]
-    Double_t threshold = 3000.;                  // Threshold [e]
-    Double_t timeResol = fPitchTns / sqrt(12.);  // Time resolution [ns]
-    Double_t deadTime  = 800.;                   // Channel dead time [ns]
-    Double_t noiseRms  = 1000.;                  // RMS of noise [e]
-    Double_t znr       = 3.9789e-3;              // Zero-crossing noise rate [1/ns]
+    UShort_t nAdc      = 32;              // Number of ADC channels (5 bit)
+    Double_t dynRange  = 75000.;          // Dynamic range [e]
+    Double_t threshold = 3000.;           // Threshold [e]
+    Double_t timeResol = fResolutionTns;  // Time resolution [ns]
+    Double_t deadTime  = 800.;            // Channel dead time [ns]
+    Double_t noiseRms  = 1000.;           // RMS of noise [e]
+    Double_t znr       = 3.9789e-3;       // Zero-crossing noise rate [1/ns]
     CbmStsParAsic userParAsic(nAsicChannels, nAdc, dynRange, threshold, timeResol, deadTime, noiseRms, znr);
 
     CbmStsParModule userParModule(nChannels, nAsicChannels);

sim/detectors/sts/CbmStsDigitizePixel.h

@@ -54,7 +54,7 @@ class CbmStsDigitizePixel : public CbmDigitize<CbmStsDigi> {
 
 public:
   /** Constructor **/
-  CbmStsDigitizePixel(Double_t pitchXcm = 0.0010, Double_t pitchYcm = 0.0010, Double_t pitchTns = 20);
+  CbmStsDigitizePixel(Double_t pitchXcm = 0.0010, Double_t pitchYcm = 0.0010, Double_t resolutionTns = 5.);
 
 
   /** Destructor **/
@@ -92,8 +92,8 @@ public:
   /// set pitch Y [cm]
   void SetPitchY(double pitchYcm) { fPitchYcm = pitchYcm; }
 
-  /// set pitch Time [ns]
-  void SetPitchT(double pitchTns) { fPitchTns = pitchTns; }
+  /// set resolution Time [ns]
+  void SetResolutionTime(double resolutionTns) { fResolutionTns = resolutionTns; }
 
 private:
   Bool_t fIsInitialised;  ///< kTRUE if Init() was called
@@ -112,7 +112,7 @@ private:
 
   Double_t fPitchXcm {0.0058};  // pitch in X [cm]
   Double_t fPitchYcm {0.0058};  // pitch in Y [cm]
-  Double_t fPitchTns {20};      // pitch in time [ns]
+  Double_t fResolutionTns {5.};  // resolution in time [ns]
 
   /** @brief Initialise the parameters **/
   void InitParams();