PPSFastLocalSimulation::Distribution Struct Reference

Public Types
enum	Type { dtBox, dtGauss, dtGaussLimit }

Public Member Functions
	Distribution (const edm::ParameterSet &)
void	Generate (CLHEP::HepRandomEngine &rndEng, double &x, double &y)

Public Attributes
enum PPSFastLocalSimulation::Distribution::Type	type_
double	x_max_
double	x_mean_
double	x_min_
double	x_width_
double	y_max_
double	y_mean_
double	y_min_
double	y_width_

Detailed Description

Definition at line 80 of file PPSFastLocalSimulation.cc.

Member Enumeration Documentation

Type

enum PPSFastLocalSimulation::Distribution::Type

Enumerator
dtBox
dtGauss
dtGaussLimit

Definition at line 81 of file PPSFastLocalSimulation.cc.

{ dtBox, dtGauss, dtGaussLimit } type_;

Constructor & Destructor Documentation

Distribution()

PPSFastLocalSimulation::Distribution::Distribution

(

const edm::ParameterSet &

ps

)

Definition at line 125 of file PPSFastLocalSimulation.cc.

References Exception, edm::ParameterSet::getParameter(), and dqmiodumpindices::typeName.

                                                                          {
  // get type
  string typeName = ps.getParameter<string>("type");
  if (!typeName.compare("box"))
    type_ = dtBox;
  else if (!typeName.compare("gauss"))
    type_ = dtGauss;
  else if (!typeName.compare("gauss-limit"))
    type_ = dtGaussLimit;
  else
    throw cms::Exception("PPS") << "Unknown distribution type `" << typeName << "'.";

  x_mean_ = ps.getParameter<double>("x_mean");
  x_width_ = ps.getParameter<double>("x_width");
  x_min_ = ps.getParameter<double>("x_min");
  x_max_ = ps.getParameter<double>("x_max");

  y_mean_ = ps.getParameter<double>("y_mean");
  y_width_ = ps.getParameter<double>("y_width");
  y_min_ = ps.getParameter<double>("y_min");
  y_max_ = ps.getParameter<double>("y_max");
}

Member Function Documentation

Generate()

void PPSFastLocalSimulation::Distribution::Generate	(	CLHEP::HepRandomEngine &	rndEng,
		double &	x,
		double &	y
	)

Definition at line 150 of file PPSFastLocalSimulation.cc.

References mathSSE::sqrt(), and x.

Referenced by PPSFastLocalSimulation::GenerateTrack().

                                                                                                    {
  switch (type_) {
    case dtBox:
      x = x_mean_ + x_width_ * (rndEng.flat() - 0.5);
      y = y_mean_ + y_width_ * (rndEng.flat() - 0.5);
      break;

    case dtGauss:
      x = x_mean_ + RandGauss::shoot(&rndEng) * x_width_;
      y = y_mean_ + RandGauss::shoot(&rndEng) * y_width_;
      break;

    case dtGaussLimit: {
      const double u_x = rndEng.flat(), u_y = rndEng.flat();

      const double cdf_x_min = (1. + TMath::Erf((x_min_ - x_mean_) / x_width_ / sqrt(2.))) / 2.;
      const double cdf_x_max = (1. + TMath::Erf((x_max_ - x_mean_) / x_width_ / sqrt(2.))) / 2.;
      const double a_x = cdf_x_max - cdf_x_min, b_x = cdf_x_min;

      const double cdf_y_min = (1. + TMath::Erf((y_min_ - y_mean_) / y_width_ / sqrt(2.))) / 2.;
      const double cdf_y_max = (1. + TMath::Erf((y_max_ - y_mean_) / y_width_ / sqrt(2.))) / 2.;
      const double a_y = cdf_y_max - cdf_y_min, b_y = cdf_y_min;

      x = x_mean_ + x_width_ * sqrt(2.) * TMath::ErfInverse(2. * (a_x * u_x + b_x) - 1.);
      y = y_mean_ + y_width_ * sqrt(2.) * TMath::ErfInverse(2. * (a_y * u_y + b_y) - 1.);
    }

    break;

    default:
      x = y = 0.;
  }
}

Member Data Documentation

type_

enum PPSFastLocalSimulation::Distribution::Type PPSFastLocalSimulation::Distribution::type_

Referenced by Modules.Service::__str__(), Modules.Service::_placeImpl(), Modules.ESSource::_placeImpl(), Modules.ESProducer::_placeImpl(), Modules.ESPrefer::copy(), Mixins._TypedParameterizable::dumpPython(), Modules.ESPrefer::dumpPythonAs(), Modules.Service::dumpSequencePython(), RandomRunSource.RandomRunSource::insertInto(), Modules.Service::insertInto(), ExternalGeneratorFilter.ExternalGeneratorFilter::insertInto(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::moduleLabel_(), Modules.ESProducer::moduleLabel_(), Modules.ESSource::nameInProcessDesc_(), Modules.ESProducer::nameInProcessDesc_(), and Modules.ESPrefer::nameInProcessDesc_().

x_max_

double PPSFastLocalSimulation::Distribution::x_max_

Definition at line 82 of file PPSFastLocalSimulation.cc.

x_mean_

double PPSFastLocalSimulation::Distribution::x_mean_

Definition at line 82 of file PPSFastLocalSimulation.cc.

x_min_

double PPSFastLocalSimulation::Distribution::x_min_

Definition at line 82 of file PPSFastLocalSimulation.cc.

x_width_

double PPSFastLocalSimulation::Distribution::x_width_

Definition at line 82 of file PPSFastLocalSimulation.cc.

y_max_

double PPSFastLocalSimulation::Distribution::y_max_

Definition at line 83 of file PPSFastLocalSimulation.cc.

y_mean_

double PPSFastLocalSimulation::Distribution::y_mean_

Definition at line 83 of file PPSFastLocalSimulation.cc.

y_min_

double PPSFastLocalSimulation::Distribution::y_min_

Definition at line 83 of file PPSFastLocalSimulation.cc.

y_width_

double PPSFastLocalSimulation::Distribution::y_width_

Definition at line 83 of file PPSFastLocalSimulation.cc.