- Timestamp:
- May 28, 2009, 4:26:57 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/coherent_elastic/include/G4DiffuseElastic.hh
r1007 r1055 25 25 // 26 26 // 27 // $Id: G4DiffuseElastic.hh,v 1.1 3 2007/11/06 17:01:20grichine Exp $28 // GEANT4 tag $Name: geant4-09-0 2$27 // $Id: G4DiffuseElastic.hh,v 1.15 2009/02/20 15:38:07 grichine Exp $ 28 // GEANT4 tag $Name: geant4-09-03-beta-cand-01 $ 29 29 // 30 30 // … … 164 164 G4double tmass, G4double thetaLab); 165 165 166 166 void TestAngleTable(const G4ParticleDefinition* theParticle, G4double partMom, 167 G4double Z, G4double A); 167 168 168 169 … … 175 176 G4double GetDiffElasticProb(G4double theta); 176 177 G4double GetDiffElasticSumProb(G4double theta); 178 G4double GetDiffElasticSumProbA(G4double alpha); 177 179 G4double GetIntegrandFunction(G4double theta); 178 180 … … 443 445 G4double r0; 444 446 445 if( A < 50.)446 { 447 if( A > 10.) r0 = 1.16*( 1 - std::pow(A, -2./3.) )*fermi; // 1.08*fermi;448 else r0 = 1.1*fermi;447 if( A < 50. ) 448 { 449 if( A > 10. ) r0 = 1.16*( 1 - std::pow(A, -2./3.) )*fermi; // 1.08*fermi; 450 else r0 = 1.1*fermi; 449 451 450 452 fNuclearRadius = r0*std::pow(A, 1./3.); … … 452 454 else 453 455 { 454 r0 = 1.7*fermi; 455 456 fNuclearRadius = r0*std::pow(A, 0.27); 456 r0 = 1.7*fermi; // 1.7*fermi; 457 458 fNuclearRadius = r0*std::pow(A, 0.27); // 0.27); 457 459 } 458 460 return fNuclearRadius;
Note: See TracChangeset
for help on using the changeset viewer.