|
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
135
136
137
138
139
140
141
142
143
144
145
146
|
// COCOA class implementation file
// Id: Measurement.C
// CAT: Model
// ---------------------------------------------------------------------------
// History: v1.0
// Authors:
// Pedro Arce
#include "Alignment/CocoaModel/interface/MeasurementCOPS.h"
#include "Alignment/CocoaModel/interface/LightRay.h"
#include "Alignment/CocoaModel/interface/OpticalObject.h"
#ifdef COCOA_VIS
#include "Alignment/CocoaVisMgr/interface/ALIVRMLMgr.h"
#include "Alignment/IgCocoaFileWriter/interface/IgCocoaFileMgr.h"
#include "Alignment/IgCocoaFileWriter/interface/ALIVisLightPath.h"
#endif
#include <iostream>
#include <iomanip>
#include <cstdlib>
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ calculate the simulated value propagating the light ray through the OptO that take part in the Measurement
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void MeasurementCOPS::calculateSimulatedValue(ALIbool firstTime) {
if (ALIUtils::debug >= 2)
printStartCalculateSimulatedValue(this); // important for Examples/FakeMeas
//---------- Create light ray
LightRay* lightray = new LightRay;
//---------- Define types of OptO that may take part in the Measurement
ALIuint isec = 0; //security variable to check OptOList().size()
//---------- Loop list of OptO that take part in measurement
std::vector<OpticalObject*>::const_iterator vocite = OptOList().begin();
if (ALIUtils::debug >= 5)
std::cout << "OptOList size" << OptOList().size() << std::endl;
//----- Check that first object is 'Xlaser'
if ((*vocite)->type() != "Xlaser") {
std::cerr << "!!ERROR MeasurementCOPS: first Optical object should be 'Xlaser'" << std::endl;
DumpBadOrderOptOs();
exit(1);
}
//---------- Check that last object is a COPS Sensor (that makes measuremnt and kill the lightray)
if ((*(OptOList().end() - 1))->type() != "COPS") {
std::cerr << "!!ERROR MeasurementCOPS: last Optical object should be 'COPS'" << std::endl;
DumpBadOrderOptOs();
exit(1);
}
#ifdef COCOA_VIS
ALIVisLightPath* vispath = 0;
if (ALIUtils::getFirstTime()) {
GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
if (gomgr->GlobalOptions()["VisWriteIguana"] > 1) {
vispath = IgCocoaFileMgr::getInstance().newLightPath(name());
}
}
#endif
std::cout << (vocite == OptOList().end()) << " vocite " << (*vocite)->name() << std::endl;
while (vocite != OptOList().end()) {
if (ALIUtils::debug >= -2)
std::cout << std::endl << "@@@@ LR:OBJECT " << (*vocite)->name() << std::endl;
isec++;
//---------- Get the behaviour of the object w.r.t the measurement (if it reflects the light, let it traverse it, ...)
ALIstring behav = getMeasuringBehaviour(vocite);
if (lightray) {
(*vocite)->participateInMeasurement(*lightray, *this, behav);
#ifdef COCOA_VIS
if (ALIUtils::getFirstTime()) {
GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
if (gomgr->GlobalOptions()["VisWriteVRML"] > 1) {
ALIVRMLMgr::getInstance().addLightPoint(lightray->point());
// std::cout << "ALIVRMLMg addLightPoint " << lightray->point() << (*vocite)->name() << std::endl;
}
if (gomgr->GlobalOptions()["VisWriteIguana"] > 1) {
vispath->addLightPoint(lightray->point(), *vocite);
}
}
#endif
} else {
std::cerr << "!! Last object is not Sensor 2D in measurement " << name() << std::endl;
DumpBadOrderOptOs();
exit(1);
}
++vocite;
if (isec > OptOList().size()) {
std::cerr << "ERROR DE PROGRAMACION EN GetSimulatedValue" << std::endl;
exit(5);
}
}
if (ALIUtils::debug >= 9)
std::cout << "end calculateSimulatedValue" << std::endl;
delete lightray;
}
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@ You input 8 numbers after 'TILMETER':
//@@
//@@ set the conversion factor from mV to mrad and the pedestal
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void MeasurementCOPS::setConversionFactor(const std::vector<ALIstring>& wordlist) {
//--------- Set it to 0
ALIuint ii;
for (ii = 0; ii < dim(); ii++) {
theDisplace[ii] = 0.;
}
//--------- Check that the format is OK
if (wordlist.size() == 1)
return;
if (wordlist.size() != 3 || !ALIUtils::IsNumber(wordlist[1]) || !ALIUtils::IsNumber(wordlist[2]) ||
!ALIUtils::IsNumber(wordlist[3]) || !ALIUtils::IsNumber(wordlist[4])) {
std::cerr << "!! SensorCOPS Measurement setConversionFactor: WRONG FORMAT " << std::endl
<< "It should be: SENSOR2D displace_U displace_D displace_L displace_R " << std::endl
<< "It is: ";
ALIUtils::dumpVS(wordlist, " ", std::cerr);
exit(1);
}
for (ii = 0; ii < dim(); ii++) {
theDisplace[ii] = atof(wordlist[ii + 1].c_str()) * valueDimensionFactor();
}
}
void MeasurementCOPS::correctValueAndSigma() {
//---------- Make displacement
ALIuint ii;
for (ii = 0; ii < dim(); ii++) {
ALIdouble val = value()[ii];
val += theDisplace[ii];
if (ALIUtils::debug >= 9)
std::cout << "MeasurementCOPS::correctValueAndSigma: old value X " << value()[ii] << " new " << val << " +- "
<< std::endl;
setValue(ii, val);
}
}
|