Message: Adding hadronic models Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Adding hadronic models 

Forum: Hadronic Processes
Date: 19 Aug, 2010
From: Neoh <Neoh>

Hi,

I am learning to simulate deuteron beam on beryllium target at about 20MeV. For hadronic process I use physics list QGSP_BIC_HP.

I noticed G4LEDeuteronInelastic is already included in G4IonBinaryCascadePhysics.cc in QGSP_BIC_HP. But if I register again the model G4LEDeuteronInelastic in my physics list, there is a big difference in neutron production. It gives ~50% more neutrons in one test simulation of 5M deuteron beams, compared with using QGSP_BIC_HP alone (for hadronic process).

I am a little confused why is this so. Will the order of registering hadronic processes in the code makes any difference?

My physics list is similar to hadronic example h01. This is what I add in end of PhysicsList::ConstructProcess():

  // for low-energy deuteron physics
  G4ParticleDefinition *theDeuteron = G4Deuteron::DeuteronDefinition();
  G4ProcessManager *deuteronProcMan = theDeuteron->GetProcessManager();
  G4DeuteronInelasticProcess* deuteronInelasticProcess = new G4DeuteronInelasticProcess();
  G4LEDeuteronInelastic* model = new G4LEDeuteronInelastic();
  deuteronInelasticProcess->RegisterMe(model);
  deuteronProcMan->AddDiscreteProcess(deuteronInelasticProcess);

Thanks for your advice.

Neoh

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

1 Idea: Re: Adding hadronic models   (Vladimir Ivanchenko - 28 Aug, 2010)
(_ None: Re: Adding hadronic models   (Neoh - 30 Aug, 2010)
 Add Message Add Message
to: "Adding hadronic models"

 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 ]