-
Notifications
You must be signed in to change notification settings - Fork 1
/
HiRASiCalibration.cpp
134 lines (111 loc) · 4.73 KB
/
HiRASiCalibration.cpp
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
127
128
129
130
131
132
133
134
#include "include/HiRASiCalibration.h"
//____________________________________________________
HiRASiCalibration::HiRASiCalibration(int TotNumTel, int TotNumStripf, int TotNumStripb) :
fNumTel(TotNumTel),
fNumStripf(TotNumStripf),
fNumStripb(TotNumStripb),
fIntercept(new double[fNumTel*fNumStripf+fNumTel*fNumStripb]),
fSlope(new double[fNumTel*fNumStripf+fNumTel*fNumStripb]),
fMatchingIntercept(new double[fNumTel*fNumStripf+fNumTel*fNumStripb]),
fMatchingSlope(new double[fNumTel*fNumStripf+fNumTel*fNumStripb]),
fMatchingSaturationPoint(new double[fNumTel*fNumStripf+fNumTel*fNumStripb]),
fCalibrationLoaded(false),
fMatchingLoaded(false)
{}
//____________________________________________________
HiRASiCalibration::~HiRASiCalibration()
{
delete [] fIntercept;
delete [] fSlope;
delete [] fMatchingIntercept;
delete [] fMatchingSlope;
delete [] fMatchingSaturationPoint;
}
//____________________________________________________
int HiRASiCalibration::LoadCalibration(const char * file_name)
{
std::ifstream FileIn(file_name);
if(!FileIn.is_open()) {
return -1;
}
int NRead=0;
while (!FileIn.eof())
{
std::string LineRead;
std::getline(FileIn, LineRead);
if(LineRead.empty()) continue;
if(LineRead.find('*')==0) continue;
if(LineRead.find_first_not_of(' ') == std::string::npos) continue;
std::istringstream LineStream(LineRead);
int numtel;
int numstrip;
bool FB;
LineStream >> FB >> numtel >> numstrip;
LineStream >> fIntercept[FB*fNumTel*fNumStripf+numtel*(FB?fNumStripb:fNumStripf)+numstrip]
>> fSlope[FB*fNumTel*fNumStripf+numtel*(FB?fNumStripb:fNumStripf)+numstrip];
NRead++;
}
fCalibrationLoaded=true;
return NRead;
}
//____________________________________________________
int HiRASiCalibration::LoadMatching(const char * file_name)
{
std::ifstream FileIn(file_name);
if(!FileIn.is_open()) {
return -1;
}
int NRead=0;
while (!FileIn.eof())
{
std::string LineRead;
std::getline(FileIn, LineRead);
if(LineRead.empty()) continue;
if(LineRead.find('*')==0) continue;
if(LineRead.find_first_not_of(' ') == std::string::npos) continue;
std::istringstream LineStream(LineRead);
int numtel;
int numstrip;
bool FB;
LineStream >> numtel >> numstrip >> FB;
LineStream >> fMatchingIntercept[FB*fNumTel*fNumStripf+numtel*(FB?fNumStripb:fNumStripf)+numstrip]
>> fMatchingSlope[FB*fNumTel*fNumStripf+numtel*(FB?fNumStripb:fNumStripf)+numstrip]>>fMatchingSaturationPoint[FB*fNumTel*fNumStripf+numtel*(FB?fNumStripb:fNumStripf)+numstrip];
NRead++;
}
fMatchingLoaded=true;
return NRead;
}
//____________________________________________________
double HiRASiCalibration::GetEnergy(double ch, int num_tel, int num_strip, bool FrontBack) const
{
return fCalibrationLoaded ? ch*fSlope[FrontBack*fNumTel*fNumStripf+num_tel*(FrontBack?fNumStripb:fNumStripf)+num_strip]
+fIntercept[FrontBack*fNumTel*fNumStripf+num_tel*(FrontBack?fNumStripb:fNumStripf)+num_strip] : -9999;
}
//____________________________________________________
double HiRASiCalibration::GetSlope(int num_tel, int num_strip, bool FrontBack) const
{
return fCalibrationLoaded ? fSlope[FrontBack*fNumTel*fNumStripf+num_tel*(FrontBack?fNumStripb:fNumStripf)+num_strip] : -9999;
}
//____________________________________________________
double HiRASiCalibration::GetIntercept(int num_tel, int num_strip, bool FrontBack) const
{
return fCalibrationLoaded ? fIntercept[FrontBack*fNumTel*fNumStripf+num_tel*(FrontBack?fNumStripb:fNumStripf)+num_strip] : -9999;
}
//____________________________________________________
double HiRASiCalibration::GetHiLowMatched(double chHi, double chLow, int num_tel, int num_strip, bool FrontBack) const
{
if(!fMatchingLoaded) return -9999;
if(chLow<fMatchingSaturationPoint[FrontBack*fNumTel*fNumStripf+num_tel*(FrontBack?fNumStripb:fNumStripf)+num_strip]) {
return chHi;
}
else {
return chLow*fMatchingSlope[FrontBack*fNumTel*fNumStripf+num_tel*(FrontBack?fNumStripb:fNumStripf)+num_strip]
+fMatchingIntercept[FrontBack*fNumTel*fNumStripf+num_tel*(FrontBack?fNumStripb:fNumStripf)+num_strip];
}
}
//____________________________________________________
double HiRASiCalibration::GetEnergyHiLowMatched(double chHi, double chLow, int num_tel, int num_strip, bool FrontBack) const
{
return fMatchingLoaded ? GetHiLowMatched(chHi,chLow,num_tel,num_strip,FrontBack)*fSlope[FrontBack*fNumTel*fNumStripf+num_tel*(FrontBack?fNumStripb:fNumStripf)+num_strip]
+fIntercept[FrontBack*fNumTel*fNumStripf+num_tel*(FrontBack?fNumStripb:fNumStripf)+num_strip] : -9999;
}