def TICL_iterations(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.985 # ~10 degrees ) process.MultiClustersFromTrackstersMIP = multiClustersFromTrackstersProducer.clone( label = "MIPMultiClustersFromTracksterByCA", Tracksters = "TrackstersMIP" ) process.FilteredLayerClusters = filteredLayerClustersProducer.clone( algo_number = 8, iteration_label = "algo8" ) 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.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.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(process): process.FEVTDEBUGHLTEventContent.outputCommands.extend(['keep *_MultiClustersFromTracksters*_*_*']) 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.985 # ~10 degrees ) process.MultiClustersFromTrackstersMIP = multiClustersFromTrackstersProducer.clone( label = "MIPMultiClustersFromTracksterByCA", Tracksters = "TrackstersMIP" ) process.FilteredLayerClusters = filteredLayerClustersProducer.clone( algo_number = 8, iteration_label = "algo8" ) 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.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.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
clusterFilter="ClusterFilterByAlgoAndSize", min_cluster_size=3, # inclusive algo_number=8, iteration_label="HAD", LayerClustersInputMask="ticlTrackstersTrk") # CA - PATTERN RECOGNITION ticlTrackstersHAD = _trackstersProducer.clone( filtered_mask="filteredLayerClustersHAD:HAD", original_mask='ticlTrackstersTrk', seeding_regions="ticlSeedingGlobal", # For the moment we mask everything w/o requirements since we are last # filter_on_categories = [5], # filter neutral hadrons # pid_threshold = 0.7, skip_layers=1, min_layers_per_trackster=12, min_cos_theta=0.866, # ~30 degrees min_cos_pointing=0.819, # ~35 degrees max_delta_time=-1, itername="HADRONIC") # MULTICLUSTERS ticlMultiClustersFromTrackstersHAD = _multiClustersFromTrackstersProducer.clone( Tracksters="ticlTrackstersHAD") ticlHADStepTask = cms.Task(ticlSeedingGlobal, filteredLayerClustersHAD, ticlTrackstersHAD, ticlMultiClustersFromTrackstersHAD)
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
import FWCore.ParameterSet.Config as cms from RecoHGCal.TICL.MIPStep_cff import * 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,
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
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", skip_layers = 3, min_layers_per_trackster = 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", skip_layers = 2, min_layers_per_trackster = 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
algo_number=8, LayerClustersInputMask='ticlTrackstersEM', iteration_label="Trk") # CA - PATTERN RECOGNITION ticlTrackstersTrk = _trackstersProducer.clone( filtered_mask="filteredLayerClustersTrk:Trk", seeding_regions="ticlSeedingTrk", original_mask='ticlTrackstersEM', filter_on_categories=[2, 4], # filter muons and charged hadrons pid_threshold=0.0, skip_layers=3, min_layers_per_trackster=10, min_cos_theta=0.866, # ~30 degrees min_cos_pointing=0.798, # ~ 37 degrees max_delta_time=-1., algo_verbosity=2, oneTracksterPerTrackSeed=True, promoteEmptyRegionToTrackster=True, itername="TRK") # MULTICLUSTERS ticlMultiClustersFromTrackstersTrk = _multiClustersFromTrackstersProducer.clone( Tracksters="ticlTrackstersTrk") ticlTrkStepTask = cms.Task(ticlSeedingTrk, filteredLayerClustersTrk, ticlTrackstersTrk, ticlMultiClustersFromTrackstersTrk)
# CA - PATTERN RECOGNITION ticlTrackstersMIP = _trackstersProducer.clone( filtered_mask="filteredLayerClustersMIP:MIP", seeding_regions="ticlSeedingGlobal", missing_layers=3, min_clusters_per_ntuplet=10, min_cos_theta=0.99, # ~10 degrees min_cos_pointing=0.5, out_in_dfs=False, itername="MIP", max_delta_time=-1) # MULTICLUSTERS ticlMultiClustersFromTrackstersMIP = _multiClustersFromTrackstersProducer.clone( Tracksters="ticlTrackstersMIP") ticlMIPStepTask = cms.Task(ticlSeedingGlobal, filteredLayerClustersMIP, ticlTrackstersMIP, ticlMultiClustersFromTrackstersMIP) filteredLayerClustersHFNoseMIP = filteredLayerClustersMIP.clone( LayerClusters='hgcalLayerClustersHFNose', LayerClustersInputMask="hgcalLayerClustersHFNose:InitialLayerClustersMask", iteration_label="MIPn", algo_number=9) ticlTrackstersHFNoseMIP = ticlTrackstersMIP.clone( detector="HFNose", layer_clusters="hgcalLayerClustersHFNose", layer_clusters_hfnose_tiles="ticlLayerTileHFNose",
import FWCore.ParameterSet.Config as cms from RecoHGCal.TICL.trackstersFromSimClustersProducer_cfi import trackstersFromSimClustersProducer as _trackstersFromSimClustersProducer from RecoHGCal.TICL.multiClustersFromTrackstersProducer_cfi import multiClustersFromTrackstersProducer as _multiClustersFromTrackstersProducer from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer # CA - PATTERN RECOGNITION filteredLayerClustersSimTracksters = _filteredLayerClustersProducer.clone( clusterFilter="ClusterFilterByAlgoAndSize", algo_number=8, min_cluster_size=0, # inclusive iteration_label="ticlSimTracksters") ticlSimTracksters = _trackstersFromSimClustersProducer.clone() from Configuration.ProcessModifiers.premix_stage2_cff import premix_stage2 premix_stage2.toModify(ticlSimTracksters, simclusters="mixData:MergedCaloTruth") ticlMultiClustersFromSimTracksters = _multiClustersFromTrackstersProducer.clone( Tracksters="ticlSimTracksters") ticlSimTrackstersTask = cms.Task(filteredLayerClustersSimTracksters, ticlSimTracksters, ticlMultiClustersFromSimTracksters)