00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040 #include "G4QNeutronHPBuilder.hh"
00041 #include "G4SystemOfUnits.hh"
00042
00043 G4QNeutronHPBuilder::G4QNeutronHPBuilder():
00044 theNeutrons(0)
00045 , theNeutronFission(0)
00046 , theNeutronCapture(0)
00047 , theCHIPSNGamma(0)
00048 , theHPNeutron(0)
00049 , wasActivated(false)
00050 {
00051 theNeutronInelastic = new G4NeutronInelasticProcess;
00052 theCHIPSInelastic = new G4QInelastic;
00053 const G4ParticleDefinition* proj = G4Neutron::Neutron();
00054 const G4String& INprocessName = "MixedHPQNeutronInelasticProcess";
00055 theInProcessMixer= new G4QDiscProcessMixer(INprocessName, proj);
00056 const G4String& NGprocessName = "MixedHPQNeutronGammaProcess";
00057 theNgProcessMixer= new G4QDiscProcessMixer(NGprocessName, proj);
00058 const G4String& FIprocessName = "MixedHPQNeutronFissionProcess";
00059 theFiProcessMixer= new G4QDiscProcessMixer(FIprocessName, proj);
00060 }
00061
00062 G4QNeutronHPBuilder::~G4QNeutronHPBuilder()
00063 {
00064 delete theCHIPSInelastic;
00065 delete theCHIPSNGamma;
00066 delete theNeutronInelastic;
00067 delete theHPNeutron;
00068 delete theInProcessMixer;
00069 delete theNgProcessMixer;
00070 delete theFiProcessMixer;
00071 }
00072
00073 void G4QNeutronHPBuilder::Build()
00074 {
00075 wasActivated = true;
00076
00077 theNeutrons = new G4NeutronBuilder;
00078 theNeutrons->RegisterMe( theHPNeutron = new G4NeutronHPBuilder );
00079
00080 std::vector<G4VNeutronBuilder *>::iterator i;
00081 for(i = theModelCollections.begin(); i != theModelCollections.end(); i++)
00082 {
00083 (*i)->Build(theNeutronInelastic);
00084 (*i)->Build(theNeutronCapture);
00085 (*i)->Build(theNeutronFission);
00086 }
00087 G4ProcessManager * theProcMan = G4Neutron::Neutron()->GetProcessManager();
00088
00089 theInProcessMixer->AddDiscreteProcess(theCHIPSInelastic, 1.E8*megaelectronvolt);
00090 theInProcessMixer->AddDiscreteProcess(theNeutronInelastic, 19.9*megaelectronvolt);
00091
00092 theNgProcessMixer->AddDiscreteProcess(theCHIPSNGamma, 1.E8*megaelectronvolt);
00093 theNgProcessMixer->AddDiscreteProcess(theNeutronCapture, 19.9*megaelectronvolt);
00094
00095
00096
00097
00098 theProcMan->AddDiscreteProcess(theInProcessMixer);
00099 theProcMan->AddDiscreteProcess(theNgProcessMixer);
00100 theProcMan->AddDiscreteProcess(theNeutronFission);
00101
00102 }
00103