Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 // 26 // 27 // ------------------------------------------- 27 // ------------------------------------------------------------------- 28 // 28 // 29 // GEANT4 Class file 29 // GEANT4 Class file 30 // 30 // 31 // 31 // 32 // File name: G4PhotoElectricAngularGenera 32 // File name: G4PhotoElectricAngularGeneratorSauterGavrila 33 // 33 // 34 // Creation date: 10 May 2004 34 // Creation date: 10 May 2004 35 // 35 // 36 // Modifications: 36 // Modifications: 37 // 10 May 2003 P. Rodrigues First imple 37 // 10 May 2003 P. Rodrigues First implementation acording with new design 38 // 38 // 39 // Class Description: 39 // Class Description: 40 // 40 // 41 // Concrete class for PhotoElectric Electron A 41 // Concrete class for PhotoElectric Electron Angular Distribution Generation 42 // This model is a re-implementation of the Ph 42 // This model is a re-implementation of the Photolectric angular distribution 43 // developed my M. Maire for the Standard EM P 43 // developed my M. Maire for the Standard EM Physics G4PhotoElectricEffect 44 // 44 // 45 // Class Description: End 45 // Class Description: End 46 // 46 // 47 // ------------------------------------------- 47 // ------------------------------------------------------------------- 48 // 48 // >> 49 // 49 50 50 #include "G4PhotoElectricAngularGeneratorSaute 51 #include "G4PhotoElectricAngularGeneratorSauterGavrila.hh" 51 #include "G4PhysicalConstants.hh" << 52 #include "Randomize.hh" 52 #include "Randomize.hh" 53 53 54 // ------------------------------------------- << 54 // 55 G4PhotoElectricAngularGeneratorSauterGavrila:: << 56 G4VEmAngularDistribution("AngularGenSauterGa << 57 {} << 58 55 59 // ------------------------------------------- << 56 G4PhotoElectricAngularGeneratorSauterGavrila::G4PhotoElectricAngularGeneratorSauterGavrila(const G4String& name):G4VPhotoElectricAngularDistribution(name) >> 57 {;} >> 58 >> 59 // 60 60 61 G4PhotoElectricAngularGeneratorSauterGavrila:: 61 G4PhotoElectricAngularGeneratorSauterGavrila::~G4PhotoElectricAngularGeneratorSauterGavrila() 62 {} << 62 {;} 63 63 64 // ------------------------------------------- << 64 // 65 65 66 G4ThreeVector& << 66 G4ThreeVector G4PhotoElectricAngularGeneratorSauterGavrila::GetPhotoElectronDirection(const G4ThreeVector& direction, const G4double eKineticEnergy, const G4ThreeVector&, const G4int) const 67 G4PhotoElectricAngularGeneratorSauterGavrila:: << 68 const G4DynamicPartic << 69 G4double, G4int, cons << 70 { 67 { 71 68 72 // Compute Theta distribution of the emitted 69 // Compute Theta distribution of the emitted electron, with respect to the 73 // incident Gamma. 70 // incident Gamma. 74 // The Sauter-Gavrila distribution for the K << 71 // The Sauter-Gavrila distribution for the K-shell is used. (adapted from G4PhotoElectricEffect) >> 72 75 G4double costeta = 1.; 73 G4double costeta = 1.; 76 G4double Phi = twopi * G4UniformRand(); 74 G4double Phi = twopi * G4UniformRand(); 77 G4double cosphi = std::cos(Phi); 75 G4double cosphi = std::cos(Phi); 78 G4double sinphi = std::sin(Phi); 76 G4double sinphi = std::sin(Phi); 79 G4double sinteta = 0; 77 G4double sinteta = 0; 80 G4double gamma = 1. + dp->GetKineticEnergy << 78 G4double gamma = 1. + eKineticEnergy/electron_mass_c2; 81 79 82 if (gamma > 5.) { 80 if (gamma > 5.) { 83 fLocalDirection = dp->GetMomentumDirection << 81 G4ThreeVector direction (sinteta*cosphi, sinteta*sinphi, costeta); 84 return fLocalDirection; << 82 return costeta; 85 // Bugzilla 1120 << 86 // SI on 05/09/2010 as suggested by JG 04/ << 87 } 83 } 88 84 89 G4double beta = std::sqrt((gamma - 1)*(gamm << 85 G4double beta = std::sqrt(gamma*gamma-1.)/gamma; 90 G4double b = 0.5*gamma*(gamma - 1)*(gamm << 86 G4double b = 0.5*gamma*(gamma-1.)*(gamma-2); 91 87 92 G4double rndm,term,greject,grejsup; 88 G4double rndm,term,greject,grejsup; 93 if (gamma < 2.) grejsup = gamma*gamma*(1.+b- 89 if (gamma < 2.) grejsup = gamma*gamma*(1.+b-beta*b); 94 else grejsup = gamma*gamma*(1.+b+ 90 else grejsup = gamma*gamma*(1.+b+beta*b); 95 91 96 do { rndm = 1.-2*G4UniformRand(); 92 do { rndm = 1.-2*G4UniformRand(); 97 costeta = (rndm+beta)/(rndm*beta+1.); 93 costeta = (rndm+beta)/(rndm*beta+1.); 98 term = 1.-beta*costeta; 94 term = 1.-beta*costeta; 99 greject = (1.-costeta*costeta)*(1.+b*te 95 greject = (1.-costeta*costeta)*(1.+b*term)/(term*term); 100 } while(greject < G4UniformRand()*grejsup); 96 } while(greject < G4UniformRand()*grejsup); 101 97 102 sinteta = std::sqrt((1 - costeta)*(1 + coste << 98 103 fLocalDirection.set(sinteta*cosphi, sinteta* << 99 sinteta = std::sqrt(1.-costeta*costeta); 104 fLocalDirection.rotateUz(dp->GetMomentumDire << 100 G4ThreeVector photoelectrondirection (sinteta*cosphi, sinteta*sinphi, costeta); 105 return fLocalDirection; << 101 photoelectrondirection.rotateUz(direction); >> 102 return photoelectrondirection; 106 } 103 } 107 104 108 // ------------------------------------------- << 105 // 109 106 110 void G4PhotoElectricAngularGeneratorSauterGavr 107 void G4PhotoElectricAngularGeneratorSauterGavrila::PrintGeneratorInformation() const 111 { 108 { 112 G4cout << "\n" << G4endl; 109 G4cout << "\n" << G4endl; 113 G4cout << "" << G4endl; 110 G4cout << "" << G4endl; 114 G4cout << "Re-implementation of the photolec 111 G4cout << "Re-implementation of the photolectric angular distribution" << G4endl; 115 G4cout << "developed my M. Maire for the Sta 112 G4cout << "developed my M. Maire for the Standard EM Physics G4PhotoElectricEffect" << G4endl; 116 G4cout << "It computes the theta distributio 113 G4cout << "It computes the theta distribution of the emitted electron, with respect to the" << G4endl; 117 G4cout << "incident Gamma, using the Sauter- 114 G4cout << "incident Gamma, using the Sauter-Gavrila distribution for the K-shell\n" << G4endl; 118 } 115 } 119 116