Geant4 Cross Reference |
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 /// \file hadronic/Hadr01/src/PhysicsList.cc 26 /// \file hadronic/Hadr01/src/PhysicsList.cc 27 /// \brief Implementation of the PhysicsList c 27 /// \brief Implementation of the PhysicsList class 28 // 28 // 29 // 29 // >> 30 // $Id: PhysicsList.cc 92421 2015-09-01 07:38:57Z gcosmo $ 30 // 31 // 31 ////////////////////////////////////////////// 32 ///////////////////////////////////////////////////////////////////////// 32 // 33 // 33 // PhysicsList 34 // PhysicsList 34 // 35 // 35 // Created: 31.04.2006 V.Ivanchenko 36 // Created: 31.04.2006 V.Ivanchenko 36 // 37 // 37 // Modified: 38 // Modified: 38 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchen 39 // 04.06.2006 Adoptation of Hadr01 (V.Ivanchenko) 39 // 26.04.2007 Physics according to 8.3 Physics 40 // 26.04.2007 Physics according to 8.3 Physics List (V.Ivanchenko) 40 // 16.10.2012 Renamed used classes (A.Ribon) 41 // 16.10.2012 Renamed used classes (A.Ribon) 41 // 42 // 42 ////////////////////////////////////////////// 43 //////////////////////////////////////////////////////////////////////// 43 // << 44 // 44 45 45 #include "PhysicsList.hh" 46 #include "PhysicsList.hh" 46 << 47 #include "PhysicsListMessenger.hh" 47 #include "PhysicsListMessenger.hh" 48 48 49 #include "G4DecayPhysics.hh" 49 #include "G4DecayPhysics.hh" 50 #include "G4Electron.hh" << 51 #include "G4EmExtraPhysics.hh" << 52 #include "G4EmLivermorePhysics.hh" << 53 #include "G4EmParameters.hh" << 54 #include "G4EmPenelopePhysics.hh" << 55 #include "G4EmStandardPhysics.hh" 50 #include "G4EmStandardPhysics.hh" 56 #include "G4EmStandardPhysicsGS.hh" << 57 #include "G4EmStandardPhysics_option1.hh" 51 #include "G4EmStandardPhysics_option1.hh" 58 #include "G4EmStandardPhysics_option2.hh" 52 #include "G4EmStandardPhysics_option2.hh" 59 #include "G4EmStandardPhysics_option3.hh" 53 #include "G4EmStandardPhysics_option3.hh" 60 #include "G4EmStandardPhysics_option4.hh" 54 #include "G4EmStandardPhysics_option4.hh" 61 #include "G4Gamma.hh" << 55 #include "G4EmStandardPhysicsGS.hh" >> 56 #include "G4EmLivermorePhysics.hh" >> 57 #include "G4EmPenelopePhysics.hh" 62 #include "G4HadronElasticPhysics.hh" 58 #include "G4HadronElasticPhysics.hh" 63 #include "G4HadronElasticPhysicsHP.hh" << 64 #include "G4HadronElasticPhysicsXS.hh" 59 #include "G4HadronElasticPhysicsXS.hh" >> 60 #include "G4HadronElasticPhysicsHP.hh" 65 #include "G4HadronHElasticPhysics.hh" 61 #include "G4HadronHElasticPhysics.hh" 66 #include "G4HadronInelasticQBBC.hh" << 62 #include "G4ChargeExchangePhysics.hh" >> 63 #include "G4NeutronTrackingCut.hh" >> 64 #include "G4NeutronCrossSectionXS.hh" >> 65 #include "G4StoppingPhysics.hh" >> 66 #include "G4IonBinaryCascadePhysics.hh" >> 67 #include "G4IonPhysics.hh" >> 68 #include "G4EmExtraPhysics.hh" >> 69 #include "G4EmProcessOptions.hh" >> 70 67 #include "G4HadronPhysicsFTFP_BERT.hh" 71 #include "G4HadronPhysicsFTFP_BERT.hh" 68 #include "G4HadronPhysicsFTFP_BERT_HP.hh" 72 #include "G4HadronPhysicsFTFP_BERT_HP.hh" 69 #include "G4HadronPhysicsFTFP_BERT_TRV.hh" 73 #include "G4HadronPhysicsFTFP_BERT_TRV.hh" 70 #include "G4HadronPhysicsFTF_BIC.hh" 74 #include "G4HadronPhysicsFTF_BIC.hh" >> 75 #include "G4HadronInelasticQBBC.hh" 71 #include "G4HadronPhysicsQGSP_BERT.hh" 76 #include "G4HadronPhysicsQGSP_BERT.hh" 72 #include "G4HadronPhysicsQGSP_BERT_HP.hh" 77 #include "G4HadronPhysicsQGSP_BERT_HP.hh" 73 #include "G4HadronPhysicsQGSP_BIC.hh" 78 #include "G4HadronPhysicsQGSP_BIC.hh" 74 #include "G4HadronPhysicsQGSP_BIC_HP.hh" 79 #include "G4HadronPhysicsQGSP_BIC_HP.hh" 75 #include "G4HadronPhysicsQGSP_FTFP_BERT.hh" 80 #include "G4HadronPhysicsQGSP_FTFP_BERT.hh" 76 #include "G4HadronPhysicsQGS_BIC.hh" 81 #include "G4HadronPhysicsQGS_BIC.hh" 77 #include "G4IonBinaryCascadePhysics.hh" << 82 78 #include "G4IonElasticPhysics.hh" << 83 #include "G4LossTableManager.hh" 79 #include "G4IonPhysics.hh" << 84 80 #include "G4IonPhysicsXS.hh" << 85 #include "G4ProcessManager.hh" 81 #include "G4NeutronCrossSectionXS.hh" << 82 #include "G4NeutronTrackingCut.hh" << 83 #include "G4ParticleTable.hh" << 84 #include "G4ParticleTypes.hh" 86 #include "G4ParticleTypes.hh" 85 #include "G4PhysListFactoryMessenger.hh" << 87 #include "G4ParticleTable.hh" >> 88 #include "G4Gamma.hh" >> 89 #include "G4Electron.hh" 86 #include "G4Positron.hh" 90 #include "G4Positron.hh" 87 #include "G4ProcessManager.hh" << 88 #include "G4Proton.hh" 91 #include "G4Proton.hh" 89 #include "G4RadioactiveDecayPhysics.hh" << 92 90 #include "G4StoppingPhysics.hh" << 91 #include "G4SystemOfUnits.hh" 93 #include "G4SystemOfUnits.hh" 92 94 93 //....oooOO0OOooo........oooOO0OOooo........oo 95 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 94 96 95 PhysicsList::PhysicsList() : G4VModularPhysics << 97 PhysicsList::PhysicsList() >> 98 : G4VModularPhysicsList(), >> 99 fEmPhysicsList(0), fParticleList(0), fMessenger(0) 96 { 100 { 97 SetDefaultCutValue(0.7 * CLHEP::mm); << 101 G4LossTableManager::Instance(); 98 verboseLevel = 1; << 102 defaultCutValue = 0.7*mm; >> 103 fCutForGamma = defaultCutValue; >> 104 fCutForElectron = defaultCutValue; >> 105 fCutForPositron = defaultCutValue; >> 106 fCutForProton = defaultCutValue; >> 107 verboseLevel = 1; 99 108 100 // Messengers << 101 fMessenger = new PhysicsListMessenger(this); 109 fMessenger = new PhysicsListMessenger(this); 102 fFactMessenger = new G4PhysListFactoryMessen << 103 110 104 // Particles 111 // Particles 105 fParticleList = new G4DecayPhysics(verboseLe << 112 fParticleList = new G4DecayPhysics("decays"); 106 113 107 // EM physics 114 // EM physics 108 fEmPhysicsList = new G4EmStandardPhysics(ver << 115 fEmPhysicsList = new G4EmStandardPhysics(); 109 } 116 } 110 117 111 //....oooOO0OOooo........oooOO0OOooo........oo 118 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 112 119 113 PhysicsList::~PhysicsList() 120 PhysicsList::~PhysicsList() 114 { 121 { 115 delete fMessenger; 122 delete fMessenger; 116 delete fParticleList; 123 delete fParticleList; 117 delete fEmPhysicsList; 124 delete fEmPhysicsList; 118 for (std::size_t i = 0; i < fHadronPhys.size << 125 for(size_t i=0; i<fHadronPhys.size(); i++) { 119 delete fHadronPhys[i]; 126 delete fHadronPhys[i]; 120 } 127 } 121 } 128 } 122 129 123 //....oooOO0OOooo........oooOO0OOooo........oo 130 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 124 131 125 void PhysicsList::ConstructParticle() 132 void PhysicsList::ConstructParticle() 126 { 133 { 127 fParticleList->ConstructParticle(); 134 fParticleList->ConstructParticle(); 128 } 135 } 129 136 130 //....oooOO0OOooo........oooOO0OOooo........oo 137 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 131 138 132 void PhysicsList::ConstructProcess() 139 void PhysicsList::ConstructProcess() 133 { 140 { 134 AddTransportation(); 141 AddTransportation(); 135 fEmPhysicsList->ConstructProcess(); 142 fEmPhysicsList->ConstructProcess(); 136 fParticleList->ConstructProcess(); 143 fParticleList->ConstructProcess(); 137 for (std::size_t i = 0; i < fHadronPhys.size << 144 for(size_t i=0; i<fHadronPhys.size(); i++) { 138 fHadronPhys[i]->ConstructProcess(); 145 fHadronPhys[i]->ConstructProcess(); 139 } 146 } 140 } 147 } 141 148 142 //....oooOO0OOooo........oooOO0OOooo........oo 149 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 143 150 144 void PhysicsList::AddPhysicsList(const G4Strin 151 void PhysicsList::AddPhysicsList(const G4String& name) 145 { 152 { 146 if (verboseLevel > 0) { << 153 if (verboseLevel>0) { 147 G4cout << "PhysicsList::AddPhysicsList: <" 154 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl; 148 } 155 } 149 if (name == "emstandard_opt0") { 156 if (name == "emstandard_opt0") { >> 157 150 delete fEmPhysicsList; 158 delete fEmPhysicsList; 151 fEmPhysicsList = new G4EmStandardPhysics(v << 159 fEmPhysicsList = new G4EmStandardPhysics(); 152 } << 160 153 else if (name == "emstandard_opt1") { << 161 } else if (name == "emstandard_opt1") { >> 162 154 delete fEmPhysicsList; 163 delete fEmPhysicsList; 155 fEmPhysicsList = new G4EmStandardPhysics_o << 164 fEmPhysicsList = new G4EmStandardPhysics_option1(); 156 } << 165 157 else if (name == "emstandard_opt2") { << 166 } else if (name == "emstandard_opt2") { >> 167 158 delete fEmPhysicsList; 168 delete fEmPhysicsList; 159 fEmPhysicsList = new G4EmStandardPhysics_o << 169 fEmPhysicsList = new G4EmStandardPhysics_option2(); 160 } << 170 161 else if (name == "emstandard_opt3") { << 171 } else if (name == "emstandard_opt3") { >> 172 162 delete fEmPhysicsList; 173 delete fEmPhysicsList; 163 fEmPhysicsList = new G4EmStandardPhysics_o << 174 fEmPhysicsList = new G4EmStandardPhysics_option3(); 164 } << 175 165 else if (name == "emstandard_opt4") { << 176 } else if (name == "emstandard_opt4") { >> 177 166 delete fEmPhysicsList; 178 delete fEmPhysicsList; 167 fEmPhysicsList = new G4EmStandardPhysics_o << 179 fEmPhysicsList = new G4EmStandardPhysics_option4(); 168 } << 180 169 else if (name == "emstandardGS") { << 181 } else if (name == "emstandardGS") { >> 182 170 delete fEmPhysicsList; 183 delete fEmPhysicsList; 171 fEmPhysicsList = new G4EmStandardPhysicsGS << 184 fEmPhysicsList = new G4EmStandardPhysicsGS(); 172 } << 185 173 else if (name == "FTFP_BERT_EMV") { << 186 } else if (name == "FTFP_BERT_EMV") { 174 AddPhysicsList("FTFP_BERT"); << 187 175 AddPhysicsList("emstandard_opt1"); 188 AddPhysicsList("emstandard_opt1"); 176 } << 177 else if (name == "FTFP_BERT_EMX") { << 178 AddPhysicsList("FTFP_BERT"); 189 AddPhysicsList("FTFP_BERT"); >> 190 >> 191 } else if (name == "FTFP_BERT_EMX") { >> 192 179 AddPhysicsList("emstandard_opt2"); 193 AddPhysicsList("emstandard_opt2"); 180 } << 181 else if (name == "FTFP_BERT_EMY") { << 182 AddPhysicsList("FTFP_BERT"); 194 AddPhysicsList("FTFP_BERT"); >> 195 >> 196 } else if (name == "FTFP_BERT_EMY") { >> 197 183 AddPhysicsList("emstandard_opt3"); 198 AddPhysicsList("emstandard_opt3"); 184 } << 185 else if (name == "FTFP_BERT_EMZ") { << 186 AddPhysicsList("FTFP_BERT"); 199 AddPhysicsList("FTFP_BERT"); >> 200 >> 201 } else if (name == "FTFP_BERT_EMZ") { >> 202 187 AddPhysicsList("emstandard_opt4"); 203 AddPhysicsList("emstandard_opt4"); 188 } << 204 AddPhysicsList("FTFP_BERT"); 189 else if (name == "FTFP_BERT") { << 205 >> 206 } else if (name == "FTFP_BERT") { >> 207 190 SetBuilderList0(false); 208 SetBuilderList0(false); 191 fHadronPhys.push_back(new G4HadronPhysicsF << 209 fHadronPhys.push_back( new G4HadronPhysicsFTFP_BERT()); 192 } << 210 193 else if (name == "FTFP_BERT_TRV") { << 211 } else if (name == "FTFP_BERT_TRV") { >> 212 194 AddPhysicsList("emstandardGS"); 213 AddPhysicsList("emstandardGS"); 195 SetBuilderList1(false); 214 SetBuilderList1(false); 196 fHadronPhys.push_back(new G4HadronPhysicsF << 215 fHadronPhys.push_back( new G4HadronPhysicsFTFP_BERT_TRV()); 197 } << 216 198 else if (name == "FTF_BIC") { << 217 } else if (name == "FTF_BIC") { >> 218 199 SetBuilderList0(false); 219 SetBuilderList0(false); 200 fHadronPhys.push_back(new G4HadronPhysicsF << 220 fHadronPhys.push_back( new G4HadronPhysicsFTF_BIC()); 201 } << 221 202 else if (name == "QBBC") { << 222 } else if (name == "QBBC") { >> 223 203 AddPhysicsList("emstandard_opt0"); 224 AddPhysicsList("emstandard_opt0"); 204 SetBuilderList2(); 225 SetBuilderList2(); 205 fHadronPhys.push_back(new G4HadronInelasti << 226 fHadronPhys.push_back( new G4HadronInelasticQBBC()); 206 } << 227 207 else if (name == "QGSP_BERT") { << 228 } else if (name == "QGSP_BERT") { >> 229 208 SetBuilderList0(false); 230 SetBuilderList0(false); 209 fHadronPhys.push_back(new G4HadronPhysicsQ << 231 fHadronPhys.push_back( new G4HadronPhysicsQGSP_BERT()); 210 } << 232 211 else if (name == "QGSP_FTFP_BERT") { << 233 } else if (name == "QGSP_FTFP_BERT") { >> 234 212 SetBuilderList0(false); 235 SetBuilderList0(false); 213 fHadronPhys.push_back(new G4HadronPhysicsQ << 236 fHadronPhys.push_back( new G4HadronPhysicsQGSP_FTFP_BERT()); 214 } << 237 215 else if (name == "QGSP_FTFP_BERT_EMV") { << 238 } else if (name == "QGSP_FTFP_BERT_EMV") { 216 AddPhysicsList("QGSP_FTFP_BERT"); << 239 217 AddPhysicsList("emstandard_opt1"); 240 AddPhysicsList("emstandard_opt1"); 218 } << 241 AddPhysicsList("QGSP_FTFP_BERT"); 219 else if (name == "QGSP_BERT_EMV") { << 242 220 AddPhysicsList("QGSP_BERT"); << 243 } else if (name == "QGSP_BERT_EMV") { >> 244 221 AddPhysicsList("emstandard_opt1"); 245 AddPhysicsList("emstandard_opt1"); 222 } << 223 else if (name == "QGSP_BERT_EMX") { << 224 AddPhysicsList("QGSP_BERT"); 246 AddPhysicsList("QGSP_BERT"); >> 247 >> 248 } else if (name == "QGSP_BERT_EMX") { >> 249 225 AddPhysicsList("emstandard_opt2"); 250 AddPhysicsList("emstandard_opt2"); 226 } << 251 AddPhysicsList("QGSP_BERT"); 227 else if (name == "QGSP_BERT_HP") { << 252 >> 253 } else if (name == "QGSP_BERT_HP") { >> 254 228 SetBuilderList0(true); 255 SetBuilderList0(true); 229 fHadronPhys.push_back(new G4HadronPhysicsQ << 256 fHadronPhys.push_back( new G4HadronPhysicsQGSP_BERT_HP()); 230 AddPhysicsList("RadioactiveDecay"); << 257 231 } << 258 } else if (name == "QGSP_BIC") { 232 else if (name == "QGSP_BIC") { << 259 233 SetBuilderList0(false); 260 SetBuilderList0(false); 234 fHadronPhys.push_back(new G4HadronPhysicsQ << 261 fHadronPhys.push_back( new G4HadronPhysicsQGSP_BIC()); 235 fHadronPhys.push_back(new G4IonElasticPhys << 262 236 } << 263 } else if (name == "QGSP_BIC_EMY") { 237 else if (name == "QGSP_BIC_EMY") { << 264 238 AddPhysicsList("QGSP_BIC"); << 239 AddPhysicsList("emstandard_opt3"); 265 AddPhysicsList("emstandard_opt3"); 240 } << 241 else if (name == "QGS_BIC") { << 242 SetBuilderList0(false); 266 SetBuilderList0(false); 243 fHadronPhys.push_back(new G4HadronPhysicsQ << 267 fHadronPhys.push_back( new G4HadronPhysicsQGSP_BIC()); 244 } << 268 245 else if (name == "QGSP_BIC_HP") { << 269 } else if (name == "QGS_BIC") { >> 270 >> 271 SetBuilderList0(false); >> 272 fHadronPhys.push_back( new G4HadronPhysicsQGS_BIC()); >> 273 >> 274 } else if (name == "QGSP_BIC_HP") { >> 275 246 SetBuilderList0(true); 276 SetBuilderList0(true); 247 fHadronPhys.push_back(new G4HadronPhysicsQ << 277 fHadronPhys.push_back( new G4HadronPhysicsQGSP_BIC_HP()); 248 AddPhysicsList("RadioactiveDecay"); << 278 249 } << 279 } else { 250 else if (name == "RadioactiveDecay") { << 280 251 fHadronPhys.push_back(new G4RadioactiveDec << 252 } << 253 else { << 254 G4cout << "PhysicsList::AddPhysicsList: <" 281 G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" 255 << " is not defined" << G4endl; << 282 << " is not defined" >> 283 << G4endl; 256 } 284 } 257 } 285 } 258 286 259 //....oooOO0OOooo........oooOO0OOooo........oo 287 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 260 288 261 void PhysicsList::SetBuilderList0(G4bool flagH 289 void PhysicsList::SetBuilderList0(G4bool flagHP) 262 { 290 { 263 fHadronPhys.push_back(new G4EmExtraPhysics(v << 291 fHadronPhys.push_back( new G4EmExtraPhysics(verboseLevel)); 264 if (flagHP) { << 292 if(flagHP) { 265 fHadronPhys.push_back(new G4HadronElasticP << 293 fHadronPhys.push_back( new G4HadronElasticPhysicsHP(verboseLevel) ); 266 } << 294 } else { 267 else { << 295 fHadronPhys.push_back( new G4HadronElasticPhysics(verboseLevel) ); 268 fHadronPhys.push_back(new G4HadronElasticP << 296 } 269 } << 297 fHadronPhys.push_back( new G4StoppingPhysics(verboseLevel)); 270 fHadronPhys.push_back(new G4StoppingPhysics( << 298 fHadronPhys.push_back( new G4IonPhysics(verboseLevel)); 271 fHadronPhys.push_back(new G4IonPhysics(verbo << 299 fHadronPhys.push_back( new G4NeutronTrackingCut(verboseLevel)); 272 if (!flagHP) { << 273 fHadronPhys.push_back(new G4NeutronTrackin << 274 } << 275 } 300 } 276 301 277 //....oooOO0OOooo........oooOO0OOooo........oo 302 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 278 303 279 void PhysicsList::SetBuilderList1(G4bool flagH 304 void PhysicsList::SetBuilderList1(G4bool flagHP) 280 { 305 { 281 fHadronPhys.push_back(new G4EmExtraPhysics(v << 306 fHadronPhys.push_back( new G4EmExtraPhysics(verboseLevel)); 282 if (flagHP) { << 307 if(flagHP) { 283 fHadronPhys.push_back(new G4HadronElasticP << 308 fHadronPhys.push_back( new G4HadronElasticPhysicsHP(verboseLevel) ); 284 } << 309 } else { 285 else { << 310 fHadronPhys.push_back( new G4HadronHElasticPhysics(verboseLevel) ); 286 fHadronPhys.push_back(new G4HadronHElastic << 311 } 287 } << 312 fHadronPhys.push_back( new G4StoppingPhysics(verboseLevel)); 288 fHadronPhys.push_back(new G4StoppingPhysics( << 313 fHadronPhys.push_back( new G4IonPhysics(verboseLevel)); 289 fHadronPhys.push_back(new G4IonPhysics(verbo << 314 fHadronPhys.push_back( new G4NeutronTrackingCut(verboseLevel)); 290 if (!flagHP) { << 291 fHadronPhys.push_back(new G4NeutronTrackin << 292 } << 293 } 315 } 294 316 295 //....oooOO0OOooo........oooOO0OOooo........oo 317 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 296 318 297 void PhysicsList::SetBuilderList2() 319 void PhysicsList::SetBuilderList2() 298 { 320 { 299 fHadronPhys.push_back(new G4EmExtraPhysics(v << 321 fHadronPhys.push_back( new G4EmExtraPhysics(verboseLevel)); 300 fHadronPhys.push_back(new G4HadronElasticPhy << 322 fHadronPhys.push_back( new G4HadronElasticPhysicsXS(verboseLevel) ); 301 fHadronPhys.push_back(new G4StoppingPhysics( << 323 fHadronPhys.push_back( new G4StoppingPhysics(verboseLevel)); 302 fHadronPhys.push_back(new G4IonPhysicsXS(ver << 324 fHadronPhys.push_back( new G4IonPhysics(verboseLevel)); 303 fHadronPhys.push_back(new G4IonElasticPhysic << 325 fHadronPhys.push_back( new G4NeutronTrackingCut(verboseLevel)); 304 fHadronPhys.push_back(new G4NeutronTrackingC << 305 } 326 } 306 327 307 //....oooOO0OOooo........oooOO0OOooo........oo 328 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 308 329 >> 330 void PhysicsList::SetCuts() >> 331 { >> 332 >> 333 if (verboseLevel >0){ >> 334 G4cout << "PhysicsList::SetCuts:"; >> 335 G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl; >> 336 } >> 337 >> 338 // set cut values for gamma at first and for e- second and next for e+, >> 339 // because some processes for e+/e- need cut values for gamma >> 340 SetCutValue(fCutForGamma, "gamma"); >> 341 SetCutValue(fCutForElectron, "e-"); >> 342 SetCutValue(fCutForPositron, "e+"); >> 343 SetCutValue(fCutForProton, "proton"); >> 344 >> 345 if (verboseLevel>0) { DumpCutValuesTable(); } >> 346 } >> 347 >> 348 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... >> 349 >> 350 void PhysicsList::SetCutForGamma(G4double cut) >> 351 { >> 352 fCutForGamma = cut; >> 353 SetParticleCuts(fCutForGamma, G4Gamma::Gamma()); >> 354 } >> 355 >> 356 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 357 >> 358 void PhysicsList::SetCutForElectron(G4double cut) >> 359 { >> 360 fCutForElectron = cut; >> 361 SetParticleCuts(fCutForElectron, G4Electron::Electron()); >> 362 } >> 363 >> 364 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 365 >> 366 void PhysicsList::SetCutForPositron(G4double cut) >> 367 { >> 368 fCutForPositron = cut; >> 369 SetParticleCuts(fCutForPositron, G4Positron::Positron()); >> 370 } >> 371 >> 372 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 373 >> 374 void PhysicsList::SetCutForProton(G4double cut) >> 375 { >> 376 fCutForProton = cut; >> 377 SetParticleCuts(fCutForProton, G4Proton::Proton()); >> 378 } >> 379 309 void PhysicsList::List() 380 void PhysicsList::List() 310 { 381 { 311 G4cout << "### PhysicsLists available: FTFP_ 382 G4cout << "### PhysicsLists available: FTFP_BERT FTFP_BERT_EMV " 312 << "FTFP_BERT_EMX FTFP_BERT_EMZ FTFP_ << 383 << "FTFP_BERT_EMX FTFP_BERT_EMZ FTFP_BERT_TRV" >> 384 << G4endl; 313 G4cout << " FTF_B 385 G4cout << " FTF_BIC QBBC QGSP_BERT " 314 << "QGSP_BERT_EMV QGSP_BERT_EMX" << G << 386 << "QGSP_BERT_EMV QGSP_BERT_EMX" >> 387 << G4endl; 315 G4cout << " QGSP_ 388 G4cout << " QGSP_BERT_HP QGSP_FTFP_BERT " 316 << "QGSP_FTFP_BERT_EMV" << G4endl; << 389 << "QGSP_FTFP_BERT_EMV" >> 390 << G4endl; 317 G4cout << " QGS_B 391 G4cout << " QGS_BIC QGSP_BIC QGSP_BIC_EMY " 318 << "QGSP_BIC_HP" << G4endl; << 392 << "QGSP_BIC_HP" >> 393 << G4endl; 319 } 394 } 320 395 321 //....oooOO0OOooo........oooOO0OOooo........oo 396 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 397 322 398