TMDlib_UnintegratedGluonKS.h

Go to the documentation of this file.

//------------------------------------------------------------------------------
//
// Copyright (c) 2012, Sebastian Sapeta
//
//  This is part of KSgluon package.
//
//  KSgluon is free software; you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation; either version 2 of the License, or
//  (at your option) any later version.
//
//  KSgluon is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  You should have received a copy of the GNU General Public License
//  along with FastJet; if not, write to the Free Software
//  Foundation, Inc.:
//      59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
//------------------------------------------------------------------------------
#ifndef __TMDlib_UNINTEGRATEDGLUONKS_H__
#define __TMDlib_UNINTEGRATEDGLUONKS_H__

#include <string>
#include <vector>
#include "TMDlib_InterpolationKS.h"
using namespace std;

//------------------------------------------------------------------------------
enum UnintegratedGluonName {
  grid2d,  // gluon interpolated from 2d grid
  grid3d   // gluon interpolated from 3d grid
};

//------------------------------------------------------------------------------
class UnintegratedGluon {
  public:
  UnintegratedGluon ();
  UnintegratedGluon (string filename, gsl_interp_type int_type);
  //pair<unsigned int,unsigned int> grid_size() const;
  vector<int> grid_size() const;
  vector<pair <double, double> > grid_limits() const {
    return _updf->grid_limits();}
  double xg(double logx, double logkt2, double logmu2 = 0.0);
  virtual ~UnintegratedGluon ();
  private:
  // function determining the size of the 2d grid
  void _get_2dgrid_dimensions(string filename);
  // function determining the size of the 3d grid
  void _get_3dgrid_dimensions(string filename);
  void _set_from_2dgrid(string filename);
  void _set_from_3dgrid(string filename);
  // functions setting gbw pdf
  void _set_gbw();
  UnintegratedGluonName _determine_grid_type(string filename);
  unsigned int _n1, _n2, _n3;
  ofstream *_errorfile;
  static double _norm_F(const vector<double> &x);
  BaseInterpolation *_updf;
  UnintegratedGluonName _name;
  // pointer to the normalization function for pdf
  double (*_norm)(const vector<double> &x);
  gsl_interp_type _int_type;
  double _some_eps;
};


#endif // __TMDlib_UNINTEGRATEDGLUONKS_H__