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/Hadr02/src/IonUrQMDPhysics. 26 /// \file hadronic/Hadr02/src/IonUrQMDPhysics.cc 27 /// \brief Implementation of the IonUrQMDPhysi 27 /// \brief Implementation of the IonUrQMDPhysics class 28 // 28 // >> 29 // $Id: IonUrQMDPhysics.cc 77519 2013-11-25 10:54:57Z gcosmo $ 29 // 30 // 30 //-------------------------------------------- 31 //--------------------------------------------------------------------------- 31 // 32 // 32 // Class: IonUrQMDPhysics 33 // Class: IonUrQMDPhysics 33 // 34 // 34 // Author: 2012 Andrea Dotti 35 // Author: 2012 Andrea Dotti 35 // 36 // 36 // 37 // 37 // Modified: 38 // Modified: 38 // 39 // 39 // ------------------------------------------- 40 // ------------------------------------------------------------ 40 // << 41 // 41 #ifdef G4_USE_URQMD 42 #ifdef G4_USE_URQMD 42 # include "IonUrQMDPhysics.hh" << 43 #include "IonUrQMDPhysics.hh" >> 44 #include "G4ParticleDefinition.hh" >> 45 #include "G4ProcessManager.hh" >> 46 #include "G4Deuteron.hh" >> 47 #include "G4Triton.hh" >> 48 #include "G4He3.hh" >> 49 #include "G4Alpha.hh" >> 50 #include "G4GenericIon.hh" >> 51 >> 52 #include "G4HadronInelasticProcess.hh" >> 53 #include "G4TripathiCrossSection.hh" >> 54 #include "G4TripathiLightCrossSection.hh" >> 55 #include "G4IonsShenCrossSection.hh" >> 56 #include "G4IonProtonCrossSection.hh" >> 57 >> 58 #include "G4UrQMD1_3Model.hh" 43 59 44 # include "G4Alpha.hh" << 60 #include "G4BuilderType.hh" 45 # include "G4BuilderType.hh" << 61 #include "G4SystemOfUnits.hh" 46 # include "G4ComponentGGNuclNuclXsc.hh" << 47 # include "G4CrossSectionInelastic.hh" << 48 # include "G4Deuteron.hh" << 49 # include "G4GenericIon.hh" << 50 # include "G4HadronInelasticProcess.hh" << 51 # include "G4HadronicParameters.hh" << 52 # include "G4He3.hh" << 53 # include "G4ParticleDefinition.hh" << 54 # include "G4ProcessManager.hh" << 55 # include "G4SystemOfUnits.hh" << 56 # include "G4Triton.hh" << 57 # include "G4UrQMD1_3Model.hh" << 58 62 59 using namespace std; 63 using namespace std; 60 64 61 //....oooOO0OOooo........oooOO0OOooo........oo 65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 62 66 63 IonUrQMDPhysics::IonUrQMDPhysics(G4int ver) 67 IonUrQMDPhysics::IonUrQMDPhysics(G4int ver) 64 : G4VHadronPhysics("ionInelasticUrQMD"), ver << 68 : G4VHadronPhysics("ionInelasticUrQMD"),verbose(ver), >> 69 fWasActivated(false) 65 { 70 { 66 fIonXS = nullptr; << 71 fTripathi = fTripathiLight = fShen = fIonH = 0; 67 fModel = nullptr; << 72 fModel = 0; 68 SetPhysicsType(bIons); 73 SetPhysicsType(bIons); 69 if (fVerbose > 1) { << 74 if(fVerbose > 1) { G4cout << "### IonUrQMDPhysics" << G4endl; } 70 G4cout << "### IonUrQMDPhysics" << G4endl; << 71 } << 72 } 75 } 73 76 74 //....oooOO0OOooo........oooOO0OOooo........oo 77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 75 78 76 IonUrQMDPhysics::~IonUrQMDPhysics() {} << 79 IonUrQMDPhysics::~IonUrQMDPhysics() >> 80 {} 77 81 78 //....oooOO0OOooo........oooOO0OOooo........oo 82 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 79 void IonUrQMDPhysics::ConstructProcess() 83 void IonUrQMDPhysics::ConstructProcess() 80 { 84 { 81 if (fWasActivated) { << 85 if(fWasActivated) { return; } 82 return; << 83 } << 84 fWasActivated = true; 86 fWasActivated = true; 85 87 86 G4double emin = 0. * MeV; << 88 G4double emin = 0.*MeV; 87 G4double emax = G4HadronicParameters::Instan << 89 G4double emax = 100.*TeV; 88 90 89 fModel = new G4UrQMD1_3Model(); 91 fModel = new G4UrQMD1_3Model(); 90 fModel->SetMinEnergy(emin); << 92 fModel->SetMinEnergy( emin ); 91 fModel->SetMaxEnergy(emax); << 93 fModel->SetMaxEnergy( emax ); 92 << 93 fIonXS = new G4CrossSectionInelastic(new G4C << 94 << 95 AddProcess("dInelastic", G4Deuteron::Deutero << 96 AddProcess("tInelastic", G4Triton::Triton()) << 97 AddProcess("He3Inelastic", G4He3::He3()); << 98 AddProcess("alphaInelastic", G4Alpha::Alpha( << 99 AddProcess("ionInelastic", G4GenericIon::Gen << 100 94 101 if (fVerbose > 1) { << 95 fShen = new G4IonsShenCrossSection(); 102 G4cout << "IonUrQMDPhysics::ConstructProce << 96 fTripathi = new G4TripathiCrossSection(); >> 97 fTripathiLight = new G4TripathiLightCrossSection(); >> 98 fIonH = new G4IonProtonCrossSection(); >> 99 fShen->SetMaxKinEnergy( emax ); >> 100 fTripathi->SetMaxKinEnergy( emax ); >> 101 fTripathiLight->SetMaxKinEnergy( emax ); >> 102 fIonH->SetMaxKinEnergy( emax ); >> 103 >> 104 >> 105 AddProcess("dInelastic", G4Deuteron::Deuteron(),false); >> 106 AddProcess("tInelastic",G4Triton::Triton(),false); >> 107 AddProcess("He3Inelastic",G4He3::He3(),true); >> 108 AddProcess("alphaInelastic", G4Alpha::Alpha(),true); >> 109 AddProcess("ionInelastic",G4GenericIon::GenericIon(),true); >> 110 >> 111 if(fVerbose > 1) { >> 112 G4cout << "IonUrQMDPhysics::ConstructProcess done! " >> 113 << G4endl; 103 } 114 } 104 } 115 } 105 116 106 //....oooOO0OOooo........oooOO0OOooo........oo 117 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 107 118 108 void IonUrQMDPhysics::AddProcess(const G4Strin << 119 void IonUrQMDPhysics::AddProcess(const G4String& name, >> 120 G4ParticleDefinition* part, >> 121 G4bool isIon) 109 { 122 { 110 G4HadronInelasticProcess* hadi = new G4Hadro 123 G4HadronInelasticProcess* hadi = new G4HadronInelasticProcess(name, part); 111 G4ProcessManager* pManager = part->GetProces 124 G4ProcessManager* pManager = part->GetProcessManager(); 112 pManager->AddDiscreteProcess(hadi); 125 pManager->AddDiscreteProcess(hadi); 113 hadi->AddDataSet(fIonXS); << 126 hadi->AddDataSet(fShen); 114 hadi->RegisterMe(fModel); << 127 // hadi->AddDataSet(fTripathi); >> 128 // hadi->AddDataSet(fTripathiLight); >> 129 if(isIon) { hadi->AddDataSet(fIonH); } >> 130 hadi->RegisterMe( fModel ); 115 } 131 } 116 132 117 //....oooOO0OOooo........oooOO0OOooo........oo 133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 118 #endif // URQMD << 134 #endif //URQMD 119 135