NestAlg.h

Go to the documentation of this file.

#ifndef NESTALG_H
#define NESTALG_H

#include "Geant4/globals.hh"
#include "Geant4/G4Step.hh"
#include "Geant4/G4Track.hh"
#include "Geant4/G4VParticleChange.hh"

// forward declarations
namespace CLHEP { class HepRandomEngine; }

class NestAlg {

 public :

  NestAlg(CLHEP::HepRandomEngine& engine);
  NestAlg(double yieldFactor, CLHEP::HepRandomEngine& engine);

  const G4VParticleChange& CalculateIonizationAndScintillation(G4Track const& aTrack,
                                                               G4Step  const& aStep);

  void   SetScintillationYieldFactor(double const& yf)     { fYieldFactor = yf;       } 
  void   SetScintillationExcitationRatio(double const& er) { fExcitationRatio = er;   } 
  int    NumberScintillationPhotons() const                { return fNumScintPhotons; }
  int    NumberIonizationElectrons()  const                { return fNumIonElectrons; }
  double EnergyDeposition()           const                { return fEnergyDep;       }

 private:

  G4double GetGasElectronDriftSpeed(G4double efieldinput,G4double density);
  G4double GetLiquidElectronDriftSpeed(double T, double F, G4bool M, G4int Z);
  G4double CalculateElectronLET ( G4double E, G4int Z );
  G4double UnivScreenFunc ( G4double E, G4double Z, G4double A );
  G4int BinomFluct(G4int N0, G4double prob); //function for doing fluctuations
  void InitMatPropValues ( G4MaterialPropertiesTable* nobleElementMat, int z );

  double             fYieldFactor;     
  double             fExcitationRatio; 
  int                fNumScintPhotons; 
  int                fNumIonElectrons; 
  double             fEnergyDep;       
  G4VParticleChange  fParticleChange;  
  std::map<int,bool> fElementPropInit; 
  CLHEP::HepRandomEngine& fEngine;     
};

#endif //NESTALG_H