/********************************************************************\

Name: analyzer.c

Created by: Roger Fearick

Contents: System part of Analyzer code for K600

\********************************************************************/

/* standard includes */

#include <stdio.h>

#include <string.h>

#include <time.h>

/* midas includes */

#include "midas.h"

#include "experim.h"

#include "analyzer.h"

/* cernlib includes */

#ifdef OS_WINNT

#define VISUAL_CPLUSPLUS

#endif

#ifdef __linux__

#define f2cFortran

#endif

/* root includes */

#include <TTree.h>

#include <TFile.h>

/* home-made includes */

#include "Parameters.h"

/* // what is this!?! I thought we are working with ROOT?

#ifdef HAVE_HBOOK

#include <cfortran.h>

#include <hbook.h>

PAWC_DEFINE(2000000);

#endif // HAVE_HBOOK

*/

/*-- Globals -------------------------------------------------------*/

// ODB structures : what is this doing?

RUNINFO runinfo;

EXP_PARAM exp_param;

//GLOBAL_PARAM globals; //these are useless here? I need them in f-plane.c, which is outside this scope

// The analyzer name (client name) as seen by other MIDAS clients

char *analyzer_name = "Analyzer";

// analyzer_loop is called with this interval in ms (0 to disable)

INT analyzer_loop_period = 0;

// default ODB size

INT odb_size = DEFAULT_ODB_SIZE;

// counters used in f-plane to look at bad tdc events: either empty or too large

int toolarge_tdc_event_counter=0;

int empty_tdc_counter=0;

int qdc_counter1=0, qdc_counter2=0;

int adc_counter1=0, adc_counter2=0;

int tdc_counter=0;

int scaler_counter=0;

int beamline_counter=0;

int trailer_TDCERR_counter=0;

int trailer_triglost_counter=0;

int trailer_bufoverflow_counter=0;

//int skipcounter=0;

//int runnr=0;

/*-- Module declarations -------------------------------------------*/

extern ANA_MODULE scaler_accum_module;

extern ANA_MODULE main_module;

extern ANA_MODULE qdc_module;

extern ANA_MODULE adc_module;

extern ANA_MODULE beamline_module;

ANA_MODULE *beam_module[] = {

&beamline_module,

NULL

};

ANA_MODULE *scaler_module[] = {

&scaler_accum_module,

NULL

};

ANA_MODULE *trigger_module[] = {

&adc_module,

&qdc_module,

&main_module,

NULL

};

/*-- Bank definitions ----------------------------------------------*/

BANK_LIST ana_trigger_bank_list[] = {

// online banks

{"ADC0", TID_DWORD, N_ADC, NULL},

{"QDC0", TID_DWORD, N_QDC, NULL},

{"TDC0", TID_DWORD, N_TDC, NULL},

// calculated banks

//{"CADC", TID_FLOAT, N_ADC, NULL},

//{"ASUM", TID_STRUCT, sizeof(ASUM_BANK), asum_bank_str},

{""},

};

BANK_LIST ana_scaler_bank_list[] = {

// online banks

{"SCLR", TID_DWORD, N_SCLR, NULL},

// calculated banks

// {"ACUM", TID_DOUBLE, N_ADC, NULL},

{""},

};

BANK_LIST ana_beam_bank_list[] = {

// online banks

{"MSRD", TID_FLOAT, N_MSRD, NULL},

{""},

};

/*-- Event request list --------------------------------------------*/

ANALYZE_REQUEST analyze_request[] = {

{"Trigger", /* equipment name */

{1, /* event ID */

TRIGGER_ALL, /* trigger mask */

GET_SOME, /* get some events */

"SYSTEM", /* event buffer */

TRUE, /* enabled */

"", "",}

,

NULL, /* analyzer routine */

trigger_module, /* module list */

ana_trigger_bank_list, /* bank list */

1000, /* RWNT buffer size */

TRUE, /* Use tests for this event */

}

,

{"Scaler", /* equipment name */

{2, /* event ID */

TRIGGER_ALL, /* trigger mask */

GET_ALL, /* get all events */

"SYSTEM", /* event buffer */

TRUE, /* enabled */

"", "",}

,

NULL, /* analyzer routine */

scaler_module, /* module list */

ana_scaler_bank_list, /* bank list */

100, /* RWNT buffer size */

}

,

{"Beamline", /* equipment name */

{3, /* event ID */

TRIGGER_ALL, /* trigger mask */

GET_ALL, /* get all events */

"SYSTEM", /* event buffer */

TRUE, /* enabled */

"", "",}

,

NULL, /* analyzer routine */

beam_module, /* module list */

ana_beam_bank_list, /* bank list */

100, /* RWNT buffer size */

}

,

{""}

,

};

/*-- Analyzer Init -------------------------------------------------*/

INT analyzer_init()

{

HNDLE hDB, hKey;

char str[80];

RUNINFO_STR(runinfo_str); // RUNINFO_STR in experim.h rn: not true !?!?!

EXP_PARAM_STR(exp_param_str); // EXP_PARAM_STR in experim.h rn: not true !!?!?!

//GLOBAL_PARAM_STR(global_param_str); // GLOBAL_PARAM_STR in experim.h

/* open ODB structures */

cm_get_experiment_database(&hDB, NULL);

db_create_record(hDB, 0, "/Runinfo", strcomb(runinfo_str));

db_find_key(hDB, 0, "/Runinfo", &hKey);

if (db_open_record(hDB, hKey, &runinfo, sizeof(runinfo), MODE_READ, NULL, NULL) !=

DB_SUCCESS) {

cm_msg(MERROR, "analyzer_init", "Cannot open \"/Runinfo\" tree in ODB");

return 0;

}

db_create_record(hDB, 0, "/Experiment/Run Parameters", strcomb(exp_param_str));

db_find_key(hDB, 0, "/Experiment/Run Parameters", &hKey);

if (db_open_record(hDB, hKey, &exp_param, sizeof(exp_param), MODE_READ, NULL, NULL) !=

DB_SUCCESS) {

cm_msg(MERROR, "analyzer_init",

"Cannot open \"/Experiment/Run Parameters\" tree in ODB");

return 0;

}

//printf("\n---testing testing 123----------------------------------------------- %i \n\n\n\n ",runinfo.run_number);

//runnr=runinfo.run_number;

ParameterInit();

return SUCCESS;

}

/*-- Analyzer Exit -------------------------------------------------*/

INT analyzer_exit()

{

return CM_SUCCESS;

}

/*-- Begin of Run --------------------------------------------------*/

INT ana_begin_of_run(INT run_number, char *error)

{

return CM_SUCCESS;

}

/*-- End of Run ----------------------------------------------------*/

INT ana_end_of_run(INT run_number, char *error)

{

FILE *f;

time_t now;

char str[256];

int size;

double n;

HNDLE hDB;

BOOL flag;

cm_get_experiment_database(&hDB, NULL);

/* update run log if run was written and running online */

size = sizeof(flag);

db_get_value(hDB, 0, "/Logger/Write data", &flag, &size, TID_BOOL, TRUE);

if (flag && runinfo.online_mode == 1) {

/* update run log */

size = sizeof(str);

str[0] = 0;

db_get_value(hDB, 0, "/Logger/Data Dir", str, &size, TID_STRING, TRUE);

if (str[0] != 0)

if (str[strlen(str) - 1] != DIR_SEPARATOR)

strcat(str, DIR_SEPARATOR_STR);

strcat(str, "runlog.txt");

f = fopen(str, "a");

time(&now);

strcpy(str, ctime(&now));

str[10] = 0;

fprintf(f, "%s\t%3d\t", str, runinfo.run_number);

strcpy(str, runinfo.start_time);

str[19] = 0;

fprintf(f, "%s\t", str + 11);

strcpy(str, ctime(&now));

str[19] = 0;

fprintf(f, "%s\t", str + 11);

size = sizeof(n);

db_get_value(hDB, 0, "/Equipment/Trigger/Statistics/Events sent", &n, &size,

TID_DOUBLE, TRUE);

fprintf(f, "%5.1lfk\t", n / 1000);

fprintf(f, "%s\n", exp_param.comment);

fclose(f);

}

printf("\n------------------------------------------------------");

printf("\n| |");

printf("\n| TDC Events: %i ",tdc_counter);

printf("\n| QDC Events: %i ",qdc_counter1);

printf("\n| ADC Events: %i ",adc_counter1);

printf("\n| Scaler Events %i ",scaler_counter);

// printf("\n| ");

// printf("\n| Events with skips %i ",skipcounter);

printf("\n| ");

printf("\n| Events with no TDC datawords : %i ",empty_tdc_counter);

printf("\n| Events with too many TDC datawords : %i ",toolarge_tdc_event_counter);

printf("\n| Events with no QDC datawords: %i ",qdc_counter2);

printf("\n| Events with no ADC datawords: %i ",adc_counter2);

printf("\n| ");

printf("\n| TDC chip ERRORS: %i ",trailer_TDCERR_counter); // counter in TDC trailer

printf("\n| TDC triggers lost: %i ",trailer_triglost_counter); // counter in TDC trailer

printf("\n| TDC buffer overflow: %i ",trailer_bufoverflow_counter); // counter in TDC trailer

printf("\n| |");

printf("\n------------------------------------------------------\n");

// close the root file (that containst the trees) created in f-plane.c

extern TFile *f1; // declared in main.c

extern TTree *t1; // declared in main.c

f1->cd();

f1->Write();

f1->Close();

return CM_SUCCESS;

}

/*-- Pause Run -----------------------------------------------------*/

INT ana_pause_run(INT run_number, char *error)

{

return CM_SUCCESS;

}

/*-- Resume Run ----------------------------------------------------*/

INT ana_resume_run(INT run_number, char *error)

{

return CM_SUCCESS;

}

/*-- Analyzer Loop -------------------------------------------------*/

INT analyzer_loop()

{

return CM_SUCCESS;

}

/*------------------------------------------------------------------*/