4 from basf2
import conditions
as b2conditions
7 import simulation
as sim
9 import modularAnalysis
as ma
32 bkgFiles = glob.glob(
'/sw/belle2/bkg/*.root')
34 simulateJitter =
False
37 eventinfosetter = register_module(
'EventInfoSetter')
38 eventinfosetter.param(
'expList', expList)
39 eventinfosetter.param(
'runList', [0])
40 eventinfosetter.param(
'evtNumList', [numEvents])
41 main.add_module(eventinfosetter)
43 main.add_module(
'EvtGenInput')
48 forceSetPXDDataReduction=
True,
49 usePXDDataReduction=ROIfinding,
50 simulateT0jitter=simulateJitter)
54 b2conditions.override_globaltags()
55 b2conditions.globaltags = [
"online"]
58 main.add_module(
'RootInput', branchNames=[
'RawPXDs',
'RawSVDs',
'RawCDCs'])
59 raw.add_unpackers(main, components=[
'PXD',
'SVD',
'CDC'])
68 trk.add_tracking_reconstruction(
70 mcTrackFinding=MCTracking,
71 trackFitHypotheses=[211])
76 ma.applyEventCuts("nTracks ==2", path=main)
78 mySelection = 'pt>1.0 and abs(dz)<0.5 and dr<0.4'
79 ma.fillParticleList('mu+:DQM', mySelection, path=main)
80 ma.reconstructDecay('Upsilon(4S):IPDQM -> mu+:DQM mu-:DQM', '10<M<11', path=main)
82 skimfilter = register_module('SkimFilter')
83 skimfilter.set_name('SkimFilter_MUMU')
84 skimfilter.param('particleLists', ['Upsilon(4S):IPDQM'])
85 main.add_module(skimfilter)
86 filter_path = create_path()
87 skimfilter.if_value('=1', filter_path, AfterConditionPath.CONTINUE)
91 main.add_module(
'SVDPerformanceTTree', outputFileName=
"SVDPerformanceTree"+str(tag)+
".root")
94 main.add_module(
'OverlapResiduals', ExpertLevel=
True)
97 main.add_module(
'Progress')