Message: Re: <14MeV neutron setup Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: <14MeV neutron setup 

Keywords: neutron interaction hadronic process
Forum: Particles
Re: Question <14MeV neutron setup (Andrew)
Date: 25 Apr, 2005
From: Dennis H. Wright <Dennis H. Wright>

Hi Andrew,

   The following code should do the trick:

  ************************************************************

  G4ProcessManager* pManager = G4Neutron::Neutron()->GetProcessManager();

  // Neutron elastic process, models and cross sections

  G4HadronElasticProcess* elasticProcess = new G4HadronElasticProcess();
  G4NeutronHPElastic* theHPElasticModel = new G4NeutronHPElastic();
  G4NeutronHPElasticData* theHPElasticData = new G4NeutronHPElasticData();
  elasticProcess->RegisterMe(theHPElasticModel);
  elasticProcess->AddDataSet(theHPElasticData);
  pManager->AddDiscreteProcess(elasticProcess);

  // Neutron inelastic process, models and cross sections

  G4NeutronInelasticProcess* ninelProc = new G4NeutronInelasticProcess();
  G4NeutronHPInelastic* theHPInelasticModel = new G4NeutronHPInelastic();
  G4NeutronHPInelasticData* theHPInelasticData = new G4NeutronHPInelasticData();
  ninelProc->RegisterMe(theHPInelasticModel);
  ninelProc->AddDataSet(theHPInelasticData);
  pManager->AddDiscreteProcess(ninelProc);

  // Neutron-induced fission process, models and cross sections

  G4HadronFissionProcess* neutronFission = new G4HadronFissionProcess();
  G4NeutronHPFission* theHPFissionModel = new G4NeutronHPFission();
  G4NeutronHPFissionData* theHPFissionData=new G4NeutronHPFissionData();
  neutronFission->RegisterMe(theHPFissionModel);
  neutronFission->AddDataSet(theHPFissionData);
  pManager->AddDiscreteProcess(neutronFission);

  // Neutron capture process, models and cross sections

  G4HadronCaptureProcess* neutronCapture = new G4HadronCaptureProcess();
  G4NeutronHPCapture* theHPCaptureModel = new G4NeutronHPCapture();
  G4NeutronHPCaptureData* theHPCaptureData = new G4NeutronHPCaptureData();
  neutronCapture->RegisterMe(theHPCaptureModel);
  neutronCapture->AddDataSet(theHPCaptureData);
  pManager->AddDiscreteProcess(neutronCapture);

  ********************************************************************

You will also need to set an environment variable:

setenv NeutronHPCrossSections somedir/G4NDL3.7

where somedir is where you have put your copy of the G4 Neutron
Data Library.


Dennis


Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

 Add Message Add Message
to: "Re: <14MeV neutron setup"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]