def __configure(self):
recoSeq = GaudiSequencer("Hlt2JetRecoSequence", ModeOR = True, ShortCircuit = False)
from HltTracking.Hlt2TrackingConfigurations import Hlt2BiKalmanFittedForwardTracking
from HltTracking.Hlt2TrackingConfigurations import Hlt2BiKalmanFittedDownstreamTracking
# Long charged proto particles
longTracking = Hlt2BiKalmanFittedForwardTracking()
longChargedProtos = longTracking.hlt2ChargedNoPIDsProtos()
longSeq = GaudiSequencer("Hlt2JetLongRecoSequence", Members = longChargedProtos.members())
# Downstream charged proto particles
downTracking = Hlt2BiKalmanFittedDownstreamTracking()
downChargedProtos = downTracking.hlt2ChargedNoPIDsProtos()
downSeq = GaudiSequencer("Hlt2JetDownRecoSequence", Members = downChargedProtos.members())
from Hlt2SharedParticles.TrackFittedBasicParticles import BiKalmanFittedPhotons as Photons
from Hlt2SharedParticles.Pi0 import ResolvedPi0s
from Hlt2SharedParticles.Pi0 import MergedPi0s
from Hlt2SharedParticles.Ks import KsLLTF as KsLL
from Hlt2SharedParticles.Ks import KsLLTF as KsDD
from Hlt2SharedParticles.Lambda import LambdaLLTrackFitted as LambdaLL
from Hlt2SharedParticles.Lambda import LambdaDDTrackFitted as LambdaDD
self.__caloProcessor = longTracking.caloProcessor()
ecalSeq = self.__caloProcessor.clusters()
inputs = [
['Particle', 'particle', self.__sharedParticleInput(Photons, "Photons")],
['Particle', 'particle', self.__sharedParticleInput(ResolvedPi0s, "ResolvedPi0s")],
['Particle', 'particle', self.__sharedParticleInput(MergedPi0s, "MergedPi0s")],
['Particle', 'particle', self.__sharedParticleInput(KsLL, "KsLL")],
['Particle', 'particle', self.__sharedParticleInput(KsDD, "KsDD")],
['Particle', 'particle', self.__sharedParticleInput(LambdaLL, "LambdaLL")],
['Particle', 'particle', self.__sharedParticleInput(LambdaDD, "LambdaDD")],
['ProtoParticle', 'best', (longSeq, longChargedProtos.outputSelection())],
['ProtoParticle', 'best', (downSeq, downChargedProtos.outputSelection())],
['CaloCluster', 'gamma', (None, self.__findCaloLocation(Photons.members(), CellularAutomatonAlg, "OutputData"))],
['CaloCluster', 'gamma', self.__hcalClusters()],
['IClusTrTable2D', 'ecal', (None, self.__findCaloLocation(Photons.members(), PhotonMatchAlg, "Output"))],
['IClusTrTable2D', 'hcal', self.__hcal2Track(longTracking, Photons)]
]
pfInputs = []
for inputClass, inputType, (seq, loc) in inputs:
if seq and seq not in recoSeq.Members:
recoSeq.Members += [seq]
pfInputs.append([inputClass, inputType, loc])
from Configurables import HltParticleFlow
particleFlow = HltParticleFlow("HltParticleFlow", **self.getProp('ParticleFlowOptions'))
particleFlow.Inputs = pfInputs
particleFlow.Output = "Hlt2/ParticleFlow/Particles"
from HltLine.HltLine import bindMembers
self.__particleFlow = bindMembers("Hlt2Jet", [recoSeq, particleFlow])
genJB.Inputs = [genPF.Output]
genJB.Output = 'Phys/JB/MCParticles'
# Create the reconstructed jets.
from Configurables import HltParticleFlow, HltJetBuilder
from StandardParticles import (StdLooseKsDD, StdLooseKsLL, StdLooseKsLD,
StdLooseLambdaDD, StdLooseLambdaLL,
StdLooseLambdaLD)
recPF = HltParticleFlow('recPF')
recPF.Inputs = [
#['Particle', 'daughters', 'Phys/Z02MuMuLine/Particles'],
['Particle', 'daughters', Zs.outputLocation()],
['ProtoParticle', 'best', 'Rec/ProtoP/Charged'],
['ProtoParticle', 'gamma', 'Rec/ProtoP/Neutrals']
]
recPF.Output = 'Phys/PF/Particles'
recPF.ProBestNames = ['mu+', 'e+', 'p+', 'K+', 'pi+']
recPF.ProBestKeys = [701, 700, 704, 703, 702]
recPF.ProBestMins = [0.5, 0.5, 0.5, 0.5, 0.5]
recPF.EcalBest = True
recPF.SprRecover = False
recPF.TrkLnErrMax = 10
recPF.TrkUpErrMax = 10
recPF.TrkDnErrMax = 10
recJB = HltJetBuilder('recJB')
recJB.JetEcPath = ''
recJB.Inputs = [recPF.Output]
recJB.Output = 'Phys/JB/Particles'
recJB.JetPtMin = JetPtMin
#recSVs_seq = SelectionSequence('recSVs_Seq', TopSelection=recSVs)
pf.Inputs += [['ProtoParticle', 'best', 'Rec/ProtoP/Charged'],
['ProtoParticle', 'gamma', 'Rec/ProtoP/Neutrals']]
pf.Inputs += [
['Particle', 'daughters',
convLLSeq.outputLocation()],
['Particle', 'daughters',
convDDSeq.outputLocation()],
['Particle', 'daughters',
caloSeq.outputLocation()],
#['Particle', 'daughters', mergedpi0Seq.outputLocation()],
['Particle', 'daughters',
dimuSeq.outputLocation()]
]
pf.Output = 'Phys/PF/Particles'
pf.ProBestNames = ['mu+', 'e+', 'p+', 'K+', 'pi+']
pf.ProBestKeys = [701, 700, 704, 703, 702]
pf.ProBestMins = [0.5, 0.5, 0.5, 0.5, 0.5]
pf.EcalBest = True
pf.SprRecover = False
pf.TrkLnErrMax = 10
pf.TrkUpErrMax = 10
pf.TrkDnErrMax = 10
jb = HltJetBuilder('jb')
jb.JetEcPath = ''
jb.Inputs = [pf.Output]
jb.Output = 'Phys/JB/Particles'
jb.JetPtMin = 0
DaVinci().appendToMainSequence([pf, jb])