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