geometryDiff Namespace Reference

Functions
def	loopover (which)
def	matrixmult (a, b)
def	rotFromEuler (g)
def	rotFromPhi (g)
def	transpose (a)

Variables
	args
int	epsilon = 1e-6
	file
	geom1 = MuonGeometry(file(args[0]))
	geom2 = MuonGeometry(file(args[1]))
	opts = dict(opts)
	sqrtepsilon = sqrt(epsilon)
string	usage = "Usage: geometryDiff.py [-h|--help] [-e|--epsilon epsilon] geometry1.xml geometry2.xml"

Function Documentation

loopover()

def geometryDiff.loopover

(

which

)

Definition at line 62 of file geometryDiff.py.

def loopover(which):
  if which == "DT":
    keys = geom1.dt.keys()
    keys.sort(dtorder)
 
  elif which == "CSC":
    keys = geom1.csc.keys()
    keys.sort(cscorder)
 
  else: raise Exception
 
  for key in keys:
    if which == "DT":
      g1 = geom1.dt[key]
      g2 = geom2.dt[key]
    else:
      g1 = geom1.csc[key]
      g2 = geom2.csc[key]
 
    if g1.relativeto != g2.relativeto:
      print("%s %s relativeto=\"%s\" versus relativeto=\"%s\"" % (which, str(key), g1.relativeto, g2.relativeto))
 
    if abs(g1.x - g2.x) > epsilon or abs(g1.y - g2.y) > epsilon or abs(g1.z - g2.z) > epsilon:
      print("%s %s position difference: (%g, %g, %g) - (%g, %g, %g) = (%g, %g, %g)" % \
            (which, str(key), g1.x, g1.y, g1.z, g2.x, g2.y, g2.z, g1.x - g2.x, g1.y - g2.y, g1.z - g2.z))
 
    if "phix" in g1.__dict__:
      g1type = "phi"
      g1a, g1b, g1c = g1.phix, g1.phiy, g1.phiz
      g1rot = rotFromPhi(g1)
    else:
      g1type = "euler"
      g1a, g1b, g1c = g1.alpha, g1.beta, g1.gamma
      g1rot = rotFromEuler(g1)
 
    if "phix" in g2.__dict__:
      g2type = "phi"
      g2a, g2b, g2c = g2.phix, g2.phiy, g2.phiz
      g2rot = rotFromPhi(g2)
    else:
      g2type = "euler"
      g2a, g2b, g2c = g2.alpha, g2.beta, g2.gamma
      g2rot = rotFromEuler(g2)
    
    diff = matrixmult(g1rot, transpose(g2rot))
    if abs(diff[0][0] - 1.) > sqrtepsilon or abs(diff[1][1] - 1.) > sqrtepsilon or abs(diff[2][2] - 1.) > sqrtepsilon or \
       abs(diff[0][1]) > epsilon or abs(diff[0][2]) > epsilon or abs(diff[1][2]) > epsilon:
      print("%s %s rotation difference: %s(%g, %g, %g) - %s(%g, %g, %g) = %s" % \
            (which, str(key), g1type, g1a, g1b, g1c, g2type, g2a, g2b, g2c, str(diff)))
 
loopover("DT")
loopover("CSC")
 
 

References funct.abs(), matrixmult(), print(), rotFromEuler(), rotFromPhi(), str, and transpose().

matrixmult()

def geometryDiff.matrixmult	(		a,
			b
	)

Definition at line 36 of file geometryDiff.py.

def matrixmult(a, b):
  return [[sum([i*j for i, j in zip(row, col)]) for col in zip(*b)] for row in a]
 

References ComparisonHelper.zip().

Referenced by loopover(), and rotFromPhi().

rotFromEuler()

def geometryDiff.rotFromEuler

(

g

)

Definition at line 55 of file geometryDiff.py.

def rotFromEuler(g):
  s1, s2, s3 = sin(g.alpha), sin(g.beta), sin(g.gamma)
  c1, c2, c3 = cos(g.alpha), cos(g.beta), cos(g.gamma)
  return [[c2 * c3,    c1 * s3 + s1 * s2 * c3,   s1 * s3 - c1 * s2 * c3,],
          [-c2 * s3,   c1 * c3 - s1 * s2 * s3,   s1 * c3 + c1 * s2 * s3,],
          [s2,        -s1 * c2,                  c1 * c2,               ]]
 

References funct.cos(), and funct.sin().

Referenced by loopover().

rotFromPhi()

def geometryDiff.rotFromPhi

(

g

)

Definition at line 42 of file geometryDiff.py.

def rotFromPhi(g):
  phix, phiy, phiz = g.phix, g.phiy, g.phiz
  rotX = [[1.,         0.,         0.,       ],
          [0.,         cos(phix),  sin(phix),],
          [0.,        -sin(phix),  cos(phix),]]
  rotY = [[cos(phiy),  0.,        -sin(phiy),],
          [0.,         1.,         0.,       ],
          [sin(phiy),  0.,         cos(phiy),]]
  rotZ = [[cos(phiz),  sin(phiz),  0.,       ],
          [-sin(phiz), cos(phiz),  0.,       ],
          [0.,         0.,         1.,       ]]
  return matrixmult(rotX, matrixmult(rotY, rotZ))
 

References funct.cos(), matrixmult(), and funct.sin().

Referenced by loopover().

transpose()

def geometryDiff.transpose

(

a

)

Definition at line 39 of file geometryDiff.py.

def transpose(a):
  return [[a[j][i] for j in range(len(a[i]))] for i in range(len(a))]
 

References FastTimerService_cff.range.

Referenced by AlignmentCorrelationsStore.correlations(), AlignmentExtendedCorrelationsStore.correlations(), AlignmentCorrelationsStore.correlationsAvailable(), AlignmentExtendedCorrelationsStore.correlationsAvailable(), AlignmentCorrelationsStore.fillCorrelationsTable(), AlignmentExtendedCorrelationsStore.fillCorrelationsTable(), AlignmentExtendedCorrelationsStore.getCorrelations(), loopover(), AlignmentCorrelationsStore.setCorrelations(), AlignmentExtendedCorrelationsStore.setCorrelations(), and npstat::ArrayND< Num1, Len1, Dim1 >.transpose().

Variable Documentation

args

geometryDiff.args

Definition at line 11 of file geometryDiff.py.

epsilon

geometryDiff.epsilon = 1e-6

Definition at line 26 of file geometryDiff.py.

Referenced by ReferenceTrajectory.addMaterialEffectsCurvlinGbl(), ReferenceTrajectory.addMaterialEffectsLocalGbl(), trklet::VarAdjustK.adjust(), almost_equal(), angle_units::operators.almostEqual(), GlobalTrackerMuonAlignment.analyzeTrackTrack(), GlobalTrackerMuonAlignment.analyzeTrackTrajectory(), trklet::TrackletCalculatorDisplaced.approxprojdisk(), GeneralBinFinderInPhi< T >.binIndex(), magneticfield::bSlab.bSlab(), TrackWithHistory.checkAtEnd(), TauDiscriminantCutMultiplexerT< TauType, TauTypeRef, ParentClass >.discriminate(), CMSCGEN.dNdEmudEnu(), IntegralOverPhiFunction.DoEval(), IntegrandThetaFunction.DoEval(), DTSegment4DQuality.dqmAnalyze(), TMultiDimFet.EvalControl(), KinematicPerigeeConversions.extendedPerigeeFromKinematicParameters(), DAClusterizerInZ_vect.fill(), DAClusterizerInZT_vect.fill(), Thrust.finalAxis(), PerigeeConversions.ftsToPerigeeParameters(), funct::GaussLegendreIntegrator.GaussLegendreIntegrator(), GEMEfficiencyAnalyzer.GEMEfficiencyAnalyzer(), PixelInactiveAreaFinder.getPhiSpanBarrel(), PixelInactiveAreaFinder.getPhiSpanEndcap(), ConversionProducer.getTrackImpactPosition(), AlignmentParameterStore.hierarchyConstraints(), SmartPropagator.initTkVolume(), HiPuRhoProducer.inputTowers(), VirtualJetProducer.inputTowers(), CSCLayerGeometry.inside(), SimpleNavigationSchool.linkNextLargerLayer(), SimpleNavigationSchool.linkNextLayerInGroup(), angle0to2pi.make0To2pi(), DTTopology.onWhichBorder(), DTTopology.onWhichBorder_old(), CMSInsideOutAlgorithm::ListIteratorLesserByDeltaR.operator()(), NumericSafeLessByEt< T >.operator()(), NumericSafeLessByPt< T >.operator()(), NumericSafeGreaterByEt< T >.operator()(), NumericSafeGreaterByPt< T >.operator()(), DTSegmentCand.operator==(), SimpleNavigationSchool.outerRadiusIncrease(), placeAssembly(), RecHitPropagator.propagate(), SteppingHelixPropagator.propagate(), tauImpactParameter::ErrorMatrixPropagator.propagateError(), SiLinearChargeDivider.readPulseShape(), FixedFlatOccupancy1d.set_epsilon_max(), CSC01.set_epsilon_max(), FixedFlatOccupancy1d.set_epsilon_min(), NAEquationSolver< Function >.SetTolerance(), HLLHCEvtVtxGenerator.sigma(), SmartPropagator.SmartPropagator(), DAClusterizerInZ.split(), DAClusterizerInZ_vect.split(), DAClusterizerInZ.splitAll(), muonisolation::NominalEfficiencyThresholds.thresholdValueForEfficiency(), trklet::VarAdjustK.VarAdjustK(), and trklet::VarAdjustKR.VarAdjustKR().

file

geometryDiff.file

Definition at line 13 of file geometryDiff.py.

Referenced by cms::cuda.abortOnCudaError(), ConfigBuilder.ConfigBuilder.addCustomise(), addFilesWithFork(), CmsShowNavigator.addFilter(), SiStripGainFromData.algoEndJob(), DumpFWRecoGeometry.analyze(), DumpFWTGeoRecoGeometry.analyze(), CSCFileDumper.analyze(), L1GctPrintLuts.analyze(), FFTJetTreeDump.analyze(), L1GtVhdlTemplateFile.append(), CmsShowNavigator.appendFile(), BeamProfileVtxGenerator.BeamProfileVtxGenerator(), ExternalLHEProducer.beginRunProduce(), ESPedestalTask.bookHistograms(), fwlite::BranchMapReader.BranchMapReader(), fwlite::internal::BranchMapReaderStrategyV1.BranchMapReaderStrategyV1(), fwlite::internal::BranchMapReaderStrategyV11.BranchMapReaderStrategyV11(), fwlite::internal::BranchMapReaderStrategyV17.BranchMapReaderStrategyV17(), fwlite::internal::BranchMapReaderStrategyV7.BranchMapReaderStrategyV7(), fwlite::internal::BranchMapReaderStrategyV8.BranchMapReaderStrategyV8(), CTPPSCompositeESSource.buildDirectSimuData(), SiStripPsuDetIdMap.BuildMap(), CmsShowNavigator.canEditFiltersExternally(), FWTEveViewer.CaptureAndSaveImage(), eostools.cat(), CDFRunInfo.CDFRunInfo(), CmsShowNavigator.changeFilter(), check_file(), check_files(), dqm::DQMFileSaverOnline.checkError(), QualityTester.configureTests(), SiPixelLorentzAngleCalibration.createFromTree(), SiStripLorentzAngleCalibration.createFromTree(), SiStripBackplaneCalibration.createFromTree(), cms::cuda.cudaCheck_(), hitfit::Defaults_Textrep.Defaults_Textrep(), TrackerDpgAnalysis.delay(), XrdAdaptor::Source.determineHostExcludeString(), FWLiteJobMetadataManager.doUpdate(), DTDigiReader.DTDigiReader(), LASBarrelAlgorithm.Dump(), ProvenanceDumper.dumpEventFilteringParameterSets_(), GeometryConfiguration.dumpFileList(), FIPConfiguration.dumpFileList(), CastorDumpConditions.dumpIt(), RunManagerMTWorker.DumpMagneticField(), edm::DuplicateTreeSentry.dup(), sistrip::EnsembleCalibrationLA.endJob(), QualityCutsAnalyzer.endJob(), ESDataFormatterV4.ESDataFormatterV4(), ESDigiToRaw.ESDigiToRaw(), ESElectronicsMapper.ESElectronicsMapper(), ESIntegrityClient.ESIntegrityClient(), ESIntegrityTask.ESIntegrityTask(), ESPedestalClient.ESPedestalClient(), ESUnpacker.ESUnpacker(), DTExpandMap.expandSteering(), FWGUIManager.exportAllViews(), reco::ExpressionEvaluator.ExpressionEvaluator(), JsonOutputProducer::JsonConfigurationBlock.extractFileBase(), spf::SherpackFetcher.Fetch(), fftjetcms.fftjet_Grid2d_parser(), fftjetcms.fftjet_LinearInterpolator2d_parser(), fftjetcms.fftjet_MembershipFunction_parser(), fftjetcms.fftjet_PeakSelector_parser(), CmsShowMain.fileChangedSlot(), dqm::DQMFileSaverBase.filename(), edm::FileRandomKEThetaGunProducer.FileRandomKEThetaGunProducer(), edm::FileRandomMultiParticlePGunProducer.FileRandomMultiParticlePGunProducer(), SiPixelCPEGenericErrorParm.fillCPEGenericErrorParm(), TkHistoMap.fillFromAscii(), HGCalTriggerGeometryV9Imp1.fillNeighborMap(), HGCalTriggerGeometryHexLayerBasedImp1.fillNeighborMaps(), AdaptiveVertexFitter.fit(), PedeSteerer.fixParameter(), trklet::FPGAWord.FPGAWord(), EcalLaserCondTools.from_hdf_to_db(), getRunInfo.Func_FillInfoHlt(), getRunInfo.Func_FillInfoMagnetHistory(), getRunInfo.Func_FillInfoRunSummary(), ExternalLHEProducer.generateRunInfo(), getCentralityFromFile(), PedeSteererWeakModeConstraints.getFile(), fwlite::BranchMapReader.getFileVersion(), getHists(), L1MuBMTEtaPatternLut.getIgnoredLines(), L1MuBMTQualPatternLut.getIgnoredLines(), production_tasks.MonitorJobs.getjobid(), findQualityFiles.getJSONGoodRuns(), CmsShowNavigator.getNSelectedEvents(), CmsShowNavigator.getNTotalEvents(), getplot(), ZIterativeAlgorithmWithFit.getStatWeights(), CmsShowNavigator.goToRunEvent(), FWConfigurationManager.guessAndReadFromFile(), GzInputStream.GzInputStream(), h2DSegm.h2DSegm(), h4DSegm.h4DSegm(), HCalSD.HCalSD(), HDelta.HDelta(), PedeSteerer.hierarchyConstraint(), HLikelihoodVSPart.HLikelihoodVSPart(), HMassVSPart.HMassVSPart(), HMassVSPartProfile.HMassVSPartProfile(), HResolutionVSPart.HResolutionVSPart(), hTMaxCell.hTMaxCell(), PixelToLNKAssociateFromAscii.init(), PixelToFEDAssociateFromAscii.init(), fit::RootMinuitCommands< Function >.init(), HGCalTriggerClusterIdentificationBDT.initialize(), ElectronEnergyRegressionEvaluate.initialize(), EGammaMvaEleEstimatorCSA14.initialize(), EGammaMvaEleEstimator.initialize(), CSCOverlapsAlignmentAlgorithm.initialize(), XrdAdaptor::RequestManager.initialize(), dqm::DQMFileSaverPB.initRun(), XrdAdaptor::Source.isDCachePool(), CmsShowNavigator.isFirstEvent(), CmsShowNavigator.isLastEvent(), JetCalibrationParameterSetTauJet.JetCalibrationParameterSetTauJet(), JetPartonNamespace::JetPartonCalibrationParameterSet.JetPartonCalibrationParameterSet(), L1MuDTPtaLut.load(), L1MuDTPhiLut.load(), L1MuDTEtaPatternLut.load(), L1MuBMTEtaPatternLut.load(), L1MuDTQualPatternLut.load(), L1MuBMTQualPatternLut.load(), L1MuDTExtLut.load(), svgfig.load(), L1TMuonBarrelParamsHelper.load_ext(), L1TMuonBarrelParamsHelper.load_phi(), L1TMuonBarrelParamsHelper.load_pt(), cmsHarvester.load_ref_hist_mappings(), FWGeometry.loadMap(), PPSDirectSimulationData.loadObject(), ClusterShapeHitFilter.loadPixelLimits(), edm::LogDebug_.LogDebug_(), logger.logger.logJobs(), motionPolicyChamber.loopover(), ls_cert_type(), main(), BeamSpotWorkflow.main(), jpt::Map.Map(), sistrip::MeasureLA.MeasureLA(), mergeLUTs(), multibsvspvplots(), MVAVariableManager< reco::GsfElectron >.MVAVariableManager(), fwlite::BranchMapReader.newStrategy(), nlumis(), TFWLiteSelectorBasic.Notify(), DavixStorageMaker.open(), LocalStorageMaker.open(), XrdFile.open(), DCacheStorageMaker.open(), XrdStorageMaker.open(), StormStorageMaker.open(), StormLcgGtStorageMaker.open(), CmsShowNavigator.openFile(), trklet.openfile(), l1t::L1TCaloStage1LutWriter.openOutputFile(), gen::Pythia6Service.openPYUPDA(), gen::Pythia6Service.openSLHA(), perftools::EdmEventSize.parseFile(), CTPPSPixelDAQMappingESSourceXML.ParseXML(), TotemDAQMappingESSourceXML.ParseXML(), IODConfig.populateClob(), PPSPixelDigiAnalyzer.PPSPixelDigiAnalyzer(), CTPPSRPAlignmentCorrectionsDataESSourceXMLCommon.PrepareSequence(), print_rates(), printSOF(), processDirectory(), FFTJetTreeDump.processTreeData(), SiPixelFakeTemplateDBObjectESSource.produce(), SiPixelFakeGenErrorDBObjectESSource.produce(), perftools::EdmEventSize.produceHistos(), edmIntegrityCheck.PublishToFileSystem.publish(), SiPixelGenError.pushfile(), SiPixelTemplate2D.pushfile(), SiPixelTemplate.pushfile(), Histos.put(), edmIntegrityCheck.PublishToFileSystem.read(), MuonResidualsTwoBin.read(), MuonResidualsFitter.read(), hitfit::Defaults_Textrep.read_file(), edmplugin.readCacheFile(), PhysicsTools::MVAComputer.readCalibration(), DTSurvey.ReadChambers(), DTCalibrationMap.readConsts(), HcalLaserEventFilter2012.readEventListFile(), DQMRootSource.readFile_(), L1MuCSCPtLut.readFromDBS(), l1tpf_impl.readFromFile(), reco::details.readGzipFile(), evf::EvFDaqDirector.readLastLSEntry(), CSCTFConfigProducer.readLUT(), l1tpf_impl.readManyFromFile(), LASAlignmentTubeAlgorithm.ReadMisalignmentFromFile(), LASBarrelAlgorithm.ReadMisalignmentFromFile(), FedRawDataInputSource.readNextChunkIntoBuffer(), MuScleFitBase.readProbabilityDistributionsFromFile(), EnergyScaleCorrection.readScalesFromFile(), EnergyScaleCorrection.readSmearingsFromFile(), LASBarrelAlgorithm.ReadStartParametersFromFile(), MuonAlignmentFromReference.readTmpFiles(), RootTreeHandler.readTree(), FedRawDataInputSource.readWorker(), CmsShowNavigator.removeFilter(), Phase1PixelMaps.retrieveCorners(), AdaptiveVertexFitter.reWeightTracks(), production_tasks.SourceCFG.run(), production_tasks.FullCFG.run(), production_tasks.RunTestEvents.run(), production_tasks.ExpandConfig.run(), production_tasks.MonitorJobs.run(), production_tasks.CheckJobStatus.run(), production_tasks.WriteJobReport.run(), LegacyIOHelper.save(), SiPixelPhase1Analyzer.SaveDetectorData(), FWTEveViewer.SaveJpg(), FWTEveViewer.SavePng(), edm::service::RandomNumberGeneratorService.saveStatesToFile(), JME::JetResolutionObject.saveToFile(), SendMonitoringInfo(), trklet::FPGAWord.set(), EcalUncalibRecHitWorkerFixedAlphaBetaFit.setAlphaBeta(), CommonAnalyzer.setFile(), GenericBenchmark.setfile(), CmsShowNavigator.setFrom(), FWGeometryTableViewManager.setGeoManagerFromFile(), gen::Pythia6Service.setPYUPDAParams(), gen::Pythia6Service.setSLHAParams(), GenericBenchmark.setup(), cond::CredentialStore.setUpForService(), alcaDQMUpload.startUpload(), edm::LogDebug_.stripLeadingDirectoryTree(), MillePedeAlignmentAlgorithm.terminate(), edm::roottree.trainCache(), HcaluLUTTPGCoder.update(), fwlite::internal::Strategy.updateFile(), fwlite::BranchMapReader.updateFile(), fwlite::internal::BranchMapReaderStrategyV1.updateFile(), fwlite::internal::BranchMapReaderStrategyV8.updateFile(), fwlite::internal::BranchMapReaderStrategyV11.updateFile(), fwlite::internal::BranchMapReaderStrategyV17.updateFile(), CmsShowNavigator.updateFileFilters(), MultiVertexFitter.updateSeeds(), CmsShowNavigator.updateSelectorsInfo(), fileCollector2.uploadFile(), MultiVertexFitter.vertices(), AdaptiveVertexFitter.weightTracks(), XrdFile.write(), MuonResidualsTwoBin.write(), MuonResidualsFitter.write(), sistrip::EnsembleCalibrationLA.write_calibrations(), MCScenario_CRAFT1_22X.write_check_cfg(), MCScenario_CRAFT1_22X.write_conversion_cfg(), cmsHarvester.write_crab_config(), sistrip::EnsembleCalibrationLA.write_ensembles_plots(), sistrip::EnsembleCalibrationLA.write_ensembles_text(), cmsHarvester.write_harvesting_config(), cmsHarvester.write_me_extraction_config(), cmsHarvester.write_multicrab_config(), sistrip::MeasureLA.write_report_plots(), sistrip::MeasureLA.write_report_text(), sistrip::MeasureLA.write_report_text_ms(), sistrip::EnsembleCalibrationLA.write_samples_plots(), MCScenario_CRAFT1_22X.write_xml(), createJobs.writeAlignCfg(), PhysicsTools::MVAComputer.writeCalibration(), FWPartialConfigSaveGUI.WriteConfig(), createJobs.writeGatherCfg(), HLTriggerJSONMonitoring.writeIniFile(), L1TriggerJSONMonitoring.writeIniFile(), l1t::L1TCaloStage1LutWriter.writeIsoTauLut(), HLTriggerJSONMonitoring.writeJsdFile(), L1TriggerJSONMonitoring.writeJsdFile(), l1tpf_impl.writeManyToFile(), MuonAlignmentFromReference.writeTmpFiles(), l1tpf_impl.writeToFile(), FWConfigurationManager.writeToFile(), Fp420AnalysisHistManager.WriteToFile(), BscAnalysisHistManager.WriteToFile(), edmtest::HcalDumpConditions.writeToFile(), SiPixelLorentzAngleCalibration.writeTree(), SiStripLorentzAngleCalibration.writeTree(), SiStripBackplaneCalibration.writeTree(), createJobs.writeValidationCfg(), EcalDccWeightBuilder.writeWeightToAsciiFile(), EcalDccWeightBuilder.writeWeightToRootFile(), ZSPJPTJetCorrector.ZSPJPTJetCorrector(), AnalysisTasksAnalyzerJEC.~AnalysisTasksAnalyzerJEC(), DQMRootSource.~DQMRootSource(), and DTDigiReader.~DTDigiReader().

geom1

geometryDiff.geom1 = MuonGeometry(file(args[0]))

Definition at line 30 of file geometryDiff.py.

geom2

geometryDiff.geom2 = MuonGeometry(file(args[1]))

Definition at line 31 of file geometryDiff.py.

opts

geometryDiff.opts = dict(opts)

Definition at line 11 of file geometryDiff.py.

Referenced by edm::StreamSchedule.initializeEarlyDelete(), edm::MallocOpts.operator!=(), edm.operator<<(), edm::MallocOpts.operator==(), edm::StreamSchedule.StreamSchedule(), and edm::service::TriggerNamesService.TriggerNamesService().

sqrtepsilon

geometryDiff.sqrtepsilon = sqrt(epsilon)

Definition at line 34 of file geometryDiff.py.

usage

geometryDiff.usage = "Usage: geometryDiff.py [-h|--help] [-e|--epsilon epsilon] geometry1.xml geometry2.xml"

Definition at line 8 of file geometryDiff.py.