NestAlg.h

Go to the documentation of this file.

#ifndef NESTALG_H
#define NESTALG_H

#include "Geant4/G4Types.hh"
#include "Geant4/G4VParticleChange.hh"
#include <map>

class G4MaterialPropertiesTable;
class G4Step;
class G4Track;

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