00001 // 00002 // ******************************************************************** 00003 // * License and Disclaimer * 00004 // * * 00005 // * The Geant4 software is copyright of the Copyright Holders of * 00006 // * the Geant4 Collaboration. It is provided under the terms and * 00007 // * conditions of the Geant4 Software License, included in the file * 00008 // * LICENSE and available at http://cern.ch/geant4/license . These * 00009 // * include a list of copyright holders. * 00010 // * * 00011 // * Neither the authors of this software system, nor their employing * 00012 // * institutes,nor the agencies providing financial support for this * 00013 // * work make any representation or warranty, express or implied, * 00014 // * regarding this software system or assume any liability for its * 00015 // * use. Please see the license in the file LICENSE and URL above * 00016 // * for the full disclaimer and the limitation of liability. * 00017 // * * 00018 // * This code implementation is the result of the scientific and * 00019 // * technical work of the GEANT4 collaboration. * 00020 // * By using, copying, modifying or distributing the software (or * 00021 // * any work based on the software) you agree to acknowledge its * 00022 // * use in resulting scientific publications, and indicate your * 00023 // * acceptance of all terms of the Geant4 Software license. * 00024 // ******************************************************************** 00025 // 00026 // 00027 // ------------------------------------------------------------------- 00028 // $Id$ 00029 // 00030 // Author: L.Pandola 00031 // 00032 // History: 00033 // ----------- 00034 // 23 Nov 2010 L. Pandola 1st implementation 00035 // 24 May 2011 L. Pandola Renamed (make default Penelope) 00036 // 13 Mar 2012 L. Pandola Made a derived class of G4VEmAngularDistribution 00037 // and update the interface accordingly 00038 // 18 Jul 2012 L. Pandola Migrate to the new interface of G4VEmAngularDistribution 00039 // 00040 // Class description: 00041 // Calculation of angular distribution for Penelope Bremsstrahlung 00042 // version 2008 00043 // -------------------------------------------------------------- 00044 00045 00046 #ifndef G4PENELOPEBREMSSTRAHLUNGANGULAR_HH 00047 #define G4PENELOPEBREMSSTRAHLUNGANGULAR_HH 1 00048 #include "globals.hh" 00049 #include <map> 00050 #include "G4VEmAngularDistribution.hh" 00051 00052 class G4PhysicsTable; 00053 class G4Material; 00054 00055 class G4PenelopeBremsstrahlungAngular : public G4VEmAngularDistribution 00056 { 00057 00058 public: 00059 G4PenelopeBremsstrahlungAngular(); 00060 ~G4PenelopeBremsstrahlungAngular(); 00061 00063 void Initialize(); 00064 00067 G4double PolarAngle(const G4double initial_energy, 00068 const G4double final_energy, 00069 const G4int Z); 00070 00073 G4ThreeVector& SampleDirection(const G4DynamicParticle* dp, 00074 G4double out_energy, 00075 G4int Z, 00076 const G4Material* mat = 0); 00077 00079 void SetVerbosityLevel(G4int vl){verbosityLevel = vl;}; 00080 G4int GetVerbosityLevel(){return verbosityLevel;}; 00081 00082 private: 00083 void PrepareInterpolationTables(G4double Zeq); 00084 void ClearTables(); 00085 00086 G4double GetEffectiveZ(const G4Material* material); 00087 std::map<const G4Material*,G4double> *theEffectiveZSq; 00088 00089 //Tables containing the Lorentz sampling coefficients 00090 //The key is the effective Z of the material 00091 std::map<G4double,G4PhysicsTable*> *theLorentzTables1; 00092 std::map<G4double,G4PhysicsTable*> *theLorentzTables2; 00093 00094 void ReadDataFile(); 00095 G4bool dataRead; 00096 00097 static const G4int NumberofZPoints=6; 00098 static const G4int NumberofEPoints=6; 00099 static const G4int NumberofKPoints=4; 00100 00101 G4double QQ1[NumberofZPoints][NumberofEPoints][NumberofKPoints]; 00102 G4double QQ2[NumberofZPoints][NumberofEPoints][NumberofKPoints]; 00103 00104 G4int verbosityLevel; 00105 00106 }; 00107 00108 00109 00110 #endif