edgeset.cxx Source File

00001 // -*- C++ -*-  $Id: edgeset.cxx,v 1.1 2000/09/01 20:20:19 burnett Exp $
00002 //
00003 // This file is part of Gismo 2
00004 //
00005 #include "Brempr.h"
00006 /*
00007 "                                K-EDGE VERSION -- 27 JUL 1988/2400"
00008 "******************************************************************"
00009 " SUBPROGRAM TO SET UP PARAMETERS FOR K-FLUORESCENCE TRANSPORT     "
00010 " IN ELEMENTS WITH Z=1 THROUGH 100                                 "
00011 "******************************************************************"
00012 " Programmers:  W. R. Nelson and T. M. Jenkins (SLAC)              "
00013 "               W. S. Lockman (U. C. Santa Cruz)                   "
00014 "               Keith Weaver (UCSF) made all elements available    "
00015 "               D.W.O.Rogers (NRC) removed IVALID check, made for  "
00016 "                      a sequence of regions plus cosmetics        "
00017 "******************************************************************"
00018 "                                                                  "
00019 "  This routine must be called after HATCH, after IEDGEFL has been "
00020 "  filled and before the first call to SHOWER                      "
00021 "                                                                  "
00022 " The flag IEDGFL(IR) contains the atomic number of the medium for "
00023 " each region (IR), or it contains zero in order to disable K-edge "
00024 " fluorescence for that region.  The cross sections and energies   "
00025 " are from:  E. Storm and H. I. Israel, At. Data and Nucl. Data    "
00026 "            Tables 7 (1970) 565.                                  "
00027 " K-edge fluorescent yields are from:  Lederer et al, Table of the "
00028 "                                      Isotopes (6th Edition).     "
00029 "******************************************************************"
00030 
00031  Adapted for Gismo by T. Burnett (U of W).
00032  Here it is a member function for KEdge, just to set the data.
00033 
00034 */
00035 //"------------------------------------------------------------------"
00036 //"Average photon energies (EALF) (L -> K electron transitions).     "
00037 //"                  (Table 5 of Storm and Israel)                   "
00038 //"------------------------------------------------------------------"
00039 
00040 float ealf[]={
00041           0.   , 0.   , 0.054, 0.109, 0.184, 0.297, 0.393, 0.524,
00042           0.675, 0.849, 1.041, 1.255, 1.487, 1.739, 2.014, 2.307,
00043           2.622, 2.957, 3.312, 3.690, 4.088, 4.508, 4.949, 5.411,
00044           5.895, 6.400, 6.925, 7.472, 8.041, 8.631, 9.243, 9.876,
00045          10.532,11.210,11.907,12.630,13.375,14.142,14.933,15.746,
00046          16.584,17.443,18.327,19.235,20.167,21.122,22.103,23.108,
00047          24.138,25.192,26.272,27.378,28.510,29.667,30.851,32.062,
00048          33.297,34.564,35.858,37.179,38.528,39.906,41.313,42.750,
00049          44.218,45.714,47.242,48.801,50.392,52.104,53.670,55.356,
00050          57.078,58.832,60.620,62.443,64.303,66.200,68.133,70.103,
00051          72.113,74.161,76.246,78.378,80.547,82.757,85.018,87.321,
00052          89.662,92.050,94.491,96.977,99.516,102.10,104.74,107.44,
00053          110.20,113.03,115.93,118.89};
00054 //"------------------------------------------------------------------"
00055 //"Average photon energies (EBET) (M -> K electron transitions).     "
00056 //"                  (Table 5 of Storm and Israel)                   "
00057 //"------------------------------------------------------------------"
00058 float ebet[]={
00059           0.   , 0.   , 0.   , 0.   , 0.   , 0.   , 0.   , 0.   ,
00060           0.   , 0.   , 0.   , 0.   , 0.   , 1.838, 2.142, 2.468,
00061           2.817, 3.191, 3.589, 4.012, 4.459, 4.931, 5.427, 5.947,
00062           6.492, 7.059, 7.649, 8.265, 8.907, 9.572,10.263,10.983,
00063          11.730,12.503,13.300,14.126,14.980,15.859,16.767,17.700,
00064          18.661,19.648,20.785,21.705,22.778,23.878,25.008,26.166,
00065          27.354,28.573,29.821,31.103,32.414,33.757,35.131,36.535,
00066          37.966,39.431,40.930,42.460,44.024,45.622,47.253,48.918,
00067          50.618,52.352,54.123,55.930,57.772,59.652,61.572,63.531,
00068          65.529,67.564,69.642,71.759,73.919,76.121,78.367,80.656,
00069          82.991,85.370,87.796,90.273,92.794,95.365,97.989,100.66,
00070          103.39,106.17,109.01,111.90,114.85,117.86,120.93,124.08,
00071          127.29,130.58,133.96,137.41};
00072 //"------------------------------------------------------------------"
00073 //"Fluorescent yields (OMEG) (probability to get fluorescent photon  "
00074 //" per K-shell vacancy. (Table 10 of Lederer et al)                 "
00075 //"------------------------------------------------------------------"
00076 float omeg[]={
00077           0.    ,0.    ,0.    ,0.    ,0.    ,0.    ,0.    ,0.    ,
00078           0.    ,0.    ,0.    ,0.    ,0.0357,0.0470,0.0604,0.0761,
00079           0.0942,0.115 ,0.138 ,0.163 ,0.190 ,0.219 ,0.250 ,0.282 ,
00080           0.314 ,0.347 ,0.381 ,0.414 ,0.445 ,0.479 ,0.510 ,0.540 ,
00081           0.567 ,0.596 ,0.622 ,0.646 ,0.669 ,0.691 ,0.711 ,0.730 ,
00082           0.748 ,0.764 ,0.779 ,0.793 ,0.807 ,0.819 ,0.830 ,0.840 ,
00083           0.850 ,0.859 ,0.867 ,0.875 ,0.882 ,0.889 ,0.895 ,0.901 ,
00084           0.906 ,0.911 ,0.915 ,0.920 ,0.924 ,0.928 ,0.931 ,0.934 ,
00085           0.937 ,0.940 ,0.943 ,0.945 ,0.948 ,0.950 ,0.952 ,0.954 ,
00086           0.956 ,0.957 ,0.959 ,0.961 ,0.962 ,0.963 ,0.964 ,0.966 ,
00087           0.967 ,0.968 ,0.969 ,0.970 ,0.971 ,0.972 ,0.972 ,0.973 ,
00088           0.974 ,0.975 ,0.975 ,0.976 ,0.977 ,0.977 ,0.978 ,0.978 ,
00089           0.979 ,0.979 ,0.980 ,0.980 };
00090 //"------------------------------------------------------------------"
00091 //"Probability for the removal of a K-electron by PE effect (PHOTOK) "
00092 //" =(K-shell PE cross section)/(total PE cross section).            "
00093 //"                  (Table 8 of Storm and Israel)                   "
00094 //"------------------------------------------------------------------"
00095 float photok[]={
00096             1.   , 1.   , 1.   , 1.   , 1.   , 1.   , 1.   , 1.   ,
00097             1.   , 1.   , 0.928, 0.922, 0.916, 0.911, 0.907, 0.903,
00098             0.900, 0.896, 0.893, 0.890, 0.888, 0.885, 0.883, 0.880,
00099             0.878, 0.876, 0.874, 0.872, 0.871, 0.870, 0.869, 0.868,
00100             0.866, 0.865, 0.863, 0.861, 0.860, 0.858, 0.856, 0.855,
00101             0.853, 0.852, 0.850, 0.848, 0.846, 0.845, 0.843, 0.841,
00102             0.840, 0.839, 0.838, 0.836, 0.835, 0.833, 0.832, 0.830,
00103             0.829, 0.827, 0.825, 0.824, 0.823, 0.821, 0.819, 0.818,
00104             0.816, 0.815, 0.814, 0.813, 0.812, 0.811, 0.809, 0.807,
00105             0.805, 0.803, 0.802, 0.800, 0.798, 0.796, 0.794, 0.792,
00106             0.790, 0.788, 0.786, 0.784, 0.782, 0.780, 0.778, 0.776,
00107             0.774, 0.772, 0.770, 0.767, 0.765, 0.762, 0.759, 0.756,
00108             0.754, 0.751, 0.749, 0.746};
00109 //"------------------------------------------------------------------"
00110 //"Relative probability of a K-alpha photon emission (PKA)           "
00111 //" =(K-alpha intensity)/(K-alpha intensity + K-beta intensity).     "
00112 //"                  (Table 6 of Storm and Israel)                   "
00113 //"------------------------------------------------------------------"
00114 float pka[]={
00115          1.   , 1.   , 1.   , 1.   , 1.   , 1.   , 1.   , 1.   ,
00116          1.   , 1.   , 1.   , 1.   , 1.   , 1.   , 1.   , 0.955,
00117          0.940, 0.924, 0.914, 0.904, 0.901, 0.898, 0.897, 0.895,
00118          0.894, 0.893, 0.8925,0.892, 0.891, 0.890, 0.887, 0.884,
00119          0.879, 0.875, 0.870, 0.866, 0.862, 0.859, 0.855, 0.852,
00120          0.849, 0.847, 0.844, 0.842, 0.840, 0.838, 0.835, 0.833,
00121          0.831, 0.829, 0.826, 0.824, 0.821, 0.819, 0.817, 0.815,
00122          0.813, 0.812, 0.811, 0.810, 0.809, 0.808, 0.807, 0.806,
00123          0.805, 0.804, 0.803, 0.802, 0.801, 0.800, 0.799, 0.798,
00124          0.797, 0.795, 0.794, 0.794, 0.793, 0.792, 0.791, 0.790,
00125          0.789, 0.788, 0.787, 0.786, 0.785, 0.784, 0.783, 0.782,
00126          0.781, 0.781, 0.780, 0.779, 0.778, 0.778, 0.777, 0.776,
00127          0.775, 0.774, 0.773, 0.773 };
00128 
00129 
00130 void
00131 Kedge::edgeSetup(int iz)
00132 {
00133   iz--;       // first entry is Z=1
00134   iedgfl = 1;
00135   bkr1 = omeg[iz] * photok[iz];
00136   bkr2 = bkr1 * pka[iz];
00137   ekalph = ealf[iz] * 1.e-3;
00138   ekbeta = ebet[iz] * 1.e-3;
00139 }
00140