diff --git a/Alignment/MuonAlignment/test/BuildFile.xml b/Alignment/MuonAlignment/test/BuildFile.xml deleted file mode 100644 index 49ff1177511b4..0000000000000 --- a/Alignment/MuonAlignment/test/BuildFile.xml +++ /dev/null @@ -1,12 +0,0 @@ - - - - - - - - - - - - diff --git a/Alignment/MuonAlignment/test/ComputeTransformation.cpp b/Alignment/MuonAlignment/test/ComputeTransformation.cpp deleted file mode 100644 index 06a11fe27f594..0000000000000 --- a/Alignment/MuonAlignment/test/ComputeTransformation.cpp +++ /dev/null @@ -1,86 +0,0 @@ -#include - -#include "FWCore/Framework/interface/EDAnalyzer.h" -#include "FWCore/Framework/interface/EventSetup.h" -#include "FWCore/Framework/interface/ESHandle.h" -#include "FWCore/Framework/interface/MakerMacros.h" - -#include "Geometry/DTGeometry/interface/DTGeometry.h" -#include "Geometry/CSCGeometry/interface/CSCGeometry.h" -#include "Geometry/Records/interface/MuonGeometryRecord.h" -#include "DataFormats/MuonDetId/interface/DTChamberId.h" -#include "DataFormats/MuonDetId/interface/CSCDetId.h" - -// -// -// class declaration -// - -class ComputeTransformation : public edm::EDAnalyzer { -public: - explicit ComputeTransformation(const edm::ParameterSet&); - ~ComputeTransformation(); - - virtual void analyze(const edm::Event&, const edm::EventSetup&); - -private: - // ----------member data --------------------------- - std::string m_fileName; -}; - -// -// constructors and destructor -// -ComputeTransformation::ComputeTransformation(const edm::ParameterSet& iConfig) - : m_fileName(iConfig.getParameter("fileName")) {} - -ComputeTransformation::~ComputeTransformation() {} - -void ComputeTransformation::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) { - edm::ESHandle dtGeometry; - iSetup.get().get(dtGeometry); - - edm::ESHandle cscGeometry; - iSetup.get().get(cscGeometry); - - std::ofstream output; - output.open(m_fileName.c_str()); - output << "rotation = {}" << std::endl; - - for (int wheel = -2; wheel <= 2; wheel++) { - for (int station = 1; station <= 4; station++) { - for (int sector = 1; sector <= 14; sector++) { - if (station != 4 and sector > 12) - continue; - DTChamberId id(wheel, station, sector); - // globalcoords = rot * localcoords - Surface::RotationType rot = dtGeometry->idToDet(id)->surface().rotation(); - output << "rotation[\"DT\", " << wheel << ", " << station << ", 0, " << sector << "] = [[" << rot.xx() << ", " - << rot.xy() << ", " << rot.xz() << "], [" << rot.yx() << ", " << rot.yy() << ", " << rot.yz() << "], [" - << rot.zx() << ", " << rot.zy() << ", " << rot.zz() << "]]" << std::endl; - } - } - } - - for (int endcap = 1; endcap <= 2; endcap++) { - for (int station = 1; station <= 4; station++) { - for (int ring = 1; ring <= 3; ring++) { - if (station > 1 and ring == 3) - continue; - for (int sector = 1; sector <= 36; sector++) { - if (station > 1 && ring == 1 && sector > 18) - continue; - CSCDetId id(endcap, station, ring, sector); - // globalcoords = rot * localcoords - Surface::RotationType rot = cscGeometry->idToDet(id)->surface().rotation(); - output << "rotation[\"CSC\", " << endcap << ", " << station << ", " << ring << ", " << sector << "] = [[" - << rot.xx() << ", " << rot.xy() << ", " << rot.xz() << "], [" << rot.yx() << ", " << rot.yy() << ", " - << rot.yz() << "], [" << rot.zx() << ", " << rot.zy() << ", " << rot.zz() << "]]" << std::endl; - } - } - } - } -} - -//define this as a plug-in -DEFINE_FWK_MODULE(ComputeTransformation); diff --git a/Alignment/MuonAlignment/test/MuonGeometryArrange.py b/Alignment/MuonAlignment/test/MuonGeometryArrange.py deleted file mode 100644 index c2395db06fede..0000000000000 --- a/Alignment/MuonAlignment/test/MuonGeometryArrange.py +++ /dev/null @@ -1,106 +0,0 @@ -import FWCore.ParameterSet.Config as cms -process =cms.Process("TEST") - -#Ideal geometry -process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi") -process.load('Configuration.Geometry.GeometryExtended2021_cff') -process.load("Geometry.MuonNumbering.muonNumberingInitialization_cfi") -process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff") - -from Configuration.AlCa.GlobalTag import GlobalTag -process.GlobalTag = GlobalTag(process.GlobalTag, "auto:phase1_2021_design") - -import Geometry.DTGeometryBuilder.dtGeometryDB_cfi -import Geometry.CSCGeometryBuilder.cscGeometryDB_cfi -import Geometry.GEMGeometryBuilder.gemGeometryDB_cfi - -process.DTGeometryIdeal = Geometry.DTGeometryBuilder.dtGeometryDB_cfi.DTGeometryESModule.clone() -process.DTGeometryIdeal.appendToDataLabel = 'MuonGeometryArrangeGeomIdeal' -process.DTGeometryIdeal.applyAlignment = cms.bool(False) -process.CSCGeometryIdeal = Geometry.CSCGeometryBuilder.cscGeometryDB_cfi.CSCGeometryESModule.clone() -process.CSCGeometryIdeal.appendToDataLabel = 'MuonGeometryArrangeGeomIdeal' -process.CSCGeometryIdeal.applyAlignment = cms.bool(False) -process.GEMGeometryIdeal = Geometry.GEMGeometryBuilder.gemGeometryDB_cfi.GEMGeometryESModule.clone() -process.GEMGeometryIdeal.appendToDataLabel = 'MuonGeometryArrangeGeomIdeal' -process.GEMGeometryIdeal.applyAlignment = cms.bool(False) - -process.DTGeometryMuonGeometryArrange1 = Geometry.DTGeometryBuilder.dtGeometryDB_cfi.DTGeometryESModule.clone() -process.DTGeometryMuonGeometryArrange1.appendToDataLabel = 'MuonGeometryArrangeLabel1' -process.DTGeometryMuonGeometryArrange1.applyAlignment = cms.bool(False) -process.CSCGeometryMuonGeometryArrange1 = Geometry.CSCGeometryBuilder.cscGeometryDB_cfi.CSCGeometryESModule.clone() -process.CSCGeometryMuonGeometryArrange1.appendToDataLabel = 'MuonGeometryArrangeLabel1' -process.CSCGeometryMuonGeometryArrange1.applyAlignment = cms.bool(False) -process.GEMGeometryMuonGeometryArrange1 = Geometry.GEMGeometryBuilder.gemGeometryDB_cfi.GEMGeometryESModule.clone() -process.GEMGeometryMuonGeometryArrange1.appendToDataLabel = 'MuonGeometryArrangeLabel1' -process.GEMGeometryMuonGeometryArrange1.applyAlignment = cms.bool(False) - -process.DTGeometryMuonGeometryArrange2 = Geometry.DTGeometryBuilder.dtGeometryDB_cfi.DTGeometryESModule.clone() -process.DTGeometryMuonGeometryArrange2.appendToDataLabel = 'MuonGeometryArrangeLabel2' -process.DTGeometryMuonGeometryArrange2.applyAlignment = cms.bool(False) -process.CSCGeometryMuonGeometryArrange2 = Geometry.CSCGeometryBuilder.cscGeometryDB_cfi.CSCGeometryESModule.clone() -process.CSCGeometryMuonGeometryArrange2.appendToDataLabel = 'MuonGeometryArrangeLabel2' -process.CSCGeometryMuonGeometryArrange2.applyAlignment = cms.bool(False) -process.GEMGeometryMuonGeometryArrange2 = Geometry.GEMGeometryBuilder.gemGeometryDB_cfi.GEMGeometryESModule.clone() -process.GEMGeometryMuonGeometryArrange2.appendToDataLabel = 'MuonGeometryArrangeLabel2' -process.GEMGeometryMuonGeometryArrange2.applyAlignment = cms.bool(False) - -process.DTGeometryMuonGeometryArrange2a = Geometry.DTGeometryBuilder.dtGeometryDB_cfi.DTGeometryESModule.clone() -process.DTGeometryMuonGeometryArrange2a.appendToDataLabel = 'MuonGeometryArrangeLabel2a' -process.DTGeometryMuonGeometryArrange2a.applyAlignment = cms.bool(False) -process.CSCGeometryMuonGeometryArrange2a = Geometry.CSCGeometryBuilder.cscGeometryDB_cfi.CSCGeometryESModule.clone() -process.CSCGeometryMuonGeometryArrange2a.appendToDataLabel = 'MuonGeometryArrangeLabel2a' -process.CSCGeometryMuonGeometryArrange2a.applyAlignment = cms.bool(False) -process.GEMGeometryMuonGeometryArrange2a = Geometry.GEMGeometryBuilder.gemGeometryDB_cfi.GEMGeometryESModule.clone() -process.GEMGeometryMuonGeometryArrange2a.appendToDataLabel = 'MuonGeometryArrangeLabel2a' -process.GEMGeometryMuonGeometryArrange2a.applyAlignment = cms.bool(False) - -process.MessageLogger = cms.Service("MessageLogger", - cerr = cms.untracked.PSet( - threshold = cms.untracked.string('INFO') - ), - files = cms.untracked.PSet( - info_txt = cms.untracked.PSet( - threshold = cms.untracked.string('INFO') - ) - ) -) - - -process.source = cms.Source("EmptySource") -process.maxEvents = cms.untracked.PSet( - input = cms.untracked.int32(1) - ) - - -# Full configuration for Muon Geometry Comparison Tool -#process.MuonGeometryCompare = cms.EDFilter("MuonGeometryArrange", -process.MuonGeometryCompare = cms.EDAnalyzer("MuonGeometryArrange", - outputFile = cms.untracked.string('output.root'), - - detIdFlag = cms.untracked.bool(False), - detIdFlagFile = cms.untracked.string('blah.txt'), - weightById = cms.untracked.bool(False), - levels = cms.untracked.vstring('Det'), - weightBy = cms.untracked.string('SELF'), - weightByIdFile = cms.untracked.string('blah2.txt'), - treeName = cms.untracked.string('alignTree'), -# Root input files are not used yet. - inputROOTFile1 = cms.untracked.string('IDEAL'), - inputROOTFile2 = cms.untracked.string('idealmuon2.root'), -# Geometries are read from xml files -# inputXMLCurrent = cms.untracked.string('B.xml'), -# inputXMLCurrent = cms.untracked.string('A.xml'), -# inputXMLCurrent = cms.untracked.string('moveRing.xml'), -# inputXMLCurrent = cms.untracked.string('movedRing.xml'), -# inputXMLCurrent = cms.untracked.string('fiddleMuon.xml'), -# inputXMLCurrent = cms.untracked.string('fiddle2Muon.xml'), - inputXMLCurrent = cms.untracked.string('fiddle3Muon.xml'), - inputXMLReference = cms.untracked.string('idealMuon.xml'), -# A few defaults. You pick. - endcapNumber = cms.untracked.int32(2), - stationNumber = cms.untracked.int32(3), - ringNumber = cms.untracked.int32(2) -) - -process.p = cms.Path( process.MuonGeometryCompare ) - diff --git a/Alignment/MuonAlignment/test/TestMisalign.cpp b/Alignment/MuonAlignment/test/TestMisalign.cpp deleted file mode 100644 index 8b96e0605940b..0000000000000 --- a/Alignment/MuonAlignment/test/TestMisalign.cpp +++ /dev/null @@ -1,130 +0,0 @@ -// -*- C++ -*- -// -// Package: TestMisalign -// Class: TestMisalign -// -// -// Description: Module to test the Alignment software -// -// - -// system include files -#include -#include -#include -#include - -// user include files -#include "FWCore/MessageLogger/interface/MessageLogger.h" -#include "FWCore/Framework/interface/EDAnalyzer.h" - -#include "FWCore/Framework/interface/EventSetup.h" -#include "FWCore/Framework/interface/ESHandle.h" -#include "FWCore/Framework/interface/MakerMacros.h" - -#include "FWCore/ParameterSet/interface/ParameterSet.h" - -#include - -#include "Geometry/DTGeometry/interface/DTGeometry.h" -#include "Geometry/CSCGeometry/interface/CSCGeometry.h" - -// -// -// class declaration -// - -class TestMisalign : public edm::EDAnalyzer { -public: - explicit TestMisalign(const edm::ParameterSet&); - ~TestMisalign(); - - virtual void analyze(const edm::Event&, const edm::EventSetup&); - -private: - void fillTree(const GeomDet* geomDet); - - // ----------member data --------------------------- - TTree* theTree; - TFile* theFile; - float x, y, z, phi, theta, length, thick, width; - int Id_; - TRotMatrix* rot; -}; - -// -// constructors and destructor -// -TestMisalign::TestMisalign(const edm::ParameterSet& iConfig) { - // Open root file and define tree - std::string fileName = iConfig.getUntrackedParameter("fileName", "test.root"); - theFile = new TFile(fileName.c_str(), "RECREATE"); - theTree = new TTree("theTree", "Detector units positions"); - - theTree->Branch("Id", &Id_, "Id/I"); - theTree->Branch("x", &x, "x/F"); - theTree->Branch("y", &y, "y/F"); - theTree->Branch("z", &z, "z/F"); - theTree->Branch("phi", &phi, "phi/F"); - theTree->Branch("theta", &theta, "theta/F"); - theTree->Branch("length", &length, "length/F"); - theTree->Branch("width", &width, "width/F"); - theTree->Branch("thick", &thick, "thick/F"); - rot = 0; - theTree->Branch("rot", "TRotMatrix", &rot); -} - -TestMisalign::~TestMisalign() { - theTree->Write(); - theFile->Close(); -} - -void TestMisalign::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) { - edm::LogInfo("MuonAlignment") << "Starting!"; - - // - // Retrieve muon geometry from event setup - // - edm::ESHandle pDT; - edm::ESHandle pCSC; - - iSetup.get().get(pDT); - iSetup.get().get(pCSC); - - // Now loop on detector units, and store position and orientation - for (auto iGeomDet = pDT->dets().begin(); iGeomDet != pDT->dets().end(); iGeomDet++) - this->fillTree(*iGeomDet); - for (auto iGeomDet = pCSC->dets().begin(); iGeomDet != pCSC->dets().end(); iGeomDet++) - this->fillTree(*iGeomDet); - - edm::LogInfo("MuonAlignment") << "Done!"; -} - -//__________________________________________________________________________________________________ -void TestMisalign::fillTree(const GeomDet* geomDet) { - Id_ = geomDet->geographicalId().rawId(); - x = geomDet->position().x(); - y = geomDet->position().y(); - z = geomDet->position().z(); - phi = geomDet->surface().normalVector().phi(); - theta = geomDet->surface().normalVector().theta(); - length = geomDet->surface().bounds().length(); - width = geomDet->surface().bounds().width(); - thick = geomDet->surface().bounds().thickness(); - - double matrix[9] = {geomDet->rotation().xx(), - geomDet->rotation().xy(), - geomDet->rotation().xz(), - geomDet->rotation().yx(), - geomDet->rotation().yy(), - geomDet->rotation().yz(), - geomDet->rotation().zx(), - geomDet->rotation().zy(), - geomDet->rotation().zz()}; - rot = new TRotMatrix("rot", "rot", matrix); - - theTree->Fill(); -} - -//define this as a plug-in -DEFINE_FWK_MODULE(TestMisalign); diff --git a/Alignment/MuonAlignment/test/TestMuonAnalyzer.cpp b/Alignment/MuonAlignment/test/TestMuonAnalyzer.cpp deleted file mode 100644 index 674fbe1cc5cbf..0000000000000 --- a/Alignment/MuonAlignment/test/TestMuonAnalyzer.cpp +++ /dev/null @@ -1,137 +0,0 @@ -// -*- C++ -*- -// -// Package: TestMuonAnalyzer -// Class: TestMuonAnalyzer -// -// -// Description: Module to test the Alignment software -// -// - -// system include files -#include -#include -#include -#include - -// user include files -#include "FWCore/Framework/interface/Frameworkfwd.h" -#include "FWCore/MessageLogger/interface/MessageLogger.h" -#include "FWCore/Framework/interface/EDAnalyzer.h" - -#include "FWCore/Framework/interface/Event.h" -#include "FWCore/Framework/interface/EventSetup.h" -#include "FWCore/Framework/interface/ESHandle.h" -#include "FWCore/Framework/interface/MakerMacros.h" - -#include "FWCore/ParameterSet/interface/ParameterSet.h" - -#include "Alignment/MuonAlignment/interface/AlignableMuon.h" -#include "Alignment/MuonAlignment/interface/AlignableDTChamber.h" -#include "Alignment/MuonAlignment/interface/AlignableCSCChamber.h" - -#include "Geometry/DTGeometry/interface/DTGeometry.h" -#include "Geometry/DTGeometry/interface/DTChamber.h" -#include "Geometry/CSCGeometry/interface/CSCGeometry.h" -#include "Geometry/CSCGeometry/interface/CSCChamber.h" -#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h" -#include "Geometry/Records/interface/MuonGeometryRecord.h" -// -// -// class declaration -// - -class TestMuonAnalyzer : public edm::EDAnalyzer { -public: - explicit TestMuonAnalyzer(const edm::ParameterSet&); - ~TestMuonAnalyzer(); - - virtual void analyze(const edm::Event&, const edm::EventSetup&); - -private: - void fillTree(const GeomDet* geomDet); - - // ----------member data --------------------------- - TTree* theTree; - TFile* theFile; - float x, y, z, phi, theta, length, thick, width; - int Id_; - TRotMatrix* rot; -}; - -// -// constructors and destructor -// -TestMuonAnalyzer::TestMuonAnalyzer(const edm::ParameterSet& iConfig) { - // Open root file and define tree - std::string fileName = iConfig.getUntrackedParameter("fileName", "test.root"); - theFile = new TFile(fileName.c_str(), "RECREATE"); - theTree = new TTree("theTree", "Detector units positions"); - - theTree->Branch("Id", &Id_, "Id/I"); - theTree->Branch("x", &x, "x/F"); - theTree->Branch("y", &y, "y/F"); - theTree->Branch("z", &z, "z/F"); - theTree->Branch("phi", &phi, "phi/F"); - theTree->Branch("theta", &theta, "theta/F"); - theTree->Branch("length", &length, "length/F"); - theTree->Branch("width", &width, "width/F"); - theTree->Branch("thick", &thick, "thick/F"); - rot = 0; - theTree->Branch("rot", "TRotMatrix", &rot); -} - -TestMuonAnalyzer::~TestMuonAnalyzer() { - theTree->Write(); - theFile->Close(); -} - -void TestMuonAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) { - edm::LogInfo("MuonAlignment") << "Starting!"; - - // - // Retrieve muon geometry from event setup - // - edm::ESHandle pDT; - edm::ESHandle pCSC; - - iSetup.get().get(pDT); - iSetup.get().get(pCSC); - - // Now loop on detector units, and store position and orientation - for (auto iGeomDet = pDT->dets().begin(); iGeomDet != pDT->dets().end(); iGeomDet++) - this->fillTree(*iGeomDet); - for (auto iGeomDet = pCSC->dets().begin(); iGeomDet != pCSC->dets().end(); iGeomDet++) - this->fillTree(*iGeomDet); - - edm::LogInfo("MuonAlignment") << "Done!"; -} - -//__________________________________________________________________________________________________ -void TestMuonAnalyzer::fillTree(const GeomDet* geomDet) { - Id_ = geomDet->geographicalId().rawId(); - x = geomDet->position().x(); - y = geomDet->position().y(); - z = geomDet->position().z(); - phi = geomDet->surface().normalVector().phi(); - theta = geomDet->surface().normalVector().theta(); - length = geomDet->surface().bounds().length(); - width = geomDet->surface().bounds().width(); - thick = geomDet->surface().bounds().thickness(); - - double matrix[9] = {geomDet->rotation().xx(), - geomDet->rotation().xy(), - geomDet->rotation().xz(), - geomDet->rotation().yx(), - geomDet->rotation().yy(), - geomDet->rotation().yz(), - geomDet->rotation().zx(), - geomDet->rotation().zy(), - geomDet->rotation().zz()}; - rot = new TRotMatrix("rot", "rot", matrix); // mem. leak? - - theTree->Fill(); -} - -//define this as a plug-in -DEFINE_FWK_MODULE(TestMuonAnalyzer); diff --git a/Alignment/MuonAlignment/test/TestMuonHierarchy.cpp b/Alignment/MuonAlignment/test/TestMuonHierarchy.cpp deleted file mode 100644 index 6e2360a419706..0000000000000 --- a/Alignment/MuonAlignment/test/TestMuonHierarchy.cpp +++ /dev/null @@ -1,114 +0,0 @@ -// -*- C++ -*- -// -// Package: TestMuonHierarchy -// Class: TestMuonHierarchy -// -// -// Description: Module to test the Alignment software -// -// -// Original Author: Frederic Ronga -// Created: March 16, 2006 -// - -// system include files -#include -#include -#include - -// user include files -#include "FWCore/Framework/interface/EDAnalyzer.h" -#include "FWCore/Framework/interface/MakerMacros.h" -#include "FWCore/Framework/interface/EventSetup.h" -#include "FWCore/Framework/interface/ESHandle.h" -#include "FWCore/MessageLogger/interface/MessageLogger.h" - -#include "Geometry/Records/interface/MuonGeometryRecord.h" - -#include "Alignment/MuonAlignment/interface/AlignableMuon.h" -#include "Alignment/CommonAlignment/interface/AlignableObjectId.h" - -static const int kLEAD_WIDTH = 40; // First field width - -// -// class declaration -// -class TestMuonHierarchy : public edm::EDAnalyzer { -public: - explicit TestMuonHierarchy(const edm::ParameterSet&) {} - - virtual void analyze(const edm::Event&, const edm::EventSetup&); - -private: - // ----------member data --------------------------- - void dumpAlignable(const Alignable*, unsigned int, unsigned int); - void printInfo(const Alignable*, unsigned int); - - std::unique_ptr alignableMuon_; - std::string leaders_, blank_, filled_; - - edm::ESGetToken esTokenDT_; - edm::ESGetToken esTokenCSC_; - edm::ESGetToken esTokenGEM_; -}; - -void TestMuonHierarchy::analyze(const edm::Event&, const edm::EventSetup& setup) { - edm::LogInfo("MuonHierarchy") << "Starting!"; - edm::ESHandle dtGeometry = setup.getHandle(esTokenDT_); - edm::ESHandle cscGeometry = setup.getHandle(esTokenCSC_); - edm::ESHandle gemGeometry = setup.getHandle(esTokenGEM_); - - alignableMuon_ = std::make_unique(&(*dtGeometry), &(*cscGeometry), &(*gemGeometry)); - - leaders_ = ""; - blank_ = " "; // These two... - filled_ = "| "; // ... must have the same length - - // Now dump mother of each alignable - //const Alignable* alignable = (&(*theAlignableMuon))->pixelHalfBarrels()[0]; - dumpAlignable(alignableMuon_.get(), 1, 1); - - edm::LogInfo("MuonAlignment") << "Done!"; -} - -//__________________________________________________________________________________________________ -// Recursive loop on alignable hierarchy -void TestMuonHierarchy::dumpAlignable(const Alignable* alignable, unsigned int idau, unsigned int ndau) { - printInfo(alignable, idau); - - if (ndau != idau) - leaders_ += filled_; - else - leaders_ += blank_; - - const align::Alignables& comps = alignable->components(); - if (unsigned int ndau_ = comps.size()) { - unsigned int idau_ = 0; - for (align::Alignables::const_iterator iter = comps.begin(); iter != comps.end(); ++iter) - dumpAlignable(*iter, ++idau_, ndau_); - } - - leaders_ = leaders_.substr(0, leaders_.length() - blank_.length()); -} - -//__________________________________________________________________________________________________ -// Do the actual printout -void TestMuonHierarchy::printInfo(const Alignable* alignable, unsigned int idau) { - int width = kLEAD_WIDTH - leaders_.length(); - - std::ostringstream name, pos, rot; - - name << alignableMuon_->objectIdProvider().idToString(alignable->alignableObjectId()) << idau; - - // Position - pos.setf(std::ios::fixed); - pos << "(" << std::right << std::setw(8) << std::setprecision(4) << alignable->globalPosition().x() << "," - << std::setw(8) << std::setprecision(4) << alignable->globalPosition().y() << "," << std::setw(8) - << std::setprecision(4) << alignable->globalPosition().z() << ")"; - - edm::LogVerbatim("DumpAlignable") << leaders_ << "+-> " << std::setw(width) << std::left << name.str() << " | " - << std::setw(3) << std::left << alignable->components().size() << " | " - << std::setw(11) << std::left << alignable->id() << " | " << pos.str(); -} -//define this as a plug-in -DEFINE_FWK_MODULE(TestMuonHierarchy); diff --git a/Alignment/MuonAlignment/test/TestReader.cpp b/Alignment/MuonAlignment/test/TestReader.cpp deleted file mode 100644 index ba4b9a43633c2..0000000000000 --- a/Alignment/MuonAlignment/test/TestReader.cpp +++ /dev/null @@ -1,263 +0,0 @@ -// -*- C++ -*- -// -// -// Description: Module to test the Alignment software -// -// - -// system include files -#include -#include -#include -#include - -// user include files -#include "FWCore/MessageLogger/interface/MessageLogger.h" -#include "FWCore/Framework/interface/EDAnalyzer.h" - -#include "FWCore/Framework/interface/EventSetup.h" -#include "FWCore/Framework/interface/ESHandle.h" -#include "FWCore/Framework/interface/ESTransientHandle.h" -#include "FWCore/Framework/interface/MakerMacros.h" - -#include "CondFormats/Alignment/interface/Alignments.h" -#include "CondFormats/Alignment/interface/AlignTransform.h" -#include "CondFormats/AlignmentRecord/interface/DTAlignmentRcd.h" -#include "CondFormats/AlignmentRecord/interface/CSCAlignmentRcd.h" -#include "CondFormats/Alignment/interface/AlignmentErrorsExtended.h" -#include "CondFormats/Alignment/interface/AlignTransformErrorExtended.h" -#include "CondFormats/AlignmentRecord/interface/DTAlignmentErrorExtendedRcd.h" -#include "CondFormats/AlignmentRecord/interface/CSCAlignmentErrorExtendedRcd.h" - -#include "Geometry/DTGeometry/interface/DTGeometry.h" -#include "Geometry/CSCGeometry/interface/CSCGeometry.h" -#include "Geometry/Records/interface/MuonGeometryRecord.h" - -#include "Alignment/MuonAlignment/interface/AlignableMuon.h" -#include "Alignment/CommonAlignment/interface/Utilities.h" -#include "Alignment/CommonAlignmentAlgorithm/interface/AlignmentParameterBuilder.h" - -#include "DataFormats/MuonDetId/interface/CSCDetId.h" - -// -// -// class declaration -// - -class TestMuonReader : public edm::EDAnalyzer { -public: - explicit TestMuonReader(const edm::ParameterSet&); - ~TestMuonReader(); - - void recursiveGetMuChambers(const align::Alignables& composite, align::Alignables& chambers, int kind); - align::EulerAngles toPhiXYZ(const align::RotationType&); - virtual void analyze(const edm::Event&, const edm::EventSetup&); - -private: - // ----------member data --------------------------- - TTree* theTree; - TFile* theFile; - float x, y, z, phi, theta, length, thick, width; - TRotMatrix* rot; - - edm::ESGetToken esTokenDT_; - edm::ESGetToken esTokenCSC_; - edm::ESGetToken esTokenGEM_; -}; - -// -// constructors and destructor -// -TestMuonReader::TestMuonReader(const edm::ParameterSet& iConfig) - : theTree(0), - theFile(0), - x(0.), - y(0.), - z(0.), - phi(0.), - theta(0.), - length(0.), - thick(0.), - width(0.), - rot(0), - esTokenDT_(esConsumes(edm::ESInputTag("", "idealForTestReader"))), - esTokenCSC_(esConsumes(edm::ESInputTag("", "idealForTestReader"))), - esTokenGEM_(esConsumes(edm::ESInputTag("", "idealForTestReader"))) {} - -TestMuonReader::~TestMuonReader() {} - -void TestMuonReader::recursiveGetMuChambers(const align::Alignables& composites, - align::Alignables& chambers, - int kind) { - for (const auto& cit : composites) { - if (cit->alignableObjectId() == kind) { - chambers.push_back(cit); - continue; - } else { - const auto& components = cit->components(); - recursiveGetMuChambers(components, chambers, kind); - } - } -} - -align::EulerAngles TestMuonReader::toPhiXYZ(const align::RotationType& rot) { - align::EulerAngles angles(3); - angles(1) = std::atan2(rot.yz(), rot.zz()); - angles(2) = std::asin(-rot.xz()); - angles(3) = std::atan2(rot.xy(), rot.xx()); - return angles; -} - -void TestMuonReader::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) { - // first, get chamber alignables from ideal geometry: - edm::ESHandle dtGeometry = iSetup.getHandle(esTokenDT_); - edm::ESHandle cscGeometry = iSetup.getHandle(esTokenCSC_); - edm::ESHandle gemGeometry = iSetup.getHandle(esTokenGEM_); - - AlignableMuon ideal_alignableMuon(&(*dtGeometry), &(*cscGeometry), &(*gemGeometry)); - - const auto& ideal_barrels = ideal_alignableMuon.DTBarrel(); - const auto& ideal_endcaps = ideal_alignableMuon.CSCEndcaps(); - - align::Alignables ideal_mb_chambers, ideal_me_chambers; - recursiveGetMuChambers(ideal_barrels, ideal_mb_chambers, align::AlignableDTChamber); - recursiveGetMuChambers(ideal_endcaps, ideal_me_chambers, align::AlignableCSCChamber); - //std::cout<<" #ideals dts="< dtAlignments; - iSetup.get().get( dtAlignments ); - edm::ESHandle dtAlignmentErrorsExtended; - iSetup.get().get( dtAlignmentErrorsExtended ); - - for ( std::vector::const_iterator it = dtAlignments->m_align.begin(); - it != dtAlignments->m_align.end(); it++ ) - { - CLHEP::HepRotation rot( (*it).rotation() ); - align::RotationType rotation( rot.xx(), rot.xy(), rot.xz(), - rot.yx(), rot.yy(), rot.yz(), - rot.zx(), rot.zy(), rot.zz() ); - - std::cout << (*it).rawId() - << " " << (*it).translation().x() - << " " << (*it).translation().y() - << " " << (*it).translation().z() - << " " << rotation.xx() << " " << rotation.xy() << " " << rotation.xz() - << " " << rotation.yx() << " " << rotation.yy() << " " << rotation.yz() - << " " << rotation.zx() << " " << rotation.zy() << " " << rotation.zz() - << std::endl; - - } - std::cout << std::endl << "----------------------" << std::endl; - - for ( std::vector::const_iterator it = dtAlignmentErrorsExtended->m_alignError.begin(); - it != dtAlignmentErrorsExtended->m_alignError.end(); it++ ) - { - CLHEP::HepSymMatrix error = (*it).matrix(); - std::cout << (*it).rawId() << " "; - for ( int i=0; i cscAlignments; - iSetup.get().get(cscAlignments); - //edm::ESHandle cscAlignmentErrorsExtended; - //iSetup.get().get( cscAlignmentErrorsExtended ); - - //align::Alignables::const_iterator csc_ideal = ideal_endcaps.begin(); - std::cout << std::setprecision(3) << std::fixed; - //std::cout<<" lens : "<m_align.size()<::const_iterator it = cscAlignments->m_align.begin(); - it != cscAlignments->m_align.end(); - it++) { - CSCDetId id((*it).rawId()); - if (id.layer() > 0) - continue; // look at chambers only, skip layers - - if (id.station() == 1 && id.ring() == 4) - continue; // not interested in duplicated ME1/4 - - char nme[100]; - sprintf(nme, "%d/%d/%02d", id.station(), id.ring(), id.chamber()); - std::string me = "ME+"; - if (id.endcap() == 2) - me = "ME-"; - me += nme; - - // find this chamber in ideal geometry - const Alignable* ideal{nullptr}; - for (const auto& cideal : ideal_me_chambers) - if (cideal->geomDetId().rawId() == (*it).rawId()) { - ideal = cideal; - break; - } - if (ideal == nullptr) { - std::cout << " no ideal chamber for " << id << std::endl; - continue; - } - - //if (ideal->geomDetId().rawId() != (*it).rawId()) std::cout<<" badid : "<<(*csc_ideal)->geomDetId().rawId()<<" "<<(*it).rawId()<globalPosition(); - align::RotationType idealRotation = ideal->globalRotation(); - //align::EulerAngles abg = align::toAngles(idealRotation); - //std::cout << me <<" "<< (*it).rawId()<<" "<0.00001 || fabs(rabg[1]-rxyz[1])>0.00001 ||fabs(rabg[1]-rxyz[1])>0.00001) - // std::cout << me <<" large angle diff = "<::const_iterator it = cscAlignmentErrorsExtended->m_alignError.begin(); - it != cscAlignmentErrorsExtended->m_alignError.end(); it++ ) - { - CLHEP::HepSymMatrix error = (*it).matrix(); - std::cout << (*it).rawId() << " "; - for ( int i=0; i -#include -#include -#include - -// user include files -#include "FWCore/MessageLogger/interface/MessageLogger.h" -#include "FWCore/Framework/interface/EDAnalyzer.h" - -#include "FWCore/Framework/interface/EventSetup.h" -#include "FWCore/Framework/interface/ESTransientHandle.h" -#include "FWCore/Framework/interface/ESHandle.h" -#include "FWCore/Framework/interface/MakerMacros.h" - -#include "FWCore/ParameterSet/interface/ParameterSet.h" - -#include "Alignment/MuonAlignment/interface/AlignableMuon.h" -#include - -#include "Geometry/DTGeometry/interface/DTGeometry.h" - -#include "DataFormats/GeometrySurface/interface/Surface.h" - -#include - -// -// -// class declaration -// - -class TestRotation : public edm::EDAnalyzer { -public: - explicit TestRotation(const edm::ParameterSet&); - ~TestRotation(); - - virtual void analyze(const edm::Event&, const edm::EventSetup&); - -private: - // ----------member data --------------------------- - TTree* theTree; - TFile* theFile; - float x, y, z, phi, theta, length, thick, width; - TRotMatrix* dir; - - edm::ESGetToken esTokenDT_; - edm::ESGetToken esTokenCSC_; - edm::ESGetToken esTokenGEM_; -}; - -// -// constructors and destructor -// -TestRotation::TestRotation(const edm::ParameterSet& iConfig) { - // Open root file and define tree - std::string fileName = iConfig.getUntrackedParameter("fileName", "test.root"); - theFile = new TFile(fileName.c_str(), "RECREATE"); - theTree = new TTree("theTree", "Detector units positions"); - - theTree->Branch("x", &x, "x/F"); - theTree->Branch("y", &y, "y/F"); - theTree->Branch("z", &z, "z/F"); - theTree->Branch("phi", &phi, "phi/F"); - theTree->Branch("theta", &theta, "theta/F"); - theTree->Branch("length", &length, "length/F"); - theTree->Branch("width", &width, "width/F"); - theTree->Branch("thick", &thick, "thick/F"); - dir = 0; - theTree->Branch("dir", "TRotMatrix", &dir); -} - -TestRotation::~TestRotation() { - theTree->Write(); - theFile->Close(); -} - -void TestRotation::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) { - edm::LogInfo("MuonAlignment") << "Starting!"; - - /* - // - // Build alignable muon geometry from event setup - // - - edm::ESTransientHandle cpv; - iRecord.getRecord().get( cpv ); - - DTGeometryBuilderFromDDD DTGeometryBuilder; - CSCGeometryBuilderFromDDD CSCGeometryBuilder; - - theDTGeometry = std::shared_ptr( DTGeometryBuilder.build( &(*cpv) ) ); - theCSCGeometry = std::shared_ptr( CSCGeometryBuilder.build( &(*cpv) ) ); - - AlignableMuon* theAlignableMuon = new AlignableMuon( &(*theDTGeometry) , &(*theCSCGeometry) ); -*/ - - // - // Retrieve ideal geometry from event setup - // - edm::ESHandle dtGeometry = iSetup.getHandle(esTokenDT_); - edm::ESHandle cscGeometry = iSetup.getHandle(esTokenCSC_); - edm::ESHandle gemGeometry = iSetup.getHandle(esTokenGEM_); - - AlignableMuon* theAlignableMuon = new AlignableMuon(&(*dtGeometry), &(*cscGeometry), &(*gemGeometry)); - - // Apply alignment - - const auto& theDTWheels = theAlignableMuon->DTWheels(); - std::cout << "Number of wheels=" << theDTWheels.size() << std::endl; - - for (const auto& iter : theDTWheels) { - std::cout << "------------------------" << std::endl << " BEFORE ROTATION " << std::endl; - - align::GlobalPoint pos_i = iter->globalPosition(); - align::RotationType dir_i = iter->globalRotation(); - - std::cout << "x=" << pos_i.x() << ", y=" << pos_i.y() << ", z=" << pos_i.z() << std::endl; - std::cout << "xx=" << dir_i.xx() << ", yx=" << dir_i.yx() << ", zx=" << dir_i.zx() << std::endl; - std::cout << "xy=" << dir_i.xy() << ", yy=" << dir_i.yy() << ", zy=" << dir_i.zy() << std::endl; - std::cout << "xz=" << dir_i.xz() << ", yz=" << dir_i.yz() << ", zz=" << dir_i.zz() << std::endl; - - // x = iter->surface().position().x(); - // y = iter->surface().position().y(); - // z = iter->surface().position().z(); - // std::cout << "X=" << x << ", Y= " << y << ", Z=" << z << std::endl ; - - double deltaPhi = 3.1415926 / 180 * 45; - - iter->rotateAroundGlobalZ(deltaPhi); - - std::cout << "------------------------" << std::endl << " AFTER ROTATION " << std::endl; - - align::GlobalPoint pos_f = iter->globalPosition(); - align::RotationType dir_f = iter->globalRotation(); - - std::cout << "x=" << pos_f.x() << ", y=" << pos_f.y() << ", z=" << pos_f.z() << std::endl; - std::cout << "xx=" << dir_f.xx() << ", yx=" << dir_f.yx() << ", zx=" << dir_f.zx() << std::endl; - std::cout << "xy=" << dir_f.xy() << ", yy=" << dir_f.yy() << ", zy=" << dir_f.zy() << std::endl; - std::cout << "xz=" << dir_f.xz() << ", yz=" << dir_f.yz() << ", zz=" << dir_f.zz() << std::endl; - std::cout << "------------------------" << std::endl; - } - - // delete theAlignableMuon ; - - edm::LogInfo("MuonAlignment") << "Done!"; -} - -//define this as a plug-in -DEFINE_FWK_MODULE(TestRotation); diff --git a/Alignment/MuonAlignment/test/TestTranslation.cpp b/Alignment/MuonAlignment/test/TestTranslation.cpp deleted file mode 100644 index 758d445308924..0000000000000 --- a/Alignment/MuonAlignment/test/TestTranslation.cpp +++ /dev/null @@ -1,155 +0,0 @@ -// -*- C++ -*- -// -// Package: TestTranslation -// Class: TestTranslation -// -// -// Description: Module to test the Alignment software -// -// - -// system include files -#include -#include -#include -#include - -// user include files -#include "FWCore/MessageLogger/interface/MessageLogger.h" -#include "FWCore/Framework/interface/EDAnalyzer.h" - -#include "FWCore/Framework/interface/EventSetup.h" -#include "FWCore/Framework/interface/ESTransientHandle.h" -#include "FWCore/Framework/interface/ESHandle.h" -#include "FWCore/Framework/interface/MakerMacros.h" - -#include "FWCore/ParameterSet/interface/ParameterSet.h" - -#include "Alignment/MuonAlignment/interface/AlignableMuon.h" -#include - -#include "Geometry/DTGeometry/interface/DTGeometry.h" - -#include "DataFormats/GeometrySurface/interface/Surface.h" - -#include - -// -// -// class declaration -// - -class TestTranslation : public edm::EDAnalyzer { -public: - explicit TestTranslation(const edm::ParameterSet&); - ~TestTranslation(); - - virtual void analyze(const edm::Event&, const edm::EventSetup&); - -private: - // ----------member data --------------------------- - TTree* theTree; - TFile* theFile; - float x, y, z, phi, theta, length, thick, width; - TRotMatrix* dir; - - typedef Surface::RotationType RotationType; - typedef Surface::PositionType PositionType; - - void apply(Alignable*); -}; - -// -// constructors and destructor -// -TestTranslation::TestTranslation(const edm::ParameterSet& iConfig) { - // Open root file and define tree - std::string fileName = iConfig.getUntrackedParameter("fileName", "test.root"); - theFile = new TFile(fileName.c_str(), "RECREATE"); - theTree = new TTree("theTree", "Detector units positions"); - - theTree->Branch("x", &x, "x/F"); - theTree->Branch("y", &y, "y/F"); - theTree->Branch("z", &z, "z/F"); - theTree->Branch("phi", &phi, "phi/F"); - theTree->Branch("theta", &theta, "theta/F"); - theTree->Branch("length", &length, "length/F"); - theTree->Branch("width", &width, "width/F"); - theTree->Branch("thick", &thick, "thick/F"); - dir = 0; - theTree->Branch("dir", "TRotMatrix", &dir); -} - -TestTranslation::~TestTranslation() { - theTree->Write(); - theFile->Close(); -} - -void TestTranslation::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) { - edm::LogInfo("MuonAlignment") << "Starting!"; - - /* - // - // Build alignable muon geometry from event setup - // - edm::ESTransientHandle cpv; - iRecord.getRecord().get( cpv ); - - DTGeometryBuilderFromDDD DTGeometryBuilder; - CSCGeometryBuilderFromDDD CSCGeometryBuilder; - - theDTGeometry = std::shared_ptr( DTGeometryBuilder.build( &(*cpv) ) ); - theCSCGeometry = std::shared_ptr( CSCGeometryBuilder.build( &(*cpv) ) ); - - AlignableMuon* theAlignableMuon = new AlignableMuon( &(*theDTGeometry) , &(*theCSCGeometry) ); -*/ - - // - // Retrieve ideal geometry from event setup - // - edm::ESHandle dtGeometry; - edm::ESHandle cscGeometry; - edm::ESHandle gemGeometry; - iSetup.get().get(dtGeometry); - iSetup.get().get(cscGeometry); - iSetup.get().get(gemGeometry); - - AlignableMuon* theAlignableMuon = new AlignableMuon(&(*dtGeometry), &(*cscGeometry), &(*gemGeometry)); - - // Apply alignment - for (const auto& iter : theAlignableMuon->DTChambers()) - apply(iter); - for (const auto& iter : theAlignableMuon->CSCEndcaps()) - apply(iter); - - edm::LogInfo("MuonAlignment") << "Done!"; -} - -void TestTranslation::apply(Alignable* it) { - std::cout << "------------------------" << std::endl << " BEFORE TRANSLATION " << std::endl; - - align::GlobalPoint pos_i = (it)->globalPosition(); - // RotationType dir_i = (it)->globalRotation(); - - std::cout << "x=" << pos_i.x() << ", y=" << pos_i.y() << ", z=" << pos_i.z() << std::endl; - - double dx = 1.0; - double dy = 2.0; - double dz = 3.0; - align::GlobalVector dr(dx, dy, dz); - it->move(dr); - - std::cout << "------------------------" << std::endl << " AFTER TRANSLATION " << std::endl; - - align::GlobalPoint pos_f = (it)->globalPosition(); - // RotationType dir_f = (it)->globalRotation(); - - std::cout << "x=" << pos_f.x() << ", y=" << pos_f.y() << ", z=" << pos_f.z() << std::endl; - - std::cout << "------------------------" << std::endl; -} - -// delete theAlignableMuon ; - -//define this as a plug-in -DEFINE_FWK_MODULE(TestTranslation); diff --git a/Alignment/MuonAlignment/test/computeIdealTransformation_cfg.py b/Alignment/MuonAlignment/test/computeIdealTransformation_cfg.py deleted file mode 100644 index ea4f76282d831..0000000000000 --- a/Alignment/MuonAlignment/test/computeIdealTransformation_cfg.py +++ /dev/null @@ -1,35 +0,0 @@ -import FWCore.ParameterSet.Config as cms - -process = cms.Process("COMPUTEIDEAL") -process.source = cms.Source("EmptySource") -process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(1)) - -process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi") -process.load("Geometry.MuonNumbering.muonNumberingInitialization_cfi") -process.load("Geometry.DTGeometry.dtGeometry_cfi") -process.load("Geometry.CSCGeometry.cscGeometry_cfi") - -# Ideal records -# frontier://FrontierProd/CMS_COND_31X_FROM21X DTAlignmentRcd DTIdealGeometry200_mc -# frontier://FrontierProd/CMS_COND_31X_FROM21X DTAlignmentErrorExtendedRcd DTIdealGeometryErrors200_mc -# frontier://FrontierProd/CMS_COND_31X_ALIGNMENT CSCAlignmentRcd CSCIdealGeometry310me42_mc -# frontier://FrontierProd/CMS_COND_31X_ALIGNMENT CSCAlignmentErrorExtendedRcd CSCIdealGeometryErrors310me42_mc -# frontier://FrontierProd/CMS_COND_31X_FROM21X GlobalPositionRcd IdealGeometry -# frontier://FrontierProd/CMS_COND_31X_FROM21X TrackerAlignmentRcd TrackerIdealGeometry210_mc -# frontier://FrontierProd/CMS_COND_31X_FROM21X TrackerAlignmentErrorExtendedRcd TrackerIdealGeometryErrors210_mc - -process.load("CondCore.DBCommon.CondDBSetup_cfi") -process.ideal31Xfrom21X = cms.ESSource("PoolDBESSource", - process.CondDBSetup, - connect = cms.string("frontier://FrontierProd/CMS_COND_31X_FROM21X"), - toGet = cms.VPSet(cms.PSet(record = cms.string("DTAlignmentRcd"), tag = cms.string("DTIdealGeometry200_mc")), - cms.PSet(record = cms.string("DTAlignmentErrorExtendedRcd"), tag = cms.string("DTIdealGeometryErrors200_mc")), - cms.PSet(record = cms.string("GlobalPositionRcd"), tag = cms.string("IdealGeometry")))) -process.ideal31X = cms.ESSource("PoolDBESSource", - process.CondDBSetup, - connect = cms.string("frontier://FrontierProd/CMS_COND_31X_ALIGNMENT"), - toGet = cms.VPSet(cms.PSet(record = cms.string("CSCAlignmentRcd"), tag = cms.string("CSCIdealGeometry310me42_mc")), - cms.PSet(record = cms.string("CSCAlignmentErrorExtendedRcd"), tag = cms.string("CSCIdealGeometryErrors310me42_mc")))) - -process.ComputeTransformation = cms.EDAnalyzer("ComputeTransformation", fileName = cms.string("Alignment/MuonAlignment/data/idealTransformation.py")) -process.Path = cms.Path(process.ComputeTransformation) diff --git a/Alignment/MuonAlignment/test/test-db-read_cfg.py b/Alignment/MuonAlignment/test/test-db-read_cfg.py deleted file mode 100644 index 48981077b91fe..0000000000000 --- a/Alignment/MuonAlignment/test/test-db-read_cfg.py +++ /dev/null @@ -1,67 +0,0 @@ -import FWCore.ParameterSet.Config as cms - -process = cms.Process("TEST") -# Message logger service -process.load("FWCore.MessageService.MessageLogger_cfi") - -# DT and CSC geometry -process.load("Geometry.MuonCommonData.muonIdealGeometryXML_cfi") -process.load("Geometry.MuonNumbering.muonNumberingInitialization_cfi") - -# Reading misalignments from DB -process.load("CondCore.DBCommon.CondDBSetup_cfi") - -process.maxEvents = cms.untracked.PSet( - input = cms.untracked.int32(1) -) -process.source = cms.Source("EmptySource") - - -process.DTGeometryMisalignedMuonProducer = cms.ESProducer("DTGeometryESModule", - appendToDataLabel = cms.string('idealForTestReader'), - applyAlignment = cms.bool(False), - alignmentsLabel = cms.string(''), - fromDDD = cms.bool(True) -) - -process.CSCGeometryMisalignedMuonProducer = cms.ESProducer("CSCGeometryESModule", - appendToDataLabel = cms.string('idealForTestReader'), - debugV = cms.untracked.bool(False), - useGangedStripsInME1a = cms.bool(False), - alignmentsLabel = cms.string(''), - useOnlyWiresInME1a = cms.bool(False), - useRealWireGeometry = cms.bool(True), - useCentreTIOffsets = cms.bool(False), - applyAlignment = cms.bool(False), - fromDDD = cms.bool(True), - fromDD4hep = cms.bool(False) -) - - -from CondCore.DBCommon.CondDBSetup_cfi import CondDBSetup -process.PoolDBESSource = cms.ESSource("PoolDBESSource", - process.CondDBSetup, - toGet = cms.VPSet(cms.PSet( - record = cms.string('DTAlignmentRcd'), - tag = cms.string('DT100InversepbScenario') - ), - cms.PSet( - record = cms.string('DTAlignmentErrorExtendedRcd'), - tag = cms.string('DT100InversepbScenarioErrors') - ), - cms.PSet( - record = cms.string('CSCAlignmentRcd'), - tag = cms.string('CSC100InversepbScenario') - ), - cms.PSet( - record = cms.string('CSCAlignmentErrorExtendedRcd'), - tag = cms.string('CSC100InversepbScenarioErrors') - )), - connect = cms.string('sqlite_file:Alignments.db') -) - -process.prod = cms.EDAnalyzer("TestMuonReader") - -process.p1 = cms.Path(process.prod) - - diff --git a/Alignment/MuonAlignment/test/test-hierarchy_cfg.py b/Alignment/MuonAlignment/test/test-hierarchy_cfg.py deleted file mode 100644 index bac775160ae86..0000000000000 --- a/Alignment/MuonAlignment/test/test-hierarchy_cfg.py +++ /dev/null @@ -1,64 +0,0 @@ -import FWCore.ParameterSet.Config as cms - -process = cms.Process("TEST") -# Message logger service -process.load("FWCore.MessageService.MessageLogger_cfi") - -# Ideal geometry producer -process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi") -process.load("Geometry.MuonNumbering.muonNumberingInitialization_cfi") - -# testing ideal geometry from Fake... -#process.load("Geometry.DTGeometry.dtGeometry_cfi") -#process.load("Geometry.CSCGeometry.cscGeometry_cfi") -#process.load("CalibMuon.Configuration.Muon_FakeAlignment_cff") - -## ...or testing geometry from DB/sqlite... -#process.load("Geometry.DTGeometry.dtGeometry_cfi") -#process.load("Geometry.CSCGeometry.cscGeometry_cfi") -#process.load("Alignment.CommonAlignmentProducer.GlobalPosition_Frontier_cff") -# -#from CondCore.DBCommon.CondDBSetup_cfi import * -#process.PoolDBESSource = cms.ESSource( -# "PoolDBESSource", -# CondDBSetup, -# #connect = cms.string('frontier://FrontierProd/CMS_COND_31X_ALIGNMENT'), -# connect = cms.string('sqlite_file:Alignments.db'), -# toGet = cms.VPSet(cms.PSet(record = cms.string('DTAlignmentRcd'), -# tag = cms.string('DT100InversepbScenario') -# ), -# cms.PSet(record = cms.string('DTAlignmentErrorExtendedRcd'), -# tag = cms.string('DT100InversepbScenarioErrors') -# ), -# cms.PSet(record = cms.string('CSCAlignmentRcd'), -# tag = cms.string('CSC100InversepbScenario') -# ), -# cms.PSet(record = cms.string('CSCAlignmentErrorExtendedRcd'), -# tag = cms.string('CSC100InversepbScenarioErrors') -# ) -# ) -# ) - -## ... or testing any of the scenarios from cff -process.load("Alignment.MuonAlignment.MisalignedMuon_cfi") -import Alignment.MuonAlignment.Scenarios_cff as _Scenarios -#process.MisalignedMuon.scenario = _Scenarios.MuonNoMovementsScenario -#process.MisalignedMuon.scenario = _Scenarios.Muon10InversepbScenario -process.MisalignedMuon.scenario = _Scenarios.Muon10inversePbScenario2008 - -process.maxEvents = cms.untracked.PSet( - input = cms.untracked.int32(1) -) -process.source = cms.Source("EmptySource") - -process.prod = cms.EDAnalyzer("TestMuonHierarchy") - -process.p1 = cms.Path(process.prod) -process.MessageLogger.cout = cms.untracked.PSet( - threshold = cms.untracked.string('INFO'), - default = cms.untracked.PSet( - limit = cms.untracked.int32(10000000) - ) -) - - diff --git a/Alignment/MuonAlignment/test/test-rotation_cfg.py b/Alignment/MuonAlignment/test/test-rotation_cfg.py deleted file mode 100644 index d258cc38b7fc3..0000000000000 --- a/Alignment/MuonAlignment/test/test-rotation_cfg.py +++ /dev/null @@ -1,41 +0,0 @@ -import FWCore.ParameterSet.Config as cms - -process = cms.Process("TEST") -process.load("Geometry.MuonCommonData.muonIdealGeometryXML_cfi") - -process.load("Geometry.MuonNumbering.muonNumberingInitialization_cfi") - -# DT geometry -# include "Geometry/DTGeometry/data/dtGeometry.cfi" -# CSC geometry -# include "Geometry/CSCGeometry/data/cscGeometry.cfi" -# Misalignment example scenario producer -process.load("Alignment.MuonAlignment.MisalignedMuon_cfi") - -process.MessageLogger = cms.Service("MessageLogger", - info_txt = cms.untracked.PSet( - threshold = cms.untracked.string('INFO') - ), - cerr = cms.untracked.PSet( - threshold = cms.untracked.string('ERROR') - ), - destinations = cms.untracked.vstring('info_txt', - 'cerr') -) - -process.maxEvents = cms.untracked.PSet( - input = cms.untracked.int32(1) -) -process.source = cms.Source("EmptySource") - -process.MisalignedTracker = cms.ESProducer("MisalignedTrackerESProducer", - process.MuonNoMovementsScenario -) - -process.prod = cms.EDAnalyzer("TestRotation", - fileName = cms.untracked.string('misaligned-2.root') -) - -process.p1 = cms.Path(process.prod) - - diff --git a/Alignment/MuonAlignment/test/test-translation_cfg.py b/Alignment/MuonAlignment/test/test-translation_cfg.py deleted file mode 100644 index 78cf7b7ae491e..0000000000000 --- a/Alignment/MuonAlignment/test/test-translation_cfg.py +++ /dev/null @@ -1,41 +0,0 @@ -import FWCore.ParameterSet.Config as cms - -process = cms.Process("TEST") -process.load("Geometry.MuonCommonData.muonIdealGeometryXML_cfi") - -process.load("Geometry.MuonNumbering.muonNumberingInitialization_cfi") - -# DT geometry -# include "Geometry/DTGeometry/data/dtGeometry.cfi" -# -# CSC geometry -# include "Geometry/CSCGeometry/data/cscGeometry.cfi" -process.load("Alignment.MuonAlignment.MisalignedMuon_cfi") - -process.MessageLogger = cms.Service("MessageLogger", - info_txt = cms.untracked.PSet( - threshold = cms.untracked.string('INFO') - ), - cerr = cms.untracked.PSet( - threshold = cms.untracked.string('ERROR') - ), - destinations = cms.untracked.vstring('info_txt', - 'cerr') -) - -process.maxEvents = cms.untracked.PSet( - input = cms.untracked.int32(1) -) -process.source = cms.Source("EmptySource") - -process.MisalignedTracker = cms.ESProducer("MisalignedTrackerESProducer", - process.MuonNoMovementsScenario -) - -process.prod = cms.EDAnalyzer("TestTranslation", - fileName = cms.untracked.string('misaligned-2.root') -) - -process.p1 = cms.Path(process.prod) - -