---

GRBpulse.h

#ifndef GRB_PULSE_H
#define GRB_PULSE_H

#include <vector>

class HepRandomEngine;



class GRBpulse
{
public:
    
    GRBpulse();
    
    ~GRBpulse();

    

    // Accessors
    const std::vector<double> &tmax()  const            { return m_tmax; }
    const std::vector< std::vector<double> > &pulse()   { return m_pulse; }
    const std::vector<long>   &nphotpul() const         { return m_nphotpul; }
    double univFWHM() const                             { return m_univFWHM; }
    


    long data(HepRandomEngine *engine, const double ethres, const long nphoton, 
        const int npuls, const double duration);
    
private:
    
    // Accessors
    const std::vector<double> &amplitude() const        { return m_amplitude; }
    const std::vector<double> &tdiff() const            { return m_tdiff; }
    const std::vector<double> &sigma()  const           { return m_sigma; }

    
    long createSigmaTdiff(HepRandomEngine *engine);
    
    void fillVector(std::vector<long> &ngtwid) const;
    
    double getAlpha(const double value) const;
    
    void getAmplitude(HepRandomEngine *engine, const int npuls);
    
    void getNphotpul(const long nphoton, const int npuls);
    
    void getPulse(const int npuls);
    
    void getTmax(HepRandomEngine *engine, const int npuls, const double duration);
    
    long index(HepRandomEngine *engine, const long diff, const long minval, 
        const std::vector<long> &in) const;
    
    void pickWidth(HepRandomEngine *engine, const double ethres, 
        const double duration);
    
    void universalWidth(HepRandomEngine *engine, const double ethres, 
        const double duration, const long diff, const long minval, 
        const std::vector<long> &in, const std::vector<double> &v);
    
    
    // Data Members
    std::vector<double>                 m_amplitude;
    std::vector<double>                 m_tdiff;
    std::vector<double>                 m_tmax;
    std::vector<double>                 m_sigma;
    std::vector< std::vector<double> >  m_pulse;
    std::vector<long>                   m_nphotpul;
    double                              m_univFWHM;
};

#endif // GRB_PULSE_H

---