Go to the source code of this file.
Defines | |
| #define | DOUBLECHOOZ1 0 |
| #define | DOUBLECHOOZ2 1 |
| #define | DAYA_BAY 2 |
| #define | T2K 3 |
| #define | NOVALE 4 |
| #define | NOVAME 5 |
| #define | LAR5LE 6 |
| #define | LAR5ME 7 |
| #define | LAR5FEHO 8 |
| #define | LAR100FEHO 9 |
Functions | |
| const char * | experiment_name (int which) |
| void | set_run_name (const char *name) |
| void | set_test_hierachy (int which) |
| void | set_test_dcp (double dcp) |
| void | set_test_sstt13 (double sstt13) |
| void | set_test_theta23 (double th23) |
| const char * | run_name () |
| double | test_ldms () |
| double | test_sdms () |
| double | test_hierarchy () |
| double | test_theta12 () |
| double | test_theta23 () |
| double | test_theta13 () |
| double | test_dcp () |
| void | dump_rates () |
| void | th13dms_scan () |
| void | ssttdcp_scan () |
| void | hier_sens () |
|
|
Definition at line 6 of file nova-globes.h. Referenced by dump_rates(), experiment_name(), and th13dms_scan(). |
|
|
Definition at line 4 of file nova-globes.h. Referenced by dump_rates(), and experiment_name(). |
|
|
Definition at line 5 of file nova-globes.h. Referenced by dump_rates(), and experiment_name(). |
|
|
Definition at line 13 of file nova-globes.h. Referenced by dump_rates(), and experiment_name(). |
|
|
Definition at line 12 of file nova-globes.h. Referenced by dump_rates(), and experiment_name(). |
|
|
Definition at line 10 of file nova-globes.h. Referenced by dump_rates(), and experiment_name(). |
|
|
Definition at line 11 of file nova-globes.h. Referenced by dump_rates(), and experiment_name(). |
|
|
Definition at line 8 of file nova-globes.h. Referenced by dump_rates(), and experiment_name(). |
|
|
Definition at line 9 of file nova-globes.h. Referenced by dump_rates(), experiment_name(), and th13dms_scan(). |
|
|
Definition at line 7 of file nova-globes.h. Referenced by dump_rates(), experiment_name(), and th13dms_scan(). |
|
|
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 }
|
|
|
Definition at line 140 of file hier_sens.c. References chisqr(), gs_test_values, hier_sens_best_chisqr(), run_name(), SSTTHIGH, SSTTLOW, test_dcp(), test_hierarchy(), test_ldms(), test_sdms(), test_theta12(), test_theta13(), and test_theta23(). Referenced by main(). 00141 {
00142 int i, j, k;
00143 double theta12 = test_theta12();
00144 double theta23 = test_theta23();
00145 double theta13 = test_theta13();
00146 double deltacp = test_dcp();
00147 double sdms = test_sdms();
00148 double ldms = test_ldms();
00149 double sstt13;
00150 double dcp;
00151 double chisqr;
00152 double ssttmin;
00153 double dcpmin;
00154 double th23min;
00155 double chisqr_best;
00156 double sstt_best;
00157 double dcp_best;
00158
00159 glb_params true_values = glbAllocParams();
00160 if (test_hierarchy()>0.0) {
00161 glbDefineParams(true_values, theta12, theta13, theta23, deltacp,
00162 sdms, sdms+ldms);
00163 }
00164 else {
00165 glbDefineParams(true_values, theta12, theta13, theta23, deltacp,
00166 sdms, -ldms);
00167 }
00168 glbSetDensityParams(true_values, 1.0, GLB_ALL);
00169
00170 char fname[256];
00171 sprintf(fname,"%s_hier_sens.txt",run_name());
00172 FILE* fp = fopen(fname,"w");
00173
00174 /* Go to every point in sstt13-dcp space, compute new "true" rates
00175 * and then find the best solution using the opposite choice for the
00176 * mass hierarchy
00177 */
00178 gs_test_values = glbAllocParams();
00179 for (i=0; i<NSAMPLEX; ++i) {
00180 sstt13 = SSTTLOW+(SSTTHIGH-SSTTLOW)*((float)i+0.5)/NSAMPLEX;
00181 theta13 = 0.5*asin(sqrt(sstt13));
00182 for (j=0; j<NSAMPLEY; ++j) {
00183 dcp = 2.0*M_PI*((float)j+0.5)/NSAMPLEY;
00184
00185 glbSetOscParams(true_values, theta13, GLB_THETA_13);
00186 glbSetOscParams(true_values, dcp, GLB_DELTA_CP);
00187 if (test_hierarchy()>0.0) {
00188 glbSetOscParams(true_values, sdms+ldms, GLB_DM_31);
00189 }
00190 else {
00191 glbSetOscParams(true_values, -ldms, GLB_DM_31);
00192 }
00193 glbSetOscillationParameters(true_values);
00194 glbSetRates();
00195
00196 chisqr = hier_sens_best_chisqr(theta12,
00197 theta13,
00198 theta23,
00199 dcp,
00200 sdms,
00201 ldms,
00202 -test_hierarchy(),
00203 &ssttmin,
00204 &dcpmin,
00205 &th23min);
00206 fprintf(fp, "%f %f %f %f %f %f\n",
00207 sstt13, dcp/M_PI, chisqr, ssttmin, dcpmin, th23min);
00208 } /* j values of dcp */
00209 } /* i values of sstt13 */
00210 fclose(fp);
00211 }
|
|
|
Definition at line 40 of file test_points.c. Referenced by cp_sens(), hier_sens(), octant_sens(), and ssttdcp_scan(). 00040 { return gs_run_name; }
|
|
|
Definition at line 16 of file test_points.c. References gs_run_name. Referenced by main(). 00016 {
00017 if (strlen(name)>255) {
00018 strncpy(gs_run_name, name, 255);
00019 }
00020 else {
00021 strcpy(gs_run_name, name);
00022 }
00023 }
|
|
|
Definition at line 34 of file test_points.c. References gs_test_dcp. Referenced by main(). 00034 {
00035 gs_test_dcp = dcp;
00036 }
|
|
|
|
|
|
Definition at line 37 of file test_points.c. References gs_test_th13. Referenced by main(). 00037 {
00038 gs_test_th13 = 0.5*asin(sqrt(sstt13));
00039 }
|
|
|
Definition at line 31 of file test_points.c. References gs_test_th23. Referenced by main(). 00031 {
00032 gs_test_th23 = th23;
00033 }
|
|
|
Definition at line 14 of file ssttdcp_scan.c. References run_name(), SSTTHIGH, SSTTLOW, test_dcp(), test_hierarchy(), test_ldms(), test_sdms(), test_theta12(), test_theta13(), and test_theta23(). Referenced by main(). 00015 {
00016 int i, j, k;
00017 double sstt13;
00018 double thetheta13;
00019 double dcp;
00020 double resnh, resih;
00021 double resnhmin, resihmin;
00022 /* Define standard oscillation parameters */
00023 double theta12 = test_theta12();
00024 double theta23 = test_theta23();
00025 double theta13 = test_theta13();
00026 double deltacp = test_dcp();
00027 double sdm = test_sdms();
00028 double ldm = test_ldms();
00029 double sig_theta12 = 0.0;
00030 double sig_theta23 = 0.0;
00031 double sig_theta13 = 0.0;
00032 double sig_deltacp = 0.0;
00033 double sig_sdm = 0.0;
00034 double sig_ldm = 0.0;
00035
00036 /* Initialize parameter vector(s) */
00037 glb_params true_values = glbAllocParams();
00038 glb_params test_values = glbAllocParams();
00039 glb_params input_errors = glbAllocParams();
00040 glb_params minimum = glbAllocParams();
00041
00042 /* Define the true parameters with their uncertainties */
00043 if (test_hierarchy()>0.0) {
00044 glbDefineParams(true_values,
00045 theta12, theta13, theta23, deltacp,
00046 sdm, sdm+ldm);
00047 }
00048 else {
00049 glbDefineParams(true_values,
00050 theta12, theta13, theta23, deltacp,
00051 sdm, -ldm);
00052 }
00053 glbSetDensityParams(true_values,1.0,GLB_ALL);
00054
00055 /* Define errors on input parameters */
00056 glbDefineParams(input_errors,
00057 sig_theta12, sig_theta13, sig_theta23, sig_deltacp,
00058 sig_sdm, sqrt(sig_sdm*sig_sdm+sig_ldm*sig_ldm));
00059 glbSetDensityParams(input_errors, 0, GLB_ALL);
00060
00061 /* The simulated "true" data rates are computed */
00062 glbSetOscillationParameters(true_values);
00063 glbSetInputErrors(input_errors);
00064 glbSetRates();
00065
00066 /* Set up the test values */
00067 glbDefineParams(test_values,theta12,theta13,theta23,deltacp,sdm,sdm+ldm);
00068 glbSetDensityParams(test_values,1.0,GLB_ALL);
00069
00070 char fname[256];
00071 sprintf(fname,"%s_ssttdcp_scan.txt",run_name());
00072 FILE* fp = fopen(fname,"w");
00073
00074 /* Iteration over all values to be computed */
00075 for(i=0; i<NSAMPLEX; ++i) {
00076 sstt13 = SSTTLOW+(SSTTHIGH-SSTTLOW)*((float)i+0.5)/NSAMPLEX;
00077 thetheta13 = 0.5*asin(sqrt(sstt13));
00078 for(j=0; j<NSAMPLEY; ++j) {
00079 dcp = 2.0*M_PI*((float)j+0.5)/NSAMPLEY;
00080 resihmin = 1.E30;
00081 resnhmin = 1.E30;
00082 for (k=0; k<=30; ++k) {
00083 theta23 = (30.0+k)*M_PI/180.0;
00084
00085 glbSetOscParams(test_values,thetheta13,GLB_THETA_13);
00086 glbSetOscParams(test_values,dcp, GLB_DELTA_CP);
00087 glbSetOscParams(test_values,theta23, GLB_THETA_23);
00088
00089 /* Test normal hierarchy */
00090 glbSetOscParams(test_values,sdm+ldm, GLB_DM_31);
00091 resnh = glbChiSys(test_values, GLB_ALL, GLB_ALL);
00092
00093 glbSetOscParams(test_values, -ldm, GLB_DM_31);
00094 resih = glbChiSys(test_values, GLB_ALL, GLB_ALL);
00095
00096 /* Keep track of best solution */
00097 if (resih<resihmin) resihmin = resih;
00098 if (resnh<resnhmin) resnhmin = resnh;
00099 }
00100 fprintf(fp,"%f %f %f %f\n", sstt13,dcp/M_PI, resnhmin, resihmin);
00101 }
00102 }
00103 close(fp);
00104 }
|
|
|
Definition at line 47 of file test_points.c. Referenced by dump_rates(), hier_sens(), main(), octant_sens(), and ssttdcp_scan(). 00047 { return gs_test_dcp; }
|
|
|
Definition at line 41 of file test_points.c. Referenced by cp_sens(), dump_rates(), hier_sens(), main(), octant_sens(), and ssttdcp_scan(). 00041 { return gs_test_hierarchy; }
|
|
|
Definition at line 43 of file test_points.c. Referenced by cp_sens(), dump_rates(), hier_sens(), main(), octant_sens(), and ssttdcp_scan(). 00043 { return gs_test_ldms; }
|
|
|
Definition at line 42 of file test_points.c. Referenced by cp_sens(), dump_rates(), hier_sens(), octant_sens(), and ssttdcp_scan(). 00042 { return gs_test_sdms; }
|
|
|
Definition at line 44 of file test_points.c. Referenced by cp_sens(), dump_rates(), hier_sens(), octant_sens(), and ssttdcp_scan(). 00044 { return gs_test_th12; }
|
|
|
Definition at line 46 of file test_points.c. Referenced by cp_sens(), dump_rates(), hier_sens(), main(), octant_sens(), and ssttdcp_scan(). 00046 { return gs_test_th13; }
|
|
|
Definition at line 45 of file test_points.c. Referenced by cp_sens(), dump_rates(), hier_sens(), main(), octant_sens(), and ssttdcp_scan(). 00045 { return gs_test_th23; }
|
|
|
Definition at line 17 of file th13dms_scan.c. References DAYA_BAY, NOVAME, and T2K. 00018 {
00019 /* Define standard oscillation parameters */
00020 double theta12 = 35.0*M_PI/180.0;
00021 double theta23 = 45.0*M_PI/180.0;
00022 double theta13 = asin(sqrt(0.0))/2;
00023 double deltacp = 0.0;
00024 double sdm = 7.6e-5;
00025 double ldm = 2.5e-3;
00026
00027 /* Initialize parameter vector(s) */
00028 glb_params true_values = glbAllocParams();
00029 glb_params test_values = glbAllocParams();
00030
00031 glbDefineParams(true_values,theta12,theta13,theta23,deltacp,sdm,sdm+ldm);
00032 glbSetDensityParams(true_values,1.0,GLB_ALL);
00033
00034 glbDefineParams(test_values,theta12,theta13,theta23,deltacp,sdm,sdm+ldm);
00035 glbSetDensityParams(test_values,1.0,GLB_ALL);
00036
00037 /* The simulated data are computed */
00038 glbSetOscillationParameters(true_values);
00039 glbSetRates();
00040
00041 int i,j;
00042 double sstt13,thetheta13,dcp;
00043 double log10sstt;
00044 double resnh0, resnh1, resnh2;
00045 double resih0, resih1, resih2;
00046 double res0, res1, res2;
00047 double res01, res012;
00048 double best0, best1, best2;
00049 double best01, best012;
00050
00051 FILE* fp = fopen("th13dms_scan.txt","w");
00052 for (i=0; i<NSAMPLE; ++i) {
00053 log10sstt = -3.0+3.0*((float)i+0.5)/NSAMPLE;
00054 theta13 = 0.5*asin(sqrt(pow(10.0,log10sstt)));
00055 for (j=0; j<NSAMPLE; ++j) {
00056 ldm = 0.0+5.0E-3*((float)j+0.5)/NSAMPLE;
00057
00058 /* Set rates for this assumption of ldm */
00059 glbDefineParams(true_values,theta12,0.0,theta23,deltacp,sdm,sdm+ldm);
00060 glbSetDensityParams(true_values,1.0,GLB_ALL);
00061 glbSetOscillationParameters(true_values);
00062 glbSetRates();
00063
00064 /* Set vector of test values */
00065 glbSetOscParams(test_values,theta13,GLB_THETA_13);
00066 glbSetOscParams(test_values,sdm+ldm, GLB_DM_31);
00067
00068 resnh0 = glbChiSys(test_values,DAYA_BAY,GLB_ALL); // Reactor
00069 resnh1 = glbChiSys(test_values,T2K, GLB_ALL); // T2K
00070 resnh2 = glbChiSys(test_values,NOVAME, GLB_ALL); // NOvA-ME
00071
00072 /* Set rates for this assumption of ldm */
00073 glbDefineParams(true_values,theta12,0.0,theta23,deltacp,sdm,sdm-ldm);
00074 glbSetDensityParams(true_values,1.0,GLB_ALL);
00075 glbSetOscillationParameters(true_values);
00076 glbSetRates();
00077
00078 glbSetOscParams(test_values,theta13,GLB_THETA_13);
00079 glbSetOscParams(test_values,sdm-ldm, GLB_DM_31);
00080
00081 resih0 = glbChiSys(test_values,DAYA_BAY,GLB_ALL); // Reactor
00082 resih1 = glbChiSys(test_values,T2K, GLB_ALL); // T2K
00083 resih2 = glbChiSys(test_values,NOVAME, GLB_ALL); // NOvA-ME
00084
00085 fprintf(fp,"%f %f %f %f %f %f %f %f %f %f %f %f %f\n",
00086 log10sstt,ldm,
00087 resnh0, resnh1, resnh2,
00088 resih0, resih1, resih2,
00089 res0, res1, res2,
00090 res01, res012);
00091 }
00092 }
00093 close(fp);
00094 }
|
1.3.9.1