Dear Masoomeh,
I'm not sure if I've understood your question, but maybe I could help a little about plotting USRBIN Error:
When you use a USRBIN card (for example, using unit 82 BIN), and plot the result using Flair/Plot, it will generate an output file named ...82_plot.dat (and not "... 82_plot.geo.dat"). In such file, there's DOSE and Error related to that interval you plotted (defined in "Projection & Limits")
This file present 4 columns: first and second with bin coordinates; third with DOSE, DOSE-EQ (type of particle/information you are recording at USRBIN) and fourth column present the error.
You can plot the error, for example, using MATLAB:
clc; clear; close all; fclose('all');
addpath('C:\[address of the file]')
MyData=load('My_file_plot_dat.txt'); %Import file (I use to convert to txt, but perhaps it's not necessary)
%% Input about number of bins (in my case, I'd plotted in X plane)
Ybin=301; %number of bins in Y
Zbin=201; %number of bins in Z
%%
Z=zeros(Zbin,1);
Y=zeros(Ybin,1);
for i=1:Zbin
for j=1:Ybin
Z(i)=MyData((i-1)*Ybin+j,1); %position Z
Y(j)=MyData((i-1)*Ybin+j,2); %position Y
DOSE(i,j)=MyData((i-1)*Ybin+j,3); %dose-eq in pSv/primary weight
DOSE(i,j)=FS*DOSE(i,j); %dose-Eq em uSv/h, where FS is my conversion factor
ERROR(i,j)=MyData((i-1)*Ybin+j,4);
end
end
figure('color','white');
surf(Y,Z,ERROR,'EdgeColor','none'); %plot ERROR
set(gca, 'ZScale', 'log');
h=colorbar;
-----Mensagem original-----
De: owner-fluka-discuss_at_mi.infn.it [mailto:owner-fluka-discuss_at_mi.infn.it] Em nome de Masoomeh Yarmohammadi
Enviada em: segunda-feira, 7 de março de 2016 06:49
Para: Vittorio Boccone <dr.vittorio.boccone_at_ieee.org>
Cc: fluka-discuss <fluka-discuss_at_fluka.org>
Assunto: Re: [fluka-discuss]: Statistical error
Dear Vittorio,
Thanks again,
True, I mentioned the error a bit general by writing "Less that 10%". while I was interested to the error in the position of x= -100 +_ 0.5 cm, y= 137 +_ 0.5 cm, z= 549 +_ 0.5 cm, that is the area that I set for VALUE Mesh 1*1*1.
Could you please let me know, if the error and dose data, point by point or bin by bin, is available in a file that I can open and plot it?
I will be grateful if you give your suggestion about the post that I had with the title of "Biasing and Redusing the statistical error".
Best regards
Masoomeh Yarmohammadi Satri
----- Original Message -----
From: "Vittorio Boccone" <dr.vittorio.boccone_at_ieee.org>
To: "Masoomeh Yarmohammadi" <myarmohammadi_at_ipm.ir>
Cc: "fluka-discuss" <fluka-discuss_at_fluka.org>
Sent: Sunday, March 6, 2016 5:19:59 PM
Subject: Re: [fluka-discuss]: Statistical error
Dear Masoomeh,
as you compare two different bins with two different bin sizes you will have different errors.
In the las mail you mentioned 2% vs 24%, now 10% vs 24. Your comparison is very qualitative and in order to be more precise you should tell me at which position you are doing the comparison, as the first USRBIN does not contain the second one (5 cm downstream) and the statistic is very poor (just plot a slice of your USRBIN to see the effects…). Still bin will have different sizes and are from different places.
Cartesian binning n. 1 "DEq-Be-WF " , generalized particle n. 240
X coordinate: from -1.6000E+02 to 5.0000E+01 cm, 100 bins ( 2.1000E+00 cm wide)
Y coordinate: from 0.0000E+00 to 2.5000E+02 cm, 100 bins ( 2.5000E+00 cm wide)
Z coordinate: from 5.3500E+02 to 5.4500E+02 cm, 100 bins ( 1.0000E-01 cm wide)
Cartesian binning n. 2 "Value " , generalized particle n. 240
X coordinate: from -1.0050E+02 to -9.9500E+01 cm, 1 bins ( 1.0000E+00 cm wide)
Y coordinate: from 1.3650E+02 to 1.3750E+02 cm, 1 bins ( 1.0000E+00 cm wide)
Z coordinate: from 5.4850E+02 to 5.4950E+02 cm, 1 bins ( 1.0000E+00 cm wide)
Data follow in a matrix A(ix,iy,iz), format (1(5x,1p,10(1x,e11.4)))
From my side there’s absolutely no mismatch, just different bin sizes with different statistics in different places.
Best,
V.
On 06 Mar 2016, at 13:17, Masoomeh Yarmohammadi < myarmohammadi_at_ipm.ir > wrote:
It seems this averaging for flair plotting and reading about 1 bin are not in agreement about the error, less that 10% and 24% respectively. But they are in a good agreement about Does Eq. (Photon scoring) values. Does it correct?!?
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id�c_info
(Ƨځ맲r*'~&x%Zm 0nߖi'iw
Received on Tue Mar 08 2016 - 12:48:31 CET