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/HadronPhysicsHIJ 26 /// \file hadronic/Hadr02/src/HadronPhysicsHIJING.cc 27 /// \brief Implementation of the HadronPhysics 27 /// \brief Implementation of the HadronPhysicsHIJING class 28 // 28 // 29 // 29 // 30 //-------------------------------------------- 30 //--------------------------------------------------------------------------- 31 // 31 // 32 // ClassName: << 32 // ClassName: 33 // 33 // 34 // Author: 2012 A. Dotti 34 // Author: 2012 A. Dotti 35 // created from HadronPhysicsHIJING 35 // created from HadronPhysicsHIJING 36 // 36 // 37 // Modified: 37 // Modified: 38 // 38 // 39 //-------------------------------------------- 39 //---------------------------------------------------------------------------- 40 // 40 // 41 #ifdef G4_USE_HIJING 41 #ifdef G4_USE_HIJING 42 # include "HadronPhysicsHIJING.hh" << 42 #include "HadronPhysicsHIJING.hh" 43 43 44 # include "G4AntiBarionBuilder.hh" << 44 #include "globals.hh" 45 # include "G4BaryonConstructor.hh" << 45 #include "G4ios.hh" 46 # include "G4BertiniNeutronBuilder.hh" << 46 #include <iomanip> 47 # include "G4BertiniPiKBuilder.hh" << 47 #include "G4ParticleDefinition.hh" 48 # include "G4BertiniProtonBuilder.hh" << 48 #include "G4ParticleTable.hh" 49 # include "G4ChipsKaonMinusInelasticXS.hh" << 49 50 # include "G4ChipsKaonPlusInelasticXS.hh" << 50 #include "G4MesonConstructor.hh" 51 # include "G4ChipsKaonZeroInelasticXS.hh" << 51 #include "G4BaryonConstructor.hh" 52 # include "G4CrossSectionDataSetRegistry.hh" << 52 #include "G4ShortLivedConstructor.hh" 53 # include "G4FTFPAntiBarionBuilder.hh" << 53 #include "G4ChipsKaonMinusInelasticXS.hh" 54 # include "G4FTFPNeutronBuilder.hh" << 54 #include "G4ChipsKaonPlusInelasticXS.hh" 55 # include "G4FTFPPiKBuilder.hh" << 55 #include "G4ChipsKaonZeroInelasticXS.hh" 56 # include "G4FTFPProtonBuilder.hh" << 56 #include "G4CrossSectionDataSetRegistry.hh" 57 # include "G4HyperonFTFPBuilder.hh" << 57 58 # include "G4MesonConstructor.hh" << 58 #include "G4ProcessManager.hh" 59 # include "G4NeutronBuilder.hh" << 59 #include "G4PhysListUtil.hh" 60 # include "G4ParticleDefinition.hh" << 60 61 # include "G4ParticleTable.hh" << 61 #include "G4PiKBuilder.hh" 62 # include "G4PhysListUtil.hh" << 62 #include "G4FTFPPiKBuilder.hh" 63 # include "G4PiKBuilder.hh" << 63 #include "G4QGSPPiKBuilder.hh" 64 # include "G4ProcessManager.hh" << 64 #include "G4BertiniPiKBuilder.hh" 65 # include "G4ProtonBuilder.hh" << 65 66 # include "G4QGSPNeutronBuilder.hh" << 66 #include "G4ProtonBuilder.hh" 67 # include "G4QGSPPiKBuilder.hh" << 67 #include "G4FTFPProtonBuilder.hh" 68 # include "G4QGSPProtonBuilder.hh" << 68 #include "G4QGSPProtonBuilder.hh" 69 # include "G4ShortLivedConstructor.hh" << 69 #include "G4BertiniProtonBuilder.hh" 70 # include "G4ios.hh" << 70 71 # include "globals.hh" << 71 #include "G4NeutronBuilder.hh" >> 72 #include "G4FTFPNeutronBuilder.hh" >> 73 #include "G4QGSPNeutronBuilder.hh" >> 74 #include "G4BertiniNeutronBuilder.hh" >> 75 >> 76 #include "G4HyperonFTFPBuilder.hh" >> 77 #include "G4AntiBarionBuilder.hh" >> 78 #include "G4FTFPAntiBarionBuilder.hh" 72 79 73 # include <iomanip> << 80 HadronPhysicsHIJING::HadronPhysicsHIJING(G4int) 74 << 81 : G4VPhysicsConstructor("hInelastic HIJING") 75 HadronPhysicsHIJING::HadronPhysicsHIJING(G4int << 76 { 82 { 77 fNeutrons = 0; 83 fNeutrons = 0; 78 fHIJINGNeutron = 0; 84 fHIJINGNeutron = 0; 79 fPiK = 0; 85 fPiK = 0; 80 fHIJINGPiK = 0; 86 fHIJINGPiK = 0; 81 fPro = 0; 87 fPro = 0; 82 fHIJINGPro = 0; << 88 fHIJINGPro = 0; 83 fHyperon = 0; 89 fHyperon = 0; 84 fFTFPHyperon = 0; 90 fFTFPHyperon = 0; 85 fAntiBaryon = 0; 91 fAntiBaryon = 0; 86 fHIJINGAntiBaryon = 0; 92 fHIJINGAntiBaryon = 0; 87 fCHIPSInelastic = 0; 93 fCHIPSInelastic = 0; 88 } 94 } 89 95 90 //....oooOO0OOooo........oooOO0OOooo........oo 96 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 91 97 92 void HadronPhysicsHIJING::CreateModels() 98 void HadronPhysicsHIJING::CreateModels() 93 { 99 { 94 G4bool quasiElasFTF = false; 100 G4bool quasiElasFTF = false; 95 G4bool quasiElasQGS = true; 101 G4bool quasiElasQGS = true; 96 G4double minFTFP = 6 * GeV; << 102 G4double minFTFP = 6*GeV; 97 G4double maxBERT = 8 * GeV; << 103 G4double maxBERT = 8*GeV; 98 G4double minHIJING = 12 * GeV; << 104 G4double minHIJING = 12*GeV; 99 G4double minQGSP = 12 * GeV; << 105 G4double minQGSP = 12*GeV; 100 G4double maxFTFP = 25 * GeV; << 106 G4double maxFTFP = 25*GeV; 101 << 107 102 // Proton << 108 //Proton 103 fPro = new G4ProtonBuilder; << 109 fPro=new G4ProtonBuilder; 104 fHIJINGPro = new HIJINGProtonBuilder(); << 110 fHIJINGPro=new HIJINGProtonBuilder(); 105 fHIJINGPro->SetMinEnergy(minHIJING); 111 fHIJINGPro->SetMinEnergy(minHIJING); 106 fPro->RegisterMe(fHIJINGPro); 112 fPro->RegisterMe(fHIJINGPro); 107 G4FTFPProtonBuilder* FTFPPro = new G4FTFPPro << 113 G4FTFPProtonBuilder* FTFPPro=new G4FTFPProtonBuilder(quasiElasFTF); 108 FTFPPro->SetMinEnergy(minFTFP); << 114 FTFPPro->SetMinEnergy( minFTFP ); 109 FTFPPro->SetMaxEnergy(maxFTFP); << 115 FTFPPro->SetMaxEnergy( maxFTFP ); 110 fPro->RegisterMe(FTFPPro); << 116 fPro->RegisterMe( FTFPPro ); 111 G4BertiniProtonBuilder* BertPro = new G4Bert 117 G4BertiniProtonBuilder* BertPro = new G4BertiniProtonBuilder(); 112 BertPro->SetMaxEnergy(maxBERT); << 118 BertPro->SetMaxEnergy( maxBERT ); 113 119 114 fNeutrons = new G4NeutronBuilder; << 120 fNeutrons=new G4NeutronBuilder; 115 fHIJINGNeutron = new HIJINGNeutronBuilder(); << 121 fHIJINGNeutron=new HIJINGNeutronBuilder(); 116 fHIJINGNeutron->SetMinEnergy(minHIJING); 122 fHIJINGNeutron->SetMinEnergy(minHIJING); 117 fNeutrons->RegisterMe(fHIJINGNeutron); 123 fNeutrons->RegisterMe(fHIJINGNeutron); 118 // G4QGSPNeutronBuilder* QGSPNeu = new G4QGS << 124 //G4QGSPNeutronBuilder* QGSPNeu = new G4QGSPNeutronBuilder(quasiElasQGS); 119 // QGSPNeu->SetMinEnergy(minQGSP); << 125 //QGSPNeu->SetMinEnergy(minQGSP); 120 // fNeutrons->RegisterMe(QGSPNeu); << 126 //fNeutrons->RegisterMe(QGSPNeu); 121 G4FTFPNeutronBuilder* FTFPNeu = new G4FTFPNe 127 G4FTFPNeutronBuilder* FTFPNeu = new G4FTFPNeutronBuilder(quasiElasFTF); 122 FTFPNeu->SetMinEnergy(minFTFP); << 128 FTFPNeu->SetMinEnergy( minFTFP ); 123 FTFPNeu->SetMaxEnergy(maxFTFP); << 129 FTFPNeu->SetMaxEnergy( maxFTFP ); 124 fNeutrons->RegisterMe(FTFPNeu); << 130 fNeutrons->RegisterMe( FTFPNeu ); 125 G4BertiniNeutronBuilder* BertNeu = new G4Ber 131 G4BertiniNeutronBuilder* BertNeu = new G4BertiniNeutronBuilder(); 126 BertNeu->SetMaxEnergy(maxBERT); << 132 BertNeu->SetMaxEnergy( maxBERT ); 127 133 128 fPiK = new G4PiKBuilder; << 134 fPiK=new G4PiKBuilder; 129 G4QGSPPiKBuilder* QGSPPiK = new G4QGSPPiKBui << 135 G4QGSPPiKBuilder* QGSPPiK=new G4QGSPPiKBuilder(quasiElasQGS); 130 fPiK->RegisterMe(QGSPPiK); 136 fPiK->RegisterMe(QGSPPiK); 131 QGSPPiK->SetMinEnergy(minQGSP); 137 QGSPPiK->SetMinEnergy(minQGSP); 132 G4FTFPPiKBuilder* FTFPPiK = new G4FTFPPiKBui << 138 G4FTFPPiKBuilder* FTFPPiK=new G4FTFPPiKBuilder(quasiElasFTF); 133 fPiK->RegisterMe(FTFPPiK); 139 fPiK->RegisterMe(FTFPPiK); 134 FTFPPiK->SetMaxEnergy(maxFTFP); 140 FTFPPiK->SetMaxEnergy(maxFTFP); 135 FTFPPiK->SetMinEnergy(minFTFP); << 141 FTFPPiK->SetMinEnergy(minFTFP); 136 G4BertiniPiKBuilder* BertiniPiK = new G4Bert << 142 G4BertiniPiKBuilder* BertiniPiK=new G4BertiniPiKBuilder(); 137 fPiK->RegisterMe(BertiniPiK); 143 fPiK->RegisterMe(BertiniPiK); 138 BertiniPiK->SetMaxEnergy(maxBERT); 144 BertiniPiK->SetMaxEnergy(maxBERT); 139 145 140 // For Hyperons use FTF model << 146 //For Hyperons use FTF model 141 fHyperon = new G4HyperonBuilder; << 147 fHyperon=new G4HyperonBuilder; 142 fFTFPHyperon = new G4HyperonFTFPBuilder; << 148 fFTFPHyperon=new G4HyperonFTFPBuilder; 143 fHyperon->RegisterMe(fFTFPHyperon); << 149 fHyperon->RegisterMe( fFTFPHyperon ); 144 << 150 145 fAntiBaryon = new G4AntiBarionBuilder; << 151 fAntiBaryon=new G4AntiBarionBuilder; 146 // FTFPAntiBaryon=new G4FTFPAntiBarionBuilde << 152 //FTFPAntiBaryon=new G4FTFPAntiBarionBuilder(quasiElasFTF); 147 G4FTFPAntiBarionBuilder* FTFPAB = new G4FTFP 153 G4FTFPAntiBarionBuilder* FTFPAB = new G4FTFPAntiBarionBuilder(quasiElasFTF); 148 fAntiBaryon->RegisterMe(FTFPAB); << 154 fAntiBaryon->RegisterMe( FTFPAB ); 149 } 155 } 150 156 151 //....oooOO0OOooo........oooOO0OOooo........oo 157 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 152 158 153 HadronPhysicsHIJING::~HadronPhysicsHIJING() 159 HadronPhysicsHIJING::~HadronPhysicsHIJING() 154 { 160 { 155 delete fNeutrons; 161 delete fNeutrons; 156 delete fHIJINGNeutron; 162 delete fHIJINGNeutron; 157 163 158 delete fPiK; 164 delete fPiK; 159 delete fHIJINGPiK; 165 delete fHIJINGPiK; 160 << 166 161 delete fPro; 167 delete fPro; 162 delete fHIJINGPro; << 168 delete fHIJINGPro; 163 << 169 164 delete fHyperon; 170 delete fHyperon; 165 delete fFTFPHyperon; 171 delete fFTFPHyperon; 166 delete fAntiBaryon; 172 delete fAntiBaryon; 167 delete fHIJINGAntiBaryon; 173 delete fHIJINGAntiBaryon; 168 << 174 169 delete fCHIPSInelastic; 175 delete fCHIPSInelastic; 170 } 176 } 171 177 172 //....oooOO0OOooo........oooOO0OOooo........oo 178 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 173 179 174 void HadronPhysicsHIJING::ConstructParticle() 180 void HadronPhysicsHIJING::ConstructParticle() 175 { 181 { 176 G4MesonConstructor pMesonConstructor; 182 G4MesonConstructor pMesonConstructor; 177 pMesonConstructor.ConstructParticle(); 183 pMesonConstructor.ConstructParticle(); 178 184 179 G4BaryonConstructor pBaryonConstructor; 185 G4BaryonConstructor pBaryonConstructor; 180 pBaryonConstructor.ConstructParticle(); 186 pBaryonConstructor.ConstructParticle(); 181 187 182 G4ShortLivedConstructor pShortLivedConstruct 188 G4ShortLivedConstructor pShortLivedConstructor; 183 pShortLivedConstructor.ConstructParticle(); << 189 pShortLivedConstructor.ConstructParticle(); 184 } 190 } 185 191 186 //....oooOO0OOooo........oooOO0OOooo........oo 192 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 187 193 188 void HadronPhysicsHIJING::ConstructProcess() 194 void HadronPhysicsHIJING::ConstructProcess() 189 { 195 { 190 CreateModels(); 196 CreateModels(); 191 fNeutrons->Build(); 197 fNeutrons->Build(); 192 fPro->Build(); 198 fPro->Build(); 193 fPiK->Build(); 199 fPiK->Build(); 194 200 195 // use CHIPS cross sections also for Kaons << 201 // use CHIPS cross sections also for Kaons 196 G4VCrossSectionDataSet* ChipsKaonMinus = << 202 G4VCrossSectionDataSet* ChipsKaonMinus = 197 G4CrossSectionDataSetRegistry::Instance()- << 203 G4CrossSectionDataSetRegistry::Instance()-> 198 G4ChipsKaonMinusInelasticXS::Default_Nam << 204 GetCrossSectionDataSet(G4ChipsKaonMinusInelasticXS::Default_Name()); 199 G4VCrossSectionDataSet* ChipsKaonPlus = << 205 G4VCrossSectionDataSet* ChipsKaonPlus = 200 G4CrossSectionDataSetRegistry::Instance()- << 206 G4CrossSectionDataSetRegistry::Instance()-> 201 G4ChipsKaonPlusInelasticXS::Default_Name << 207 GetCrossSectionDataSet(G4ChipsKaonPlusInelasticXS::Default_Name()); 202 G4VCrossSectionDataSet* ChipsKaonZero = << 208 G4VCrossSectionDataSet* ChipsKaonZero = 203 G4CrossSectionDataSetRegistry::Instance()- << 209 G4CrossSectionDataSetRegistry::Instance()-> 204 G4ChipsKaonZeroInelasticXS::Default_Name << 210 GetCrossSectionDataSet(G4ChipsKaonZeroInelasticXS::Default_Name()); 205 // << 211 // 206 << 212 207 G4PhysListUtil::FindInelasticProcess(G4KaonM << 213 G4PhysListUtil::FindInelasticProcess(G4KaonMinus::KaonMinus())-> 208 G4PhysListUtil::FindInelasticProcess(G4KaonP << 214 AddDataSet(ChipsKaonMinus); 209 G4PhysListUtil::FindInelasticProcess(G4KaonZ << 215 G4PhysListUtil::FindInelasticProcess(G4KaonPlus::KaonPlus())-> 210 G4PhysListUtil::FindInelasticProcess(G4KaonZ << 216 AddDataSet(ChipsKaonPlus); 211 << 217 G4PhysListUtil::FindInelasticProcess(G4KaonZeroShort::KaonZeroShort())-> >> 218 AddDataSet(ChipsKaonZero ); >> 219 G4PhysListUtil::FindInelasticProcess(G4KaonZeroLong::KaonZeroLong())-> >> 220 AddDataSet(ChipsKaonZero ); >> 221 212 fHyperon->Build(); 222 fHyperon->Build(); 213 fAntiBaryon->Build(); 223 fAntiBaryon->Build(); 214 } 224 } 215 225 216 //....oooOO0OOooo........oooOO0OOooo........oo 226 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 217 227 218 G4HadronicProcess* HadronPhysicsHIJING::FindIn << 228 G4HadronicProcess* >> 229 HadronPhysicsHIJING::FindInelasticProcess(const G4ParticleDefinition* p) 219 { 230 { 220 G4HadronicProcess* had = 0; 231 G4HadronicProcess* had = 0; 221 if (p) { << 232 if(p) { 222 G4ProcessVector* pvec = p->GetProcessManag << 233 G4ProcessVector* pvec = p->GetProcessManager()->GetProcessList(); 223 size_t n = pvec->size(); << 234 size_t n = pvec->size(); 224 if (0 < n) { << 235 if(0 < n) { 225 for (size_t i = 0; i < n; ++i) { << 236 for(size_t i=0; i<n; ++i) { 226 if (fHadronInelastic == ((*pvec)[i])-> << 237 if(fHadronInelastic == ((*pvec)[i])->GetProcessSubType()) { 227 had = static_cast<G4HadronicProcess* << 238 had = static_cast<G4HadronicProcess*>((*pvec)[i]); 228 << 239 229 break; << 240 break; 230 } << 241 } 231 } << 242 } 232 } << 243 } 233 } 244 } 234 return had; 245 return had; 235 } 246 } 236 247 237 #endif // G4_USE_HIJING << 248 #endif //G4_USE_HIJING 238 249