---

GRBConstants.h

#ifndef GRBCONSTANTS_HH
#define GRBCONSTANTS_HH 1
#include <string>
#include <cmath>
#include "CLHEP/Random/RandomEngine.h"

namespace cst
{
  // Universal constants :
  const double mpc2      = 938.2; 
  const double mec2      = 0.510999;       //MeV
  const double erg2MeV   = 624151.0;
  const double mpc2cm    = 3.0857e+24;
  // Conversion fron \e Gauss to \e MeV
  //  const double G2MeV     = 815.78;
  
  const double st        = 6.65225e-25;
  const double mu0       = 4.0*M_PI*1.0e-7; 
  const double   BQ        = 4.413e+13; 
  const double c         = 2.98e+10;
  const double c2        = c*c;
  const double hplanck   = 4.13567e-15;
  const double pi        = 3.1415926535897932385;
  const double Hubble    = 6.5e+1; 
  const double wzel      = 0.0;
  const double csi       = 1.0; 
  const double alphae    = .5;
  const double alphab    = .1;
  const double p         = 2.5;
  const double viscosity = 0.0;
  const int    nstep        = 300; 
  const int    enstep       =  50;
  const double enmin     = 1.0e+3;
  const double enmax     =1.0e+12;
  // Minimal Temporal separation between 2 photons
  // const double DeadTime  = 1.0e-5; //sec
  
  const float flagIC     = 1;
  const bool  flagQG     = false;
  const bool savef=false;
  const std::string paramFile= "GRBdata.txt";
  const std::string pulse_shape="agauss"; 
}

class GRBConstants 
{ 
 public:

  GRBConstants();

  ~GRBConstants() { }

  static void InitializeRandom(long seed = 0);
  
  static HepRandomEngine* GetTheRandomEngine(long seed = 0);
  
  int ReadParam();
  
  void Print();
  
  void Save(bool flag=false);


  inline void setEngineType(int value){m_engineType = value;}
  
  inline int EngineType(){return m_engineType;}
    

  //  inline int Nshell() {return m_nshell;}

  //  inline void setNshell(int value=10){m_nshell = value;}
    

  inline int Nshock() {return m_nshock;}

  inline void setNshock(int value=10){m_nshock = value;}

  inline double Duration(){return m_duration;}

  inline void setDuration(double value=0.4){m_duration = value;}

  inline double RiseTime(){return m_rt;}

  inline void setRiseTime(double value=0.4){m_rt = value;}
 
  inline double DecayTime(){return m_dt;}
 
  inline void setDecayTime(double value=0.4){m_dt = value;}
 
  inline double PeakEnergy(){return m_peak;}
 
  inline void setPeakEnergy(double value=0.4){m_peak = value;}
  inline void setShellType(int value){m_shellType = value;}
  
  inline int ShellType(){return m_shellType;}
  

  inline double ShellRadius(){return m_d0;}
  
  inline void setShellRadius(double value){m_d0 = value;}

 
  inline double JetRadius(){return m_r0;}
  
  inline void setJetRadius(double value){m_r0 = value;}
  
  inline double JetAngle(){return m_angle;}
  
  inline void setJetAngle(double value){m_angle = M_PI/180. * value;}
  
 
  inline double Thickness(){return m_t0;}

  inline void setThickness(double value){m_t0 = value;}
   
  inline double Etot(){return m_etot;}

  inline void setEtot(double value){m_etot = value;}

  inline double GammaMin(){return m_g0;}

  inline void setGammaMin(double value=100.0){m_g0 = value;}

  inline double GammaMax(){return m_g1;}

  inline void setGammaMax(double value=1000.0){m_g1 = value;}

  inline double Redshift() {return m_redshift;}

  inline void setRedshift(double value=1.0){m_redshift = value;}
  inline void setEnergyPh(double value=25.0e+3){m_enph = value;}
  
  inline double EnergyPh(){return m_enph;}
  
  double Distance();
  
  inline std::string GetParamFile() {return  m_paramFile;}
 
 private:
  char* m_burst_type;
  //  int m_nshell;
  double m_d0;
  int m_nshock;
  double m_r0;
  double m_angle;
  double m_redshift;
  double m_etot;
  double m_g0,m_g1;
  double m_t0;
  double m_duration;
  double m_rt,m_dt,m_peak;
  double m_enph;
  std::string m_paramFile;
  int m_engineType;
  int m_shellType;
};

#endif

---