|
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
|
// COCOA class header file
//Id: ALIUtils.h
//CAT: Model
//
// Class with some utility function
//
// History: v1.0
// Pedro Arce
#ifndef CocoaUtils_HH
#define CocoaUtils_HH
#include "Alignment/CocoaUtilities/interface/CocoaGlobals.h"
#include "CLHEP/Vector/ThreeVector.h"
#include "CLHEP/Vector/Rotation.h"
#include <vector>
#include <ctime>
#include <fstream>
#include <iostream>
#include <cmath>
class ALIUtils {
public:
ALIUtils() {}
~ALIUtils() {}
static int IsNumber(const ALIstring& str);
static void dump3v(const CLHEP::Hep3Vector& vec, const std::string& msg);
static void dumprm(const CLHEP::HepRotation& rm, const std::string& msg, std::ostream& out = std::cout);
// public static DATA MEMBERS
static ALIint report;
static ALIint debug;
static ALIdouble deg;
static void setReportVerbosity(ALIint val) { report = val; }
static void setDebugVerbosity(ALIint val) { debug = val; }
static time_t time_now() { return _time_now; }
static void set_time_now(time_t now) { _time_now = now; }
//! Convert a string to an float, checking that it is really a number
static double getFloat(const ALIstring& str);
//! Convert a string to an integer, checking that it is really an integer
static int getInt(const ALIstring& str);
//! Convert a bool to an integer, checking that it is really a bool
static bool getBool(const ALIstring& str);
//! dumps a vector of strings with a message to outs
static void dumpVS(const std::vector<ALIstring>& wl, const std::string& msg, std::ostream& outs = std::cout);
//---------- Dimension factors
static void SetLengthDimensionFactors();
static void SetAngleDimensionFactors();
static void SetOutputLengthDimensionFactors();
static void SetOutputAngleDimensionFactors();
static ALIdouble CalculateLengthDimensionFactorFromInt(ALIint ad);
static ALIdouble CalculateAngleDimensionFactorFromInt(ALIint ad);
static ALIdouble CalculateLengthDimensionFactorFromString(ALIstring dimstr);
static ALIdouble CalculateAngleDimensionFactorFromString(ALIstring dimstr);
static void dumpDimensions(std::ofstream& fout);
static ALIdouble LengthValueDimensionFactor() { return _LengthValueDimensionFactor; }
static ALIdouble LengthSigmaDimensionFactor() { return _LengthSigmaDimensionFactor; }
static ALIdouble AngleValueDimensionFactor() { return _AngleValueDimensionFactor; }
static ALIdouble AngleSigmaDimensionFactor() { return _AngleSigmaDimensionFactor; }
static ALIdouble OutputLengthValueDimensionFactor() { return _OutputLengthValueDimensionFactor; }
static ALIdouble OutputLengthSigmaDimensionFactor() { return _OutputLengthSigmaDimensionFactor; }
static ALIdouble OutputAngleValueDimensionFactor() { return _OutputAngleValueDimensionFactor; }
static ALIdouble OutputAngleSigmaDimensionFactor() { return _OutputAngleSigmaDimensionFactor; }
static ALIdouble val0(ALIdouble val) {
//-std::cout << val << " val " << ( (val <= 1.E-9) ? 0. : val) << std::endl;
// return (abs(val) <= 1.E-9) ? 0. : val; }
if (std::fabs(val) <= 1.E-9) {
return 0.;
} else {
return val;
};
}
static ALIstring subQuotes(const ALIstring& str);
static ALIdouble getDimensionValue(const ALIstring& dim, const ALIstring& dimType);
static std::string changeName(const std::string& oldName, const std::string& subsstr1, const std::string& subsstr2);
static ALIbool getFirstTime() { return firstTime; }
static void setFirstTime(ALIbool val) { firstTime = val; }
static ALIdouble getMaximumDeviationDerivative() { return maximum_deviation_derivative; }
static void setMaximumDeviationDerivative(ALIdouble val) { maximum_deviation_derivative = val; }
static std::vector<double> getRotationAnglesFromMatrix(const CLHEP::HepRotation& rmLocal,
double origAngleX,
double origAngleY,
double origAngleZ);
static double diff2pi(double ang1, double ang2);
static bool eq2ang(double ang1, double ang2);
static double approxTo0(double val);
static double addPii(double val);
static int checkMatrixEquations(double angleX, double angleY, double angleZ, const CLHEP::HepRotation& rot);
private:
static ALIdouble _LengthValueDimensionFactor;
static ALIdouble _LengthSigmaDimensionFactor;
static ALIdouble _AngleValueDimensionFactor;
static ALIdouble _AngleSigmaDimensionFactor;
static ALIdouble _OutputLengthValueDimensionFactor;
static ALIdouble _OutputLengthSigmaDimensionFactor;
static ALIdouble _OutputAngleValueDimensionFactor;
static ALIdouble _OutputAngleSigmaDimensionFactor;
static time_t _time_now;
static ALIbool firstTime;
static ALIdouble maximum_deviation_derivative;
};
/*
template<class T>
ALIuint FindItemInVector( const T* item, const std::vector<T*>& item_vector )
{
std::vector<T*>::const_iterator vtcite;
ALIuint vfound = 1;
for( vtcite = item_vector.begin(); vtcite != item_vector.end(); vtcite++) {
if( (*vtcite) == item ) {
}
}
}
*/
//std::ostream& operator << (std::ostream& os, const CLHEP::HepRotation& c);
#endif
|