|
|
|
#include <Alignment/MuonAlignment/interface/MuonAlignmentOutputXML.h>
Public Member Functions | |
| MuonAlignmentOutputXML (const edm::ParameterSet &iConfig) | |
| void | write (AlignableMuon *alignableMuon, const edm::EventSetup &iSetup) const |
| virtual | ~MuonAlignmentOutputXML () |
Private Member Functions | |
| MuonAlignmentOutputXML (const MuonAlignmentOutputXML &) | |
| const MuonAlignmentOutputXML & | operator= (const MuonAlignmentOutputXML &) |
| void | writeComponents (align::Alignables &alignables, align::Alignables &ideals, std::map< align::ID, CLHEP::HepSymMatrix > &errors, std::ofstream &outputFile, bool DT) const |
Private Attributes | |
| bool | m_eulerAngles |
| std::string | m_fileName |
| int | m_precision |
| bool | m_rawIds |
| int | m_relativeto |
| bool | m_suppressCSCChambers |
| bool | m_suppressCSCEndcaps |
| bool | m_suppressCSCLayers |
| bool | m_suppressCSCRings |
| bool | m_suppressCSCStations |
| bool | m_suppressDTBarrel |
| bool | m_suppressDTChambers |
| bool | m_suppressDTLayers |
| bool | m_suppressDTStations |
| bool | m_suppressDTSuperLayers |
| bool | m_suppressDTWheels |
| bool | m_survey |
Description: <one line="" class="" summary>="">
Usage: <usage>
Definition at line 34 of file MuonAlignmentOutputXML.h.
| MuonAlignmentOutputXML::MuonAlignmentOutputXML | ( | const edm::ParameterSet & | iConfig | ) |
Definition at line 44 of file MuonAlignmentOutputXML.cc.
References Exception, edm::ParameterSet::getParameter(), m_relativeto, and AlCaHLTBitMon_QueryRunRegistry::string.
: m_fileName(iConfig.getParameter<std::string>("fileName")) , m_survey(iConfig.getParameter<bool>("survey")) , m_rawIds(iConfig.getParameter<bool>("rawIds")) , m_eulerAngles(iConfig.getParameter<bool>("eulerAngles")) , m_precision(iConfig.getParameter<int>("precision")) , m_suppressDTBarrel(iConfig.getUntrackedParameter<bool>("suppressDTBarrel", false)) , m_suppressDTWheels(iConfig.getUntrackedParameter<bool>("suppressDTWheels", false)) , m_suppressDTStations(iConfig.getUntrackedParameter<bool>("suppressDTStations", false)) , m_suppressDTChambers(iConfig.getUntrackedParameter<bool>("suppressDTChambers", false)) , m_suppressDTSuperLayers(iConfig.getUntrackedParameter<bool>("suppressDTSuperLayers", false)) , m_suppressDTLayers(iConfig.getUntrackedParameter<bool>("suppressDTLayers", false)) , m_suppressCSCEndcaps(iConfig.getUntrackedParameter<bool>("suppressCSCEndcaps", false)) , m_suppressCSCStations(iConfig.getUntrackedParameter<bool>("suppressCSCStations", false)) , m_suppressCSCRings(iConfig.getUntrackedParameter<bool>("suppressCSCRings", false)) , m_suppressCSCChambers(iConfig.getUntrackedParameter<bool>("suppressCSCChambers", false)) , m_suppressCSCLayers(iConfig.getUntrackedParameter<bool>("suppressCSCLayers", false)) { std::string str_relativeto = iConfig.getParameter<std::string>("relativeto"); if (str_relativeto == std::string("none")) { m_relativeto = 0; } else if (str_relativeto == std::string("ideal")) { m_relativeto = 1; } else if (str_relativeto == std::string("container")) { m_relativeto = 2; } else { throw cms::Exception("BadConfig") << "relativeto must be \"none\", \"ideal\", or \"container\"" << std::endl; } }
| MuonAlignmentOutputXML::~MuonAlignmentOutputXML | ( | ) | [virtual] |
Definition at line 83 of file MuonAlignmentOutputXML.cc.
{
}
| MuonAlignmentOutputXML::MuonAlignmentOutputXML | ( | const MuonAlignmentOutputXML & | ) | [private] |
| const MuonAlignmentOutputXML& MuonAlignmentOutputXML::operator= | ( | const MuonAlignmentOutputXML & | ) | [private] |
| void MuonAlignmentOutputXML::write | ( | AlignableMuon * | alignableMuon, |
| const edm::EventSetup & | iSetup | ||
| ) | const |
Definition at line 103 of file MuonAlignmentOutputXML.cc.
References DTGeometryBuilderFromDDD::build(), CSCGeometryBuilderFromDDD::build(), AlignableMuon::cscAlignmentErrors(), AlignableMuon::CSCEndcaps(), AlignableMuon::dtAlignmentErrors(), AlignableMuon::DTBarrel(), benchmark_cfg::errors, edm::EventSetup::get(), AlignmentErrors::m_alignError, m_fileName, m_precision, m_relativeto, download_sqlite_cfg::outputFile, and writeComponents().
Referenced by MuonAlignment::writeXML().
{
std::ofstream outputFile(m_fileName.c_str());
outputFile << std::setprecision(m_precision) << std::fixed;
outputFile << "<?xml version=\"1.0\" encoding=\"UTF-8\"?>" << std::endl;
outputFile << "<?xml-stylesheet type=\"text/xml\" href=\"MuonAlignment.xsl\"?>" << std::endl;
outputFile << "<MuonAlignment>" << std::endl << std::endl;
std::map<align::ID, CLHEP::HepSymMatrix> errors;
AlignmentErrors *dtErrors = alignableMuon->dtAlignmentErrors();
AlignmentErrors *cscErrors = alignableMuon->cscAlignmentErrors();
for (std::vector<AlignTransformError>::const_iterator dtError = dtErrors->m_alignError.begin(); dtError != dtErrors->m_alignError.end(); ++dtError) {
errors[dtError->rawId()] = dtError->matrix();
}
for (std::vector<AlignTransformError>::const_iterator cscError = cscErrors->m_alignError.begin(); cscError != cscErrors->m_alignError.end(); ++cscError) {
errors[cscError->rawId()] = cscError->matrix();
}
align::Alignables barrels = alignableMuon->DTBarrel();
align::Alignables endcaps = alignableMuon->CSCEndcaps();
if (m_relativeto == 1) {
edm::ESTransientHandle<DDCompactView> cpv;
iSetup.get<IdealGeometryRecord>().get(cpv);
edm::ESHandle<MuonDDDConstants> mdc;
iSetup.get<MuonNumberingRecord>().get(mdc);
DTGeometryBuilderFromDDD DTGeometryBuilder;
CSCGeometryBuilderFromDDD CSCGeometryBuilder;
boost::shared_ptr<DTGeometry> dtGeometry(new DTGeometry );
DTGeometryBuilder.build(dtGeometry, &(*cpv), *mdc);
boost::shared_ptr<CSCGeometry> boost_cscGeometry(new CSCGeometry);
CSCGeometryBuilder.build(boost_cscGeometry, &(*cpv), *mdc);
AlignableMuon ideal_alignableMuon(&(*dtGeometry), &(*boost_cscGeometry));
align::Alignables ideal_barrels = ideal_alignableMuon.DTBarrel();
align::Alignables ideal_endcaps = ideal_alignableMuon.CSCEndcaps();
writeComponents(barrels, ideal_barrels, errors, outputFile, true);
writeComponents(endcaps, ideal_endcaps, errors, outputFile, false);
}
else {
align::Alignables empty1, empty2;
writeComponents(barrels, empty1, errors, outputFile, true);
writeComponents(endcaps, empty2, errors, outputFile, false);
}
outputFile << "</MuonAlignment>" << std::endl;
}
| void MuonAlignmentOutputXML::writeComponents | ( | align::Alignables & | alignables, |
| align::Alignables & | ideals, | ||
| std::map< align::ID, CLHEP::HepSymMatrix > & | errors, | ||
| std::ofstream & | outputFile, | ||
| bool | DT | ||
| ) | const [private] |
Definition at line 157 of file MuonAlignmentOutputXML.cc.
References align::AlignableCSCChamber, align::AlignableCSCEndcap, align::AlignableCSCRing, align::AlignableCSCStation, align::AlignableDetUnit, align::AlignableDTBarrel, align::AlignableDTChamber, align::AlignableDTStation, align::AlignableDTSuperLayer, align::AlignableDTWheel, PV3DBase< T, PVType, FrameType >::basicVector(), makeMuonMisalignmentScenario::components, gather_cfg::cout, Reference_intrackfit_cff::endcap, CSCDetId::endcap(), Exception, AlignableObjectId::idToString(), m_eulerAngles, m_rawIds, m_relativeto, m_suppressCSCChambers, m_suppressCSCEndcaps, m_suppressCSCLayers, m_suppressCSCRings, m_suppressCSCStations, m_suppressDTBarrel, m_suppressDTChambers, m_suppressDTLayers, m_suppressDTStations, m_suppressDTSuperLayers, m_suppressDTWheels, m_survey, NULL, pos, relativeConstraints::ring, makeMuonMisalignmentScenario::rot, relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, align::toAngles(), TkRotation< T >::transposed(), DTChamberId::wheel(), PV3DBase< T, PVType, FrameType >::x(), TkRotation< T >::xx(), TkRotation< T >::xy(), TkRotation< T >::xz(), PV3DBase< T, PVType, FrameType >::y(), TkRotation< T >::yz(), PV3DBase< T, PVType, FrameType >::z(), and TkRotation< T >::zz().
Referenced by write().
{
align::Alignables::const_iterator ideal = ideals.begin();
for (align::Alignables::const_iterator alignable = alignables.begin(); alignable != alignables.end(); ++alignable) {
if (m_survey && (*alignable)->survey() == NULL) {
throw cms::Exception("Alignment") << "SurveyDets must all be defined when writing to XML" << std::endl;
} // now I can assume it's okay everywhere
align::StructureType alignableObjectId = (*alignable)->alignableObjectId();
if ((alignableObjectId == align::AlignableDTBarrel && !m_suppressDTBarrel) ||
(alignableObjectId == align::AlignableDTWheel && !m_suppressDTWheels) ||
(alignableObjectId == align::AlignableDTStation && !m_suppressDTStations) ||
(alignableObjectId == align::AlignableDTChamber && !m_suppressDTChambers) ||
(DT && alignableObjectId == align::AlignableDTSuperLayer && !m_suppressDTSuperLayers) ||
(DT && alignableObjectId == align::AlignableDetUnit && !m_suppressDTLayers) ||
(alignableObjectId == align::AlignableCSCEndcap && !m_suppressCSCEndcaps) ||
(alignableObjectId == align::AlignableCSCStation && !m_suppressCSCStations) ||
(alignableObjectId == align::AlignableCSCRing && !m_suppressCSCRings) ||
(alignableObjectId == align::AlignableCSCChamber && !m_suppressCSCChambers) ||
(!DT && alignableObjectId == align::AlignableDetUnit && !m_suppressCSCLayers)) {
unsigned int rawId = (*alignable)->geomDetId().rawId();
outputFile << "<operation>" << std::endl;
if (DT) {
if (m_rawIds && rawId != 0) {
std::string typeName = AlignableObjectId::idToString(alignableObjectId);
if (alignableObjectId == align::AlignableDTSuperLayer) typeName = std::string("DTSuperLayer");
if (alignableObjectId == align::AlignableDetUnit) typeName = std::string("DTLayer");
outputFile << " <" << typeName << " rawId=\"" << rawId << "\" />" << std::endl;
}
else {
if (alignableObjectId == align::AlignableDetUnit) {
DTLayerId id(rawId);
outputFile << " <DTLayer wheel=\"" << id.wheel() << "\" station=\"" << id.station() << "\" sector=\"" << id.sector() << "\" superlayer=\"" << id.superlayer() << "\" layer=\"" << id.layer() << "\" />" << std::endl;
}
else if (alignableObjectId == align::AlignableDTSuperLayer) {
DTSuperLayerId id(rawId);
outputFile << " <DTSuperLayer wheel=\"" << id.wheel() << "\" station=\"" << id.station() << "\" sector=\"" << id.sector() << "\" superlayer=\"" << id.superlayer() << "\" />" << std::endl;
}
else if (alignableObjectId == align::AlignableDTChamber) {
DTChamberId id(rawId);
outputFile << " <DTChamber wheel=\"" << id.wheel() << "\" station=\"" << id.station() << "\" sector=\"" << id.sector() << "\" />" << std::endl;
}
else {
DTChamberId id((*alignable)->id());
if (alignableObjectId == align::AlignableDTStation) {
outputFile << " <DTStation wheel=\"" << id.wheel() << "\" station=\"" << id.station() << "\" />" << std::endl;
}
else if (alignableObjectId == align::AlignableDTWheel) {
outputFile << " <DTWheel wheel=\"" << id.wheel() << "\" />" << std::endl;
}
else if (alignableObjectId == align::AlignableDTBarrel) {
outputFile << " <DTBarrel />" << std::endl;
}
else throw cms::Exception("Alignment") << "Unknown DT Alignable StructureType" << std::endl;
}
} // end if not rawId
} // end if DT
else { // CSC
if (m_rawIds && rawId != 0) {
std::string typeName = AlignableObjectId::idToString(alignableObjectId);
if (alignableObjectId == align::AlignableDetUnit) typeName = std::string("CSCLayer");
outputFile << " <" << typeName << " rawId=\"" << rawId << "\" />" << std::endl;
}
else {
if (alignableObjectId == align::AlignableDetUnit) {
CSCDetId id(rawId);
outputFile << " <CSCLayer endcap=\"" << id.endcap() << "\" station=\"" << id.station() << "\" ring=\"" << id.ring() << "\" chamber=\"" << id.chamber() << "\" layer=\"" << id.layer() << "\" />" << std::endl;
}
else if (alignableObjectId == align::AlignableCSCChamber) {
CSCDetId id(rawId);
outputFile << " <CSCChamber endcap=\"" << id.endcap() << "\" station=\"" << id.station() << "\" ring=\"" << id.ring() << "\" chamber=\"" << id.chamber() << "\" />" << std::endl;
}
else {
CSCDetId id((*alignable)->id());
if (alignableObjectId == align::AlignableCSCRing) {
outputFile << " <CSCRing endcap=\"" << id.endcap() << "\" station=\"" << id.station() << "\" ring=\"" << id.ring() << "\" />" << std::endl;
}
else if (alignableObjectId == align::AlignableCSCStation) {
outputFile << " <CSCStation endcap=\"" << id.endcap() << "\" station=\"" << id.station() << "\" />" << std::endl;
}
else if (alignableObjectId == align::AlignableCSCEndcap) {
outputFile << " <CSCEndcap endcap=\"" << id.endcap() << "\" />" << std::endl;
}
else throw cms::Exception("Alignment") << "Unknown CSC Alignable StructureType" << std::endl;
}
} // end if not rawId
} // end if CSC
align::PositionType pos = (*alignable)->globalPosition();
align::RotationType rot = (*alignable)->globalRotation();
if (m_survey) {
pos = (*alignable)->survey()->position();
rot = (*alignable)->survey()->rotation();
}
std::string str_relativeto;
if (m_relativeto == 0) {
str_relativeto = std::string("none");
}
else if (m_relativeto == 1) {
if (ideal == ideals.end() || (*ideal)->alignableObjectId() != alignableObjectId || (*ideal)->id() != (*alignable)->id()) {
throw cms::Exception("Alignment") << "AlignableMuon and ideal_AlignableMuon are out of sync!" << std::endl;
}
align::PositionType idealPosition = (*ideal)->globalPosition();
align::RotationType idealRotation = (*ideal)->globalRotation();
pos = align::PositionType(idealRotation * (pos.basicVector() - idealPosition.basicVector()));
rot = rot * idealRotation.transposed();
str_relativeto = std::string("ideal");
bool csc_debug=0;
if (csc_debug && !DT) {
CSCDetId id(rawId);
if(id.endcap()==1 && id.station()==1 && id.ring()==1 && id.chamber()==33 ){
std::cout<<" investigating "<<id<<std::endl<<(*alignable)->globalRotation()<<std::endl<<std::endl
<<idealRotation.transposed()<<std::endl<<std::endl<<rot<<std::endl<<std::endl;
double phix = atan2(rot.yz(), rot.zz());
double phiy = asin(-rot.xz());
double phiz = atan2(rot.xy(), rot.xx());
std::cout << "phix=\"" << phix << "\" phiy=\"" << phiy << "\" phiz=\"" << phiz << std::endl;
align::EulerAngles eulerAngles = align::toAngles((*alignable)->globalRotation());
std::cout << "alpha=\"" << eulerAngles(1) << "\" beta=\"" << eulerAngles(2) << "\" gamma=\"" << eulerAngles(3) << std::endl;
eulerAngles = align::toAngles(idealRotation);
std::cout << "alpha=\"" << eulerAngles(1) << "\" beta=\"" << eulerAngles(2) << "\" gamma=\"" << eulerAngles(3) << std::endl;
eulerAngles = align::toAngles(rot);
std::cout << "alpha=\"" << eulerAngles(1) << "\" beta=\"" << eulerAngles(2) << "\" gamma=\"" << eulerAngles(3) << std::endl;
}
}
}
else if (m_relativeto == 2 && (*alignable)->mother() != NULL) {
align::PositionType globalPosition = (*alignable)->mother()->globalPosition();
align::RotationType globalRotation = (*alignable)->mother()->globalRotation();
pos = align::PositionType(globalRotation * (pos.basicVector() - globalPosition.basicVector()));
rot = rot * globalRotation.transposed();
str_relativeto = std::string("container");
}
else assert(false); // can't happen: see constructor
outputFile << " <setposition relativeto=\"" << str_relativeto << "\" "
<< "x=\"" << pos.x() << "\" y=\"" << pos.y() << "\" z=\"" << pos.z() << "\" ";
if (m_eulerAngles) {
align::EulerAngles eulerAngles = align::toAngles(rot);
outputFile << "alpha=\"" << eulerAngles(1) << "\" beta=\"" << eulerAngles(2) << "\" gamma=\"" << eulerAngles(3) << "\" />" << std::endl;
}
else {
// the angle convention originally used in alignment, also known as "non-standard Euler angles with a Z-Y-X convention"
// this also gets the sign convention right
double phix = atan2(rot.yz(), rot.zz());
double phiy = asin(-rot.xz());
double phiz = atan2(rot.xy(), rot.xx());
outputFile << "phix=\"" << phix << "\" phiy=\"" << phiy << "\" phiz=\"" << phiz << "\" />" << std::endl;
}
if (m_survey) {
align::ErrorMatrix err = (*alignable)->survey()->errors();
outputFile << " <setsurveyerr"
<< " xx=\"" << err(0,0) << "\" xy=\"" << err(0,1) << "\" xz=\"" << err(0,2) << "\" xa=\"" << err(0,3) << "\" xb=\"" << err(0,4) << "\" xc=\"" << err(0,5)
<< "\" yy=\"" << err(1,1) << "\" yz=\"" << err(1,2) << "\" ya=\"" << err(1,3) << "\" yb=\"" << err(1,4) << "\" yc=\"" << err(1,5)
<< "\" zz=\"" << err(2,2) << "\" za=\"" << err(2,3) << "\" zb=\"" << err(2,4) << "\" zc=\"" << err(2,5)
<< "\" aa=\"" << err(3,3) << "\" ab=\"" << err(3,4) << "\" ac=\"" << err(3,5)
<< "\" bb=\"" << err(4,4) << "\" bc=\"" << err(4,5)
<< "\" cc=\"" << err(5,5) << "\" />" << std::endl;
}
else if (rawId != 0) {
CLHEP::HepSymMatrix err = errors[(*alignable)->id()];
outputFile << " <setape xx=\"" << err(1,1) << "\" xy=\"" << err(1,2) << "\" xz=\"" << err(1,3)
<< "\" yy=\"" << err(2,2) << "\" yz=\"" << err(2,3) << "\" zz=\"" << err(3,3) << "\" />" << std::endl;
}
outputFile << "</operation>" << std::endl << std::endl;
} // end if not suppressed
// write superstructures before substructures: this is important because <setape> overwrites all substructures' APEs
if (ideal != ideals.end()) {
align::Alignables components = (*alignable)->components();
align::Alignables ideal_components = (*ideal)->components();
writeComponents(components, ideal_components, errors, outputFile, DT);
++ideal; // important for synchronization in the "for" loop!
}
else {
align::Alignables components = (*alignable)->components();
align::Alignables dummy;
writeComponents(components, dummy, errors, outputFile, DT);
}
} // end loop over alignables
}
bool MuonAlignmentOutputXML::m_eulerAngles [private] |
Definition at line 59 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
std::string MuonAlignmentOutputXML::m_fileName [private] |
Definition at line 57 of file MuonAlignmentOutputXML.h.
Referenced by write().
int MuonAlignmentOutputXML::m_precision [private] |
Definition at line 60 of file MuonAlignmentOutputXML.h.
Referenced by write().
bool MuonAlignmentOutputXML::m_rawIds [private] |
Definition at line 59 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
int MuonAlignmentOutputXML::m_relativeto [private] |
Definition at line 58 of file MuonAlignmentOutputXML.h.
Referenced by MuonAlignmentOutputXML(), write(), and writeComponents().
bool MuonAlignmentOutputXML::m_suppressCSCChambers [private] |
Definition at line 62 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_suppressCSCEndcaps [private] |
Definition at line 62 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_suppressCSCLayers [private] |
Definition at line 62 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_suppressCSCRings [private] |
Definition at line 62 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_suppressCSCStations [private] |
Definition at line 62 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_suppressDTBarrel [private] |
Definition at line 61 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_suppressDTChambers [private] |
Definition at line 61 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_suppressDTLayers [private] |
Definition at line 61 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_suppressDTStations [private] |
Definition at line 61 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_suppressDTSuperLayers [private] |
Definition at line 61 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_suppressDTWheels [private] |
Definition at line 61 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
bool MuonAlignmentOutputXML::m_survey [private] |
Definition at line 59 of file MuonAlignmentOutputXML.h.
Referenced by writeComponents().
1.7.3