---

GRBSim.h

#include "GRBShock.h"
#include "GRBConstants.h"
#include "GRBengine.h"
#include "SpectObj.h"
#include "GRBSynchrotron.h"
#include "GRBICompton.h"
#include <vector>

#ifndef GRBSIM_H
#define GRBSIM_H 1

class GRBSim 
{
 public:
  
  GRBSim(long seed=0); 
  GRBSim(const std::string& params);
  ~GRBSim();
  
  GRBConstants *myParam;
  
  void MakeGRB(double time_offset = 0.);
  
  std::vector<double> ComputeFlux(const double time);
  inline void setSpectrum(std::vector<double> myspectrum) 
    {m_spectrum=myspectrum;}
  inline std::pair<float,float> GRBdir()       {return m_direction;}
  inline std::vector<double>    Spectrum()     {return m_spectrum;}
  inline std::vector<double>    Energy()       {return m_energy;}
  inline std::vector<double>    DeltaE()       {return m_de;}
  
  inline double                 Tmax()         {return m_duration;}
  
  inline double                 Area()         {return m_area;}
  
  inline double                 Flux(int en)   {return m_spectrum[en];}
  
  inline double                 Energy(int en) {return m_energy[en];}
  
  double              IFlux( std::vector<double> spctrmVec,
                             double enmin=cst::enmin,
                             double enmax=cst::enmax);
  
  double              IRate(std::vector<double> spctrmVec,
                            double enmin=cst::enmin,
                            double enmax=cst::enmax);
  double                DrawPhotonFromSpectrum(std::vector<double> spctrmVec, 
                                               float u=0.0, 
                                               double enmin=cst::enmin);
  inline double EnergyPh(){return m_enph;}
  long parseParamList(std::string input, int index);
 private:
  //data member
  std::vector<GRBShell> theShells;
  std::vector<GRBShock> theShocks;
  std::vector<double>    m_energy, m_de, m_spectrum;
  std::vector<std::vector<double> > m_Fvt;
  std::pair<float,float> m_direction;
  
  double m_duration;
  double m_ftot;
  double m_phtot;
  double m_distance;
  double m_area;
  double m_DeadTime;
  double m_enph;
  double m_jetangle;
  long m_seed;
  HepRandomEngine *m_engine;
  GRBSynchrotron m_synchrotron;
  GRBICompton m_icompton;
  SpectObj m_spectobj;
};

    #endif

---