Hadr04/src/NeutronHPphysics.cc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 /// \file NeutronHPphysics.cc
0027 /// \brief Implementation of the NeutronHPphysics class
0028 //
0029 //
0030 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0031
0032 #include "NeutronHPphysics.hh"
0033
0034 #include "G4GenericMessenger.hh"
0035 #include "G4ParticleDefinition.hh"
0036 #include "G4ProcessManager.hh"
0037 #include "G4ProcessTable.hh"
0038
0039 // Processes
0040
0041 #include "G4HadronElasticProcess.hh"
0042 #include "G4HadronInelasticProcess.hh"
0043 #include "G4NeutronCaptureProcess.hh"
0044 #include "G4NeutronFissionProcess.hh"
0045 #include "G4ParticleHPCapture.hh"
0046 #include "G4ParticleHPCaptureData.hh"
0047 #include "G4ParticleHPElastic.hh"
0048 #include "G4ParticleHPElasticData.hh"
0049 #include "G4ParticleHPFission.hh"
0050 #include "G4ParticleHPFissionData.hh"
0051 #include "G4ParticleHPInelastic.hh"
0052 #include "G4ParticleHPInelasticData.hh"
0053 #include "G4ParticleHPThermalScattering.hh"
0054 #include "G4ParticleHPThermalScatteringData.hh"
0055 #include "G4SystemOfUnits.hh"
0056
0057 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0058
0059 NeutronHPphysics::NeutronHPphysics(const G4String& name) : G4VPhysicsConstructor(name)
0060 {
0061   // define commands for this class
0062   DefineCommands();
0063 }
0064
0065 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0066
0067 NeutronHPphysics::~NeutronHPphysics()
0068 {
0069   delete fMessenger;
0070 }
0071
0072 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0073
0074 void NeutronHPphysics::ConstructProcess()
0075 {
0076   G4ParticleDefinition* neutron = G4Neutron::Neutron();
0077   G4ProcessManager* pManager = neutron->GetProcessManager();
0078
0079   // delete all neutron processes if already registered
0080   //
0081   G4VProcess* process = 0;
0082   process = pManager->GetProcess("hadElastic");
0083   if (process) pManager->RemoveProcess(process);
0084   //
0085   process = pManager->GetProcess("neutronInelastic");
0086   if (process) pManager->RemoveProcess(process);
0087   //
0088   process = pManager->GetProcess("nCapture");
0089   if (process) pManager->RemoveProcess(process);
0090   //
0091   process = pManager->GetProcess("nFission");
0092   if (process) pManager->RemoveProcess(process);
0093
0094   // (re) create process: elastic
0095   //
0096   G4HadronElasticProcess* process1 = new G4HadronElasticProcess();
0097   pManager->AddDiscreteProcess(process1);
0098   //
0099   // model1a
0100   G4ParticleHPElastic* model1a = new G4ParticleHPElastic();
0101   process1->RegisterMe(model1a);
0102   process1->AddDataSet(new G4ParticleHPElasticData());
0103   //
0104   // model1b
0105   if (fThermal) {
0106     model1a->SetMinEnergy(4 * eV);
0107     G4ParticleHPThermalScattering* model1b = new G4ParticleHPThermalScattering();
0108     process1->RegisterMe(model1b);
0109     process1->AddDataSet(new G4ParticleHPThermalScatteringData());
0110   }
0111
0112   // (re) create process: inelastic
0113   //
0114   G4HadronInelasticProcess* process2 =
0115     new G4HadronInelasticProcess("neutronInelastic", G4Neutron::Definition());
0116   pManager->AddDiscreteProcess(process2);
0117   //
0118   // cross section data set
0119   G4ParticleHPInelasticData* dataSet2 = new G4ParticleHPInelasticData();
0120   process2->AddDataSet(dataSet2);
0121   //
0122   // models
0123   G4ParticleHPInelastic* model2 = new G4ParticleHPInelastic();
0124   process2->RegisterMe(model2);
0125
0126   // (re) create process: nCapture
0127   //
0128   G4NeutronCaptureProcess* process3 = new G4NeutronCaptureProcess();
0129   pManager->AddDiscreteProcess(process3);
0130   //
0131   // cross section data set
0132   G4ParticleHPCaptureData* dataSet3 = new G4ParticleHPCaptureData();
0133   process3->AddDataSet(dataSet3);
0134   //
0135   // models
0136   G4ParticleHPCapture* model3 = new G4ParticleHPCapture();
0137   process3->RegisterMe(model3);
0138
0139   // (re) create process: nFission
0140   //
0141   G4NeutronFissionProcess* process4 = new G4NeutronFissionProcess();
0142   pManager->AddDiscreteProcess(process4);
0143   //
0144   // cross section data set
0145   G4ParticleHPFissionData* dataSet4 = new G4ParticleHPFissionData();
0146   process4->AddDataSet(dataSet4);
0147   //
0148   // models
0149   G4ParticleHPFission* model4 = new G4ParticleHPFission();
0150   process4->RegisterMe(model4);
0151 }
0152
0153 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0154
0155 void NeutronHPphysics::DefineCommands()
0156 {
0157   // Define /testhadr/phys command directory using generic messenger class
0158   fMessenger = new G4GenericMessenger(this, "/testhadr/phys/", "physics list commands");
0159
0160   // thermal scattering command
0161   auto& thermalCmd = fMessenger->DeclareProperty("thermalScattering", fThermal);
0162
0163   thermalCmd.SetGuidance("set thermal scattering model");
0164   thermalCmd.SetParameterName("thermal", false);
0165   thermalCmd.SetStates(G4State_PreInit);
0166 }
0167
0168 //..oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......