RE: [fluka-discuss]: K alpha of Copper and Bragg peak

From: Philippe Schoofs <philippe.schoofs_at_cern.ch>
Date: Wed, 14 Jun 2017 13:07:01 +0000

x-rays of course, sorry for the typo

Philippe

From: owner-fluka-discuss_at_mi.infn.it [mailto:owner-fluka-discuss_at_mi.infn.it] On Behalf Of Philippe Schoofs
Sent: Wednesday, June 14, 2017 2:33 PM
To: Javier Lozano Díaz <jlozano_at_clpu.es>; fluka-discuss_at_fluka.org
Subject: RE: [fluka-discuss]: K alpha of Copper and Bragg peak

Hello Javier,

The reason you don’t see the r-rays you are looking for is because they are not explicitly produced in the code. Their emission effect on energy loss is instead accounted for continuously throughout the particle travel in matter. If you look for PIXE (particle-induced X-ray emission) you will find more information about how to simulate that. I’m no specialist though, I can’t really say much more than that.

Although you might not need them anymore, let me answer the rest of your questions :

· The type of USRBDX allows the user to specify if he or she wants to account for particles crossing in one or both directions, and if the binnings (energy and solid angle) should be logarithmic or linear.

· If you don’t specify a min energy and a max energy for your binning, then the default values are taken. In this case, you would have had a binning from 1MeV to 10 MeV (log) or from 0 to 10 MeV (linear). Note that in the case of a log scale, you would have missed the photons, even had they been explicitly produced in the code.

· In EMFCUT you specify the energy down to which particles are transported/produced (depending on SDUM)

For your Bragg peak your question is essentially : what is a USRBIN and what numbers does it return ?
You can find lots of info in the manual, under USRBIN. I will highlight here the following :
“5. Energy deposition will be expressed in GeV per cm3 per unit primary weight.”
There is also lots of info in the scoring lesson online
https://indico.cern.ch/event/604018/contributions/2435691/attachments/1395454/2242757/08_Scoring_2017.pdf


So the number you see on the y-axis if a measure in GeV per cm3 of the energy deposited into your bin(s). In order to represent it on a plot, you have to choose a way to “project” the data. You have multiple options you can choose from in flair (1D max, 1D projection,…), in the drop down list on the right hand side of the plot window and the numerical result on the plot will reflect your choice.


Hope this helps
Have a good day
Philippe



From: owner-fluka-discuss_at_mi.infn.it<mailto:owner-fluka-discuss_at_mi.infn.it> [mailto:owner-fluka-discuss_at_mi.infn.it] On Behalf Of Javier Lozano Díaz
Sent: Tuesday, June 13, 2017 1:06 PM
To: fluka-discuss_at_fluka.org<mailto:fluka-discuss_at_fluka.org>
Subject: [fluka-discuss]: K alpha of Copper and Bragg peak

Dear Fluka experts.

I have a problem with my simulations concerning the K alpha experiment in copeer and Bragg peak. They are the following:


K alpha of Copper: (see proton k alpha.flair and inp)

We send a monoenergetic proton beam of 10 MeV on a target of Copper. We want to see the k alpha (8.04 KeV). For to see the spectrum, we want to use the card USRBDX and the cards to activate the photoelectric effect: EMF, EMFFLUO, EMFCUT prod, EMFCUT transport. We want to put the threshold at 1 KeV to be sure to see the k alpha at 8.04 KeV and how many K alpha there will be.
1. Concerning URSBDX option (see attached file):
-We have a doubt about wich "type" to choose and the meaning of them.
-We choose the photon in the option "part".
-To use the card URSBDX, we create around the copper target 2 thin layers of vacuum in order to give the direction of view in the option "reg... to reg"
-We don't specify the min and max energy

2. Concerning EMFCUT:

- I have a doubt about the cut off, if it means the minimum or the maximum energy to be transported?
3. About the run:

It looks that we can see photon, but around 10 MeV and not 8 KeV.


Also, if we want to simulate a proton beam of 1 MeV, we don't see photon with URSBDX.


Bragg peak: (see BP10.flair and inp)

We want to see the bragg peak of a proton beam stopping into a scintillator. We want to see dE/dX according to dX. For this, we use the cards URSBIN, in 1D.
1. We are not sure about what means the Y axis (see picture BP10.png). More we decrease the area of the detector and more this value increase.In the picture, the peak goes until 0.07 in the Y axis, what it means ?


Thank you very much.

Kind regards.


--
D. Javier Lozano Díaz
Técnico de apoyo en Radioprotección
Centro de láseres Pulsados Ultracortos (CLPU)
Calle Adaja nº8 C.P. 37185 Villamayor, Salamanca.
Edificio M5, Parque Científico USAL.
Tel: 923 338 121 Ext: 26
E-mail: jlozano_at_clpu.es<mailto:jlozano_at_clpu.es>
www.clpu.es<http://www.clpu.es>

__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Wed Jun 14 2017 - 16:15:14 CEST

This archive was generated by hypermail 2.3.0 : Wed Jun 14 2017 - 16:15:18 CEST