Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/binary_cascade/src/G4PionMinusField.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/G4PionMinusField.cc (Version 11.3.0) and /processes/hadronic/models/binary_cascade/src/G4PionMinusField.cc (Version 8.1)


  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:     G4PionMinusField.cc         32 //      File name:     G4PionMinusField.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 "G4PionMinusField.hh"                     38 #include "G4PionMinusField.hh"
 40 #include "G4PhysicalConstants.hh"              <<  39 #include "G4NucleiPropertiesTable.hh"
 41 #include "G4SystemOfUnits.hh"                  << 
 42 #include "G4NucleiProperties.hh"               << 
 43 #include "G4VNuclearDensity.hh"                    40 #include "G4VNuclearDensity.hh"
 44 #include "G4FermiMomentum.hh"                      41 #include "G4FermiMomentum.hh"
 45 #include "G4PionMinus.hh"                          42 #include "G4PionMinus.hh"
 46 #include "G4HadTmpUtil.hh"                         43 #include "G4HadTmpUtil.hh"
 47 #include "G4Pow.hh"                            << 
 48                                                    44 
 49 G4PionMinusField::G4PionMinusField(G4V3DNucleu     45 G4PionMinusField::G4PionMinusField(G4V3DNucleus * nucleus, G4double coeff)
 50   : G4VNuclearField(nucleus)                       46   : G4VNuclearField(nucleus)
 51 {                                                  47 {
 52   theCoeff = coeff;                                48   theCoeff = coeff;
 53 }                                                  49 }
 54                                                    50 
 55                                                    51 
 56 G4PionMinusField::~G4PionMinusField()              52 G4PionMinusField::~G4PionMinusField()
 57 { }                                                53 { }
 58                                                    54 
                                                   >>  55 
                                                   >>  56 const G4PionMinusField & G4PionMinusField::operator=(const G4PionMinusField &)
                                                   >>  57 {
                                                   >>  58   throw G4HadronicException(__FILE__, __LINE__, "G4PionMinusField::operator= meant not to be accessible");
                                                   >>  59   return *this;
                                                   >>  60 }
                                                   >>  61 
                                                   >>  62 
                                                   >>  63 G4int G4PionMinusField::operator==(const G4PionMinusField &) const
                                                   >>  64 {
                                                   >>  65   throw G4HadronicException(__FILE__, __LINE__, "G4PionMinusField::operator== meant not to be accessible");
                                                   >>  66   return 0;
                                                   >>  67 }
                                                   >>  68 
                                                   >>  69 
                                                   >>  70 G4int G4PionMinusField::operator!=(const G4PionMinusField &) const
                                                   >>  71 {
                                                   >>  72   throw G4HadronicException(__FILE__, __LINE__, "G4PionMinusField::operator!= meant not to be accessible");
                                                   >>  73   return 1;
                                                   >>  74 }
                                                   >>  75 
                                                   >>  76 
 59 G4double G4PionMinusField::GetField(const G4Th     77 G4double G4PionMinusField::GetField(const G4ThreeVector & aPosition)
 60 {                                                  78 {
 61 // Field is 0 out of the nucleus!                  79 // Field is 0 out of the nucleus!
 62   if(aPosition.mag() >= radius) return 0.0;        80   if(aPosition.mag() >= radius) return 0.0;
 63                                                    81 
 64   G4int A = theNucleus->GetMassNumber();       <<  82   G4double A = theNucleus->GetMassNumber();
 65   G4int Z = theNucleus->GetCharge();           <<  83   G4double Z = theNucleus->GetCharge();
 66   G4double pionMinusMass = G4PionMinus::PionMi     84   G4double pionMinusMass = G4PionMinus::PionMinus()->GetPDGMass();
 67   G4double bindingEnergy = G4NucleiProperties: <<  85   G4double bindingEnergy = G4NucleiPropertiesTable::GetBindingEnergy(G4lrint(Z), G4lrint(A));
 68   G4double nucleusMass = Z*proton_mass_c2+(A-Z     86   G4double nucleusMass = Z*proton_mass_c2+(A-Z)*neutron_mass_c2+bindingEnergy;
 69   G4double reducedMass = pionMinusMass*nucleus     87   G4double reducedMass = pionMinusMass*nucleusMass/(pionMinusMass+nucleusMass);
 70                                                    88 
 71   G4double density = A*theNucleus->GetNuclearD     89   G4double density = A*theNucleus->GetNuclearDensity()->GetDensity(aPosition);
 72   G4double nucleonMass = (proton_mass_c2+neutr     90   G4double nucleonMass = (proton_mass_c2+neutron_mass_c2)/2;
 73                                                    91 
 74   return 2.*pi*hbarc*hbarc/reducedMass*(1+pion     92   return 2.*pi*hbarc*hbarc/reducedMass*(1+pionMinusMass/nucleonMass)*theCoeff*density + GetBarrier();
 75 }                                                  93 }
 76                                                    94 
 77                                                    95 
 78 G4double G4PionMinusField::GetBarrier()            96 G4double G4PionMinusField::GetBarrier()
 79 {                                                  97 {
 80   G4int A = theNucleus->GetMassNumber();       <<  98   G4double A = theNucleus->GetMassNumber();
 81   G4int Z = theNucleus->GetCharge();           <<  99   G4double Z = theNucleus->GetCharge();
 82   G4double coulombBarrier = (1.44/1.14) * MeV  << 100   G4double coulombBarrier = (1.44/1.14) * MeV * Z / (1.0 + std::pow(A,1./3.));
 83   return -coulombBarrier;                         101   return -coulombBarrier;
 84 }                                                 102 }
 85                                                   103 
 86                                                   104 
 87                                                   105 
 88                                                   106 
 89                                                   107 
 90                                                   108 
 91                                                   109 
 92                                                   110 
 93                                                   111 
 94                                                   112 
 95                                                   113 
 96                                                   114 
 97                                                   115 
 98                                                   116 
 99                                                   117 
100                                                   118 
101                                                   119