Message: Re: Poisson Noise Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Poisson Noise 

Forum: Particles
Re: Question Poisson Noise (A. Hod)
Date: 23 Dec, 2017
From: Tim Fulcher <Tim Fulcher>

Hi A.Hod

if you use the /gps source as opposed to the gun source in your macro file you'll get a great deal more flexibility. The /gps source will, supposedly*, allow you to determine the energy spectra manually. So you can generate an energy spectra with noise in, say, Matlab and punch it into your .mac file.

/gps/ene/type User
/gps/hist/type energy
/gps/hist/point 242.5 0.001129
/gps/hist/point 243.5 0.006832
/gps/hist/point 244.5 0.03135
/gps/hist/point 245.5 0.109
/gps/hist/point 246.5 0.2875
/gps/hist/point 247.5 0.5746
/gps/hist/point 248.5 0.8707
/gps/hist/point 259.5 1.0000
/gps/hist/point 250.5 0.8707
/gps/hist/point 251.5 0.5746
/gps/hist/point 252.5 0.2875
/gps/hist/point 253.5 0.109
/gps/hist/point 254.5 0.03135
/gps/hist/point 255.5 0.006832
/gps/hist/point 256.5 0.001129

If your baseline code uses /gun it isn't too difficult to cut over to /gps. Just a few hacks here and there.

*I'm inserting a disclaimer here. I'm having a helluva time trying to get a manually generated Gaussian energy distribution to match with the default Geant4 Gaussian distribution:

/gps/ene/type Gauss
/gps/ene/mono 250 MeV
/gps/ene/sigma 12.5 MeV

Regards

Tim

 Add Message Add Message
to: "Re: Poisson Noise"

 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 ]