#include <PtYDistributor.h>

Public Member Functions
double	firePt ()

double	firePt (double ptmin, double ptmax)

double	fireY ()

double	fireY (double ymin, double ymax)

	PtYDistributor ()

	PtYDistributor (edm::FileInPath fip, CLHEP::HepRandomEngine &fRandomEngine, 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 15 of file PtYDistributor.h.

Constructor & Destructor Documentation

gen::PtYDistributor::PtYDistributor ( )

inline

Definition at line 17 of file PtYDistributor.h.

17 {};

PtYDistributor::PtYDistributor	(	edm::FileInPath	fip,
		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 gather_cfg::cout, edm::hlt::Exception, f, fPtGenerator, edm::FileInPath::fullPath(), fYGenerator, i, edm::errors::NullPointerError, ptbins_, ptmax_, ptmin_, create_public_lumi_plots::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(fRandomEngine,aProbFunc1,ybins_);
   fPtGenerator = new CLHEP::RandGeneral(fRandomEngine,aProbFunc2,ptbins_);
 
   f.Close();
   
 } // from file

virtual gen::PtYDistributor::~PtYDistributor ( )

inlinevirtual

Definition at line 22 of file PtYDistributor.h.

22 {};

Member Function Documentation

double PtYDistributor::firePt ( )

Definition at line 68 of file PtYDistributor.cc.

References ptmax_, and ptmin_.

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

                       {
    return firePt(ptmin_,ptmax_);
 }

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

Definition at line 88 of file PtYDistributor.cc.

References edm::hlt::Exception, 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::fireY ( )

Definition at line 63 of file PtYDistributor.cc.

References ymax_, and ymin_.

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

                      {
    return fireY(ymin_,ymax_);
 }

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

Definition at line 73 of file PtYDistributor.cc.

References edm::hlt::Exception, fYGenerator, edm::errors::NullPointerError, detailsBasic3DVector::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;
 }