SM: Executing: make clean make -> creates replay (executing program) >replay >asking for seed <6000 > RANLUX LUXURY LEVEL SET BY RLUXGO : 3 P= 223 > RANLUX INITIALIZED BY RLUXGO FROM SEEDS 6000 0 0 > asking for Event number (how many events you want to process) <2000 //for the first time put some value No. of corrupted events = 0 / 0.% [micherdz@pc-micherdz SM]$ After this execution, you will find: sm.rz hist.dat fort.35 Make a look on fort.35 (in this case this will consists of 47230 lines: 8 1 0.677106351 28. 20 1 3.05867087 28. 21 2 1.55084824 28. 22 2 4.77554626 28. 23 2 4.58409405 28. 32 1 2.11453874 28. ... 499962 1 4.17611797 28. 499974 2 4.6806356 28. 499978 2 0.474533817 28. 499980 1 4.26317842 28. 499985 1 3.7901247 28. Where: first column - event number (i) second column - how many times neutrons scatters througt SM (j) third column - position in SM (x) fourth column - position in SM (z) (here allways 28 - end of SM in z dir.) if you make now look on host.dat, this will be probably empty, or by some luck you will find something here. With open fort.35 find some numnber of scattering neutrons on which you want to make closer look. Find for example 13 (maximal number of times neutrons catch SM plates). S now you know event number, for us this like 384486. Execute again replay: put the same seed: 6000 and as event number: 384486 After this hits.dat includes 14 lines, with x, y, z (?): 0. 1.8018069 0.013026433 2.1251425 1.82949149 -6.27246155 4.42777729 1.85418532 0.00842089762 6.73063768 1.87357793 -6.27706744 9.03371147 1.88766801 0.00381466283 11.3369864 1.89645427 -6.28167402 13.6404504 1.89993546 -0.000792271349 15.9440911 1.89811035 -0.0030960007 18.2478963 1.89097774 -0.00539990491 20.5518538 1.87853647 -0.00770398395 22.8559514 1.86078539 -0.0100082378 25.1601768 1.83772338 -0.0123126666 27.4645178 1.80934937 -0.0146172702 You can now execute pltest.kumac to see how this look like on SM plates. PAW>./pltest 14 And open sm.rz in PAW to see various things: PAW> h/file 1 sm.rz 8190 PAW> nt/pri 10 you will see: ****************************************************************** * Ntuple ID = 10 Entries = 47230 smnt ****************************************************************** * Var numb * Type * Packing * Range * Block * Name * ****************************************************************** * 1 * I*4 * * * SM * i * 2 * I*4 * * * SM * j * 3 * R*4 * * * SM * nlam * 4 * R*4 * * * SM * nx0 * 5 * R*4 * * * SM * ny0 * 6 * R*4 * * * SM * nz0 * 7 * R*4 * * * SM * nthx0 * 8 * R*4 * * * SM * nthy0 * 9 * R*4 * * * SM * nx * 10 * R*4 * * * SM * ny * 11 * R*4 * * * SM * nz * 12 * R*4 * * * SM * nthx * 13 * R*4 * * * SM * nthy ****************************************************************** * Block * Entries * Unpacked * Packed * Packing Factor * ****************************************************************** * SM * 47230 * 52 * 52 * 1.000 * * Total * --- * 52 * 52 * 1.000 * ****************************************************************** * Blocks = 1 Variables = 13 Columns = 13 * ****************************************************************** Now you can see given variable by: PAW> nt/plo 10.i (to see 1dim) PAW> nt/plo 10.nx%ny (to see 2dim) (see pictres attached) %%%%% In program %%%%%%% ranlux(vect,4) genetates vec(4) where: vec(1) vec(2) vec(3) vec(4) 0.615414619 0.657832861 0.299835086 0.039093554 0.696331561 0.784972012 0.867112398 0.662244439 0.171229184 0.167748809 0.307739496 0.460097671 0.924537897 0.225444436 0.775110841 0.682788491 0.784990609 0.054633677 0.508425713 0.0144108534 0.00689584017 0.327261806 0.954927027 0.665641129 0.74920541 0.235763729 0.856392801 0.226807892 0.678548872 0.584852457 0.387437761 0.711544752 0.336177468 0.723112166 0.312846541 0.374962628 0.8417238 0.186721981 0.310952723 0.466303468 0.759673297 0.997933149 0.280667365 0.0697531104