Skip to content

Commit

Permalink
Merge pull request #41589 from OlivieFranklova/split_layerClusterProd…
Browse files Browse the repository at this point in the history
…ucer

Split HGCal layer cluster producer
  • Loading branch information
cmsbuild authored Jun 1, 2023
2 parents 282ef9e + 651e4d6 commit 2c948f2
Show file tree
Hide file tree
Showing 96 changed files with 1,365 additions and 951 deletions.
1 change: 0 additions & 1 deletion DataFormats/CaloRecHit/interface/CaloCluster.h
Original file line number Diff line number Diff line change
Expand Up @@ -39,7 +39,6 @@ namespace reco {
particleFlow = 5,
hgcal_em = 6,
hgcal_had = 7,
hgcal_mixed = 8,
hfnose = 9,
undefined = 1000
};
Expand Down
2 changes: 1 addition & 1 deletion DataFormats/HGCalReco/interface/Common.h
Original file line number Diff line number Diff line change
Expand Up @@ -29,7 +29,7 @@ namespace ticl {
} // namespace ticl

namespace ticl {
typedef std::vector<std::pair<unsigned int, float> > HgcalClusterFilterMask;
typedef std::vector<std::pair<unsigned int, float> > TICLClusterFilterMask;
} // namespace ticl

#endif // DataFormats_HGCalReco_Common_h
3 changes: 2 additions & 1 deletion Fireworks/Calo/plugins/FWCaloClusterProxyBuilder.cc
Original file line number Diff line number Diff line change
Expand Up @@ -106,7 +106,8 @@ void FWCaloClusterProxyBuilder::build(const reco::CaloCluster &iData,
continue;

// HGCal
if (iData.algo() == 8 || (type >= 8 && type <= 10)) {
if (iData.algo() == reco::CaloCluster::hgcal_em || iData.algo() == reco::CaloCluster::hgcal_had ||
(type >= 8 && type <= 10)) {
if (heatmap && hitmap->find(it->first) == hitmap->end())
continue;

Expand Down
Original file line number Diff line number Diff line change
@@ -1,9 +1,8 @@
import FWCore.ParameterSet.Config as cms

filteredLayerClustersCLUE3DHighL1Seeded = cms.EDProducer("FilteredLayerClustersProducer",
LayerClusters = cms.InputTag("hgcalLayerClustersL1Seeded"),
LayerClustersInputMask = cms.InputTag("hgcalLayerClustersL1Seeded","InitialLayerClustersMask"),
algo_number = cms.int32(8),
LayerClusters = cms.InputTag("hgcalMergeLayerClustersL1Seeded"),
LayerClustersInputMask = cms.InputTag("hgcalMergeLayerClustersL1Seeded","InitialLayerClustersMask"),
clusterFilter = cms.string('ClusterFilterByAlgoAndSize'),
iteration_label = cms.string('CLUE3DHigh'),
max_cluster_size = cms.int32(9999),
Expand Down
Original file line number Diff line number Diff line change
@@ -1,9 +1,8 @@
import FWCore.ParameterSet.Config as cms

filteredLayerClustersCLUE3DHigh = cms.EDProducer("FilteredLayerClustersProducer",
LayerClusters = cms.InputTag("hgcalLayerClusters"),
LayerClustersInputMask = cms.InputTag("hgcalLayerClusters","InitialLayerClustersMask"),
algo_number = cms.int32(8),
LayerClusters = cms.InputTag("hgcalMergeLayerClusters"),
LayerClustersInputMask = cms.InputTag("hgcalMergeLayerClusters","InitialLayerClustersMask"),
clusterFilter = cms.string('ClusterFilterByAlgoAndSize'),
iteration_label = cms.string('CLUE3DHigh'),
max_cluster_size = cms.int32(9999),
Expand Down
Original file line number Diff line number Diff line change
@@ -1,9 +1,8 @@
import FWCore.ParameterSet.Config as cms

filteredLayerClustersEM = cms.EDProducer("FilteredLayerClustersProducer",
LayerClusters = cms.InputTag("hgcalLayerClusters"),
LayerClusters = cms.InputTag("hgcalMergeLayerClusters"),
LayerClustersInputMask = cms.InputTag("ticlTrackstersTrkEM"),
algo_number = cms.int32(8),
clusterFilter = cms.string('ClusterFilterByAlgoAndSizeAndLayerRange'),
iteration_label = cms.string('EM'),
max_cluster_size = cms.int32(9999),
Expand Down
Original file line number Diff line number Diff line change
@@ -1,9 +1,8 @@
import FWCore.ParameterSet.Config as cms

filteredLayerClustersHAD = cms.EDProducer("FilteredLayerClustersProducer",
LayerClusters = cms.InputTag("hgcalLayerClusters"),
LayerClusters = cms.InputTag("hgcalMergeLayerClusters"),
LayerClustersInputMask = cms.InputTag("ticlTrackstersTrk"),
algo_number = cms.int32(8),
clusterFilter = cms.string('ClusterFilterByAlgoAndSize'),
iteration_label = cms.string('HAD'),
max_cluster_size = cms.int32(9999),
Expand Down
Original file line number Diff line number Diff line change
@@ -1,9 +1,8 @@
import FWCore.ParameterSet.Config as cms

filteredLayerClustersTrkEM = cms.EDProducer("FilteredLayerClustersProducer",
LayerClusters = cms.InputTag("hgcalLayerClusters"),
LayerClustersInputMask = cms.InputTag("hgcalLayerClusters","InitialLayerClustersMask"),
algo_number = cms.int32(8),
LayerClusters = cms.InputTag("hgcalMergeLayerClusters"),
LayerClustersInputMask = cms.InputTag("hgcalMergeLayerClusters","InitialLayerClustersMask"),
clusterFilter = cms.string('ClusterFilterByAlgoAndSizeAndLayerRange'),
iteration_label = cms.string('TrkEM'),
max_cluster_size = cms.int32(9999),
Expand Down
Original file line number Diff line number Diff line change
@@ -1,9 +1,8 @@
import FWCore.ParameterSet.Config as cms

filteredLayerClustersTrk = cms.EDProducer("FilteredLayerClustersProducer",
LayerClusters = cms.InputTag("hgcalLayerClusters"),
LayerClusters = cms.InputTag("hgcalMergeLayerClusters"),
LayerClustersInputMask = cms.InputTag("ticlTrackstersEM"),
algo_number = cms.int32(8),
clusterFilter = cms.string('ClusterFilterByAlgoAndSize'),
iteration_label = cms.string('Trk'),
max_cluster_size = cms.int32(9999),
Expand Down
Original file line number Diff line number Diff line change
@@ -0,0 +1,61 @@
import FWCore.ParameterSet.Config as cms

hgcalLayerClustersEEL1Seeded = cms.EDProducer('HGCalLayerClusterProducer',
plugin = cms.PSet(
thresholdW0 = cms.vdouble(
2.9,
2.9,
2.9
),
positionDeltaRho2 = cms.double(1.69),
deltac = cms.vdouble(
1.3,
1.3,
1.3,
0.0315
),
dependSensor = cms.bool(True),
ecut = cms.double(3),
kappa = cms.double(9),
verbosity = cms.untracked.uint32(3),
dEdXweights = cms.vdouble(
0.0, 8.894541, 10.937907, 10.937907, 10.937907,
10.937907, 10.937907, 10.937907, 10.937907, 10.937907,
10.932882, 10.932882, 10.937907, 10.937907, 10.938169,
10.938169, 10.938169, 10.938169, 10.938169, 10.938169,
10.938169, 10.938169, 10.938169, 10.938169, 10.938169,
10.938169, 10.938169, 10.938169, 32.332097, 51.574301,
51.444192, 51.444192, 51.444192, 51.444192, 51.444192,
51.444192, 51.444192, 51.444192, 51.444192, 51.444192,
69.513118, 87.582044, 87.582044, 87.582044, 87.582044,
87.582044, 87.214571, 86.888309, 86.92952, 86.92952,
86.92952),
thicknessCorrection = cms.vdouble(
0.77, 0.77, 0.77, 0.84, 0.84,
0.84
),
sciThicknessCorrection = cms.double(0.9),
deltasi_index_regemfac = cms.int32(3),
maxNumberOfThickIndices = cms.uint32(6),
fcPerMip = cms.vdouble( 2.06, 3.43, 5.15, 2.06, 3.43, 5.15),
fcPerEle = cms.double(0),
noises = cms.vdouble( 2000.0, 2400.0, 2000.0, 2000.0, 2400.0,
2000.0),
noiseMip = cms.PSet(
scaleByDose = cms.bool(False),
scaleByDoseAlgo = cms.uint32(0),
scaleByDoseFactor = cms.double(1),
referenceIdark = cms.double(-1),
referenceXtalk = cms.double(-1),
noise_MIP = cms.double(0.01)
),
use2x2 = cms.bool(True),
type = cms.string('SiCLUE')

),
detector = cms.string('EE'),
recHits = cms.InputTag('hltRechitInRegionsHGCAL', 'HGCEERecHits'),
timeClname = cms.string('timeLayerCluster'),
nHitsTime = cms.uint32(3),
mightGet = cms.optional.untracked.vstring
)
Original file line number Diff line number Diff line change
@@ -0,0 +1,62 @@
import FWCore.ParameterSet.Config as cms

hgcalLayerClustersEE = cms.EDProducer('HGCalLayerClusterProducer',
plugin = cms.PSet(
thresholdW0 = cms.vdouble(
2.9,
2.9,
2.9
),
positionDeltaRho2 = cms.double(1.69),
deltac = cms.vdouble(
1.3,
1.3,
1.3,
0.0315
),
dependSensor = cms.bool(True),
ecut = cms.double(3),
kappa = cms.double(9),
verbosity = cms.untracked.uint32(3),
dEdXweights = cms.vdouble(
0.0, 8.894541, 10.937907, 10.937907, 10.937907,
10.937907, 10.937907, 10.937907, 10.937907, 10.937907,
10.932882, 10.932882, 10.937907, 10.937907, 10.938169,
10.938169, 10.938169, 10.938169, 10.938169, 10.938169,
10.938169, 10.938169, 10.938169, 10.938169, 10.938169,
10.938169, 10.938169, 10.938169, 32.332097, 51.574301,
51.444192, 51.444192, 51.444192, 51.444192, 51.444192,
51.444192, 51.444192, 51.444192, 51.444192, 51.444192,
69.513118, 87.582044, 87.582044, 87.582044, 87.582044,
87.582044, 87.214571, 86.888309, 86.92952, 86.92952,
86.92952),
thicknessCorrection = cms.vdouble(
0.77, 0.77, 0.77, 0.84, 0.84,
0.84
),
sciThicknessCorrection = cms.double(0.9),
deltasi_index_regemfac = cms.int32(3),
maxNumberOfThickIndices = cms.uint32(6),
fcPerMip = cms.vdouble( 2.06, 3.43, 5.15, 2.06, 3.43, 5.15),
fcPerEle = cms.double(0),
noises = cms.vdouble( 2000.0, 2400.0, 2000.0, 2000.0, 2400.0,
2000.0),
noiseMip = cms.PSet(
scaleByDose = cms.bool(False),
scaleByDoseAlgo = cms.uint32(0),
scaleByDoseFactor = cms.double(1),
referenceIdark = cms.double(-1),
referenceXtalk = cms.double(-1),
noise_MIP = cms.double(0.01)
),
use2x2 = cms.bool(True),
type = cms.string('SiCLUE')

),
detector = cms.string('EE'),
recHits = cms.InputTag('HGCalRecHit', 'HGCEERecHits'),
timeClname = cms.string('timeLayerCluster'),
nHitsTime = cms.uint32(3),
mightGet = cms.optional.untracked.vstring
)

Original file line number Diff line number Diff line change
@@ -0,0 +1,61 @@
import FWCore.ParameterSet.Config as cms

hgcalLayerClustersHSciL1Seeded = cms.EDProducer('HGCalLayerClusterProducer',
plugin = cms.PSet(
thresholdW0 = cms.vdouble(
2.9,
2.9,
2.9
),
positionDeltaRho2 = cms.double(1.69),
deltac = cms.vdouble(
1.3,
1.3,
1.3,
0.0315
),
dependSensor = cms.bool(True),
ecut = cms.double(3),
kappa = cms.double(9),
verbosity = cms.untracked.uint32(3),
dEdXweights = cms.vdouble(
0.0, 8.894541, 10.937907, 10.937907, 10.937907,
10.937907, 10.937907, 10.937907, 10.937907, 10.937907,
10.932882, 10.932882, 10.937907, 10.937907, 10.938169,
10.938169, 10.938169, 10.938169, 10.938169, 10.938169,
10.938169, 10.938169, 10.938169, 10.938169, 10.938169,
10.938169, 10.938169, 10.938169, 32.332097, 51.574301,
51.444192, 51.444192, 51.444192, 51.444192, 51.444192,
51.444192, 51.444192, 51.444192, 51.444192, 51.444192,
69.513118, 87.582044, 87.582044, 87.582044, 87.582044,
87.582044, 87.214571, 86.888309, 86.92952, 86.92952,
86.92952),
thicknessCorrection = cms.vdouble(
0.77, 0.77, 0.77, 0.84, 0.84,
0.84
),
sciThicknessCorrection = cms.double(0.9),
deltasi_index_regemfac = cms.int32(3),
maxNumberOfThickIndices = cms.uint32(6),
fcPerMip = cms.vdouble( 2.06, 3.43, 5.15, 2.06, 3.43, 5.15),
fcPerEle = cms.double(0),
noises = cms.vdouble( 2000.0, 2400.0, 2000.0, 2000.0, 2400.0,
2000.0),
noiseMip = cms.PSet(
scaleByDose = cms.bool(False),
scaleByDoseAlgo = cms.uint32(0),
scaleByDoseFactor = cms.double(1),
referenceIdark = cms.double(-1),
referenceXtalk = cms.double(-1),
noise_MIP = cms.double(0.01)
),
use2x2 = cms.bool(True),
type = cms.string('SciCLUE')

),
detector = cms.string('BH'),
recHits = cms.InputTag('hltRechitInRegionsHGCAL', 'HGCHEBRecHits'),
timeClname = cms.string('timeLayerCluster'),
nHitsTime = cms.uint32(3),
mightGet = cms.optional.untracked.vstring
)
Original file line number Diff line number Diff line change
@@ -0,0 +1,62 @@
import FWCore.ParameterSet.Config as cms

hgcalLayerClustersHSci = cms.EDProducer('HGCalLayerClusterProducer',
plugin = cms.PSet(
thresholdW0 = cms.vdouble(
2.9,
2.9,
2.9
),
positionDeltaRho2 = cms.double(1.69),
deltac = cms.vdouble(
1.3,
1.3,
1.3,
0.0315
),
dependSensor = cms.bool(True),
ecut = cms.double(3),
kappa = cms.double(9),
verbosity = cms.untracked.uint32(3),
dEdXweights = cms.vdouble(
0.0, 8.894541, 10.937907, 10.937907, 10.937907,
10.937907, 10.937907, 10.937907, 10.937907, 10.937907,
10.932882, 10.932882, 10.937907, 10.937907, 10.938169,
10.938169, 10.938169, 10.938169, 10.938169, 10.938169,
10.938169, 10.938169, 10.938169, 10.938169, 10.938169,
10.938169, 10.938169, 10.938169, 32.332097, 51.574301,
51.444192, 51.444192, 51.444192, 51.444192, 51.444192,
51.444192, 51.444192, 51.444192, 51.444192, 51.444192,
69.513118, 87.582044, 87.582044, 87.582044, 87.582044,
87.582044, 87.214571, 86.888309, 86.92952, 86.92952,
86.92952),
thicknessCorrection = cms.vdouble(
0.77, 0.77, 0.77, 0.84, 0.84,
0.84
),
sciThicknessCorrection = cms.double(0.9),
deltasi_index_regemfac = cms.int32(3),
maxNumberOfThickIndices = cms.uint32(6),
fcPerMip = cms.vdouble( 2.06, 3.43, 5.15, 2.06, 3.43, 5.15),
fcPerEle = cms.double(0),
noises = cms.vdouble( 2000.0, 2400.0, 2000.0, 2000.0, 2400.0,
2000.0),
noiseMip = cms.PSet(
scaleByDose = cms.bool(False),
scaleByDoseAlgo = cms.uint32(0),
scaleByDoseFactor = cms.double(1),
referenceIdark = cms.double(-1),
referenceXtalk = cms.double(-1),
noise_MIP = cms.double(0.01)
),
use2x2 = cms.bool(True),
type = cms.string('SciCLUE')

),
detector = cms.string('BH'),
recHits = cms.InputTag('HGCalRecHit', 'HGCHEBRecHits'),
timeClname = cms.string('timeLayerCluster'),
nHitsTime = cms.uint32(3),
mightGet = cms.optional.untracked.vstring
)

Loading

0 comments on commit 2c948f2

Please sign in to comment.