from AthenaCommon.AlgSequence import AlgSequence
theJob = AlgSequence()
from FastTrackSimWrap.FastTrackSimWrapConf import FTKRegionalWrapper
if hasattr(runArgs, "outputNTUP_FTKIPFile"):
OutputNTUP_FTKIPFile = runArgs.outputNTUP_FTKIPFile
else:
OutputNTUP_FTKIPFile = "ftksim_32Towers_wrap.root"
from AthenaCommon.AppMgr import ToolSvc
from TrigFTKSim.TrigFTKSimConf import FTK_SGHitInput
FTKSGInput = FTK_SGHitInput(maxEta=3.2, minPt=0.8 * GeV)
FTKSGInput.OutputLevel = DEBUG
FTKSGInput.ReadTruthTracks = True
ToolSvc += FTKSGInput
print "Output file", OutputNTUP_FTKIPFile
wrapper = FTKRegionalWrapper(OutputLevel=DEBUG,
PMapPath=pmap_path,
RMapPath=rmap_path,
OutFileName=OutputNTUP_FTKIPFile)
wrapper.IBLMode = 2
wrapper.FixEndcapL0 = False
wrapper.HitInputTool = FTKSGInput
theJob += wrapper
print theJob
###############################################################
# FTKRegionalWrapper job options file
#
#==============================================================
# Helper function from transforms
from PyJobTransforms.trfUtils import findFile
pmap_path = findFile(os.environ['DATAPATH'], 'ftk_configuration/map_files/raw_11L.pmap')
print "Using PMAP:", pmap_path
rmap_path = findFile(os.environ['DATAPATH'], 'ftk_configuration/map_files/raw_11L.tmap')
print "Using RMAP:", rmap_path
from AthenaCommon.AlgSequence import AlgSequence
theJob = AlgSequence()
from FastTrackSimWrap.FastTrackSimWrapConf import FTKRegionalWrapper
if hasattr(runArgs,"outputNTUP_FTKIPFile") :
OutputNTUP_FTKIPFile = runArgs.outputNTUP_FTKIPFile
else :
OutputNTUP_FTKIPFile = 'ftksim_newwrap_raw.root'
print "Output file", OutputNTUP_FTKIPFile
theJob += FTKRegionalWrapper(OutputLevel = DEBUG,
PMapPath = pmap_path,
RMapPath = rmap_path,
OutFileName = "ftksim_regionalwrap.root")
print theJob
###############################################################
pmap_path = findFile(os.environ['DATAPATH'],
'ftk_configuration/map_files/raw_11L.pmap')
print "Using PMAP:", pmap_path
rmap_path = findFile(os.environ['DATAPATH'],
'ftk_configuration/map_files/raw_11L.rmap')
print "Using RMAP:", rmap_path
from AthenaCommon.AlgSequence import AlgSequence
theJob = AlgSequence()
from FastTrackSimWrap.FastTrackSimWrapConf import FTKRegionalWrapper
if hasattr(runArgs, "outputNTUP_FTKIPFile"):
OutputNTUP_FTKIPFile = runArgs.outputNTUP_FTKIPFile
else:
OutputNTUP_FTKIPFile = "ftksim_8Reg_wrap.root"
from TrigFTKSim.TrigFTKSimConf import FTK_SGHitInput
FTKSGInput = FTK_SGHitInput()
FTKSGInput.ReadTruthTracks = True
ToolSvc += FTKSGInput
print "Output file", OutputNTUP_FTKIPFile
theJob += FTKRegionalWrapper(OutputLevel=DEBUG,
PMapPath=pmap_path,
RMapPath=rmap_path,
OutFileName=OutputNTUP_FTKIPFile)
print theJob
###############################################################
FTKSGInput = FTK_SGHitInput(maxEta=3.2, minPt=0.8 * GeV)
FTKSGInput.OutputLevel = DEBUG
FTKSGInput.ReadTruthTracks = True
ToolSvc += FTKSGInput
from TrigFTKTrackConverter.TrigFTKTrackConverterConf import TrigFTKClusterConverterTool
ClusterConvTool = TrigFTKClusterConverterTool(UsePixelCalibSvc=False)
ToolSvc += ClusterConvTool
print "Output file", OutputNTUP_FTKIPFile
wrapper = FTKRegionalWrapper(OutputLevel=DEBUG,
PMapPath=pmap_path,
RMapPath=rmap_path,
OutFileName=OutputNTUP_FTKIPFile,
WriteClustersToESD=True,
SaveHits=True,
Clustering=True,
TrigFTKClusterConverterTool=ClusterConvTool)
wrapper.IBLMode = 1
wrapper.HitInputTool = FTKSGInput
wrapper.Clustering = True
theJob += wrapper
from RecExConfig.ObjKeyStore import objKeyStore
objKeyStore.addStreamESD("InDet::SCT_ClusterContainer", "FTK_SCT_Cluster")
objKeyStore.addStreamESD("InDet::PixelClusterContainer", "FTK_Pixel_Clusters")
objKeyStore.addStreamESD("Trk::PRD_MultiTruthCollection",
"PRD_MultiTruthSCT_FTK")
objKeyStore.addStreamESD("PRD_MultiTruthCollection", "PRD_MultiTruthSCT_FTK")
objKeyStore.addStreamESD("PRD_MultiTruthCollection", "PRD_MultiTruthPixel_FTK")
theJob = AlgSequence()
from FastTrackSimWrap.FastTrackSimWrapConf import FTKRegionalWrapper
if hasattr(runArgs,"outputNTUP_FTKIPFile") :
OutputNTUP_FTKIPFile = runArgs.outputNTUP_FTKIPFile
else :
OutputNTUP_FTKIPFile = "ftksim_ITk_wrap.root"
from AthenaCommon.AppMgr import ToolSvc
from TrigFTKSim.TrigFTKSimConf import FTK_SGHitInput
FTKSGInput = FTK_SGHitInput(maxEta=3.0, minPt=4*GeV)
FTKSGInput.OutputLevel = DEBUG
FTKSGInput.ReadTruthTracks = True
FTKSGInput.DoOutFileRawHits = True
ToolSvc += FTKSGInput
print "Output file", OutputNTUP_FTKIPFile
wrapper = FTKRegionalWrapper(OutputLevel = DEBUG,
PMapPath = pmap_path,
RMapPath = rmap_path,
OutFileName = OutputNTUP_FTKIPFile)
wrapper.IBLMode = 0
wrapper.ITkMode = True
wrapper.HitInputTool = FTKSGInput
theJob += wrapper
print theJob
###############################################################
OutputNTUP_FTKIPFile = runArgs.outputNTUP_FTKIPFile
else:
OutputNTUP_FTKIPFile = "ftksim_32Towers_wrap.root"
from AthenaCommon.AppMgr import ToolSvc
from TrigFTKSim.TrigFTKSimConf import FTK_SGHitInput
FTKSGInput = FTK_SGHitInput(maxEta=3.2, minPt=0.8 * GeV)
FTKSGInput.OutputLevel = DEBUG
FTKSGInput.ReadTruthTracks = False
ToolSvc += FTKSGInput
print "Output file", OutputNTUP_FTKIPFile
wrapper = FTKRegionalWrapper(OutputLevel=DEBUG,
PMapPath=pmap_path,
RMapPath=rmap_path,
OutFileName=OutputNTUP_FTKIPFile)
wrapper.IBLMode = 2
wrapper.HitInputTool = FTKSGInput
wrapper.Clustering = True
wrapper.PixelClusteringMode = 101
wrapper.Ibl3DRealistic = False
wrapper.SavePerPlane = True
wrapper.EmulateDF = True
wrapper.SaveRawHits = True
wrapper.SaveHits = True
wrapper.DuplicateGanged = False
wrapper.SctClustering = True
print "FTKRegionalWrapper setting ROBS"
wrapper.pixRodIds = [
OutputNTUP_FTKIPFile = runArgs.outputNTUP_FTKIPFile
else:
OutputNTUP_FTKIPFile = "ftksim_64Towers_wrap.root"
from AthenaCommon.AppMgr import ToolSvc
from TrigFTKSim.TrigFTKSimConf import FTK_SGHitInput
FTKSGInput = FTK_SGHitInput(maxEta=3.2, minPt=0.8 * GeV)
FTKSGInput.OutputLevel = DEBUG
FTKSGInput.ReadTruthTracks = False
ToolSvc += FTKSGInput
print "Output file", OutputNTUP_FTKIPFile
wrapper = FTKRegionalWrapper(OutputLevel=INFO,
PMapPath=pmap_path,
RMapPath=rmap_path,
OutFileName=OutputNTUP_FTKIPFile)
wrapper.IBLMode = 1
wrapper.HitInputTool = FTKSGInput
wrapper.Clustering = True
wrapper.PixelClusteringMode = 100
wrapper.Ibl3DRealistic = False
wrapper.SavePerPlane = True
wrapper.EmulateDF = True
wrapper.SaveRawHits = True
wrapper.SaveHits = True
#wrapper.pixRodIds = {0x130011, 0x111510, 0x111508, 0x112414, 0x130015, 0x111716, 0x112416}
#wrapper.sctRodIds = {0x220005, 0x210005, 0x220007}
#wrapper.pixRodIds = {130011, 111510, 111508, 112414, 130015, 111716, 112416}
OutputNTUP_FTKIPFile = runArgs.outputNTUP_FTKIPFile
else:
OutputNTUP_FTKIPFile = "ftksim_64Towers_wrap.root"
from AthenaCommon.AppMgr import ToolSvc
from TrigFTKSim.TrigFTKSimConf import FTK_SGHitInput
FTKSGInput = FTK_SGHitInput(maxEta=3.2, minPt=0.8 * GeV)
FTKSGInput.OutputLevel = DEBUG
FTKSGInput.ReadTruthTracks = True
ToolSvc += FTKSGInput
print "Output file", OutputNTUP_FTKIPFile
wrapper = FTKRegionalWrapper(OutputLevel=DEBUG,
PMapPath=pmap_path,
RMapPath=rmap_path,
OutFileName=OutputNTUP_FTKIPFile)
wrapper.IBLMode = 2
wrapper.FixEndcapL0 = True
wrapper.HitInputTool = FTKSGInput
wrapper.PixelClusteringMode = 101
wrapper.Ibl3DRealistic = False
wrapper.SctClustering = True
wrapper.Clustering = True
####wrapper.GetOffline = True
wrapper.GetOffline = False
theJob += wrapper
print theJob