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

None neutron_HP_Physics and NIST_Material 

Keywords: HP-neutron-physics ; NIST-material
Forum: Hadronic Processes
Date: 08 May, 2006
From: Sylvia Studeny <sylvia.studeny@gsf.de>

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

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   (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   (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: "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 ]