Message: Re: neutron_HP_Physics and NIST_Material Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: neutron_HP_Physics and NIST_Material 

Forum: Hadronic Processes
Re: None neutron_HP_Physics and NIST_Material (Sylvia Studeny)
Date: 08 May, 2006
From: Vladimir Ivantchenko <vnivanch@mail.cern.ch>

On Mon, 8 May 2006, Sylvia Studeny wrote:

> *** Discussion title: Hadronic Processes
> Email replies to PublicHyperNews@slac.stanford.edu must include:
>   In-Reply-To: <"/hadronprocess/491"@geant4-hn.slac.stanford.edu>
>   Subject: ...change this to be about your reply.
> 
> Hi,
> 
> I'm doing calculations with low energy neutrons. I started with example
> /runandevent/Re02 and changed a few things. All was running fine so far.
> Now I wanted to include physics which is better adapted to my problem.
> So I changed the neutron physics into:
> 
> // neutron Physics
>    pManager = G4Neutron::Neutron()->GetProcessManager();
> 
>    // add process
>    G4HadronElasticProcess* thenElasticProcess 
>                          = new G4HadronElasticProcess();
>    thenElasticProcess->AddDataSet(new G4NeutronHPElasticData); 
> 
>    G4LElastic* thenElasticModel = new G4LElastic();
>    thenElasticModel->SetMinEnergy(20.*MeV);
>    G4NeutronHPElastic* theNeutronHPElasticModel= new G4NeutronHPElastic();
>    theNeutronHPElasticModel ->SetMaxEnergy(20.*MeV);
>    thenElasticProcess->RegisterMe(thenElasticModel);
>    thenElasticProcess->RegisterMe(theNeutronHPElasticModel);
> 
>    pManager->AddDiscreteProcess(thenElasticProcess);
> 
>    G4NeutronInelasticProcess* theNeutronInelasticProcess 
>                          = new G4NeutronInelasticProcess();
>    theNeutronInelasticProcess->AddDataSet(new G4NeutronHPInelasticData);  
> 
>    G4NeutronHPInelastic* theNeutronHPInelasticModel= new G4NeutronHPInelastic();
>    theNeutronHPInelasticModel ->SetMaxEnergy(20.*MeV);
>    G4BinaryCascade* theNeutronBCModel = new G4BinaryCascade();
>    theNeutronBCModel->SetMinEnergy(20.*MeV);
>    theNeutronBCModel->SetMaxEnergy(6.*GeV);
>    G4LENeutronInelastic* theNeutronLEPModel = new G4LENeutronInelastic();
>    theNeutronLEPModel->SetMinEnergy(4.*GeV);
>    G4HENeutronInelastic* theNeutronHEPModel = new G4HENeutronInelastic();
>    theNeutronInelasticProcess->RegisterMe(theNeutronHPInelasticModel);
>    theNeutronInelasticProcess->RegisterMe(theNeutronBCModel);
>    theNeutronInelasticProcess->RegisterMe(theNeutronLEPModel);
>    theNeutronInelasticProcess->RegisterMe(theNeutronHEPModel);
>    pManager->AddDiscreteProcess(theNeutronInelasticProcess);
> 
>    G4HadronFissionProcess* thenFission
>                          = new G4HadronFissionProcess();
>    thenFission->AddDataSet(new G4NeutronHPFissionData);  
> 
>    G4LFission* thenFissionModel = new G4LFission();
>    thenFissionModel->SetMinEnergy(20.*MeV);
>    G4NeutronHPFission* theNeutronHPFissionModel= new G4NeutronHPFission();
>    theNeutronHPFissionModel ->SetMaxEnergy(20.*MeV);
>    thenFission->RegisterMe(thenFissionModel);
>    thenFission->RegisterMe(theNeutronHPFissionModel);
>    pManager->AddDiscreteProcess(thenFission);
> 
>    G4HadronCaptureProcess* thenCapture
>                          = new G4HadronCaptureProcess();
>    thenCapture->AddDataSet(new G4NeutronHPCaptureData);
> 
>    G4LCapture* thenCaptureModel = new G4LCapture();
>    thenCaptureModel->SetMinEnergy(20.*MeV);
>    G4NeutronHPCapture* theNeutronHPCaptureModel= new G4NeutronHPCapture();
>    theNeutronHPCaptureModel ->SetMaxEnergy(20.*MeV);
>    thenCapture->RegisterMe(thenCaptureModel);
>    thenCapture->RegisterMe(theNeutronHPCaptureModel);
>    pManager->AddDiscreteProcess(thenCapture);
> 
> (Basically I added the HP-Physics). Now I have the Problem, that as long
> as I use G4_Water as material all looks good. But when I change the
> material to G4_TISSUE_SOFT_ICRP (or some other ICRP-material), the
> tracks start looking very strange, the neutrons get reflected at the
> surface and there are a lot of electrons produced.
> 
> My question is: Did I make a mistake in the neutron-physics list? Or do
> I have to add something for the NIST-materials? In the end I wanted to
> use several materials out of the NIST-list, so it would be very
> convenient if there is a way to use it!
> 
> Thanks a lot, Sylvia
> 

Hello,

Material definition via NIST does not do something magic - you can dump 
the material list and see what is defined - it is equvivalent what you 
need. So, if you have problems, they are in physics.

I do not see any obvious problem in your physics list. Probably, you need 
to study few events and identify, which process provides artifitial 
reflections.

VI

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

1 None: Re: neutron_HP_Physics and NIST_Material   (Sylvia Studeny - 08 May, 2006)
(_ None: Re: neutron_HP_Physics and NIST_Material   (Vladimir Ivantchenko - 08 May, 2006)
(_ None: Re: neutron_HP_Physics and NIST_Material   (Sylvia Studeny - 08 May, 2006)
(_ None: Re: neutron_HP_Physics and NIST_Material   (Vladimir Ivantchenko - 08 May, 2006)
(_ None: Re: neutron_HP_Physics and NIST_Material   (Sylvia Studeny - 08 May, 2006)
(_ None: Re: neutron_HP_Physics and NIST_Material   (Sylvia Studeny - 08 May, 2006)
(_ None: Re: neutron_HP_Physics and NIST_Material   (Vladimir Ivantchenko - 08 May, 2006)
(_ None: Re: neutron_HP_Physics and NIST_Material   (Sylvia Studeny - 09 May, 2006)
 Add Message Add Message
to: "Re: neutron_HP_Physics and NIST_Material"

 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 ]