def __init__(self, name):
code = ("(MM > %(MinMass)s) & (PT> %(Pt)s)" +
" & (MINTREE('mu-' == ABSID, PT) > %(MuPt)s)" +
" & (VFASPF(VCHI2PDOF)<%(VertexChi2)s )" +
" & (MAXTREE('mu-' == ABSID, TRCHI2DOF) < %(TrChi2Tight)s)")
inputs = [TrackFittedDiMuon]
Hlt2ParticleFilter.__init__(self, name, code, inputs, shared = True)
def __init__(self, name):
code = ("(MM > %(MinMass)s) & (PT > %(Pt)s) " +
"& ~in_range(%(MinExcMass)s, MM, %(MaxExcMass)s) " +
"& (MAXTREE('mu-' == ABSID, TRCHI2DOF) < %(TkChi2)s) " +
"& (MINTREE('mu-' == ABSID, PT) > %(MuPt)s)")
inputs = [TrackFittedDiMuon]
Hlt2ParticleFilter.__init__(self, name, code, inputs)
def __init__(self, name, inputs, nickname = None):
cut = ("(MIPCHI2DV(PRIMARY) < %(BPVIPCHI2_MAX)s)")
if nickname is None:
nickname = name
Hlt2ParticleFilter.__init__(self, 'PID' + name + 'PromptIPChi2Filter', cut,
inputs, shared = True, nickname = nickname,
UseP2PVRelations = False)
def __init__(self, name):
code = ("(MINTREE('mu-' == ABSID, MIPCHI2DV(PRIMARY)) > %(IPChi2)s)" +
" & (BPVDLS > %(DLS)s)")
inputs = [DiMuonFilter('DiMuon')]
Hlt2ParticleFilter.__init__(self, name, code, inputs,
dependencies = [PV3D('Hlt2')],
UseP2PVRelations = False)
def __init__(self, tag, n, inputs, varmap, props):
pids = ['K+', 'KS0', 'Lambda0', 'Lambda~0']
if tag == 'RAD' and n == 3: pids.append('gamma')
varmap = PrepVarMap(pids, varmap)
params = ('$PARAMFILESROOT/data/Hlt2Topo%dBody_BDTParams_%s.txt' %
(n, props['RUN1_BDT_%iBODY_PARAMS' %
(2 if tag == 'RAD' else n)]))
bdttool = self.__classifier(params, varmap, "TrgBBDT")
pc = ("(VALUE('%s/%s') > %%(RUN1_BDT_%iBODY%s_MIN)s)" %
(bdttool.Type.getType(), bdttool.Name, n, tag))
if tag == 'MU': pc += " & INTREE(HASPROTO & HASMUON & ISMUON)"
elif tag == 'E':
pc += (" & INTREE(HASPROTO " "& (PIDe > %(RUN1_PIDE_MIN)s))")
from HltTracking.HltPVs import PV3D
if tag == 'RAD' and props['RUN1_RAD_TOS']:
Hlt2ParticleFilter.__init__(self,
'BDT',
pc,
inputs,
tistos='RUN1_TisTosSpec',
shared=False,
tools=[bdttool],
dependencies=[PV3D('Hlt2')])
else:
Hlt2ParticleFilter.__init__(self,
'BDT',
pc,
inputs,
shared=False,
tools=[bdttool],
dependencies=[PV3D('Hlt2')])
def __init__(self, name, inputs):
from HltTracking.HltPVs import PV3D
cut = "in_range(%(Sig_M_MIN_HHmumu)s, M, %(Sig_M_MAX_HHmumu)s)"
Hlt2ParticleFilter.__init__(self,
name,
cut,
inputs,
dependencies=[PV3D('Hlt2')])
def __init__(self, name, inputs):
code = ("(MINTREE('mu+'==ABSID,BPVIPCHI2()) > %(MuIPChi2)s)"
"& (MINTREE('mu+'==ABSID,PROBNNmu) > %(MuProbNNmu)s)"
"& (PT > %(PT)s)"
"& (HASVERTEX)"
"& (BPVVDCHI2 > %(FDChi2)s)")
Hlt2ParticleFilter.__init__(self, name, code, inputs, dependencies = [PV3D('Hlt2')])
def __init__(self, name, inputs, vetoes = { }, windows = { }, misid_windows = { }):
cuts = ["( ( WM(" + ids + ") > %(MVETO_" + hypo + "_GT)s ) | ( WM(" + ids + ") < %(MVETO_" + hypo + "_LT)s ) )" for hypo, ids in vetoes.iteritems()]
cuts += ["( WM(" + ids + ") > %(MVETO_" + hypo + "_LOW)s )" for hypo, ids in misid_windows.iteritems()]
cuts += ["(in_range(%(M_" + comb + "_LOW)s, MASS(" + var + "), %(M_" + comb + "_HIGH)s))" for comb, var in windows.iteritems()]
cut = (" & ".join(cuts))
Hlt2ParticleFilter.__init__(self, 'PID' + name + 'Veto', cut, inputs,
shared = True, nickname = name,
UseP2PVRelations = False)
def __init__(self, name):
code = ("(PT > %(Pt)s)")
inputs = [Pions, DownPions, Photons]
Hlt2ParticleFilter.__init__(self,
name,
code,
inputs,
dependencies=[PV3D('Hlt2')])
def __init__(self, name):
code = "BPVDLS > %(DLS)s"
from HltLine.Hlt2Monitoring import Hlt2Monitor, Hlt2MonitorMinMax
args = {'UseP2PVRelations' : False,
'PostMonitor' : Hlt2Monitor("M", "M(#mu#mu)", 3686, 200, 'M_out', nbins = 25)}
inputs = [Psi2SFilter('Psi2S')]
Hlt2ParticleFilter.__init__(self, name, code, inputs,
dependencies = [PV3D('Hlt2')],**args)
def __init__(self, name):
code = "(MM>%(MinMass)s) & (PT>%(Pt)s)"
from HltLine.Hlt2Monitoring import Hlt2Monitor, Hlt2MonitorMinMax
args = {'PreMonitor' : Hlt2Monitor("M", "M(#mu#mu)", 3097, 200, 'M_in', nbins = 25),
'PostMonitor' : Hlt2Monitor("M", "M(#mu#mu)", 90000, 10000, 'M_out', nbins=25)}
inputs = [TrackFittedDiMuon]
Hlt2ParticleFilter.__init__(self, name, code, inputs, **args)
def __init__(self, name):
cut = "(PT > %(mu_PTmin)s)"
from Inputs import Hlt2LooseMuons
inputs = [Hlt2LooseMuons]
Hlt2ParticleFilter.__init__(self,
name,
cut,
inputs,
dependencies=[TrackGEC(name)])
def __init__(self, name, inputs, pidVar = None, pidULim = False):
if pidVar:
pidCut = " & (%s %s %%(PID_LIM)s)" % (pidVar, ('<' if pidULim else '>'))
else:
pidCut = ""
cut = ("(TRGHOSTPROB < %(Trk_GhostProb_MAX)s) & (PT > %(Trk_PT_MIN)s)"
+ " & (MIPCHI2DV(PRIMARY) > %(Trk_MinIPChi2)s)" + pidCut)
Hlt2ParticleFilter.__init__(self, 'PID' + name + 'PIDGhostFilter', cut, inputs,
shared = True, nickname = name, UseP2PVRelations = False)
def __init__(self, name):
cut = "(PT > %(pi0_PTmin)s)"
from Inputs import MergedPi0s, ResolvedPi0s
inputs = [MergedPi0s, ResolvedPi0s]
Hlt2ParticleFilter.__init__(self,
name,
cut,
inputs,
dependencies=[TrackGEC(name)])
def __init__(self, name, inputs):
from HltLine.Hlt2Monitoring import Hlt2MonitorMinMax
cut = "(DTF_FUN ( M , True , 'D0' , %(Dst_MaxDTFCHI2NDF)s , %(Dst_MaxDTFM)s+1.0 ) < %(Dst_MaxDTFM)s)"
args = {'PreMonitor' : (Hlt2MonitorMinMax("DTF_FUN(M,True,'D0')","D^{*} M_{DTF}(D^0,PV) [MeV]", 2000,2040,"Dst_DTF_M_in",nbins=50)
+"&"+Hlt2MonitorMinMax("DTF_CHI2NDOF(True,'D0')","DecayTreeFitter D^{*} #chi^{2}/ndf", 0,20,"Dst_DTF_CHI2_in",nbins=50)),
'PostMonitor' : (Hlt2MonitorMinMax("DTF_FUN(M,True,'D0')","DecayTreeFitter D^{*} mass [MeV]", 2000,2040,"Dst_DTF_M_out",nbins=50)
+"&"+Hlt2MonitorMinMax("DTF_CHI2NDOF(True,'D0')","DecayTreeFitter D^{*} #chi^{2}/ndf", 0,20,"Dst_DTF_CHI2_out",nbins=50))
}
Hlt2ParticleFilter.__init__(self, name, cut, inputs,**args)
def __init__(self, name):
code = ("(MM > %(MinMass)s) " +
"& (MAXTREE('mu-' == ABSID, TRCHI2DOF) < %(TrChi2Tight)s) " +
"& (VFASPF(VCHI2PDOF) < %(VertexChi2)s) ")
from HltLine.Hlt2Monitoring import Hlt2Monitor, Hlt2MonitorMinMax
args = {'PreMonitor' : Hlt2Monitor("M", "M(#mu#mu)", 3097, 200, 'M_in', nbins = 25),
'PostMonitor' : Hlt2Monitor("M", "M(#mu#mu)", 5300, 700, 'M_out', nbins=25)}
inputs = [TrackFittedDiMuon]
Hlt2ParticleFilter.__init__(self, name, code, inputs, **args)
def __init__(self, name):
code = ("(MIPCHI2DV(PRIMARY)>%(RareIPChi2)s)"+
"& (PT> %(RarePtMin)s) & (PT < %(RarePtMax)s) & ( P < %(RarePMax)s )"+
"& (TRCHI2DOF<%(RareTrChi2Tight)s)"+
"& (PIDmu >%(muID)s)")
inputs = [ SingleMuonFilter("SharedSingleMuon") ]
Hlt2ParticleFilter.__init__(self, name, code, inputs,
dependencies = [PV3D('Hlt2')],
UseP2PVRelations = False)
def __init__(self, mode, nickname):
inputs = [LongAssocMuons]
cut = "(Q>0)" if mode == 1 else "(Q<0)"
Hlt2ParticleFilter.__init__(self,
'TrackEff' + nickname + str(mode) + 'Long',
cut,
inputs,
dependencies=[PV3D('Hlt2')],
nickname=nickname,
shared=True)
def __init__(self, name, inputs, pidVar=None, pidULim=False): # {
cut = (" (TRCHI2DOF < %(Trk_ALL_TRCHI2DOF_MAX)s )" +
"& (PT > %(Trk_ALL_PT_MIN)s)" + "& (P > %(Trk_ALL_P_MIN)s)")
if pidVar: # {
pidCut = "& (%s %s %%(PID_LIM)s)" % (pidVar,
('<' if pidULim else '>'))
cut = cut + pidCut
# }
Hlt2ParticleFilter.__init__(self, name, cut, inputs, shared=True)
def __init__(self, name, input, ismuon = False):
ismuonstr = " & (CHILDCUT( ISMUONLOOSE , 1) )" if ismuon else ""
cut = ("( ADWM( 'KS0' , WM( 'pi+' , 'pi-') ) > %(LambdaKsVeto)s ) & " +
"(ADMASS('Lambda0') < %(LambdaMWindow)s ) & " +
"(CHILDCUT( ( P > %(LambdaProtonP)s ) & (PT > %(LambdaProtonPT)s), 1) ) & " +
"(MIPCHI2DV(PRIMARY) < %(LambdaIPChi2Max)s )" +
ismuonstr) # ASSUME proton is 1st child
Hlt2ParticleFilter.__init__(self, 'PID' + name, cut, [ input ],
shared = True,
nickname = name,
UseP2PVRelations = False)
def __init__(self, name, inputs, highmass = False):
code = ''
if highmass: code = '(M > %(M)s) &'
code += ("(MINTREE('mu+'==ABSID,PT) > %(MuPT)s) & (MINTREE('mu+'==ABSID,P) > %(MuP)s)"
"& (MINTREE('mu+'==ABSID,BPVIPCHI2()) < %(MuIPChi2)s)"
"& (MINTREE('mu+'==ABSID,PROBNNmu) > %(MuProbNNmu)s)"
"& (PT > %(PT)s)"
"& (HASVERTEX)"
"& (BPVVDCHI2 < %(FDChi2)s)")
Hlt2ParticleFilter.__init__(self, name, code, inputs, dependencies = [PV3D('Hlt2')])
def __init__(self, name):
code = ("(MM > %(MinMass)s) & (PT> %(Pt)s)" +
" & (MINTREE('mu-' == ABSID, PT) > %(MuPt)s)" +
" & (VFASPF(VCHI2PDOF)<%(VertexChi2)s )" +
" & (MAXTREE('mu-' == ABSID, TRCHI2DOF) < %(TrChi2Tight)s)" +
" & (MINTREE('mu-' == ABSID, MIPCHI2DV(PRIMARY)) > %(IPChi2)s)" +
" & (BPVDLS > %(DLS)s)")
inputs = [TrackFittedDiMuon]
Hlt2ParticleFilter.__init__(self, name, code, inputs,
dependencies = [PV3D('Hlt2')],
UseP2PVRelations = False)
def __init__(self):
inputs = [Hlt2Muons]
code = ("(PT > %(MU_PT)s) & (PROBNNmu > %(MU_PROBNNMU)s)"
"& (TRGHOSTPROB < %(GHOSTPROB)s)")
Hlt2ParticleFilter.__init__(self,
'JetsMuon',
code,
inputs,
shared=True,
dependencies=[PV3D('Hlt2')])
def __init__(self, name):
cut = ("(PT > %(KS_PT_MIN)s)" + " & (BPVVDCHI2 > %(KS_FD_MIN)s)" +
" & CHILDCUT((TRGHOSTPROB < %(GHOSTPROB_MAX)s),1)" +
" & CHILDCUT((TRGHOSTPROB < %(GHOSTPROB_MAX)s),2)" +
" & (ADMASS('KS0') < %(KS_MASS_WINDOW)s)" +
" & (VFASPF(VCHI2PDOF) < %(KS_VCHI2NDOF_MAX)s)")
Hlt2ParticleFilter.__init__(self,
name,
cut,
inputs=[KsLLTF, KsDD],
shared=True)
def __init__(self, name, nickname, inputs,hadron):
dc = "(MIPCHI2DV(PRIMARY) > %(DTrk_ALL_MIPCHI2DV_MIN)s)"\
"& (PT > %(DTrk_ALL_PT_MIN)s)" \
"& (P > %(DTrk_ALL_P_MIN)s)"
if (hadron=='pi'):
dc = dc + "& (PIDK < %(Pi_PIDK)s)"
if (hadron=='K'):
dc = dc + "& (PIDK > %(K_PIDK)s)"
if (hadron=='p'):
dc = dc + "& (PIDp > %(P_PIDp)s)"
from HltTracking.HltPVs import PV3D
Hlt2ParticleFilter.__init__(self, name, dc, inputs,nickname = nickname, dependencies=[PV3D('Hlt2')],shared=True)
def __init__(self, inputs, gec):
pc = ("(PT > %(TRK_PT_MIN)s) & (P > %(TRK_P_MIN)s) "
"& (TRCHI2DOF < %(TRK_CHI2_MAX)s) "
"& (MIPCHI2DV(PRIMARY) > 16)")
from HltTracking.HltPVs import PV3D
deps = [FilterEvent(), PV3D('Hlt2')] if gec else [PV3D('Hlt2')]
Hlt2ParticleFilter.__init__(self,
'TopoInput4',
pc,
inputs,
shared=True,
dependencies=deps)
def __init__(self, name, inputs):
from HltTracking.HltPVs import PV3D
cut = "(TRCHI2DOF< %(TrkChi2_HLL)s )" \
"& (PT> %(TrkPt_HLL)s)" \
"& (P> %(TrkP_HLL)s)" \
"& (MIPCHI2DV(PRIMARY)> %(TrkPVIPChi2_HLL)s )"
Hlt2ParticleFilter.__init__(self,
name,
cut,
inputs,
dependencies=[PV3D('Hlt2')],
shared=True)
def __init__(self, name):
code = ("(ADMASS(3686.09) < %(MassWindow)s) " +
"& (MAXTREE('mu-' == ABSID, TRCHI2DOF) < %(TrChi2Tight)s) " +
"& (PT < %(PtMax)s) " +
"& (MINTREE('mu-' == ABSID, PT) > %(MuPt)s) " +
"& (VFASPF(VCHI2PDOF) < %(VertexChi2)s )")
from HltLine.Hlt2Monitoring import Hlt2Monitor, Hlt2MonitorMinMax
args = {'PreMonitor' : Hlt2Monitor("M", "M(#mu#mu)", 3097, 200, 'M_in', nbins = 25),
'PostMonitor' : Hlt2Monitor("M", "M(#mu#mu)", 3686, 200, 'M_out', nbins = 25)}
inputs = [TrackFittedDiMuon]
Hlt2ParticleFilter.__init__(self, name, code, inputs, shared = True, **args)
def __init__(self, name, inputs, nickname=None, shared=True):
pidp_cut = "(PIDp > %(PIDp_MIN)s )"
deltaPID_cut = "( (PIDp-PIDK) > %(DeltaPID_MIN)s )"
momentum_cut = "( P > %(P_MIN)s )"
cut = pidp_cut + " & " + deltaPID_cut + " & " + momentum_cut
nickname = name if nickname == None else nickname
name = name if not shared else 'B2HHHH%sTighterProton' % name
Hlt2ParticleFilter.__init__(self,
name,
cut,
inputs,
nickname=nickname,
shared=shared)
def __init__(self, tag, m, n, inputs):
from GaudiKernel.SystemOfUnits import MeV
if m == 2:
pc = ("(M < %f) & (VFASPF(VCHI2)"
" < %%(RUN1_V2BODYCHI2_MAX)s)" %
(5000 * MeV if n == 4 else 6000 * MeV))
else:
pc = '(M < %f)' % (6000 * MeV)
Hlt2ParticleFilter.__init__(self,
'Run1Topo%s%iFor%i' % (tag, m, n),
pc,
inputs,
shared=True)