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