---

GRBobsUtilities.h

#ifndef GRB_OBS_UTILITIES_H
#define GRB_OBS_UTILITIES_H

#include <vector>
#include <deque>
#include <string>
#include <functional>  // for unary_function and plus

#include "GRBobsConstants.h"

using namespace grbobstypes;

class HepRandomEngine;



class GRBobsUtilities
{
public:
    static double result(HepRandomEngine *engine, const double lo, 
        const double hi, const double p);
    
    static void sortVector(const long index, const std::vector<double> &in, 
        const std::vector<double> &sorted, std::vector<double> &out);
     
    struct multiplier : public std::unary_function<double, double>
    {
        multiplier(double value) : m_value(value) {}
        double operator() (double x);
        
        double m_value;
    };
    
    
    struct randGen : public std::unary_function<double, double>
    {
        randGen(HepRandomEngine *engine) : m_engine(engine) {}
        double operator() (double x); 
        
        HepRandomEngine *m_engine;
    };
    
    
    template<class T>
        static void getSorter(const std::vector<T> &in, 
        std::vector<long> &sorter)
    {
        std::vector<T> sorted(in);
        std::sort(sorted.begin(), sorted.end());
        
        LongSize sz = in.size();
        
        std::deque<bool> qused(sz, 0);
        
        for (LongSize i=0; i<sz; ++i)
        {
            for (LongSize j=0; j<sz; ++j)
            {
                if ((!qused[j]) && (sorted[i] == in[j]))
                {
                    sorter[i] = j;
                    qused[j] = 1;
                    break;
                }
            }
        }
    }
    
    
    template<class S, class T>
        static void cumulativeSum(const std::vector<S> &in, std::vector<T> &out)
    {
        out.clear();
        
        out.reserve(in.size());
        
        std::vector<S>::const_iterator it = in.begin();
        T value = *it;
        out.push_back(value);
        
        ++it;
        
        while (it != in.end())
        {
            value += *it;
            out.push_back(value);
            ++it;
        }
    }
};

#endif

---