19def BPHY28Kernel(flags):
20
21
22 augList = []
23 skimList = []
24 thinList = []
25 thinTrkVtxList = []
26
27
28 from DerivationFrameworkBPhys.commonBPHYMethodsCfg import (
29 BPHY_V0ToolCfg, BPHY_InDetDetailedTrackSelectorToolCfg,
30 BPHY_VertexPointEstimatorCfg, BPHY_TrkVKalVrtFitterCfg,
31 AugOriginalCountsCfg)
32 from JpsiUpsilonTools.JpsiUpsilonToolsConfig import PrimaryVertexRefittingToolCfg
33 acc = ComponentAccumulator()
34 isSimulation = flags.Input.isMC
35 V0Tools = acc.popToolsAndMerge(BPHY_V0ToolCfg(flags, BPHYDerivationName))
36 vkalvrt = acc.popToolsAndMerge(BPHY_TrkVKalVrtFitterCfg(flags, BPHYDerivationName))
37 acc.addPublicTool(vkalvrt)
38 acc.addPublicTool(V0Tools)
39 trackselect = acc.popToolsAndMerge(BPHY_InDetDetailedTrackSelectorToolCfg(flags, BPHYDerivationName))
40 acc.addPublicTool(trackselect)
41 vpest = acc.popToolsAndMerge(BPHY_VertexPointEstimatorCfg(flags, BPHYDerivationName))
42 acc.addPublicTool(vpest)
43 PVrefit = acc.popToolsAndMerge(PrimaryVertexRefittingToolCfg(flags))
44 acc.addPublicTool(PVrefit)
45
46
47 doLRT = flags.Tracking.doLargeD0
48 if not doLRT :
print(
"BPHY28: LRT tracks disabled")
49 mainMuonInput = "StdWithLRTMuons" if doLRT else "Muons"
50 mainIDInput = "InDetWithLRTTrackParticles" if doLRT else "InDetTrackParticles"
51 if doLRT:
52 from DerivationFrameworkLLP.LLPToolsConfig import LRTMuonMergerAlg
53 from AthenaConfiguration.Enums import LHCPeriod
54 acc.merge(LRTMuonMergerAlg( flags,
55 PromptMuonLocation = "Muons",
56 LRTMuonLocation = "MuonsLRT",
57 OutputMuonLocation = mainMuonInput,
58 CreateViewCollection = True,
59 UseRun3WP = flags.GeoModel.Run == LHCPeriod.Run3))
60 from DerivationFrameworkInDet.InDetToolsConfig import InDetLRTMergeCfg
61 acc.merge(InDetLRTMergeCfg(flags))
62 toRelink = ["InDetTrackParticles", "InDetLargeD0TrackParticles"] if doLRT else []
63 MuonReLink = [ "Muons", "MuonsLRT" ] if doLRT else []
64
65 BPHY28_AugOriginalCounts = acc.popToolsAndMerge(
66 AugOriginalCountsCfg(flags, name = "BPHY28_AugOriginalCounts"))
67 augList += [ BPHY28_AugOriginalCounts ]
68
69 BPHY28MuMuFinder = CompFactory.Analysis.JpsiFinder(
70 name = "BPHY28MuMuFinder",
71 muAndMu = True,
72 muAndTrack = False,
73 TrackAndTrack = False,
74 assumeDiMuons = True,
75 invMassUpper = 100000.0,
76 invMassLower = 0.0,
77 Chi2Cut = 100.,
78 oppChargesOnly = True,
79 combOnly = True,
80 atLeastOneComb = False,
81 useCombinedMeasurement = False,
82 muonCollectionKey = mainMuonInput,
83 TrackParticleCollection = mainIDInput,
84 useV0Fitter = False,
85 TrkVertexFitterTool = vkalvrt,
86 V0VertexFitterTool = None,
87 TrackSelectorTool = trackselect,
88 VertexPointEstimator = vpest,
89 useMCPCuts = False )
90
91 BPHY28MuMuSelectAndWrite = CompFactory.DerivationFramework.Reco_Vertex(name = "BPHY28MuMuSelectAndWrite",
92 VertexSearchTool = BPHY28MuMuFinder,
93 OutputVtxContainerName = MuMuContainerName,
94 PVContainerName = "PrimaryVertices",
95 V0Tools = V0Tools,
96 PVRefitter = PVrefit,
97 RefPVContainerName = "SHOULDNOTBEUSED",
98 RelinkTracks = toRelink,
99 RelinkMuons = MuonReLink,
100 DoVertexType = 7)
101 augList += [ BPHY28MuMuSelectAndWrite ]
102 thinTrkVtxList += [ MuMuContainerName ]
103
104
105 BPHY28_Select_Jpsi2mumu = CompFactory.DerivationFramework.Select_onia2mumu(
106 name = "BPHY28_Select_Jpsi2mumu",
107 HypothesisName = "Jpsi",
108 InputVtxContainerName = MuMuContainerName,
109 V0Tools = V0Tools,
110 VtxMassHypo = 3096.916,
111 MassMin = 2000.0,
112 MassMax = 3600.0,
113 Chi2Max = 100, Do3d = False,
114 DoVertexType = 7)
115 augList += [ BPHY28_Select_Jpsi2mumu ]
116
117
118 BPHY28BsKKMuMu = CompFactory.Analysis.JpsiPlus2Tracks(name = "BPHY28BsKKMuMu",
119 kaonkaonHypothesis = True,
120 pionpionHypothesis = False,
121 kaonpionHypothesis = False,
122 trkThresholdPt = 500.0,
123 trkMaxEta = 3.0,
124 BMassUpper = 5900.0,
125 BMassLower = 4900.0,
126 DiTrackMassUpper = 1220,
127 DiTrackMassLower = 820,
128 Chi2Cut = 100.0,
129 TrkQuadrupletMassUpper = 6000.0,
130 TrkQuadrupletMassLower = 4800.0,
131 JpsiContainerKey = MuMuContainerName,
132 TrackParticleCollection = mainIDInput,
133 MuonsUsedInJpsi = mainMuonInput,
134 TrkVertexFitterTool = vkalvrt,
135 TrackSelectorTool = trackselect,
136 UseMassConstraint = False)
137
138
139 BPHY28BsKKSelectAndWrite = CompFactory.DerivationFramework.Reco_Vertex(name = "BPHY28BsKKSelectAndWrite",
140 VertexSearchTool = BPHY28BsKKMuMu,
141 OutputVtxContainerName = BsPhiMuMuContainerName,
142 PVContainerName = "PrimaryVertices",
143 V0Tools = V0Tools,
144 PVRefitter = PVrefit,
145 RefPVContainerName = "BPHY28RefittedPrimaryVertices",
146 RefitPV = True, Do3d = False,
147 RelinkTracks = toRelink,
148 MaxPVrefit = 10000, DoVertexType = 7)
149 augList += [ BPHY28BsKKSelectAndWrite ]
150 thinTrkVtxList += [ BsPhiMuMuContainerName ]
151
152
153 BPHY28_Select_Bs2KKMuMu = CompFactory.DerivationFramework.Select_onia2mumu(
154 name = "BPHY28_Select_Bs2KKMuMu",
155 HypothesisName = "Bs",
156 InputVtxContainerName = BsPhiMuMuContainerName,
157 V0Tools = V0Tools,
158 TrkMasses = [105.658, 105.658, 493.677, 493.677],
159 VtxMassHypo = 5366.3,
160 MassMin = 4900.0,
161 MassMax = 5900.0, Do3d = False,
162 Chi2Max = 100)
163 augList += [ BPHY28_Select_Bs2KKMuMu ]
164
165
166 if not isSimulation:
167 BPHY28_SelectBsKKMuMuEvent = CompFactory.DerivationFramework.xAODStringSkimmingTool(
168 name = "BPHY28_SelectBsKKMuMuEvent",
169 expression = "count(BPHY28BsKKMuMuCandidates.passed_Bs) > 0")
170
171 skimList += [ BPHY28_SelectBsKKMuMuEvent ]
172
173
174
176 BPHY28_Thin_VtxTracks = CompFactory.DerivationFramework.Thin_vtxTrk( name = "BPHY28_Thin_VtxTracks",
177 StreamName = streamName,
178 TrackParticleContainerName = "InDetTrackParticles",
179 VertexContainerNames = thinTrkVtxList,
180 IgnoreFlags = True )
181 thinList += [ BPHY28_Thin_VtxTracks ]
182
183
184 if doLRT:
185 BPHY28_Thin_VtxTracks_LRT = CompFactory.DerivationFramework.Thin_vtxTrk( name = "BPHY28_Thin_VtxTracks_LRT",
186 StreamName = streamName,
187 TrackParticleContainerName = "InDetLargeD0TrackParticles",
188 VertexContainerNames = thinTrkVtxList,
189 IgnoreFlags = True )
190 thinList += [ BPHY28_Thin_VtxTracks_LRT ]
191
192
193 for t in augList + skimList + thinList : acc.addPublicTool(t)
194 acc.addEventAlgo(CompFactory.DerivationFramework.DerivationKernel("BPHY28Kernel",
195 AugmentationTools = augList,
196
197 SkimmingTools = skimList,
198 ThinningTools = thinList))
199
200 return acc
201
202
void print(char *figname, TCanvas *c1)