CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
MuonAlignmentInputXML Class Reference

#include <Alignment/MuonAlignment/interface/MuonAlignmentInputXML.h>

Inheritance diagram for MuonAlignmentInputXML:
MuonAlignmentInputMethod

Public Member Functions

 MuonAlignmentInputXML (std::string fileName)
 
virtual AlignableMuonnewAlignableMuon (const edm::EventSetup &iSetup) const
 
virtual ~MuonAlignmentInputXML ()
 
- Public Member Functions inherited from MuonAlignmentInputMethod
 MuonAlignmentInputMethod ()
 
virtual ~MuonAlignmentInputMethod ()
 

Private Member Functions

void do_moveglobal (const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
 
void do_movelocal (const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
 
void do_rotatebeamline (const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
 
void do_rotateglobalaxis (const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
 
void do_rotatelocal (const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
 
void do_setape (const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
 
void do_setposition (const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
 
void do_setsurveyerr (const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
 
void fillAliToIdeal (std::map< Alignable *, Alignable * > &alitoideal, const align::Alignables &alignables, const align::Alignables &ideals) const
 
AlignablegetCSCnode (align::StructureType structureType, std::map< unsigned int, Alignable * > &alignableNavigator, const XERCES_CPP_NAMESPACE::DOMElement *node) const
 
AlignablegetDTnode (align::StructureType structureType, std::map< unsigned int, Alignable * > &alignableNavigator, const XERCES_CPP_NAMESPACE::DOMElement *node) const
 
AlignablegetNode (std::map< unsigned int, Alignable * > &alignableNavigator, const XERCES_CPP_NAMESPACE::DOMElement *node) const
 
 MuonAlignmentInputXML (const MuonAlignmentInputXML &)
 
const MuonAlignmentInputXMLoperator= (const MuonAlignmentInputXML &)
 
double parseDouble (const XMLCh *str, const char *attribute) const
 
void recursiveGetId (std::map< unsigned int, Alignable * > &alignableNavigator, const align::Alignables &alignables) const
 
void set_one_position (Alignable *ali, const align::PositionType &pos, const align::RotationType &rot) const
 

Private Attributes

std::string m_fileName
 
XMLCh * str_aa
 
XMLCh * str_ab
 
XMLCh * str_ac
 
XMLCh * str_alpha
 
XMLCh * str_angle
 
XMLCh * str_axisx
 
XMLCh * str_axisy
 
XMLCh * str_axisz
 
XMLCh * str_bb
 
XMLCh * str_bc
 
XMLCh * str_beta
 
XMLCh * str_cc
 
XMLCh * str_chamber
 
XMLCh * str_collection
 
XMLCh * str_container
 
XMLCh * str_CSCChamber
 
XMLCh * str_CSCEndcap
 
XMLCh * str_CSCLayer
 
XMLCh * str_CSCRing
 
XMLCh * str_CSCStation
 
XMLCh * str_DTBarrel
 
XMLCh * str_DTChamber
 
XMLCh * str_DTLayer
 
XMLCh * str_DTStation
 
XMLCh * str_DTSuperLayer
 
XMLCh * str_DTWheel
 
XMLCh * str_endcap
 
XMLCh * str_gamma
 
XMLCh * str_ideal
 
XMLCh * str_layer
 
XMLCh * str_moveglobal
 
XMLCh * str_movelocal
 
XMLCh * str_name
 
XMLCh * str_none
 
XMLCh * str_operation
 
XMLCh * str_phi
 
XMLCh * str_phix
 
XMLCh * str_phiy
 
XMLCh * str_phiz
 
XMLCh * str_rawId
 
XMLCh * str_relativeto
 
XMLCh * str_ring
 
XMLCh * str_rotatebeamline
 
XMLCh * str_rotateglobalaxis
 
XMLCh * str_rotatelocal
 
XMLCh * str_rphi
 
XMLCh * str_sector
 
XMLCh * str_setape
 
XMLCh * str_setposition
 
XMLCh * str_setsurveyerr
 
XMLCh * str_station
 
XMLCh * str_superlayer
 
XMLCh * str_wheel
 
XMLCh * str_x
 
XMLCh * str_xa
 
XMLCh * str_xb
 
XMLCh * str_xc
 
XMLCh * str_xx
 
XMLCh * str_xy
 
XMLCh * str_xz
 
XMLCh * str_y
 
XMLCh * str_ya
 
XMLCh * str_yb
 
XMLCh * str_yc
 
XMLCh * str_yy
 
XMLCh * str_yz
 
XMLCh * str_z
 
XMLCh * str_za
 
XMLCh * str_zb
 
XMLCh * str_zc
 
XMLCh * str_zz
 

Additional Inherited Members

- Protected Member Functions inherited from MuonAlignmentInputMethod
boost::shared_ptr< CSCGeometryidealCSCGeometry (const edm::EventSetup &iSetup) const
 
boost::shared_ptr< DTGeometryidealDTGeometry (const edm::EventSetup &iSetup) const
 

Detailed Description

Description: <one line="" class="" summary>="">

Usage: <usage>

Definition at line 34 of file MuonAlignmentInputXML.h.

Constructor & Destructor Documentation

XERCES_CPP_NAMESPACE_USE MuonAlignmentInputXML::MuonAlignmentInputXML ( std::string  fileName)

Definition at line 46 of file MuonAlignmentInputXML.cc.

References str_aa, str_ab, str_ac, str_alpha, str_angle, str_axisx, str_axisy, str_axisz, str_bb, str_bc, str_beta, str_cc, str_chamber, str_collection, str_container, str_CSCChamber, str_CSCEndcap, str_CSCLayer, str_CSCRing, str_CSCStation, str_DTBarrel, str_DTChamber, str_DTLayer, str_DTStation, str_DTSuperLayer, str_DTWheel, str_endcap, str_gamma, str_ideal, str_layer, str_moveglobal, str_movelocal, str_name, str_none, str_operation, str_phi, str_phix, str_phiy, str_phiz, str_rawId, str_relativeto, str_ring, str_rotatebeamline, str_rotateglobalaxis, str_rotatelocal, str_rphi, str_sector, str_setape, str_setposition, str_setsurveyerr, str_station, str_superlayer, str_wheel, str_x, str_xa, str_xb, str_xc, str_xx, str_xy, str_xz, str_y, str_ya, str_yb, str_yc, str_yy, str_yz, str_z, str_za, str_zb, str_zc, and str_zz.

48 {
49  str_operation = XMLString::transcode("operation");
50  str_collection = XMLString::transcode("collection");
51  str_name = XMLString::transcode("name");
52  str_DTBarrel = XMLString::transcode("DTBarrel");
53  str_DTWheel = XMLString::transcode("DTWheel");
54  str_DTStation = XMLString::transcode("DTStation");
55  str_DTChamber = XMLString::transcode("DTChamber");
56  str_DTSuperLayer = XMLString::transcode("DTSuperLayer");
57  str_DTLayer = XMLString::transcode("DTLayer");
58  str_CSCEndcap = XMLString::transcode("CSCEndcap");
59  str_CSCStation = XMLString::transcode("CSCStation");
60  str_CSCRing = XMLString::transcode("CSCRing");
61  str_CSCChamber = XMLString::transcode("CSCChamber");
62  str_CSCLayer = XMLString::transcode("CSCLayer");
63  str_setposition = XMLString::transcode("setposition");
64  str_setape = XMLString::transcode("setape");
65  str_setsurveyerr = XMLString::transcode("setsurveyerr");
66  str_moveglobal = XMLString::transcode("moveglobal");
67  str_movelocal = XMLString::transcode("movelocal");
68  str_rotatelocal = XMLString::transcode("rotatelocal");
69  str_rotatebeamline = XMLString::transcode("rotatebeamline");
70  str_rotateglobalaxis = XMLString::transcode("rotateglobalaxis");
71  str_relativeto = XMLString::transcode("relativeto");
72  str_rawId = XMLString::transcode("rawId");
73  str_wheel = XMLString::transcode("wheel");
74  str_station = XMLString::transcode("station");
75  str_sector = XMLString::transcode("sector");
76  str_superlayer = XMLString::transcode("superlayer");
77  str_layer = XMLString::transcode("layer");
78  str_endcap = XMLString::transcode("endcap");
79  str_ring = XMLString::transcode("ring");
80  str_chamber = XMLString::transcode("chamber");
81  str_axisx = XMLString::transcode("axisx");
82  str_axisy = XMLString::transcode("axisy");
83  str_axisz = XMLString::transcode("axisz");
84  str_angle = XMLString::transcode("angle");
85  str_x = XMLString::transcode("x");
86  str_y = XMLString::transcode("y");
87  str_z = XMLString::transcode("z");
88  str_phix = XMLString::transcode("phix");
89  str_phiy = XMLString::transcode("phiy");
90  str_phiz = XMLString::transcode("phiz");
91  str_alpha = XMLString::transcode("alpha");
92  str_beta = XMLString::transcode("beta");
93  str_gamma = XMLString::transcode("gamma");
94  str_rphi = XMLString::transcode("rphi");
95  str_phi = XMLString::transcode("phi");
96  str_xx = XMLString::transcode("xx");
97  str_xy = XMLString::transcode("xy");
98  str_xz = XMLString::transcode("xz");
99  str_xa = XMLString::transcode("xa");
100  str_xb = XMLString::transcode("xb");
101  str_xc = XMLString::transcode("xc");
102  str_yy = XMLString::transcode("yy");
103  str_yz = XMLString::transcode("yz");
104  str_ya = XMLString::transcode("ya");
105  str_yb = XMLString::transcode("yb");
106  str_yc = XMLString::transcode("yc");
107  str_zz = XMLString::transcode("zz");
108  str_za = XMLString::transcode("za");
109  str_zb = XMLString::transcode("zb");
110  str_zc = XMLString::transcode("zc");
111  str_aa = XMLString::transcode("aa");
112  str_ab = XMLString::transcode("ab");
113  str_ac = XMLString::transcode("ac");
114  str_bb = XMLString::transcode("bb");
115  str_bc = XMLString::transcode("bc");
116  str_cc = XMLString::transcode("cc");
117  str_none = XMLString::transcode("none");
118  str_ideal = XMLString::transcode("ideal");
119  str_container = XMLString::transcode("container");
120 }
MuonAlignmentInputXML::~MuonAlignmentInputXML ( )
virtual

Definition at line 127 of file MuonAlignmentInputXML.cc.

References fetchall_from_DQM_v2::release, str_aa, str_ab, str_ac, str_alpha, str_angle, str_axisx, str_axisy, str_axisz, str_bb, str_bc, str_beta, str_cc, str_chamber, str_collection, str_container, str_CSCChamber, str_CSCEndcap, str_CSCLayer, str_CSCRing, str_CSCStation, str_DTBarrel, str_DTChamber, str_DTLayer, str_DTStation, str_DTSuperLayer, str_DTWheel, str_endcap, str_gamma, str_ideal, str_layer, str_moveglobal, str_movelocal, str_name, str_none, str_operation, str_phi, str_phix, str_phiy, str_phiz, str_rawId, str_relativeto, str_ring, str_rotatebeamline, str_rotateglobalaxis, str_rotatelocal, str_rphi, str_sector, str_setape, str_setposition, str_setsurveyerr, str_station, str_superlayer, str_wheel, str_x, str_xa, str_xb, str_xc, str_xx, str_xy, str_xz, str_y, str_ya, str_yb, str_yc, str_yy, str_yz, str_z, str_za, str_zb, str_zc, and str_zz.

127  {
199 }
MuonAlignmentInputXML::MuonAlignmentInputXML ( const MuonAlignmentInputXML )
private

Member Function Documentation

void MuonAlignmentInputXML::do_moveglobal ( const XERCES_CPP_NAMESPACE::DOMElement *  node,
std::map< Alignable *, bool > &  aliset,
std::map< Alignable *, Alignable * > &  alitoideal 
) const
private

Definition at line 1015 of file MuonAlignmentInputXML.cc.

References SurveyDet::errors(), edm::hlt::Exception, Alignable::move(), NULL, parseDouble(), Alignable::setSurvey(), str_x, str_y, str_z, Alignable::surface(), Alignable::survey(), x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

Referenced by newAlignableMuon().

1015  {
1016  DOMAttr *node_x = node->getAttributeNode(str_x);
1017  DOMAttr *node_y = node->getAttributeNode(str_y);
1018  DOMAttr *node_z = node->getAttributeNode(str_z);
1019  if (node_x == NULL) throw cms::Exception("XMLException") << "<moveglobal> is missing required \"x\" attribute" << std::endl;
1020  if (node_y == NULL) throw cms::Exception("XMLException") << "<moveglobal> is missing required \"y\" attribute" << std::endl;
1021  if (node_z == NULL) throw cms::Exception("XMLException") << "<moveglobal> is missing required \"z\" attribute" << std::endl;
1022 
1023  double x = parseDouble(node_x->getValue(), "x");
1024  double y = parseDouble(node_y->getValue(), "y");
1025  double z = parseDouble(node_z->getValue(), "z");
1026  align::GlobalVector vect(x, y, z);
1027 
1028  for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1029  Alignable *ali = aliiter->first;
1030 
1031  ali->move(vect);
1032 
1033  align::ErrorMatrix matrix6x6 = ROOT::Math::SMatrixIdentity();
1034  matrix6x6 *= 1000.; // initial assumption: infinitely weak constraint
1035 
1036  const SurveyDet *survey = ali->survey();
1037  if (survey != NULL) {
1038  matrix6x6 = survey->errors(); // save the constraint information
1039  }
1040  ali->setSurvey(new SurveyDet(ali->surface(), matrix6x6));
1041  } // end loop over alignables
1042 }
double parseDouble(const XMLCh *str, const char *attribute) const
const SurveyDet * survey() const
Return survey info.
Definition: Alignable.h:205
#define NULL
Definition: scimark2.h:8
virtual void move(const GlobalVector &displacement)=0
Movement with respect to the global reference frame.
const align::ErrorMatrix & errors() const
Definition: SurveyDet.h:73
float float float z
tuple node
Definition: Node.py:50
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:126
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:268
math::Error< 6 >::type ErrorMatrix
Definition: Definitions.h:39
Definition: DDAxes.h:10
void MuonAlignmentInputXML::do_movelocal ( const XERCES_CPP_NAMESPACE::DOMElement *  node,
std::map< Alignable *, bool > &  aliset,
std::map< Alignable *, Alignable * > &  alitoideal 
) const
private

Definition at line 1044 of file MuonAlignmentInputXML.cc.

References SurveyDet::errors(), edm::hlt::Exception, Alignable::move(), NULL, parseDouble(), Alignable::setSurvey(), str_x, str_y, str_z, Alignable::surface(), Alignable::survey(), AlignableSurface::toGlobal(), x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

Referenced by newAlignableMuon().

1044  {
1045  DOMAttr *node_x = node->getAttributeNode(str_x);
1046  DOMAttr *node_y = node->getAttributeNode(str_y);
1047  DOMAttr *node_z = node->getAttributeNode(str_z);
1048  if (node_x == NULL) throw cms::Exception("XMLException") << "<movelocal> is missing required \"x\" attribute" << std::endl;
1049  if (node_y == NULL) throw cms::Exception("XMLException") << "<movelocal> is missing required \"y\" attribute" << std::endl;
1050  if (node_z == NULL) throw cms::Exception("XMLException") << "<movelocal> is missing required \"z\" attribute" << std::endl;
1051 
1052  double x = parseDouble(node_x->getValue(), "x");
1053  double y = parseDouble(node_y->getValue(), "y");
1054  double z = parseDouble(node_z->getValue(), "z");
1055  align::LocalVector vect(x, y, z);
1056 
1057  for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1058  Alignable *ali = aliiter->first;
1059 
1060  align::GlobalVector globalVector = ali->surface().toGlobal(vect);
1061  ali->move(globalVector);
1062 
1063  align::ErrorMatrix matrix6x6 = ROOT::Math::SMatrixIdentity();
1064  matrix6x6 *= 1000.; // initial assumption: infinitely weak constraint
1065 
1066  const SurveyDet *survey = ali->survey();
1067  if (survey != NULL) {
1068  matrix6x6 = survey->errors(); // save the constraint information
1069  }
1070  ali->setSurvey(new SurveyDet(ali->surface(), matrix6x6));
1071  } // end loop over alignables
1072 }
double parseDouble(const XMLCh *str, const char *attribute) const
const SurveyDet * survey() const
Return survey info.
Definition: Alignable.h:205
#define NULL
Definition: scimark2.h:8
virtual void move(const GlobalVector &displacement)=0
Movement with respect to the global reference frame.
const align::ErrorMatrix & errors() const
Definition: SurveyDet.h:73
float float float z
tuple node
Definition: Node.py:50
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:126
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:268
align::GlobalPoints toGlobal(const align::LocalPoints &) const
Return in global coord given a set of local points.
math::Error< 6 >::type ErrorMatrix
Definition: Definitions.h:39
Definition: DDAxes.h:10
void MuonAlignmentInputXML::do_rotatebeamline ( const XERCES_CPP_NAMESPACE::DOMElement *  node,
std::map< Alignable *, bool > &  aliset,
std::map< Alignable *, Alignable * > &  alitoideal 
) const
private

Definition at line 1106 of file MuonAlignmentInputXML.cc.

References funct::cos(), SurveyDet::errors(), edm::hlt::Exception, Alignable::move(), NULL, parseDouble(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), CosmicsPD_Skims::radius, Alignable::rotateAroundGlobalZ(), Alignable::setSurvey(), funct::sin(), str_phi, str_rphi, Alignable::surface(), Alignable::survey(), AlignableSurface::toGlobal(), and relativeConstraints::value.

Referenced by newAlignableMuon().

1106  {
1107  DOMAttr *node_rphi = node->getAttributeNode(str_rphi);
1108  DOMAttr *node_phi = node->getAttributeNode(str_phi);
1109  if (node_rphi == NULL && node_phi == NULL) throw cms::Exception("XMLException") << "<rotatebeamline> is missing required \"*phi\" attribute" << std::endl;
1110  if (node_rphi != NULL && node_phi != NULL) throw cms::Exception("XMLException") << "<rotatebeamline> can't have both an \"rphi\" and a \"phi\" attribute" << std::endl;
1111 
1112  double value;
1113  if (node_rphi != NULL) {
1114  value = parseDouble(node_rphi->getValue(), "rphi");
1115  }
1116  else {
1117  value = parseDouble(node_phi->getValue(), "phi");
1118  }
1119 
1120  for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1121  Alignable *ali = aliiter->first;
1122 
1124 
1125  double radius = pos.perp();
1126  double phi0 = pos.phi();
1127  double deltaphi = value;
1128  if (node_rphi != NULL) deltaphi = value / radius;
1129 
1130  ali->rotateAroundGlobalZ(deltaphi);
1131  ali->move(align::GlobalVector(radius * (cos(phi0 + deltaphi) - cos(phi0)),
1132  radius * (sin(phi0 + deltaphi) - sin(phi0)),
1133  0.));
1134 
1135  align::ErrorMatrix matrix6x6 = ROOT::Math::SMatrixIdentity();
1136  matrix6x6 *= 1000.; // initial assumption: infinitely weak constraint
1137 
1138  const SurveyDet *survey = ali->survey();
1139  if (survey != NULL) {
1140  matrix6x6 = survey->errors(); // save the constraint information
1141  }
1142  ali->setSurvey(new SurveyDet(ali->surface(), matrix6x6));
1143  } // end loop over alignables
1144 }
T perp() const
Definition: PV3DBase.h:72
double parseDouble(const XMLCh *str, const char *attribute) const
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
const SurveyDet * survey() const
Return survey info.
Definition: Alignable.h:205
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
#define NULL
Definition: scimark2.h:8
virtual void move(const GlobalVector &displacement)=0
Movement with respect to the global reference frame.
const align::ErrorMatrix & errors() const
Definition: SurveyDet.h:73
tuple node
Definition: Node.py:50
virtual void rotateAroundGlobalZ(Scalar radians)
Rotation around global z-axis.
Definition: Alignable.cc:174
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:126
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:268
align::GlobalPoints toGlobal(const align::LocalPoints &) const
Return in global coord given a set of local points.
math::Error< 6 >::type ErrorMatrix
Definition: Definitions.h:39
void MuonAlignmentInputXML::do_rotateglobalaxis ( const XERCES_CPP_NAMESPACE::DOMElement *  node,
std::map< Alignable *, bool > &  aliset,
std::map< Alignable *, Alignable * > &  alitoideal 
) const
private

Definition at line 1146 of file MuonAlignmentInputXML.cc.

References angle(), funct::cos(), SurveyDet::errors(), edm::hlt::Exception, Alignable::move(), NULL, parseDouble(), SiStripMonitorClusterAlca_cfi::q0, q1, q2, Alignable::rotateAroundGlobalAxis(), Alignable::setSurvey(), funct::sin(), mathSSE::sqrt(), str_angle, str_x, str_y, str_z, Alignable::surface(), Alignable::survey(), AlignableSurface::toGlobal(), x, PV3DBase< T, PVType, FrameType >::x(), detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by newAlignableMuon().

1146  {
1147  DOMAttr *node_x = node->getAttributeNode(str_x);
1148  DOMAttr *node_y = node->getAttributeNode(str_y);
1149  DOMAttr *node_z = node->getAttributeNode(str_z);
1150  DOMAttr *node_angle = node->getAttributeNode(str_angle);
1151  if (node_x == NULL) throw cms::Exception("XMLException") << "<rotateglobalaxis> is missing required \"x\" attribute" << std::endl;
1152  if (node_y == NULL) throw cms::Exception("XMLException") << "<rotateglobalaxis> is missing required \"y\" attribute" << std::endl;
1153  if (node_z == NULL) throw cms::Exception("XMLException") << "<rotateglobalaxis> is missing required \"z\" attribute" << std::endl;
1154  if (node_angle == NULL) throw cms::Exception("XMLException") << "<rotateglobalaxis> is missing required \"angle\" attribute" << std::endl;
1155 
1156  double x = parseDouble(node_x->getValue(), "x");
1157  double y = parseDouble(node_y->getValue(), "y");
1158  double z = parseDouble(node_z->getValue(), "z");
1159  double angle = parseDouble(node_angle->getValue(), "angle");
1160 
1161  for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1162  Alignable *ali = aliiter->first;
1164 
1165  ali->rotateAroundGlobalAxis(align::GlobalVector(x, y, z), angle);
1166 
1167  double aprime = x/sqrt(x*x + y*y + z*z);
1168  double bprime = y/sqrt(x*x + y*y + z*z);
1169  double cprime = z/sqrt(x*x + y*y + z*z);
1170  double q0 = cos(angle/2.);
1171  double q1 = sin(angle/2.) * aprime;
1172  double q2 = sin(angle/2.) * bprime;
1173  double q3 = sin(angle/2.) * cprime;
1174 
1175  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();
1176  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();
1177  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();
1178 
1179  double movex = pos2x - pos.x();
1180  double movey = pos2y - pos.y();
1181  double movez = pos2z - pos.z();
1182  ali->move(align::GlobalVector(movex, movey, movez));
1183 
1184  align::ErrorMatrix matrix6x6 = ROOT::Math::SMatrixIdentity();
1185  matrix6x6 *= 1000.; // initial assumption: infinitely weak constraint
1186 
1187  const SurveyDet *survey = ali->survey();
1188  if (survey != NULL) {
1189  matrix6x6 = survey->errors(); // save the constraint information
1190  }
1191  ali->setSurvey(new SurveyDet(ali->surface(), matrix6x6));
1192  } // end loop over alignables
1193 }
double parseDouble(const XMLCh *str, const char *attribute) const
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
const SurveyDet * survey() const
Return survey info.
Definition: Alignable.h:205
T y() const
Definition: PV3DBase.h:63
#define NULL
Definition: scimark2.h:8
virtual void move(const GlobalVector &displacement)=0
Movement with respect to the global reference frame.
const align::ErrorMatrix & errors() const
Definition: SurveyDet.h:73
virtual void rotateAroundGlobalAxis(const GlobalVector &axis, Scalar radians)
Rotation around arbitratry global axis.
Definition: Alignable.cc:104
float float float z
double q2[4]
Definition: TauolaWrapper.h:88
tuple node
Definition: Node.py:50
T sqrt(T t)
Definition: SSEVec.h:48
T z() const
Definition: PV3DBase.h:64
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:126
double q1[4]
Definition: TauolaWrapper.h:87
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:268
align::GlobalPoints toGlobal(const align::LocalPoints &) const
Return in global coord given a set of local points.
math::Error< 6 >::type ErrorMatrix
Definition: Definitions.h:39
Definition: DDAxes.h:10
T x() const
Definition: PV3DBase.h:62
T angle(T x1, T y1, T z1, T x2, T y2, T z2)
Definition: angle.h:11
void MuonAlignmentInputXML::do_rotatelocal ( const XERCES_CPP_NAMESPACE::DOMElement *  node,
std::map< Alignable *, bool > &  aliset,
std::map< Alignable *, Alignable * > &  alitoideal 
) const
private

Definition at line 1074 of file MuonAlignmentInputXML.cc.

References angle(), SurveyDet::errors(), edm::hlt::Exception, NULL, parseDouble(), Alignable::rotateAroundLocalAxis(), Alignable::setSurvey(), str_angle, str_axisx, str_axisy, str_axisz, Alignable::surface(), Alignable::survey(), x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

Referenced by newAlignableMuon().

1074  {
1075  DOMAttr *node_axisx = node->getAttributeNode(str_axisx);
1076  DOMAttr *node_axisy = node->getAttributeNode(str_axisy);
1077  DOMAttr *node_axisz = node->getAttributeNode(str_axisz);
1078  DOMAttr *node_angle = node->getAttributeNode(str_angle);
1079  if (node_axisx == NULL) throw cms::Exception("XMLException") << "<rotatelocal> is missing required \"axisx\" attribute" << std::endl;
1080  if (node_axisy == NULL) throw cms::Exception("XMLException") << "<rotatelocal> is missing required \"axisy\" attribute" << std::endl;
1081  if (node_axisz == NULL) throw cms::Exception("XMLException") << "<rotatelocal> is missing required \"axisz\" attribute" << std::endl;
1082  if (node_angle == NULL) throw cms::Exception("XMLException") << "<rotatelocal> is missing required \"angle\" attribute" << std::endl;
1083 
1084  double x = parseDouble(node_axisx->getValue(), "x");
1085  double y = parseDouble(node_axisy->getValue(), "y");
1086  double z = parseDouble(node_axisz->getValue(), "z");
1087  double angle = parseDouble(node_angle->getValue(), "angle");
1088  align::LocalVector vect(x, y, z);
1089 
1090  for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1091  Alignable *ali = aliiter->first;
1092 
1093  ali->rotateAroundLocalAxis(vect, angle);
1094 
1095  align::ErrorMatrix matrix6x6 = ROOT::Math::SMatrixIdentity();
1096  matrix6x6 *= 1000.; // initial assumption: infinitely weak constraint
1097 
1098  const SurveyDet *survey = ali->survey();
1099  if (survey != NULL) {
1100  matrix6x6 = survey->errors(); // save the constraint information
1101  }
1102  ali->setSurvey(new SurveyDet(ali->surface(), matrix6x6));
1103  } // end loop over alignables
1104 }
virtual void rotateAroundLocalAxis(const LocalVector &axis, Scalar radians)
Rotation around arbitratry local axis.
Definition: Alignable.cc:113
double parseDouble(const XMLCh *str, const char *attribute) const
const SurveyDet * survey() const
Return survey info.
Definition: Alignable.h:205
#define NULL
Definition: scimark2.h:8
const align::ErrorMatrix & errors() const
Definition: SurveyDet.h:73
float float float z
tuple node
Definition: Node.py:50
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:126
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:268
math::Error< 6 >::type ErrorMatrix
Definition: Definitions.h:39
Definition: DDAxes.h:10
T angle(T x1, T y1, T z1, T x2, T y2, T z2)
Definition: angle.h:11
void MuonAlignmentInputXML::do_setape ( const XERCES_CPP_NAMESPACE::DOMElement *  node,
std::map< Alignable *, bool > &  aliset,
std::map< Alignable *, Alignable * > &  alitoideal 
) const
private

Definition at line 912 of file MuonAlignmentInputXML.cc.

References edm::hlt::Exception, NULL, parseDouble(), str_xx, str_xy, str_xz, str_yy, str_yz, and str_zz.

Referenced by newAlignableMuon().

912  {
913  DOMAttr *node_xx = node->getAttributeNode(str_xx);
914  DOMAttr *node_xy = node->getAttributeNode(str_xy);
915  DOMAttr *node_xz = node->getAttributeNode(str_xz);
916  DOMAttr *node_yy = node->getAttributeNode(str_yy);
917  DOMAttr *node_yz = node->getAttributeNode(str_yz);
918  DOMAttr *node_zz = node->getAttributeNode(str_zz);
919 
920  if (node_xx == NULL) throw cms::Exception("XMLException") << "<setape> is missing required \"xx\" attribute" << std::endl;
921  if (node_xy == NULL) throw cms::Exception("XMLException") << "<setape> is missing required \"xy\" attribute" << std::endl;
922  if (node_xz == NULL) throw cms::Exception("XMLException") << "<setape> is missing required \"xz\" attribute" << std::endl;
923  if (node_yy == NULL) throw cms::Exception("XMLException") << "<setape> is missing required \"yy\" attribute" << std::endl;
924  if (node_yz == NULL) throw cms::Exception("XMLException") << "<setape> is missing required \"yz\" attribute" << std::endl;
925  if (node_zz == NULL) throw cms::Exception("XMLException") << "<setape> is missing required \"zz\" attribute" << std::endl;
926 
927  AlgebraicSymMatrix33 matrix3x3;
928  matrix3x3(0,0) = parseDouble(node_xx->getValue(), "xx");
929  matrix3x3(0,1) = parseDouble(node_xy->getValue(), "xy");
930  matrix3x3(0,2) = parseDouble(node_xz->getValue(), "xz");
931  matrix3x3(1,1) = parseDouble(node_yy->getValue(), "yy");
932  matrix3x3(1,2) = parseDouble(node_yz->getValue(), "yz");
933  matrix3x3(2,2) = parseDouble(node_zz->getValue(), "zz");
934 
935  for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
936  // this sets APEs at this level and (since 2nd argument is true) all lower levels
937  aliiter->first->setAlignmentPositionError(AlignmentPositionError(matrix3x3), true);
938  }
939 }
double parseDouble(const XMLCh *str, const char *attribute) const
#define NULL
Definition: scimark2.h:8
tuple node
Definition: Node.py:50
ROOT::Math::SMatrix< double, 3, 3, ROOT::Math::MatRepSym< double, 3 > > AlgebraicSymMatrix33
void MuonAlignmentInputXML::do_setposition ( const XERCES_CPP_NAMESPACE::DOMElement *  node,
std::map< Alignable *, bool > &  aliset,
std::map< Alignable *, Alignable * > &  alitoideal 
) const
private

Definition at line 719 of file MuonAlignmentInputXML.cc.

References alpha, PV3DBase< T, PVType, FrameType >::basicVector(), beta, funct::cos(), edm::hlt::Exception, Alignable::globalPosition(), Alignable::globalRotation(), python.rootplot.argparse::message, Alignable::mother(), NULL, parseDouble(), makeMuonMisalignmentScenario::rot, set_one_position(), funct::sin(), str_alpha, str_beta, str_container, str_gamma, str_ideal, str_none, str_phix, str_phiy, str_phiz, str_relativeto, str_x, str_y, str_z, align::toAngles(), align::toMatrix(), TkRotation< T >::transposed(), x, PV3DBase< T, PVType, FrameType >::x(), TkRotation< T >::xx(), TkRotation< T >::xy(), TkRotation< T >::xz(), detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::y(), TkRotation< T >::yz(), detailsBasic3DVector::z, PV3DBase< T, PVType, FrameType >::z(), and TkRotation< T >::zz().

Referenced by newAlignableMuon().

719  {
720  DOMAttr *node_relativeto = node->getAttributeNode(str_relativeto);
721  if (node_relativeto == NULL) throw cms::Exception("XMLException") << "<setposition> is missing required \"relativeto\" attribute" << std::endl;
722 
723  int relativeto = 0;
724  if (XMLString::equals(node_relativeto->getValue(), str_none)) {
725  relativeto = 0;
726  }
727  else if (XMLString::equals(node_relativeto->getValue(), str_ideal)) {
728  relativeto = 1;
729  }
730  else if (XMLString::equals(node_relativeto->getValue(), str_container)) {
731  relativeto = 2;
732  }
733  else {
734  char *message = XMLString::transcode(node_relativeto->getValue());
735  throw cms::Exception("XMLException") << "relativeto must be \"none\", \"ideal\", or \"container\", not \"" << message << "\"" << std::endl;
736  }
737 
738  for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
739  // first reconstruct the old position: how it would look in this coordinate system
740 
741  Alignable *ali = aliiter->first;
742  Alignable *ideal = alitoideal[ali];
743 
744  align::PositionType oldpos = ali->globalPosition();
745  align::RotationType oldrot = ali->globalRotation();
746 
747  if (relativeto == 0) {}
748 
749  else if (relativeto == 1) {
750  align::PositionType idealPosition = ideal->globalPosition();
751  align::RotationType idealRotation = ideal->globalRotation();
752 
753  oldpos = align::PositionType(idealRotation * (oldpos.basicVector() - idealPosition.basicVector()));
754  oldrot = oldrot * idealRotation.transposed();
755  }
756 
757  else if (relativeto == 2 && ali->mother() != NULL) {
758  align::PositionType globalPosition = ali->mother()->globalPosition();
759  align::RotationType globalRotation = ali->mother()->globalRotation();
760 
761  oldpos = align::PositionType(globalRotation * (oldpos.basicVector() - globalPosition.basicVector()));
762  oldrot = oldrot * globalRotation.transposed();
763  }
764 
765  double x = oldpos.x();
766  double y = oldpos.y();
767  double z = oldpos.z();
768 
769  double phix = atan2(oldrot.yz(), oldrot.zz());
770  double phiy = asin(-oldrot.xz());
771  double phiz = atan2(oldrot.xy(), oldrot.xx());
772 
773  align::EulerAngles oldEulerAngles = align::toAngles(oldrot);
774  double alpha = oldEulerAngles(1);
775  double beta = oldEulerAngles(2);
776  double gamma = oldEulerAngles(3);
777 
778  // now get the new information; if it's incomplete, use the old position for those coordinates
779 
780  DOMAttr *node_x = node->getAttributeNode(str_x);
781  DOMAttr *node_y = node->getAttributeNode(str_y);
782  DOMAttr *node_z = node->getAttributeNode(str_z);
783 
784  if (node_x != NULL) x = parseDouble(node_x->getValue(), "x");
785  if (node_y != NULL) y = parseDouble(node_y->getValue(), "y");
786  if (node_z != NULL) z = parseDouble(node_z->getValue(), "z");
787  align::PositionType pos(x, y, z);
788 
789  DOMAttr *node_phix = node->getAttributeNode(str_phix);
790  DOMAttr *node_phiy = node->getAttributeNode(str_phiy);
791  DOMAttr *node_phiz = node->getAttributeNode(str_phiz);
792  DOMAttr *node_alpha = node->getAttributeNode(str_alpha);
793  DOMAttr *node_beta = node->getAttributeNode(str_beta);
794  DOMAttr *node_gamma = node->getAttributeNode(str_gamma);
796 
797  bool phixyz = (node_phix != NULL || node_phiy != NULL || node_phiz != NULL);
798  bool alphabetagamma = (node_alpha != NULL || node_beta != NULL || node_gamma != NULL);
799  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;
800  if (!phixyz && !alphabetagamma) alphabetagamma = true;
801 
802  if (phixyz) {
803  if (node_phix != NULL) phix = parseDouble(node_phix->getValue(), "phix");
804  if (node_phiy != NULL) phiy = parseDouble(node_phiy->getValue(), "phiy");
805  if (node_phiz != NULL) phiz = parseDouble(node_phiz->getValue(), "phiz");
806 
807  // the angle convention originally used in alignment, also known as "non-standard Euler angles with a Z-Y-X convention"
808  // this also gets the sign convention right
809  align::RotationType rotX( 1., 0., 0.,
810  0., cos(phix), sin(phix),
811  0., -sin(phix), cos(phix));
812  align::RotationType rotY( cos(phiy), 0., -sin(phiy),
813  0., 1., 0.,
814  sin(phiy), 0., cos(phiy));
815  align::RotationType rotZ( cos(phiz), sin(phiz), 0.,
816  -sin(phiz), cos(phiz), 0.,
817  0., 0., 1.);
818 
819  rot = rotX * rotY * rotZ;
820  }
821 
822  else if (alphabetagamma) {
823  if (node_alpha != NULL) alpha = parseDouble(node_alpha->getValue(), "alpha");
824  if (node_beta != NULL) beta = parseDouble(node_beta->getValue(), "beta");
825  if (node_gamma != NULL) gamma = parseDouble(node_gamma->getValue(), "gamma");
826 
827  // standard Euler angles (how they're internally stored in the database)
828  align::EulerAngles eulerAngles(3);
829  eulerAngles(1) = alpha;
830  eulerAngles(2) = beta;
831  eulerAngles(3) = gamma;
832  rot = align::RotationType(align::toMatrix(eulerAngles));
833  }
834 
835  else assert(false); // see above
836 
837  if (relativeto == 0) {
838  set_one_position(aliiter->first, pos, rot);
839  } // end relativeto="none"
840 
841  else if (relativeto == 1) {
842  Alignable *ali = aliiter->first;
843  Alignable *ideal = alitoideal[ali];
844 
845  align::PositionType idealPosition = ideal->globalPosition();
846  align::RotationType idealRotation = ideal->globalRotation();
847  align::PositionType newpos = align::PositionType(idealRotation.transposed() * pos.basicVector() + idealPosition.basicVector());
848  align::RotationType newrot = rot * idealRotation;
849 
850  set_one_position(ali, newpos, newrot);
851  } // end relativeto="ideal"
852 
853  else if (relativeto == 2) {
854  Alignable *ali = aliiter->first;
855  Alignable *container = ali->mother();
856 
857  if (container != NULL) {
858  align::PositionType globalPosition = container->globalPosition();
859  align::RotationType globalRotation = container->globalRotation();
860  align::PositionType newpos = align::PositionType(globalRotation.transposed() * pos.basicVector() + globalPosition.basicVector());
861  align::RotationType newrot = rot * globalRotation;
862  set_one_position(ali, newpos, newrot);
863  }
864  else {
865  set_one_position(ali, pos, rot);
866  }
867  } // end relativeto="container"
868 
869  } // end loop over alignables
870 }
T xx() const
const double beta
TkRotation< Scalar > RotationType
Definition: Definitions.h:29
float alpha
Definition: AMPTWrapper.h:95
double parseDouble(const XMLCh *str, const char *attribute) const
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
T y() const
Definition: PV3DBase.h:63
#define NULL
Definition: scimark2.h:8
const RotationType & globalRotation() const
Return the global orientation of the object.
Definition: Alignable.h:132
float float float z
tuple node
Definition: Node.py:50
Point3DBase< Scalar, GlobalTag > PositionType
Definition: Definitions.h:30
T xy() const
T zz() const
T z() const
Definition: PV3DBase.h:64
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
EulerAngles toAngles(const RotationType &)
Convert rotation matrix to angles about x-, y-, z-axes (frame rotation).
Definition: Utilities.cc:7
AlgebraicVector EulerAngles
Definition: Definitions.h:36
T xz() const
TkRotation transposed() const
RotationType toMatrix(const EulerAngles &)
Convert rotation angles about x-, y-, z-axes to matrix.
Definition: Utilities.cc:40
const PositionType & globalPosition() const
Return the global position of the object.
Definition: Alignable.h:129
Definition: DDAxes.h:10
void set_one_position(Alignable *ali, const align::PositionType &pos, const align::RotationType &rot) const
T x() const
Definition: PV3DBase.h:62
const BasicVectorType & basicVector() const
Definition: PV3DBase.h:56
Alignable * mother() const
Return pointer to container alignable (if any)
Definition: Alignable.h:85
T yz() const
void MuonAlignmentInputXML::do_setsurveyerr ( const XERCES_CPP_NAMESPACE::DOMElement *  node,
std::map< Alignable *, bool > &  aliset,
std::map< Alignable *, Alignable * > &  alitoideal 
) const
private

Definition at line 941 of file MuonAlignmentInputXML.cc.

References edm::hlt::Exception, NULL, parseDouble(), Alignable::setSurvey(), str_aa, str_ab, str_ac, str_bb, str_bc, str_cc, str_xa, str_xb, str_xc, str_xx, str_xy, str_xz, str_ya, str_yb, str_yc, str_yy, str_yz, str_za, str_zb, str_zc, str_zz, and Alignable::surface().

Referenced by newAlignableMuon().

941  {
942  DOMAttr *node_xx = node->getAttributeNode(str_xx);
943  DOMAttr *node_xy = node->getAttributeNode(str_xy);
944  DOMAttr *node_xz = node->getAttributeNode(str_xz);
945  DOMAttr *node_xa = node->getAttributeNode(str_xa);
946  DOMAttr *node_xb = node->getAttributeNode(str_xb);
947  DOMAttr *node_xc = node->getAttributeNode(str_xc);
948  DOMAttr *node_yy = node->getAttributeNode(str_yy);
949  DOMAttr *node_yz = node->getAttributeNode(str_yz);
950  DOMAttr *node_ya = node->getAttributeNode(str_ya);
951  DOMAttr *node_yb = node->getAttributeNode(str_yb);
952  DOMAttr *node_yc = node->getAttributeNode(str_yc);
953  DOMAttr *node_zz = node->getAttributeNode(str_zz);
954  DOMAttr *node_za = node->getAttributeNode(str_za);
955  DOMAttr *node_zb = node->getAttributeNode(str_zb);
956  DOMAttr *node_zc = node->getAttributeNode(str_zc);
957  DOMAttr *node_aa = node->getAttributeNode(str_aa);
958  DOMAttr *node_ab = node->getAttributeNode(str_ab);
959  DOMAttr *node_ac = node->getAttributeNode(str_ac);
960  DOMAttr *node_bb = node->getAttributeNode(str_bb);
961  DOMAttr *node_bc = node->getAttributeNode(str_bc);
962  DOMAttr *node_cc = node->getAttributeNode(str_cc);
963 
964  if (node_xx == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"xx\" attribute" << std::endl;
965  if (node_xy == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"xy\" attribute" << std::endl;
966  if (node_xz == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"xz\" attribute" << std::endl;
967  if (node_xa == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"xa\" attribute" << std::endl;
968  if (node_xb == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"xb\" attribute" << std::endl;
969  if (node_xc == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"xc\" attribute" << std::endl;
970  if (node_yy == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"yy\" attribute" << std::endl;
971  if (node_yz == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"yz\" attribute" << std::endl;
972  if (node_ya == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"ya\" attribute" << std::endl;
973  if (node_yb == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"yb\" attribute" << std::endl;
974  if (node_yc == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"yc\" attribute" << std::endl;
975  if (node_zz == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"zz\" attribute" << std::endl;
976  if (node_za == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"za\" attribute" << std::endl;
977  if (node_zb == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"zb\" attribute" << std::endl;
978  if (node_zc == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"zc\" attribute" << std::endl;
979  if (node_aa == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"aa\" attribute" << std::endl;
980  if (node_ab == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"ab\" attribute" << std::endl;
981  if (node_ac == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"ac\" attribute" << std::endl;
982  if (node_bb == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"bb\" attribute" << std::endl;
983  if (node_bc == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"bc\" attribute" << std::endl;
984  if (node_cc == NULL) throw cms::Exception("XMLException") << "<setsurveyerr> is missing required \"cc\" attribute" << std::endl;
985 
986  align::ErrorMatrix matrix6x6;
987  matrix6x6(0,0) = parseDouble(node_xx->getValue(), "xx");
988  matrix6x6(0,1) = parseDouble(node_xy->getValue(), "xy");
989  matrix6x6(0,2) = parseDouble(node_xz->getValue(), "xz");
990  matrix6x6(0,3) = parseDouble(node_xa->getValue(), "xa");
991  matrix6x6(0,4) = parseDouble(node_xb->getValue(), "xb");
992  matrix6x6(0,5) = parseDouble(node_xc->getValue(), "xc");
993  matrix6x6(1,1) = parseDouble(node_yy->getValue(), "yy");
994  matrix6x6(1,2) = parseDouble(node_yz->getValue(), "yz");
995  matrix6x6(1,3) = parseDouble(node_ya->getValue(), "ya");
996  matrix6x6(1,4) = parseDouble(node_yb->getValue(), "yb");
997  matrix6x6(1,5) = parseDouble(node_yc->getValue(), "yc");
998  matrix6x6(2,2) = parseDouble(node_zz->getValue(), "zz");
999  matrix6x6(2,3) = parseDouble(node_za->getValue(), "za");
1000  matrix6x6(2,4) = parseDouble(node_zb->getValue(), "zb");
1001  matrix6x6(2,5) = parseDouble(node_zc->getValue(), "zc");
1002  matrix6x6(3,3) = parseDouble(node_aa->getValue(), "aa");
1003  matrix6x6(3,4) = parseDouble(node_ab->getValue(), "ab");
1004  matrix6x6(3,5) = parseDouble(node_ac->getValue(), "ac");
1005  matrix6x6(4,4) = parseDouble(node_bb->getValue(), "bb");
1006  matrix6x6(4,5) = parseDouble(node_bc->getValue(), "bc");
1007  matrix6x6(5,5) = parseDouble(node_cc->getValue(), "cc");
1008 
1009  for (std::map<Alignable*, bool>::const_iterator aliiter = aliset.begin(); aliiter != aliset.end(); ++aliiter) {
1010  Alignable *ali = aliiter->first;
1011  ali->setSurvey(new SurveyDet(ali->surface(), matrix6x6));
1012  }
1013 }
double parseDouble(const XMLCh *str, const char *attribute) const
#define NULL
Definition: scimark2.h:8
tuple node
Definition: Node.py:50
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:126
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:268
math::Error< 6 >::type ErrorMatrix
Definition: Definitions.h:39
void MuonAlignmentInputXML::fillAliToIdeal ( std::map< Alignable *, Alignable * > &  alitoideal,
const align::Alignables alignables,
const align::Alignables ideals 
) const
private

Definition at line 228 of file MuonAlignmentInputXML.cc.

References edm::hlt::Exception.

Referenced by newAlignableMuon().

228  {
229  align::Alignables::const_iterator alignable = alignables.begin();
230  align::Alignables::const_iterator ideal = ideals.begin();
231 
232  while (alignable != alignables.end() && ideal != ideals.end()) {
233  alitoideal[*alignable] = *ideal;
234 
235  fillAliToIdeal(alitoideal, (*alignable)->components(), (*ideal)->components());
236 
237  ++alignable;
238  ++ideal;
239  }
240 
241  if (alignable != alignables.end() || ideal != ideals.end()) {
242  throw cms::Exception("Alignment") << "alignable and ideal-alignable trees are out of sync (this should never happen)";
243  }
244 }
void fillAliToIdeal(std::map< Alignable *, Alignable * > &alitoideal, const align::Alignables &alignables, const align::Alignables &ideals) const
Alignable * MuonAlignmentInputXML::getCSCnode ( align::StructureType  structureType,
std::map< unsigned int, Alignable * > &  alignableNavigator,
const XERCES_CPP_NAMESPACE::DOMElement *  node 
) const
private

Definition at line 550 of file MuonAlignmentInputXML.cc.

References align::AlignableCSCChamber, align::AlignableCSCEndcap, align::AlignableCSCRing, align::AlignableCSCStation, align::AlignableDetUnit, Alignable::alignableObjectId(), Reference_intrackfit_cff::endcap, edm::hlt::Exception, Alignable::mother(), NULL, DetId::rawId(), relativeConstraints::ring, relativeConstraints::station, str_chamber, str_endcap, str_layer, str_rawId, str_ring, and str_station.

Referenced by getNode().

550  {
551  unsigned int rawId;
552 
553  DOMAttr *node_rawId = node->getAttributeNode(str_rawId);
554  if (node_rawId != NULL) {
555  try {
556  rawId = XMLString::parseInt(node_rawId->getValue());
557  }
558  catch (const XMLException &toCatch) {
559  throw cms::Exception("XMLException") << "Value of \"rawId\" must be an integer" << std::endl;
560  }
561  }
562  else {
563  int endcap, station, ring, chamber, layer;
564  endcap = station = ring = chamber = layer = 1;
565 
566  DOMAttr *node_endcap = node->getAttributeNode(str_endcap);
567  if (node_endcap == NULL) throw cms::Exception("XMLException") << "CSC node is missing required \"endcap\" attribute" << std::endl;
568  try {
569  endcap = XMLString::parseInt(node_endcap->getValue());
570  }
571  catch (const XMLException &toCatch) {
572  throw cms::Exception("XMLException") << "Value of \"endcap\" must be an integer" << std::endl;
573  }
574  if (endcap == -1) endcap = 2;
575 
576  if (structureType != align::AlignableCSCEndcap) {
577  DOMAttr *node_station = node->getAttributeNode(str_station);
578  if (node_station == NULL) throw cms::Exception("XMLException") << "CSC node is missing required \"station\" attribute" << std::endl;
579  try {
580  station = XMLString::parseInt(node_station->getValue());
581  }
582  catch (const XMLException &toCatch) {
583  throw cms::Exception("XMLException") << "Value of \"station\" must be an integer" << std::endl;
584  }
585 
586  if (structureType != align::AlignableCSCStation) {
587  DOMAttr *node_ring = node->getAttributeNode(str_ring);
588  if (node_ring == NULL) throw cms::Exception("XMLException") << "CSC node is missing required \"ring\" attribute" << std::endl;
589  try {
590  ring = XMLString::parseInt(node_ring->getValue());
591  }
592  catch (const XMLException &toCatch) {
593  throw cms::Exception("XMLException") << "Value of \"ring\" must be an integer" << std::endl;
594  }
595 
596  if (structureType != align::AlignableCSCRing) {
597  DOMAttr *node_chamber = node->getAttributeNode(str_chamber);
598  if (node_chamber == NULL) throw cms::Exception("XMLException") << "CSC node is missing required \"chamber\" attribute" << std::endl;
599  try {
600  chamber = XMLString::parseInt(node_chamber->getValue());
601  }
602  catch (const XMLException &toCatch) {
603  throw cms::Exception("XMLException") << "Value of \"chamber\" must be an integer" << std::endl;
604  }
605 
606  if (structureType != align::AlignableCSCChamber) {
607  DOMAttr *node_layer = node->getAttributeNode(str_layer);
608  if (node_layer == NULL) throw cms::Exception("XMLException") << "CSC node is missing required \"layer\" attribute" << std::endl;
609  try {
610  layer = XMLString::parseInt(node_layer->getValue());
611  }
612  catch (const XMLException &toCatch) {
613  throw cms::Exception("XMLException") << "Value of \"layer\" must be an integer" << std::endl;
614  }
615 
616  } // end if we need a layer number
617  } // end if we need a chamber number
618  } // end if we need a ring number
619  } // end if we need a station number
620 
621  CSCDetId layerId(endcap, station, ring, chamber, layer);
622  rawId = layerId.rawId();
623  } // end if it's specified by endcap, station, ring, chamber, layer
624 
625  Alignable *ali = alignableNavigator[rawId];
626  if (ali == NULL) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not recognized" << std::endl;
627 
628  while (ali->alignableObjectId() != structureType) {
629  ali = ali->mother();
630 
631  if (ali == NULL) {
632  if (structureType == align::AlignableCSCEndcap) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a CSCEndcap" << std::endl;
633  else if (structureType == align::AlignableCSCStation) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a CSCStation" << std::endl;
634  else if (structureType == align::AlignableCSCRing) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a CSCRing" << std::endl;
635  else if (structureType == align::AlignableCSCChamber) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a CSCChamber" << std::endl;
636  else if (structureType == align::AlignableDetUnit) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a CSCLayer" << std::endl;
637  }
638  }
639  return ali;
640 }
#define NULL
Definition: scimark2.h:8
tuple node
Definition: Node.py:50
virtual StructureType alignableObjectId() const =0
Return the alignable type identifier.
Alignable * mother() const
Return pointer to container alignable (if any)
Definition: Alignable.h:85
Alignable * MuonAlignmentInputXML::getDTnode ( align::StructureType  structureType,
std::map< unsigned int, Alignable * > &  alignableNavigator,
const XERCES_CPP_NAMESPACE::DOMElement *  node 
) const
private

Definition at line 456 of file MuonAlignmentInputXML.cc.

References align::AlignableDetUnit, align::AlignableDTBarrel, align::AlignableDTChamber, align::AlignableDTStation, align::AlignableDTSuperLayer, align::AlignableDTWheel, Alignable::alignableObjectId(), edm::hlt::Exception, Alignable::mother(), NULL, DetId::rawId(), relativeConstraints::station, str_layer, str_rawId, str_sector, str_station, str_superlayer, and str_wheel.

Referenced by getNode().

456  {
457  unsigned int rawId = 0;
458 
459  DOMAttr *node_rawId = node->getAttributeNode(str_rawId);
460  if (node_rawId != NULL) {
461  try {
462  rawId = XMLString::parseInt(node_rawId->getValue());
463  }
464  catch (const XMLException &toCatch) {
465  throw cms::Exception("XMLException") << "Value of \"rawId\" must be an integer" << std::endl;
466  }
467  }
468  else {
469  int wheel, station, sector, superlayer, layer;
470  wheel = station = sector = superlayer = layer = 1;
471 
472  if (structureType != align::AlignableDTBarrel) {
473  DOMAttr *node_wheel = node->getAttributeNode(str_wheel);
474  if (node_wheel == NULL) throw cms::Exception("XMLException") << "DT node is missing required \"wheel\" attribute" << std::endl;
475  try {
476  wheel = XMLString::parseInt(node_wheel->getValue());
477  }
478  catch (const XMLException &toCatch) {
479  throw cms::Exception("XMLException") << "Value of \"wheel\" must be an integer" << std::endl;
480  }
481 
482  if (structureType != align::AlignableDTWheel) {
483  DOMAttr *node_station = node->getAttributeNode(str_station);
484  if (node_station == NULL) throw cms::Exception("XMLException") << "DT node is missing required \"station\" attribute" << std::endl;
485  try {
486  station = XMLString::parseInt(node_station->getValue());
487  }
488  catch (const XMLException &toCatch) {
489  throw cms::Exception("XMLException") << "Value of \"station\" must be an integer" << std::endl;
490  }
491 
492  if (structureType != align::AlignableDTStation) {
493  DOMAttr *node_sector = node->getAttributeNode(str_sector);
494  if (node_sector == NULL) throw cms::Exception("XMLException") << "DT node is missing required \"sector\" attribute" << std::endl;
495  try {
496  sector = XMLString::parseInt(node_sector->getValue());
497  }
498  catch (const XMLException &toCatch) {
499  throw cms::Exception("XMLException") << "Value of \"sector\" must be an integer" << std::endl;
500  }
501 
502  if (structureType != align::AlignableDTChamber) {
503  DOMAttr *node_superlayer = node->getAttributeNode(str_superlayer);
504  if (node_superlayer == NULL) throw cms::Exception("XMLException") << "DT node is missing required \"superlayer\" attribute" << std::endl;
505  try {
506  superlayer = XMLString::parseInt(node_superlayer->getValue());
507  }
508  catch (const XMLException &toCatch) {
509  throw cms::Exception("XMLException") << "Value of \"superlayer\" must be an integer" << std::endl;
510  }
511 
512  if (structureType != align::AlignableDTSuperLayer) {
513  DOMAttr *node_layer = node->getAttributeNode(str_layer);
514  if (node_layer == NULL) throw cms::Exception("XMLException") << "DT node is missing required \"layer\" attribute" << std::endl;
515  try {
516  layer = XMLString::parseInt(node_layer->getValue());
517  }
518  catch (const XMLException &toCatch) {
519  throw cms::Exception("XMLException") << "Value of \"layer\" must be an integer" << std::endl;
520  }
521 
522  } // end if we need a layer number
523  } // end if we need a superlayer number
524  } // end if we need a sector number
525  } // end if we need a station number
526  } // end if we need a wheel number
527 
528  DTLayerId layerId(wheel, station, sector, superlayer, layer);
529  rawId = layerId.rawId();
530  } // end if it's specified by wheel, station, sector, superlayer, layer
531 
532  Alignable *ali = alignableNavigator[rawId];
533  if (ali == NULL) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not recognized" << std::endl;
534 
535  while (ali->alignableObjectId() != structureType) {
536  ali = ali->mother();
537 
538  if (ali == NULL) {
539  if (structureType == align::AlignableDTBarrel) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a DTBarrel" << std::endl;
540  else if (structureType == align::AlignableDTWheel) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a DTWheel" << std::endl;
541  else if (structureType == align::AlignableDTStation) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a DTStation" << std::endl;
542  else if (structureType == align::AlignableDTChamber) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a DTChamber" << std::endl;
543  else if (structureType == align::AlignableDTSuperLayer) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a DTSuperLayer" << std::endl;
544  else if (structureType == align::AlignableDetUnit) throw cms::Exception("XMLException") << "rawId \"" << rawId << "\" is not a DTLayer" << std::endl;
545  }
546  }
547  return ali;
548 }
#define NULL
Definition: scimark2.h:8
tuple node
Definition: Node.py:50
virtual StructureType alignableObjectId() const =0
Return the alignable type identifier.
Alignable * mother() const
Return pointer to container alignable (if any)
Definition: Alignable.h:85
Alignable * MuonAlignmentInputXML::getNode ( std::map< unsigned int, Alignable * > &  alignableNavigator,
const XERCES_CPP_NAMESPACE::DOMElement *  node 
) const
private

Definition at line 441 of file MuonAlignmentInputXML.cc.

References align::AlignableCSCChamber, align::AlignableCSCEndcap, align::AlignableCSCRing, align::AlignableCSCStation, align::AlignableDetUnit, align::AlignableDTBarrel, align::AlignableDTChamber, align::AlignableDTStation, align::AlignableDTSuperLayer, align::AlignableDTWheel, getCSCnode(), getDTnode(), NULL, str_CSCChamber, str_CSCEndcap, str_CSCLayer, str_CSCRing, str_CSCStation, str_DTBarrel, str_DTChamber, str_DTLayer, str_DTStation, str_DTSuperLayer, and str_DTWheel.

Referenced by python.TagTree.tagTree::deleteNode(), python.TagTree.tagTree::deleteSubtree(), python.TagTree.tagTree::getPath(), python.TagTree.tagTree::getSubtree(), python.TagTree.tagTree::insertNode(), python.TagTree.tagTree::nChildren(), and newAlignableMuon().

441  {
442  if (XMLString::equals(node->getNodeName(), str_DTBarrel)) return getDTnode(align::AlignableDTBarrel, alignableNavigator, node);
443  else if (XMLString::equals(node->getNodeName(), str_DTWheel)) return getDTnode(align::AlignableDTWheel, alignableNavigator, node);
444  else if (XMLString::equals(node->getNodeName(), str_DTStation)) return getDTnode(align::AlignableDTStation, alignableNavigator, node);
445  else if (XMLString::equals(node->getNodeName(), str_DTChamber)) return getDTnode(align::AlignableDTChamber, alignableNavigator, node);
446  else if (XMLString::equals(node->getNodeName(), str_DTSuperLayer)) return getDTnode(align::AlignableDTSuperLayer, alignableNavigator, node);
447  else if (XMLString::equals(node->getNodeName(), str_DTLayer)) return getDTnode(align::AlignableDetUnit, alignableNavigator, node);
448  else if (XMLString::equals(node->getNodeName(), str_CSCEndcap)) return getCSCnode(align::AlignableCSCEndcap, alignableNavigator, node);
449  else if (XMLString::equals(node->getNodeName(), str_CSCStation)) return getCSCnode(align::AlignableCSCStation, alignableNavigator, node);
450  else if (XMLString::equals(node->getNodeName(), str_CSCRing)) return getCSCnode(align::AlignableCSCRing, alignableNavigator, node);
451  else if (XMLString::equals(node->getNodeName(), str_CSCChamber)) return getCSCnode(align::AlignableCSCChamber, alignableNavigator, node);
452  else if (XMLString::equals(node->getNodeName(), str_CSCLayer)) return getCSCnode(align::AlignableDetUnit, alignableNavigator, node);
453  else return NULL;
454 }
#define NULL
Definition: scimark2.h:8
tuple node
Definition: Node.py:50
Alignable * getDTnode(align::StructureType structureType, std::map< unsigned int, Alignable * > &alignableNavigator, const XERCES_CPP_NAMESPACE::DOMElement *node) const
Alignable * getCSCnode(align::StructureType structureType, std::map< unsigned int, Alignable * > &alignableNavigator, const XERCES_CPP_NAMESPACE::DOMElement *node) const
AlignableMuon * MuonAlignmentInputXML::newAlignableMuon ( const edm::EventSetup iSetup) const
virtual

Reimplemented from MuonAlignmentInputMethod.

Definition at line 246 of file MuonAlignmentInputXML.cc.

References align::AlignableMuon, runEdmFileComparison::collection, AlignableMuon::CSCEndcaps(), do_moveglobal(), do_movelocal(), do_rotatebeamline(), do_rotateglobalaxis(), do_rotatelocal(), do_setape(), do_setposition(), do_setsurveyerr(), asciidump::doc, AlignableMuon::DTBarrel(), edm::hlt::Exception, fillAliToIdeal(), getNode(), i, MuonAlignmentInputMethod::idealCSCGeometry(), MuonAlignmentInputMethod::idealDTGeometry(), j, m_fileName, python.rootplot.argparse::message, mergeVDriftHistosByStation::name, python.Node::node, NULL, geometryXMLtoCSV::parser, recursiveGetId(), fetchall_from_DQM_v2::release, str_collection, str_moveglobal, str_movelocal, str_name, str_operation, str_rotatebeamline, str_rotateglobalaxis, str_rotatelocal, str_setape, str_setposition, str_setsurveyerr, AlCaHLTBitMon_QueryRunRegistry::string, cms::concurrency::xercesInitialize(), and cms::concurrency::xercesTerminate().

246  {
247  boost::shared_ptr<DTGeometry> dtGeometry = idealDTGeometry(iSetup);
248  boost::shared_ptr<CSCGeometry> cscGeometry = idealCSCGeometry(iSetup);
249 
250  AlignableMuon *alignableMuon = new AlignableMuon(&(*dtGeometry), &(*cscGeometry));
251  std::map<unsigned int, Alignable*> alignableNavigator; // real AlignableNavigators don't have const methods
252  recursiveGetId(alignableNavigator, alignableMuon->DTBarrel());
253  recursiveGetId(alignableNavigator, alignableMuon->CSCEndcaps());
254 
255  AlignableMuon *ideal_alignableMuon = new AlignableMuon(&(*dtGeometry), &(*cscGeometry));
256  std::map<unsigned int, Alignable*> ideal_alignableNavigator; // real AlignableNavigators don't have const methods
257  recursiveGetId(ideal_alignableNavigator, ideal_alignableMuon->DTBarrel());
258  recursiveGetId(ideal_alignableNavigator, ideal_alignableMuon->CSCEndcaps());
259 
260  try {
262  }
263  catch (const XMLException &toCatch) {
264  throw cms::Exception("XMLException") << "Xerces XML parser threw an exception on initialization." << std::endl;
265  }
266 
267  XercesDOMParser *parser = new XercesDOMParser();
268  parser->setValidationScheme(XercesDOMParser::Val_Always);
269 
270  ErrorHandler *errHandler = (ErrorHandler*)(new HandlerBase());
271  parser->setErrorHandler(errHandler);
272 
273  try {
274  parser->parse(m_fileName.c_str());
275  }
276  catch (const XMLException &toCatch) {
277  char *message = XMLString::transcode(toCatch.getMessage());
278  throw cms::Exception("XMLException") << "Xerces XML parser threw this exception: " << message << std::endl;
279  }
280  catch (const DOMException &toCatch) {
281  char *message = XMLString::transcode(toCatch.msg);
282  throw cms::Exception("XMLException") << "Xerces XML parser threw this exception: " << message << std::endl;
283  }
284  catch (const SAXException &toCatch) {
285  char *message = XMLString::transcode(toCatch.getMessage());
286  throw cms::Exception("XMLException") << "Xerces XML parser threw this exception: " << message << std::endl;
287  }
288  catch (...) {
289  throw cms::Exception("XMLException") << "Xerces XML parser threw an unknown exception" << std::endl;
290  }
291 
292  DOMDocument *doc = parser->getDocument();
293  DOMElement *node_MuonAlignment = doc->getDocumentElement();
294  DOMNodeList *collections = doc->getElementsByTagName(str_collection);
295  DOMNodeList *operations = doc->getElementsByTagName(str_operation);
296 
297  std::map<Alignable*, Alignable*> alitoideal;
298  fillAliToIdeal(alitoideal, alignableMuon->DTBarrel(), ideal_alignableMuon->DTBarrel());
299  fillAliToIdeal(alitoideal, alignableMuon->CSCEndcaps(), ideal_alignableMuon->CSCEndcaps());
300 
301  std::map<std::string, std::map<Alignable*, bool> > alicollections;
302  for (unsigned int i = 0; i < collections->getLength(); i++) {
303  DOMElement *collection = (DOMElement*)(collections->item(i));
304  if (collection->getParentNode() == node_MuonAlignment) {
305  DOMNodeList *children = collection->getChildNodes();
306 
307  DOMAttr *node_name = collection->getAttributeNode(str_name);
308  if (node_name == NULL) {
309  throw cms::Exception("XMLException") << "<collection> requires a name attribute" << std::endl;
310  }
311  char *ascii_name = XMLString::transcode(node_name->getValue());
312  std::string name(ascii_name);
313  XMLString::release(&ascii_name);
314 
315  std::map<Alignable*, bool> aliset;
316  for (unsigned int j = 0; j < children->getLength(); j++) {
317  DOMNode *node = children->item(j);
318 
319  if (node->getNodeType() == DOMNode::ELEMENT_NODE) {
320  Alignable *ali = getNode(alignableNavigator, (DOMElement*)(node));
321  if (ali == NULL) {
322  throw cms::Exception("XMLException") << "<collection> must contain only alignables" << std::endl;
323  }
324 
325  aliset[ali] = true;
326  } // end if this node is an element
327  } // end loop over collection's children
328 
329  alicollections[name] = aliset;
330  } // end if this is a top-level collection
331  } // end loop over collections
332 
333  for (unsigned int i = 0; i < operations->getLength(); i++) {
334  DOMElement *operation = (DOMElement*)(operations->item(i));
335  if (operation->getParentNode() != node_MuonAlignment) {
336  throw cms::Exception("XMLException") << "All operations must be top-level elements" << std::endl;
337  }
338 
339  DOMNodeList *children = operation->getChildNodes();
340 
341  std::map<Alignable*, bool> aliset;
342  std::vector<DOMNode*> nodesToRemove;
343  for (unsigned int j = 0; j < children->getLength(); j++) {
344  DOMNode *node = children->item(j);
345 
346  if (node->getNodeType() == DOMNode::ELEMENT_NODE) {
347  Alignable *ali = getNode(alignableNavigator, (DOMElement*)(node));
348  if (ali != NULL) {
349  aliset[ali] = true;
350  nodesToRemove.push_back(node);
351  } // end if this node is an alignable
352 
353  else if (XMLString::equals(node->getNodeName(), str_collection)) {
354  DOMAttr *node_name = ((DOMElement*)(node))->getAttributeNode(str_name);
355  if (node_name == NULL) {
356  throw cms::Exception("XMLException") << "<collection> requires a name attribute" << std::endl;
357  }
358  char *ascii_name = XMLString::transcode(node_name->getValue());
359  std::string name(ascii_name);
360  XMLString::release(&ascii_name);
361 
362  std::map<std::string, std::map<Alignable*, bool> >::const_iterator alicollections_iter = alicollections.find(name);
363  if (alicollections_iter == alicollections.end()) {
364  throw cms::Exception("XMLException") << "<collection name=\"" << name << "\"> hasn't been defined" << std::endl;
365  }
366 
367  for (std::map<Alignable*, bool>::const_iterator aliiter = alicollections_iter->second.begin();
368  aliiter != alicollections_iter->second.end();
369  ++aliiter) {
370  aliset[aliiter->first] = true;
371  } // end loop over alignables in this collection
372 
373  nodesToRemove.push_back(node);
374  } // end if this node is a collection
375 
376  else {} // anything else? assume it's a position/rotation directive
377 
378  } // end if node is node is an element
379  } // end first loop over operation's children
380 
381  // from now on, we only want to see position/rotation directives
382  for (std::vector<DOMNode*>::const_iterator node = nodesToRemove.begin(); node != nodesToRemove.end(); ++node) {
383  operation->removeChild(*node);
384  }
385  children = operation->getChildNodes();
386 
387  for (unsigned int j = 0; j < children->getLength(); j++) {
388  DOMNode *node = children->item(j);
389  if (node->getNodeType() == DOMNode::ELEMENT_NODE) {
390 
391  if (XMLString::equals(node->getNodeName(), str_setposition)) {
392  do_setposition((DOMElement*)(node), aliset, alitoideal);
393  }
394 
395  else if (XMLString::equals(node->getNodeName(), str_setape)) {
396  do_setape((DOMElement*)(node), aliset, alitoideal);
397  }
398 
399  else if (XMLString::equals(node->getNodeName(), str_setsurveyerr)) {
400  do_setsurveyerr((DOMElement*)(node), aliset, alitoideal);
401  }
402 
403  else if (XMLString::equals(node->getNodeName(), str_moveglobal)) {
404  do_moveglobal((DOMElement*)(node), aliset, alitoideal);
405  }
406 
407  else if (XMLString::equals(node->getNodeName(), str_movelocal)) {
408  do_movelocal((DOMElement*)(node), aliset, alitoideal);
409  }
410 
411  else if (XMLString::equals(node->getNodeName(), str_rotatelocal)) {
412  do_rotatelocal((DOMElement*)(node), aliset, alitoideal);
413  }
414 
415  else if (XMLString::equals(node->getNodeName(), str_rotatebeamline)) {
416  do_rotatebeamline((DOMElement*)(node), aliset, alitoideal);
417  }
418 
419  else if (XMLString::equals(node->getNodeName(), str_rotateglobalaxis)) {
420  do_rotateglobalaxis((DOMElement*)(node), aliset, alitoideal);
421  }
422 
423  else {
424  char *message = XMLString::transcode(node->getNodeName());
425  throw cms::Exception("XMLException") << "Unrecognized operation: \"" << message << "\"" << std::endl;
426  }
427 
428  } // end if node is an element
429  } // end second loop over operation's children
430  } // end loop over operations
431 
432  delete parser;
433  delete errHandler;
434 
436 
437  delete ideal_alignableMuon;
438  return alignableMuon;
439 }
void do_setsurveyerr(const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
int i
Definition: DBlmapReader.cc:9
void do_rotatelocal(const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
void xercesTerminate()
Definition: Xerces.cc:22
align::Alignables DTBarrel()
#define NULL
Definition: scimark2.h:8
void do_movelocal(const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
void do_moveglobal(const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
void xercesInitialize()
Definition: Xerces.cc:17
tuple node
Definition: Node.py:50
void do_rotatebeamline(const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
align::Alignables CSCEndcaps()
void recursiveGetId(std::map< unsigned int, Alignable * > &alignableNavigator, const align::Alignables &alignables) const
Alignable * getNode(std::map< unsigned int, Alignable * > &alignableNavigator, const XERCES_CPP_NAMESPACE::DOMElement *node) const
int j
Definition: DBlmapReader.cc:9
void do_setape(const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
tuple doc
Definition: asciidump.py:381
boost::shared_ptr< CSCGeometry > idealCSCGeometry(const edm::EventSetup &iSetup) const
void do_rotateglobalaxis(const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
void fillAliToIdeal(std::map< Alignable *, Alignable * > &alitoideal, const align::Alignables &alignables, const align::Alignables &ideals) const
void do_setposition(const XERCES_CPP_NAMESPACE::DOMElement *node, std::map< Alignable *, bool > &aliset, std::map< Alignable *, Alignable * > &alitoideal) const
Constructor of the full muon geometry.
Definition: AlignableMuon.h:36
boost::shared_ptr< DTGeometry > idealDTGeometry(const edm::EventSetup &iSetup) const
const MuonAlignmentInputXML& MuonAlignmentInputXML::operator= ( const MuonAlignmentInputXML )
private
double MuonAlignmentInputXML::parseDouble ( const XMLCh *  str,
const char *  attribute 
) const
private

Definition at line 642 of file MuonAlignmentInputXML.cc.

References edm::hlt::Exception, i, reco::HaloData::minus, convertSQLitetoXML_cfg::output, funct::pow(), and fetchall_from_DQM_v2::release.

Referenced by do_moveglobal(), do_movelocal(), do_rotatebeamline(), do_rotateglobalaxis(), do_rotatelocal(), do_setape(), do_setposition(), and do_setsurveyerr().

642  {
643  unsigned int len = XMLString::stringLen(str);
644  char *cstr = XMLString::transcode(str);
645  std::stringstream errmessage;
646  errmessage << "Value of \"" << attribute << "\" must be a double, not \"" << cstr << "\"" << std::endl;
647 
648  unsigned int i = 0;
649 
650  bool minus = false;
651  if (cstr[i] == '-') {
652  minus = true;
653  i++;
654  }
655  else if (cstr[i] == '+') i++;
656 
657  double output = 0.;
658 
659  while (cstr[i] != '.' && cstr[i] != 'e' && cstr[i] != 'E' && i < len) {
660  if (cstr[i] < '0' || cstr[i] > '9') {
661  XMLString::release(&cstr);
662  throw cms::Exception("XMLException") << errmessage.str();
663  }
664 
665  output *= 10;
666  output += cstr[i] - '0';
667  i++;
668  }
669 
670  if (cstr[i] == '.') {
671  double place = 0.1;
672  i++;
673 
674  while (cstr[i] != 'e' && cstr[i] != 'E' && i < len) {
675  if (cstr[i] < '0' || cstr[i] > '9') {
676  XMLString::release(&cstr);
677  throw cms::Exception("XMLException") << errmessage.str();
678  }
679 
680  output += (cstr[i] - '0') * place;
681  place /= 10.;
682  i++;
683  }
684  }
685 
686  if (cstr[i] == 'e' || cstr[i] == 'E') {
687  i++;
688 
689  int exponent = 0;
690  bool expminus = false;
691  if (cstr[i] == '-') {
692  expminus = true;
693  i++;
694  }
695  else if (cstr[i] == '+') i++;
696 
697  while (i < len) {
698  if (cstr[i] < '0' || cstr[i] > '9') {
699  XMLString::release(&cstr);
700  throw cms::Exception("XMLException") << errmessage.str();
701  }
702 
703  exponent *= 10;
704  exponent += cstr[i] - '0';
705  i++;
706  }
707 
708  if (expminus) exponent *= -1;
709 
710  output *= pow(10., exponent);
711  }
712 
713  if (minus) output *= -1.;
714 
715  XMLString::release(&cstr);
716  return output;
717 }
int i
Definition: DBlmapReader.cc:9
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
void MuonAlignmentInputXML::recursiveGetId ( std::map< unsigned int, Alignable * > &  alignableNavigator,
const align::Alignables alignables 
) const
private

Definition at line 217 of file MuonAlignmentInputXML.cc.

References align::AlignableCSCChamber, align::AlignableCSCLayer, align::AlignableDet, align::AlignableDetUnit, align::AlignableDTChamber, align::AlignableDTLayer, and align::AlignableDTSuperLayer.

Referenced by newAlignableMuon().

217  {
218  for (align::Alignables::const_iterator ali = alignables.begin(); ali != alignables.end(); ++ali) {
219  if ((*ali)->alignableObjectId() == align::AlignableDetUnit || (*ali)->alignableObjectId() == align::AlignableDet ||
220  (*ali)->alignableObjectId() == align::AlignableDTChamber || (*ali)->alignableObjectId() == align::AlignableDTSuperLayer || (*ali)->alignableObjectId() == align::AlignableDTLayer ||
221  (*ali)->alignableObjectId() == align::AlignableCSCChamber || (*ali)->alignableObjectId() == align::AlignableCSCLayer) {
222  alignableNavigator[(*ali)->geomDetId().rawId()] = *ali;
223  }
224  recursiveGetId(alignableNavigator, (*ali)->components());
225  }
226 }
void recursiveGetId(std::map< unsigned int, Alignable * > &alignableNavigator, const align::Alignables &alignables) const
void MuonAlignmentInputXML::set_one_position ( Alignable ali,
const align::PositionType pos,
const align::RotationType rot 
) const
private

Definition at line 872 of file MuonAlignmentInputXML.cc.

References SurveyDet::errors(), Alignable::globalPosition(), Alignable::globalRotation(), Alignable::move(), TkRotation< T >::multiplyInverse(), NULL, align::rectify(), Alignable::rotateInGlobalFrame(), Alignable::setSurvey(), Alignable::surface(), and Alignable::survey().

Referenced by do_setposition().

872  {
873  const align::PositionType& oldpos = ali->globalPosition();
874  const align::RotationType& oldrot = ali->globalRotation();
875 
876  // shift needed to move from current to new position
877  align::GlobalVector posDiff = pos - oldpos;
878  align::RotationType rotDiff = oldrot.multiplyInverse(rot);
879  align::rectify(rotDiff); // correct for rounding errors
880  ali->move(posDiff);
881  ali->rotateInGlobalFrame(rotDiff);
882 
883 // // check for consistency
884 // const align::PositionType& newpos = ali->globalPosition();
885 // const align::RotationType& newrot = ali->globalRotation();
886 // align::GlobalVector posDiff2 = pos - newpos;
887 // align::RotationType rotDiff2 = newrot.multiplyInverse(rot);
888 // align::rectify(rotDiff2); // correct for rounding errors
889 
890 // if (fabs(posDiff2.x()) > 1e-6 || fabs(posDiff2.y()) > 1e-6 || fabs(posDiff2.z()) > 1e-6) {
891 // std::cout << "zeropos " << posDiff2 << std::endl;
892 // }
893 // if (fabs(rotDiff2.xx() - 1.) > 1e-4 ||
894 // fabs(rotDiff2.yy() - 1.) > 1e-4 ||
895 // fabs(rotDiff2.zz() - 1.) > 1e-4 ||
896 // fabs(rotDiff2.xy()) > 1e-8 ||
897 // fabs(rotDiff2.xz()) > 1e-8 ||
898 // fabs(rotDiff2.yz()) > 1e-8) {
899 // std::cout << "zerorot " << rotDiff2 << std::endl;
900 // }
901 
902  align::ErrorMatrix matrix6x6 = ROOT::Math::SMatrixIdentity();
903  matrix6x6 *= 1000.; // initial assumption: infinitely weak constraint
904 
905  const SurveyDet *survey = ali->survey();
906  if (survey != NULL) {
907  matrix6x6 = survey->errors(); // save the constraint information
908  }
909  ali->setSurvey(new SurveyDet(ali->surface(), matrix6x6));
910 }
const SurveyDet * survey() const
Return survey info.
Definition: Alignable.h:205
#define NULL
Definition: scimark2.h:8
const RotationType & globalRotation() const
Return the global orientation of the object.
Definition: Alignable.h:132
virtual void move(const GlobalVector &displacement)=0
Movement with respect to the global reference frame.
const align::ErrorMatrix & errors() const
Definition: SurveyDet.h:73
void rectify(RotationType &)
Correct a rotation matrix for rounding errors.
Definition: Utilities.cc:196
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:126
virtual void rotateInGlobalFrame(const RotationType &rotation)=0
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:268
const PositionType & globalPosition() const
Return the global position of the object.
Definition: Alignable.h:129
math::Error< 6 >::type ErrorMatrix
Definition: Definitions.h:39
Basic3DVector< T > multiplyInverse(const Basic3DVector< T > &v) const

Member Data Documentation

std::string MuonAlignmentInputXML::m_fileName
private

Definition at line 73 of file MuonAlignmentInputXML.h.

Referenced by newAlignableMuon().

XMLCh* MuonAlignmentInputXML::str_aa
private
XMLCh* MuonAlignmentInputXML::str_ab
private
XMLCh* MuonAlignmentInputXML::str_ac
private
XMLCh* MuonAlignmentInputXML::str_alpha
private
XMLCh* MuonAlignmentInputXML::str_angle
private
XMLCh* MuonAlignmentInputXML::str_axisx
private
XMLCh* MuonAlignmentInputXML::str_axisy
private
XMLCh* MuonAlignmentInputXML::str_axisz
private
XMLCh* MuonAlignmentInputXML::str_bb
private
XMLCh* MuonAlignmentInputXML::str_bc
private
XMLCh* MuonAlignmentInputXML::str_beta
private
XMLCh* MuonAlignmentInputXML::str_cc
private
XMLCh* MuonAlignmentInputXML::str_chamber
private
XMLCh* MuonAlignmentInputXML::str_collection
private
XMLCh* MuonAlignmentInputXML::str_container
private
XMLCh* MuonAlignmentInputXML::str_CSCChamber
private

Definition at line 87 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_CSCEndcap
private

Definition at line 84 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_CSCLayer
private

Definition at line 88 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_CSCRing
private

Definition at line 86 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_CSCStation
private

Definition at line 85 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_DTBarrel
private

Definition at line 78 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_DTChamber
private

Definition at line 81 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_DTLayer
private

Definition at line 83 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_DTStation
private

Definition at line 80 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_DTSuperLayer
private

Definition at line 82 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_DTWheel
private

Definition at line 79 of file MuonAlignmentInputXML.h.

Referenced by getNode(), MuonAlignmentInputXML(), and ~MuonAlignmentInputXML().

XMLCh* MuonAlignmentInputXML::str_endcap
private
XMLCh* MuonAlignmentInputXML::str_gamma
private
XMLCh* MuonAlignmentInputXML::str_ideal
private
XMLCh* MuonAlignmentInputXML::str_layer
private
XMLCh* MuonAlignmentInputXML::str_moveglobal
private
XMLCh* MuonAlignmentInputXML::str_movelocal
private
XMLCh* MuonAlignmentInputXML::str_name
private
XMLCh* MuonAlignmentInputXML::str_none
private
XMLCh* MuonAlignmentInputXML::str_operation
private
XMLCh* MuonAlignmentInputXML::str_phi
private
XMLCh* MuonAlignmentInputXML::str_phix
private
XMLCh* MuonAlignmentInputXML::str_phiy
private
XMLCh* MuonAlignmentInputXML::str_phiz
private
XMLCh* MuonAlignmentInputXML::str_rawId
private
XMLCh* MuonAlignmentInputXML::str_relativeto
private
XMLCh* MuonAlignmentInputXML::str_ring
private
XMLCh* MuonAlignmentInputXML::str_rotatebeamline
private
XMLCh* MuonAlignmentInputXML::str_rotateglobalaxis
private
XMLCh* MuonAlignmentInputXML::str_rotatelocal
private
XMLCh* MuonAlignmentInputXML::str_rphi
private
XMLCh* MuonAlignmentInputXML::str_sector
private
XMLCh* MuonAlignmentInputXML::str_setape
private
XMLCh* MuonAlignmentInputXML::str_setposition
private
XMLCh* MuonAlignmentInputXML::str_setsurveyerr
private
XMLCh* MuonAlignmentInputXML::str_station
private
XMLCh* MuonAlignmentInputXML::str_superlayer
private
XMLCh* MuonAlignmentInputXML::str_wheel
private
XMLCh* MuonAlignmentInputXML::str_x
private
XMLCh* MuonAlignmentInputXML::str_xa
private
XMLCh* MuonAlignmentInputXML::str_xb
private
XMLCh* MuonAlignmentInputXML::str_xc
private
XMLCh* MuonAlignmentInputXML::str_xx
private
XMLCh* MuonAlignmentInputXML::str_xy
private
XMLCh* MuonAlignmentInputXML::str_xz
private
XMLCh* MuonAlignmentInputXML::str_y
private
XMLCh* MuonAlignmentInputXML::str_ya
private
XMLCh* MuonAlignmentInputXML::str_yb
private
XMLCh* MuonAlignmentInputXML::str_yc
private
XMLCh* MuonAlignmentInputXML::str_yy
private
XMLCh* MuonAlignmentInputXML::str_yz
private
XMLCh* MuonAlignmentInputXML::str_z
private
XMLCh* MuonAlignmentInputXML::str_za
private
XMLCh* MuonAlignmentInputXML::str_zb
private
XMLCh* MuonAlignmentInputXML::str_zc
private
XMLCh* MuonAlignmentInputXML::str_zz
private