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.

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

 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")
 
 

◆ matrixmult()

def geometryDiff.matrixmult	(	a,
		b
	)

Definition at line 36 of file geometryDiff.py.

References ComparisonHelper.zip().

Referenced by loopover(), and rotFromPhi().

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

◆ rotFromEuler()

def geometryDiff.rotFromEuler ( g )

Definition at line 55 of file geometryDiff.py.

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

Referenced by loopover().

 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,               ]]
 

◆ rotFromPhi()

def geometryDiff.rotFromPhi ( g )

Definition at line 42 of file geometryDiff.py.

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

Referenced by loopover().

 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))
 

◆ transpose()

def geometryDiff.transpose ( a )

Definition at line 39 of file geometryDiff.py.

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< Numeric >.transpose().

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

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(), BoostedDoubleSVProducer.calcNsubjettiness(), 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(), 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(), TemplatedSecondaryVertexProducer< IPTI, VTX >.produce(), RecHitPropagator.propagate(), SteppingHelixPropagator.propagate(), SiLinearChargeDivider.readPulseShape(), scaleGraphByLuminosity(), FixedFlatOccupancy1d.set_epsilon_max(), CSC01.set_epsilon_max(), FixedFlatOccupancy1d.set_epsilon_min(), SiPixelRecHitQuality::Packing.setProbabilityQ(), SiPixelRecHitQuality::Packing.setProbabilityXY(), 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.

◆ 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::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.

Functions

Variables