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*_*_*',
'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.ticlSeedingRegionProducer_cfi import ticlSeedingRegionProducer as _ticlSeedingRegionProducer # SEEDING REGION ticlSeedingGlobal = _ticlSeedingRegionProducer.clone(algoId=2) ticlSeedingTrk = _ticlSeedingRegionProducer.clone(algoId=1) ticlSeedingGlobalHFNose = _ticlSeedingRegionProducer.clone(algoId=2)
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
import FWCore.ParameterSet.Config as cms
from RecoHGCal.TICL.ticlSeedingRegionProducer_cfi import ticlSeedingRegionProducer as _ticlSeedingRegionProducer
# SEEDING REGION
ticlSeedingGlobal = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionGlobal"))
ticlSeedingTrk = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionByTracks"))
ticlSeedingGlobalHFNose = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionGlobal"))
import FWCore.ParameterSet.Config as cms
from RecoHGCal.TICL.ticlSeedingRegionProducer_cfi import ticlSeedingRegionProducer as _ticlSeedingRegionProducer
# SEEDING REGION
ticlSeedingGlobal = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionGlobal"))
ticlSeedingTrk = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionByTracks"))
ticlSeedingGlobalHFNose = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionGlobal"))
# fixme: empty seed, need to work on propagator
ticlSeedingTrkHFNose = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionByTracks",
cutTk=cms.string(
'3. < abs(eta) < 4. && pt > 1. && quality("highPurity") && hitPattern().numberOfLostHits("MISSING_OUTER_HITS") < 5'
),
detector=cms.string("HFNose")))
ticlSeedingByHFHFNose = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionByHF"))
import FWCore.ParameterSet.Config as cms
from RecoHGCal.TICL.ticlSeedingRegionProducer_cfi import ticlSeedingRegionProducer as _ticlSeedingRegionProducer
# SEEDING REGION
ticlSeedingGlobal = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionGlobal"))
ticlSeedingTrk = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionByTracks"))
ticlSeedingGlobalHFNose = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionGlobal"))
# fixme: empty seed, need to work on propagator
ticlSeedingTrkHFNose = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(
type="SeedingRegionByTracks",
cutTk=cms.string(
'3. < abs(eta) < 4. && pt > 2. &&' +
'quality("highPurity") && numberOfValidHits() > 12 && normalizedChi2() > 0.7 &&'
+ 'hitPattern().numberOfLostHits("MISSING_OUTER_HITS") < 5'),
detector=cms.string("HFNose"),
propagator=cms.string("RungeKuttaTrackerPropagator")))
ticlSeedingByHFHFNose = _ticlSeedingRegionProducer.clone(
seedingPSet=_ticlSeedingRegionProducer.seedingPSet.clone(