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 // G4RToEConvForProton class implementation 26 // G4RToEConvForProton class implementation 27 // 27 // 28 // Author: H.Kurashige, 05 October 2002 - Firs 28 // Author: H.Kurashige, 05 October 2002 - First implementation 29 // ------------------------------------------- 29 // -------------------------------------------------------------------- 30 30 31 #include "G4RToEConvForProton.hh" 31 #include "G4RToEConvForProton.hh" 32 #include "G4ParticleTable.hh" 32 #include "G4ParticleTable.hh" >> 33 #include "G4Material.hh" >> 34 #include "G4PhysicsLogVector.hh" 33 35 >> 36 #include "G4ios.hh" 34 #include "G4PhysicalConstants.hh" 37 #include "G4PhysicalConstants.hh" 35 #include "G4SystemOfUnits.hh" 38 #include "G4SystemOfUnits.hh" 36 39 37 // ------------------------------------------- 40 // -------------------------------------------------------------------- 38 G4RToEConvForProton::G4RToEConvForProton() 41 G4RToEConvForProton::G4RToEConvForProton() 39 : G4VRangeToEnergyConverter() 42 : G4VRangeToEnergyConverter() 40 { 43 { 41 theParticle = G4ParticleTable::GetParticleTa 44 theParticle = G4ParticleTable::GetParticleTable()->FindParticle("proton"); 42 if (theParticle == nullptr) 45 if (theParticle == nullptr) 43 { 46 { 44 #ifdef G4VERBOSE 47 #ifdef G4VERBOSE 45 if (GetVerboseLevel()>0) 48 if (GetVerboseLevel()>0) 46 { 49 { 47 G4cout << "G4RToEConvForProton::G4RToECo 50 G4cout << "G4RToEConvForProton::G4RToEConvForProton() - "; 48 G4cout << "Proton is not defined !!" << 51 G4cout << "Proton is not defined !!" << G4endl; 49 } 52 } 50 #endif 53 #endif 51 } 54 } 52 else << 55 else 53 { 56 { 54 fPDG = theParticle->GetPDGEncoding(); << 57 Mass = theParticle->GetPDGMass(); 55 } << 58 } 56 } 59 } 57 60 58 // ------------------------------------------- 61 // -------------------------------------------------------------------- 59 G4RToEConvForProton::~G4RToEConvForProton() 62 G4RToEConvForProton::~G4RToEConvForProton() 60 {} << 63 { >> 64 } 61 65 62 // ------------------------------------------- 66 // -------------------------------------------------------------------- 63 G4double G4RToEConvForProton::Convert(const G4 << 67 G4double G4RToEConvForProton::Convert(G4double rangeCut, const G4Material* ) 64 const G4 << 65 { 68 { 66 #ifdef G4VERBOSE << 67 if (GetVerboseLevel()>3) << 68 { << 69 G4cout << "G4RToEConvForProton::Convert() << 70 G4cout << " with Range Cut " << rangeCut/m << 71 } << 72 #endif << 73 // Simple formula - range = Ekin/(100*keV)*( 69 // Simple formula - range = Ekin/(100*keV)*(1*mm); 74 return (rangeCut/(1.0*CLHEP::mm)) * (100.0*C << 70 >> 71 return (rangeCut/(1.0*mm)) * (100.0*keV); 75 } 72 } 76 73 77 // ------------------------------------------- << 74 // ********************************************************************** 78 G4double G4RToEConvForProton::ComputeValue(con << 75 // ************************* ComputeLoss ******************************** >> 76 // ********************************************************************** >> 77 G4double G4RToEConvForProton::ComputeLoss(G4double AtomicNumber, >> 78 G4double KineticEnergy) 79 { 79 { 80 return 0.0; << 80 // calculate dE/dx >> 81 const G4double z2Particle = 1.0; >> 82 >> 83 if( std::fabs(AtomicNumber-Z)>0.1 ) >> 84 { >> 85 // recalculate constants >> 86 Z = AtomicNumber; >> 87 G4double Z13 = std::exp(std::log(Z)/3.); >> 88 tau0 = 0.1*Z13*MeV/proton_mass_c2; >> 89 taum = 0.035*Z13*MeV/proton_mass_c2; >> 90 taul = 2.*MeV/proton_mass_c2; >> 91 ionpot = 1.6e-5*MeV*std::exp(0.9*std::log(Z)); >> 92 cc = (taul+1.)*(taul+1.)*std::log(2.*electron_mass_c2*taul*(taul+2.)/ionpot) >> 93 / (taul*(taul+2.))-1.; >> 94 cc = 2.*twopi_mc2_rcl2*Z*cc*std::sqrt(taul); >> 95 ca = cc/((1.-0.5*std::sqrt(tau0/taum))*tau0); >> 96 cba = -0.5/std::sqrt(taum); >> 97 } >> 98 >> 99 G4double tau = KineticEnergy/Mass; >> 100 G4double dEdx; >> 101 if ( tau <= tau0 ) >> 102 { >> 103 dEdx = ca*(std::sqrt(tau)+cba*tau); >> 104 } >> 105 else >> 106 { >> 107 if( tau <= taul ) >> 108 { >> 109 dEdx = cc/std::sqrt(tau); >> 110 } >> 111 else >> 112 { >> 113 dEdx = (tau+1.)*(tau+1.) >> 114 * std::log(2.*electron_mass_c2*tau*(tau+2.)/ionpot) >> 115 / (tau*(tau+2.))-1.; >> 116 dEdx = 2.*twopi_mc2_rcl2*Z*dEdx; >> 117 } >> 118 } >> 119 return dEdx*z2Particle ; 81 } 120 } 82 121 83 // ------------------------------------------- << 122 // ********************************************************************** >> 123 // ************************* Reset ******************************** >> 124 // ********************************************************************** >> 125 void G4RToEConvForProton::Reset() >> 126 { >> 127 // do nothing because loss tables and range vectors are not used >> 128 >> 129 return; >> 130 } 84 131