Belle II Software development
analysisDQM.py
1# !/usr/bin/env python3
2
3
10
11"""
12This module defines functions to add analysis DQM modules.
13"""
14
15import basf2 as b2
16from stdPi0s import stdPi0s
17from stdCharged import stdPi, stdK
18import modularAnalysis as ma
19import stdV0s
20import vertex
21from variables import variables as vm
22
23
24def add_analysis_dqm(path):
25 """Add the analysis DQM modules to the ``path``.
26 Builds a list of muons, pi0's, Kshorts and the Fox Wolfram event shape variables.
27 Also M(mumu) for nominal beam energy monitoring.
28
29 Parameters:
30 path (basf2.Path): modules are loaded onto this path
31 """
32
33 # muons, Kshorts and pi0s
34 ma.fillParticleList('mu+:KLMDQM', 'p>1.5', path=path)
35 ma.fillParticleList('mu+:KLMDQM2', 'p>1.5 and abs(d0) < 2 and abs(z0) < 4', path=path)
36 ma.fillParticleList('gamma:physDQM', 'E > 0.15', path=path)
37 ma.fillParticleList('mu+:physDQM', 'pt>2. and abs(d0) < 2 and abs(z0) < 4', path=path)
38 ma.reconstructDecay('pi0:physDQM -> gamma:physDQM gamma:physDQM', '0.10 < M < 0.15', 1, writeOut=False, path=path)
39 # std Kshorts-TreeFit
40 stdV0s.stdKshorts(path=path, fitter='TreeFit', updateAllDaughters=True, writeOut=False, addSuffix=True)
41 ma.reconstructDecay('Upsilon:physDQM -> mu-:physDQM mu+:physDQM', '9 < M < 12', 1, writeOut=False, path=path)
42 # bhabha,hadrons
43 ma.fillParticleList('e+:physDQM', 'pt>0.2 and abs(d0) < 2 and abs(z0) < 4 and thetaInCDCAcceptance', path=path)
44 ma.reconstructDecay('Upsilon:ephysDQM -> e-:physDQM e+:physDQM', '4 < M < 12', 1, writeOut=False, path=path)
45 ma.fillParticleList('pi+:hadbphysDQM', 'p>0.1 and abs(d0) < 2 and abs(z0) < 4 and thetaInCDCAcceptance', path=path)
46
47 # have to manually create "all" lists of pi+ and photons to use inside buildEventShape
48 # to avoid loading the photons' beamBackgroundMVA variable on the DQM
49 ma.fillParticleList('pi+:evtshape', '', path=path)
50 ma.fillParticleList('gamma:evtshape', '', path=path)
51 ma.buildEventShape(
52 path=path,
53 inputListNames=['pi+:evtshape', 'gamma:evtshape'],
54 default_cleanup=False, # do not want any clean up
55 foxWolfram=True,
56 cleoCones=False,
57 collisionAxis=False,
58 harmonicMoments=False,
59 jets=False,
60 sphericity=False,
61 thrust=False)
62
63 dqm = b2.register_module('PhysicsObjectsDQM')
64 dqm.param('PI0PListName', 'pi0:physDQM')
65 dqm.param('KS0PListName', 'K_S0:merged_TreeFit')
66 dqm.param('UpsPListName', 'Upsilon:physDQM')
67 # bhabha,hadrons
68 dqm.param('UpsBhabhaPListName', 'Upsilon:ephysDQM')
69 dqm.param('UpsHadPListName', 'pi+:hadbphysDQM')
70
71 path.add_module(dqm)
72
73
74def add_mirabelle_dqm(path):
75 """Add the mirabelle DQM modules to the ``path``.
76 Runs on conditional paths depending on the software trigger results.
77 Building D*'s or dimuons on the conditional paths.
78
79 Parameters:
80 path (basf2.Path): modules are loaded onto this path
81 """
82 # Software Trigger to divert the path
83 MiraBelleMumu_path = b2.create_path()
84 MiraBelleZ0_path = b2.create_path()
85 MiraBelleDst1_path = b2.create_path()
86 MiraBelleNotDst1_path = b2.create_path()
87 MiraBelleDst2_path = b2.create_path()
88 MiraBelleTau_path = b2.create_path()
89 # bhabha,hadrons
90 MiraBelleBhabha_path = b2.create_path()
91 MiraBellehadronb2_path = b2.create_path()
92
93 trigger_skim_mumutight = path.add_module(
94 "TriggerSkim",
95 triggerLines=["software_trigger_cut&skim&accept_mumutight"],
96 resultOnMissing=0,
97 )
98 trigger_skim_mumutight.if_value("==1", MiraBelleMumu_path, b2.AfterConditionPath.CONTINUE)
99
100 trigger_skim_singlemuon = path.add_module(
101 "TriggerSkim",
102 triggerLines=["software_trigger_cut&filter&single_muon"],
103 resultOnMissing=0,
104 )
105 trigger_skim_singlemuon.if_value("==1", MiraBelleZ0_path, b2.AfterConditionPath.CONTINUE)
106
107 trigger_skim_dstar_1 = path.add_module(
108 "TriggerSkim",
109 triggerLines=["software_trigger_cut&skim&accept_dstar_1"],
110 resultOnMissing=0,
111 )
112 trigger_skim_dstar_1.if_value("==1", MiraBelleDst1_path, b2.AfterConditionPath.CONTINUE)
113
114 trigger_skim_not_dstar_1 = path.add_module(
115 "TriggerSkim",
116 triggerLines=["software_trigger_cut&skim&accept_dstar_1"],
117 expectedResult=0,
118 resultOnMissing=0,
119 )
120 trigger_skim_not_dstar_1.if_value("==1", MiraBelleNotDst1_path, b2.AfterConditionPath.CONTINUE)
121 trigger_skim_dstar_2 = MiraBelleNotDst1_path.add_module(
122 "TriggerSkim",
123 triggerLines=["software_trigger_cut&skim&accept_dstar_2"],
124 resultOnMissing=0,
125 )
126 trigger_skim_dstar_2.if_value("==1", MiraBelleDst2_path, b2.AfterConditionPath.CONTINUE)
127 trigger_skim_tau = path.add_module(
128 "TriggerSkim",
129 triggerLines=["software_trigger_cut&skim&accept_tau_2trk", "software_trigger_cut&skim&accept_tau_Ntrk"],
130 expectedResult=1,
131 logicMode="or",
132 resultOnMissing=0,
133 )
134 trigger_skim_tau.if_value("==1", MiraBelleTau_path, b2.AfterConditionPath.CONTINUE)
135 # bhabha,hadrons
136 trigger_skim_bhabhaall = path.add_module(
137 "TriggerSkim",
138 triggerLines=["software_trigger_cut&skim&accept_bhabha_trk_ecl"],
139 resultOnMissing=0,
140 )
141 trigger_skim_bhabhaall.if_value("==1", MiraBelleBhabha_path, b2.AfterConditionPath.CONTINUE)
142 trigger_skim_hadronb2 = path.add_module(
143 "TriggerSkim",
144 triggerLines=["software_trigger_cut&skim&accept_hadronb2"],
145 resultOnMissing=0,
146 )
147 trigger_skim_hadronb2.if_value("==1", MiraBellehadronb2_path, b2.AfterConditionPath.CONTINUE)
148
149 # MiraBelle di-muon path
150 ma.fillParticleList('mu+:physMiraBelle', '', path=MiraBelleMumu_path)
151 ma.reconstructDecay('Upsilon:physMiraBelle -> mu+:physMiraBelle mu-:physMiraBelle', '9 < M < 12', path=MiraBelleMumu_path)
152 MiraBelleMumu = b2.register_module('PhysicsObjectsMiraBelle')
153 MiraBelleMumu.param('MuPListName', 'mu+:physMiraBelle')
154 MiraBelleMumu.param('MuMuPListName', 'Upsilon:physMiraBelle')
155 MiraBelleMumu_path.add_module(MiraBelleMumu)
156
157 # MiraBelle Z0 path
158 ma.fillParticleList('mu+:physMiraBelleZ0', 'abs(dr) < 2 and abs(dz) < 5', path=MiraBelleZ0_path)
159 ma.reconstructDecay('Z0:physMiraBelle -> mu+:physMiraBelleZ0 mu-:physMiraBelleZ0',
160 'nCleanedTracks(abs(dr) < 2 and abs(dz) < 5) and 9 < M and M < 12', path=MiraBelleZ0_path)
161 MiraBelleModule = b2.register_module('PhysicsObjectsMiraBelle')
162 MiraBelleModule.param('Z0PListName', 'Z0:physMiraBelle')
163 MiraBelleZ0_path.add_module(MiraBelleModule)
164
165 # MiraBelle D* (followed by D0 -> K pi) path
166 ma.fillParticleList('pi+:MiraBelleDst1', 'abs(d0)<0.5 and abs(z0)<3', path=MiraBelleDst1_path)
167 ma.fillParticleList('K+:MiraBelleDst1', 'abs(d0)<0.5 and abs(z0)<3', path=MiraBelleDst1_path)
168 ma.reconstructDecay('D0:MiraBelleDst1_kpi -> K-:MiraBelleDst1 pi+:MiraBelleDst1', '1.7 < M < 2.1', path=MiraBelleDst1_path)
169 ma.reconstructDecay('D*+:MiraBelleDst1_kpi -> D0:MiraBelleDst1_kpi pi+:MiraBelleDst1',
170 'useCMSFrame(p) > 2.5 and massDifference(0) < 0.16', path=MiraBelleDst1_path)
171 MiraBelleDst1 = b2.register_module('PhysicsObjectsMiraBelleDst')
172 MiraBelleDst1.param('DstListName', 'D*+:MiraBelleDst1_kpi')
173 MiraBelleDst1_path.add_module(MiraBelleDst1)
174
175 # MiraBelle D* (followed by D0 -> K pi pi0) path
176 ma.fillParticleList('pi+:MiraBelleDst2', 'abs(d0)<0.5 and abs(z0)<3', path=MiraBelleDst2_path)
177 ma.fillParticleList('K+:MiraBelleDst2', 'abs(d0)<0.5 and abs(z0)<3', path=MiraBelleDst2_path)
178 stdPi0s(listtype='eff60_May2020', writeOut=False, path=MiraBelleDst2_path)
179 ma.reconstructDecay(
180 'D0:MiraBelleDst2_kpipi0 -> K-:MiraBelleDst2 pi+:MiraBelleDst2 pi0:eff60_May2020',
181 '1.7 < M < 2.1',
182 path=MiraBelleDst2_path)
183 ma.reconstructDecay('D*+:MiraBelleDst2_kpipi0 -> D0:MiraBelleDst2_kpipi0 pi+:MiraBelleDst2',
184 'useCMSFrame(p) > 2.5 and massDifference(0) < 0.16', path=MiraBelleDst2_path)
185 MiraBelleDst2 = b2.register_module('PhysicsObjectsMiraBelleDst2')
186 MiraBelleDst2.param('DstListName', 'D*+:MiraBelleDst2_kpipi0')
187 MiraBelleDst2_path.add_module(MiraBelleDst2)
188
189 # MiraBelle tau path
190 ma.fillParticleList('pi+:MiraBelleTau', 'abs(d0)<1 and abs(z0)<3', path=MiraBelleTau_path)
191 ma.fillParticleList('gamma:MiraBelleTau', '0.1<E and 1.5<clusterNHits and thetaInCDCAcceptance', path=MiraBelleTau_path)
192 ParticleListForEvent = ['pi+:MiraBelleTau', 'gamma:MiraBelleTau']
193 ma.buildEventShape(ParticleListForEvent,
194 allMoments=False, foxWolfram=False, cleoCones=False,
195 sphericity=False, jets=False, path=MiraBelleTau_path)
196 ma.buildEventKinematics(ParticleListForEvent, path=MiraBelleTau_path)
197
198 vm.addAlias("nGoodTracks", "nParticlesInList(pi+:MiraBelleTau)")
199 vm.addAlias("maxPt", "maxPtInList(pi+:MiraBelleTau)")
200 vm.addAlias("E_ECLtrk", "formula(totalECLEnergyOfParticlesInList(pi+:MiraBelleTau))")
201 vm.addAlias("maxOp", "useCMSFrame(maxOpeningAngleInList(pi+:MiraBelleTau))")
202
203 ma.cutAndCopyLists('pi+:MiraBelleTau_para', 'pi+:MiraBelleTau', 'cosToThrustOfEvent>0', path=MiraBelleTau_path)
204 ma.cutAndCopyLists('pi+:MiraBelleTau_anti', 'pi+:MiraBelleTau', 'cosToThrustOfEvent<0', path=MiraBelleTau_path)
205 vm.addAlias('pionInParaThrust', 'countInList(pi+:MiraBelleTau_para)')
206 vm.addAlias('pionInAntiThrust', 'countInList(pi+:MiraBelleTau_anti)')
207
208 Tau1x1Cuts = '[nGoodTracks == 2] and [3 < visibleEnergyOfEventCMS < 10]'
209 Tau1x1Cuts += 'and [missingMomentumOfEvent_theta < 2.618]'
210 Tau1x1Cuts += 'and [E_ECLtrk < 6] and [maxPt > 1] and [maxOp < 3.106686]'
211 Tau1x1Cuts += 'and [pionInParaThrust == 1 and pionInAntiThrust == 1]'
212 ma.reconstructDecay('Z0:physMiraBelleTau1x1 -> pi+:MiraBelleTau pi-:MiraBelleTau',
213 f'{Tau1x1Cuts}', path=MiraBelleTau_path)
214
215 Tau1x3Cuts = '[nGoodTracks == 4] and [3 < visibleEnergyOfEventCMS < 10.5]'
216 Tau1x3Cuts += 'and [E_ECLtrk < 6] and [maxPt > 1] and [maxOp < 3.106686]'
217 Tau1x3Cuts += 'and [[pionInParaThrust == 1 and pionInAntiThrust == 3] or [pionInParaThrust == 3 and pionInAntiThrust == 1]]'
218 ma.reconstructDecay('Z0:physMiraBelleTau1x3 -> pi+:MiraBelleTau pi-:MiraBelleTau pi+:MiraBelleTau pi-:MiraBelleTau',
219 f'{Tau1x3Cuts}', path=MiraBelleTau_path)
220
221 MiraBelleTau = b2.register_module('PhysicsObjectsMiraBelleTau')
222 MiraBelleTau.param('tautau1x1PListName', 'Z0:physMiraBelleTau1x1')
223 MiraBelleTau.param('tautau1x3PListName', 'Z0:physMiraBelleTau1x3')
224 MiraBelleTau_path.add_module(MiraBelleTau)
225
226 # bhabha,hadrons
227 ma.fillParticleList(
228 'e+:physMiraBelle',
229 'pt>0.2 and abs(d0) < 2 and abs(z0) < 4 and thetaInCDCAcceptance',
230 path=MiraBelleBhabha_path)
231 ma.reconstructDecay('Upsilon:ephysMiraBelle -> e+:physMiraBelle e-:physMiraBelle', '4 < M < 12', path=MiraBelleBhabha_path)
232 MiraBelleBhabha = b2.register_module('PhysicsObjectsMiraBelleBhabha')
233 MiraBelleBhabha.param('ePListName', 'e+:physMiraBelle')
234 MiraBelleBhabha.param('bhabhaPListName', 'Upsilon:ephysMiraBelle')
235 MiraBelleBhabha_path.add_module(MiraBelleBhabha)
236 ma.fillParticleList(
237 'pi+:hadb2physMiraBelle',
238 'p>0.1 and abs(d0) < 2 and abs(z0) < 4 and thetaInCDCAcceptance',
239 path=MiraBellehadronb2_path)
240 MiraBellehadronb = b2.register_module('PhysicsObjectsMiraBelleHadron')
241 MiraBellehadronb.param('hadronb2piPListName', 'pi+:hadb2physMiraBelle')
242 MiraBellehadronb2_path.add_module(MiraBellehadronb)
243
244
245# Selection for the EcmsBB analysis
246def get_hadB_path(path):
247 """ Selects the hadronic B decays, function returns corresponding path """
248
249 # module to be run prior the collector
250 path_hadB = b2.create_path()
251 trigger_skim_BB = path.add_module(
252 "TriggerSkim",
253 triggerLines=["software_trigger_cut&skim&accept_btocharm"],
254 resultOnMissing=0,
255 )
256 trigger_skim_BB.if_value("==1", path_hadB, b2.AfterConditionPath.CONTINUE)
257
258 stdPi(listtype='loose', path=path_hadB)
259 stdK(listtype='good', path=path_hadB)
260 stdPi0s(listtype='eff40_May2020', path=path_hadB)
261
262 ma.cutAndCopyList("pi+:hadB", "pi+:loose", "[abs(dz)<2.0] and [abs(dr)<0.5]", path=path_hadB)
263 ma.cutAndCopyList("K+:hadB", "K+:good", "[abs(dz)<2.0] and [abs(dr)<0.5]", path=path_hadB)
264
265 ma.cutAndCopyList("pi0:hadB", "pi0:eff40_May2020", "", path=path_hadB)
266
267
270
271 DcutLoose = '1.7 < M < 2.1'
272 Dcut = '1.830 < M < 1.894'
273 # Reconstructs D0s and sets decay mode identifiers
274 ma.reconstructDecay(decayString='D0:hadB_Kpi -> K-:hadB pi+:hadB', cut=DcutLoose, dmID=1, path=path_hadB)
275 ma.reconstructDecay(decayString='D0:hadB_Kpipi0 -> K-:hadB pi+:hadB pi0:hadB',
276 cut=DcutLoose, dmID=2, path=path_hadB)
277 ma.reconstructDecay(decayString='D0:hadB_Kpipipi -> K-:hadB pi+:hadB pi-:hadB pi+:hadB',
278 cut=DcutLoose, dmID=3, path=path_hadB)
279
280 # Performs mass constrained fit for all D0 candidates
281 vertex.kFit(list_name='D0:hadB_Kpi', conf_level=0.0, fit_type='mass', path=path_hadB)
282 # vertex.kFit(list_name='D0:hadB_Kpipi0', conf_level=0.0, fit_type='mass', path=path_hadB)
283 vertex.kFit(list_name='D0:hadB_Kpipipi', conf_level=0.0, fit_type='mass', path=path_hadB)
284
285 ma.applyCuts("D0:hadB_Kpi", Dcut, path=path_hadB)
286 ma.applyCuts("D0:hadB_Kpipi0", Dcut, path=path_hadB)
287 ma.applyCuts("D0:hadB_Kpipipi", Dcut, path=path_hadB)
288
289 DStarcutLoose = 'massDifference(0) < 0.16'
290
291 # Reconstructs D*-s and sets decay mode identifiers
292 ma.reconstructDecay(decayString='D*+:hadB_D0pi_Kpi -> D0:hadB_Kpi pi+:hadB', cut=DStarcutLoose, dmID=1, path=path_hadB)
293 ma.reconstructDecay(decayString='D*+:hadB_D0pi_Kpipi0 -> D0:hadB_Kpipi0 pi+:hadB',
294 cut=DStarcutLoose, dmID=2, path=path_hadB)
295 ma.reconstructDecay(decayString='D*+:hadB_D0pi_Kpipipi -> D0:hadB_Kpipipi pi+:hadB',
296 cut=DStarcutLoose, dmID=3, path=path_hadB)
297
298 BcutLoose = '[ useCMSFrame(p) < 1.6 ] and [abs(dM) < 0.25]'
299 Bcut = '[ useCMSFrame(p) < 1.2 ] and [abs(dM) < 0.05]'
300
301 # Reconstructs the signal B0 candidates from Dstar
302 ma.reconstructDecay(decayString='B0:hadB_Dstpi_D0pi_Kpi -> D*-:hadB_D0pi_Kpi pi+:hadB',
303 cut=BcutLoose,
304 dmID=1, path=path_hadB)
305 ma.reconstructDecay(decayString='B0:hadB_Dstpi_D0pi_Kpipi0 -> D*-:hadB_D0pi_Kpipi0 pi+:hadB',
306 cut=BcutLoose,
307 dmID=2, path=path_hadB)
308 ma.reconstructDecay(decayString='B0:hadB_Dstpi_D0pi_Kpipipi -> D*-:hadB_D0pi_Kpipipi pi+:hadB',
309 cut=BcutLoose,
310 dmID=3, path=path_hadB)
311
312 vertex.treeFit('B0:hadB_Dstpi_D0pi_Kpi', updateAllDaughters=True, ipConstraint=True, path=path_hadB)
313 vertex.treeFit('B0:hadB_Dstpi_D0pi_Kpipi0', updateAllDaughters=True, ipConstraint=True, path=path_hadB)
314 vertex.treeFit('B0:hadB_Dstpi_D0pi_Kpipipi', updateAllDaughters=True, ipConstraint=True, path=path_hadB)
315
316
319
320 # Reconstructs charged D mesons and sets decay mode identifiers
321 ma.reconstructDecay(decayString='D-:hadB_Kpipi -> K+:hadB pi-:hadB pi-:hadB',
322 cut=DcutLoose, dmID=4, path=path_hadB)
323
324 vertex.kFit(list_name='D-:hadB_Kpipi', conf_level=0.0, fit_type='mass', path=path_hadB)
325 ma.applyCuts("D-:hadB_Kpipi", '1.844 < M < 1.894', path=path_hadB)
326
327 # Reconstructs the signal B candidates
328 ma.reconstructDecay(decayString='B0:hadB_Dpi_Kpipi -> D-:hadB_Kpipi pi+:hadB',
329 cut=BcutLoose, dmID=4, path=path_hadB)
330
331
334
335 # Reconstructs the signal B- candidates
336 ma.reconstructDecay(decayString='B-:hadB_D0pi_Kpi -> D0:hadB_Kpi pi-:hadB',
337 cut=BcutLoose,
338 dmID=5, path=path_hadB)
339 ma.reconstructDecay(decayString='B-:hadB_D0pi_Kpipi0 -> D0:hadB_Kpipi0 pi-:hadB',
340 cut=BcutLoose,
341 dmID=6, path=path_hadB)
342 ma.reconstructDecay(decayString='B-:hadB_D0pi_Kpipipi -> D0:hadB_Kpipipi pi-:hadB',
343 cut=BcutLoose,
344 dmID=7, path=path_hadB)
345
346 vertex.treeFit('B-:hadB_D0pi_Kpi', updateAllDaughters=True, ipConstraint=True, path=path_hadB)
347 vertex.treeFit('B-:hadB_D0pi_Kpipi0', updateAllDaughters=True, ipConstraint=True, path=path_hadB)
348 vertex.treeFit('B-:hadB_D0pi_Kpipipi', updateAllDaughters=True, ipConstraint=True, path=path_hadB)
349
350 ma.copyLists(
351 outputListName='B0:hadB_combined',
352 inputListNames=[
353 'B0:hadB_Dstpi_D0pi_Kpi',
354 'B0:hadB_Dstpi_D0pi_Kpipi0',
355 'B0:hadB_Dstpi_D0pi_Kpipipi',
356 'B0:hadB_Dpi_Kpipi'
357 ],
358 path=path_hadB)
359
360 ma.copyLists(
361 outputListName='B-:hadB_combined',
362 inputListNames=[
363 'B-:hadB_D0pi_Kpi',
364 'B-:hadB_D0pi_Kpipi0',
365 'B-:hadB_D0pi_Kpipipi',
366 ],
367 path=path_hadB)
368
369 # Builds the rest of event object, which contains all particles not used in the reconstruction of B0 candidates.
370 ma.buildRestOfEvent(target_list_name='B0:hadB_combined', path=path_hadB)
371
372 # CleanMask is identical for B0 and B+
373 cleanMask = ('cleanMask', 'nCDCHits > 0 and useCMSFrame(p)<=3.2', 'p >= 0.05 and useCMSFrame(p)<=3.2')
374
375 # Calculates the continuum suppression variables
376 ma.appendROEMasks(list_name='B0:hadB_combined', mask_tuples=[cleanMask], path=path_hadB)
377 ma.buildContinuumSuppression(list_name='B0:hadB_combined', roe_mask='cleanMask', path=path_hadB)
378
379 # Builds the rest of event object, which contains all particles not used in the reconstruction of B- candidates.
380 ma.buildRestOfEvent(target_list_name='B-:hadB_combined', path=path_hadB)
381
382 # Calculates the continuum suppression variables
383 ma.appendROEMasks(list_name='B-:hadB_combined', mask_tuples=[cleanMask], path=path_hadB)
384 ma.buildContinuumSuppression(list_name='B-:hadB_combined', roe_mask='cleanMask', path=path_hadB)
385
386 ma.applyCuts("B0:hadB_combined", "[R2 < 0.3] and " + Bcut, path=path_hadB)
387 ma.applyCuts("B-:hadB_combined", "[R2 < 0.3] and " + Bcut, path=path_hadB)
388
389 MiraBelleEcmsBB = b2.register_module('PhysicsObjectsMiraBelleEcmsBB')
390 MiraBelleEcmsBB.param('B0ListName', 'B0:hadB_combined')
391 MiraBelleEcmsBB.param('BmListName', 'B-:hadB_combined')
392 path_hadB.add_module(MiraBelleEcmsBB)
393
394 return path_hadB
stdV0s.stdKshorts
stdKshorts(prioritiseV0=True, fitter="TreeFit", path=None, updateAllDaughters=False, writeOut=False, addSuffix=False)
Definition stdV0s.py:25
vertex.treeFit
treeFit(list_name, conf_level=0.001, massConstraint=[], ipConstraint=False, updateAllDaughters=False, massConstraintDecayString='', massConstraintMassValues=[], customOriginConstraint=False, customOriginVertex=[0.001, 0, 0.0116], customOriginCovariance=[0.0048, 0, 0, 0, 0.003567, 0, 0, 0, 0.0400], originDimension=3, treatAsInvisible='', ignoreFromVertexFit='', path=None)
Definition vertex.py:237
vertex.kFit
kFit(list_name, conf_level, fit_type='vertex', constraint='', daughtersUpdate=False, decay_string='', massConstraint=[], recoilMass=0, smearing=0, path=None)
Definition vertex.py:126