19 #include "xercesc/parsers/XercesDOMParser.hpp"
20 #include "xercesc/dom/DOM.hpp"
21 #include "xercesc/sax/HandlerBase.hpp"
22 #include "xercesc/util/XMLString.hpp"
23 #include "xercesc/util/PlatformUtils.hpp"
24 #include "xercesc/util/XercesDefs.hpp"
25 XERCES_CPP_NAMESPACE_USE
48 : m_fileName(fileName)
52 str_name = XMLString::transcode(
"name");
73 str_rawId = XMLString::transcode(
"rawId");
74 str_wheel = XMLString::transcode(
"wheel");
78 str_layer = XMLString::transcode(
"layer");
80 str_ring = XMLString::transcode(
"ring");
82 str_axisx = XMLString::transcode(
"axisx");
83 str_axisy = XMLString::transcode(
"axisy");
84 str_axisz = XMLString::transcode(
"axisz");
85 str_angle = XMLString::transcode(
"angle");
86 str_x = XMLString::transcode(
"x");
87 str_y = XMLString::transcode(
"y");
88 str_z = XMLString::transcode(
"z");
89 str_phix = XMLString::transcode(
"phix");
90 str_phiy = XMLString::transcode(
"phiy");
91 str_phiz = XMLString::transcode(
"phiz");
92 str_alpha = XMLString::transcode(
"alpha");
93 str_beta = XMLString::transcode(
"beta");
94 str_gamma = XMLString::transcode(
"gamma");
95 str_rphi = XMLString::transcode(
"rphi");
96 str_phi = XMLString::transcode(
"phi");
97 str_xx = XMLString::transcode(
"xx");
98 str_xy = XMLString::transcode(
"xy");
99 str_xz = XMLString::transcode(
"xz");
100 str_xa = XMLString::transcode(
"xa");
101 str_xb = XMLString::transcode(
"xb");
102 str_xc = XMLString::transcode(
"xc");
103 str_yy = XMLString::transcode(
"yy");
104 str_yz = XMLString::transcode(
"yz");
105 str_ya = XMLString::transcode(
"ya");
106 str_yb = XMLString::transcode(
"yb");
107 str_yc = XMLString::transcode(
"yc");
108 str_zz = XMLString::transcode(
"zz");
109 str_za = XMLString::transcode(
"za");
110 str_zb = XMLString::transcode(
"zb");
111 str_zc = XMLString::transcode(
"zc");
112 str_aa = XMLString::transcode(
"aa");
113 str_ab = XMLString::transcode(
"ab");
114 str_ac = XMLString::transcode(
"ac");
115 str_bb = XMLString::transcode(
"bb");
116 str_bc = XMLString::transcode(
"bc");
117 str_cc = XMLString::transcode(
"cc");
118 str_none = XMLString::transcode(
"none");
119 str_ideal = XMLString::transcode(
"ideal");
219 for (align::Alignables::const_iterator ali = alignables.begin(); ali != alignables.end(); ++ali) {
223 alignableNavigator[(*ali)->geomDetId().rawId()] = *ali;
230 align::Alignables::const_iterator alignable = alignables.begin();
231 align::Alignables::const_iterator ideal = ideals.begin();
233 while (alignable != alignables.end() && ideal != ideals.end()) {
234 alitoideal[*alignable] = *ideal;
236 fillAliToIdeal(alitoideal, (*alignable)->components(), (*ideal)->components());
242 if (alignable != alignables.end() || ideal != ideals.end()) {
243 throw cms::Exception(
"Alignment") <<
"alignable and ideal-alignable trees are out of sync (this should never happen)";
252 std::map<unsigned int, Alignable*> alignableNavigator;
257 std::map<unsigned int, Alignable*> ideal_alignableNavigator;
258 recursiveGetId(ideal_alignableNavigator, ideal_alignableMuon->DTBarrel());
259 recursiveGetId(ideal_alignableNavigator, ideal_alignableMuon->CSCEndcaps());
264 catch (
const XMLException &toCatch) {
265 throw cms::Exception(
"XMLException") <<
"Xerces XML parser threw an exception on initialization." << std::endl;
268 XercesDOMParser *
parser =
new XercesDOMParser();
269 parser->setValidationScheme(XercesDOMParser::Val_Always);
271 XERCES_CPP_NAMESPACE::ErrorHandler *errHandler = (XERCES_CPP_NAMESPACE::ErrorHandler*)(
new HandlerBase());
272 parser->setErrorHandler(errHandler);
277 catch (
const XMLException &toCatch) {
278 char *
message = XMLString::transcode(toCatch.getMessage());
279 throw cms::Exception(
"XMLException") <<
"Xerces XML parser threw this exception: " << message << std::endl;
281 catch (
const DOMException &toCatch) {
282 char *
message = XMLString::transcode(toCatch.msg);
283 throw cms::Exception(
"XMLException") <<
"Xerces XML parser threw this exception: " << message << std::endl;
285 catch (
const SAXException &toCatch) {
286 char *
message = XMLString::transcode(toCatch.getMessage());
287 throw cms::Exception(
"XMLException") <<
"Xerces XML parser threw this exception: " << message << std::endl;
290 throw cms::Exception(
"XMLException") <<
"Xerces XML parser threw an unknown exception" << std::endl;
293 DOMDocument *doc = parser->getDocument();
294 DOMElement *node_MuonAlignment = doc->getDocumentElement();
295 DOMNodeList *collections = doc->getElementsByTagName(
str_collection);
296 DOMNodeList *operations = doc->getElementsByTagName(
str_operation);
298 std::map<Alignable*, Alignable*> alitoideal;
302 std::map<std::string, std::map<Alignable*, bool> > alicollections;
303 for (
unsigned int i = 0;
i < collections->getLength();
i++) {
304 DOMElement *
collection = (DOMElement*)(collections->item(
i));
305 if (collection->getParentNode() == node_MuonAlignment) {
306 DOMNodeList *children = collection->getChildNodes();
308 DOMAttr *node_name = collection->getAttributeNode(
str_name);
309 if (node_name ==
NULL) {
310 throw cms::Exception(
"XMLException") <<
"<collection> requires a name attribute" << std::endl;
312 char *ascii_name = XMLString::transcode(node_name->getValue());
316 std::map<Alignable*, bool> aliset;
317 for (
unsigned int j = 0;
j < children->getLength();
j++) {
318 DOMNode *node = children->item(
j);
320 if (node->getNodeType() == DOMNode::ELEMENT_NODE) {
323 throw cms::Exception(
"XMLException") <<
"<collection> must contain only alignables" << std::endl;
330 alicollections[
name] = aliset;
334 for (
unsigned int i = 0;
i < operations->getLength();
i++) {
335 DOMElement *operation = (DOMElement*)(operations->item(
i));
336 if (operation->getParentNode() != node_MuonAlignment) {
337 throw cms::Exception(
"XMLException") <<
"All operations must be top-level elements" << std::endl;
340 DOMNodeList *children = operation->getChildNodes();
342 std::map<Alignable*, bool> aliset;
343 std::vector<DOMNode*> nodesToRemove;
344 for (
unsigned int j = 0;
j < children->getLength();
j++) {
345 DOMNode *node = children->item(
j);
347 if (node->getNodeType() == DOMNode::ELEMENT_NODE) {
351 nodesToRemove.push_back(node);
354 else if (XMLString::equals(node->getNodeName(),
str_collection)) {
355 DOMAttr *node_name = ((DOMElement*)(node))->getAttributeNode(
str_name);
356 if (node_name ==
NULL) {
357 throw cms::Exception(
"XMLException") <<
"<collection> requires a name attribute" << std::endl;
359 char *ascii_name = XMLString::transcode(node_name->getValue());
363 std::map<std::string, std::map<Alignable*, bool> >::const_iterator alicollections_iter = alicollections.find(name);
364 if (alicollections_iter == alicollections.end()) {
365 throw cms::Exception(
"XMLException") <<
"<collection name=\"" << name <<
"\"> hasn't been defined" << std::endl;
368 for (std::map<Alignable*, bool>::const_iterator aliiter = alicollections_iter->second.begin();
369 aliiter != alicollections_iter->second.end();
371 aliset[aliiter->first] =
true;
374 nodesToRemove.push_back(node);
383 for (std::vector<DOMNode*>::const_iterator node = nodesToRemove.begin(); node != nodesToRemove.end(); ++node) {
384 operation->removeChild(*node);
386 children = operation->getChildNodes();
388 for (
unsigned int j = 0;
j < children->getLength();
j++) {
389 DOMNode *node = children->item(
j);
390 if (node->getNodeType() == DOMNode::ELEMENT_NODE) {
396 else if (XMLString::equals(node->getNodeName(),
str_setape)) {
397 do_setape((DOMElement*)(node), aliset, alitoideal);
404 else if (XMLString::equals(node->getNodeName(),
str_moveglobal)) {
408 else if (XMLString::equals(node->getNodeName(),
str_movelocal)) {
425 char *
message = XMLString::transcode(node->getNodeName());
426 throw cms::Exception(
"XMLException") <<
"Unrecognized operation: \"" << message <<
"\"" << std::endl;
438 delete ideal_alignableMuon;
439 return alignableMuon;
458 unsigned int rawId = 0;
460 DOMAttr *node_rawId = node->getAttributeNode(
str_rawId);
461 if (node_rawId !=
NULL) {
463 rawId = XMLString::parseInt(node_rawId->getValue());
465 catch (
const XMLException &toCatch) {
466 throw cms::Exception(
"XMLException") <<
"Value of \"rawId\" must be an integer" << std::endl;
470 int wheel,
station, sector, superlayer, layer;
471 wheel = station = sector = superlayer = layer = 1;
474 DOMAttr *node_wheel = node->getAttributeNode(
str_wheel);
475 if (node_wheel ==
NULL)
throw cms::Exception(
"XMLException") <<
"DT node is missing required \"wheel\" attribute" << std::endl;
477 wheel = XMLString::parseInt(node_wheel->getValue());
479 catch (
const XMLException &toCatch) {
480 throw cms::Exception(
"XMLException") <<
"Value of \"wheel\" must be an integer" << std::endl;
484 DOMAttr *node_station = node->getAttributeNode(
str_station);
485 if (node_station ==
NULL)
throw cms::Exception(
"XMLException") <<
"DT node is missing required \"station\" attribute" << std::endl;
487 station = XMLString::parseInt(node_station->getValue());
489 catch (
const XMLException &toCatch) {
490 throw cms::Exception(
"XMLException") <<
"Value of \"station\" must be an integer" << std::endl;
494 DOMAttr *node_sector = node->getAttributeNode(
str_sector);
495 if (node_sector ==
NULL)
throw cms::Exception(
"XMLException") <<
"DT node is missing required \"sector\" attribute" << std::endl;
497 sector = XMLString::parseInt(node_sector->getValue());
499 catch (
const XMLException &toCatch) {
500 throw cms::Exception(
"XMLException") <<
"Value of \"sector\" must be an integer" << std::endl;
504 DOMAttr *node_superlayer = node->getAttributeNode(
str_superlayer);
505 if (node_superlayer ==
NULL)
throw cms::Exception(
"XMLException") <<
"DT node is missing required \"superlayer\" attribute" << std::endl;
507 superlayer = XMLString::parseInt(node_superlayer->getValue());
509 catch (
const XMLException &toCatch) {
510 throw cms::Exception(
"XMLException") <<
"Value of \"superlayer\" must be an integer" << std::endl;
514 DOMAttr *node_layer = node->getAttributeNode(
str_layer);
515 if (node_layer ==
NULL)
throw cms::Exception(
"XMLException") <<
"DT node is missing required \"layer\" attribute" << std::endl;
517 layer = XMLString::parseInt(node_layer->getValue());
519 catch (
const XMLException &toCatch) {
520 throw cms::Exception(
"XMLException") <<
"Value of \"layer\" must be an integer" << std::endl;
529 DTLayerId layerId(wheel, station, sector, superlayer, layer);
530 rawId = layerId.
rawId();
533 Alignable *ali = alignableNavigator[rawId];
534 if (ali ==
NULL)
throw cms::Exception(
"XMLException") <<
"rawId \"" << rawId <<
"\" is not recognized" << std::endl;
541 <<
"rawId \"" << rawId <<
"\" is not a "
551 DOMAttr *node_rawId = node->getAttributeNode(
str_rawId);
552 if (node_rawId !=
NULL) {
554 rawId = XMLString::parseInt(node_rawId->getValue());
556 catch (
const XMLException &toCatch) {
557 throw cms::Exception(
"XMLException") <<
"Value of \"rawId\" must be an integer" << std::endl;
562 endcap = station = ring = chamber = layer = 1;
564 DOMAttr *node_endcap = node->getAttributeNode(
str_endcap);
565 if (node_endcap ==
NULL)
throw cms::Exception(
"XMLException") <<
"CSC node is missing required \"endcap\" attribute" << std::endl;
567 endcap = XMLString::parseInt(node_endcap->getValue());
569 catch (
const XMLException &toCatch) {
570 throw cms::Exception(
"XMLException") <<
"Value of \"endcap\" must be an integer" << std::endl;
572 if (endcap == -1) endcap = 2;
575 DOMAttr *node_station = node->getAttributeNode(
str_station);
576 if (node_station ==
NULL)
throw cms::Exception(
"XMLException") <<
"CSC node is missing required \"station\" attribute" << std::endl;
578 station = XMLString::parseInt(node_station->getValue());
580 catch (
const XMLException &toCatch) {
581 throw cms::Exception(
"XMLException") <<
"Value of \"station\" must be an integer" << std::endl;
585 DOMAttr *node_ring = node->getAttributeNode(
str_ring);
586 if (node_ring ==
NULL)
throw cms::Exception(
"XMLException") <<
"CSC node is missing required \"ring\" attribute" << std::endl;
588 ring = XMLString::parseInt(node_ring->getValue());
590 catch (
const XMLException &toCatch) {
591 throw cms::Exception(
"XMLException") <<
"Value of \"ring\" must be an integer" << std::endl;
595 DOMAttr *node_chamber = node->getAttributeNode(
str_chamber);
596 if (node_chamber ==
NULL)
throw cms::Exception(
"XMLException") <<
"CSC node is missing required \"chamber\" attribute" << std::endl;
598 chamber = XMLString::parseInt(node_chamber->getValue());
600 catch (
const XMLException &toCatch) {
601 throw cms::Exception(
"XMLException") <<
"Value of \"chamber\" must be an integer" << std::endl;
605 DOMAttr *node_layer = node->getAttributeNode(
str_layer);
606 if (node_layer ==
NULL)
throw cms::Exception(
"XMLException") <<
"CSC node is missing required \"layer\" attribute" << std::endl;
608 layer = XMLString::parseInt(node_layer->getValue());
610 catch (
const XMLException &toCatch) {
611 throw cms::Exception(
"XMLException") <<
"Value of \"layer\" must be an integer" << std::endl;
619 CSCDetId layerId(endcap, station, ring, chamber, layer);
620 rawId = layerId.
rawId();
623 Alignable *ali = alignableNavigator[rawId];
624 if (ali ==
NULL)
throw cms::Exception(
"XMLException") <<
"rawId \"" << rawId <<
"\" is not recognized" << std::endl;
631 <<
"rawId \"" << rawId <<
"\" is not a "
639 unsigned int len = XMLString::stringLen(str);
640 char *cstr = XMLString::transcode(str);
641 std::stringstream errmessage;
642 errmessage <<
"Value of \"" << attribute <<
"\" must be a double, not \"" << cstr <<
"\"" << std::endl;
647 if (cstr[i] ==
'-') {
651 else if (cstr[i] ==
'+') i++;
655 while (cstr[i] !=
'.' && cstr[i] !=
'e' && cstr[i] !=
'E' && i < len) {
656 if (cstr[i] <
'0' || cstr[i] >
'9') {
662 output += cstr[
i] -
'0';
666 if (cstr[i] ==
'.') {
670 while (cstr[i] !=
'e' && cstr[i] !=
'E' && i < len) {
671 if (cstr[i] <
'0' || cstr[i] >
'9') {
676 output += (cstr[
i] -
'0') * place;
682 if (cstr[i] ==
'e' || cstr[i] ==
'E') {
686 bool expminus =
false;
687 if (cstr[i] ==
'-') {
691 else if (cstr[i] ==
'+') i++;
694 if (cstr[i] <
'0' || cstr[i] >
'9') {
700 exponent += cstr[
i] -
'0';
704 if (expminus) exponent *= -1;
706 output *=
pow(10., exponent);
709 if (minus) output *= -1.;
716 DOMAttr *node_relativeto = node->getAttributeNode(
str_relativeto);
717 if (node_relativeto ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setposition> is missing required \"relativeto\" attribute" << std::endl;
720 if (XMLString::equals(node_relativeto->getValue(),
str_none)) {
723 else if (XMLString::equals(node_relativeto->getValue(),
str_ideal)) {
726 else if (XMLString::equals(node_relativeto->getValue(),
str_container)) {
730 char *
message = XMLString::transcode(node_relativeto->getValue());
731 throw cms::Exception(
"XMLException") <<
"relativeto must be \"none\", \"ideal\", or \"container\", not \"" << message <<
"\"" << std::endl;
734 for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
743 if (relativeto == 0) {}
745 else if (relativeto == 1) {
753 else if (relativeto == 2 && ali->
mother() !=
NULL) {
758 oldrot = oldrot * globalRotation.
transposed();
761 double x = oldpos.
x();
762 double y = oldpos.
y();
763 double z = oldpos.
z();
765 double phix = atan2(oldrot.
yz(), oldrot.
zz());
766 double phiy = asin(-oldrot.
xz());
767 double phiz = atan2(oldrot.
xy(), oldrot.
xx());
770 double alpha = oldEulerAngles(1);
771 double beta = oldEulerAngles(2);
772 double gamma = oldEulerAngles(3);
776 DOMAttr *node_x = node->getAttributeNode(
str_x);
777 DOMAttr *node_y = node->getAttributeNode(
str_y);
778 DOMAttr *node_z = node->getAttributeNode(
str_z);
785 DOMAttr *node_phix = node->getAttributeNode(
str_phix);
786 DOMAttr *node_phiy = node->getAttributeNode(
str_phiy);
787 DOMAttr *node_phiz = node->getAttributeNode(
str_phiz);
788 DOMAttr *node_alpha = node->getAttributeNode(
str_alpha);
789 DOMAttr *node_beta = node->getAttributeNode(
str_beta);
790 DOMAttr *node_gamma = node->getAttributeNode(
str_gamma);
793 bool phixyz = (node_phix !=
NULL || node_phiy !=
NULL || node_phiz !=
NULL);
794 bool alphabetagamma = (node_alpha !=
NULL || node_beta !=
NULL || node_gamma !=
NULL);
795 if (phixyz && alphabetagamma)
throw cms::Exception(
"XMLException") <<
"<setposition> must either have phix, phiy, and phiz or alpha, beta, and gamma, but not both" << std::endl;
796 if (!phixyz && !alphabetagamma) alphabetagamma =
true;
799 if (node_phix !=
NULL) phix =
parseDouble(node_phix->getValue(),
"phix");
800 if (node_phiy !=
NULL) phiy =
parseDouble(node_phiy->getValue(),
"phiy");
801 if (node_phiz !=
NULL) phiz =
parseDouble(node_phiz->getValue(),
"phiz");
807 0., -
sin(phix),
cos(phix));
810 sin(phiy), 0.,
cos(phiy));
812 -
sin(phiz),
cos(phiz), 0.,
815 rot = rotX * rotY * rotZ;
818 else if (alphabetagamma) {
819 if (node_alpha !=
NULL) alpha =
parseDouble(node_alpha->getValue(),
"alpha");
820 if (node_beta !=
NULL) beta =
parseDouble(node_beta->getValue(),
"beta");
821 if (node_gamma !=
NULL) gamma =
parseDouble(node_gamma->getValue(),
"gamma");
825 eulerAngles(1) =
alpha;
826 eulerAngles(2) =
beta;
827 eulerAngles(3) = gamma;
833 if (relativeto == 0) {
837 else if (relativeto == 1) {
849 else if (relativeto == 2) {
853 if (container !=
NULL) {
902 if (survey !=
NULL) {
903 matrix6x6 = survey->
errors();
908 void MuonAlignmentInputXML::do_setape(
const XERCES_CPP_NAMESPACE::DOMElement *node, std::map<Alignable*, bool> &aliset, std::map<Alignable*, Alignable*> &alitoideal)
const {
909 DOMAttr *node_xx = node->getAttributeNode(
str_xx);
910 DOMAttr *node_xy = node->getAttributeNode(
str_xy);
911 DOMAttr *node_xz = node->getAttributeNode(
str_xz);
912 DOMAttr *node_xa = node->getAttributeNode(
str_xa);
913 DOMAttr *node_xb = node->getAttributeNode(
str_xb);
914 DOMAttr *node_xc = node->getAttributeNode(
str_xc);
915 DOMAttr *node_yy = node->getAttributeNode(
str_yy);
916 DOMAttr *node_yz = node->getAttributeNode(
str_yz);
917 DOMAttr *node_ya = node->getAttributeNode(
str_ya);
918 DOMAttr *node_yb = node->getAttributeNode(
str_yb);
919 DOMAttr *node_yc = node->getAttributeNode(
str_yc);
920 DOMAttr *node_zz = node->getAttributeNode(
str_zz);
921 DOMAttr *node_za = node->getAttributeNode(
str_za);
922 DOMAttr *node_zb = node->getAttributeNode(
str_zb);
923 DOMAttr *node_zc = node->getAttributeNode(
str_zc);
924 DOMAttr *node_aa = node->getAttributeNode(
str_aa);
925 DOMAttr *node_ab = node->getAttributeNode(
str_ab);
926 DOMAttr *node_ac = node->getAttributeNode(
str_ac);
927 DOMAttr *node_bb = node->getAttributeNode(
str_bb);
928 DOMAttr *node_bc = node->getAttributeNode(
str_bc);
929 DOMAttr *node_cc = node->getAttributeNode(
str_cc);
931 if (node_xx ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"xx\" attribute" << std::endl;
932 if (node_xy ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"xy\" attribute" << std::endl;
933 if (node_xz ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"xz\" attribute" << std::endl;
934 if (node_xa ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"xa\" attribute" << std::endl;
935 if (node_xb ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"xb\" attribute" << std::endl;
936 if (node_xc ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"xc\" attribute" << std::endl;
937 if (node_yy ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"yy\" attribute" << std::endl;
938 if (node_yz ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"yz\" attribute" << std::endl;
939 if (node_ya ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"ya\" attribute" << std::endl;
940 if (node_yb ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"yb\" attribute" << std::endl;
941 if (node_yc ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"yc\" attribute" << std::endl;
942 if (node_zz ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"zz\" attribute" << std::endl;
943 if (node_za ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"za\" attribute" << std::endl;
944 if (node_zb ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"zb\" attribute" << std::endl;
945 if (node_zc ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"zc\" attribute" << std::endl;
946 if (node_aa ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"aa\" attribute" << std::endl;
947 if (node_ab ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"ab\" attribute" << std::endl;
948 if (node_ac ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"ac\" attribute" << std::endl;
949 if (node_bb ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"bb\" attribute" << std::endl;
950 if (node_bc ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"bc\" attribute" << std::endl;
951 if (node_cc ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setape> is missing required \"cc\" attribute" << std::endl;
954 matrix6x6(0,0) =
parseDouble(node_xx->getValue(),
"xx");
955 matrix6x6(0,1) =
parseDouble(node_xy->getValue(),
"xy");
956 matrix6x6(0,2) =
parseDouble(node_xz->getValue(),
"xz");
957 matrix6x6(0,3) =
parseDouble(node_xa->getValue(),
"xa");
958 matrix6x6(0,4) =
parseDouble(node_xb->getValue(),
"xb");
959 matrix6x6(0,5) =
parseDouble(node_xc->getValue(),
"xc");
960 matrix6x6(1,1) =
parseDouble(node_yy->getValue(),
"yy");
961 matrix6x6(1,2) =
parseDouble(node_yz->getValue(),
"yz");
962 matrix6x6(1,3) =
parseDouble(node_ya->getValue(),
"ya");
963 matrix6x6(1,4) =
parseDouble(node_yb->getValue(),
"yb");
964 matrix6x6(1,5) =
parseDouble(node_yc->getValue(),
"yc");
965 matrix6x6(2,2) =
parseDouble(node_zz->getValue(),
"zz");
966 matrix6x6(2,3) =
parseDouble(node_za->getValue(),
"za");
967 matrix6x6(2,4) =
parseDouble(node_zb->getValue(),
"zb");
968 matrix6x6(2,5) =
parseDouble(node_zc->getValue(),
"zc");
969 matrix6x6(3,3) =
parseDouble(node_aa->getValue(),
"aa");
970 matrix6x6(3,4) =
parseDouble(node_ab->getValue(),
"ab");
971 matrix6x6(3,5) =
parseDouble(node_ac->getValue(),
"ac");
972 matrix6x6(4,4) =
parseDouble(node_bb->getValue(),
"bb");
973 matrix6x6(4,5) =
parseDouble(node_bc->getValue(),
"bc");
974 matrix6x6(5,5) =
parseDouble(node_cc->getValue(),
"cc");
976 for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
983 DOMAttr *node_xx = node->getAttributeNode(
str_xx);
984 DOMAttr *node_xy = node->getAttributeNode(
str_xy);
985 DOMAttr *node_xz = node->getAttributeNode(
str_xz);
986 DOMAttr *node_xa = node->getAttributeNode(
str_xa);
987 DOMAttr *node_xb = node->getAttributeNode(
str_xb);
988 DOMAttr *node_xc = node->getAttributeNode(
str_xc);
989 DOMAttr *node_yy = node->getAttributeNode(
str_yy);
990 DOMAttr *node_yz = node->getAttributeNode(
str_yz);
991 DOMAttr *node_ya = node->getAttributeNode(
str_ya);
992 DOMAttr *node_yb = node->getAttributeNode(
str_yb);
993 DOMAttr *node_yc = node->getAttributeNode(
str_yc);
994 DOMAttr *node_zz = node->getAttributeNode(
str_zz);
995 DOMAttr *node_za = node->getAttributeNode(
str_za);
996 DOMAttr *node_zb = node->getAttributeNode(
str_zb);
997 DOMAttr *node_zc = node->getAttributeNode(
str_zc);
998 DOMAttr *node_aa = node->getAttributeNode(
str_aa);
999 DOMAttr *node_ab = node->getAttributeNode(
str_ab);
1000 DOMAttr *node_ac = node->getAttributeNode(
str_ac);
1001 DOMAttr *node_bb = node->getAttributeNode(
str_bb);
1002 DOMAttr *node_bc = node->getAttributeNode(
str_bc);
1003 DOMAttr *node_cc = node->getAttributeNode(
str_cc);
1005 if (node_xx ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"xx\" attribute" << std::endl;
1006 if (node_xy ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"xy\" attribute" << std::endl;
1007 if (node_xz ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"xz\" attribute" << std::endl;
1008 if (node_xa ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"xa\" attribute" << std::endl;
1009 if (node_xb ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"xb\" attribute" << std::endl;
1010 if (node_xc ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"xc\" attribute" << std::endl;
1011 if (node_yy ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"yy\" attribute" << std::endl;
1012 if (node_yz ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"yz\" attribute" << std::endl;
1013 if (node_ya ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"ya\" attribute" << std::endl;
1014 if (node_yb ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"yb\" attribute" << std::endl;
1015 if (node_yc ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"yc\" attribute" << std::endl;
1016 if (node_zz ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"zz\" attribute" << std::endl;
1017 if (node_za ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"za\" attribute" << std::endl;
1018 if (node_zb ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"zb\" attribute" << std::endl;
1019 if (node_zc ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"zc\" attribute" << std::endl;
1020 if (node_aa ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"aa\" attribute" << std::endl;
1021 if (node_ab ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"ab\" attribute" << std::endl;
1022 if (node_ac ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"ac\" attribute" << std::endl;
1023 if (node_bb ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"bb\" attribute" << std::endl;
1024 if (node_bc ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"bc\" attribute" << std::endl;
1025 if (node_cc ==
NULL)
throw cms::Exception(
"XMLException") <<
"<setsurveyerr> is missing required \"cc\" attribute" << std::endl;
1028 matrix6x6(0,0) =
parseDouble(node_xx->getValue(),
"xx");
1029 matrix6x6(0,1) =
parseDouble(node_xy->getValue(),
"xy");
1030 matrix6x6(0,2) =
parseDouble(node_xz->getValue(),
"xz");
1031 matrix6x6(0,3) =
parseDouble(node_xa->getValue(),
"xa");
1032 matrix6x6(0,4) =
parseDouble(node_xb->getValue(),
"xb");
1033 matrix6x6(0,5) =
parseDouble(node_xc->getValue(),
"xc");
1034 matrix6x6(1,1) =
parseDouble(node_yy->getValue(),
"yy");
1035 matrix6x6(1,2) =
parseDouble(node_yz->getValue(),
"yz");
1036 matrix6x6(1,3) =
parseDouble(node_ya->getValue(),
"ya");
1037 matrix6x6(1,4) =
parseDouble(node_yb->getValue(),
"yb");
1038 matrix6x6(1,5) =
parseDouble(node_yc->getValue(),
"yc");
1039 matrix6x6(2,2) =
parseDouble(node_zz->getValue(),
"zz");
1040 matrix6x6(2,3) =
parseDouble(node_za->getValue(),
"za");
1041 matrix6x6(2,4) =
parseDouble(node_zb->getValue(),
"zb");
1042 matrix6x6(2,5) =
parseDouble(node_zc->getValue(),
"zc");
1043 matrix6x6(3,3) =
parseDouble(node_aa->getValue(),
"aa");
1044 matrix6x6(3,4) =
parseDouble(node_ab->getValue(),
"ab");
1045 matrix6x6(3,5) =
parseDouble(node_ac->getValue(),
"ac");
1046 matrix6x6(4,4) =
parseDouble(node_bb->getValue(),
"bb");
1047 matrix6x6(4,5) =
parseDouble(node_bc->getValue(),
"bc");
1048 matrix6x6(5,5) =
parseDouble(node_cc->getValue(),
"cc");
1050 for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1057 DOMAttr *node_x = node->getAttributeNode(
str_x);
1058 DOMAttr *node_y = node->getAttributeNode(
str_y);
1059 DOMAttr *node_z = node->getAttributeNode(
str_z);
1060 if (node_x ==
NULL)
throw cms::Exception(
"XMLException") <<
"<moveglobal> is missing required \"x\" attribute" << std::endl;
1061 if (node_y ==
NULL)
throw cms::Exception(
"XMLException") <<
"<moveglobal> is missing required \"y\" attribute" << std::endl;
1062 if (node_z ==
NULL)
throw cms::Exception(
"XMLException") <<
"<moveglobal> is missing required \"z\" attribute" << std::endl;
1069 for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1078 if (survey !=
NULL) {
1079 matrix6x6 = survey->
errors();
1086 DOMAttr *node_x = node->getAttributeNode(
str_x);
1087 DOMAttr *node_y = node->getAttributeNode(
str_y);
1088 DOMAttr *node_z = node->getAttributeNode(
str_z);
1089 if (node_x ==
NULL)
throw cms::Exception(
"XMLException") <<
"<movelocal> is missing required \"x\" attribute" << std::endl;
1090 if (node_y ==
NULL)
throw cms::Exception(
"XMLException") <<
"<movelocal> is missing required \"y\" attribute" << std::endl;
1091 if (node_z ==
NULL)
throw cms::Exception(
"XMLException") <<
"<movelocal> is missing required \"z\" attribute" << std::endl;
1098 for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1102 ali->
move(globalVector);
1108 if (survey !=
NULL) {
1109 matrix6x6 = survey->
errors();
1116 DOMAttr *node_axisx = node->getAttributeNode(
str_axisx);
1117 DOMAttr *node_axisy = node->getAttributeNode(
str_axisy);
1118 DOMAttr *node_axisz = node->getAttributeNode(
str_axisz);
1119 DOMAttr *node_angle = node->getAttributeNode(
str_angle);
1120 if (node_axisx ==
NULL)
throw cms::Exception(
"XMLException") <<
"<rotatelocal> is missing required \"axisx\" attribute" << std::endl;
1121 if (node_axisy ==
NULL)
throw cms::Exception(
"XMLException") <<
"<rotatelocal> is missing required \"axisy\" attribute" << std::endl;
1122 if (node_axisz ==
NULL)
throw cms::Exception(
"XMLException") <<
"<rotatelocal> is missing required \"axisz\" attribute" << std::endl;
1123 if (node_angle ==
NULL)
throw cms::Exception(
"XMLException") <<
"<rotatelocal> is missing required \"angle\" attribute" << std::endl;
1131 for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1140 if (survey !=
NULL) {
1141 matrix6x6 = survey->
errors();
1148 DOMAttr *node_rphi = node->getAttributeNode(
str_rphi);
1149 DOMAttr *node_phi = node->getAttributeNode(
str_phi);
1150 if (node_rphi ==
NULL && node_phi ==
NULL)
throw cms::Exception(
"XMLException") <<
"<rotatebeamline> is missing required \"*phi\" attribute" << std::endl;
1151 if (node_rphi !=
NULL && node_phi !=
NULL)
throw cms::Exception(
"XMLException") <<
"<rotatebeamline> can't have both an \"rphi\" and a \"phi\" attribute" << std::endl;
1154 if (node_rphi !=
NULL) {
1155 value =
parseDouble(node_rphi->getValue(),
"rphi");
1161 for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1167 double phi0 = pos.
phi();
1168 double deltaphi =
value;
1169 if (node_rphi !=
NULL) deltaphi = value /
radius;
1173 radius * (
sin(phi0 + deltaphi) -
sin(phi0)),
1180 if (survey !=
NULL) {
1181 matrix6x6 = survey->
errors();
1188 DOMAttr *node_x = node->getAttributeNode(
str_x);
1189 DOMAttr *node_y = node->getAttributeNode(
str_y);
1190 DOMAttr *node_z = node->getAttributeNode(
str_z);
1191 DOMAttr *node_angle = node->getAttributeNode(
str_angle);
1192 if (node_x ==
NULL)
throw cms::Exception(
"XMLException") <<
"<rotateglobalaxis> is missing required \"x\" attribute" << std::endl;
1193 if (node_y ==
NULL)
throw cms::Exception(
"XMLException") <<
"<rotateglobalaxis> is missing required \"y\" attribute" << std::endl;
1194 if (node_z ==
NULL)
throw cms::Exception(
"XMLException") <<
"<rotateglobalaxis> is missing required \"z\" attribute" << std::endl;
1195 if (node_angle ==
NULL)
throw cms::Exception(
"XMLException") <<
"<rotateglobalaxis> is missing required \"angle\" attribute" << std::endl;
1202 for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1208 double aprime = x/
sqrt(x*x + y*y + z*z);
1209 double bprime = y/
sqrt(x*x + y*y + z*z);
1210 double cprime = z/
sqrt(x*x + y*y + z*z);
1211 double q0 =
cos(angle/2.);
1212 double q1 =
sin(angle/2.) * aprime;
1213 double q2 =
sin(angle/2.) * bprime;
1214 double q3 =
sin(angle/2.) * cprime;
1216 double pos2x = (q0*q0 + q1*q1 - q2*q2 - q3*q3) * pos.
x() + 2.*(q1*q2 - q0*q3) * pos.
y() + 2.*(q1*q3 + q0*
q2) * pos.
z();
1217 double pos2y = 2.*(q2*q1 + q0*q3) * pos.
x() + (q0*q0 - q1*q1 + q2*q2 - q3*q3) * pos.
y() + 2.*(q2*q3 - q0*
q1) * pos.
z();
1218 double pos2z = 2.*(q3*q1 - q0*
q2) * pos.
x() + 2.*(q3*q2 + q0*
q1) * pos.
y() + (q0*q0 - q1*q1 - q2*q2 + q3*q3) * pos.
z();
1220 double movex = pos2x - pos.
x();
1221 double movey = pos2y - pos.
y();
1222 double movez = pos2z - pos.
z();
1229 if (survey !=
NULL) {
1230 matrix6x6 = survey->
errors();
TkRotation< Scalar > RotationType
virtual void rotateAroundLocalAxis(const LocalVector &axis, Scalar radians)
Rotation around arbitratry local axis.
Sin< T >::type sin(const T &t)
align::Alignables DTBarrel()
const SurveyDet * survey() const
Return survey info.
Geom::Phi< T > phi() const
const RotationType & globalRotation() const
Return the global orientation of the object.
virtual void move(const GlobalVector &displacement)=0
Movement with respect to the global reference frame.
const align::ErrorMatrix & errors() const
virtual void rotateAroundGlobalAxis(const GlobalVector &axis, Scalar radians)
Rotation around arbitratry global axis.
virtual void rotateAroundGlobalZ(Scalar radians)
Rotation around global z-axis.
uint32_t rawId() const
get the raw id
Point3DBase< Scalar, GlobalTag > PositionType
align::Alignables CSCEndcaps()
void rectify(RotationType &)
Correct a rotation matrix for rounding errors.
Cos< T >::type cos(const T &t)
virtual StructureType alignableObjectId() const =0
Return the alignable type identifier.
EulerAngles toAngles(const RotationType &)
Convert rotation matrix to angles about x-, y-, z-axes (frame rotation).
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
AlgebraicVector EulerAngles
virtual void rotateInGlobalFrame(const RotationType &rotation)=0
void setSurvey(const SurveyDet *)
Set survey info.
std::vector< Alignable * > Alignables
TkRotation transposed() const
align::GlobalPoints toGlobal(const align::LocalPoints &) const
Return in global coord given a set of local points.
RotationType toMatrix(const EulerAngles &)
Convert rotation angles about x-, y-, z-axes to matrix.
const PositionType & globalPosition() const
Return the global position of the object.
math::Error< 6 >::type ErrorMatrix
static const char * idToString(align::StructureType type)
Basic3DVector< T > multiplyInverse(const Basic3DVector< T > &v) const
Constructor of the full muon geometry.
const BasicVectorType & basicVector() const
Alignable * mother() const
Return pointer to container alignable (if any)
Power< A, B >::type pow(const A &a, const B &b)
T angle(T x1, T y1, T z1, T x2, T y2, T z2)