def TICL_iterations(process):
process.FEVTDEBUGHLTEventContent.outputCommands.extend(
['keep *_multiClustersFromTracksters*_*_*'])
process.ticlLayerTileProducer = ticlLayerTileProducer.clone()
process.ticlSeedingGlobal = ticlSeedingRegionProducer.clone(algoId=2)
process.filteredLayerClustersMIP = filteredLayerClustersProducer.clone(
clusterFilter="ClusterFilterBySize",
algo_number=8,
max_cluster_size=2, # inclusive
iteration_label="MIP")
process.trackstersMIP = trackstersProducer.clone(
filtered_mask="filteredLayerClustersMIP:MIP",
seeding_regions="ticlSeedingGlobal",
missing_layers=3,
min_clusters_per_ntuplet=15,
min_cos_theta=0.99, # ~10 degrees
)
process.multiClustersFromTrackstersMIP = multiClustersFromTrackstersProducer.clone(
label="MIPMultiClustersFromTracksterByCA", Tracksters="trackstersMIP")
process.filteredLayerClusters = filteredLayerClustersProducer.clone(
clusterFilter="ClusterFilterByAlgoAndSize",
min_cluster_size=2,
algo_number=8,
iteration_label="algo8")
process.tracksters = trackstersProducer.clone(
original_mask="trackstersMIP",
filtered_mask="filteredLayerClusters:algo8",
seeding_regions="ticlSeedingGlobal",
missing_layers=2,
min_clusters_per_ntuplet=15,
min_cos_theta=0.94, # ~20 degrees
min_cos_pointing=0.7)
process.multiClustersFromTracksters = multiClustersFromTrackstersProducer.clone(
Tracksters="tracksters")
process.HGCalUncalibRecHit = HGCalUncalibRecHit
process.HGCalRecHit = HGCalRecHit
process.hgcalLayerClusters = hgcalLayerClusters
process.hgcalMultiClusters = hgcalMultiClusters
process.TICL_Task = cms.Task(
process.HGCalUncalibRecHit, process.HGCalRecHit,
process.hgcalLayerClusters, process.filteredLayerClustersMIP,
process.ticlLayerTileProducer, process.ticlSeedingGlobal,
process.trackstersMIP, process.multiClustersFromTrackstersMIP,
process.filteredLayerClusters, process.tracksters,
process.multiClustersFromTracksters, process.hgcalMultiClusters)
process.schedule = cms.Schedule(process.raw2digi_step,
process.FEVTDEBUGHLToutput_step)
process.schedule.associate(process.TICL_Task)
return process
def TICL_iterations_withReco(process):
process.FEVTDEBUGHLTEventContent.outputCommands.extend(
['keep *_MultiClustersFromTracksters*_*_*'])
process.TICLLayerTileProducer = ticlLayerTileProducer.clone()
process.FilteredLayerClustersMIP = filteredLayerClustersProducer.clone(
clusterFilter="ClusterFilterBySize",
algo_number=8,
max_cluster_size=2, # inclusive
iteration_label="MIP")
process.TrackstersMIP = trackstersProducer.clone(
filtered_mask=cms.InputTag("FilteredLayerClustersMIP", "MIP"),
missing_layers=3,
min_clusters_per_ntuplet=15,
min_cos_theta=0.99, # ~10 degrees
min_cos_pointing=0.9)
process.MultiClustersFromTrackstersMIP = multiClustersFromTrackstersProducer.clone(
label="MIPMultiClustersFromTracksterByCA", Tracksters="TrackstersMIP")
process.FilteredLayerClusters = filteredLayerClustersProducer.clone(
clusterFilter="ClusterFilterByAlgoAndSize",
min_cluster_size=2,
algo_number=8,
iteration_label="algo8",
LayerClustersInputMask="TrackstersMIP")
process.Tracksters = trackstersProducer.clone(
original_mask="TrackstersMIP",
filtered_mask=cms.InputTag("FilteredLayerClusters", "algo8"),
missing_layers=2,
min_clusters_per_ntuplet=15,
min_cos_theta=0.94, # ~20 degrees
min_cos_pointing=0.7)
process.MultiClustersFromTracksters = multiClustersFromTrackstersProducer.clone(
Tracksters="Tracksters")
process.hgcalMultiClusters = hgcalMultiClusters
process.TICL_Task = cms.Task(
process.TICLLayerTileProducer, process.FilteredLayerClustersMIP,
process.TrackstersMIP, process.MultiClustersFromTrackstersMIP,
process.FilteredLayerClusters, process.Tracksters,
process.MultiClustersFromTracksters)
process.schedule.associate(process.TICL_Task)
return process
def TICL_iterations_withReco(process):
process.FEVTDEBUGHLTEventContent.outputCommands.extend([
'keep *_multiClustersFromTracksters*_*_*',
'keep *_ticlCandidateFromTrackstersProducer*_*_*',
'keep *_pfTICLProducer*_*_*'
])
process.ticlLayerTileProducer = ticlLayerTileProducer.clone()
process.ticlSeedingTrk = ticlSeedingRegionProducer.clone(algoId=1)
process.filteredLayerClustersTrk = filteredLayerClustersProducer.clone(
clusterFilter="ClusterFilterByAlgo",
algo_number=8,
iteration_label="Trk")
process.trackstersTrk = trackstersProducer.clone(
filtered_mask="filteredLayerClustersTrk:Trk",
seeding_regions="ticlSeedingTrk",
missing_layers=3,
min_clusters_per_ntuplet=5,
min_cos_theta=
0.99, # ~10 degrees
min_cos_pointing=0.9)
process.multiClustersFromTrackstersTrk = multiClustersFromTrackstersProducer.clone(
label="TrkMultiClustersFromTracksterByCA", Tracksters="trackstersTrk")
process.ticlSeedingGlobal = ticlSeedingRegionProducer.clone(algoId=2)
process.filteredLayerClustersMIP = filteredLayerClustersProducer.clone(
clusterFilter="ClusterFilterBySize",
algo_number=8,
max_cluster_size=2, # inclusive
iteration_label="MIP")
process.trackstersMIP = trackstersProducer.clone(
filtered_mask="filteredLayerClustersMIP:MIP",
seeding_regions="ticlSeedingGlobal",
missing_layers=3,
min_clusters_per_ntuplet=15,
min_cos_theta=0.99, # ~10 degrees
min_cos_pointing=0.9,
out_in_dfs=False,
)
process.multiClustersFromTrackstersMIP = multiClustersFromTrackstersProducer.clone(
label="MIPMultiClustersFromTracksterByCA", Tracksters="trackstersMIP")
process.filteredLayerClusters = filteredLayerClustersProducer.clone(
clusterFilter="ClusterFilterByAlgoAndSize",
min_cluster_size=2,
algo_number=8,
iteration_label="algo8",
LayerClustersInputMask="trackstersMIP")
process.trackstersEM = trackstersProducer.clone(
max_out_in_hops=4,
original_mask="trackstersMIP",
filtered_mask="filteredLayerClusters:algo8",
seeding_regions="ticlSeedingGlobal",
missing_layers=1,
min_clusters_per_ntuplet=10,
min_cos_theta=0.984, # ~10 degrees
min_cos_pointing=0.9 # ~26 degrees
)
process.multiClustersFromTrackstersEM = multiClustersFromTrackstersProducer.clone(
Tracksters="trackstersEM")
process.trackstersHAD = trackstersProducer.clone(
filtered_mask="filteredLayerClusters:algo8",
seeding_regions="ticlSeedingGlobal",
missing_layers=2,
min_clusters_per_ntuplet=10,
min_cos_theta=0.8,
min_cos_pointing=0.7)
process.multiClustersFromTrackstersHAD = multiClustersFromTrackstersProducer.clone(
Tracksters="trackstersHAD")
process.ticlCandidateFromTrackstersProducer = ticlCandidateFromTrackstersProducer.clone(
)
process.pfTICLProducer = pfTICLProducer.clone()
process.hgcalMultiClusters = hgcalMultiClusters
process.TICL_Task = cms.Task(
process.ticlLayerTileProducer, process.ticlSeedingTrk,
process.filteredLayerClustersTrk, process.trackstersTrk,
process.multiClustersFromTrackstersTrk, process.ticlSeedingGlobal,
process.filteredLayerClustersMIP, process.trackstersMIP,
process.multiClustersFromTrackstersMIP, process.filteredLayerClusters,
process.trackstersEM, process.multiClustersFromTrackstersEM,
process.trackstersHAD, process.multiClustersFromTrackstersHAD,
process.ticlCandidateFromTrackstersProducer, process.pfTICLProducer)
process.schedule.associate(process.TICL_Task)
process.ticlPFValidation = ticlPFValidation
process.hgcalValidation.insert(-1, process.ticlPFValidation)
if getattr(process, 'hgcalValidator'):
process.hgcalValidator.label_lcl = "hgcalLayerClusters"
process.hgcalValidator.label_mcl = [
"multiClustersFromTrackstersEM:MultiClustersFromTracksterByCA",
"multiClustersFromTrackstersHAD:MultiClustersFromTracksterByCA"
]
process.hgcalValidator.domulticlustersPlots = True
return process
from RecoHGCal.TICL.TrkEMStep_cff import *
from RecoHGCal.TICL.TrkStep_cff import *
from RecoHGCal.TICL.EMStep_cff import *
from RecoHGCal.TICL.HADStep_cff import *
from RecoHGCal.TICL.ticlLayerTileProducer_cfi import ticlLayerTileProducer
from RecoHGCal.TICL.pfTICLProducer_cfi import pfTICLProducer as _pfTICLProducer
from RecoHGCal.TICL.trackstersMergeProducer_cfi import trackstersMergeProducer as _trackstersMergeProducer
from RecoHGCal.TICL.multiClustersFromTrackstersProducer_cfi import multiClustersFromTrackstersProducer as _multiClustersFromTrackstersProducer
ticlLayerTileTask = cms.Task(ticlLayerTileProducer)
ticlTrackstersMerge = _trackstersMergeProducer.clone()
ticlMultiClustersFromTrackstersMerge = _multiClustersFromTrackstersProducer.clone(
Tracksters="ticlTrackstersMerge")
ticlTracksterMergeTask = cms.Task(ticlTrackstersMerge,
ticlMultiClustersFromTrackstersMerge)
pfTICL = _pfTICLProducer.clone()
ticlPFTask = cms.Task(pfTICL)
iterTICLTask = cms.Task(ticlLayerTileTask, ticlTrkEMStepTask, ticlEMStepTask,
ticlTrkStepTask, ticlHADStepTask,
ticlTracksterMergeTask, ticlPFTask)
ticlLayerTileHFNose = ticlLayerTileProducer.clone(detector='HFNose')
ticlLayerTileHFNoseTask = cms.Task(ticlLayerTileHFNose)
iterHFNoseTICLTask = cms.Task(ticlLayerTileHFNoseTask, ticlHFNoseMIPStepTask,
ticlHFNoseEMStepTask)
def TICL_iterations_withReco(process):
process.FEVTDEBUGHLTEventContent.outputCommands.extend(['keep *_multiClustersFromTracksters*_*_*'])
process.ticlLayerTileProducer = ticlLayerTileProducer.clone()
process.ticlSeedingTrk = ticlSeedingRegionProducer.clone(
algoId = 1
)
process.filteredLayerClustersTrk = filteredLayerClustersProducer.clone(
clusterFilter = "ClusterFilterByAlgo",
algo_number = 8,
iteration_label = "Trk"
)
process.trackstersTrk = trackstersProducer.clone(
filtered_mask = cms.InputTag("filteredLayerClustersTrk", "Trk"),
seeding_regions = "ticlSeedingTrk",
missing_layers = 3,
min_clusters_per_ntuplet = 5,
min_cos_theta = 0.99, # ~10 degrees
min_cos_pointing = 0.9
)
process.multiClustersFromTrackstersTrk = multiClustersFromTrackstersProducer.clone(
label = "TrkMultiClustersFromTracksterByCA",
Tracksters = "trackstersTrk"
)
process.ticlSeedingGlobal = ticlSeedingRegionProducer.clone(
algoId = 2
)
process.filteredLayerClustersMIP = filteredLayerClustersProducer.clone(
clusterFilter = "ClusterFilterBySize",
algo_number = 8,
max_cluster_size = 2, # inclusive
iteration_label = "MIP"
)
process.trackstersMIP = trackstersProducer.clone(
filtered_mask = cms.InputTag("filteredLayerClustersMIP", "MIP"),
seeding_regions = "ticlSeedingGlobal",
missing_layers = 3,
min_clusters_per_ntuplet = 15,
min_cos_theta = 0.99, # ~10 degrees
min_cos_pointing = 0.9,
out_in_dfs = False,
)
process.multiClustersFromTrackstersMIP = multiClustersFromTrackstersProducer.clone(
label = "MIPMultiClustersFromTracksterByCA",
Tracksters = "trackstersMIP"
)
process.filteredLayerClusters = filteredLayerClustersProducer.clone(
clusterFilter = "ClusterFilterByAlgoAndSize",
min_cluster_size = 2,
algo_number = 8,
iteration_label = "algo8",
LayerClustersInputMask = "trackstersMIP"
)
process.trackstersEM = trackstersProducer.clone(
original_mask = "trackstersMIP",
filtered_mask = cms.InputTag("filteredLayerClusters", "algo8"),
seeding_regions = "ticlSeedingGlobal",
missing_layers = 2,
min_clusters_per_ntuplet = 10,
min_cos_theta = 0.94, # ~20 degrees
min_cos_pointing = 0.7
)
process.multiClustersFromTrackstersEM = multiClustersFromTrackstersProducer.clone(
Tracksters = "trackstersEM"
)
process.trackstersHAD = trackstersProducer.clone(
filtered_mask = cms.InputTag("filteredLayerClusters", "algo8"),
seeding_regions = "ticlSeedingGlobal",
missing_layers = 2,
min_clusters_per_ntuplet = 10,
min_cos_theta = 0.8,
min_cos_pointing = 0.7
)
process.multiClustersFromTrackstersHAD = multiClustersFromTrackstersProducer.clone(
Tracksters = "trackstersHAD"
)
process.hgcalMultiClusters = hgcalMultiClusters
process.TICL_Task = cms.Task(
process.ticlLayerTileProducer,
process.ticlSeedingTrk,
process.filteredLayerClustersTrk,
process.trackstersTrk,
process.multiClustersFromTrackstersTrk,
process.ticlSeedingGlobal,
process.filteredLayerClustersMIP,
process.trackstersMIP,
process.multiClustersFromTrackstersMIP,
process.filteredLayerClusters,
process.trackstersEM,
process.multiClustersFromTrackstersEM,
process.trackstersHAD,
process.multiClustersFromTrackstersHAD)
process.schedule.associate(process.TICL_Task)
return process
