def configureEarlyDataAlignment(fixQOverPBias=True):
TAlignment().WriteCondSubDetList += ['Velo', 'TT', 'IT', 'OT', 'Muon']
# define the alignment elements
elements = Alignables()
elements.Velo("None")
elements.VeloRight("Tx")
elements.VeloLeft("Tx")
elements.IT("None")
elements.ITBoxes("TxTzRz")
elements.ITLayers("TxTz")
elements.OT("None")
elements.OTCFrames("Tx")
elements.OTCFrameLayers("Tz")
elements.TT("None")
elements.TTLayers("None")
elements.TTModules("TxRz")
elements.Tracker("None")
TAlignment().ElementsToAlign = list(elements)
# make sure that the velo stays where it was
TAlignment().Constraints = constraints = []
constraints.append("VeloHalfAverage : Velo/Velo(Left|Right) : Tx ")
# fix the q/p scale by not moving T in X. note that you do not
# want to do this if you use D0 in the alignment
if fixQOverPBias:
constraints.append("OT3X : OT/T3X1U.Side : Tx")
# tweak the survey a little bit to fix the z-scale to survey
surveyconstraints = SurveyConstraints()
surveyconstraints.All()
# make sure we fix the z-scale
surveyconstraints.XmlUncertainties += [
"OT/T3X1U : 0.5 0.5 0.00001 0.0001 0.0001 0.0001"
]
def configureVeloHalfAlignment():
TAlignment().WriteCondSubDetList += ['Velo']
# define the alignment elements
elements = Alignables()
elements.Velo("None")
elements.VeloRight("TxTyTzRxRyRz")
elements.VeloLeft("TxTyTzRxRyRz")
TAlignment().ElementsToAlign = list(elements)
# make sure that the velo stays where it was
TAlignment().Constraints = constraints = []
constraints.append(
"VeloHalfAverage : Velo/Velo(Left|Right) : Tx Ty Tz Rx Ry Rz")
# tweak the survey a little bit to fix the z-scale to survey
surveyconstraints = SurveyConstraints()
surveyconstraints.All()
def configureVeloModuleAlignment():
'''
Align Modues only for only the main degrees of freedom Tx Ty Rz
Constrain two modules in each half
'''
TAlignment().WriteCondSubDetList += ['Velo']
elements = Alignables()
elements.Velo("None")
elements.VeloRight("None")
elements.VeloLeft("None")
elements.VeloModules("TxTyRz")
elements.VeloPhiSensors("None")
elements.VeloRSensors("None")
TAlignment().ElementsToAlign += list(elements)
# Constraints
surveyconstraints = SurveyConstraints()
surveyconstraints.All()
constraints = []
# Fix 2 modules
constraints.append("VeloFixModule10 : Velo/VeloLeft/Module10: Tx Ty Rz")
constraints.append("VeloFixModule11 : Velo/VeloRight/Module11: Tx Ty Rz")
constraints.append("VeloFixModule32 : Velo/VeloLeft/Module32: Tx Ty Rz")
constraints.append("VeloFixModule33 : Velo/VeloRight/Module33: Tx Ty Rz")
constraints.append(
"VeloFixSensors10 : Velo/VeloLeft/Module10/RPhiPair10/.*: Tx Ty Rz")
constraints.append(
"VeloFixSensors11 : Velo/VeloRight/Module11/RPhiPair11/.*: Tx Ty Rz")
constraints.append(
"VeloFixSensors32 : Velo/VeloLeft/Module32/RPhiPair32/.*: Tx Ty Rz")
constraints.append(
"VeloFixSensors33 : Velo/VeloRight/Module33/RPhiPair33/.*: Tx Ty Rz")
# Constrain average
# constraints.append("VeloAInternal : Velo/VeloRight/Module..: Tx Ty Rz Szx Szy")
# constraints.append("VeloCInternal : Velo/VeloLeft/Module..: Tx Ty Rz Szx Szy")
TAlignment().Constraints = constraints
def configureVeloSensorAlignment():
TAlignment().WriteCondSubDetList += ['Velo']
# To align at the sensor level, we align 'modules' and
#'phi-sensors'. The R-sensor is always equal to the modules. (We
#keep its 'deltas' equal to zero.) However, because of a feature
#in the way we keep track of mother-daughter relations, it is
#important to create alignables for the R-sensors as well: use
#"None" to specify an alignable without active degrees of
#freedom.
# Note also that it doesn't make sense to align R modules for
# rotations around the z-axis. Rather than aligning the phi sensor
# for this rotation, we put the rotation at the module level.
elements = Alignables()
elements.Velo("None")
elements.VeloRight("TxTyTzRxRyRz")
elements.VeloLeft("TxTyTzRxRyRz")
elements.VeloModules("TxTyTzRxRyRz")
elements.VeloPhiSensors("TxTy")
elements.VeloRSensors("None")
TAlignment().ElementsToAlign += list(elements)
# make sure that the velo stays where it was. Important note: the
# dofs here must match the ones that we actually align for. If you
# specify too many, things will go rather wrong.
TAlignment().Constraints.append(
"VeloHalfAverage : Velo/Velo(Left|Right) : Tx Ty Tz Rx Ry Rz")
# now constrain the total shearing and z-scale in each Velo half
# to get the z-axis nominal. we'll do this by tightening the
# survey errors.
# we need to align the relative twist of the two halves. the
# C-side (minus X, Right) is known to have a corkscrew ('twist')
# in the base plate. so, for the right half we don't fix the Rz rotation:
surveyconstraints = SurveyConstraints()
surveyconstraints.All()
def configureMuonAlignment(surveyVer='latest'):
TAlignment().WriteCondSubDetList += ['Muon']
# change convergence criteria
updatetool = Al__AlignUpdateTool("Al::AlignUpdateTool")
updatetool.MaxDeltaChi2PDofForConvergence = 2
updatetool.MaxModeDeltaChi2ForConvergence = 10
# define the alignment elements
elements = Alignables()
elements.MuonHalfStations("TxTy")
TAlignment().ElementsToAlign = list(elements)
# make sure that the velo stays where it was
constraints = []
TAlignment().Constraints = constraints
if surveyVer == None: return
# tweak the survey a little bit to fix the z-scale to survey
surveyconstraints = SurveyConstraints()
surveyconstraints.All(surveyVer)
"""
def configureTrackerAlignment(fixQOverPBias=False, surveyVer='latest'):
TAlignment().WriteCondSubDetList += ['TT', 'IT', 'OT']
# define the alignment elements
elements = Alignables()
elements.Velo("None")
elements.VeloRight("None")
elements.VeloLeft("None")
elements.IT("None")
elements.ITBoxes("TxTzRz")
elements.ITLayers("None")
elements.OT("None")
elements.OTCFrames("TxRz")
elements.OTCFrameLayers("Tz")
elements.OTModules("TxRz")
elements.TT("None")
elements.TTLayers("Tz")
elements.TTModules("TxRz")
elements.Tracker("None")
TAlignment().ElementsToAlign = list(elements)
# make sure that the velo stays where it was
constraints = []
# fix the q/p scale by not moving T in X. note that you do not
# want to do this if you use D0 in the alignment
if fixQOverPBias:
constraints.append("OT3X : OT/T3X1U.Side : Tx")
TAlignment().Constraints = constraints
if surveyVer == None: return
# tweak the survey a little bit to fix the z-scale to survey
surveyconstraints = SurveyConstraints()
surveyconstraints.All(surveyVer)
# make sure we fix the z-scale
surveyconstraints.XmlUncertainties += [
"OT/T3X1U : 0.5 0.5 0.00001 0.0001 0.0001 0.0001"
]
#elements.TTBoxes("Ry")
#elements.TTHalfLayers("Tx")
#elements.TTHalfModules("TxRz")
elements.TTModules("TxRz")
elements.Tracker("None")
# let's try this:
# elements.Tracker("TxTzRxRz") This doesn't solve anything ... we need to look somewhere else
# to get Fred's y alignment from survey
elements.OTStations("Ty")
elements.TT("Ty")
elements.ITBoxes("Ty")
elements.ITLayers("Ty")
surveyconstraints = SurveyConstraints()
surveyconstraints.All()
constraints = []
# make sure that the velo stays where it was
constraints.append(
"VeloHalfAverage : Velo/Velo(Left|Right) : Tx Ty Tz Rx Ry : total")
# we fix the OT x, until we have seen magnet-off data
# fix the X-position of OT in T3 two C-frames:
# release the constraint now that we have J/psi
# constraints.append("OT3X : OT/T3X1U.Side : Tx")
# now fix OT2 because this was the only one that wasn't moved
constraints.append("OT2X : OT/T2X1U.Side : Tx")
# fix two C-frame layers in Z
def configureVeloAlignment(fixMeanInHalf=True):
'''
This should be the default alignment for the Automatic alignment procedure
Align 2-halves for all degree of freedom and
Modues only for only the main degrees of freedom Tx Ty Rz
Constrain the global Velo position and two modules in each half
'''
TAlignment().WriteCondSubDetList += ['Velo']
elements = Alignables()
elements.Velo("None")
elements.VeloRight("TxTyTzRxRyRz")
elements.VeloLeft("TxTyTzRxRyRz")
elements.VeloModules("TxTyRz")
elements.VeloPhiSensors("None")
elements.VeloRSensors("None")
TAlignment().ElementsToAlign += list(elements)
# Constraints
surveyconstraints = SurveyConstraints()
surveyconstraints.All()
constraints = []
surveyconstraints.Constraints += [
"Velo : 0 0 0 -0.0001 0 -0.0001 : 0.2 0.2 0.2 0.0001 0.0001 0.001",
"Velo/Velo(Right|Left) : 0 0 0 0 0 0 : 10 1 0.4 0.01 0.01 0.001"
]
# make sure that the velo stays where it was. Important note: the
# dofs here must match the ones that we actually align for. If you
# specify too many, things will go rather wrong.
constraints.append(
"VeloHalfAverage : Velo/Velo(Left|Right) : Tx Ty Tz Rx Ry Rz")
if fixMeanInHalf:
# Constrain average
constraints.append(
"VeloAInternal : Velo/VeloRight/Module..: Tx Ty Rz Szx Szy")
constraints.append(
"VeloCInternal : Velo/VeloLeft/Module..: Tx Ty Rz Szx Szy")
else:
# Fix 2 modules
constraints.append(
"VeloFixModule10 : Velo/VeloLeft/Module10: Tx Ty Rz")
constraints.append(
"VeloFixModule11 : Velo/VeloRight/Module11: Tx Ty Rz")
constraints.append(
"VeloFixModule32 : Velo/VeloLeft/Module32: Tx Ty Rz")
constraints.append(
"VeloFixModule33 : Velo/VeloRight/Module33: Tx Ty Rz")
constraints.append(
"VeloFixSensors10 : Velo/VeloLeft/Module10/RPhiPair10/.*: Tx Ty Rz"
)
constraints.append(
"VeloFixSensors11 : Velo/VeloRight/Module11/RPhiPair11/.*: Tx Ty Rz"
)
constraints.append(
"VeloFixSensors32 : Velo/VeloLeft/Module32/RPhiPair32/.*: Tx Ty Rz"
)
constraints.append(
"VeloFixSensors33 : Velo/VeloRight/Module33/RPhiPair33/.*: Tx Ty Rz"
)
TAlignment().Constraints = constraints