Ramesh Subedi 07/13/2005 This program runs in any linux machine. To compile this code just hit make It creates an executable vaxeff. TO run the compiled code issue the command ./vaxeff file.out where the file.dat is the input file and file.out is the output file. For example if the file.dat is 10by40.dat then file.out could be 10by40.out and the command to run the code would be ./vaxeff <10by40.dat> 10by40.out Your result is in the file.out file. For more information please see the file: letter.lis. To understand the meaning of each line of the file.dat file, please see the file instr.lis. Mandatory files to be present in the directory "eff" are: 1. vaxeff.f 2. mt19937.f 3. crosec.in 4. Makefile 5. file.dat *************************************************************************** The overall statistical accuracy of the efficiency data is about 3 to 5%. Better accuracy for high energy neutrons and lesser for low energy ones. *************************************************************************** How to read efficiency and kinetic energy from the output? See the following. Lets take a look of one small sample of output data. ############################################ a sample of output file ############################### 40 BY 4 BY 4 NE102 DETECTOR SENSE SWITCH SETTINGS AND NUMBER OF EVENTS TO DUMP (6I1,I4):NSW,MDUMP 122222 0. BIAS(MEVEE),TIMING(NSEC)BIN WIDTH,POSITION(IN) FOR STUDY(3F8.3):TBIAS,DTIM,DPOS 0.200 0.100 0.200 TWO VALUES OF ONE-PE LEVEL(MEVEE) FOR RES. SMEARING (2F6.2):ELOPE(2),ELOPE(3) 0.002 0.020 SEED FOR RANDOM NUMBER GENERATOR (I5):KS1,KS2 2223 0 DENSITY (GM/CM**3) AND RATIO OF H TO C ATOMS (2F8.3):RHO, COMP 1.032 1.104 LIGHT OUTPUT COEF. FOR PROTONS (6F12.5): PA,PB,PC,PR,PS,PT -8.00000000 -0.10000000 0.89999998 0.94999999 0.00000000 0.00000000 LIGHT OUTPUT COEF. FOR ALPHAS (6F12.5): AH,AO,AP,AQ,AU -5.90000010 -0.06500000 1.00999999 0.41000000 0.00000000 0.00000000 SCINTILLATOR BOUNDARIES(XB,YB,ZB all in INCHES),IGEO,IRANP (3F6.2,2I2) 20.00 2.00 4.00 0 1 POSITIONS DIRECTION COSINES ENERGY DELTA E COUNTS BINWIDTH -*-(XI(I),I=1,3) (CI(I),I=1,3) E1 DE NEV BINW 40" X 4" X 4" DEEP NE-102 DETECTOR RECTANGULAR SCINTILLATOR (INCHES) 4.000 DEEP BY 40.000 BY 4.000 40" X 4" X 4" DEEP NE-102 DETECTOR DENSITY= 1.032 GM/CM**3 H/C RATIO= 1.104 0.00 0.00 0.00 0.0000 0.0000 1.0000 40.00 0.000 100000. 0.100 ** RANDOM INITIAL POSITION ** INTEGRAL EFFICIENCY,ONE P.E.= 0.002 MEV BIN WIDTH 0.100 MEV NOTE FIRST BIN CORRESPONDS TO ZERO BIAS (i.e., no threshold). Energy = 40.00 with Bin Width= 0.100 output Number of Counts 310 42.067 23.852 22.870 22.569 22.359 22.123 21.942 21.729 21.547 21.361 21.181 21.012 20.839 20.656 20.496 20.334 20.144 19.957 19.761 19.532 19.334 19.149 18.946 18.722 18.508 18.318 18.109 17.909 17.683 17.440 17.227 17.026 16.827 16.601 16.393 16.165 15.946 15.730 15.514 15.297 15.074 14.850 14.621 14.418 14.218 14.022 13.821 13.624 13.445 13.284 13.114 12.909 12.772 12.610 12.421 12.274 12.125 11.955 11.810 11.675 11.549 11.399 11.262 11.121 10.993 10.863 10.767 10.634 10.513 10.398 10.277 10.167 10.055 9.954 9.842 9.736 9.633 9.540 9.436 9.337 9.241 9.166 9.092 9.020 8.950 8.879 8.804 8.732 8.638 8.576 8.495 8.398 8.320 8.246 8.167 8.093 8.025 7.954 7.873 7.790 7.719 7.644 7.556 7.481 7.418 7.341 7.287 7.211 7.147 7.074 6.981 6.906 6.841 6.757 6.675 6.604 6.538 6.470 6.390 6.327 6.256 6.190 6.121 6.050 5.985 5.920 5.850 5.784 5.724 5.664 5.606 5.552 5.485 5.430 5.377 5.312 5.245 5.201 5.154 5.101 5.054 5.009 4.973 4.917 4.878 4.842 4.792 4.754 4.710 4.678 4.645 4.613 4.577 4.549 4.525 4.497 4.464 4.428 4.395 4.371 4.334 4.299 4.281 4.257 4.236 4.220 4.197 4.179 4.155 4.135 4.109 4.091 4.075 4.053 4.033 4.016 3.985 3.953 3.920 3.904 3.880 3.856 3.838 3.812 3.791 3.773 3.757 3.733 3.706 3.687 3.664 3.647 3.619 3.591 3.571 3.547 3.512 3.498 3.472 3.455 3.431 3.413 3.391 3.361 3.339 3.316 3.296 3.273 3.251 3.223 3.197 3.161 3.135 3.105 3.085 3.058 3.030 3.001 2.971 2.940 2.907 2.880 2.853 2.819 2.802 2.777 2.756 2.741 2.719 2.703 2.678 2.653 2.620 2.592 2.564 2.530 2.502 2.478 2.441 2.406 2.377 2.335 2.301 2.269 2.243 2.212 2.192 2.157 2.116 2.084 2.053 2.028 2.003 1.967 1.940 1.908 1.870 1.837 1.813 1.792 1.756 1.729 1.690 1.653 1.600 1.560 1.527 1.484 1.460 1.424 1.390 1.342 1.308 1.279 1.250 1.205 1.157 1.128 1.089 1.054 1.013 0.970 0.934 0.895 0.859 0.816 0.788 0.745 0.704 0.664 0.606 0.548 0.510 0.474 0.439 0.398 0.362 0.314 0.281 0.246 0.213 0.170 0.135 0.106 0.073 0.044 0.023 0.010 0.004 0.003 ############################################################################### This portion of output data is for an incident neutron of kinetic energy 40.0 MeV. Here is how you read any chunk of data: In the INTEGRAL EFFICIENCY part of the data, the first entry (42.067) corresponds to zero bias (that means no threshold). The next entry right to this entry ,which is 23.852, is for 0.1 MeVee bias (threshold) since the BIN WIDTH is 0.100 MeV. Each next entry to the right is for another 0.1 extra MeVee threshold. Once the line of 10 entry is finished, the same process repeats for the next line.This scheme (left to right and below) continues until the last entry, which is almost a zero value of efficiency, is read. Each entry for efficiency is read in units of percent. All chunk of data for light output is also read in the same fashion. ************************************************************************************** More explanation is here: In this example,the efficiency for a threshold of 4.0 MeVee would be 15.074 which is the first entry of the row number 4 below. The efficiency for a thrsshold of 4.3 MeVee would be 14.418 which is the column number 3 of the same row. Each row after the first row increments the threshold by 1 MeVee.Each column after the first column increments the threshold by 0.1 MeVee. Column-> 0 1 2 3 4 5 6 7 8 9 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0 42.067 23.852 22.870 22.569 22.359 22.123 21.942 21.729 21.547 21.361 1 21.181 21.012 20.839 20.656 20.496 20.334 20.144 19.957 19.761 19.532 2 19.334 19.149 18.946 18.722 18.508 18.318 18.109 17.909 17.683 17.440 3 17.227 17.026 16.827 16.601 16.393 16.165 15.946 15.730 15.514 15.297 4 15.074 14.850 14.621 14.418 14.218 14.022 13.821 13.624 13.445 13.284 5 13.114 12.909 12.772 12.610 12.421 12.274 12.125 11.955 11.810 11.675 6 11.549 11.399 11.262 11.121 10.993 10.863 10.767 10.634 10.513 10.398 7 10.277 10.167 10.055 9.954 9.842 9.736 9.633 9.540 9.436 9.337 8 9.241 9.166 9.092 9.020 8.950 8.879 8.804 8.732 8.638 8.576 9 8.495 8.398 8.320 8.246 8.167 8.093 8.025 7.954 7.873 7.790 10 7.719 7.644 7.556 7.481 7.418 7.341 7.287 7.211 7.147 7.074 11 6.981 6.906 6.841 6.757 6.675 6.604 6.538 6.470 6.390 6.327 12 6.256 6.190 6.121 6.050 5.985 5.920 5.850 5.784 5.724 5.664 13 5.606 5.552 5.485 5.430 5.377 5.312 5.245 5.201 5.154 5.101 14 5.054 5.009 4.973 4.917 4.878 4.842 4.792 4.754 4.710 4.678 15 4.645 4.613 4.577 4.549 4.525 4.497 4.464 4.428 4.395 4.371 16 4.334 4.299 4.281 4.257 4.236 4.220 4.197 4.179 4.155 4.135 17 4.109 4.091 4.075 4.053 4.033 4.016 3.985 3.953 3.920 3.904 18 3.880 3.856 3.838 3.812 3.791 3.773 3.757 3.733 3.706 3.687 19 3.664 3.647 3.619 3.591 3.571 3.547 3.512 3.498 3.472 3.455 20 3.431 3.413 3.391 3.361 3.339 3.316 3.296 3.273 3.251 3.223 21 3.197 3.161 3.135 3.105 3.085 3.058 3.030 3.001 2.971 2.940 22 2.907 2.880 2.853 2.819 2.802 2.777 2.756 2.741 2.719 2.703 23 2.678 2.653 2.620 2.592 2.564 2.530 2.502 2.478 2.441 2.406 24 2.377 2.335 2.301 2.269 2.243 2.212 2.192 2.157 2.116 2.084 25 2.053 2.028 2.003 1.967 1.940 1.908 1.870 1.837 1.813 1.792 26 1.756 1.729 1.690 1.653 1.600 1.560 1.527 1.484 1.460 1.424 27 1.390 1.342 1.308 1.279 1.250 1.205 1.157 1.128 1.089 1.054 28 1.013 0.970 0.934 0.895 0.859 0.816 0.788 0.745 0.704 0.664 29 0.606 0.548 0.510 0.474 0.439 0.398 0.362 0.314 0.281 0.246 30 0.213 0.170 0.135 0.106 0.073 0.044 0.023 0.010 0.004 0.003 ***************************************************************************************************