|
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
|
// COCOA class implementation file
//Id: FittedEntry.cc
//CAT: Model
//
// History: v1.0
// Pedro Arce
#include "Alignment/CocoaFit/interface/FittedEntry.h"
#include "Alignment/CocoaModel/interface/OpticalObject.h"
#include "Alignment/CocoaModel/interface/Entry.h"
#include "Alignment/CocoaModel/interface/Model.h"
#include "Alignment/CocoaUtilities/interface/ALIUtils.h"
#include <iostream>
#include <iomanip>
#include <cstdlib>
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@2
FittedEntry::FittedEntry(Entry* entry, ALIint order, ALIdouble sigma) {
theEntry = entry;
theOrder = order;
theOptOName = entry->OptOCurrent()->longName();
if (ALIUtils::debug >= 5)
std::cout << " creating FittedEntry " << theOptOName << std::endl;
theEntryName = entry->name();
BuildName();
//------ store values and sigmas in dimensions indicated by global options
ALIdouble dimv = entry->OutputValueDimensionFactor();
ALIdouble dims = entry->OutputSigmaDimensionFactor();
theValue = (entry->value() + entry->valueDisplacementByFitting()) / dimv;
theSigma = sigma / dims;
theOrigValue = entry->value() / dimv;
theOrigSigma = entry->sigma() / dims;
theQuality = entry->quality();
//-std::cout << this << " FE value" << this->theValue << "sigma" << this->theSigma << std::endl;
}
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@2
FittedEntry::FittedEntry(ALIstring name, float value, float sigma) {
//ar.lass1.laser.centre_X
theOrder = 0;
theOptOName = "s";
ALIint point = -1;
ALIint pointold = 0;
for (;;) {
point = name.find('.', point + 1);
if (point == -1)
break;
theOptOName += "/" + name.substr(pointold, point - pointold);
pointold = point + 1;
}
theEntryName = name.substr(pointold, name.size());
// std::cout << " building theEntryName " << theEntryName << " " << pointold << " " << name << std::endl;
Entry* entry = Model::getEntryByName(theOptOName, theEntryName);
theEntry = nullptr;
//------ store values and sigmas in dimensions indicated by global options
ALIdouble dimv = entry->OutputValueDimensionFactor();
ALIdouble dims = entry->OutputSigmaDimensionFactor();
theValue = value * dimv;
theSigma = sigma * dims;
theOrigValue = value * dimv;
theOrigSigma = sigma * dims;
theQuality = 2;
//-std::cout << this << " FE value" << this->theValue << "sigma" << this->theSigma << std::endl;
}
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@2
FittedEntry::FittedEntry(const std::vector<FittedEntry*>& _vFEntry) {
//----- Average the entries
std::vector<FittedEntry*> vFEntry = _vFEntry;
std::vector<FittedEntry*>::iterator ite;
//--- First check that all entries are from the same OptO and Entry
theOptOName = (vFEntry[0]->getOptOName());
theEntryName = (vFEntry[0]->getEntryName());
theOrder = (vFEntry[0]->getOrder());
theEntry = (vFEntry[0]->getEntry());
theQuality = (vFEntry[0]->getQuality());
theValue = 0.;
theSigma = 0.;
theOrigValue = 0.;
theOrigSigma = 0.;
for (ite = vFEntry.begin(); ite != vFEntry.end(); ++ite) {
if ((*ite)->getOptOName() != theOptOName || (*ite)->getEntryName() != theEntryName) {
std::cerr << "!!! FATAL ERROR FittedEntry::FittedEntry one entry in list has different opto or entry names : "
<< (*ite)->getOptOName() << " != " << theOptOName << " " << (*ite)->getEntryName()
<< " != " << theEntryName << std::endl;
exit(1);
}
theValue += (*ite)->getValue();
theSigma += (*ite)->getSigma();
theOrigValue += (*ite)->getOrigValue();
theOrigSigma += (*ite)->getOrigSigma();
}
ALIint siz = vFEntry.size();
theValue /= siz;
theSigma /= siz;
theOrigValue /= siz;
theOrigSigma /= siz;
}
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@2
void FittedEntry::BuildName() {
//----- substitute '/' by '.' in opto name
theName = theOptOName.substr(2, theOptOName.size()); // skip the first 's/'
ALIint slash = -1;
for (;;) {
slash = theName.find('/', slash + 1);
if (slash == -1)
break;
theName[slash] = '.';
}
//----- Check if there is a ' ' in entry (should not happen now)
ALIint space = theEntryName.rfind(' ');
theName.append(".");
ALIstring en = theEntryName;
if (space != -1)
en[space] = '_';
//----- Merge opto and entry names
// now it is not used as theName theName.append( en + ".out");
theName.append(en);
}
|