4import modularAnalysis
as ma
5from variables
import variables
as vm
10 filename=b2.find_file(
"starterkit/2021/fei_skimmed_xulnu.udst.root",
"examples"),
15 "dr < 0.5 and abs(dz) < 2 and nCDCHits > 20 and thetaInCDCAcceptance"
18ma.fillParticleList(
"mu-",
"muonID > 0.9 and " + good_track, path=main)
19ma.fillParticleList(
"pi-",
"pionID > 0.5 and " + good_track, path=main)
21ma.reconstructDecay(
"B0:signal -> pi- mu+ ?nu", cut=
"", path=main)
24 "Upsilon(4S):opposite_cp -> B0:generic anti-B0:signal", cut=
"", path=main
27 decayString=
"Upsilon(4S):same_cp -> B0:generic B0:signal",
33 outputListName=
"Upsilon(4S)",
34 inputListNames=[
"Upsilon(4S):opposite_cp",
"Upsilon(4S):same_cp"],
38ma.buildRestOfEvent(
"Upsilon(4S)", path=main)
39track_based_cuts =
"thetaInCDCAcceptance and pt > 0.075 and dr < 2 and abs(dz) < 4"
40ecl_based_cuts =
"thetaInCDCAcceptance and E > 0.05"
41roe_mask = (
"my_mask", track_based_cuts, ecl_based_cuts)
42ma.appendROEMasks(
"Upsilon(4S)", [roe_mask], path=main)
44ma.matchMCTruth(list_name=
"Upsilon(4S)", path=main)
46vm.addAlias(
"Btag_SigProb",
"daughter(0, extraInfo(SignalProbability))")
47vm.addAlias(
"Btag_decayModeID",
"daughter(0, extraInfo(decayModeID))")
48vm.addAlias(
"Btag_Mbc",
"daughter(0, Mbc)")
50vm.addAlias(
"Bsig_isSignal",
51 "daughter(1, isSignal)")
52vm.addAlias(
"nCharged",
"nROE_Charged(my_mask)")
65 filename=
'Upsilon4S.root',
69b2.process(main, calculateStatistics=
True)