Skip to content
Snippets Groups Projects

Lines for Dark Matter searches (from Lambda decays)

Merged Lines for Dark Matter searches (from Lambda decays)
StrippingBDecaysDM_QEE into 2017-patches
Merged Saul Lopez Solino requested to merge StrippingBDecaysDM_QEE into 2017-patches
Compare
and
1 file
+ 279
− 0
Compare changes
  • Side-by-side
  • Inline
Phys/StrippingSelections/python/StrippingSelections/StrippingQEE/StrippingLambdaDecaysDM.py 0 → 100644
+ 279
− 0
###############################################################################
# (c) Copyright 2000-2021 CERN for the benefit of the LHCb Collaboration #
# #
# This software is distributed under the terms of the GNU General Public #
# Licence version 3 (GPL Version 3), copied verbatim in the file "COPYING". #
# #
# In applying this licence, CERN does not waive the privileges and immunities #
# granted to it by virtue of its status as an Intergovernmental Organization #
# or submit itself to any jurisdiction. #
###############################################################################
"""
Lines for searches on DM candidates produced from Lambda decays (arXiv:2101.02706).
"""
### Report using 2017 data (163k events):
# StrippingReport INFO Event 163109, Good event 163109
# | *Decision name*|*Rate,%*|*Accepted*| *Mult*|*ms/evt*|
# |_StrippingGlobal_ | 0.1600| 261| | 8.450|
# |_StrippingSequenceStreamBhadronCompleteEvent_ | 0.1600| 261| | 8.428|
# |!StrippingLambdaDecaysDMLambda1520Line | 0.0417| 68| 1.206| 5.608|
# |!StrippingLambdaDecaysDMLambda1520ControlLine | 0.0313| 51| 1.118| 0.069|
# |!StrippingLambdaDecaysDMLambda2595Line | 0.0215| 35| 1.171| 0.884|
# |!StrippingLambdaDecaysDMLambda2595ControlLine | 0.0454| 74| 1.432| 0.217|
# |!StrippingLambdaDecaysDMLambda1520Line_TIMING | 0.0417| 68| 1.206| 0.268|
# |!StrippingLambdaDecaysDMLambda1520ControlLine_TIMING | 0.0325| 53| 1.189| 0.062|
# |!StrippingLambdaDecaysDMLambda2595Line_TIMING | 0.0215| 35| 1.171| 0.062|
# |!StrippingLambdaDecaysDMLambda2595ControlLine_TIMING | 0.0484| 79| 1.544| 0.120|
from Gaudi.Configuration import *
from LHCbKernel.Configuration import *
from CommonParticles.Utils import *
from CommonParticles import StdAllNoPIDsPions
from StandardParticles import StdAllLooseKaons, StdAllLooseProtons, StdAllLoosePions
from GaudiConfUtils.ConfigurableGenerators import CombineParticles
from PhysSelPython.Wrappers import Selection
from StrippingConf.StrippingLine import StrippingLine
from StrippingUtils.Utils import LineBuilder
from GaudiKernel.SystemOfUnits import MeV, picosecond, millimeter
__author__ = ['Xabier Cid Vidal', 'Carlos Vazquez Sierra', 'Saul Lopez Solino']
__date__ = '18/03/2021'
__all__ = 'LambdaDecaysDMConf', 'default_config'
default_config = {
'NAME' : 'LambdaDecaysDM',
'BUILDERTYPE' : 'LambdaDecaysDMConf',
'WGs' : [ 'QEE' ],
'STREAMS' : [ 'BhadronCompleteEvent' ],
'CONFIG' : {
'Common': {
'checkPV' : False
},
'StrippingDMLambda1520Line': {
'Prescale' : 1.0,
'Postscale' : 1.0,
'Prescale_control' : 0.1,
'Postscale_control' : 1.0,
'Lambda1520_PT' : 1000*MeV,
'Lambda1520_iso' : 0.85,
'Lambda1520_ProbNNK' : 0.5,
'Lambda1520_ProbNNp' : 0.8,
'Lambda1520_Daug_PT' : 500*MeV,
'Lambda1520_GhostProb': 0.3,
'Lambda1520_IPChi2' : 25,
'Lambda1520_VChi2' : 2,
'Lambda1520_M_Min' : 1460,
'Lambda1520_M_Max' : 1580,
'Lambda1520_FDChi2' : 45,
'Lambda1520_MaxDoca' : 0.1*millimeter,
},
'StrippingDMLambda2595Line': {
'Prescale' : 1.0,
'Postscale' : 1.0,
'Prescale_control' : 0.6,
'Postscale_control' : 1.0,
'K_ProbNNghost' : 0.1,
'K_ProbNNk' : 0.5,
'K_PT' : 250*MeV,
'p_ProbNNghost' : 0.1,
'p_ProbNNp' : 0.5,
'p_PT' : 250*MeV,
'pi_lambda2595_IPchi2' : 0.5,
'pi_lambda2595_ProbNNghost' : 0.1,
'pi_lambda2595_PT' : 100*MeV,
'pi_lambdac_ProbNNghost' : 0.1,
'pi_lambdac_PT' : 250*MeV,
'Lambda_DOCA' : 1.0*millimeter,
'Lambda_FDCHI2' : 50,
'Lambda_M_Min' : 2236*MeV,
'Lambda_M_Max' : 2336*MeV,
'Lambda_PT' : 1500*MeV,
'Lambda_VCHI2' : 4,
'Lambda2595_DOCA' : 1.0*millimeter,
'Lambda2595_FDCHI2' : 5,
'Lambda2595_ISO' : 0.6,
'Lambda2595_M_Min' : 2545*MeV,
'Lambda2595_M_Max' : 2645*MeV,
'Lambda2595_PT' : 1500*MeV,
'Lambda2595_VCHI2' : 4
},
}
}
class LambdaDecaysDMConf(LineBuilder) :
__configuration_keys__ = default_config['CONFIG'].keys()
def __init__(self, name, config):
LineBuilder.__init__(self, name, config)
sel_lambda1520 = self.combineLambda1520(name, config['StrippingDMLambda1520Line'])
lambda1520line = StrippingLine(name + 'Lambda1520Line',
prescale = config['StrippingDMLambda1520Line']['Prescale'],
postscale = config['StrippingDMLambda1520Line']['Postscale'],
checkPV = config['Common']['checkPV'],
selection = sel_lambda1520)
self.registerLine(lambda1520line)
sel_lambda1520_control = self.combineLambda1520(name+'_Control', config['StrippingDMLambda1520Line'], isolation = False)
lambda1520line_control = StrippingLine(name + 'Lambda1520ControlLine',
prescale = config['StrippingDMLambda1520Line']['Prescale_control'],
postscale = config['StrippingDMLambda1520Line']['Postscale_control'],
checkPV = config['Common']['checkPV'],
selection = sel_lambda1520_control)
self.registerLine(lambda1520line_control)
sel_lambda2595 = self.combineLambda2595(name, config['StrippingDMLambda2595Line'], isolation = True)
lambda2595line = StrippingLine(name + 'Lambda2595Line',
prescale = config['StrippingDMLambda2595Line']['Prescale'],
postscale = config['StrippingDMLambda2595Line']['Postscale'],
checkPV = config['Common']['checkPV'],
selection = sel_lambda2595)
self.registerLine(lambda2595line)
sel_lambda2595_control = self.combineLambda2595(name+'_control', config['StrippingDMLambda2595Line'], isolation = False)
lambda2595line_control = StrippingLine(name + 'Lambda2595ControlLine',
prescale = config['StrippingDMLambda2595Line']['Prescale_control'],
postscale = config['StrippingDMLambda2595Line']['Postscale_control'],
checkPV = config['Common']['checkPV'],
selection = sel_lambda2595_control)
self.registerLine(lambda2595line_control)
def combineLambda1520(self, name, config, isolation=True):
daugh_cut={'K': ("(PT > {Lambda1520_Daug_PT}) "
"& (MIPCHI2DV(PRIMARY) > {Lambda1520_IPChi2}) "
"& (TRGHOSTPROB < {Lambda1520_GhostProb}) "
"& (PROBNNK > {Lambda1520_ProbNNK}) ").format(**config),
'p': ("(PT > {Lambda1520_Daug_PT}) "
"& (MIPCHI2DV(PRIMARY) > {Lambda1520_IPChi2}) "
"& (TRGHOSTPROB < {Lambda1520_GhostProb}) "
"& (PROBNNp > {Lambda1520_ProbNNp}) ").format(**config)}
comb_cut = ("(APT > {Lambda1520_PT}) &"
"(AM > {Lambda1520_M_Min}) & (AM < {Lambda1520_M_Max}) &"
"(AMAXDOCA('') < {Lambda1520_MaxDoca}) ".format(**config))
if not isolation:
mother_cut = ("(PT > {Lambda1520_PT})"
"& (M > {Lambda1520_M_Min}) & (M < {Lambda1520_M_Max})"
"& (HASVERTEX)"
"& (VFASPF(VCHI2PDOF) < {Lambda1520_VChi2}) "
"& (BPVIPCHI2() > {Lambda1520_IPChi2})"
"& (BPVVDCHI2 > {Lambda1520_FDChi2})".format(**config))
else:
mother_cut = ("(PT > {Lambda1520_PT})"
"& (M > {Lambda1520_M_Min}) & (M < {Lambda1520_M_Max})"
"& (HASVERTEX)"
"& (VFASPF(VCHI2PDOF) < {Lambda1520_VChi2}) "
"& (BPVIPCHI2() > {Lambda1520_IPChi2})"
"& (BPVVDCHI2 > {Lambda1520_FDChi2})"
"& ((PT/(PT+PTCONE)) > {Lambda1520_iso})".format(**config))
decay = "[Lambda(1520)0 -> K- p+]cc"
namesel = name + "Lambda1520Sel"
_combination = CombineParticles( DecayDescriptor = decay,
CombinationCut = comb_cut,
DaughtersCuts = { 'K-' : daugh_cut["K"],
'p+' : daugh_cut["p"] },
MotherCut = mother_cut,
Preambulo = ["PTCONE = SUMCONE ( 0.6**2 , PT , '/Event/Phys/StdAllNoPIDsPions/Particles')"]
)
return Selection ( namesel,
Algorithm = _combination,
RequiredSelections = [ StdAllLooseProtons, StdAllLooseKaons ]
)
def combineLambda_c(self, name, config):
daugh_cut={'K': ("(PT > {K_PT} )"
"& (TRGHOSTPROB < {K_ProbNNghost})"
"& (PROBNNK > {K_ProbNNk})").format(**config),
'p': ("(PT > {p_PT} )"
"& (TRGHOSTPROB < {p_ProbNNghost})"
"& (PROBNNp > {p_ProbNNp})").format(**config),
'pi': ("(PT > {pi_lambdac_PT} )"
"& (TRGHOSTPROB < {pi_lambdac_ProbNNghost})").format(**config)}
comb_cut = ("(APT > {Lambda_PT}) &"
"(AM > {Lambda_M_Min}) & (AM < {Lambda_M_Max}) &"
"(AMAXDOCA('') < {Lambda_DOCA})".format(**config))
mother_cut = ("(PT > {Lambda_PT} ) &"
"( M > {Lambda_M_Min}) & (M < {Lambda_M_Max}) &"
"(VFASPF(VCHI2PDOF) < {Lambda_VCHI2}) &"
"(BPVVDCHI2 > {Lambda_FDCHI2})".format(**config))
decay = "[Lambda_c+ -> K- p+ pi+]cc"
namesel = name+"Lambdac_Lambda2595Sel"
_combination = CombineParticles( DecayDescriptor = decay,
CombinationCut = comb_cut,
DaughtersCuts = { 'K-' : daugh_cut["K"],
'p+' : daugh_cut["p"],
'pi+': daugh_cut["pi"] },
MotherCut = mother_cut
)
return Selection ( namesel,
Algorithm = _combination,
RequiredSelections = [ StdAllLooseKaons , StdAllLooseProtons, StdAllLoosePions ]
)
def combineLambda2595(self, name, config, isolation = True):
SelLambda_c = self.combineLambda_c("Lambda_c"+name, config)
decay = "[Lambda_c(2595)+ -> Lambda_c+ pi+ pi-]cc"
daugh_cut={'pi1': ("(PT > {pi_lambda2595_PT} )"
"& (MIPCHI2DV(PRIMARY) > {pi_lambda2595_IPchi2})"
"& (TRGHOSTPROB < {pi_lambda2595_ProbNNghost})").format(**config),
'pi2': ("(PT > {pi_lambda2595_PT} )"
"& (MIPCHI2DV(PRIMARY) > {pi_lambda2595_IPchi2})"
"& (TRGHOSTPROB < {pi_lambda2595_ProbNNghost})").format(**config)}
comb_cut = ("(APT > {Lambda2595_PT} ) &"
"(AM > {Lambda2595_M_Min}) & (AM < {Lambda2595_M_Max}) &"
"(AMAXDOCA('') < {Lambda2595_DOCA})".format(**config))
if not isolation:
mother_cut = ("(PT > {Lambda2595_PT} ) &"
"( M > {Lambda2595_M_Min}) & (M < {Lambda2595_M_Max}) &"
" (VFASPF(VCHI2PDOF) < {Lambda2595_VCHI2}) & "
" (BPVVDCHI2 > {Lambda2595_FDCHI2})".format(**config))
else:
mother_cut = ("(PT > {Lambda2595_PT} ) &"
"( M > {Lambda2595_M_Min}) & (M < {Lambda2595_M_Max}) &"
" (VFASPF(VCHI2PDOF) < {Lambda2595_VCHI2}) & "
" (BPVVDCHI2 > {Lambda2595_FDCHI2}) & "
" ((PT/(PT+PTCONE)) > {Lambda2595_ISO})".format(**config))
Lambda2595 = CombineParticles( DecayDescriptor = decay,
CombinationCut = comb_cut,
DaughtersCuts = { 'pi+' : daugh_cut["pi1"],
'pi-' : daugh_cut["pi2"]},
MotherCut = mother_cut,
Preambulo = ["PTCONE = SUMCONE ( 0.6**2 , PT , '/Event/Phys/StdAllNoPIDsPions/Particles')"]
)
return Selection("SelLambda2595"+name,
Algorithm = Lambda2595,
RequiredSelections = [SelLambda_c, StdAllLoosePions]
)