LogXGaussSmearing.h

Go to the documentation of this file.

#ifndef _CachedXShowerRegeneratorAG_LogXGaussSmearing_h_
#define _CachedXShowerRegeneratorAG_LogXGaussSmearing_h_

#include <CLHEP/Random/RandGauss.h>

#include <utl/Point.h>
#include <utl/Vector.h>
#include <utl/CoordinateSystemPtr.h>
#include <utl/TimeStamp.h>
#include <utl/TimeInterval.h>
#include <utl/RandomEngine.h>


namespace CachedXShowerRegeneratorAG {

  class LogXGaussSmearing {
  public:
    LogXGaussSmearing(const double width, utl::RandomEngine::RandomEngineType* const engine)
      : fWidth(width), fRandomEngine(engine) { }

    double
    GetSmearedTime(const double planeFrontTime, const double particleTime)
      const
    {
      if (particleTime > planeFrontTime) {
        const double fact = exp(fWidth * CLHEP::RandGauss::shoot(fRandomEngine));
        return planeFrontTime * (1 - fact) + particleTime * fact;
      }
      return particleTime;
    }

  private:
    const double fWidth;
    utl::RandomEngine::RandomEngineType* const fRandomEngine;
  };

}


#endif