xsec_vs_energy Class Reference

#include <neugen3_dummy.h>

List of all members.

Public Member Functions
	xsec_vs_energy ()
	xsec_vs_energy (int nentries, double E[], double xsec[])
	xsec_vs_energy (const xsec_vs_energy *xs)
	~xsec_vs_energy ()
void	load_from_file (const char *filename)
double	xsec (double E) const
TGraph *	get_as_graph (int npoints=100, bool scale_with_e=false) const
void	build_spline (int nentries, double E[], double xsec[])
double	E_min (void) const
double	E_max (void) const
Private Member Functions
void	init (void)
Private Attributes
TSpline3 *	_interpolator
pair< double, double >	_validity_range

---

Detailed Description

____________________________________________________________________________

NuValidator/generators/nugen

Costas Andreopoulos (Rutherford Lab.) <C.V.Andreopoulos@rl.ac.uk> Hugh Gallagher (Tufts University) <gallag@minos.phy.tufts.edu> _______________________________________________________________________________

Definition at line 30 of file neugen3_dummy.h.

---

Constructor & Destructor Documentation

xsec_vs_energy::xsec_vs_energy (   )

Definition at line 21 of file xsec_vs_energy.cxx. References init(). { init(); }

xsec_vs_energy::xsec_vs_energy (  int  nentries, double  E[], double  xsec[] )

Definition at line 26 of file xsec_vs_energy.cxx. References build_spline(), and init(). { init(); build_spline(nentries, E, xsec); }

xsec_vs_energy::xsec_vs_energy (  const xsec_vs_energy *  xs  )

Definition at line 33 of file xsec_vs_energy.cxx. References _interpolator, _validity_range, and xsec(). { _validity_range.first = xs->_validity_range.first; _validity_range.second = xs->_validity_range.second; // TSpline3 copy ctor was made protected (was dangerous anyway // since it did a shallow copy and thus shared underlying data) //_interpolator = new TSpline3( *(xs->_interpolator) ); // effectively do the copy ctor int nentries = xs->_interpolator->GetNpx(); double* E = new double[nentries]; double* xsec = new double[nentries]; for (int i=0; i<nentries; ++i) xs->_interpolator->GetKnot(i,E[i],xsec[i]); _interpolator = new TSpline3("spl3",E,xsec,nentries,"0", _validity_range.first, _validity_range.second); delete [] E; delete [] xsec; }

xsec_vs_energy::~xsec_vs_energy (   )

Definition at line 54 of file xsec_vs_energy.cxx. { if(_interpolator) delete _interpolator; }

---

Member Function Documentation

void xsec_vs_energy::build_spline (  int  nentries, double  E[], double  xsec[] )

Definition at line 74 of file xsec_vs_energy.cxx. References _interpolator, _validity_range, E_max(), and E_min(). Referenced by load_from_file(), and xsec_vs_energy(). { double E_min = E[ TMath::LocMin(nentries, E) ]; // minimum energy in spline double E_max = E[ TMath::LocMax(nentries, E) ]; // maximum energy in spline _validity_range.first = E_min; _validity_range.second = E_max; if(_interpolator) delete _interpolator; _interpolator = new TSpline3("spl3", E, xsec, nentries, "0", E_min, E_max); }

double xsec_vs_energy::E_max (  void   )  const [inline]

Definition at line 38 of file xsec_vs_energy.h. References _validity_range. Referenced by build_spline(), and get_as_graph(). { return _validity_range.second; }

double xsec_vs_energy::E_min (  void   )  const [inline]

Definition at line 37 of file xsec_vs_energy.h. References _validity_range. Referenced by build_spline(), and get_as_graph(). { return _validity_range.first; }

TGraph * xsec_vs_energy::get_as_graph (  int  npoints = 100, bool  scale_with_e = false )  const

Definition at line 96 of file xsec_vs_energy.cxx. References _validity_range, E_max(), E_min(), and xsec(). { double E_min = _validity_range.first; double E_max = _validity_range.second; if(npoints < 2) npoints = 2; double dE = (E_max - E_min) / (npoints-1); double * E = new double[npoints]; double * xs = new double[npoints]; for(int i=0; i<npoints; i++) { E[i] = E_min + i*dE; xs[i] = xsec( E[i] ); if (scale_with_e) xs[i] /= E[i]; } TGraph * graph = new TGraph(npoints, E, xs); //delete [] E; //delete [] xs; return graph; }

void xsec_vs_energy::init (  void   )  [private]

Definition at line 124 of file xsec_vs_energy.cxx. References _interpolator, and _validity_range. Referenced by xsec_vs_energy(). { _validity_range.first = 0; _validity_range.second = 0; _interpolator = 0; }

void xsec_vs_energy::load_from_file (  const char *  filename  )

Definition at line 59 of file xsec_vs_energy.cxx. References build_spline(), and xsec(). Referenced by nulook_wrapper::run(). { TNtuple nt("ntuple","xsec vs energy","energy:xsec"); nt.ReadFile(filename); nt.Draw("energy:xsec","","GOFF"); // select all entries - graphics off int nentries = nt.GetSelectedRows(); double * energy = nt.GetV1(); double * xsec = nt.GetV2(); build_spline(nentries, energy, xsec); }

double xsec_vs_energy::xsec (  double  E  )  const

Definition at line 87 of file xsec_vs_energy.cxx. References _interpolator, and _validity_range. Referenced by get_as_graph(), load_from_file(), and xsec_vs_energy(). { // check for E_min < E < E_max bool is_in_range = (_validity_range.first < E && E < _validity_range.second); if(is_in_range) return _interpolator->Eval(E); else return -1; }

---

Member Data Documentation

TSpline3* xsec_vs_energy::_interpolator [private]

Definition at line 44 of file xsec_vs_energy.h. Referenced by build_spline(), init(), xsec(), and xsec_vs_energy().

pair<double, double> xsec_vs_energy::_validity_range [private]

Definition at line 45 of file xsec_vs_energy.h. Referenced by build_spline(), E_max(), E_min(), get_as_graph(), init(), xsec(), and xsec_vs_energy().

---

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

• neugen3_dummy.h
• xsec_vs_energy.h
• xsec_vs_energy.cxx

---