confused on the results about neutron yield and spectra produced from reaction (d,n)

Help: Forums, postings and threads

Forums, postings and threads

Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.

Layout of the posting thread(s)

The posting threads can be shown with/without the text of the postings ("Inline depth"). In addition the depth of the thread can be controlled ("Outline depth"). For both cases several values are provided, which can be toggled. Note that the displayed values are those, which can be chosen, they do not correspond to the actual display status of the thread(s).

Add message

Add a new posting thread to this Forum

Show subscribers

Click "Show subscribers" to get the List of all subscribers, who receive e-mail from postings to this Forum, if they have not explictly unsubscribed from a posting thread. This includes:

Subscribers to the forum
Explicit subscribers to the posting or any of its ancestors
The author of the posting (if he has not unsubscribed)
The authors of any ancestor posting (if they have not unsubscribed)

Some browsers display old pages from the cache

Some browsers (e.g. Firefox) display a page from the cache even though the page has been updated. As a result, one may not get the latest messages in the forum's outline or the logon ID of the previous session appears on the upper right corner. If one suspects that he or she does not get the latest content of a page, try clearing up the personal data from the browser's cache as the following:

Tools -> Clear Private Data -> Clear Private Data (from browser Firefox)

Tools -> Delete Private Data (from browser Opera)

then reflesh the current page to get the latest content.

Details on "Inline depth" and "Outline depth"

Threads and the follow-up messages are numbered. There is one line per posting, normally starting with the posting number. The URLs of a posting and its responses reflect the depth of the follow-up postings. Watch the URLs to see how this is done. Depending on the chosen display options, the thread-display may be incomplete; the number of messages at the next level is, however, displayed, e.g. "…3 message(s)". You may not see all contributors unless you choose "All".

Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.

Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.

Argument Default value Alternative Comment

1 0 1 no text / full text of posting

2 1 All text for level 1 only / text for All postings

Argument	Default value	Alternative		Comment
1	0	1		no text / full text of posting
2	1	All		text for level 1 only / text for All postings

Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.

Argument Default value Alternative Comment

1 2 1 2 levels / 1 level (original posting)

2 3 2 3 levels / 2 levels

3 3 All 3 levels / all levels (all postings)

Argument	Default value	Alternative	Comment
1	2	1	2 levels / 1 level (original posting)
2	3	2	3 levels / 2 levels
3	3	All	3 levels / all levels (all postings)

confused on the results about neutron yield and spectra produced from reaction (d,n)

Hi all, I am trying to simulate the neutron production yield and neutron spectra from the reactions: 9Be(d,n)10B and 12C(d,n)13N. Here is physics processes that I applied for deuteron:

#include "IonPhysics.hh"

#include "globals.hh"
#include "G4ios.hh"
#include <iomanip>

// Elastic scattering, multiple scattering and ionisation

IonPhysics::IonPhysics(const G4String& name)
                 :  G4VPhysicsConstructor(name)
{
}

IonPhysics::~IonPhysics() { }

#include "G4ParticleDefinition.hh"
#include "G4ParticleTable.hh"

#include "G4ProcessManager.hh"

void IonPhysics::ConstructProcess()
{
   G4ProcessManager * pManager = 0;

      // Deuteron
   pManager = G4Deuteron::Deuteron()->GetProcessManager();

   // add process
   G4HadronElasticProcess* thedueElasticProcess
                         = new G4HadronElasticProcess();
   G4LElastic* thedueElasticModel = new G4LElastic();
   thedueElasticProcess->RegisterMe(thedueElasticModel);
   pManager->AddDiscreteProcess(thedueElasticProcess);

   G4VProcess* thedueMultipleScattering = new G4MultipleScattering();
   G4VProcess* thedueIonisation        = new G4hIonisation();
   //
   pManager->AddProcess(thedueIonisation);
   pManager->AddProcess(thedueMultipleScattering);
   //
   // set ordering for AlongStepDoIt
   pManager->SetProcessOrdering(thedueMultipleScattering, idxAlongStep,1);
   pManager->SetProcessOrdering(thedueIonisation,        idxAlongStep,2);
   //
   // set ordering for PostStepDoIt
   pManager->SetProcessOrdering(thedueMultipleScattering, idxPostStep,1);
   pManager->SetProcessOrdering(thedueIonisation,        idxPostStep,2);

} // inelastic processes

#include "PhysListIonBinaryCascade.hh"
#include "G4ParticleDefinition.hh"
#include "G4ProcessManager.hh"

#include "G4Deuteron.hh"
#include "G4Triton.hh"
#include "G4Alpha.hh"
#include "G4GenericIon.hh"

#include "G4HadronInelasticProcess.hh"
#include "G4HadronElasticProcess.hh"
#include "G4BinaryLightIonReaction.hh"
#include "G4TripathiCrossSection.hh"
#include "G4TripathiLightCrossSection.hh"
#include "G4IonsShenCrossSection.hh"
#include "G4LEDeuteronInelastic.hh"
#include "G4DeuteronInelasticProcess.hh"
#include "G4LETritonInelastic.hh"
#include "G4LEAlphaInelastic.hh"
#include "G4LElastic.hh"

//EM process

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

PhysListIonBinaryCascade::PhysListIonBinaryCascade(const G4String& name)
   :  G4VPhysicsConstructor(name)
{}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

PhysListIonBinaryCascade::~PhysListIonBinaryCascade() {}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void PhysListIonBinaryCascade::ConstructProcess()
{ 
  // Binary Cascade
  G4ParticleDefinition* particle = 0;
  G4ProcessManager* pmanager = 0;

 // Cross section data sets

  // TRIPATHI CROSS SECTION
  // Implementation of formulas taken from NASA technical paper 3621 by 
  // Tripathi, et al. Cross-sections for ion ion scattering

  G4TripathiCrossSection * TripathiCrossSection= new G4TripathiCrossSection;

  // IONS SHEN CROSS SECTION
  // Implementation of formulas 
  // Shen et al. Nuc. Phys. A 491 130 (1989) 
  // Total Reaction Cross Section for Heavy-Ion Collisions

  G4IonsShenCrossSection * aShen = new G4IonsShenCrossSection;

  // Intra-nuclear transport: Binary Cascade Model
  // Binary Cascade for deuteron, triton, alpha particle

  G4BinaryLightIonReaction* theIonBC = new G4BinaryLightIonReaction();

  theIonBC->SetMinEnergy(80*MeV);
  theIonBC->SetMaxEnergy(10*GeV);

  // deuteron
  particle = G4Deuteron::Deuteron();
  pmanager = particle->GetProcessManager();

 // Final state production model for deuteron inelastic scattering below 100 MeV: Low Energy Parameterised model

  G4LEDeuteronInelastic* theDIModel = new G4LEDeuteronInelastic;

  // Set the maximum energy for LEP model

  theDIModel->SetMaxEnergy(100*MeV);

 // Activate the Tripathi and Shen Cross Section

  theIPdeuteron.AddDataSet(TripathiCrossSection);
  theIPdeuteron.AddDataSet(aShen);

 // Register the Parameterised Deuteron Inelastic Model and the Ion Binary Cascade Model

  theIPdeuteron.RegisterMe(theDIModel);
  //theIPdeuteron.RegisterMe(theTheoModel);
  theIPdeuteron.RegisterMe(theIonBC);

 // Activate inelastic scattering

  pmanager->AddDiscreteProcess(&theIPdeuteron);

Here is physics processes for neutron:

#include "HadronPhysics.hh"

#include "globals.hh"
#include "G4ios.hh"
#include <iomanip>

HadronPhysics::HadronPhysics(const G4String& name)
                   :  G4VPhysicsConstructor(name)
{
}

HadronPhysics::~HadronPhysics() { }

#include "G4ParticleDefinition.hh"
#include "G4ParticleTable.hh"

#include "G4ProcessManager.hh"
#include "G4NeutronInelasticCrossSection.hh"

void HadronPhysics::ConstructProcess()
{
   G4ProcessManager * pManager = 0;

   // neutron Physics
   pManager = G4Neutron::Neutron()->GetProcessManager();

   // add process
   G4HadronElasticProcess* thenElasticProcess 
                         = new G4HadronElasticProcess();
   G4LElastic* thenElasticModel = new G4LElastic();
   thenElasticProcess->RegisterMe(thenElasticModel);
   pManager->AddDiscreteProcess(thenElasticProcess);

   G4NeutronInelasticCrossSection* theNeutronCrossSection
						 = new G4NeutronInelasticCrossSection() ;

   G4NeutronInelasticProcess* theNeutronInelasticProcess 
                         = new G4NeutronInelasticProcess();

   G4LENeutronInelastic* theNeutronLEPModel = new G4LENeutronInelastic();
   G4HENeutronInelastic* theNeutronHEPModel = new G4HENeutronInelastic();
   //theNeutronHEPModel->SetMaxEnergy(20.*MeV);
   theNeutronInelasticProcess->AddDataSet(theNeutronCrossSection);

   theNeutronInelasticProcess->RegisterMe(theNeutronLEPModel);
   theNeutronInelasticProcess->RegisterMe(theNeutronHEPModel);
   pManager->AddDiscreteProcess(theNeutronInelasticProcess);

   //neutron + fission 
   G4HadronFissionProcess* thenFission
                         = new G4HadronFissionProcess();
   G4LFission* thenFissionModel = new G4LFission();
   thenFission->RegisterMe(thenFissionModel);
   pManager->AddDiscreteProcess(thenFission);

   //neutron + capture

   G4HadronCaptureProcess* thenCapture
                         = new G4HadronCaptureProcess();
   G4LCapture* thenCaptureModel = new G4LCapture();
   thenCapture->RegisterMe(thenCaptureModel);
   pManager->AddDiscreteProcess(thenCapture);

   // anti-neutron Physics
   pManager = G4AntiNeutron::AntiNeutron()->GetProcessManager();

   // add process
   G4HadronElasticProcess* theanElasticProcess 
                         = new G4HadronElasticProcess();
   G4LElastic* theanElasticModel = new G4LElastic();
   theanElasticProcess->RegisterMe(theanElasticModel);
   pManager->AddDiscreteProcess(theanElasticProcess);

   G4AntiNeutronInelasticProcess* theAntiNeutronInelasticProcess 
                         = new G4AntiNeutronInelasticProcess();

   G4LEAntiNeutronInelastic* theAntiNeutronLEPModel = new G4LEAntiNeutronInelastic();
   G4HEAntiNeutronInelastic* theAntiNeutronHEPModel = new G4HEAntiNeutronInelastic();
   theAntiNeutronInelasticProcess->RegisterMe(theAntiNeutronLEPModel);
   theAntiNeutronInelasticProcess->RegisterMe(theAntiNeutronHEPModel);
   pManager->AddDiscreteProcess(theAntiNeutronInelasticProcess);

   G4AntiNeutronAnnihilationAtRest* theAntiNeutronAnnihilation
                            =  new G4AntiNeutronAnnihilationAtRest();
   pManager->AddRestProcess(theAntiNeutronAnnihilation);

However, at low energy ( deuteron which has energy lower than 100 MeV) my program can not reproduce the experimental data about neutron yields and spectra. Speccially, with 160 MeV deuteron i can not see any neutron.

It would be very grateful if you can let me know what is error in my program and suggest me some ideas to solve this problem.

Thank you very much in advance!

Phan Viet Cuong

Re: confused on the results about neutron yield and spectra produced from reaction (d,n) by Phan Viet Cuong <Phan Viet Cuong>, 12 Aug, 2008 Re: confused on the results about neutron yield and spectra produced from reaction (d,n) (Phan Viet Cuong) Dear All, I do not know why my question was removed from the forum physicslist so I want to post it again. Cuong ---------------------------------------------------------------------- Hi all, I am trying to simulate the neutron production yield and neutron spectra from the reactions: 9Be(d,n)10B and 12C(d,n)13N. Here is physics processes that I applied for deuteron: #include "IonPhysics.hh" #include "globals.hh" #include "G4ios.hh" #include <iomanip> // Elastic scattering, multiple scattering and ionisation IonPhysics::IonPhysics(const G4String& name) : G4VPhysicsConstructor(name) { } IonPhysics::~IonPhysics() { } #include "G4ParticleDefinition.hh" #include "G4ParticleTable.hh" #include "G4ProcessManager.hh" void IonPhysics::ConstructProcess() { G4ProcessManager * pManager = 0; // Deuteron pManager = G4Deuteron::Deuteron()->GetProcessManager(); // add process G4HadronElasticProcess* thedueElasticProcess = new G4HadronElasticProcess(); G4LElastic* thedueElasticModel = new G4LElastic(); thedueElasticProcess->RegisterMe(thedueElasticModel); pManager->AddDiscreteProcess(thedueElasticProcess); G4VProcess* thedueMultipleScattering = new G4MultipleScattering(); G4VProcess* thedueIonisation = new G4hIonisation(); // pManager->AddProcess(thedueIonisation); pManager->AddProcess(thedueMultipleScattering); // // set ordering for AlongStepDoIt pManager->SetProcessOrdering(thedueMultipleScattering, idxAlongStep,1); pManager->SetProcessOrdering(thedueIonisation, idxAlongStep,2); // // set ordering for PostStepDoIt pManager->SetProcessOrdering(thedueMultipleScattering, idxPostStep,1); pManager->SetProcessOrdering(thedueIonisation, idxPostStep,2); } // inelastic processes #include "PhysListIonBinaryCascade.hh" #include "G4ParticleDefinition.hh" #include "G4ProcessManager.hh" #include "G4Deuteron.hh" #include "G4Triton.hh" #include "G4Alpha.hh" #include "G4GenericIon.hh" #include "G4HadronInelasticProcess.hh" #include "G4HadronElasticProcess.hh" #include "G4BinaryLightIonReaction.hh" #include "G4TripathiCrossSection.hh" #include "G4TripathiLightCrossSection.hh" #include "G4IonsShenCrossSection.hh" #include "G4LEDeuteronInelastic.hh" #include "G4DeuteronInelasticProcess.hh" #include "G4LETritonInelastic.hh" #include "G4LEAlphaInelastic.hh" #include "G4LElastic.hh" //EM process //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... PhysListIonBinaryCascade::PhysListIonBinaryCascade(const G4String& name) : G4VPhysicsConstructor(name) {} //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... PhysListIonBinaryCascade::~PhysListIonBinaryCascade() {} //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... void PhysListIonBinaryCascade::ConstructProcess() { // Binary Cascade G4ParticleDefinition* particle = 0; G4ProcessManager* pmanager = 0; // Cross section data sets // TRIPATHI CROSS SECTION // Implementation of formulas taken from NASA technical paper 3621 by // Tripathi, et al. Cross-sections for ion ion scattering G4TripathiCrossSection * TripathiCrossSection= new G4TripathiCrossSection; // IONS SHEN CROSS SECTION // Implementation of formulas // Shen et al. Nuc. Phys. A 491 130 (1989) // Total Reaction Cross Section for Heavy-Ion Collisions G4IonsShenCrossSection * aShen = new G4IonsShenCrossSection; // Intra-nuclear transport: Binary Cascade Model // Binary Cascade for deuteron, triton, alpha particle G4BinaryLightIonReaction* theIonBC = new G4BinaryLightIonReaction(); theIonBC->SetMinEnergy(80*MeV); theIonBC->SetMaxEnergy(10*GeV); // deuteron particle = G4Deuteron::Deuteron(); pmanager = particle->GetProcessManager(); // Final state production model for deuteron inelastic scattering below 100 MeV: Low Energy Parameterised model G4LEDeuteronInelastic* theDIModel = new G4LEDeuteronInelastic; // Set the maximum energy for LEP model theDIModel->SetMaxEnergy(100*MeV); // Activate the Tripathi and Shen Cross Section theIPdeuteron.AddDataSet(TripathiCrossSection); theIPdeuteron.AddDataSet(aShen); // Register the Parameterised Deuteron Inelastic Model and the Ion Binary Cascade Model theIPdeuteron.RegisterMe(theDIModel); //theIPdeuteron.RegisterMe(theTheoModel); theIPdeuteron.RegisterMe(theIonBC); // Activate inelastic scattering pmanager->AddDiscreteProcess(&theIPdeuteron); Here is physics processes for neutron: #include "HadronPhysics.hh" #include "globals.hh" #include "G4ios.hh" #include <iomanip> HadronPhysics::HadronPhysics(const G4String& name) : G4VPhysicsConstructor(name) { } HadronPhysics::~HadronPhysics() { } #include "G4ParticleDefinition.hh" #include "G4ParticleTable.hh" #include "G4ProcessManager.hh" #include "G4NeutronInelasticCrossSection.hh" void HadronPhysics::ConstructProcess() { G4ProcessManager * pManager = 0; // neutron Physics pManager = G4Neutron::Neutron()->GetProcessManager(); // add process G4HadronElasticProcess* thenElasticProcess = new G4HadronElasticProcess(); G4LElastic* thenElasticModel = new G4LElastic(); thenElasticProcess->RegisterMe(thenElasticModel); pManager->AddDiscreteProcess(thenElasticProcess); G4NeutronInelasticCrossSection* theNeutronCrossSection = new G4NeutronInelasticCrossSection() ; G4NeutronInelasticProcess* theNeutronInelasticProcess = new G4NeutronInelasticProcess(); G4LENeutronInelastic* theNeutronLEPModel = new G4LENeutronInelastic(); G4HENeutronInelastic* theNeutronHEPModel = new G4HENeutronInelastic(); //theNeutronHEPModel->SetMaxEnergy(20.*MeV); theNeutronInelasticProcess->AddDataSet(theNeutronCrossSection); theNeutronInelasticProcess->RegisterMe(theNeutronLEPModel); theNeutronInelasticProcess->RegisterMe(theNeutronHEPModel); pManager->AddDiscreteProcess(theNeutronInelasticProcess); //neutron + fission G4HadronFissionProcess* thenFission = new G4HadronFissionProcess(); G4LFission* thenFissionModel = new G4LFission(); thenFission->RegisterMe(thenFissionModel); pManager->AddDiscreteProcess(thenFission); //neutron + capture G4HadronCaptureProcess* thenCapture = new G4HadronCaptureProcess(); G4LCapture* thenCaptureModel = new G4LCapture(); thenCapture->RegisterMe(thenCaptureModel); pManager->AddDiscreteProcess(thenCapture); // anti-neutron Physics pManager = G4AntiNeutron::AntiNeutron()->GetProcessManager(); // add process G4HadronElasticProcess* theanElasticProcess = new G4HadronElasticProcess(); G4LElastic* theanElasticModel = new G4LElastic(); theanElasticProcess->RegisterMe(theanElasticModel); pManager->AddDiscreteProcess(theanElasticProcess); G4AntiNeutronInelasticProcess* theAntiNeutronInelasticProcess = new G4AntiNeutronInelasticProcess(); G4LEAntiNeutronInelastic* theAntiNeutronLEPModel = new G4LEAntiNeutronInelastic(); G4HEAntiNeutronInelastic* theAntiNeutronHEPModel = new G4HEAntiNeutronInelastic(); theAntiNeutronInelasticProcess->RegisterMe(theAntiNeutronLEPModel); theAntiNeutronInelasticProcess->RegisterMe(theAntiNeutronHEPModel); pManager->AddDiscreteProcess(theAntiNeutronInelasticProcess); G4AntiNeutronAnnihilationAtRest* theAntiNeutronAnnihilation = new G4AntiNeutronAnnihilationAtRest(); pManager->AddRestProcess(theAntiNeutronAnnihilation); } However, at low energy ( deuteron which has energy lower than 100 MeV) my program can not reproduce the experimental data about neutron yields and spectra. Speccially, with 160 MeV deuteron i can not see any neutron. It would be very grateful if you can let me know what is error in my program and suggest me some ideas to solve this problem. Thank you very much in advance! Phan Viet Cuong

Add Message

to: "confused on the results about neutron yield and spectra produced from reaction (d,n) "