#include <stdio.h>#include <stdlib.h>#include <math.h>#include <globes/globes.h>#include "nova-globes.h"Go to the source code of this file.
Functions | |
| const char * | experiment_name (int which) |
| void | dump_rates () |
|
|
Definition at line 24 of file dump_rates.c. References DAYA_BAY, LAR100FEHO, LAR5FEHO, LAR5LE, LAR5ME, NOVALE, NOVAME, T2K, test_dcp(), test_hierarchy(), test_ldms(), test_sdms(), test_theta12(), test_theta13(), and test_theta23(). Referenced by main(). 00025 {
00026 /* Define standard oscillation parameters */
00027 double theta12 = test_theta12();
00028 double sdm = test_sdms();
00029 double theta23 = test_theta23();
00030 double ldm = test_ldms();
00031 double theta13 = test_theta13();
00032 double deltacp = test_dcp();
00033
00034 /* Initialize parameter vector(s) */
00035 glb_params true_values = glbAllocParams();
00036
00037 if (test_hierarchy()>0.0) {
00038 glbDefineParams(true_values,theta12,theta13,theta23,deltacp,sdm,sdm+ldm);
00039 }
00040 else {
00041 glbDefineParams(true_values,theta12,theta13,theta23,deltacp,sdm,-ldm);
00042 }
00043 glbSetDensityParams(true_values,1.0,GLB_ALL);
00044
00045 /* The simulated data are computed */
00046 glbSetOscillationParameters(true_values);
00047 glbSetRates();
00048
00049 FILE* fp = 0;
00050
00051 /* Daya bay rates */
00052 fp = fopen("dayabay_rates_sig.txt","w");
00053 glbShowRuleRates(fp, DAYA_BAY, 0,
00054 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00055 close(fp);
00056 fp = fopen("dayabay_rates_bg.txt","w");
00057 glbShowRuleRates(fp, DAYA_BAY, 0,
00058 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00059 close(fp);
00060
00061 /* T2K rates */
00062 fp = fopen("t2k_rates_sig_0.txt","w");
00063 glbShowRuleRates(fp, T2K, 0,
00064 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00065 close(fp);
00066 fp = fopen("t2k_rates_bg_0.txt","w");
00067 glbShowRuleRates(fp, T2K, 0,
00068 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00069 close(fp);
00070 fp = fopen("t2k_rates_sig_1.txt","w");
00071 glbShowRuleRates(fp, T2K, 1,
00072 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00073 close(fp);
00074 fp = fopen("t2k_rates_bg_1.txt","w");
00075 glbShowRuleRates(fp, T2K, 1,
00076 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00077 close(fp);
00078 fp = fopen("t2k_rates_sig_2.txt","w");
00079 glbShowRuleRates(fp, T2K, 2,
00080 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00081 close(fp);
00082 fp = fopen("t2k_rates_bg_2.txt","w");
00083 glbShowRuleRates(fp, T2K, 2,
00084 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00085 close(fp);
00086 fp = fopen("t2k_rates_sig_3.txt","w");
00087 glbShowRuleRates(fp, T2K, 3,
00088 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00089 close(fp);
00090 fp = fopen("t2k_rates_bg_3.txt","w");
00091 glbShowRuleRates(fp, T2K, 3,
00092 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00093 close(fp);
00094
00095 /* NOvA in ME beam */
00096 fp = fopen("novame_rates_sig_0.txt","w");
00097 glbShowRuleRates(fp, NOVAME, 0,
00098 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00099 close(fp);
00100 fp = fopen("novame_rates_bg_0.txt","w");
00101 glbShowRuleRates(fp, NOVAME, 0,
00102 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00103 fp = fopen("novame_rates_sig_1.txt","w");
00104 glbShowRuleRates(fp, NOVAME, 1,
00105 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00106 close(fp);
00107 fp = fopen("novame_rates_bg_1.txt","w");
00108 glbShowRuleRates(fp, NOVAME, 1,
00109 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00110 fp = fopen("novame_rates_sig_2.txt","w");
00111 glbShowRuleRates(fp, NOVAME, 2,
00112 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00113 close(fp);
00114 fp = fopen("novame_rates_bg_2.txt","w");
00115 glbShowRuleRates(fp, NOVAME, 2,
00116 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00117 fp = fopen("novame_rates_sig_3.txt","w");
00118 glbShowRuleRates(fp, NOVAME, 3,
00119 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00120 close(fp);
00121 fp = fopen("novame_rates_bg_3.txt","w");
00122 glbShowRuleRates(fp, NOVAME, 3,
00123 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00124
00125 /* NOvA in LE beam */
00126 fp = fopen("novale_rates_sig_0.txt","w");
00127 glbShowRuleRates(fp, NOVALE, 0,
00128 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00129 close(fp);
00130 fp = fopen("novale_rates_bg_0.txt","w");
00131 glbShowRuleRates(fp, NOVALE, 0,
00132 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00133 fp = fopen("novale_rates_sig_1.txt","w");
00134 glbShowRuleRates(fp, NOVALE, 1,
00135 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00136 close(fp);
00137 fp = fopen("novale_rates_bg_1.txt","w");
00138 glbShowRuleRates(fp, NOVALE, 1,
00139 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00140 fp = fopen("novale_rates_sig_2.txt","w");
00141 glbShowRuleRates(fp, NOVALE, 2,
00142 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00143 close(fp);
00144 fp = fopen("novale_rates_bg_2.txt","w");
00145 glbShowRuleRates(fp, NOVALE, 2,
00146 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00147 fp = fopen("novale_rates_sig_3.txt","w");
00148 glbShowRuleRates(fp, NOVALE, 3,
00149 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00150 close(fp);
00151 fp = fopen("novale_rates_bg_3.txt","w");
00152 glbShowRuleRates(fp, NOVALE, 3,
00153 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00154
00155 /* LAr5 in LE beam */
00156 fp = fopen("lar5le_rates_sig_0.txt","w");
00157 glbShowRuleRates(fp, LAR5LE, 0,
00158 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00159 close(fp);
00160 fp = fopen("lar5le_rates_bg_0.txt","w");
00161 glbShowRuleRates(fp, LAR5LE, 0,
00162 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00163 fp = fopen("lar5le_rates_sig_1.txt","w");
00164 glbShowRuleRates(fp, LAR5LE, 1,
00165 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00166 close(fp);
00167 fp = fopen("lar5le_rates_bg_1.txt","w");
00168 glbShowRuleRates(fp, LAR5LE, 1,
00169 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00170 fp = fopen("lar5le_rates_sig_2.txt","w");
00171 glbShowRuleRates(fp, LAR5LE, 2,
00172 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00173 close(fp);
00174 fp = fopen("lar5le_rates_bg_2.txt","w");
00175 glbShowRuleRates(fp, LAR5LE, 2,
00176 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00177 fp = fopen("lar5le_rates_sig_3.txt","w");
00178 glbShowRuleRates(fp, LAR5LE, 3,
00179 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00180 close(fp);
00181 fp = fopen("lar5le_rates_bg_3.txt","w");
00182 glbShowRuleRates(fp, LAR5LE, 3,
00183 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00184
00185 /* LAr5 in ME beam */
00186 fp = fopen("lar5me_rates_sig_0.txt","w");
00187 glbShowRuleRates(fp, LAR5ME, 0,
00188 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00189 close(fp);
00190 fp = fopen("lar5me_rates_bg_0.txt","w");
00191 glbShowRuleRates(fp, LAR5ME, 0,
00192 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00193 fp = fopen("lar5me_rates_sig_1.txt","w");
00194 glbShowRuleRates(fp, LAR5ME, 1,
00195 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00196 close(fp);
00197 fp = fopen("lar5me_rates_bg_1.txt","w");
00198 glbShowRuleRates(fp, LAR5ME, 1,
00199 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00200 fp = fopen("lar5me_rates_sig_2.txt","w");
00201 glbShowRuleRates(fp, LAR5ME, 2,
00202 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00203 close(fp);
00204 fp = fopen("lar5me_rates_bg_2.txt","w");
00205 glbShowRuleRates(fp, LAR5ME, 2,
00206 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00207 fp = fopen("lar5me_rates_sig_3.txt","w");
00208 glbShowRuleRates(fp, LAR5ME, 3,
00209 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00210 close(fp);
00211 fp = fopen("lar5me_rates_bg_3.txt","w");
00212 glbShowRuleRates(fp, LAR5ME, 3,
00213 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00214 close(fp);
00215
00216 /* LAr5 in FEHO beam */
00217 fp = fopen("lar5feho_rates_sig_0.txt","w");
00218 glbShowRuleRates(fp, LAR5FEHO, 0,
00219 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00220 close(fp);
00221 fp = fopen("lar5feho_rates_bg_0.txt","w");
00222 glbShowRuleRates(fp, LAR5FEHO, 0,
00223 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00224 fp = fopen("lar5feho_rates_sig_1.txt","w");
00225 glbShowRuleRates(fp, LAR5FEHO, 1,
00226 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00227 close(fp);
00228 fp = fopen("lar5feho_rates_bg_1.txt","w");
00229 glbShowRuleRates(fp, LAR5FEHO, 1,
00230 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00231 fp = fopen("lar5feho_rates_sig_2.txt","w");
00232 glbShowRuleRates(fp, LAR5FEHO, 2,
00233 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00234 close(fp);
00235 fp = fopen("lar5feho_rates_bg_2.txt","w");
00236 glbShowRuleRates(fp, LAR5FEHO, 2,
00237 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00238 fp = fopen("lar5feho_rates_sig_3.txt","w");
00239 glbShowRuleRates(fp, LAR5FEHO, 3,
00240 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00241 close(fp);
00242 fp = fopen("lar5feho_rates_bg_3.txt","w");
00243 glbShowRuleRates(fp, LAR5FEHO, 3,
00244 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00245 close(fp);
00246
00247 /* LAr100 in FEHO beam */
00248 fp = fopen("lar100feho_rates_sig_0.txt","w");
00249 glbShowRuleRates(fp, LAR100FEHO, 0,
00250 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00251 close(fp);
00252 fp = fopen("lar100feho_rates_bg_0.txt","w");
00253 glbShowRuleRates(fp, LAR100FEHO, 0,
00254 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00255 fp = fopen("lar100feho_rates_sig_1.txt","w");
00256 glbShowRuleRates(fp, LAR100FEHO, 1,
00257 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00258 close(fp);
00259 fp = fopen("lar100feho_rates_bg_1.txt","w");
00260 glbShowRuleRates(fp, LAR100FEHO, 1,
00261 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00262 fp = fopen("lar100feho_rates_sig_2.txt","w");
00263 glbShowRuleRates(fp, LAR100FEHO, 2,
00264 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00265 close(fp);
00266 fp = fopen("lar100feho_rates_bg_2.txt","w");
00267 glbShowRuleRates(fp, LAR100FEHO, 2,
00268 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00269 fp = fopen("lar100feho_rates_sig_3.txt","w");
00270 glbShowRuleRates(fp, LAR100FEHO, 3,
00271 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_SIG);
00272 close(fp);
00273 fp = fopen("lar100feho_rates_bg_3.txt","w");
00274 glbShowRuleRates(fp, LAR100FEHO, 3,
00275 GLB_ALL, GLB_W_EFF, GLB_W_BG, GLB_W_COEFF, GLB_BG);
00276 close(fp);
00277 }
|
|
|
Definition at line 7 of file dump_rates.c. References DAYA_BAY, LAR100FEHO, LAR5FEHO, LAR5LE, LAR5ME, NOVALE, NOVAME, and T2K. 00008 {
00009 switch (which) {
00010 case DAYA_BAY: return "daya_bay";
00011 case T2K: return "t2k";
00012 case NOVALE: return "nova_le";
00013 case NOVAME: return "nova_me";
00014 case LAR5LE: return "lar5_le";
00015 case LAR5ME: return "lar5_me";
00016 case LAR5FEHO: return "lar5_feho";
00017 case LAR100FEHO: return "lar100_feho";
00018 default: return "??";
00019 }
00020 }
|
1.3.9.1