Skip to content
Snippets Groups Projects
Commit 2501f873 authored by Jan Kuechler's avatar Jan Kuechler
Browse files

Add tools/algs from TrigFTKFastSimTruth_jobOptions to FTKRec

This requires that FTKRec is scheduled earlier, thus moving it from SystemRec_config to RecExCommon before the trigger
parent c409ab36
No related branches found
No related tags found
No related merge requests found
......@@ -639,6 +639,10 @@ if recAlgs.doAtlfast():
protectedInclude ("AtlfastAlgs/Atlfast_RecExCommon_Fragment.py")
AODFix_postAtlfast()
# functionality : FTK reconstruction
if DetFlags.detdescr.FTK_on() :
protectedInclude( "FTK_RecExample/FTKRec_jobOptions.py")
pdr.flag_domain('trig')
# no trigger, if readESD _and_ doESD ! (from Simon George, #87654)
......
......@@ -24,10 +24,6 @@ pdr.flag_domain('id')
if DetFlags.detdescr.ID_on():
protectedInclude( "InDetRecExample/InDetRec_jobOptions.py" )
AODFix_postInDetRec()
# functionality : FTK reconstruction
if DetFlags.detdescr.FTK_on() :
protectedInclude( "FTK_RecExample/FTKRec_jobOptions.py")
# functionality : Muon reconstruction
......
......@@ -23,6 +23,65 @@ if rec.doFTK() and (globalflags.InputFormat() == 'bytestream' or rec.readRDO):
theFTKrandomSvc.EventReseeding = True
theFTKrandomSvc.UseOldBrokenSeeding = False
ServiceMgr += theFTKrandomSvc
elif FTKFlags.doTruthSmearing:
from PyJobTransforms.trfUtils import findFile
pmap8L_path = findFile(os.environ['DATAPATH'], 'ftk_configuration/map_files/raw_8LcIbl3D123.pmap')
pmap12L_path = findFile(os.environ['DATAPATH'], 'ftk_configuration/map_files/raw_12LiblHW3D.pmap')
print "Using PMAP 8L:", pmap8L_path
print "Using PMAP 12L:", pmap12L_path
rmap_path = findFile(os.environ['DATAPATH'], 'ftk_configuration/map_files/raw_12Libl64TmodB_3D_t13.tmap')
print "Using RMAP:", rmap_path
import glob
sectors = []
ntower = 64
for i in range(ntower):
sectors += glob.glob("/cvmfs/atlas.cern.ch/repo/sw/database/GroupData/FTK/FitConstants/ftk.64tower.2016.mc16.v2/sectors_raw_12L_reg"+str(i)+".patt.bz2")
from TrigFTKSim.TrigFTKSimConf import FTK_SGHitInput
FTKSGInput = FTK_SGHitInput( maxEta= 3.2, minPt= 0.8*GeV)
FTKSGInput.ReadTruthTracks = False
FTKSGInput.DoOutFileRawHits=False
ToolSvc += FTKSGInput
from TrigFTKFastSimTruth.TrigFTKFastSimTruthConf import TrigFTKSectorMatchTool
FTKTrajectorSelector = TrigFTKSectorMatchTool()
FTKTrajectorSelector.PMapPath8L = pmap8L_path
FTKTrajectorSelector.PMapPath12L = pmap12L_path
FTKTrajectorSelector.RMapPath = rmap_path
FTKTrajectorSelector.SectorFiles = sectors
ToolSvc+=FTKTrajectorSelector
from TrigFTKFastSimTruth.TrigFTKFastSimTruthConf import TrigFTKTrackBuilderTool
FTKTruthTrackBuilder = TrigFTKTrackBuilderTool()
FTKTruthTrackBuilder.PMapPath = pmap12L_path
ToolSvc+=FTKTruthTrackBuilder
from AthenaServices.AthenaServicesConf import AtRanluxGenSvc
theFTKrandomSvc = AtRanluxGenSvc("TrigFTKFastSimTruthRandom")
theFTKrandomSvc.EventReseeding = True
theFTKrandomSvc.UseOldBrokenSeeding = False
ServiceMgr += theFTKrandomSvc
from TrigFTKFastSimTruth.TrigFTKFastSimTruthConf import TrigFTKFastSimTruth
theTrigFTKFastSimTruth = TrigFTKFastSimTruth(OutputLevel = INFO,
PMapPath_8L = pmap8L_path,
PMapPath_12L = pmap12L_path,
SectorMatchTool = FTKTrajectorSelector,
TrackBuilderTool = FTKTruthTrackBuilder)
theTrigFTKFastSimTruth.IBLMode = 2
theTrigFTKFastSimTruth.HitInputTool = FTKSGInput
theTrigFTKFastSimTruth.Clustering = True
theTrigFTKFastSimTruth.PixelClusteringMode = 101
topSequence += theTrigFTKFastSimTruth
from TrigFTK_RawDataAlgs.TrigFTK_RawDataAlgsConf import FTK_RDO_ReaderAlgo
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment