DaughtersCuts={
'p+': '(PT > 200*MeV) & (P > 2000*MeV) & (MIPCHI2DV(PRIMARY) > 4)',
'K-': '(PT > 200*MeV) & (P > 2000*MeV) & (MIPCHI2DV(PRIMARY) > 4)',
'pi+': '(PT > 200*MeV) & (P > 2000*MeV) & (MIPCHI2DV(PRIMARY) > 4)'
},
CombinationCut=
"( (ADAMASS('Lambda_c+') < 110*MeV) | (ADAMASS('Xi_c+') < 110*MeV) )",
MotherCut=
"( (VFASPF(VCHI2/VDOF)< 10) & ( (ADMASS('Lambda_c+') < 80*MeV) | (ADMASS('Xi_c+') < 80*MeV) )"
)
from PhysConf.Selections import Selection, SelectionSequence
Lc2pKpi_sel = Selection('Lc2pKpi_sel',
Algorithm=Lc2pKpi_combiner,
RequiredSelections=[Pions, Kaons, Protons])
Lc2pKpi_seq = SelectionSequence('Lc2pKpi_seq', TopSelection=Lc2pKpi_sel)
DaVinci().UserAlgorithms += [Lc2pKpi_seq.sequence()]
tuple_Lc2pKpi.Inputs = [Lc2pKpi_sel.outputLocation()]
if (Turbo and year in ["2015", "2016"]):
tuple_Lc2pKpi.InputPrimaryVertices = '/Event/Turbo/Primary'
# (detached) B -> (Lc -> p K pi) mu nu
#line = "SelLc2PKPiforCharmFromBSemi"
#tuple_b2Lc2pKpi = DecayTreeTuple( 'tuple_b2Lc2pKpi' )
#tuple_b2Lc2pKpi.Inputs = ['Phys/{0}/Particles'.format(line)]
#tuple_b2Lc2pKpi.Decay = '[ Beauty -> ^( Lambda_c+ -> ^p ^K- ^pi+ ) ^mu- ]CC'
tuples = [tuple_Lc2pKpi]
# Define common tuple tools
dstar = CombineParticles(
'Combine_Dstar',
DecayDescriptor='[D*(2010)+ -> D0 pi+]cc',
DaughtersCuts=dstar_decay_products,
CombinationCut=dstar_comb,
MotherCut=dstar_mother
)
dstar_sel = Selection(
'Sel_Dstar',
Algorithm=dstar,
RequiredSelections=[d0_sel, Pions]
)
dstar_seq = SelectionSequence('Dstar_Seq', TopSelection=dstar_sel)
# Create an ntuple
dtt = DecayTreeTuple('TupleDstToD0pi_D0ToKK')
dtt.Inputs = dstar_seq.outputLocations()
dtt.Decay = '[D*(2010)+ -> ^D0 ^pi+]CC'
DaVinci().UserAlgorithms += [dstar_seq.sequence(), dtt]
# DaVinci configuration
DaVinci().InputType = 'DST'
DaVinci().TupleFile = 'DVntuple.root'
DaVinci().PrintFreq = 1000
DaVinci().DataType = '2016'
DaVinci().Simulation = True
# Only ask for luminosity information when not using simulated data
dstar_comb = "(ADAMASS('D*(2010)+') < 400*MeV)"
dstar_mother = (
"(abs(M-MAXTREE('D0'==ABSID,M)-145.42) < 10*MeV)"
'& (VFASPF(VCHI2/VDOF)< 9)'
)
dstar_sel = CombineSelection(
'Sel_Dstar',
[d0_sel, soft_pion_sel],
DecayDescriptor='[D*(2010)+ -> D0 pi+]cc',
CombinationCut=dstar_comb,
MotherCut=dstar_mother
)
dstar_seq = SelectionSequence('Dstar_Seq', TopSelection=dstar_sel)
# Create an ntuple
dtt = DecayTreeTuple('TupleDstToD0pi_D0ToKpi')
dtt.Inputs = dstar_seq.outputLocations()
dtt.Decay = '[D*(2010)+ -> D0 pi+]CC'
DaVinci().UserAlgorithms += [dstar_seq.sequence(), dtt]
# DaVinci configuration
DaVinci().InputType = 'DST'
DaVinci().TupleFile = 'DVntuple.root'
DaVinci().PrintFreq = 1000
DaVinci().DataType = '2012'
DaVinci().Simulation = True
# Only ask for luminosity information when not using simulated data