---

edm::PtYDistributor Class Reference

#include <GeneratorInterface/Pythia6Interface/interface/PtYDistributor.h>

List of all members.

Public Member Functions
double	firePt (double ptmin, double ptmax)
double	firePt ()
double	fireY (double ymin, double ymax)
double	fireY ()
	PtYDistributor (std::string inputfile, CLHEP::HepRandomEngine &fRandomEngine, double ptmax, double ptmin, double ymax, double ymin, int ptbins, int ybins)
	PtYDistributor ()
virtual	~PtYDistributor ()
Private Attributes
std::string	file
RandGeneral *	fPtGenerator
RandGeneral *	fYGenerator
int	ptbins_
double	ptmax_
double	ptmin_
int	ybins_
double	ymax_
double	ymin_

---

Detailed Description

Definition at line 13 of file PtYDistributor.h.

---

Constructor & Destructor Documentation

edm::PtYDistributor::PtYDistributor

(

)

[inline]

Definition at line 15 of file PtYDistributor.h.

{};

PtYDistributor::PtYDistributor	(	std::string	inputfile,
		CLHEP::HepRandomEngine &	fRandomEngine,
		double	ptmax = 100,
		double	ptmin = 0,
		double	ymax = 10,
		double	ymin = -10,
		int	ptbins = 1000,
		int	ybins = 50
	)

Definition at line 10 of file PtYDistributor.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), f, f1, fPtGenerator, edm::FileInPath::fullPath(), fYGenerator, i, edm::errors::NullPointerError, ptbins_, ptmax_, ptmin_, ybins_, ymax_, and ymin_.

                                                                                                                                                                                                 : ptmax_(ptmax),ptmin_(ptmin),ymax_(ymax),ymin_(ymin), ptbins_(ptbins), ybins_(ybins)
{  
   edm::FileInPath f1(input);
   std::string fDataFile = f1.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 RandGeneral(fRandomEngine,aProbFunc1,ybins_);
  fPtGenerator = new RandGeneral(fRandomEngine,aProbFunc2,ptbins_);

  f.Close();
  
} // from file

virtual edm::PtYDistributor::~PtYDistributor

(

)

[inline, virtual]

Definition at line 17 of file PtYDistributor.h.

{};

---

Member Function Documentation

double PtYDistributor::firePt	(	double	ptmin,
		double	ptmax
	)

Definition at line 84 of file PtYDistributor.cc.

References fPtGenerator, edm::errors::NullPointerError, ptmax_, and ptmin_.

                                                {

   double pt = -999;

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

double PtYDistributor::firePt

(

)

Definition at line 64 of file PtYDistributor.cc.

References ptmax_, and ptmin_.

Referenced by edm::PythiaSource::produce(), and edm::PythiaProducer::produce().

                      {
   return firePt(ptmin_,ptmax_);
}

double PtYDistributor::fireY	(	double	ymin,
		double	ymax
	)

Definition at line 69 of file PtYDistributor.cc.

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

                                             {

   double y = -999;

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

double PtYDistributor::fireY

(

)

Definition at line 59 of file PtYDistributor.cc.

References ymax_, and ymin_.

Referenced by edm::PythiaSource::produce(), and edm::PythiaProducer::produce().

                     {
   return fireY(ymin_,ymax_);
}

---

Member Data Documentation

std::string edm::PtYDistributor::file [private]

Definition at line 33 of file PtYDistributor.h.

RandGeneral* edm::PtYDistributor::fPtGenerator [private]

Definition at line 36 of file PtYDistributor.h.

Referenced by firePt(), and PtYDistributor().

RandGeneral* edm::PtYDistributor::fYGenerator [private]

Definition at line 35 of file PtYDistributor.h.

Referenced by fireY(), and PtYDistributor().

int edm::PtYDistributor::ptbins_ [private]

Definition at line 30 of file PtYDistributor.h.

Referenced by PtYDistributor().

double edm::PtYDistributor::ptmax_ [private]

Definition at line 25 of file PtYDistributor.h.

Referenced by firePt(), and PtYDistributor().

double edm::PtYDistributor::ptmin_ [private]

Definition at line 26 of file PtYDistributor.h.

Referenced by firePt(), and PtYDistributor().

int edm::PtYDistributor::ybins_ [private]

Definition at line 31 of file PtYDistributor.h.

Referenced by PtYDistributor().

double edm::PtYDistributor::ymax_ [private]

Definition at line 27 of file PtYDistributor.h.

Referenced by fireY(), and PtYDistributor().

double edm::PtYDistributor::ymin_ [private]

Definition at line 28 of file PtYDistributor.h.

Referenced by fireY(), and PtYDistributor().

---

The documentation for this class was generated from the following files:

• GeneratorInterface/Pythia6Interface/interface/PtYDistributor.h
• GeneratorInterface/Pythia6Interface/src/PtYDistributor.cc

---