Ignore:
Timestamp:
May 28, 2009, 4:26:57 PM (15 years ago)
Author:
garnier
Message:

maj sur la beta de geant 4.9.3

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/source/processes/hadronic/models/coherent_elastic/include/G4DiffuseElastic.hh

    r1007 r1055  
    2525//
    2626//
    27 // $Id: G4DiffuseElastic.hh,v 1.13 2007/11/06 17:01:20 grichine Exp $
    28 // GEANT4 tag $Name: geant4-09-02 $
     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 $
    2929//
    3030//
     
    164164                                G4double tmass, G4double thetaLab);
    165165
    166 
     166  void TestAngleTable(const G4ParticleDefinition* theParticle, G4double partMom,
     167                      G4double Z, G4double A);
    167168
    168169
     
    175176  G4double GetDiffElasticProb(G4double theta);
    176177  G4double GetDiffElasticSumProb(G4double theta);
     178  G4double GetDiffElasticSumProbA(G4double alpha);
    177179  G4double GetIntegrandFunction(G4double theta);
    178180
     
    443445  G4double r0;
    444446
    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;
    449451
    450452    fNuclearRadius = r0*std::pow(A, 1./3.);
     
    452454  else
    453455  {
    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);
    457459  }
    458460  return fNuclearRadius;
Note: See TracChangeset for help on using the changeset viewer.