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 /// \file hadronic/Hadr02/src/IonCRMCPhysics.c 27 /// \brief Implementation of the IonCRMCPhysic 28 // 29 // 30 //-------------------------------------------- 31 // 32 // Class: IonCRMCPhysics 33 // 34 // Author: 2018 Alberto Ribon 35 // 36 // Modified: 37 // - 18-May-2021 Alberto Ribon : Used the lat 38 // 39 //-------------------------------------------- 40 // 41 #ifdef G4_USE_CRMC 42 43 # include "IonCRMCPhysics.hh" 44 45 # include "HadronicInelasticModelCRMC.hh" 46 47 # include "G4Alpha.hh" 48 # include "G4BinaryLightIonReaction.hh" 49 # include "G4BuilderType.hh" 50 # include "G4ComponentGGNuclNuclXsc.hh" 51 # include "G4CrossSectionInelastic.hh" 52 # include "G4Deuteron.hh" 53 # include "G4ExcitationHandler.hh" 54 # include "G4FTFBuilder.hh" 55 # include "G4GenericIon.hh" 56 # include "G4HadronInelasticProcess.hh" 57 # include "G4HadronicInteraction.hh" 58 # include "G4HadronicInteractionRegistry.hh" 59 # include "G4HadronicParameters.hh" 60 # include "G4He3.hh" 61 # include "G4IonConstructor.hh" 62 # include "G4IonPhysics.hh" 63 # include "G4ParticleDefinition.hh" 64 # include "G4PreCompoundModel.hh" 65 # include "G4ProcessManager.hh" 66 # include "G4SystemOfUnits.hh" 67 # include "G4Triton.hh" 68 69 using namespace std; 70 71 // factory 72 # include "G4PhysicsConstructorFactory.hh" 73 // 74 G4_DECLARE_PHYSCONSTR_FACTORY(IonCRMCPhysics); 75 76 //....oooOO0OOooo........oooOO0OOooo........oo 77 78 const std::array<std::string, 13> IonCRMCPhysi 79 "EPOS-LHC", "EPOS-1.99", "QGSJET-01", "", 80 "", "", "QGSJETII-03", "DPM 81 82 //....oooOO0OOooo........oooOO0OOooo........oo 83 84 IonCRMCPhysics::IonCRMCPhysics(G4int ver) : G4 85 { 86 fModel = 0; //***LOOKHERE*** CRMC model: 0: 87 // 7: 88 fVerbose = ver; 89 if (fVerbose > 1) G4cout << "### IonCRMCPhys 90 SetPhysicsType(bIons); 91 } 92 93 //....oooOO0OOooo........oooOO0OOooo........oo 94 95 IonCRMCPhysics::~IonCRMCPhysics() {} 96 97 //....oooOO0OOooo........oooOO0OOooo........oo 98 99 void IonCRMCPhysics::ConstructParticle() 100 { 101 // Construct ions 102 G4IonConstructor pConstructor; 103 pConstructor.ConstructParticle(); 104 } 105 106 //....oooOO0OOooo........oooOO0OOooo........oo 107 108 void IonCRMCPhysics::ConstructProcess() 109 { 110 fModel = 0; //***LOOKHERE*** 0:EPOS-LHC, 1: 111 // 7:QGSJETII-04, 112 const G4double minCRMC = 113 100.0 114 * GeV; //***LOOKHERE*** CRMC model is app 115 const G4double maxFTFP = 116 110.0 117 * GeV; //***LOOKHERE*** FTFP model is app 118 G4HadronicInteraction* p = G4HadronicInterac 119 G4PreCompoundModel* thePreCompound = static_ 120 if (!thePreCompound) thePreCompound = new G4 121 // Binary Cascade 122 G4HadronicInteraction* theIonBC = new G4Bina 123 theIonBC->SetMinEnergy(0.0); 124 theIonBC->SetMaxEnergy(G4HadronicParameters: 125 // FTFP 126 G4FTFBuilder theBuilder("FTFP", thePreCompou 127 G4HadronicInteraction* theFTFP = theBuilder. 128 theFTFP->SetMinEnergy(G4HadronicParameters:: 129 theFTFP->SetMaxEnergy(maxFTFP); 130 // CRMC 131 G4HadronicInteraction* theCRMC = new Hadroni 132 theCRMC->SetMinEnergy(minCRMC); 133 theCRMC->SetMaxEnergy(G4HadronicParameters:: 134 // Cross section 135 G4CrossSectionInelastic* theXS = new G4Cross 136 // Processes 137 AddProcess("dInelastic", G4Deuteron::Deutero 138 AddProcess("tInelastic", G4Triton::Triton(), 139 AddProcess("He3Inelastic", G4He3::He3(), the 140 AddProcess("alphaInelastic", G4Alpha::Alpha( 141 AddProcess("ionInelastic", G4GenericIon::Gen 142 if (fVerbose > 1) G4cout << "IonCRMCPhysics: 143 } 144 145 //....oooOO0OOooo........oooOO0OOooo........oo 146 147 void IonCRMCPhysics::AddProcess(const G4String 148 G4HadronicInte 149 G4HadronicInte 150 { 151 G4HadronInelasticProcess* hadi = new G4Hadro 152 G4ProcessManager* pManager = part->GetProces 153 pManager->AddDiscreteProcess(hadi); 154 if (xs) hadi->AddDataSet(xs); 155 if (theIonBC) hadi->RegisterMe(theIonBC); 156 if (theFTFP) hadi->RegisterMe(theFTFP); 157 if (theCRMC) hadi->RegisterMe(theCRMC); 158 } 159 160 #endif // G4_USE_CRMC 161