Message: Re: wrong momentum distribution of secondary after neutron capture Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Idea Re: wrong momentum distribution of secondary after neutron capture 

Forum: Hadronic Processes
Re: None wrong momentum distribution of secondary after neutron capture (Michael Gericke)
Date: 04 Sep, 2014
From: C Mos <C Mos>

This won't fix your main problem but here are some thoughts:

A sensitive detector gets activated when a particles moves within the volume. In your case, every time the triton or proton are inside the SD, it will call the ProcessHit() function. In other words, every time either particle does a step, it call the function and you will be recording the vertex momentum and vertex position (and other parameters) more than once. I could suggest yo check the process before you store the hit. Something like this:

  const G4VProcess* process = aStep->GetPostStepPoint()->GetProcessDefinedStep();
  if(process->GetProcessName() == "nCapture"){

     [COLLECT vertex info here]

  }

Thinking about it, something is funny indeed about your plot as it is not isotropic at all. You said you started the neutrons with energies in the MeV region. I quick test would be to use cold neutrons or to implement my suggestion above to see the momentum of neutrons only at capture events.

And yes, I agree with you. Reassigning momentum directions is not a good way to do things.

Cheers,

C Mos

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

1 Sad: Re: wrong momentum distribution of secondary after neutron capture   (Michael Gericke - 04 Sep, 2014)
 Add Message Add Message
to: "Re: wrong momentum distribution of secondary after neutron capture"

 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 ]