forked from sPHENIX-Collaboration/online_distribution
/
caen_correction.cc
126 lines (108 loc) · 2.66 KB
/
caen_correction.cc
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
#include <iostream>
#include <iomanip>
#include <fstream>
#include "caen_correction.h"
#include "packet.h"
using namespace std;
caen_correction::caen_correction ( const char *calibdata)
{
int nsample[9];
_broken = 0;
int index;
ifstream IN;
int adccorr[9];
float timecorr;
IN.open(calibdata,ios_base::in);
int chip;
int i;
if ( ! IN.is_open())
{
_broken = 1;
}
else
{
for ( chip = 0; chip < 4; chip++)
{
for ( i = 0; i < 1024; i++)
{
IN >> index
>> adccorr[0]
>> adccorr[1]
>> adccorr[2]
>> adccorr[3]
>> adccorr[4]
>> adccorr[5]
>> adccorr[6]
>> adccorr[7]
>> adccorr[8]
>> timecorr
>> nsample[0]
>> nsample[1]
>> nsample[2]
>> nsample[3]
>> nsample[4]
>> nsample[5]
>> nsample[6]
>> nsample[7]
>> nsample[8];
base [index][chip*9 + 0] = adccorr[0];
base [index][chip*9 + 1]= adccorr[1];
base [index][chip*9 + 2]= adccorr[2];
base [index][chip*9 + 3]= adccorr[3];
base [index][chip*9 + 4]= adccorr[4];
base [index][chip*9 + 5]= adccorr[5];
base [index][chip*9 + 6]= adccorr[6];
base [index][chip*9 + 7]= adccorr[7];
base [index][chip*9 + 8]= adccorr[8];
timevec[index][chip] = timecorr;
}
}
IN.close();
}
}
int caen_correction::init (Packet *p)
{
int chip,c,i,idx;
for ( chip = 0; chip < 4; chip++)
{
int cell = p->iValue(chip,"INDEXCELL");
//correct time for each chip
idx = cell;
for ( i = 0; i < 1024; i++)
{
current_time[i][chip] = timevec[idx][chip];
if (idx >=1024) idx=0;
}
// the adc samples
for ( c = 0; c < 8; c++)
{
idx = cell;
for ( i = 0; i < 1024; i++)
{
current_wave[i][chip*8+c] = p->iValue(i,chip*8+c) - base[idx][chip*9+c];
idx++;
if (idx >=1024) idx=0;
}
}
// // the TR cells
// idx = cell;
//for ( i = 0; i < 1024; i++)
// {
// if (chip < 2) current_wave[i][chip*8+8] = p->iValue(i,"TR0") - base[idx][chip*9+8];
// else current_wave[i][chip*9+8] = p->iValue(i,"TR1") - base[idx][chip*9+8];
// idx++;
// if (idx >=1024) idx=0;
// }
}
return 0;
}
float caen_correction::caen_corrected(const int sample, const int channel) const
{
if ( sample < 0 || sample >1023 || channel < 0 || channel > 31) return 0;
return current_wave[sample][channel];
}
float caen_correction::caen_time(const int sample, const int channel) const
{
if ( sample < 0 || sample >1023 || channel < 0 || channel > 31) return 0;
return current_wave[sample][channel/8];
}