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