Beispiel #1
0
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
Beispiel #2
0
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
Beispiel #3
0
def TICL_iterations(process):
    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.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.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.filteredLayerClusters,
        process.tracksters, process.hgcalMultiClusters)
    process.schedule = cms.Schedule(process.raw2digi_step,
                                    process.FEVTDEBUGHLToutput_step)
    process.schedule.associate(process.TICL_Task)
    return process
Beispiel #4
0
import FWCore.ParameterSet.Config as cms

from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingGlobal
from RecoHGCal.TICL.ticlLayerTileProducer_cfi import ticlLayerTileProducer as _ticlLayerTileProducer
from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer
from RecoHGCal.TICL.multiClustersFromTrackstersProducer_cfi import multiClustersFromTrackstersProducer as _multiClustersFromTrackstersProducer

# CLUSTER FILTERING/MASKING

filteredLayerClustersHAD = _filteredLayerClustersProducer.clone(
    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,
Beispiel #5
0
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
Beispiel #6
0
import FWCore.ParameterSet.Config as cms

from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingGlobal, ticlSeedingGlobalHFNose
from RecoHGCal.TICL.ticlLayerTileProducer_cfi import ticlLayerTileProducer as _ticlLayerTileProducer
from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer

# CLUSTER FILTERING/MASKING

filteredLayerClustersHAD = _filteredLayerClustersProducer.clone(
    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,
    pluginPatternRecognitionByCA=dict(
        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),
Beispiel #7
0
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
Beispiel #8
0
import FWCore.ParameterSet.Config as cms

from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingGlobal, ticlSeedingGlobalHFNose
from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer
from RecoHGCal.TICL.multiClustersFromTrackstersProducer_cfi import multiClustersFromTrackstersProducer as _multiClustersFromTrackstersProducer

# CLUSTER FILTERING/MASKING

filteredLayerClustersEM = _filteredLayerClustersProducer.clone(
    clusterFilter="ClusterFilterByAlgoAndSize",
    min_cluster_size=2,  # inclusive
    algo_number=8,
    LayerClustersInputMask='ticlTrackstersTrk',
    iteration_label="EM")

# CA - PATTERN RECOGNITION

ticlTrackstersEM = _trackstersProducer.clone(
    filtered_mask="filteredLayerClustersEM:EM",
    original_mask='ticlTrackstersTrk',
    seeding_regions="ticlSeedingGlobal",
    filter_on_categories=[0, 1],
    pid_threshold=0.8,
    max_out_in_hops=4,
    missing_layers=1,
    min_clusters_per_ntuplet=10,
    min_cos_theta=0.978,  # ~12 degrees
    min_cos_pointing=0.9,  # ~25 degrees
    max_delta_time=3.,
    itername="EM",
Beispiel #9
0
import FWCore.ParameterSet.Config as cms

from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingGlobal, ticlSeedingGlobalHFNose
from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer
from RecoHGCal.TICL.multiClustersFromTrackstersProducer_cfi import multiClustersFromTrackstersProducer as _multiClustersFromTrackstersProducer

# CLUSTER FILTERING/MASKING

filteredLayerClustersCLUE3DHigh = _filteredLayerClustersProducer.clone(
    clusterFilter = "ClusterFilterByAlgoAndSize",
    min_cluster_size = 2, # inclusive
    algo_number = 8,
    iteration_label = "CLUE3DHigh"
)

# PATTERN RECOGNITION

ticlTrackstersCLUE3DHigh = _trackstersProducer.clone(
    filtered_mask = "filteredLayerClustersCLUE3DHigh:CLUE3DHigh",
    seeding_regions = "ticlSeedingGlobal",
    itername = "CLUE3DHigh",
    patternRecognitionBy = "CLUE3D",
    pluginPatternRecognitionByCLUE3D = dict (
        criticalEtaPhiDistance = 0.025
    )

)

ticlCLUE3DHighStepTask = cms.Task(ticlSeedingGlobal
    ,filteredLayerClustersCLUE3DHigh
Beispiel #10
0
import FWCore.ParameterSet.Config as cms

from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingTrk, ticlSeedingTrkHFNose
from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer

# CLUSTER FILTERING/MASKING

filteredLayerClustersTrkEM = _filteredLayerClustersProducer.clone(
    clusterFilter = "ClusterFilterByAlgoAndSizeAndLayerRange",
    min_cluster_size = 3, # inclusive
    max_layerId = 30, # inclusive
    algo_number = 8,
    iteration_label = "TrkEM"
)

# CA - PATTERN RECOGNITION

ticlTrackstersTrkEM = _trackstersProducer.clone(
    filtered_mask = "filteredLayerClustersTrkEM:TrkEM",
    seeding_regions = "ticlSeedingTrk",
    pluginPatternRecognitionByCA = dict(
        algo_verbosity = 0,
        filter_on_categories = [0, 1],
        pid_threshold = 0.5,
        energy_em_over_total_threshold = 0.9,
        max_longitudinal_sigmaPCA = 10,
        shower_start_max_layer = 5, #inclusive
        max_out_in_hops = 1,
        max_missing_layers_in_trackster = 2,
        skip_layers = 2,
Beispiel #11
0
import FWCore.ParameterSet.Config as cms

from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingGlobal, ticlSeedingGlobalHFNose
from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer

# CLUSTER FILTERING/MASKING

filteredLayerClustersMIP = _filteredLayerClustersProducer.clone(
    clusterFilter="ClusterFilterBySize",
    algo_number=8,
    max_cluster_size=2,  # inclusive
    iteration_label="MIP")

# CA - PATTERN RECOGNITION

ticlTrackstersMIP = _trackstersProducer.clone(
    filtered_mask="filteredLayerClustersMIP:MIP",
    seeding_regions="ticlSeedingGlobal",
    pluginPatternRecognitionByCA=dict(
        skip_layers=3,
        min_layers_per_trackster=10,
        min_cos_theta=0.99,  # ~10 degrees
        min_cos_pointing=0.5,
        out_in_dfs=False,
        max_delta_time=-1),
    itername="MIP")

ticlMIPStepTask = cms.Task(ticlSeedingGlobal, filteredLayerClustersMIP,
                           ticlTrackstersMIP)
Beispiel #12
0
import FWCore.ParameterSet.Config as cms

from RecoHGCal.TICL.trackstersFromSimClustersProducer_cfi import trackstersFromSimClustersProducer as _trackstersFromSimClustersProducer
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",
    caloparticles="mixData:MergedCaloTruth",
)

ticlSimTrackstersTask = cms.Task(filteredLayerClustersSimTracksters,
                                 ticlSimTracksters)
Beispiel #13
0
import FWCore.ParameterSet.Config as cms

from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingTrk
from RecoHGCal.TICL.ticlLayerTileProducer_cfi import ticlLayerTileProducer as _ticlLayerTileProducer
from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer
from RecoHGCal.TICL.multiClustersFromTrackstersProducer_cfi import multiClustersFromTrackstersProducer as _multiClustersFromTrackstersProducer

# CLUSTER FILTERING/MASKING

filteredLayerClustersTrk = _filteredLayerClustersProducer.clone(
    clusterFilter="ClusterFilterByAlgo", algo_number=8, iteration_label="Trk")

# CA - PATTERN RECOGNITION

ticlTrackstersTrk = _trackstersProducer.clone(
    filtered_mask="filteredLayerClustersTrk:Trk",
    seeding_regions="ticlSeedingTrk",
    filter_on_categories=[2, 4],  # filter muons and charged hadrons
    pid_threshold=0.0,
    missing_layers=3,
    min_clusters_per_ntuplet=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
Beispiel #14
0
import FWCore.ParameterSet.Config as cms

from RecoHGCal.TICL.TICLSeedingRegions_cff import ticlSeedingGlobal, ticlSeedingGlobalHFNose
from RecoHGCal.TICL.trackstersProducer_cfi import trackstersProducer as _trackstersProducer
from RecoHGCal.TICL.filteredLayerClustersProducer_cfi import filteredLayerClustersProducer as _filteredLayerClustersProducer
from RecoHGCal.TICL.multiClustersFromTrackstersProducer_cfi import multiClustersFromTrackstersProducer as _multiClustersFromTrackstersProducer

# CLUSTER FILTERING/MASKING

filteredLayerClustersFastJet = _filteredLayerClustersProducer.clone(
    clusterFilter="ClusterFilterByAlgoAndSize",
    min_cluster_size=3,  # inclusive
    algo_number=8,
    iteration_label="FastJet")

# PATTERN RECOGNITION

ticlTrackstersFastJet = _trackstersProducer.clone(
    filtered_mask="filteredLayerClustersFastJet:FastJet",
    seeding_regions="ticlSeedingGlobal",
    itername="FastJet",
    patternRecognitionBy="FastJet",
    pluginPatternRecognitionByFastJet=dict(algo_verbosity=2))

ticlFastJetStepTask = cms.Task(ticlSeedingGlobal, filteredLayerClustersFastJet,
                               ticlTrackstersFastJet)