Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/binary_cascade/src/G4NeutronField.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/hadronic/models/binary_cascade/src/G4NeutronField.cc (Version 11.3.0) and /processes/hadronic/models/binary_cascade/src/G4NeutronField.cc (Version 9.3)


  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 //      GEANT 4 class implementation file          28 //      GEANT 4 class implementation file 
 29 //                                                 29 //
 30 //      CERN, Geneva, Switzerland                  30 //      CERN, Geneva, Switzerland
 31 //                                                 31 //
 32 //      File name:     G4NeutronField.cc           32 //      File name:     G4NeutronField.cc
 33 //                                                 33 //
 34 //      Author:        Alessandro Brunengo (Al     34 //      Author:        Alessandro Brunengo (Alessandro.Brunengo@ge.infn.it)
 35 //                                                 35 // 
 36 //      Creation date: 5 June 2000                 36 //      Creation date: 5 June 2000
 37 // -------------------------------------------     37 // -------------------------------------------------------------------
 38                                                << 
 39 #include "G4NeutronField.hh"                       38 #include "G4NeutronField.hh"
 40 #include "G4PhysicalConstants.hh"              << 
 41 #include "G4SystemOfUnits.hh"                  << 
 42 #include "G4VNuclearDensity.hh"                    39 #include "G4VNuclearDensity.hh"
 43 #include "G4FermiMomentum.hh"                      40 #include "G4FermiMomentum.hh"
 44                                                    41 
                                                   >>  42 
 45 G4NeutronField::G4NeutronField(G4V3DNucleus *      43 G4NeutronField::G4NeutronField(G4V3DNucleus * aNucleus) : 
 46     G4VNuclearField(aNucleus), theDensity(theN     44     G4VNuclearField(aNucleus), theDensity(theNucleus->GetNuclearDensity())
 47 {                                                  45 { 
 48   theA = theNucleus->GetMassNumber();              46   theA = theNucleus->GetMassNumber();
 49   theZ = theNucleus->GetCharge();                  47   theZ = theNucleus->GetCharge();
 50   theFermi.Init(theA, theZ);                       48   theFermi.Init(theA, theZ);
 51   theR = 2.*theNucleus->GetOuterRadius();          49   theR = 2.*theNucleus->GetOuterRadius();
 52   for (G4double aR=0.;aR<theR; aR+=0.3*fermi)  <<  50   G4double aR=0;
                                                   >>  51   while(aR<theR)
 53   {                                                52   {
 54     G4ThreeVector aPosition(0,0,aR);               53     G4ThreeVector aPosition(0,0,aR);
 55     G4double density = GetDensity(aPosition);      54     G4double density = GetDensity(aPosition);
 56     G4double fermiMom = GetFermiMomentum(densi     55     G4double fermiMom = GetFermiMomentum(density);
 57     theFermiMomBuffer.push_back(fermiMom);         56     theFermiMomBuffer.push_back(fermiMom);
                                                   >>  57     aR+=0.3*fermi;
 58   }                                                58   }
 59   {                                                59   {
 60   G4ThreeVector aPosition(0,0,theR);               60   G4ThreeVector aPosition(0,0,theR);
 61   G4double density = GetDensity(aPosition);        61   G4double density = GetDensity(aPosition);
 62   G4double fermiMom = GetFermiMomentum(density     62   G4double fermiMom = GetFermiMomentum(density);
 63   theFermiMomBuffer.push_back(fermiMom);           63   theFermiMomBuffer.push_back(fermiMom);  
 64   }                                                64   }
 65   {                                                65   {
 66   G4ThreeVector aPosition(0,0,theR+0.001*fermi     66   G4ThreeVector aPosition(0,0,theR+0.001*fermi);
 67   theFermiMomBuffer.push_back(0);                  67   theFermiMomBuffer.push_back(0);  
 68   }                                                68   }
 69   {                                                69   {
 70   G4ThreeVector aPosition(0,0,1.*m);               70   G4ThreeVector aPosition(0,0,1.*m);
 71   theFermiMomBuffer.push_back(0);                  71   theFermiMomBuffer.push_back(0);  
 72   }                                                72   }
 73 }                                                  73 }
 74                                                    74 
 75 G4NeutronField::~G4NeutronField()                  75 G4NeutronField::~G4NeutronField()
 76 { }                                                76 { }
 77                                                    77 
                                                   >>  78 
                                                   >>  79 const G4NeutronField & G4NeutronField::operator=(const G4NeutronField &)
                                                   >>  80 {
                                                   >>  81   throw G4HadronicException(__FILE__, __LINE__, "G4NeutronField::operator= meant not to be accessible");
                                                   >>  82   return *this;
                                                   >>  83 }
                                                   >>  84 
                                                   >>  85 
                                                   >>  86 G4int G4NeutronField::operator==(const G4NeutronField &) const
                                                   >>  87 {
                                                   >>  88   throw G4HadronicException(__FILE__, __LINE__, "G4NeutronField::operator== meant not to be accessible");
                                                   >>  89   return 0;
                                                   >>  90 }
                                                   >>  91 
                                                   >>  92 
                                                   >>  93 G4int G4NeutronField::operator!=(const G4NeutronField &) const
                                                   >>  94 {
                                                   >>  95   throw G4HadronicException(__FILE__, __LINE__, "G4NeutronField::operator!= meant not to be accessible");
                                                   >>  96   return 1;
                                                   >>  97 }
                                                   >>  98 
                                                   >>  99 
 78 G4double G4NeutronField::GetField(const G4Thre    100 G4double G4NeutronField::GetField(const G4ThreeVector & aPosition)
 79 {                                                 101 {
 80   G4double x = aPosition.mag();                   102   G4double x = aPosition.mag();
 81   unsigned int index = static_cast<unsigned in << 103   G4int index = static_cast<G4int>(x/(0.3*fermi) );
 82   if( (index+2) > theFermiMomBuffer.size()) re << 104   if(index+2> static_cast<G4int>(theFermiMomBuffer.size())) return theFermiMomBuffer.back();
 83   G4double y1 = theFermiMomBuffer[index];         105   G4double y1 = theFermiMomBuffer[index];
 84   G4double y2 = theFermiMomBuffer[index+1];       106   G4double y2 = theFermiMomBuffer[index+1];
 85   G4double x1 = (0.3*fermi)*index;                107   G4double x1 = (0.3*fermi)*index;
 86   G4double x2 = (0.3*fermi)*(index+1);            108   G4double x2 = (0.3*fermi)*(index+1);
 87   G4double fermiMom = y1 + (x-x1)*(y2-y1)/(x2-    109   G4double fermiMom = y1 + (x-x1)*(y2-y1)/(x2-x1);
 88   return -1*(fermiMom*fermiMom)/(2*neutron_mas    110   return -1*(fermiMom*fermiMom)/(2*neutron_mass_c2);
 89 }                                                 111 }
                                                   >> 112 
                                                   >> 113 
 90                                                   114 
 91 G4double G4NeutronField::GetBarrier()             115 G4double G4NeutronField::GetBarrier()
 92 {                                                 116 {
 93 /*                                                117 /*
 94  *   G4double A = theNucleus->GetMassNumber();    118  *   G4double A = theNucleus->GetMassNumber();
 95  *   G4double Z = theNucleus->GetCharge();        119  *   G4double Z = theNucleus->GetCharge();
 96  *                                                120  * 
 97  *   return G4NucleiPropertiesTable::GetBindin    121  *   return G4NucleiPropertiesTable::GetBindingEnergy(Z, A)/A;
 98  */                                               122  */
 99   return 0.;                                      123   return 0.;
100 }                                                 124 }
101                                                   125 
102                                                   126 
103                                                   127 
104                                                   128 
105                                                   129 
106                                                   130 
107                                                   131