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