CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
LegacyPFClusterProducer Class Reference
Inheritance diagram for LegacyPFClusterProducer:
edm::stream::EDProducer<>

Public Member Functions

 LegacyPFClusterProducer (edm::ParameterSet const &config)
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProducer & operator= (const EDProducer &)=delete
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 

Private Member Functions

void produce (edm::Event &, const edm::EventSetup &) override
 

Private Attributes

std::unique_ptr< PFCPositionCalculatorBase > allCellsPositionCalc_
 
const bool cutsFromDB_
 
const edm::ESGetToken< HcalPFCuts, HcalPFCutsRcd > hcalCutsToken_
 
const edm::EDGetTokenT< reco::PFRecHitHostCollection > InputPFRecHitSoA_Token_
 
const edm::EDPutTokenT< reco::PFClusterCollection > legacyPfClustersToken_
 
const edm::EDGetTokenT< reco::PFClusterHostCollection > pfClusterSoAToken_
 
const edm::EDGetTokenT< reco::PFRecHitFractionHostCollection > pfRecHitFractionSoAToken_
 
std::unique_ptr< PFCPositionCalculatorBase > positionCalc_
 
const edm::EDGetTokenT< reco::PFRecHitCollection > recHitsLabel_
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDProducer<>
using CacheTypes = CacheContexts< T... >
 
using GlobalCache = typename CacheTypes::GlobalCache
 
using HasAbility = AbilityChecker< T... >
 
using InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
 
using LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
 
using LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
 
using LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
 
using RunCache = typename CacheTypes::RunCache
 
using RunContext = RunContextT< RunCache, GlobalCache >
 
using RunSummaryCache = typename CacheTypes::RunSummaryCache
 

Detailed Description

Definition at line 36 of file LegacyPFClusterProducer.cc.

Constructor & Destructor Documentation

◆ LegacyPFClusterProducer()

LegacyPFClusterProducer::LegacyPFClusterProducer ( edm::ParameterSet const &  config)
inline

Definition at line 38 of file LegacyPFClusterProducer.cc.

39  : pfClusterSoAToken_(consumes(config.getParameter<edm::InputTag>("src"))),
40  pfRecHitFractionSoAToken_(consumes(config.getParameter<edm::InputTag>("src"))),
41  InputPFRecHitSoA_Token_{consumes(config.getParameter<edm::InputTag>("PFRecHitsLabelIn"))},
42  legacyPfClustersToken_(produces()),
43  recHitsLabel_(consumes(config.getParameter<edm::InputTag>("recHitsSource"))),
44  hcalCutsToken_(esConsumes<HcalPFCuts, HcalPFCutsRcd>(edm::ESInputTag("", "withTopo"))),
45  cutsFromDB_(config.getParameter<bool>("usePFThresholdsFromDB")) {
46  edm::ConsumesCollector cc = consumesCollector();
47 
48  //setup pf cluster builder if requested
49  const edm::ParameterSet& pfcConf = config.getParameterSet("pfClusterBuilder");
50  if (!pfcConf.empty()) {
51  const auto& acConf = pfcConf.getParameterSet("positionCalc");
52  if (!acConf.empty()) {
53  const std::string& algoac = acConf.getParameter<std::string>("algoName");
54  if (!algoac.empty())
55  positionCalc_ = PFCPositionCalculatorFactory::get()->create(algoac, acConf, cc);
56  }
57 
58  const auto& acConf2 = pfcConf.getParameterSet("allCellsPositionCalc");
59  if (!acConf2.empty()) {
60  const std::string& algoac = acConf2.getParameter<std::string>("algoName");
61  if (!algoac.empty())
62  allCellsPositionCalc_ = PFCPositionCalculatorFactory::get()->create(algoac, acConf2, cc);
63  }
64  }
65  }
LegacyPFClusterProducer::positionCalc_
std::unique_ptr< PFCPositionCalculatorBase > positionCalc_
Definition: LegacyPFClusterProducer.cc:195
LegacyPFClusterProducer::pfClusterSoAToken_
const edm::EDGetTokenT< reco::PFClusterHostCollection > pfClusterSoAToken_
Definition: LegacyPFClusterProducer.cc:187
LegacyPFClusterProducer::pfRecHitFractionSoAToken_
const edm::EDGetTokenT< reco::PFRecHitFractionHostCollection > pfRecHitFractionSoAToken_
Definition: LegacyPFClusterProducer.cc:188
LegacyPFClusterProducer::InputPFRecHitSoA_Token_
const edm::EDGetTokenT< reco::PFRecHitHostCollection > InputPFRecHitSoA_Token_
Definition: LegacyPFClusterProducer.cc:189
gpuPixelDoublets::cc
uint32_t cc[maxCellsPerHit]
Definition: gpuFishbone.h:49
edm::ParameterSet::getParameterSet
ParameterSet const & getParameterSet(std::string const &) const
Definition: ParameterSet.cc:2204
config
Definition: config.py:1
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
LegacyPFClusterProducer::legacyPfClustersToken_
const edm::EDPutTokenT< reco::PFClusterCollection > legacyPfClustersToken_
Definition: LegacyPFClusterProducer.cc:190
LegacyPFClusterProducer::allCellsPositionCalc_
std::unique_ptr< PFCPositionCalculatorBase > allCellsPositionCalc_
Definition: LegacyPFClusterProducer.cc:196
edm::ParameterSet::empty
bool empty() const
Definition: ParameterSet.h:202
LegacyPFClusterProducer::hcalCutsToken_
const edm::ESGetToken< HcalPFCuts, HcalPFCutsRcd > hcalCutsToken_
Definition: LegacyPFClusterProducer.cc:192
LegacyPFClusterProducer::recHitsLabel_
const edm::EDGetTokenT< reco::PFRecHitCollection > recHitsLabel_
Definition: LegacyPFClusterProducer.cc:191
get
#define get

Member Function Documentation

◆ fillDescriptions()

static void LegacyPFClusterProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
inlinestatic

Definition at line 67 of file LegacyPFClusterProducer.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, MillePedeFileConverter_cfg::e, ProducerED_cfi::InputTag, HLT_2024v14_cff::pfClusterBuilder, and AlCaHLTBitMon_QueryRunRegistry::string.

67  {
68  edm::ParameterSetDescription desc;
69  desc.add<edm::InputTag>("src", edm::InputTag("pfClusterSoAProducer"));
70  desc.add<edm::InputTag>("PFRecHitsLabelIn", edm::InputTag("pfRecHitSoAProducerHCAL"));
71  desc.add<edm::InputTag>("recHitsSource", edm::InputTag("legacyPFRecHitProducer"));
72  desc.add<bool>("usePFThresholdsFromDB", true);
73  {
74  edm::ParameterSetDescription pfClusterBuilder;
75  pfClusterBuilder.add<unsigned int>("maxIterations", 5);
76  pfClusterBuilder.add<double>("minFracTot", 1e-20);
77  pfClusterBuilder.add<double>("minFractionToKeep", 1e-7);
78  pfClusterBuilder.add<bool>("excludeOtherSeeds", true);
79  pfClusterBuilder.add<double>("showerSigma", 10.);
80  pfClusterBuilder.add<double>("stoppingTolerance", 1e-8);
81  pfClusterBuilder.add<double>("timeSigmaEB", 10.);
82  pfClusterBuilder.add<double>("timeSigmaEE", 10.);
83  pfClusterBuilder.add<double>("maxNSigmaTime", 10.);
84  pfClusterBuilder.add<double>("minChi2Prob", 0.);
85  pfClusterBuilder.add<bool>("clusterTimeResFromSeed", false);
86  pfClusterBuilder.add<std::string>("algoName", "");
87  pfClusterBuilder.add<edm::ParameterSetDescription>("positionCalcForConvergence", {});
88  {
89  edm::ParameterSetDescription validator;
90  validator.add<std::string>("detector", "");
91  validator.add<std::vector<int>>("depths", {});
92  validator.add<std::vector<double>>("recHitEnergyNorm", {});
93  std::vector<edm::ParameterSet> vDefaults(2);
94  vDefaults[0].addParameter<std::string>("detector", "HCAL_BARREL1");
95  vDefaults[0].addParameter<std::vector<int>>("depths", {1, 2, 3, 4});
96  vDefaults[0].addParameter<std::vector<double>>("recHitEnergyNorm", {0.1, 0.2, 0.3, 0.3});
97  vDefaults[1].addParameter<std::string>("detector", "HCAL_ENDCAP");
98  vDefaults[1].addParameter<std::vector<int>>("depths", {1, 2, 3, 4, 5, 6, 7});
99  vDefaults[1].addParameter<std::vector<double>>("recHitEnergyNorm", {0.1, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2});
100  pfClusterBuilder.addVPSet("recHitEnergyNorms", validator, vDefaults);
101  }
102  {
103  edm::ParameterSetDescription bar;
104  bar.add<std::string>("algoName", "Basic2DGenericPFlowPositionCalc");
105  bar.add<double>("minFractionInCalc", 1e-9);
106  bar.add<int>("posCalcNCrystals", 5);
107  {
108  edm::ParameterSetDescription validator;
109  validator.add<std::string>("detector", "");
110  validator.add<std::vector<int>>("depths", {});
111  validator.add<std::vector<double>>("logWeightDenominator", {});
112  std::vector<edm::ParameterSet> vDefaults(2);
113  vDefaults[0].addParameter<std::string>("detector", "HCAL_BARREL1");
114  vDefaults[0].addParameter<std::vector<int>>("depths", {1, 2, 3, 4});
115  vDefaults[0].addParameter<std::vector<double>>("logWeightDenominator", {0.1, 0.2, 0.3, 0.3});
116  vDefaults[1].addParameter<std::string>("detector", "HCAL_ENDCAP");
117  vDefaults[1].addParameter<std::vector<int>>("depths", {1, 2, 3, 4, 5, 6, 7});
118  vDefaults[1].addParameter<std::vector<double>>("logWeightDenominator", {0.1, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2});
119  bar.addVPSet("logWeightDenominatorByDetector", validator, vDefaults);
120  }
121  bar.add<double>("minAllowedNormalization", 1e-9);
122  bar.add<edm::ParameterSetDescription>("timeResolutionCalcBarrel", {});
123  bar.add<edm::ParameterSetDescription>("timeResolutionCalcEndcap", {});
124  pfClusterBuilder.add("positionCalc", bar);
125  }
126  {
127  edm::ParameterSetDescription bar;
128  bar.add<std::string>("algoName", "Basic2DGenericPFlowPositionCalc");
129  bar.add<double>("minFractionInCalc", 1e-9);
130  bar.add<int>("posCalcNCrystals", -1);
131  {
132  edm::ParameterSetDescription validator;
133  validator.add<std::string>("detector", "");
134  validator.add<std::vector<int>>("depths", {});
135  validator.add<std::vector<double>>("logWeightDenominator", {});
136  std::vector<edm::ParameterSet> vDefaults(2);
137  vDefaults[0].addParameter<std::string>("detector", "HCAL_BARREL1");
138  vDefaults[0].addParameter<std::vector<int>>("depths", {1, 2, 3, 4});
139  vDefaults[0].addParameter<std::vector<double>>("logWeightDenominator", {0.1, 0.2, 0.3, 0.3});
140  vDefaults[1].addParameter<std::string>("detector", "HCAL_ENDCAP");
141  vDefaults[1].addParameter<std::vector<int>>("depths", {1, 2, 3, 4, 5, 6, 7});
142  vDefaults[1].addParameter<std::vector<double>>("logWeightDenominator", {0.1, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2});
143  bar.addVPSet("logWeightDenominatorByDetector", validator, vDefaults);
144  }
145  bar.add<double>("minAllowedNormalization", 1e-9);
146  bar.add<edm::ParameterSetDescription>("timeResolutionCalcBarrel", {});
147  bar.add<edm::ParameterSetDescription>("timeResolutionCalcEndcap", {});
148  pfClusterBuilder.add("allCellsPositionCalc", bar);
149  }
150  {
151  edm::ParameterSetDescription bar;
152  bar.add<double>("corrTermLowE", 0.);
153  bar.add<double>("threshLowE", 6.);
154  bar.add<double>("noiseTerm", 21.86);
155  bar.add<double>("constantTermLowE", 4.24);
156  bar.add<double>("noiseTermLowE", 8.);
157  bar.add<double>("threshHighE", 15.);
158  bar.add<double>("constantTerm", 2.82);
159  pfClusterBuilder.add("timeResolutionCalcBarrel", bar);
160  }
161  {
162  edm::ParameterSetDescription bar;
163  bar.add<double>("corrTermLowE", 0.);
164  bar.add<double>("threshLowE", 6.);
165  bar.add<double>("noiseTerm", 21.86);
166  bar.add<double>("constantTermLowE", 4.24);
167  bar.add<double>("noiseTermLowE", 8.);
168  bar.add<double>("threshHighE", 15.);
169  bar.add<double>("constantTerm", 2.82);
170  pfClusterBuilder.add("timeResolutionCalcEndcap", bar);
171  }
172  {
173  edm::ParameterSetDescription bar;
174  pfClusterBuilder.add("positionReCalc", bar);
175  }
176  {
177  edm::ParameterSetDescription bar;
178  pfClusterBuilder.add("energyCorrector", bar);
179  }
180  desc.add("pfClusterBuilder", pfClusterBuilder);
181  }
182  descriptions.addWithDefaultLabel(desc);
183  }
HLT_2024v14_cff.pfClusterBuilder
pfClusterBuilder
Definition: HLT_2024v14_cff.py:5446
edm::ConfigurationDescriptions::addWithDefaultLabel
void addWithDefaultLabel(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:87
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95

◆ produce()

void LegacyPFClusterProducer::produce ( edm::Event &  event,
const edm::EventSetup &  setup 
)
overrideprivate

Definition at line 199 of file LegacyPFClusterProducer.cc.

References allCellsPositionCalc_, cutsFromDB_, hcalRecHitTable_cff::detId, DivergingColor::frac, hcalCutsToken_, mps_fire::i, InputPFRecHitSoA_Token_, dqmdumpme::k, legacyPfClustersToken_, eostools::move(), dqmiodumpmetadata::n, hltrates_dqm_sourceclient-live_cfg::offset, MillePedeFileConverter_cfg::out, pfClusterSoAToken_, pfRecHitFractionSoAToken_, HLT_2024v14_cff::pfRecHits, positionCalc_, recHitsLabel_, singleTopDQM_cfi::setup, and groupFilesInBlocks::temp.

199  {
200  const reco::PFRecHitHostCollection& pfRecHits = event.get(InputPFRecHitSoA_Token_);
201 
202  HcalPFCuts const* paramPF = cutsFromDB_ ? &setup.getData(hcalCutsToken_) : nullptr;
203 
204  auto const& pfClusterSoA = event.get(pfClusterSoAToken_).const_view();
205  auto const& pfRecHitFractionSoA = event.get(pfRecHitFractionSoAToken_).const_view();
206 
207  int nRH = 0;
208  if (pfRecHits->metadata().size() != 0)
209  nRH = pfRecHits.view().size();
210  reco::PFClusterCollection out;
211  out.reserve(nRH);
212 
213  auto const rechitsHandle = event.getHandle(recHitsLabel_);
214 
215  if (nRH != 0) {
216  // Build PFClusters in legacy format
217  std::vector<int> nTopoSeeds(nRH, 0);
218 
219  for (int i = 0; i < pfClusterSoA.nSeeds(); i++) {
220  nTopoSeeds[pfClusterSoA[i].topoId()]++;
221  }
222 
223  // Looping over SoA PFClusters to produce legacy PFCluster collection
224  for (int i = 0; i < pfClusterSoA.nSeeds(); i++) {
225  unsigned int n = pfClusterSoA[i].seedRHIdx();
226  reco::PFCluster temp;
227  temp.setSeed((*rechitsHandle)[n].detId()); // Pulling the detId of this PFRecHit from the legacy format input
228  int offset = pfClusterSoA[i].rhfracOffset();
229  for (int k = offset; k < (offset + pfClusterSoA[i].rhfracSize()) && k >= 0;
230  k++) { // Looping over PFRecHits in the same topo cluster
231  if (pfRecHitFractionSoA[k].pfrhIdx() < nRH && pfRecHitFractionSoA[k].pfrhIdx() > -1 &&
232  pfRecHitFractionSoA[k].frac() > 0.0) {
233  const reco::PFRecHitRef& refhit = reco::PFRecHitRef(rechitsHandle, pfRecHitFractionSoA[k].pfrhIdx());
234  temp.addRecHitFraction(reco::PFRecHitFraction(refhit, pfRecHitFractionSoA[k].frac()));
235  }
236  }
237 
238  // Now PFRecHitFraction of this PFCluster is set. Now compute calculateAndSetPosition (energy, position etc)
239  if (nTopoSeeds[pfClusterSoA[i].topoId()] == 1 && allCellsPositionCalc_) {
240  allCellsPositionCalc_->calculateAndSetPosition(temp, paramPF);
241  } else {
242  positionCalc_->calculateAndSetPosition(temp, paramPF);
243  }
244  out.emplace_back(std::move(temp));
245  }
246  }
247 
248  event.emplace(legacyPfClustersToken_, std::move(out));
249 }
LegacyPFClusterProducer::positionCalc_
std::unique_ptr< PFCPositionCalculatorBase > positionCalc_
Definition: LegacyPFClusterProducer.cc:195
mps_fire.i
i
Definition: mps_fire.py:429
LegacyPFClusterProducer::pfClusterSoAToken_
const edm::EDGetTokenT< reco::PFClusterHostCollection > pfClusterSoAToken_
Definition: LegacyPFClusterProducer.cc:187
reco::PFCluster
Particle flow cluster, see clustering algorithm in PFClusterAlgo.
Definition: PFCluster.h:42
LegacyPFClusterProducer::pfRecHitFractionSoAToken_
const edm::EDGetTokenT< reco::PFRecHitFractionHostCollection > pfRecHitFractionSoAToken_
Definition: LegacyPFClusterProducer.cc:188
LegacyPFClusterProducer::InputPFRecHitSoA_Token_
const edm::EDGetTokenT< reco::PFRecHitHostCollection > InputPFRecHitSoA_Token_
Definition: LegacyPFClusterProducer.cc:189
reco::PFRecHitFraction
Fraction of a PFRecHit (rechits can be shared between several PFCluster&#39;s)
Definition: PFRecHitFraction.h:18
HLT_2024v14_cff.pfRecHits
pfRecHits
Definition: HLT_2024v14_cff.py:10488
LegacyPFClusterProducer::legacyPfClustersToken_
const edm::EDPutTokenT< reco::PFClusterCollection > legacyPfClustersToken_
Definition: LegacyPFClusterProducer.cc:190
LegacyPFClusterProducer::allCellsPositionCalc_
std::unique_ptr< PFCPositionCalculatorBase > allCellsPositionCalc_
Definition: LegacyPFClusterProducer.cc:196
singleTopDQM_cfi.setup
setup
SETUP
Definition: singleTopDQM_cfi.py:37
reco::PFRecHitRef
edm::Ref< PFRecHitCollection > PFRecHitRef
persistent reference to PFRecHit objects
Definition: PFRecHitFwd.h:15
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
MillePedeFileConverter_cfg.out
out
Definition: MillePedeFileConverter_cfg.py:31
LegacyPFClusterProducer::hcalCutsToken_
const edm::ESGetToken< HcalPFCuts, HcalPFCutsRcd > hcalCutsToken_
Definition: LegacyPFClusterProducer.cc:192
LegacyPFClusterProducer::recHitsLabel_
const edm::EDGetTokenT< reco::PFRecHitCollection > recHitsLabel_
Definition: LegacyPFClusterProducer.cc:191
reco::PFClusterCollection
std::vector< PFCluster > PFClusterCollection
collection of PFCluster objects
Definition: PFClusterFwd.h:9
eostools.move
def move(src, dest)
Definition: eostools.py:511
hltrates_dqm_sourceclient-live_cfg.offset
offset
Definition: hltrates_dqm_sourceclient-live_cfg.py:83

Member Data Documentation

◆ allCellsPositionCalc_

std::unique_ptr<PFCPositionCalculatorBase> LegacyPFClusterProducer::allCellsPositionCalc_
private

Definition at line 196 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ cutsFromDB_

const bool LegacyPFClusterProducer::cutsFromDB_
private

Definition at line 193 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ hcalCutsToken_

const edm::ESGetToken<HcalPFCuts, HcalPFCutsRcd> LegacyPFClusterProducer::hcalCutsToken_
private

Definition at line 192 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ InputPFRecHitSoA_Token_

const edm::EDGetTokenT<reco::PFRecHitHostCollection> LegacyPFClusterProducer::InputPFRecHitSoA_Token_
private

Definition at line 189 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ legacyPfClustersToken_

const edm::EDPutTokenT<reco::PFClusterCollection> LegacyPFClusterProducer::legacyPfClustersToken_
private

Definition at line 190 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ pfClusterSoAToken_

const edm::EDGetTokenT<reco::PFClusterHostCollection> LegacyPFClusterProducer::pfClusterSoAToken_
private

Definition at line 187 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ pfRecHitFractionSoAToken_

const edm::EDGetTokenT<reco::PFRecHitFractionHostCollection> LegacyPFClusterProducer::pfRecHitFractionSoAToken_
private

Definition at line 188 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ positionCalc_

std::unique_ptr<PFCPositionCalculatorBase> LegacyPFClusterProducer::positionCalc_
private

Definition at line 195 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ recHitsLabel_

const edm::EDGetTokenT<reco::PFRecHitCollection> LegacyPFClusterProducer::recHitsLabel_
private

Definition at line 191 of file LegacyPFClusterProducer.cc.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.14