LandauFluctuationGenerator.h

Go to the documentation of this file.

#ifndef FastSimulation_Utilities_LandauFluctuationGenerator_H
#define FastSimulation_Utilities_LandauFluctuationGenerator_H

#include "FastSimulation/Utilities/interface/BaseNumericalRandomGenerator.h"

#include <cmath>

class RandomEngineAndDistribution;

class LandauFluctuationGenerator : public BaseNumericalRandomGenerator {
public:
  LandauFluctuationGenerator() : BaseNumericalRandomGenerator(-3.5, 25.) { initialize(); }

  ~LandauFluctuationGenerator() override {}

  double landau(RandomEngineAndDistribution const* random) const { return generate(random); }

  double function(double x) override { return ersatzt(x); }

private:
  double ersatzt(double x) { return std::exp(-0.5 * (x + std::exp(-x))) / std::sqrt(2. * M_PI); }
};
#endif