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 /// \file PhysicsList.cc 27 /// \brief Implementation of the PhysicsList class 28 // 29 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 32 #include "PhysicsList.hh" 33 34 #include "ElectromagneticPhysics.hh" 35 #include "GammaNuclearPhysics.hh" 36 #include "GammaNuclearPhysicsLEND.hh" 37 #include "HadronElasticPhysicsHP.hh" 38 #include "RadioactiveDecayPhysics.hh" 39 40 #include "G4DecayPhysics.hh" 41 #include "G4EmStandardPhysics_option3.hh" 42 #include "G4HadronElasticPhysicsXS.hh" 43 #include "G4HadronInelasticQBBC.hh" 44 #include "G4HadronPhysicsFTFP_BERT_HP.hh" 45 #include "G4HadronPhysicsINCLXX.hh" 46 #include "G4HadronPhysicsQGSP_BIC_AllHP.hh" 47 #include "G4HadronPhysicsQGSP_BIC_HP.hh" 48 #include "G4HadronicInteraction.hh" 49 #include "G4IonElasticPhysics.hh" 50 #include "G4IonINCLXXPhysics.hh" 51 #include "G4IonPhysicsPHP.hh" 52 #include "G4IonPhysicsXS.hh" 53 #include "G4Neutron.hh" 54 #include "G4NuclideTable.hh" 55 #include "G4ProcessManager.hh" 56 #include "G4RadioactiveDecayPhysics.hh" 57 #include "G4StoppingPhysics.hh" 58 #include "G4SystemOfUnits.hh" 59 #include "G4UnitsTable.hh" 60 61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 62 63 PhysicsList::PhysicsList() 64 { 65 G4int verb = 0; 66 SetVerboseLevel(verb); 67 68 // add new units 69 // 70 new G4UnitDefinition("mm2/g", "mm2/g", "Surface/Mass", mm2 / g); 71 new G4UnitDefinition("um2/mg", "um2/mg", "Surface/Mass", um * um / mg); 72 73 // mandatory for G4NuclideTable 74 // 75 const G4double meanLife = 1 * nanosecond, halfLife = meanLife * std::log(2); 76 G4NuclideTable::GetInstance()->SetThresholdOfHalfLife(halfLife); 77 78 // Hadron Elastic scattering 79 fHadronElastic = new HadronElasticPhysicsHP(verb); 80 ////fHadronElastic = new G4HadronElasticPhysicsXS(verb); 81 RegisterPhysics(fHadronElastic); 82 83 // Hadron Inelastic Physics 84 ////fHadronInelastic = new G4HadronPhysicsFTFP_BERT_HP(verb); 85 fHadronInelastic = new G4HadronPhysicsQGSP_BIC_HP(verb); 86 ////fHadronInelastic = new G4HadronPhysicsQGSP_BIC_AllHP(verb); 87 ////fHadronInelastic = new G4HadronInelasticQBBC(verb); 88 ////fHadronInelastic = new G4HadronPhysicsINCLXX(verb); 89 RegisterPhysics(fHadronInelastic); 90 91 // Ion Elastic Physics 92 fIonElastic = new G4IonElasticPhysics(verb); 93 RegisterPhysics(fIonElastic); 94 95 // Ion Inelastic Physics 96 fIonInelastic = new G4IonPhysicsXS(verb); 97 ////fIonInelastic = new G4IonPhysicsPHP(verb) 98 ////fIonInelastic = new G4IonINCLXXPhysics(verb); 99 RegisterPhysics(fIonInelastic); 100 101 // stopping Particles 102 /// RegisterPhysics( new G4StoppingPhysics(verb)); 103 104 // Gamma-Nuclear Physics 105 fGammaNuclear = new GammaNuclearPhysics("gamma"); 106 ////fGammaNuclear = new GammaNuclearPhysicsLEND("gamma"); 107 RegisterPhysics(fGammaNuclear); 108 109 // EM physics 110 fElectromagnetic = new ElectromagneticPhysics(); 111 ////fElectromagnetic = new G4EmStandardPhysics_option3(); 112 RegisterPhysics(fElectromagnetic); 113 114 // Decay 115 fDecay = new G4DecayPhysics(); 116 RegisterPhysics(fDecay); 117 118 // Radioactive decay 119 fRadioactiveDecay = new RadioactiveDecayPhysics(); 120 ////fRadioactiveDecay = new G4RadioactiveDecayPhysics(); 121 RegisterPhysics(fRadioactiveDecay); 122 } 123 124 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 125 126 void PhysicsList::ConstructProcess() 127 { 128 // Transportation first (mandatory) 129 // 130 AddTransportation(); 131 132 // Physics constructors 133 // 134 fHadronElastic->ConstructProcess(); 135 fHadronInelastic->ConstructProcess(); 136 fIonElastic->ConstructProcess(); 137 fIonInelastic->ConstructProcess(); 138 fGammaNuclear->ConstructProcess(); 139 fElectromagnetic->ConstructProcess(); 140 fDecay->ConstructProcess(); 141 fRadioactiveDecay->ConstructProcess(); 142 143 // example of GetHadronicModel (due to bug in QGSP_BIC_AllHP) 144 // 145 G4ProcessManager* pManager = G4Neutron::Neutron()->GetProcessManager(); 146 G4HadronicProcess* process = dynamic_cast<G4HadronicProcess*>(pManager->GetProcess("nCapture")); 147 G4HadronicInteraction* model = process->GetHadronicModel("nRadCapture"); 148 if (model) model->SetMinEnergy(19.9 * MeV); 149 } 150 151 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 152 153 void PhysicsList::SetCuts() 154 { 155 SetCutValue(0 * mm, "proton"); 156 SetCutValue(10 * km, "e-"); 157 SetCutValue(10 * km, "e+"); 158 SetCutValue(10 * km, "gamma"); 159 } 160 161 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 162