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
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061 #include <stdio.h>
00062 #include <string.h>
00063 #include <time.h>
00064
00065
00066 #include "midas.h"
00067 #include "experim.h"
00068 #include "analyzer.h"
00069
00070
00071 #ifdef OS_WINNT
00072 #define VISUAL_CPLUSPLUS
00073 #endif
00074 #ifdef __linux__
00075 #define f2cFortran
00076 #endif
00077
00078 #ifdef HAVE_HBOOK
00079 #include <cfortran.h>
00080 #include <hbook.h>
00081
00082 PAWC_DEFINE(1000000);
00083 #endif
00084
00085
00086
00087
00088 char *analyzer_name = "Analyzer";
00089
00090
00091 INT analyzer_loop_period = 0;
00092
00093
00094 INT odb_size = DEFAULT_ODB_SIZE;
00095
00096
00097 RUNINFO runinfo;
00098 GLOBAL_PARAM global_param;
00099 EXP_PARAM exp_param;
00100 TRIGGER_SETTINGS trigger_settings;
00101
00102
00103
00104 extern ANA_MODULE scaler_accum_module;
00105 extern ANA_MODULE adc_calib_module;
00106 extern ANA_MODULE adc_summing_module;
00107
00108 ANA_MODULE *scaler_module[] = {
00109 &scaler_accum_module,
00110 NULL
00111 };
00112
00113 ANA_MODULE *trigger_module[] = {
00114 &adc_calib_module,
00115 &adc_summing_module,
00116 NULL
00117 };
00118
00119
00120
00121 ASUM_BANK_STR(asum_bank_str);
00122
00123 BANK_LIST ana_trigger_bank_list[] = {
00124
00125
00126 {"ADC0", TID_WORD, N_ADC, NULL},
00127 {"TDC0", TID_WORD, N_TDC, NULL},
00128
00129
00130 {"CADC", TID_FLOAT, N_ADC, NULL},
00131 {"ASUM", TID_STRUCT, sizeof(ASUM_BANK), asum_bank_str},
00132
00133 {""},
00134 };
00135
00136 BANK_LIST ana_scaler_bank_list[] = {
00137
00138 {"SCLR", TID_DWORD, N_ADC, NULL},
00139
00140
00141 {"ACUM", TID_DOUBLE, N_ADC, NULL},
00142 {""},
00143 };
00144
00145
00146
00147 ANALYZE_REQUEST analyze_request[] = {
00148 {"Trigger",
00149 {1,
00150 TRIGGER_ALL,
00151 GET_SOME,
00152 "SYSTEM",
00153 TRUE,
00154 "", "",}
00155 ,
00156 NULL,
00157 trigger_module,
00158 ana_trigger_bank_list,
00159 1000,
00160 TRUE,
00161 }
00162 ,
00163
00164 {"Scaler",
00165 {2,
00166 TRIGGER_ALL,
00167 GET_ALL,
00168 "SYSTEM",
00169 TRUE,
00170 "", "",}
00171 ,
00172 NULL,
00173 scaler_module,
00174 ana_scaler_bank_list,
00175 100,
00176 }
00177 ,
00178
00179 {""}
00180 ,
00181 };
00182
00183
00184
00185 INT analyzer_init()
00186 {
00187 HNDLE hDB, hKey;
00188 char str[80];
00189
00190 RUNINFO_STR(runinfo_str);
00191 EXP_PARAM_STR(exp_param_str);
00192 GLOBAL_PARAM_STR(global_param_str);
00193 TRIGGER_SETTINGS_STR(trigger_settings_str);
00194
00195
00196 cm_get_experiment_database(&hDB, NULL);
00197 db_create_record(hDB, 0, "/Runinfo", strcomb(runinfo_str));
00198 db_find_key(hDB, 0, "/Runinfo", &hKey);
00199 if (db_open_record(hDB, hKey, &runinfo, sizeof(runinfo), MODE_READ, NULL, NULL) !=
00200 DB_SUCCESS) {
00201 cm_msg(MERROR, "analyzer_init", "Cannot open \"/Runinfo\" tree in ODB");
00202 return 0;
00203 }
00204
00205 db_create_record(hDB, 0, "/Experiment/Run Parameters", strcomb(exp_param_str));
00206 db_find_key(hDB, 0, "/Experiment/Run Parameters", &hKey);
00207 if (db_open_record(hDB, hKey, &exp_param, sizeof(exp_param), MODE_READ, NULL, NULL) !=
00208 DB_SUCCESS) {
00209 cm_msg(MERROR, "analyzer_init",
00210 "Cannot open \"/Experiment/Run Parameters\" tree in ODB");
00211 return 0;
00212 }
00213
00214 sprintf(str, "/%s/Parameters/Global", analyzer_name);
00215 db_create_record(hDB, 0, str, strcomb(global_param_str));
00216 db_find_key(hDB, 0, str, &hKey);
00217 if (db_open_record
00218 (hDB, hKey, &global_param, sizeof(global_param), MODE_READ, NULL,
00219 NULL) != DB_SUCCESS) {
00220 cm_msg(MERROR, "analyzer_init", "Cannot open \"%s\" tree in ODB", str);
00221 return 0;
00222 }
00223
00224 db_create_record(hDB, 0, "/Equipment/Trigger/Settings", strcomb(trigger_settings_str));
00225 db_find_key(hDB, 0, "/Equipment/Trigger/Settings", &hKey);
00226
00227 if (db_open_record
00228 (hDB, hKey, &trigger_settings, sizeof(trigger_settings), MODE_READ, NULL,
00229 NULL) != DB_SUCCESS) {
00230 cm_msg(MERROR, "analyzer_init",
00231 "Cannot open \"/Equipment/Trigger/Settings\" tree in ODB");
00232 return 0;
00233 }
00234
00235 return SUCCESS;
00236 }
00237
00238
00239
00240 INT analyzer_exit()
00241 {
00242 return CM_SUCCESS;
00243 }
00244
00245
00246
00247 INT ana_begin_of_run(INT run_number, char *error)
00248 {
00249 return CM_SUCCESS;
00250 }
00251
00252
00253
00254 INT ana_end_of_run(INT run_number, char *error)
00255 {
00256 FILE *f;
00257 time_t now;
00258 char str[256];
00259 int size;
00260 double n;
00261 HNDLE hDB;
00262 BOOL flag;
00263
00264 cm_get_experiment_database(&hDB, NULL);
00265
00266
00267
00268 size = sizeof(flag);
00269 db_get_value(hDB, 0, "/Logger/Write data", &flag, &size, TID_BOOL, TRUE);
00270 if (flag && runinfo.online_mode == 1) {
00271
00272 size = sizeof(str);
00273 str[0] = 0;
00274 db_get_value(hDB, 0, "/Logger/Data Dir", str, &size, TID_STRING, TRUE);
00275 if (str[0] != 0)
00276 if (str[strlen(str) - 1] != DIR_SEPARATOR)
00277 strcat(str, DIR_SEPARATOR_STR);
00278 strcat(str, "runlog.txt");
00279
00280 f = fopen(str, "a");
00281
00282 time(&now);
00283 strcpy(str, ctime(&now));
00284 str[10] = 0;
00285
00286 fprintf(f, "%s\t%3d\t", str, runinfo.run_number);
00287
00288 strcpy(str, runinfo.start_time);
00289 str[19] = 0;
00290 fprintf(f, "%s\t", str + 11);
00291
00292 strcpy(str, ctime(&now));
00293 str[19] = 0;
00294 fprintf(f, "%s\t", str + 11);
00295
00296 size = sizeof(n);
00297 db_get_value(hDB, 0, "/Equipment/Trigger/Statistics/Events sent", &n, &size,
00298 TID_DOUBLE, TRUE);
00299
00300 fprintf(f, "%5.1lfk\t", n / 1000);
00301 fprintf(f, "%s\n", exp_param.comment);
00302
00303 fclose(f);
00304 }
00305
00306 return CM_SUCCESS;
00307 }
00308
00309
00310
00311 INT ana_pause_run(INT run_number, char *error)
00312 {
00313 return CM_SUCCESS;
00314 }
00315
00316
00317
00318 INT ana_resume_run(INT run_number, char *error)
00319 {
00320 return CM_SUCCESS;
00321 }
00322
00323
00324
00325 INT analyzer_loop()
00326 {
00327 return CM_SUCCESS;
00328 }
00329
00330