#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

◆ PtYDistributor() [1/2]

gen::PtYDistributor::PtYDistributor ( )

inline

Definition at line 16 of file PtYDistributor.h.

16 {};

◆ PtYDistributor() [2/2]

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.

     : 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

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

◆ ~PtYDistributor()

virtual gen::PtYDistributor::~PtYDistributor ( )

inlinevirtual

Definition at line 20 of file PtYDistributor.h.

20 {};

Member Function Documentation

◆ firePt() [1/2]

double PtYDistributor::firePt ( CLHEP::HepRandomEngine * engine )

Definition at line 72 of file PtYDistributor.cc.

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

References ptmax_, and ptmin_.

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

◆ firePt() [2/2]

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

Definition at line 87 of file PtYDistributor.cc.

                                                                                       {
   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;
 }

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

◆ fireY() [1/2]

double PtYDistributor::fireY ( CLHEP::HepRandomEngine * engine )

Definition at line 70 of file PtYDistributor.cc.

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

References ymax_, and ymin_.

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

◆ fireY() [2/2]

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

Definition at line 74 of file PtYDistributor.cc.

                                                                                    {
   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;
 }