VertexGen_Spectrum.hh

// RAT::VertexGen_Spectrum
// 18-July-2008 JR Wilson
 
#ifndef __RAT_VertexGen_Spectrum__
#define __RAT_VertexGen_Spectrum__
 
#include <G4Event.hh>
#include <G4ThreeVector.hh>
#include <RAT/GLG4VertexGen.hh>
#include <globals.hh>
#include <vector>
 
#include "RAT/DB.hh"
 
namespace RAT {
 
class VertexGen_Spectrum : public GLG4VertexGen {
 public:
  VertexGen_Spectrum(const char *arg_dbname = "SPECTRUM");
  virtual ~VertexGen_Spectrum();
  virtual void GeneratePrimaryVertex(G4Event *argEvent, G4ThreeVector &dx, G4double dt);
  virtual void SetState(G4String newValues);
  virtual G4String GetState();
 
  virtual void InitialiseSpectrum();
 
  virtual double SampleEnergy();
 
  virtual bool ELimitable() { return true; };
 
  virtual void LimitEnergies(double Elo, double Ehi);
 
  virtual double EMaximum();
 
  virtual double EMinimum();
 
 private:
  G4String _particle;           // name of the particle type
  G4ParticleDefinition *_pDef;  // particle definition
  G4String _spectrum;           // name of the spectrum to use
  DBLinkPtr _lspec;             // link to spectrum entry in database
  double _emin;                 // valid range for spectrum in MeV
  double _emax;
  std::vector<double> spec_E;       // spectrum energy values
  std::vector<double> spec_mag;     // spectrum magnitude values
  std::vector<double> spec_cummag;  // spectrum cumulative non-normalised magnitude values
  double Elim_Ulo;                  // user applied universal lower limit to energy range
  double Elim_Uhi;                  // user applied universal higher limit to energy range
  double Elim_Tlo;                  // temporary lower energy limit (applies to one event)
  double Elim_Thi;                  // temporary higher energy limit (applies to one event)
};
 
}  // namespace RAT
 
#endif