'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
tupletools = []
Пример #2
0
    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
DaVinci().Lumi = not DaVinci().Simulation
DaVinci().EvtMax = 1000

IOHelper().inputFiles([('root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2016/ALLSTREAMS.DST/00062514/0000/00062514_00000008_7.AllStreams.dst')],
                      clear=True)
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
DaVinci().Lumi = not DaVinci().Simulation
DaVinci().EvtMax = 10000

# Use the local input data
IOHelper().inputFiles([('root://eoslhcb.cern.ch/'
                        '/eos/lhcb/grid/prod/lhcb/'
                        'MC/2012/ALLSTREAMS.DST/00035742/0000/'