|
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
|
/****************************************************************************
*
* This is a part of CMS-TOTEM PPS offline software.
* Authors:
* Jan Kašpar (jan.kaspar@gmail.com)
* Helena Malbouisson
* Clemencia Mora Herrera
*
****************************************************************************/
#include "FWCore/Utilities/interface/Exception.h"
#include "CondFormats/PPSObjects/interface/CTPPSRPAlignmentCorrectionData.h"
#include <Math/RotationZYX.h>
using namespace std;
//----------------------------------------------------------------------------------------------------
CTPPSRPAlignmentCorrectionData::CTPPSRPAlignmentCorrectionData(double _sh_x,
double _sh_x_u,
double _sh_y,
double _sh_y_u,
double _sh_z,
double _sh_z_u,
double _rot_x,
double _rot_x_u,
double _rot_y,
double _rot_y_u,
double _rot_z,
double _rot_z_u)
: sh_x(_sh_x),
sh_y(_sh_y),
sh_z(_sh_z),
sh_x_unc(_sh_x_u),
sh_y_unc(_sh_y_u),
sh_z_unc(_sh_z_u),
rot_x(_rot_x),
rot_y(_rot_y),
rot_z(_rot_z),
rot_x_unc(_rot_x_u),
rot_y_unc(_rot_y_u),
rot_z_unc(_rot_z_u) {}
//----------------------------------------------------------------------------------------------------
CTPPSRPAlignmentCorrectionData::CTPPSRPAlignmentCorrectionData(
double _sh_x, double _sh_y, double _sh_z, double _rot_x, double _rot_y, double _rot_z)
: sh_x(_sh_x), sh_y(_sh_y), sh_z(_sh_z), rot_x(_rot_x), rot_y(_rot_y), rot_z(_rot_z) {}
//----------------------------------------------------------------------------------------------------
void CTPPSRPAlignmentCorrectionData::add(const CTPPSRPAlignmentCorrectionData& a,
bool sumErrors,
bool addSh,
bool addRot) {
if (addSh) {
sh_x += a.sh_x;
sh_y += a.sh_y;
sh_z += a.sh_z;
if (sumErrors) {
sh_x_unc = sqrt(sh_x_unc * sh_x_unc + a.sh_x_unc * a.sh_x_unc);
sh_y_unc = sqrt(sh_y_unc * sh_y_unc + a.sh_y_unc * a.sh_y_unc);
sh_z_unc = sqrt(sh_z_unc * sh_z_unc + a.sh_z_unc * a.sh_z_unc);
}
}
if (addRot) {
rot_x += a.rot_x;
rot_y += a.rot_y;
rot_z += a.rot_z;
if (sumErrors) {
rot_x_unc = sqrt(rot_x_unc * rot_x_unc + a.rot_x_unc * a.rot_x_unc);
rot_y_unc = sqrt(rot_y_unc * rot_y_unc + a.rot_y_unc * a.rot_y_unc);
rot_z_unc = sqrt(rot_z_unc * rot_z_unc + a.rot_z_unc * a.rot_z_unc);
}
}
}
//----------------------------------------------------------------------------------------------------
std::ostream& operator<<(std::ostream& s, const CTPPSRPAlignmentCorrectionData& corr) {
s << fixed << "shift (um) " << std::setprecision(1) << "x = " << corr.getShX() * 1E3 << " +- "
<< corr.getShXUnc() * 1E3 << ", y = " << corr.getShY() * 1E3 << " +- " << corr.getShYUnc() * 1E3
<< ", z = " << corr.getShZ() * 1E3 << " +- " << corr.getShZUnc() * 1E3 << ", rotation (mrad) "
<< std::setprecision(1) << "x = " << corr.getRotX() * 1E3 << " +- " << corr.getRotXUnc() * 1E3
<< ", y = " << corr.getRotY() * 1E3 << " +- " << corr.getRotYUnc() * 1E3 << ", z = " << corr.getRotZ() * 1E3
<< " +- " << corr.getRotZUnc() * 1E3;
return s;
}
|