Re: how to selcet weight window

Hi Liang,

I think that currently in FLUKA you can actually go up to 100000 (not 10000, which is what says the manual). Also you can start with importances lower than one (0.0001 is supposed to be minimum, I don't know if this limit has been lowered too). Remember not to set an importance ratio greater than 5 from one sheath of your shielding to the contiguous ones.

In order to help you tune your weight windows, you could edit mgdraw.f and print to a file the weight (wtrack), and kinetic energy (etrack-am(jtrack)) of a given particle type (jtrack) at the crossing of a given region boundary (you have to set the condition within 'bxdraw' entry). You could then for instance use gnuplot to plot a weight-energy cloud (set logscale xy; plot 'file.dat' us 2:1) and from there decide the weight window as a function of the particle, region and energy. By the way, in order for mgdraw to work, remember to set the USERDUMP card in your input file, and, of course, recompile the code (lfluka,...) with your modified mgdraw.f routine.

Hope you can progress on your problem with this hints.

Mario

On Aug 17, 2011, at 2:37 AM, liang wrote:

> Dear Fluka experts,
> I want to learn setting weight window with FLUKA.So i simulate a simple problem;an isotropic neutron beam of 1.6Gev,spheral iron shielding(diameter is 6m).The input file is attached.
> Refer to the FLUKA course(http://www.fluka.org/free_download/course/portugal2010/Lectures);
> Firstly,i ran the job without any biasing and printed counter of every region.
>
> Secondly,analyzed counter and assign region importance according to inverse of the attenuation(inverse ratio of total number of particles entered the regions).So the importance value of "far" region will be very large,and exceed the maximum allowable value(10000 in flair).I can not assign larger importance value to "further" region.How should i do?
>
> Finally,analyz counter,select weight windows(WW-THRES,WW-FACTO) around average weights.But the average weights get from counters have nothing to do with energy range set by WW-THRES card?i mean the counters can not give average weight of the particles for some specified energy range?is my understanding correct? So how can i get average weights of the particles for different energy range?
>
> forgive my ignorant questions and thanks for any help!
> liang
>
>
>
>
> --=-smYthzq7h1bTfvkcFxaE
> Content-Disposition: attachment; filename="impor.inp"
> Content-Type: text/plain; name="impor.inp"; charset="UTF-8"
> Content-Transfer-Encoding: 7bit
>
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> TITLE
> Practice neutron weight window
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> DEFAULTS NEW-DEFA
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BEAM -1.6 10000.0 1.0NEUTRON
> BEAMPOS 0.0 0.0 0.0 0.0 0.0
> GEOBEGIN COMBNAME
> 0 0 A
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> RPP B1 -50000.0 +50000.0 -50000.0 +50000.0 -50000.0 +50000.0
> RPP B2 -20000.0 +20000.0 -20000.0 +20000.0 -20000.0 +20000.0
> * ********************
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B11 0.0 0.0 0.0 30.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B12 0.0 0.0 0.0 60.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B13 0.0 0.0 0.0 80.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B14 0.0 0.0 0.0 100.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B15 0.0 0.0 0.0 120.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B16 0.0 0.0 0.0 150.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B17 0.0 0.0 0.0 180.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B18 0.0 0.0 0.0 210.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B19 0.0 0.0 0.0 240.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B20 0.0 0.0 0.0 255.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B21 0.0 0.0 0.0 270.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B22 0.0 0.0 0.0 300.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B23 0.0 0.0 0.0 330.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B24 0.0 0.0 0.0 345.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B25 0.0 0.0 0.0 360.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B26 0.0 0.0 0.0 390.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B27 0.0 0.0 0.0 420.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B28 0.0 0.0 0.0 450.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B29 0.0 0.0 0.0 465.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B30 0.0 0.0 0.0 480.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B31 0.0 0.0 0.0 510.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B32 0.0 0.0 0.0 540.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B33 0.0 0.0 0.0 570.0
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> SPH B34 0.0 0.0 0.0 600.0
> END
> * ********************
> * black hole
> BH1 5 +B1 -B2
> * vacuum around
> VR2 5 +B2 -B34
> * *****************
> * SHILED
> S7 5 +B11
> S8 5 +B12 -B11
> S9 5 +B13 -B12
> S10 5 +B14 -B13
> S11 5 +B15 -B14
> S12 5 +B16 -B15
> S13 5 +B17 -B16
> S14 5 +B18 -B17
> S15 5 +B19 -B18
> S16 5 +B20 -B19
> S17 5 +B21 -B20
> S18 5 +B22 -B21
> S19 5 +B23 -B22
> S20 5 +B24 -B23
> S21 5 +B25 -B24
> S22 5 +B26 -B25
> S23 5 +B27 -B26
> S24 5 +B28 -B27
> S25 5 +B29 -B28
> S26 5 +B30 -B29
> S27 5 +B31 -B30
> S28 5 +B32 -B31
> S29 5 +B33 -B32
> S30 5 +B33 -B32
> END
> GEOEND
> * ***add tangang
> MATERIAL 7.85 26.0 TANGANG
> * **************add tangang
> COMPOUND -0.50 CARBON -99.50 IRON TANGANG
> * ************************
> * ***********add tangang
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> LOW-MAT IRON 26. -8. 296. IRON
> * ********************
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * External Black Hole
> ASSIGNMA BLCKHOLE BH1
> * Vacuum
> ASSIGNMA VACUUM VR2
> * ********ADD LARGE REGION SCORING TANGANG
> ASSIGNMA TANGANG S7 S29 1.0
> * ********************importance biasing*********
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.0001 S7 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.0001 S8 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.0003 S9 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.0004 S10 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.0008 S11 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.0016 S12 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.0064 S13 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.0219 S14 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.0438 S15 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.1752 S16 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 0.3872 S17 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 1.2056 S18 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 1.5311 S19 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 6.1244 S20 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 24.4976 S21 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 69.8181 S22 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 269.4305 S23 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 296.3736 S24 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 1238.8146 S25 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 5178.3581 S26 PRINT
> **********exceed maximum value
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 10000. S27 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 10000. S28 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1. 10000. S29 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1. 10000. S30 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 1.0 10000. VR2 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING -1.0 0.0 1.0 7.0 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING -1.0 0.0 10.0 16.0 PRINT
> * *****LOW NEUTRON BIASING****************
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0001 S7 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0001 S8 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0001 S9 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0001 S10 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0001 S11 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0001 S12 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0001 S13 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0001 S14 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0003 S15 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0003 S16 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0005 S17 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0007 S18 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0016 S19 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0024 S20 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0036 S21 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0090 S22 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0423 S23 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.0519 S24 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.1587 S25 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.6359 S26 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.6359 S27 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.6359 S28 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.6359 S29 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.6359 S30 PRINT
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BIASING 3.0 0.6359 VR2 PRINT
> * *******SCORING***********
> * REGION binning of the deposited energy WHOLE
> USRBIN 12. ENERGY -65. S30 deposit
> USRBIN BH1 1.0 &
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> RANDOMIZ 1.
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> START 1000000.0 0.0
> STOP
>
> --=-smYthzq7h1bTfvkcFxaE--
>
>
>
Received on Thu Aug 25 2011 - 09:26:47 CEST