RunType Upgrade and change of TraskListMerger to TrackContainerCopy

Changed DefaultRecoSequences for the Upgrade, added RunType Upgrade in Escher Configuration. Changed TrackListMerger into TrackContainerCopy in TAlignment Configuration. Added a working version of the option file for testing the alignment (AllignVP.py).

Alignment/Escher/options/AlignVP.py

+###############################################################################
+# (c) Copyright 2000-2018 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.                                       #
+###############################################################################
+##############################################################################
+# File for running VP alignment on simulated Data
+#
+# Syntax is:
+#
+#   gaudirun.py $ESCHEROPTS/AlignVeloSensors.py $ESCHEROPTS/2008-MC-Files.py
+#   gaudirun.py $ESCHEROPTS/AlignVeloSensors.py $ESCHEROPTS/DC06-Files.py
+#
+#   $ESCHEROPTS/gaudiiter.py -n 3 -e 1000 $ESCHEROPTS/AlignVeloSensors.py $ESCHEROPTS/2008-MC-Files.py
+#
+##############################################################################
+from Configurables import Escher, TrackSys, RecSysConf, LHCbApp, DstConf, RecMoniConf
+from GaudiKernel.ProcessJobOptions import importOptions
+from Configurables import CondDB
+from Configurables import GaudiSequencer, ChargedProtoParticleMaker
+#from Configurables import PrPixelTracking, PrPixelStoreClusters, PrLHCbID2MCParticle
+from Configurables import CaloDigitConf,CaloProcessor,GlobalRecoConf
+from Configurables import TAlignment,GetElementsToBeAligned
+from TAlignment.Alignables import * #Alignables  # *
+from TAlignment.TrackSelections import TrackRefiner
+from Configurables import AlignAlgorithm
+from Configurables import TrackContainerCopy, TrackSelector, TrackMonitor
+from Configurables import VPTrackSelector
+# from TrackFitter.ConfiguredFitters import *
+from Gaudi.Configuration import appendPostConfigAction, INFO, DEBUG, VERBOSE
+from TAlignment.SurveyConstraints import * #SurveyConstraints  # *
+from GaudiConf import IOHelper
+from Configurables import CaloProcessor, GlobalRecoConf  # , CaloDigitConf
+
+
+importOptions('$STDOPTS/PreloadUnits.opts')
+
+LHCbApp().Simulation = True
+LHCbApp().DataType = 'Upgrade'
+CondDB().Upgrade = True
+CondDB().Tags = {
+    "DDDB": 'HEAD',
+    "LHCBCOND": 'default',
+    "SIMCOND": 'HEAD',#'dd4hep',
+    "ONLINE": 'fake'
+}
+
+detectors = ['VP', 'UT', 'FT', 'Rich1Pmt', 'Rich2Pmt', 'Ecal', 'Hcal', 'Muon', 'Magnet', 'Tr']
+LHCbApp().Detectors = detectors
+
+#CondDB().LoadCALIBDB = 'HLT1'
+
+Escher().InputType = "DIGI"
+Escher().Kalman = True
+Escher().Millepede = False
+Escher().EvtMax = -1
+Escher().MoniSequence = []
+#Escher().RecoSequence = ["Decoding", "VP", "VELOPIX","Tr","Vertex","CALO","PROTO"]
+Escher().RecoSequence = ["Decoding", "TrFast", "TrBest", "VP","VELOPIX","CALO","PROTO"]
+#Escher().RecoSequence = ["Decoding", "TrFast", "TrBest", "RICH", "RICHFUTURE", "CALO", "MUON", "PROTO", "SUMMARY"]
+Escher().WithMC = True
+Escher().Simulation = True
+Escher().DataType = 'Upgrade'
+Escher().Upgrade = True
+#LHCbApp().OutputLevel = VERBOSE
+# RecSysConf().Detectors = detectors
+# RecSysConf().Simulation = True
+# RecMoniConf().Detectors = detectors
+
+# Escher().ExpertTracking = ["kalmanSmoother"]
+
+# TrackSys().TrackExtraInfoAlgorithms = []
+# TrackSys().DataType = 'Upgrade'
+# TrackSys().Simulation = True
+# TrackSys().WithMC = False #True
+#TrackSys().TrackTypes   = ["Velo"]
+# TrackSys().TrackPatRecAlgorithms = ["VP", "Forward","TsaSeed", "Match", "Downstream"]
+# DstConf().Detectors = detectors
+
+#hvpseq = GaudiSequencer("VPSeq")
+#prLHCbID2mc = PrLHCbID2MCParticle()
+#prPixTr = PrPixelTracking()
+#prPixStCl = PrPixelStoreClusters()
+
+#hvpseq.Members += [prPixStCl, prPixTr, prLHCbID2mc]
+
+AlignAlgorithm("Alignment").ForcedInitialTime = 1
+
+elements = Alignables()
+dofs = "TxTyTzRxRyRz"
+elements.VPModules(dofs) # Ty Tz Rx Ry Rz')
+#elements.VPLeft("TxTy")
+# elements.VPRight("TxTyRz")
+# constraints = ["Rz","Ry","Rx","Szx","Szy","SRz"]
+constraints = ["Tx", "Ty", "Szx", "Szy", "Rz", "SRz"]
+
+# constraints = []
+# constraints.append("VPFixModule10 : VP/VPLeft/Module10WithSupport: Tx Ty Tz Rx Ry Rz")
+# constraints.append("VPFixModule11 : VP/VPRight/Module11WithSupport: Tx Ty Tz Rx Ry Rz")
+# constraints.append("VPFixModule40 : VP/VPLeft/Module40WithSupport: Tx Ty Tz Rx Ry Rz")
+# constraints.append("VPFixModule41 : VP/VPRight/Module41WithSupport: Tx Ty Tz Rx Ry Rz")
+constraints.append("VPFixModule40 : VP/VPLeft/Module40WithSupport: Tx Ty Tz Rx Ry Rz")
+constraints.append("VPFixModule41 : VP/VPRight/Module41WithSupport: Tx Ty Tz Rx Ry Rz")
+constraints.append("VPFixModule50 : VP/VPLeft/Module50WithSupport: Tx Ty Tz Rx Ry Rz")
+constraints.append("VPFixModule51 : VP/VPRight/Module51WithSupport: Tx Ty Tz Rx Ry Rz")
+# constraints = []
+print("aligning elements", elements)
+print("ok 1")
+
+trackRefitSeq = GaudiSequencer("TrackRefitSeq")
+
+# create a track list for tracks with velo hits
+velotrackselector = TrackContainerCopy("TracksWithVeloHits",
+                                       inputLocations=["Rec/Track/Best"],
+                                       outputLocation="Rec/Track/TracksWithVeloHits",
+                                       Selector=TrackSelector())
+
+# refit the tracks in that list
+#velotrackrefitter = ConfiguredEventFitter(Name="TracksWithVeloHitsFitter",
+#                                          TracksInContainer="Rec/Track/TracksWithVeloHits",
+#                                          FieldOff=True)
+#velotrackrefitter.Fitter.MeasProvider.IgnoreIT = True
+#velotrackrefitter.Fitter.MeasProvider.IgnoreOT = True
+#velotrackrefitter.Fitter.MeasProvider.IgnoreTT = True
+#velotrackrefitter.Fitter.MeasProvider.IgnoreMuon = True
+#velotrackrefitter.Fitter.MeasProvider.IgnoreUT = True
+#velotrackrefitter.Fitter.MakeNodes = True
+#velotrackrefitter.Fitter.MakeMeasurements = True
+#
+#trackRefitSeq.Members += [velotrackselector, velotrackrefitter]
+trackRefitSeq.Members += [velotrackselector]
+
+
+def redefinefilters():
+    GaudiSequencer("RecoTrSeq").Members += [trackRefitSeq]
+
+
+appendPostConfigAction(redefinefilters)
+
+
+class myVPTracksOverlap(TrackRefiner):
+    def __init__(self, Name="VPOverlapTracks", InputLocation="Rec/Track/Best", Fitted=True):
+        TrackRefiner.__init__(self, Name=Name, InputLocation=InputLocation, Fitted=Fitted)
+
+    def configureSelector(self, a):
+        a.Selector = VPTrackSelector()
+        a.Selector.MinHitsASide = 1
+        a.Selector.MinHitsCSide = 1
+        a.Selector.TrackTypes = ["Velo", "Long"]
+
+
+class myVPTracks(TrackRefiner):
+    def __init__(self, Name="VPOverlapTracks", InputLocation="Rec/Track/Best", Fitted=True):
+        TrackRefiner.__init__(self, Name=Name, InputLocation=InputLocation, Fitted=Fitted)
+
+    def configureSelector( self, a ):
+
+        a.Selector = VPTrackSelector()
+        a.Selector.TrackTypes = ["Velo","Long"]
+
+    def configureSelector(self, a):
+        a.Selector = VPTrackSelector()
+        a.Selector.TrackTypes = ["Long"]
+        a.Selector.MinHits = 5
+        a.Selector.MinPCut = 5000
+        a.Selector.MaxPCut = 200000
+        a.Selector.MinPtCut = 200
+
+        if self._fitted:
+            a.Selector.MaxChi2Cut = 5
+            a.Selector.MaxChi2PerDoFMatch = 5
+            a.Selector.MaxChi2PerDoFVelo = 5
+            a.Selector.MaxChi2PerDoFDownstream = 5
+
+
+class VPBackwards(TrackRefiner):
+    def __init__(self, Name="VPTracks", InputLocation="Rec/Track/Best", Fitted=True):
+        TrackRefiner.__init__(self, Name=Name, InputLocation=InputLocation, Fitted=Fitted)
+
+    def configureSelector(self, a):
+        a.Selector = VPTrackSelector()
+        a.Selector.TrackTypes = ["Backward"]
+
+        a.Selector.MinHits = 5
+        if self._fitted:
+            a.Selector.MaxChi2Cut = 5
+
+TAlignment().ElementsToAlign = list(elements)
+TAlignment().TrackLocation = "Rec/Track/Best" 
+TAlignment().VertexLocation = "Rec/Vertex/Primary"
+TAlignment().Constraints = constraints
+TAlignment().TrackSelections = [myVPTracks()]#, VPBackwards()]
+TAlignment().WriteCondSubDetList = ['VP']
+TAlignment().EigenValueThreshold = 100
+TAlignment().UsePreconditioning = True
+TAlignment().UseLocalFrame = True
+# include survey constraints
+
+surveyconstraints = SurveyConstraints()
+surveyconstraints.VP()
+
+
+# Convert Calo ReadoutStatus to ProcStatus
+#caloBanks = GaudiSequencer("CaloBanksHandler")
+#caloDetectors = [det for det in ['Spd', 'Prs', 'Ecal', 'Hcal'] if det in detectors]
+#CaloDigitConf(ReadoutStatusConvert=True, Sequence=caloBanks, Detectors=caloDetectors)
+#GaudiSequencer("EscherSequencer").Members += [caloBanks]
+
+CaloProcessor().setProp("NoSpdPrs", True)
+GlobalRecoConf().setProp("NoSpdPrs", True)
+
+
+def doMyChanges():
+    GaudiSequencer("HltFilterSeq").Members = []
+
+
+appendPostConfigAction(doMyChanges)
+
+IOHelper('ROOT').inputFiles([
+    '/eos/user/b/bimitres/Velo_alignment/Alignment/Boole.xdigi',
+])