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