Dear Mina
I use gnuplot to manipulate the 2 files and generate a combined plot with
the separate normalization applied to each one.
For example if you have 2 files (file1.dat and file2.dat) , each belonging
to separate runs with different projectiles with results in pSv/p, each
arranged in 4 columns of data ( coordinate 1, coordinate 2, dose in pSv/p,
%uncertainty), then in gnuplot, the command
splot '< paste file1.dat file2.dat' us 1:2: ($3*f1+$7*f2) notitle
will plot the sum of file1.dat and file2.dat with the 3rd column of the
files multiplied with f1 and f2 respectively. The 3rd column of second file
for gnuplot is the 7th data set.
f1 and f2 are constants for you, the conversion number for pSv/p to
microSv/h/unit beam current, which is a straight forward calculation.
Of course, before you issue the above command, a few more commands are
also required such as set pm3d map, set palette of your choice etc. Please
see the gnuplot manual for the syntax.
I get the file1.dat by plotting the data with flair. It then generates a
plot.dat which you rename. This is important since the data has to be in a
certain format for gnuplot to understand, and flair does this for you.
Repeat for the second data to get file2.dat. You can look up the flair
output window for more gnu commands and the color palette used, when it
creates a plot.
Please also see this link by Mario
http://www.fluka.org/web_archive/earchive/new-fluka-discuss/3119.html
You may also want to create a plot showing the combined uncertainties using
the error propagation technique.
Cheers Sunil
On Thu, Apr 2, 2015 at 2:50 AM, Mina Nozar <nozarm_at_triumf.ca> wrote:
> Hello everyone,
>
> I am working on a project that involves two targets and two beams in the
> same target station area. One beam is a 50 kW proton beam and another is a
> 500 kW electron beam. See the attached image please.
>
> This figure shows the prompt dose rates from the combined results of both
> sets of runs, one for the proton beam and one for the electron beam.
> However, as you can see, the results are shown in the raw units of
> pSv/prim. In order to show results in uSv/h ---> Sv/h, I would need to
> apply a different normalization to each set of binary files generated via
> usbsuw.
>
> Is there a way to do this mixed normalization at the usbsuw stage when
> merging the two sets of files from each run?
>
>
> Thanks in advance,
> Mina
>
>
>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Apr 02 2015 - 06:38:33 CEST