tup.addTool(MCTupleToolKinematic())
tup.MCTupleToolKinematic.Verbose = True
tup.addTool(LoKi_Photos)
tup.ToolList = [
"MCTupleToolHierarchy",
"MCTupleToolKinematic",
"LoKi::Hybrid::MCTupleTool/LoKi_Photos" # doesn't work with DaVinci v36r6
]
tup.addTool(TupleToolMCTruth, name="TruthTool")
tup.addTool(TupleToolMCBackgroundInfo, name="BackgroundInfo")
tup.ToolList += ["TupleToolMCTruth/TruthTool"]
tup.ToolList += ["TupleToolMCBackgroundInfo/BackgroundInfo"]
if OUTPUTLEVEL == DEBUG:
from Configurables import PrintMCTree, PrintMCDecayTreeTool
mctree = PrintMCTree("PrintTrue")
mctree.addTool(PrintMCDecayTreeTool)
mctree.PrintMCDecayTreeTool.Information = "Name M P Px Py Pz Pt Vx Vy Vz"
mctree.ParticleNames = ["B+", "B-"]
mctree.Depth = 3 # down to the K and mu
#########################################################################################################
# Now set up the DecayTreeTuples for the reconstructed particles
#########################################################################################################
tupletools = []
tupletools.append("TupleToolPrimaries")
tupletools.append("TupleToolKinematic")
tupletools.append("TupleToolGeometry")
tupletools.append("TupleToolTrackInfo")
tupletools.append("TupleToolPid")
tupletools.append("TupleToolRecoStats")
from Configurables import MCDecayTreeTuple
mcTuple = MCDecayTreeTuple("MCTuple")
mcTuple.Decay = "( [ [B0]nos => ^( J/psi(1S) => ^mu+ ^mu- ) ^( KS0 => ^pi+ ^pi- ) ]CC ) || ( [ [B0]os -> ^( J/psi(1S) => ^mu+ ^mu- ) ^( KS0 => ^pi+ ^pi- ) ]CC )"
mcTuple.ToolList = [
"MCTupleToolKinematic", "TupleToolEventInfo", "MCTupleToolReconstructed"
]
from Configurables import MCTupleToolReconstructed, MCReconstructed
#ToolSvc().addTool(MCReconstructed)
#ToolSvc().MCReconstructed.OutputLevel = 1
#mcTuple.OutputLevel = 1
from Configurables import PrintMCTree
pmc = PrintMCTree()
pmc.ParticleNames = ["B0", "B~0"]
########################################################################
#
# DaVinci
#
from Configurables import DaVinci
DaVinci().EvtMax = 100
DaVinci().SkipEvents = 0
DaVinci().DataType = "2010"
DaVinci().Simulation = True
#DaVinci().CondDBtag = 'MC09-20100827-vc-md100'
DaVinci().TupleFile = "DecayTreeTuple.root" # Ntuple
DaVinci().MoniSequence = [tuple, evtTuple, mcTuple] # , pmc ]
DaVinci().appendToMainSequence([sc.sequence()])
for b in bad :
del algs[b]
dod.AlgMap = algs
from Configurables import EventClockSvc, CondDB
EventClockSvc ( EventTimeDecoder = "FakeEventTime" )
CondDB ( IgnoreHeartBeat = True )
appendPostConfigAction( doIt)
mctl=[ 'MCTupleToolHierarchy', 'MCTupleToolKinematic', 'MCTupleToolPrimaries']
#mctl=[ 'MCTupleToolKinematic']
mct.ToolList=mctl
#mctB.ToolList=mctl
printMC = PrintMCTree()
printMC.ParticleNames = ["Lambda_c+","Lambda_c~-"]
printMC.ParticleNames = ["D0","D~0"]
printMC.ParticleNames = ["Xi_cc++","Xi_cc~--"]
printMC.ParticleNames = ["B0","B~0"]
printMC.ParticleNames = ["B+","B-"]
########################################################################
from Configurables import DaVinci
DaVinci().EvtMax = -1
#DaVinci().PrintFreq = 1
#DaVinci().SkipEvents = 4 # Events to skip
DaVinci().DataType = "2018"
#######################################################################
#
# Print Reconstructed Bs
#
from Configurables import PrintDecayTree, PrintDecayTreeTool
tree = PrintDecayTree("PrintFoundBs")
tree.Inputs = [SeqBs2JpsiPhi.outputLocation()]
tree.addTool(PrintDecayTreeTool, name="PrintDecay")
tree.PrintDecay.Information = "Name M P Px Py Pz Pt chi2"
DaVinci().MoniSequence = [tree] # The monitoring stuff
#######################################################################
#
# Print Print All True Bs (only on MC)
#
from Configurables import PrintMCTree, PrintMCDecayTreeTool
mctree = PrintMCTree("PrintTrueBs")
mctree.addTool(PrintMCDecayTreeTool)
mctree.PrintMCDecayTreeTool.Information = "Name M P Px Py Pz Pt chi2"
mctree.ParticleNames = ["B_s0", "B_s~0"]
mctree.Depth = 2 # down to the K and mu
DaVinci().MoniSequence += [mctree] # The monitoring stuff
#######################################################################
#
# Configure the application
#
DaVinci().UserAlgorithms = [seq] # The selection sequence
DaVinci().HistogramFile = "DVHistos_5.root" # Histogram file
DaVinci().EvtMax = 1000 # Number of events
DaVinci().DataType = "2012" #
DaVinci().Simulation = True
#
