Пример #1
0
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
        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