Skip to content
Snippets Groups Projects

AP for Lb->LcD0pi using 2016 data

Merged Chencheng Feng requested to merge fengc/Lb2LcD0pi into master
Files
2
+ 18
20
@@ -9,7 +9,7 @@ elif backend=="forAP":
###############################################################
######################## dv_opts.py ########################
######################## Lb2LcD0pi.py ########################
###############################################################
@@ -31,7 +31,7 @@ DaVinci().RootInTES = '/Event/Bhadron'
###############################################################
######################## generation.py ########################
######################## Lb2LcD0pi.py ########################
###############################################################
@@ -48,7 +48,7 @@ from Configurables import TupleToolMCTruth, TupleToolMCBackgroundInfo
from Configurables import LoKi__Hybrid__TupleTool, TupleToolVeto
from Configurables import LoKi__Hybrid__EvtTupleTool,LoKi__Hybrid__DictOfFunctors, LoKi__Hybrid__Dict2Tuple, LoKi__Hybrid__DTFDict
from Configurables import TupleToolDecayTreeFitter
#from Configurables import LoKi__Hybrid__DictofFunctors,LoKi__Hybrid__DTFDict,LoKi__Hybrid__DictValue, LoKi__Hybrid__Dict2Tuple
from DecayTreeTuple.Configuration import *
from Configurables import SubstitutePID
@@ -274,18 +274,16 @@ LcD0pi_line = 'Phys/X2LcD0PiD02KPiBeauty2CharmLine/Particles'
FilterLcD0pi = FilterDesktop("FilterLcD0pi")
DLcdaughter = { "D0" : OrderedDict({'K':1, 'pi':2}), "Lc" : OrderedDict({'p':1, 'K':2, 'pi':3}) }
FilterLcD0pi.Code = "((M > 5.1*GeV) & ( M < 6*GeV ))"
# (P>32*GeV) & (4*GeV<PT) & \
# (PT<80*GeV) & (2<ETA) & (ETA<5) &\
# (CHILD(P,1)>12*GeV) & (CHILD(P,2)>10*GeV) & (CHILD(P,3)>2.5*GeV) & (CHILD(PT,1)>1*GeV) & (CHILD(PT,2)>800*MeV) & (CHILD(PT,3)>150*MeV) &\
# (CHILD(CHILD(P,1),1)>4*GeV) & (CHILD(CHILD(P,2),1)>2.5*GeV) & (CHILD(CHILD(P,3),1)>1.5*GeV) &\
# (CHILD(CHILD(P,1),2)>2.5*GeV) & (CHILD(CHILD(P,2),2)>1.5*GeV) &\
# (CHILD(CHILD(PT,1),1)>300*MeV) & (CHILD(CHILD(PT,2),1)>150*MeV) & (CHILD(CHILD(PT,1),2)>150*MeV) &\
# (DOCA(1,2)<0.2*mm) & (DOCA(1,3)<0.25*mm) & (DOCA(2,3)<0.25*mm) &\
# (CHILD(DOCA(1,2),1)<0.3*mm) & (CHILD(DOCA(1,3),1)<0.4*mm) & (CHILD(DOCA(2,3),1)<0.4*mm) & (CHILD(DOCA(1,2),2)<0.4*mm) &\
# (CHILD(PROBNNk,3)>0.1) & (CHILD(CHILD(PROBNNk,1),2)*CHILD(CHILD(PROBNNpi,2),2)>0.01) & (BPVIPCHI2()<12) & (BPVVD>1.5*mm) &\
# (CHILD(BPVVD,1)>1.8*mm) & (CHILD(BPVVD,2)>2.0*mm) & ((CHILD(VFASPF(VZ),1) - VFASPF(VZ))>-3*mm) & ((CHILD(VFASPF(VZ),2) - VFASPF(VZ))>-2*mm)"
FilterLcD0pi.Code = "((M > 5.1*GeV) & ( M < 6*GeV ) &\
(P>32*GeV) & (4*GeV<PT) & \
(PT<80*GeV) & (2<ETA) & (ETA<5) &\
(DOCA(1,2)<0.2*mm) & (DOCA(1,3)<0.25*mm) & (DOCA(2,3)<0.25*mm) &\
(CHILD(DOCA(1,2),1)<0.4*mm) & (CHILD(DOCA(1,3),1)<0.4*mm) & (CHILD(DOCA(2,3),1)<0.4*mm) & (CHILD(DOCA(1,2),2)<0.4*mm) &\
(CHILD(PROBNNpi,3)>0.1) & (CHILD(CHILD(PROBNNk,1),2)*CHILD(CHILD(PROBNNpi,2),2)>0.01) & (BPVIPCHI2()<12) & (BPVVD>1.5*mm) &\
(CHILD(P,1)>12*GeV) & (CHILD(P,2)>10*GeV) & (CHILD(P,3)>2.5*GeV) & (CHILD(PT,1)>1*GeV) & (CHILD(PT,2)>800*MeV) & (CHILD(PT,3)>150*MeV) &\
(CHILD(CHILD(P,1),1)>4*GeV) & (CHILD(CHILD(P,2),1)>2.5*GeV) & (CHILD(CHILD(P,3),1)>1.5*GeV) &\
(CHILD(CHILD(P,1),2)>2.5*GeV) & (CHILD(CHILD(P,2),2)>1.5*GeV) &\
(CHILD(CHILD(PT,1),1)>300*MeV) & (CHILD(CHILD(PT,2),1)>150*MeV) & (CHILD(CHILD(PT,1),2)>150*MeV) )"
LcD0pi_sel = Selection("Sel_LcD0pi",
@@ -412,11 +410,11 @@ ResConfig(tuple_LcD0pi.B0)
ResConfig(tuple_LcD0pi.Lc)
ResConfig(tuple_LcD0pi.D0)
tuple_LcD0pi.pi.addTupleTool("LoKi::Hybrid::TupleTool/KLoKiTool").Variables = dau_vars()
tuple_LcD0pi.Lc_p.addTupleTool("LoKi::Hybrid::TupleTool/KLoKiTool").Variables = dau_vars()
tuple_LcD0pi.Lc_K.addTupleTool("LoKi::Hybrid::TupleTool/KLoKiTool").Variables = dau_vars()
tuple_LcD0pi.D0_K.addTupleTool("LoKi::Hybrid::TupleTool/KLoKiTool").Variables = dau_vars()
tuple_LcD0pi.D0_pi.addTupleTool("LoKi::Hybrid::TupleTool/KLoKiTool").Variables = dau_vars()
# tuple_LcD0pi.pi.addTupleTool("LoKi::Hybrid::TupleTool/KLoKiTool").Variables = dau_vars()
# tuple_LcD0pi.Lc_p.addTupleTool("LoKi::Hybrid::TupleTool/KLoKiTool").Variables = dau_vars()
# tuple_LcD0pi.Lc_K.addTupleTool("LoKi::Hybrid::TupleTool/KLoKiTool").Variables = dau_vars()
# tuple_LcD0pi.D0_K.addTupleTool("LoKi::Hybrid::TupleTool/KLoKiTool").Variables = dau_vars()
# tuple_LcD0pi.D0_pi.addTupleTool("LoKi::Hybrid::TupleTool/KLoKiTool").Variables = dau_vars()
TrackConfig(tuple_LcD0pi.pi)
TrackConfig(tuple_LcD0pi.Lc_p)
TrackConfig(tuple_LcD0pi.Lc_K)
Loading