Message: Forcing neutron to undergo inelastic scattering Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Forcing neutron to undergo inelastic scattering 

Forum: Event and Track Management
Date: 17 Oct, 2012
From: C Mos <C Mos>

Hello Neutrons from my source travel right to my target. As soon as it interacts in the target, I would like to force the neutron to undergo inelastic scattering. Since this is a local effect for testing purposes, I was thinking in modyfing the neutron reaction only inside the target (hence i did not disable processes in the PhysicsList class). Would it be allowed to change/set the process through my defined UserSteppingAction() class? Notice that I am planning to implement this behaviour only when neutrons are about the E threshold for such processes.Please see my code below:

Q2: In case I disable all processes except the inelastic processes for neutrons, what would happen after the neutron is below the energy thresold for any inelastic event? Would be the neutron be transported to the boundariesof the universe and killed?

void mySteppingAction::UserSteppingAction(const G4Step* aStep) {


//Only in target
if(aStep->GetTrack()->GetNextVolume()->GetName() == "PhysTarget") {
  G4int partPDGencoding=aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
  //Only when particle is a neutron
     //Act on neutrons only when they have at least the min. threshold energy
       if(aStep->GetTrack()->GetKineticEnergy() > 13. ){
#include "G4NeutronInelasticProcess.hh"
#include "G4LENeutronInelastic.hh"
         G4NeutronInelasticProcess* theInelasticProcess = new G4NeutronInelasticProcess("inelastic");
	 theInelasticProcess->RegisterMe(new G4LENeutronInelastic);	        



}// end of UserSteppingAction()

Thank you, C.

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

1 Warning: Re: Forcing neutron to undergo inelastic scattering   (Gumplinger Peter - 17 Oct, 2012)
 Add Message Add Message
to: "Forcing neutron to undergo inelastic scattering"

 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 ]