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