Message: G4ParticleHPInelastic can't not re-produce neutron energy spectrum from 7Li(p,n) reaction? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question G4ParticleHPInelastic can't not re-produce neutron energy spectrum from 7Li(p,n) reaction? 

Forum: Hadronic Processes
Date: 04 Apr, 2016
From: Phan Viet Cuong <Phan Viet Cuong>

Hello everybody,

I am trying to simulate the 7Li(p,n)7Be reaction which the input parameters as following: Ep=1912 keV, target LiF with the thickness of 2 mg/cm2 evaporated on 1 mm Cu (Cu backing). I used the G4ParticleHPInelastic model to simulate the interaction of proton at low energy. However, the results in neutron energy spectra obtained is not as the experimental ones.

Here is my code:

void ProtonPhysics::ConstructProcess()
{
   G4ProcessManager* pManager = G4Proton::Proton()->GetProcessManager();
   // Elastic

   G4HadronElasticProcess* pElastProc = new G4HadronElasticProcess();
   G4ParticleHPElastic* pElastModel = new G4ParticleHPElastic();
   //G4ParticleHPElasticData* pElastData = new G4ParticleHPElasticData();
   pElastProc->RegisterMe(pElastModel);
   //pElastProc->AddDataSet(pElastData);

   pManager->AddDiscreteProcess(pElastProc);

   // Inelastic
   G4ProtonInelasticProcess* inelProcess = new G4ProtonInelasticProcess();
   G4ParticleHPInelastic* thepIneModel = new G4ParticleHPInelastic(G4Proton::Proton());
   G4ParticleHPInelasticData* thepIneData = new G4ParticleHPInelasticData(G4Proton::Proton());
   thepIneData->SetMinKinEnergy(0.0);
   thepIneData->SetMaxKinEnergy(20.*MeV);
   inelProcess->RegisterMe(thepIneModel);
   inelProcess->AddDataSet(thepIneData);
   pManager->AddDiscreteProcess(inelProcess);

}

I also attaching here the comparison between the Geant4 simulation result and data. I was wondering if someone could give me some comments and advises for my case. Thank you very much!

With my best regards

Phan Viet Cuong

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/04/04/09.29-76376-SimComparison.png

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

1 None: Re: G4ParticleHPInelastic can't not re-produce neutron energy spectrum from 7Li(p,n) reaction?   (John Watts - 05 Apr, 2016)
1 Question: Re: G4ParticleHPInelastic can't not re-produce neutron energy spectrum from 7Li(p,n) reaction?   (Phan Viet Cuong - 05 Apr, 2016)
1 Feedback: Re: G4ParticleHPInelastic can't not re-produce neutron energy spectrum from 7Li(p,n) reaction?   (John Watts - 05 Apr, 2016)
... 2 Message(s)
2 Ok: Re: G4ParticleHPInelastic can't not re-produce neutron energy spectrum from 7Li(p,n) reaction?   (Phan Viet Cuong - 06 Apr, 2016)
2 None: Re: G4ParticleHPInelastic can't not re-produce neutron energy spectrum from 7Li(p,n) reaction?   (Artem Zontikov - 06 Apr, 2016)
1 Feedback: Re: G4ParticleHPInelastic can't not re-produce neutron energy spectrum from 7Li(p,n) reaction?   (John Watts - 06 Apr, 2016)
 Add Message Add Message
to: "G4ParticleHPInelastic can't not re-produce neutron energy spectrum from 7Li(p,n) reaction?"

 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 ]