Dear Fluka experts,

I’ve created an energy spectrum with USRTRACK for beta and bremsstrahlung radiation, as shown in the picture. I plotted the DX*Y quantity because as I understand, I will obtain the integrated fluence. I would like to obtain the average energy of both the radiations from this spectrum, with the relative error on the average. I’ve tried to get the average energy with a weighted mean on the fluence from the sum.lis file, but I’m missing how to obtain the average error (since I have errors for the fluences). Is it possible to take both these information directly from Fluka? If not, what is the easiest way to do it? Best regards,

Leonardo

Dear @leonardo.delpi,

Plotting the DX*Y you obtained the **differential** fluence, not the integrated one, but in the `*_sum.lis`

you could consult the integrated fluence.

FLUKA outputs the differential fluence, as requested, but does not provide directly the average value as you wish it using the `USRTRACK`

card, but using the `USRBIN`

card with a per `REGION`

scoring, where you would use the same region as for the `USRTRACK`

card. The output of the `USRBIN`

card will have the average error computed automatically.

You can nevertheless retrieve the error yourself. The errors are based on Monte Carlo statistics, meaning they follow a 1/sqrt(N) law, where N is the number of particles that have participated to the scoring (in each DX bin in your case). To retrieve the average error, you would have to:

Best,

Daniel

Dear Daniel, thank you very much for the response.

I still have some doubts. If DX*Y is the differential fluence, what is the difference between plotting Y or plotting DX*Y in a USRTRACK plot?

Can you also suggest me a program to acces the data from the sum.lis file? Usually I’ve used excel for the tab.lis file, but it seems like that it cannot properly divide the cells from the sum.lis file.

Thank you again,

Leonardo

Dear @leonardo.delpi,

I apologize for the misunderstanding. Indeed, `Y`

is the differencial fluence in GeV^{-1} \cdot cm^{-2}, while `DX*Y`

is the integral fluence over the bin energy range in cm^{-2}.

It is up to you what program you prefer to read the data files. Personally, I use coding scripts rather than programs, for example `python`

with the `pandas`

library.

Best,

Daniel

Dear Daniel,

thank you very much for your suggestion, I’ll soon try it.

Have a nice day,

Leonardo