gen::PtYDistributor Class Reference

#include <PtYDistributor.h>

Public Member Functions
double	firePt (CLHEP::HepRandomEngine *)
double	firePt (double ptmin, double ptmax, CLHEP::HepRandomEngine *)
double	fireY (CLHEP::HepRandomEngine *)
double	fireY (double ymin, double ymax, CLHEP::HepRandomEngine *)
	PtYDistributor ()
	PtYDistributor (const edm::FileInPath &fip, double ptmax, double ptmin, double ymax, double ymin, int ptbins, int ybins)
virtual	~PtYDistributor ()

Private Attributes
CLHEP::RandGeneral *	fPtGenerator
CLHEP::RandGeneral *	fYGenerator
int	ptbins_
double	ptmax_
double	ptmin_
int	ybins_
double	ymax_
double	ymin_

Detailed Description

Definition at line 14 of file PtYDistributor.h.

Constructor & Destructor Documentation

gen::PtYDistributor::PtYDistributor ( )

inline

Definition at line 16 of file PtYDistributor.h.

{};

PtYDistributor::PtYDistributor	(	const edm::FileInPath &	fip,
		double	ptmax = 100,
		double	ptmin = 0,
		double	ymax = 10,
		double	ymin = -10,
		int	ptbins = 1000,
		int	ybins = 50
	)

Definition at line 16 of file PtYDistributor.cc.

References gather_cfg::cout, Exception, validate-o2o-wbm::f, fPtGenerator, edm::FileInPath::fullPath(), fYGenerator, mps_fire::i, edm::errors::NullPointerError, ptbins_, ptmax_, ptmin_, AlCaHLTBitMon_QueryRunRegistry::string, xy(), ybins_, ymax_, and ymin_.

    : ptmax_(ptmax), ptmin_(ptmin), ymax_(ymax), ymin_(ymin), ptbins_(ptbins), ybins_(ybins) {
  // edm::FileInPath f1(input);
  std::string fDataFile = fip.fullPath();

  std::cout << " File from " << fDataFile << std::endl;
  TFile f(fDataFile.c_str(), "READ");
  TGraph* yfunc = (TGraph*)f.Get("rapidity");
  TGraph* ptfunc = (TGraph*)f.Get("pt");

  if (!yfunc)
    throw edm::Exception(edm::errors::NullPointerError, "PtYDistributor")
        << "Rapidity distribution could not be found in file " << fDataFile;

  if (!ptfunc)
    throw edm::Exception(edm::errors::NullPointerError, "PtYDistributor")
        << "Pt distribution could not be found in file " << fDataFile;

  if (ptbins_ > 100000) {
    ptbins_ = 100000;
  }

  if (ybins_ > 100000) {
    ybins_ = 100000;
  }

  double aProbFunc1[100000];
  double aProbFunc2[100000];

  for (int i = 0; i < ybins_; ++i) {
    double xy = ymin_ + i * (ymax_ - ymin_) / ybins_;
    double yy = yfunc->Eval(xy);
    aProbFunc1[i] = yy;
  }

  for (int ip = 0; ip < ptbins_; ++ip) {
    double xpt = ptmin_ + ip * (ptmax_ - ptmin_) / ptbins_;
    double ypt = ptfunc->Eval(xpt);
    aProbFunc2[ip] = ypt;
  }

  fYGenerator = new CLHEP::RandGeneral(nullptr, aProbFunc1, ybins_);
  fPtGenerator = new CLHEP::RandGeneral(nullptr, aProbFunc2, ptbins_);

  f.Close();

}  // from file

virtual gen::PtYDistributor::~PtYDistributor ( )

inlinevirtual

Definition at line 20 of file PtYDistributor.h.

{};

Member Function Documentation

double PtYDistributor::firePt

(

CLHEP::HepRandomEngine *

engine

)

Definition at line 72 of file PtYDistributor.cc.

References ptmax_, and ptmin_.

Referenced by gen::Pythia6PtYDistGun::generateEvent().

{ return firePt(ptmin_, ptmax_, engine); }

double PtYDistributor::firePt	(	double	ptmin,
		double	ptmax,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 87 of file PtYDistributor.cc.

References Exception, fPtGenerator, edm::errors::NullPointerError, DiDispStaMuonMonitor_cfi::pt, ptmax_, and ptmin_.

                                                                                      {
  double pt = -999;

  if (fPtGenerator) {
    while (pt < ptmin || pt > ptmax)
      pt = ptmin_ + (ptmax_ - ptmin_) * fPtGenerator->shoot(engine);
  } else {
    throw edm::Exception(edm::errors::NullPointerError, "PtYDistributor")
        << "Random pt requested but Random Number Generator for pt not Initialized!";
  }
  return pt;
}

double PtYDistributor::fireY

(

CLHEP::HepRandomEngine *

engine

)

Definition at line 70 of file PtYDistributor.cc.

References ymax_, and ymin_.

Referenced by gen::Pythia6PtYDistGun::generateEvent().

{ return fireY(ymin_, ymax_, engine); }

double PtYDistributor::fireY	(	double	ymin,
		double	ymax,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 74 of file PtYDistributor.cc.

References Exception, fYGenerator, edm::errors::NullPointerError, ymax_, and ymin_.

                                                                                   {
  double y = -999;

  if (fYGenerator) {
    while (y < ymin || y > ymax)
      y = ymin_ + (ymax_ - ymin_) * fYGenerator->shoot(engine);
  } else {
    throw edm::Exception(edm::errors::NullPointerError, "PtYDistributor")
        << "Random y requested but Random Number Generator for y not Initialized!";
  }
  return y;
}

Member Data Documentation

CLHEP::RandGeneral* gen::PtYDistributor::fPtGenerator

private

Definition at line 37 of file PtYDistributor.h.

Referenced by firePt(), and PtYDistributor().

CLHEP::RandGeneral* gen::PtYDistributor::fYGenerator

private

Definition at line 36 of file PtYDistributor.h.

Referenced by fireY(), and PtYDistributor().

int gen::PtYDistributor::ptbins_

private

Definition at line 33 of file PtYDistributor.h.

Referenced by PtYDistributor().

double gen::PtYDistributor::ptmax_

private

Definition at line 28 of file PtYDistributor.h.

Referenced by firePt(), and PtYDistributor().

double gen::PtYDistributor::ptmin_

private

Definition at line 29 of file PtYDistributor.h.

Referenced by firePt(), and PtYDistributor().

int gen::PtYDistributor::ybins_

private

Definition at line 34 of file PtYDistributor.h.

Referenced by PtYDistributor().

double gen::PtYDistributor::ymax_

private

Definition at line 30 of file PtYDistributor.h.

Referenced by fireY(), and PtYDistributor().

double gen::PtYDistributor::ymin_

private

Definition at line 31 of file PtYDistributor.h.

Referenced by fireY(), and PtYDistributor().