Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // 27 //--------------------------------------------------------------------------- 28 // 29 // ClassName: G4HadronElasticPhysicsHPT 30 // 31 // Author: Alberto Ribon (CERN), April 2023 32 // 33 // This class inherits from G4HadronElasticPhysicsHP and activates 34 // the special treatment of elastic scattering of thermal neutrons 35 // (i.e. with kinetic energy below 4 eV). 36 // This special treatment, called Thermal Scattering Law (TSL), is based on 37 // the S(alpha, beta) approach, which relies on both experimental measurements 38 // and molecular dynamics calculations. 39 // 40 // Modified: 41 // 42 //---------------------------------------------------------------------------- 43 // 44 45 #include "G4HadronElasticPhysicsHPT.hh" 46 #include "G4SystemOfUnits.hh" 47 #include "G4Neutron.hh" 48 #include "G4ParticleDefinition.hh" 49 #include "G4HadronicProcess.hh" 50 #include "G4ParticleHPThermalScattering.hh" 51 #include "G4ParticleHPThermalScatteringData.hh" 52 #include "G4HadronicParameters.hh" 53 #include "G4PhysListUtil.hh" 54 #include "G4PhysicsConstructorFactory.hh" 55 56 57 G4_DECLARE_PHYSCONSTR_FACTORY( G4HadronElasticPhysicsHPT ); 58 59 60 G4HadronElasticPhysicsHPT::G4HadronElasticPhysicsHPT( G4int ver ) : G4HadronElasticPhysicsHP( ver ) { 61 if ( ver > 1 ) G4cout << "### G4HadronElasticPhysicsHPT: " << GetPhysicsName() << G4endl; 62 } 63 64 65 void G4HadronElasticPhysicsHPT::ConstructProcess() { 66 G4HadronElasticPhysicsHP::ConstructProcess(); 67 const G4Neutron* neutron = G4Neutron::Neutron(); 68 G4HadronicProcess* neutronElasticProcess = G4PhysListUtil::FindElasticProcess( neutron ); 69 if ( neutronElasticProcess == nullptr ) { 70 G4cout << "### " << GetPhysicsName() 71 << " WARNING: Fail to add thermal neutron scattering" << G4endl; 72 return; 73 } 74 std::size_t sizeInteractionList = neutronElasticProcess->GetHadronicInteractionList().size(); 75 if ( sizeInteractionList < 1 ) { 76 G4cout << "### " << GetPhysicsName() 77 << " WARNING: Fail to add thermal neutron scattering ! sizeInteractionList=" 78 << sizeInteractionList << G4endl; 79 return; 80 } 81 neutronElasticProcess->GetHadronicInteractionList()[ sizeInteractionList-1 ]->SetMinEnergy( 4.0*CLHEP::eV ); 82 neutronElasticProcess->RegisterMe( new G4ParticleHPThermalScattering ); 83 neutronElasticProcess->AddDataSet( new G4ParticleHPThermalScatteringData ); 84 if ( G4HadronicParameters::Instance()->GetVerboseLevel() > 1 ) { 85 G4cout << "### HadronElasticPhysicsHPT is constructed " << G4endl; 86 } 87 } 88