|
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
|
/*
* \class TPNCor
*
* \author: Julie Malcles - CEA/Saclay
*/
#include "CalibCalorimetry/EcalLaserAnalyzer/interface/TPNCor.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <sstream>
#include <iostream>
#include <iomanip>
#include <fstream>
#include <cmath>
using namespace std;
//ClassImp(TPNCor)
// Constructor...
TPNCor::TPNCor(string filename) {
// Initialize
isFileOK = 0;
for (int i = 0; i < iSizePar; i++) {
for (int j = 0; j < iSizeGain; j++) {
corParams[j][i] = 0.0;
}
}
// Get values from file
FILE *test;
test = fopen(filename.c_str(), "r");
int gain;
double aa, bb, cc;
std::ifstream fin;
if (test) {
fclose(test);
fin.open(filename.c_str());
while (fin.peek() != EOF) {
fin >> gain >> aa >> bb >> cc;
if (gain < iSizeGain) {
corParams[gain][0] = aa;
corParams[gain][1] = bb;
corParams[gain][2] = cc;
}
}
isFileOK = 1;
fin.close();
} else {
cout << " No PN linearity corrections file found, no correction will be applied " << endl;
}
}
// Destructor
TPNCor::~TPNCor() {}
double TPNCor::getPNCorrectionFactor(double val0, int gain) {
double cor = 0;
double corr = 1.0;
double pn = val0;
double xpn = val0 / 1000.0;
if (isFileOK == 0)
return 1.0;
if (gain > iSizeGain)
cout << "Unknown gain, gain has to be lower than " << iSizeGain << endl;
if (gain < iSizeGain) {
cor = xpn * (corParams[gain][0] + xpn * (corParams[gain][1] + xpn * corParams[gain][2]));
if (pn != 0)
corr = 1.0 - cor / pn;
else
corr = 1.0;
}
return corr;
}
|