00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039 #ifndef ONE_D_DISTRIB_CLASS
00040 #define ONE_D_DISTRIB_CLASS
00041
00042 #include "distrib_base.hpp"
00043 #include "math_util.hpp"
00044 #include "messaging.hpp"
00045 #include "ENDF_file.hpp"
00046
00047
00048 class one_d_table : public distrib_base
00049
00050 {
00051 private:
00052
00053 int interp_;
00054 double b_;
00055 double x_0_;
00056 double y_0_;
00057
00058 public:
00059
00060
00061 double f(double E);
00062
00063 one_d_table(): interp_(0), b_(0.0), x_0_(0.0), y_0_(0.0) {}
00064 ~one_d_table() {}
00065
00066
00067 void expand_interp( double Max_E );
00068
00069
00070 void copy( one_d_table& list_2 );
00071
00072
00073 double line_multiple( );
00074
00075
00076 void list_interp(double e_in, one_d_table& left_list,
00077 one_d_table& right_list);
00078
00079
00080 void read_data( int NP, ENDF_file& inFile, double max_E );
00081 };
00082
00083
00084 class one_d_evap : public distrib_base
00085
00086 {
00087 private:
00088 double U_;
00089 double Theta_;
00090
00091 public:
00092
00093 double f(double E);
00094
00095 one_d_evap();
00096 ~one_d_evap() {}
00097
00098
00099 void initiate(double e_in, double Theta, double U);
00100
00101
00102 void copy( one_d_evap& list_2 );
00103
00104
00105 void extrema();
00106
00107
00108 void list_interp(double e_in, one_d_evap& left_list,
00109 one_d_evap& right_list);
00110 };
00111
00112
00113 class one_d_Maxwell : public distrib_base
00114
00115 {
00116 private:
00117 double U_;
00118 double Theta_;
00119
00120 public:
00121
00122 double f(double E);
00123
00124 one_d_Maxwell();
00125 ~one_d_Maxwell() {}
00126
00127
00128 void initiate(double e_in, double Theta, double U);
00129
00130
00131 void extrema();
00132
00133
00134 void copy( one_d_Maxwell& list_2 );
00135
00136
00137 void list_interp(double e_in,
00138 one_d_Maxwell& left_list, one_d_Maxwell& right_list);
00139 };
00140
00141
00142 class one_d_Watt : public distrib_base
00143
00144 {
00145 private:
00146 double a_;
00147 double b_;
00148 double U_;
00149
00150 public:
00151
00152
00153 int interp_type_b;
00154
00155 double f(double E);
00156
00157 one_d_Watt();
00158 ~one_d_Watt() {}
00159
00160
00161 void set_EUa( double e_in, double U, double a );
00162
00163
00164 void intermediate_b(double e_in, double U, double a, int a_interp,
00165 double b, int b_interp);
00166
00167
00168 void same_b( double b, int b_interp );
00169
00170
00171 void initiate();
00172
00173
00174
00175
00176
00177 void extrema();
00178
00179
00180 void copy( one_d_Watt& list_2 );
00181
00182
00183 void list_interp(double e_in,
00184 one_d_Watt& left_list, one_d_Watt& right_list);
00185 };
00186
00187
00188 double Watt_max(double beta, double dummy);
00189
00190
00191 class one_d_Madland : public distrib_base
00192
00193 {
00194 private:
00195 double EFL_;
00196 double EFH_;
00197 double TM_;
00198 double E_max;
00199 double g( double E_out, double E_frag );
00200 double Gamma_a;
00201
00202 public:
00203
00204 double f(double E)
00205 {
00206 return 0.5*( g(E, EFL_) + g(E, EFH_));
00207 }
00208
00209 one_d_Madland();
00210 ~one_d_Madland() {}
00211
00212
00213 void initiate(double e_in, double EFL, double EFH, double TM);
00214
00215
00216 void copy( one_d_Madland& list_2 );
00217
00218
00219 void list_interp(double e_in,
00220 one_d_Madland& left_list, one_d_Madland& right_list);
00221 };
00222
00223
00224 class one_d_Legendre : public distrib_base
00225
00226 {
00227 public:
00228 Param coef;
00229
00230 double f(double mu);
00231
00232 one_d_Legendre():coef(1,0.0) {}
00233 ~one_d_Legendre() {}
00234
00235
00236 void expand();
00237
00238
00239 void copy( one_d_Legendre& list_2 );
00240
00241
00242
00243 void list_interp(double e_in,
00244 one_d_Legendre& left_list, one_d_Legendre& right_list );
00245 };
00246
00247
00248 class one_d_isotropic : public distrib_base
00249
00250 {
00251 public:
00252
00253 double f(double E)
00254 {
00255 SevereError("one_d_isotropic::f","Write this function!");
00256 return 0.0;
00257 }
00258
00259 one_d_isotropic() {}
00260 ~one_d_isotropic() {}
00261
00262
00263 void initiate(double e_in);
00264
00265
00266 void copy( one_d_isotropic& list_2 );
00267 };
00268
00269
00270 class one_d_hat : public distrib_base
00271
00272 {
00273 private:
00274
00275 double left_scale;
00276 double top_scale;
00277 double right_scale;
00278
00279 public:
00280
00281 double f(double E_out);
00282
00283 one_d_hat() {}
00284 ~one_d_hat() {}
00285
00286
00287 void initiate(double left_end, double top, double right_end);
00288
00289
00290 void list_interp(double e_in, one_d_hat& left_list,
00291 one_d_hat& right_list)
00292 {
00293 SevereError("one_d_hats::list_interp",
00294 "Write code to interpolate between one_d_hats?");
00295 }
00296 };
00297
00298 #endif