60 def makeAlgs (self, config) :
61
62 log = logging.getLogger('MuonCalibrationConfig')
63
64
65
66 if self.calibMode == 'correctData_CB':
67 calibMode = 0
68 elif self.calibMode == 'correctData_IDMS':
69 calibMode = 1
70 elif self.calibMode == 'notCorrectData_IDMS':
71 calibMode = 2
72 elif self.calibMode == 'notCorrectData_CB':
73 calibMode = 3
74 elif self.calibMode == 'correctData_IDonly':
75 calibMode = 4
76 elif self.calibMode == 'correctData_MSonly':
77 calibMode = 5
78 else :
79 raise ValueError ("invalid calibMode: \"" + self.calibMode + "\". Allowed values are correctData_CB, correctData_IDMS, notCorrectData_IDMS, notCorrectData_CB, correctData_IDonly, correctData_MSonly")
80
81 inputContainer = "AnalysisMuons" if config.isPhyslite() else "Muons"
82 if self.inputContainer:
83 inputContainer = self.inputContainer
84 config.setSourceName (self.containerName, inputContainer)
85 config.setContainerMeta (self.containerName, 'calibMode', calibMode)
86
87
88 if config.wantCopy (self.containerName) :
89 alg = config.createAlgorithm( 'CP::AsgShallowCopyAlg', 'MuonShallowCopyAlg' )
90 alg.input = config.readName (self.containerName)
91 alg.output = config.copyName (self.containerName)
92 alg.outputType = 'xAOD::MuonContainer'
93 decorationList = ['DFCommonJetDr',
94 'DFCommonMuonPassIDCuts',
95 'DFCommonMuonPassPreselection',
96 'neflowisol20_CloseByCorr',
97 'ptvarcone30_Nonprompt_All_MaxWeightTTVA_pt1000_CloseByCorr',
98 'ptvarcone30_Nonprompt_All_MaxWeightTTVA_pt500_CloseByCorr',
99 'topoetcone20_CloseByCorr']
100 if self.addGlobalFELinksDep:
101 decorationList += ['neutralGlobalFELinks', 'chargedGlobalFELinks']
102 if config.dataType() is not DataType.Data:
103 decorationList += ['TruthLink']
104 if self.addGlobalFELinksDep or config.dataType() is not DataType.Data:
105 alg.declareDecorations = decorationList
106
107
108 alg = config.createAlgorithm( 'CP::AsgSelectionAlg',
109 'MuonEtaCutAlg' )
110 config.addPrivateTool( 'selectionTool', 'CP::AsgPtEtaSelectionTool' )
111 alg.selectionTool.maxEta = self.maxEta
112 alg.selectionDecoration = 'selectEta' + self.postfix + ',as_bits'
113 alg.particles = config.readName (self.containerName)
114 alg.preselection = config.getPreselection (self.containerName, '')
115 config.addSelection (self.containerName, '', alg.selectionDecoration)
116
117
118 alg = config.createAlgorithm( 'CP::MuonCalibrationAndSmearingAlg',
119 'MuonCalibrationAndSmearingAlg' )
120 config.addPrivateTool( 'calibrationAndSmearingTool',
121 'CP::MuonCalibTool' )
122
123 alg.calibrationAndSmearingTool.IsRun3Geo = config.geometry() >= LHCPeriod.Run3
124 alg.calibrationAndSmearingTool.calibMode = calibMode
125 if config.geometry() is LHCPeriod.Run4:
126 log.warning("Disabling NSW hits for Run4 geometry")
127 alg.calibrationAndSmearingTool.ExcludeNSWFromPrecisionLayers = True
128 else:
129 alg.calibrationAndSmearingTool.ExcludeNSWFromPrecisionLayers = self.excludeNSWFromPrecisionLayers and (config.geometry() >= LHCPeriod.Run3)
130
131
132
133 if self.useZeroPixMuons:
134 config.addPrivateTool( 'calibrationAndSmearingTool_ZeroPix',
135 'CP::MuonCalibTool' )
136 alg.calibrationAndSmearingTool_ZeroPix.IsRun3Geo = config.geometry() >= LHCPeriod.Run3
137 alg.calibrationAndSmearingTool_ZeroPix.calibMode = 5
138 alg.calibrationAndSmearingTool_ZeroPix.ExcludeNSWFromPrecisionLayers = alg.calibrationAndSmearingTool.ExcludeNSWFromPrecisionLayers
139 if self.zeroPixMuonType is not None:
140 alg.zeroPixMuonType = self.zeroPixMuonType
141
142 alg.muons = config.readName (self.containerName)
143 alg.muonsOut = config.copyName (self.containerName)
144 alg.preselection = config.getPreselection (self.containerName, '')
145 if config.isPhyslite() and not self.recalibratePhyslite :
146 alg.skipNominal = True
147
148
149 if self.minPt > 0:
150 alg = config.createAlgorithm( 'CP::AsgSelectionAlg', 'MuonPtCutAlg' )
151 alg.selectionDecoration = 'selectPt' + self.postfix + ',as_bits'
152 config.addPrivateTool( 'selectionTool', 'CP::AsgPtEtaSelectionTool' )
153 alg.particles = config.readName (self.containerName)
154 alg.selectionTool.minPt = self.minPt
155 alg.preselection = config.getPreselection (self.containerName, '')
156 config.addSelection (self.containerName, '', alg.selectionDecoration,
157 preselection = True)
158
159 alg = config.createAlgorithm( 'CP::AsgEnergyDecoratorAlg', 'EnergyDecorator' )
160 alg.particles = config.readName (self.containerName)
161
162 config.addOutputVar (self.containerName, 'pt', 'pt')
163 config.addOutputVar (self.containerName, 'eta', 'eta', noSys=True)
164 config.addOutputVar (self.containerName, 'phi', 'phi', noSys=True)
165 config.addOutputVar (self.containerName, 'e_%SYS%', 'e')
166 config.addOutputVar (self.containerName, 'charge', 'charge', noSys=True)
167
168
169 if self.decorateTruth and config.dataType() is not DataType.Data:
170 config.addOutputVar (self.containerName, "truthType", "truth_type", noSys=True)
171 config.addOutputVar (self.containerName, "truthOrigin", "truth_origin", noSys=True)
172
173 config.addOutputVar (self.containerName, 'muonType', 'muonType', noSys=True, enabled=self.writeColumnarToolVariables)
174
175