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; }

namespace gar {
  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;     ///< turns scint. on/off
    double                 fExcitationRatio; ///< N_ex/N_i, the dimensionless ratio of initial
                                         ///< excitons to ions
    int                    fNumScintPhotons; ///< number of photons produced by the step
    int                    fNumIonElectrons; ///< number of ionization electrons produced by step
    double                 fEnergyDep;       ///< energy deposited by the step
    G4VParticleChange  fParticleChange;  ///< pointer to G4VParticleChange
    std::map<int,bool> fElementPropInit; ///< map of noble element z to flag
                                         ///< for whether that element's material
                                         ///< properties table has been initialized
    CLHEP::HepRandomEngine& fEngine;     ///< random engine
  };
  
} // gar

#endif //NESTALG_H