Message: Intensity of multi-particle source Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Intensity of multi-particle source 

Forum: Particles
Date: 30 Nov, 2017
From: Rafal Lalik <Rafal Lalik>

Hi,

my first message here so welcome to everybody.

I want to create a particle source based on CR data. For that I created a pdf
distribution for each element between H and Fe, where the x axis of the
histogram is an energy and the y is a flux (see attached file). Now I am
adding these histograms to a source as follow (example for two first
elements)

/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/pos/centre 0. 3.5 0. cm 
/gps/pos/radius 7 cm
/gps/ang/type cos
/gps/particle proton
/gps/ene/type Arb
/gps/hist/type arb
/gps/hist/file cards_gps_CR/1H.txt
/gps/hist/inter Lin

/gps/source/add 1.0
/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/pos/centre 0. 3.5 0. cm 
/gps/pos/radius 7 cm
/gps/ang/type cos
/gps/particle ion
/gps/ion 2 4
/gps/ene/type Arb
/gps/hist/type arb
/gps/hist/file cards_gps_CR/4He.txt
/gps/hist/inter Lin


However, the command /gps/source/add <intensity> with a parameter 1 which
I set for every source normalizes distributions to the same intensity. As I
mentioned, my histograms already are scaled in flux so they are could be
interpreted as real PDF. Is it possible to add the source in a way that it
keeps the intensity from the histogram definition and is not renormalized? Or
the only way is to recalculate the integral of the histograms and normalize
them manually?

Regards,
Rafal

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2017/11/30/02.38-22334-CR_model.png

 Add Message Add Message
to: "Intensity of multi-particle source"

 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 ]