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 #include <stdio.h>
00049 #include <math.h>
00050
00051
00052 #include "midas.h"
00053 #include "experim.h"
00054 #include "analyzer.h"
00055
00056
00057 #include <TH1F.h>
00058 #include <TDirectory.h>
00059
00060 #ifndef PI
00061 #define PI 3.14159265359
00062 #endif
00063
00064
00065
00066 ADC_SUMMING_PARAM adc_summing_param;
00067
00068
00069
00070 INT adc_summing(EVENT_HEADER *, void *);
00071 INT adc_summing_init(void);
00072 INT adc_summing_bor(INT run_number);
00073
00074 ADC_SUMMING_PARAM_STR(adc_summing_param_str);
00075
00076 ANA_MODULE adc_summing_module = {
00077 "ADC summing",
00078 "Stefan Ritt",
00079 adc_summing,
00080 NULL,
00081 NULL,
00082 adc_summing_init,
00083 NULL,
00084 &adc_summing_param,
00085 sizeof(adc_summing_param),
00086 adc_summing_param_str,
00087 };
00088
00089
00090
00091 extern TDirectory *gManaHistsDir;
00092
00093 static TH1F *gAdcSumHist;
00094
00095
00096
00097 INT adc_summing_init(void)
00098 {
00099
00100 gAdcSumHist = (TH1F *) gManaHistsDir->GetList()->FindObject("ADCSUM");
00101
00102 if (gAdcSumHist == NULL)
00103 gAdcSumHist = new TH1F("ADCSUM", "ADC sum", 500, 0, 10000);
00104
00105 return SUCCESS;
00106 }
00107
00108
00109
00110 INT adc_summing(EVENT_HEADER * pheader, void *pevent)
00111 {
00112 INT i, j, n_adc;
00113 float *cadc;
00114 ASUM_BANK *asum;
00115
00116
00117 n_adc = bk_locate(pevent, "CADC", &cadc);
00118 if (n_adc == 0)
00119 return 1;
00120
00121
00122 bk_create(pevent, "ASUM", TID_STRUCT, &asum);
00123
00124
00125 asum->sum = 0.f;
00126 for (i = j = 0; i < n_adc; i++)
00127 if (cadc[i] > adc_summing_param.adc_threshold) {
00128 asum->sum += cadc[i];
00129 j++;
00130 }
00131
00132
00133 asum->average = j > 0 ? asum->sum / j : 0;
00134
00135
00136 gAdcSumHist->Fill(asum->sum, 1);
00137
00138
00139 bk_close(pevent, asum + 1);
00140
00141 return SUCCESS;
00142 }