Functions
def	cholesky (A)

def	mmult (a, b)
	Print out user's choices as diagnostics. More...

def	mtrans (a)

def	mvdot (m, v)

def	random6dof ()
	Generate correlated random misalignments for all chambers. More...

Variables
	action

	args

	ave_phix = sum_phix/float(len(misal))

	ave_phixphix = sum_phixphix/float(len(misal))

	ave_phixphiy = sum_phixphiy/float(len(misal))

	ave_phixphiz = sum_phixphiz/float(len(misal))

	ave_phiy = sum_phiy/float(len(misal))

	ave_phiyphiy = sum_phiyphiy/float(len(misal))

	ave_phiyphiz = sum_phiyphiz/float(len(misal))

	ave_phiz = sum_phiz/float(len(misal))

	ave_phizphiz = sum_phizphiz/float(len(misal))

	ave_x = sum_x/float(len(misal))

	ave_xphix = sum_xphix/float(len(misal))

	ave_xphiy = sum_xphiy/float(len(misal))

	ave_xphiz = sum_xphiz/float(len(misal))

	ave_xx = sum_xx/float(len(misal))

	ave_xy = sum_xy/float(len(misal))

	ave_xz = sum_xz/float(len(misal))

	ave_y = sum_y/float(len(misal))

	ave_yphix = sum_yphix/float(len(misal))

	ave_yphiy = sum_yphiy/float(len(misal))

	ave_yphiz = sum_yphiz/float(len(misal))

	ave_yy = sum_yy/float(len(misal))

	ave_yz = sum_yz/float(len(misal))

	ave_z = sum_z/float(len(misal))

	ave_zphix = sum_zphix/float(len(misal))

	ave_zphiy = sum_zphiy/float(len(misal))

	ave_zphiz = sum_zphiz/float(len(misal))

	ave_zz = sum_zz/float(len(misal))

	cfgfile = file(outputName + "_convert_cfg.py", "w")
	Convert it to an SQLite file with CMSSW. More...

	chomat = cholesky(matrix)

string	components = "xx"

	default

	dest

	endcap

	globalape = mmult(rot, mmult(localape, mtrans(rot)))

	globalxx = globalape[0][0]

	globalxy = globalape[0][1]

	globalxz = globalape[0][2]

	globalyy = globalape[1][1]

	globalyz = globalape[1][2]

	globalzz = globalape[2][2]

	help

list	localape = [[xx, xy, xz], [xy, yy, yz], [xz, yz, zz]]

list	matrix

dictionary	misal = {}

	options

	outputName = args[0]

	parser = OptionParser(usage="Usage: python %prog outputName [options] (default is unit matrix times 1e-15)")
	Get variances and covariances from the commandline. More...

	rot = rotation[system, whendcap, station, ring, sector]

int	sum_phix = 0

int	sum_phixphix = 0

int	sum_phixphiy = 0

int	sum_phixphiz = 0

int	sum_phiy = 0

int	sum_phiyphiy = 0

int	sum_phiyphiz = 0

int	sum_phiz = 0

int	sum_phizphiz = 0

int	sum_x = 0
	More diagnostics. More...

int	sum_xphix = 0

int	sum_xphiy = 0

int	sum_xphiz = 0

int	sum_xx = 0

int	sum_xy = 0

int	sum_xz = 0

int	sum_y = 0

int	sum_yphix = 0

int	sum_yphiy = 0

int	sum_yphiz = 0

int	sum_yy = 0

int	sum_yz = 0

int	sum_z = 0

int	sum_zphix = 0

int	sum_zphiy = 0

int	sum_zphiz = 0

int	sum_zz = 0

	txtfile = file(outputName + "_correlations.txt", "w")
	Delete all three files at once to make sure the user never sees stale data (e.g. More...

	wheel

	xmlfile = file(outputName + ".xml", "w")
	Make an XML representation of the misalignment. More...

Function Documentation

def makeMuonMisalignmentScenario.cholesky ( A )

Cholesky decomposition of the correlation matrix to properly normalize the transformed random deviates

Definition at line 111 of file makeMuonMisalignmentScenario.py.

References mmult(), and mathSSE.sqrt().

 def cholesky(A):
     """Cholesky decomposition of the correlation matrix to properly normalize the transformed random deviates"""
 
     # A = L * D * L^T = (L * D^0.5) * (L * D^0.5)^T where we want (L * D^0.5), the "square root" of A
     # find L and D from A using recurrence relations
     L = {}
     D = {}
     for j in range(len(A[0])):
         D[j] = A[j][j] - sum([L[j,k]**2 * D[k] for k in range(j)])
         for i in range(len(A)):
             if i > j:
                 L[i,j] = (A[i][j] - sum([L[i,k] * L[j,k] * D[k] for k in range(j)])) / D[j]
 
     L = [[    1.,     0.,     0.,     0.,     0., 0.],
          [L[1,0],     1.,     0.,     0.,     0., 0.],
          [L[2,0], L[2,1],     1.,     0.,     0., 0.],
          [L[3,0], L[3,1], L[3,2],     1.,     0., 0.],
          [L[4,0], L[4,1], L[4,2], L[4,1],     1., 0.],
          [L[5,0], L[5,1], L[5,2], L[5,1], L[5,0], 1.]]
 
     Dsqrt = [[sqrt(D[0]),         0.,          0.,         0.,         0.,         0.],
              [        0., sqrt(D[1]),          0.,         0.,         0.,         0.],
              [        0.,         0., sqrt(D[2]),          0.,         0.,         0.],
              [        0.,         0.,          0., sqrt(D[3]),         0.,         0.],
              [        0.,         0.,          0.,         0., sqrt(D[4]),         0.],
              [        0.,         0.,          0.,         0.,         0., sqrt(D[5])]]
 
     return mmult(L, Dsqrt)
 

def makeMuonMisalignmentScenario.mmult	(	a,
		b
	)

Print out user's choices as diagnostics.

Some useful mathematical transformations (why don't we have access to numpy?)

Matrix multiplication: mmult([[11, 12], [21, 22]], [[-1, 0], [0, 1]]) returns [[-11, 12], [-21, 22]]

Definition at line 99 of file makeMuonMisalignmentScenario.py.

References ComparisonHelper.zip().

Referenced by cholesky(), and mvdot().

 def mmult(a, b):
     """Matrix multiplication: mmult([[11, 12], [21, 22]], [[-1, 0], [0, 1]]) returns [[-11, 12], [-21, 22]]"""
     return [[sum([i*j for i, j in zip(row, col)]) for col in zip(*b)] for row in a]
 

def makeMuonMisalignmentScenario.mtrans ( a )

Matrix transposition: mtrans([[11, 12], [21, 22]]) returns [[11, 21], [12, 22]]

Definition at line 107 of file makeMuonMisalignmentScenario.py.

 def mtrans(a):
     """Matrix transposition: mtrans([[11, 12], [21, 22]]) returns [[11, 21], [12, 22]]"""
     return [[a[j][i] for j in range(len(a[i]))] for i in range(len(a))]
 

def makeMuonMisalignmentScenario.mvdot	(	m,
		v
	)

Applies matrix m to vector v: mvdot([[-1, 0], [0, 1]], [12, 55]) returns [-12, 55]

Definition at line 103 of file makeMuonMisalignmentScenario.py.

References mmult().

Referenced by random6dof().

 def mvdot(m, v):
     """Applies matrix m to vector v: mvdot([[-1, 0], [0, 1]], [12, 55]) returns [-12, 55]"""
     return [i[0] for i in mmult(m, [[vi] for vi in v])]
 

def makeMuonMisalignmentScenario.random6dof ( )

Generate correlated random misalignments for all chambers.

Definition at line 151 of file makeMuonMisalignmentScenario.py.

References mvdot().

 def random6dof():
     randomunit = [gauss(0., 1.), gauss(0., 1.), gauss(0., 1.), gauss(0., 1.), gauss(0., 1.), gauss(0., 1.)]
     return mvdot(chomat, randomunit)
 

Variable Documentation

makeMuonMisalignmentScenario.action

Definition at line 39 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.args

Definition at line 41 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_phix = sum_phix/float(len(misal))

Definition at line 237 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_phixphix = sum_phixphix/float(len(misal))

Definition at line 256 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_phixphiy = sum_phixphiy/float(len(misal))

Definition at line 257 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_phixphiz = sum_phixphiz/float(len(misal))

Definition at line 258 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_phiy = sum_phiy/float(len(misal))

Definition at line 238 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_phiyphiy = sum_phiyphiy/float(len(misal))

Definition at line 259 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_phiyphiz = sum_phiyphiz/float(len(misal))

Definition at line 260 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_phiz = sum_phiz/float(len(misal))

Definition at line 239 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_phizphiz = sum_phizphiz/float(len(misal))

Definition at line 261 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_x = sum_x/float(len(misal))

Definition at line 234 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_xphix = sum_xphix/float(len(misal))

Definition at line 244 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_xphiy = sum_xphiy/float(len(misal))

Definition at line 245 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_xphiz = sum_xphiz/float(len(misal))

Definition at line 246 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_xx = sum_xx/float(len(misal))

Definition at line 241 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_xy = sum_xy/float(len(misal))

Definition at line 242 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_xz = sum_xz/float(len(misal))

Definition at line 243 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_y = sum_y/float(len(misal))

Definition at line 235 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_yphix = sum_yphix/float(len(misal))

Definition at line 249 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_yphiy = sum_yphiy/float(len(misal))

Definition at line 250 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_yphiz = sum_yphiz/float(len(misal))

Definition at line 251 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_yy = sum_yy/float(len(misal))

Definition at line 247 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_yz = sum_yz/float(len(misal))

Definition at line 248 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_z = sum_z/float(len(misal))

Definition at line 236 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_zphix = sum_zphix/float(len(misal))

Definition at line 253 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_zphiy = sum_zphiy/float(len(misal))

Definition at line 254 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_zphiz = sum_zphiz/float(len(misal))

Definition at line 255 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.ave_zz = sum_zz/float(len(misal))

Definition at line 252 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.cfgfile = file(outputName + "_convert_cfg.py", "w")

Convert it to an SQLite file with CMSSW.

Definition at line 351 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.chomat = cholesky(matrix)

Definition at line 147 of file makeMuonMisalignmentScenario.py.

string makeMuonMisalignmentScenario.components = "xx"

Definition at line 57 of file makeMuonMisalignmentScenario.py.

Referenced by MultiTrajectoryStateAssembler.addState(), TrackerOfflineValidation.bookDirHists(), TrackProducerAlgorithm< reco::GsfTrack >.buildTrack(), DAFTrackProducerAlgorithm.calculateNdof(), DAFTrackProducerAlgorithm.checkHits(), ForwardDetLayer.computeSurface(), MuonSeedOrcaPatternRecognition.countHits(), DAFTrackProducerAlgorithm.countingGoodHits(), edm::eventsetup.determinePreferred(), MRHChi2MeasurementEstimator.estimate(), GsfTrackProducerBase.fillStates(), DAFTrackProducerAlgorithm.filter(), CmsTrackerDetIdBuilder.iterate(), CmsMTDDetIdBuilder.iterate(), AlignableCompositeBuilder.maxNumComponents(), MultiTrajectoryStateMode.momentumFromModeCartesian(), MultiTrajectoryStateMode.momentumFromModePPhiEta(), GaussianStateConversions.multiGaussianStateFromTSOS(), GaussianStateConversions.multiGaussianStateFromVertex(), MultiGaussianStateTransform.multiState(), MultiGaussianStateTransform.multiState1D(), PFGsfHelper.PFGsfHelper(), MultiTrajectoryStateMode.positionFromModeCartesian(), edm::eventsetup.preferEverything(), TrackerOfflineValidation.prepareSummaryHists(), MuonAlignment.recursiveCopySurveyToAlignment(), MultiTrajectoryStateTransform.stateOnSurface(), TrajAnnealing.TrajAnnealing(), GaussianStateConversions.tsosFromMultiGaussianState(), GaussianStateConversions.vertexFromMultiGaussianState(), and MuonAlignmentOutputXML.writeComponents().

makeMuonMisalignmentScenario.default

Definition at line 12 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.dest

Definition at line 12 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.endcap

Definition at line 319 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.globalape = mmult(rot, mmult(localape, mtrans(rot)))

Definition at line 323 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.globalxx = globalape[0][0]

Definition at line 324 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.globalxy = globalape[0][1]

Definition at line 325 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.globalxz = globalape[0][2]

Definition at line 326 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.globalyy = globalape[1][1]

Definition at line 327 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.globalyz = globalape[1][2]

Definition at line 328 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.globalzz = globalape[2][2]

Definition at line 329 of file makeMuonMisalignmentScenario.py.

makeMuonMisalignmentScenario.help

Definition at line 12 of file makeMuonMisalignmentScenario.py.

list makeMuonMisalignmentScenario.localape = [[xx, xy, xz], [xy, yy, yz], [xz, yz, zz]]

Definition at line 322 of file makeMuonMisalignmentScenario.py.

list makeMuonMisalignmentScenario.matrix

Initial value:

 = [[   xx,    xy,    xz,    xphix,    xphiy,    xphiz],
           [   xy,    yy,    yz,    yphix,    yphiy,    yphiz],
           [   xz,    yz,    zz,    zphix,    zphiy,    zphiz],
           [xphix, yphix, zphix, phixphix, phixphiy, phixphiz],
           [xphiy, yphiy, zphiy, phixphiy, phiyphiy, phiyphiz],
           [xphiz, yphiz, zphiz, phixphiz, phiyphiz, phizphiz]]

Definition at line 140 of file makeMuonMisalignmentScenario.py.

Referenced by MuonErrorMatrix.adjust(), TestOutliers.analyze(), AlcaBeamMonitor.beginLuminosityBlock(), clone(), reco::GsfComponent5D.covariance(), reco::BeamSpot.covariance3D(), PhysicsTools::VarProcessor.deriv(), CSCSegFit.derivativeMatrix(), MuonSegFit.derivativeMatrix(), GEMCSCSegFit.derivativeMatrix(), AlcaBeamSpotFromDB.endLuminosityBlockProduce(), CSCDbStripConditions.fetchNoisifier(), BSFitter.Fit(), BSFitter.Fit_d0phi(), BSFitter.Fit_d_likelihood(), BSFitter.Fit_d_z_likelihood(), BSFitter.Fit_dres_z_likelihood(), BSFitter.Fit_z_chi2(), BSFitter.Fit_z_likelihood(), get_transform(), DTMuonMillepede.getbcsMatrix(), OnlineMetaDataRawToDigi.getBeamSpot(), DTMuonMillepede.getbqcMatrix(), DTMuonMillepede.getbsurveyMatrix(), DTMuonMillepede.getCcsMatrix(), DTMuonMillepede.getCqcMatrix(), DTMuonMillepede.getCsurveyMatrix(), DTMuonMillepede.getLagMatrix(), matrixSaver.getMatrix(), FWGeometry.getMatrix(), matrixSaver.getMatrixVector(), reco::GsfComponent5D.GsfComponent5D(), identity(), initializeMatrix(), FWGeometry.initMap(), FWGeometry.loadMap(), PhysicsTools.loadMatrix(), CSCConfigurableStripConditions.makeNoisifier(), reco::MET.mEtCorr(), TKinFitter.nbConstraints(), MultiTrackVertexLinkKinematicConstraint.parametersDerivative(), CombinedKinematicConstraint.parametersDerivative(), MultiTrackPointingKinematicConstraint.parametersDerivative(), MultiTrackVertexLinkKinematicConstraint.positionDerivative(), CombinedKinematicConstraint.positionDerivative(), MultiTrackPointingKinematicConstraint.positionDerivative(), TKinFitter.printMatrix(), BeamSpotProducer.produce(), BeamSpotOnlineProducer.produce(), KfComponentsHolder.projection(), KFBasedPixelFitter::MyBeamSpotHit.projectionMatrix(), DTMuonSLToSL.returnbSLMatrix(), DTMuonSLToSL.returnCSLMatrix(), reco::BeamSpot.rotatedCovariance3D(), PVFitter.runBXFitter(), PVFitter.runFitter(), BeamFitter.runPVandTrkFitter(), PhysicsTools.saveMatrix(), GsfMatrixTools.trace(), SignCaloSpecificAlgo.usePreviousSignif(), ValidateGeometry.validateCSChamberGeometry(), ValidateGeometry.validateCSCLayerGeometry(), ValidateGeometry.validateDTChamberGeometry(), ValidateGeometry.validateDTLayerGeometry(), ValidateGeometry.validateRPCGeometry(), ValidateGeometry.validateTrackerGeometry(), CSCSegFit.weightMatrix(), MuonSegFit.weightMatrix(), GEMCSCSegFit.weightMatrix(), EcalTBWeightsXMLTranslator.writeChi2WeightMatrix(), EcalTBWeightsXMLTranslator.writeWeightMatrix(), and FWTGeoRecoGeometryESProducer.~FWTGeoRecoGeometryESProducer().