Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 //-------------------------------------------- 28 // 29 // ClassName: INCLXXPhysicsListHelper 30 // 31 // Author: 04.03.2013 Davide Mancusi 32 // Created by modifying and generalising the 33 // 34 //-------------------------------------------- 35 // 36 37 #include <iomanip> 38 39 #include "globals.hh" 40 #include "G4ios.hh" 41 #include "G4ProcessManager.hh" 42 #include "G4ProcessVector.hh" 43 #include "G4ParticleTypes.hh" 44 #include "G4ParticleTable.hh" 45 46 #include "G4Material.hh" 47 #include "G4MaterialTable.hh" 48 49 #include "G4DecayPhysics.hh" 50 #include "G4RadioactiveDecayPhysics.hh" 51 #include "G4EmStandardPhysics.hh" 52 #include "G4EmExtraPhysics.hh" 53 #include "G4IonINCLXXPhysics.hh" 54 #include "G4StoppingPhysics.hh" 55 #include "G4StoppingPhysicsWithINCLXX.hh" 56 #include "G4HadronElasticPhysics.hh" 57 #include "G4HadronElasticPhysicsHP.hh" 58 #include "G4NeutronTrackingCut.hh" 59 60 #include "G4HadronPhysicsINCLXX.hh" 61 62 #include "G4WarnPLStatus.hh" 63 64 template<class T, bool withNeutronHP, bool wit 65 { 66 67 if(withFTFP) 68 name = "FTFP_INCLXX"; 69 else 70 name = "QGSP_INCLXX"; 71 if(withNeutronHP) 72 name += "_HP"; 73 74 G4cout << "<<< Geant4 Physics List simulatio 75 G4cout <<G4endl; 76 77 this->defaultCutValue = 0.7*CLHEP::mm; 78 this->SetVerboseLevel(ver); 79 80 G4WarnPLStatus exp; 81 exp.Experimental(name); 82 83 // EM Physics 84 this->RegisterPhysics( new G4EmStandardPhysi 85 86 // Synchroton Radiation & GN Physics 87 this->RegisterPhysics( new G4EmExtraPhysics( 88 89 // Decays 90 this->RegisterPhysics( new G4DecayPhysics(ve 91 if(withNeutronHP) this->RegisterPhysics( new 92 93 // Hadron Elastic scattering 94 if(withNeutronHP) 95 this->RegisterPhysics( new G4HadronElastic 96 else 97 this->RegisterPhysics( new G4HadronElastic 98 99 // Hadron Physics 100 this->RegisterPhysics( new G4HadronPhysicsIN 101 102 // Stopping Physics 103 this->RegisterPhysics( new G4StoppingPhysics 104 105 // Ion Physics 106 this->RegisterPhysics( new G4IonINCLXXPhysic 107 108 // Neutron tracking cut 109 if(!withNeutronHP) 110 this->RegisterPhysics( new G4NeutronTracki 111 112 } 113 114 template<class T, bool withNeutronHP, bool wit 115 { 116 } 117 118 template<class T, bool withNeutronHP, bool wit 119 { 120 if (this->verboseLevel >1){ 121 G4cout << name << "::SetCuts:"; 122 } 123 // " G4VUserPhysicsList::SetCutsWithDefault 124 // the default cut value for all particle 125 126 this->SetCutsWithDefault(); 127 128 if(withNeutronHP) { this->SetCutValue(0.0, " 129 130 // if (this->verboseLevel >0) 131 // G4VUserPhysicsList::DumpCutValuesTable 132 } 133 134