// // ******************************************************************** // * License and Disclaimer * // * * // * The Geant4 software is copyright of the Copyright Holders of * // * the Geant4 Collaboration. It is provided under the terms and * // * conditions of the Geant4 Software License, included in the file * // * LICENSE and available at http://cern.ch/geant4/license . These * // * include a list of copyright holders. * // * * // * Neither the authors of this software system, nor their employing * // * institutes,nor the agencies providing financial support for this * // * work make any representation or warranty, express or implied, * // * regarding this software system or assume any liability for its * // * use. Please see the license in the file LICENSE and URL above * // * for the full disclaimer and the limitation of liability. * // * * // * This code implementation is the result of the scientific and * // * technical work of the GEANT4 collaboration. * // * By using, copying, modifying or distributing the software (or * // * any work based on the software) you agree to acknowledge its * // * use in resulting scientific publications, and indicate your * // * acceptance of all terms of the Geant4 Software license. * // ******************************************************************** // // $Id: G4VEmModel.hh,v 1.48 2007/10/29 08:38:58 vnivanch Exp $ // GEANT4 tag $Name: geant4-09-01-patch-02 $ // // ------------------------------------------------------------------- // // GEANT4 Class header file // // // File name: G4VEmModel // // Author: Vladimir Ivanchenko // // Creation date: 03.01.2002 // // Modifications: // // 23-12-02 V.Ivanchenko change interface before move to cut per region // 24-01-03 Cut per region (V.Ivanchenko) // 13-02-03 Add name (V.Ivanchenko) // 25-02-03 Add sample theta and displacement (V.Ivanchenko) // 23-07-03 Replace G4Material by G4MaterialCutCouple in dE/dx and CrossSection // calculation (V.Ivanchenko) // 01-03-04 L.Urban signature changed in SampleCosineTheta // 23-04-04 L.urban signature of SampleCosineTheta changed back // 17-11-04 Add method CrossSectionPerAtom (V.Ivanchenko) // 14-03-05 Reduce number of pure virtual methods and make inline part // separate (V.Ivanchenko) // 24-03-05 Remove IsInCharge and add G4VParticleChange in the constructor (VI) // 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko) // 15-04-05 optimize internal interface for msc (V.Ivanchenko) // 08-05-05 A -> N (V.Ivanchenko) // 25-07-05 Move constructor and destructor to the body (V.Ivanchenko) // 02-02-06 ComputeCrossSectionPerAtom: default value A=0. (mma) // 06-02-06 add method ComputeMeanFreePath() (mma) // 07-03-06 Optimize msc methods (V.Ivanchenko) // 29-06-06 Add member currentElement and Get/Set methods (V.Ivanchenko) // 29-10-07 Added SampleScattering (V.Ivanchenko) // // Class Description: // // Abstract interface to energy loss models // ------------------------------------------------------------------- // #ifndef G4VEmModel_h #define G4VEmModel_h 1 #include "globals.hh" #include "G4DynamicParticle.hh" #include "G4ParticleDefinition.hh" #include "G4MaterialCutsCouple.hh" #include "G4Material.hh" #include "G4Element.hh" #include "G4ElementVector.hh" #include "G4DataVector.hh" #include "G4VEmFluctuationModel.hh" #include "Randomize.hh" class G4PhysicsTable; class G4Region; class G4VParticleChange; class G4Track; class G4VEmModel { public: G4VEmModel(const G4String& nam); virtual ~G4VEmModel(); //------------------------------------------------------------------------ // Virtual methods to be implemented for the concrete model //------------------------------------------------------------------------ virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&) = 0; virtual void SampleSecondaries(std::vector*, const G4MaterialCutsCouple*, const G4DynamicParticle*, G4double tmin = 0.0, G4double tmax = DBL_MAX) = 0; //------------------------------------------------------------------------ // Methods with standard implementation; may be overwritten if needed //------------------------------------------------------------------------ virtual G4double MinEnergyCut(const G4ParticleDefinition*, const G4MaterialCutsCouple*); virtual G4double ComputeDEDX(const G4MaterialCutsCouple*, const G4ParticleDefinition*, G4double kineticEnergy, G4double cutEnergy = DBL_MAX); virtual G4double CrossSection(const G4MaterialCutsCouple*, const G4ParticleDefinition*, G4double kineticEnergy, G4double cutEnergy = 0.0, G4double maxEnergy = DBL_MAX); virtual G4double ComputeDEDXPerVolume(const G4Material*, const G4ParticleDefinition*, G4double kineticEnergy, G4double cutEnergy = DBL_MAX); virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*, G4double kinEnergy, G4double Z, G4double A = 0., G4double cutEnergy = 0.0, G4double maxEnergy = DBL_MAX); virtual G4double ComputeMeanFreePath(const G4ParticleDefinition*, G4double kineticEnergy, const G4Material*, G4double cutEnergy = 0.0, G4double maxEnergy = DBL_MAX); virtual G4double CrossSectionPerVolume(const G4Material*, const G4ParticleDefinition*, G4double kineticEnergy, G4double cutEnergy = 0.0, G4double maxEnergy = DBL_MAX); protected: virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition*, G4double kineticEnergy); //------------------------------------------------------------------------ // Methods for msc simulation which needs to be overwritten //------------------------------------------------------------------------ public: virtual void SampleScattering(const G4DynamicParticle*, G4double safety); virtual G4double ComputeTruePathLengthLimit(const G4Track& track, G4PhysicsTable* theLambdaTable, G4double currentMinimalStep); virtual G4double ComputeGeomPathLength(G4double truePathLength); virtual G4double ComputeTrueStepLength(G4double geomPathLength); virtual void DefineForRegion(const G4Region*); //------------------------------------------------------------------------ // Generic methods common to all models //------------------------------------------------------------------------ void SetParticleChange(G4VParticleChange*, G4VEmFluctuationModel*); G4VEmFluctuationModel* GetModelOfFluctuations(); G4double HighEnergyLimit(); G4double LowEnergyLimit(); void SetHighEnergyLimit(G4double); void SetLowEnergyLimit(G4double); G4double MaxSecondaryKinEnergy(const G4DynamicParticle* dynParticle); const G4Element* SelectRandomAtom(const G4Material*, const G4ParticleDefinition*, G4double kineticEnergy, G4double cutEnergy = 0.0, G4double maxEnergy = DBL_MAX); const G4String& GetName() const; protected: const G4Element* GetCurrentElement() const; void SetCurrentElement(const G4Element*); private: // hide assignment operator G4VEmModel & operator=(const G4VEmModel &right); G4VEmModel(const G4VEmModel&); G4double lowLimit; G4double highLimit; G4double xsec[40]; G4VEmFluctuationModel* fluc; const G4String name; const G4Element* currentElement; protected: G4VParticleChange* pParticleChange; }; //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::HighEnergyLimit() { return highLimit; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::LowEnergyLimit() { return lowLimit; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline void G4VEmModel::SetHighEnergyLimit(G4double val) { highLimit = val; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline void G4VEmModel::SetLowEnergyLimit(G4double val) { lowLimit = val; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline void G4VEmModel::SetParticleChange(G4VParticleChange* p, G4VEmFluctuationModel* f = 0) { if(p && pParticleChange != p) pParticleChange = p; fluc = f; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4VEmFluctuationModel* G4VEmModel::GetModelOfFluctuations() { return fluc; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::MinEnergyCut(const G4ParticleDefinition*, const G4MaterialCutsCouple*) { return 0.0; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::ComputeDEDXPerVolume( const G4Material*, const G4ParticleDefinition*, G4double, G4double) { return 0.0; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::ComputeDEDX(const G4MaterialCutsCouple* c, const G4ParticleDefinition* p, G4double kinEnergy, G4double cutEnergy) { return ComputeDEDXPerVolume(c->GetMaterial(),p,kinEnergy,cutEnergy); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::CrossSection(const G4MaterialCutsCouple* c, const G4ParticleDefinition* p, G4double kinEnergy, G4double cutEnergy, G4double maxEnergy) { return CrossSectionPerVolume(c->GetMaterial(),p,kinEnergy,cutEnergy,maxEnergy); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::ComputeCrossSectionPerAtom( const G4ParticleDefinition*, G4double, G4double, G4double, G4double, G4double) { return 0.0; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline const G4Element* G4VEmModel::SelectRandomAtom( const G4Material* material, const G4ParticleDefinition* pd, G4double kinEnergy, G4double tcut, G4double tmax) { const G4ElementVector* theElementVector = material->GetElementVector(); G4int nelm = material->GetNumberOfElements(); currentElement = (*theElementVector)[nelm-1]; if (nelm > 1) { G4double x = G4UniformRand()* CrossSectionPerVolume(material,pd,kinEnergy,tcut,tmax); for(G4int i=0; iGetDefinition(), dynParticle->GetKineticEnergy()); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::MaxSecondaryEnergy(const G4ParticleDefinition*, G4double kineticEnergy) { return kineticEnergy; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline const G4String& G4VEmModel::GetName() const { return name; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... // Methods for msc simulation //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline void G4VEmModel::SampleScattering(const G4DynamicParticle*, G4double) {} //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::ComputeTruePathLengthLimit( const G4Track&, G4PhysicsTable*, G4double) { return DBL_MAX; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::ComputeGeomPathLength(G4double truePathLength) { return truePathLength; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline G4double G4VEmModel::ComputeTrueStepLength(G4double geomPathLength) { return geomPathLength; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... inline void G4VEmModel::DefineForRegion(const G4Region*) {} //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... #endif