Plotting USRBINs in 3.4.x: maximum and integrated value

Flair
1

Dear Flair experts,

with the display of maximum and integrated value of the plotted usrbin range gone in 3.4.0, is there a convenient alternative way to obtain or plot this information that doesn’t require me to parse the .dat file manually?

Best regards
Roman

3

Dear Roman,

the new flair 3.4.0 (announcement will follow soon) contains a different utility for plotting and manipulating the 3D mesh scoring, called meshtk, allowing:

  • to process FLUKA, Phits, MCNP and vtk 3D mesh files
  • list content and statistical information like max/min,int
  • unstructured meshes
  • make simple mathematical operations between meshes + - * /
  • rebin contents
  • make projections and maximum trace
  • convert forms

If you run meshtk it will print on the screen all commands how to invoke it

4

Dear Vasilis,

thanks, this sounds promising. Is it possible to run it from within Flair or just command line yet? Is the integral value already implemented (I don’t see an option nor the output)?

Roman

5

From command line for the moment
There should be a line with the Min/Max and Integral

$ meshtk usrbin.bnn
local> meshtk ./FCC/W/fcc_w_2_5_21.bnn
MeshData: 1 neutron
        Title:
        Type: cylindrical
        Symmetries:
        Score: 8
        Weight: 200000
        R: [0 .. 20] x 100 (0.2)
        Phi: [-3.14159 .. 3.14159] x 1 (6.28319)
        Z: [0 .. 50] x 200 (0.25)
        Min: 0  +/- 0
        Max: 0.661173  +/- 1
        Integral: 67.569
        Errors: Relative
MeshData: 2 energy
        Title:
        Type: cylindrical
        Symmetries:
        Score: 208
        Weight: 200000
        R: [0 .. 20] x 100 (0.2)
        Phi: [-3.14159 .. 3.14159] x 1 (6.28319)
        Z: [0 .. 50] x 200 (0.25)
        Min: 0  +/- 0
        Max: 5.28229  +/- 1
        Integral: 2.85293
        Errors: Relative
6

Trying around, I did get an integral for some .bnn files, but not for all, e.g.

>>/usr/local/flair/meshtk run_no_wire_30.bnn
MeshData: 1 EdepW1r
	Title: 
	Type: cylindrical
	Symmetries:
	Score: 208
	Weight: 20000
	R: [0 .. 0.005] x 1 (0.005)
	Phi: [-3.14159 .. 3.14159] x 1 (6.28319)
	Z: [-1 .. 1] x 21 (0.0952381)
	Min: 0  +/- 0
	Max: 24812  +/- 0.92743
	Errors: Relative

Also I noticed that the error for the Max (0.92743) is actually just the largest of all errors, not the error of the maximum value:

>>/usr/local/flair/meshtk -az run_no_wire_30.bnn -1
# Detector n: EdepW1r 1D Projection 1
-1 -0.9047619048 0 0
-0.9047619048 -0.8095238095 0 0
-0.8095238095 -0.7142857143 8.123466492 92.74299741
-0.7142857143 -0.619047619 43.08355331 52.23546028
-0.619047619 -0.5238095238 157.6155853 36.61586642
-0.5238095238 -0.4285714286 767.6334839 11.11223921
-0.4285714286 -0.3333333333 3027.743652 7.608804107
-0.3333333333 -0.2380952381 7304.853027 4.592140019
-0.2380952381 -0.1428571429 14791.41309 2.067569271
-0.1428571429 -0.04761904762 21949.80273 2.427963912
-0.04761904762 0.04761904762 24812.00781 1.22949928
0.04761904762 0.1428571429 22343.13086 2.004500665
0.1428571429 0.2380952381 13533.95703 2.376495488
0.2380952381 0.3333333333 7450.392578 3.533706069
0.3333333333 0.4285714286 2873.555664 4.713323712
0.4285714286 0.5238095238 627.9781494 20.2527836
0.5238095238 0.619047619 296.0669861 19.46048439
0.619047619 0.7142857143 56.44048309 56.65596128
0.7142857143 0.8095238095 0 0
0.8095238095 0.9047619048 0 0
0.9047619048 1 0 0
1 1 0 0

run_no_wire_30.bnn.txt (436 Bytes)

7

Thanks @rmartin
indeed there was a typo in the meshdata module which did not report back correctly the integral.
flair was just updated to the reflect this correction. Please re-install it