GretaPolynomialReader< Log, Exp >

Description

template<double(*)(double) Log = log, double(*)(double) Exp = exp>
class DAS::GretaPolynomialReader< Log, Exp >

Wrapper to read series of 1D polynomials from a file

Expected format:

• one column showing the low edge of the rapidity bin (valid until another is defined)
• two columns showing the pt range
• one column showing the degree (= #params - 1)
• arbitrary number of columns containing the parameters (corresponding to the degree)

#include <Greta.h>

Collaboration diagram for GretaPolynomialReader< Log, Exp >:

Public Member Functions
void	read (TString fname)
double	operator() (double pt, double y)

Private Attributes
vector< double >	ymins
vector< GretaPolynomial< Log, Exp > >	pols

Member Function Documentation

operator()()

double operator() ( double  pt, double  y  )

inline

    {
        y = abs(y);
        size_t iy = ymins.size();
        do --iy;
        while (ymins.at(iy) > y);
        auto& g = pols.at(iy);
        return g(pt);
    }

read()

void read ( TString  fname )

inline

    {
        cout << "Reading " << fname << endl;
        pols.clear();
        ymins.clear();
        ifstream infile;
        infile.open(fname.Data());
        if (infile.fail()) {
            cerr << "Failure at opening the file " << fname << '\n';
            exit(EXIT_FAILURE);
        }
 
        while (infile.good()) { 
 
            // getting values
            double ymin, ptmin, ptmax;
            int Npars;
            infile >> ymin >> ptmin >> ptmax >> Npars;
            double * pars = new double[Npars];
            for (int i = 0; i < Npars; ++i) infile >> pars[i];
            
            if (infile.eof()) break;
 
            // instantiating
            if (ymins.size() == 0) { assert(abs(ymin) < deps);    }
            else                   { assert(ymin > ymins.back()); }
            ymins.push_back(ymin);
            GretaPolynomial<Log,Exp> gp(pars, Npars, ptmin, ptmax);
            pols.push_back(gp);
            //pols.insert( { ymin, gp } );
 
            // print-out
            cout << std::setw(3) << ymin
                 << std::setw(8) << (int) ptmin
                 << std::setw(8) << (int) ptmax
                 << std::setw(3) << Npars;
            for (int j = 0; j < Npars; ++j) std::cout << std::setw(15) << pars[j];
            cout << '\n';
        } 
        cout << flush;
 
        assert(ymins.size() > 0);
        assert(ymins.size() == pols.size());
 
        infile.close(); 
    }

Member Data Documentation

pols

vector<GretaPolynomial<Log,Exp> > pols

private

first index for low edge of rap bin

ymins

vector<double> ymins

private

---

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

• interface/Greta.h