Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // 27 // 28 //--------------------------------------------------------------------------- 29 // 30 // ClassName: HadronPhysicsNuBeam 31 // 32 // Author: Julia Yarba, FNAL/CD (2013) 33 // created from (molded after) HadronPhysicsFTFP_BERT 34 // 35 // Modified: 36 // 37 //---------------------------------------------------------------------------- 38 // 39 #include <iomanip> 40 41 #include "G4HadronPhysicsNuBeam.hh" 42 #include "G4QGSPLundStrFragmProtonBuilder.hh" 43 #include "G4ProtonBuilder.hh" 44 #include "G4FTFPProtonBuilder.hh" 45 #include "G4BertiniProtonBuilder.hh" 46 #include "globals.hh" 47 #include "G4ios.hh" 48 #include "G4SystemOfUnits.hh" 49 50 #include "G4PhysListUtil.hh" 51 #include "G4HadronicParameters.hh" 52 53 // factory 54 #include "G4PhysicsConstructorFactory.hh" 55 // 56 G4_DECLARE_PHYSCONSTR_FACTORY(G4HadronPhysicsNuBeam); 57 58 G4HadronPhysicsNuBeam::G4HadronPhysicsNuBeam(G4int verb) : 59 G4HadronPhysicsNuBeam("hInelasticNuBeam",false) 60 { 61 G4HadronicParameters::Instance()->SetVerboseLevel(verb); 62 } 63 64 G4HadronPhysicsNuBeam::G4HadronPhysicsNuBeam(const G4String& name, G4bool quasiElastic) 65 : G4HadronPhysicsFTFP_BERT(name,quasiElastic) 66 { 67 // specific transition energies should be defined here 68 69 //minFTFP_neutron = 4.0*GeV; 70 //maxBERT_neutron = 5.0*GeV; 71 minFTFP_proton = 3.0*GeV; 72 maxFTFP_proton = 101.0*GeV; 73 //maxBERT_proton = 3.5*GeV; 74 //minFTFP_pion = minFTFP_kaon = 3.0*GeV; 75 //maxBERT_pion = maxBERT_kaon = 3.5*GeV; 76 } 77 78 void G4HadronPhysicsNuBeam::Proton() 79 { 80 G4HadronicParameters* param = G4HadronicParameters::Instance(); 81 G4bool useFactorXS = param->ApplyFactorXS(); 82 83 auto pro = new G4ProtonBuilder; 84 AddBuilder(pro); 85 // this is the new "custom" proton builder, tentatively for NuBeam 86 // 87 // no need to set the min energy because it's set in the ProBuilder 88 // ... and theMax will be set via Build() 89 // 90 // also explicitly set quasi-elastic key ON for QGS 91 // (it should be OFF for FTF, controlled by QuasiElastic) 92 // 93 auto qgsppro = new G4QGSPLundStrFragmProtonBuilder( true ); 94 AddBuilder(qgsppro); 95 pro->RegisterMe(qgsppro); 96 // 97 // standard FTFP builder, but energy range is adjusted 98 // 99 auto ftfppro = new G4FTFPProtonBuilder(QuasiElastic); 100 AddBuilder(ftfppro); 101 pro->RegisterMe(ftfppro); 102 ftfppro->SetMinEnergy(minFTFP_proton); 103 ftfppro->SetMaxEnergy(maxFTFP_proton); 104 // 105 // standard Bertini builder 106 // 107 auto bertpro = new G4BertiniProtonBuilder; 108 AddBuilder(bertpro); 109 pro->RegisterMe(bertpro); 110 bertpro->SetMaxEnergy(maxBERT_proton); 111 pro->Build(); 112 113 const G4ParticleDefinition* proton = G4Proton::Proton(); 114 G4HadronicProcess* inel = G4PhysListUtil::FindInelasticProcess(proton); 115 if(nullptr != inel) { 116 if( useFactorXS ) inel->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() ); 117 } 118 } 119 120 void G4HadronPhysicsNuBeam::ConstructProcess() 121 { 122 if(G4Threading::IsMasterThread() && 123 G4HadronicParameters::Instance()->GetVerboseLevel() > 0) { 124 DumpBanner(); 125 } 126 CreateModels(); 127 } 128