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  if (pfcConf.exists("positionCalc")) {
52  const edm::ParameterSet& acConf = pfcConf.getParameterSet("positionCalc");
53  const std::string& algoac = acConf.getParameter<std::string>("algoName");
54  positionCalc_ = PFCPositionCalculatorFactory::get()->create(algoac, acConf, cc);
55  }
56 
57  if (pfcConf.exists("allCellsPositionCalc")) {
58  const edm::ParameterSet& acConf = pfcConf.getParameterSet("allCellsPositionCalc");
59  const std::string& algoac = acConf.getParameter<std::string>("algoName");
60  allCellsPositionCalc_ = PFCPositionCalculatorFactory::get()->create(algoac, acConf, cc);
61  }
62  }
63  }
LegacyPFClusterProducer::positionCalc_
std::unique_ptr< PFCPositionCalculatorBase > positionCalc_
Definition: LegacyPFClusterProducer.cc:188
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
LegacyPFClusterProducer::pfClusterSoAToken_
const edm::EDGetTokenT< reco::PFClusterHostCollection > pfClusterSoAToken_
Definition: LegacyPFClusterProducer.cc:180
LegacyPFClusterProducer::pfRecHitFractionSoAToken_
const edm::EDGetTokenT< reco::PFRecHitFractionHostCollection > pfRecHitFractionSoAToken_
Definition: LegacyPFClusterProducer.cc:181
LegacyPFClusterProducer::InputPFRecHitSoA_Token_
const edm::EDGetTokenT< reco::PFRecHitHostCollection > InputPFRecHitSoA_Token_
Definition: LegacyPFClusterProducer.cc:182
gpuPixelDoublets::cc
uint32_t cc[maxCellsPerHit]
Definition: gpuFishbone.h:49
edm::ParameterSet::exists
bool exists(std::string const &parameterName) const
checks if a parameter exists
Definition: ParameterSet.cc:730
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:183
LegacyPFClusterProducer::allCellsPositionCalc_
std::unique_ptr< PFCPositionCalculatorBase > allCellsPositionCalc_
Definition: LegacyPFClusterProducer.cc:189
edm::ParameterSet::empty
bool empty() const
Definition: ParameterSet.h:202
LegacyPFClusterProducer::hcalCutsToken_
const edm::ESGetToken< HcalPFCuts, HcalPFCutsRcd > hcalCutsToken_
Definition: LegacyPFClusterProducer.cc:185
LegacyPFClusterProducer::recHitsLabel_
const edm::EDGetTokenT< reco::PFRecHitCollection > recHitsLabel_
Definition: LegacyPFClusterProducer.cc:184
get
#define get

Member Function Documentation

◆ fillDescriptions()

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

Definition at line 65 of file LegacyPFClusterProducer.cc.

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

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

◆ produce()

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

Definition at line 192 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_2024v11_cff::pfRecHits, positionCalc_, recHitsLabel_, singleTopDQM_cfi::setup, and groupFilesInBlocks::temp.

192  {
193  const reco::PFRecHitHostCollection& pfRecHits = event.get(InputPFRecHitSoA_Token_);
194 
195  HcalPFCuts const* paramPF = cutsFromDB_ ? &setup.getData(hcalCutsToken_) : nullptr;
196 
197  auto const& pfClusterSoA = event.get(pfClusterSoAToken_).const_view();
198  auto const& pfRecHitFractionSoA = event.get(pfRecHitFractionSoAToken_).const_view();
199 
200  int nRH = 0;
201  if (pfRecHits->metadata().size() != 0)
202  nRH = pfRecHits.view().size();
203  reco::PFClusterCollection out;
204  out.reserve(nRH);
205 
206  auto const rechitsHandle = event.getHandle(recHitsLabel_);
207 
208  if (nRH != 0) {
209  // Build PFClusters in legacy format
210  std::vector<int> nTopoSeeds(nRH, 0);
211 
212  for (int i = 0; i < pfClusterSoA.nSeeds(); i++) {
213  nTopoSeeds[pfClusterSoA[i].topoId()]++;
214  }
215 
216  // Looping over SoA PFClusters to produce legacy PFCluster collection
217  for (int i = 0; i < pfClusterSoA.nSeeds(); i++) {
218  unsigned int n = pfClusterSoA[i].seedRHIdx();
219  reco::PFCluster temp;
220  temp.setSeed((*rechitsHandle)[n].detId()); // Pulling the detId of this PFRecHit from the legacy format input
221  int offset = pfClusterSoA[i].rhfracOffset();
222  for (int k = offset; k < (offset + pfClusterSoA[i].rhfracSize()) && k >= 0;
223  k++) { // Looping over PFRecHits in the same topo cluster
224  if (pfRecHitFractionSoA[k].pfrhIdx() < nRH && pfRecHitFractionSoA[k].pfrhIdx() > -1 &&
225  pfRecHitFractionSoA[k].frac() > 0.0) {
226  const reco::PFRecHitRef& refhit = reco::PFRecHitRef(rechitsHandle, pfRecHitFractionSoA[k].pfrhIdx());
227  temp.addRecHitFraction(reco::PFRecHitFraction(refhit, pfRecHitFractionSoA[k].frac()));
228  }
229  }
230 
231  // Now PFRecHitFraction of this PFCluster is set. Now compute calculateAndSetPosition (energy, position etc)
232  if (nTopoSeeds[pfClusterSoA[i].topoId()] == 1 && allCellsPositionCalc_) {
233  allCellsPositionCalc_->calculateAndSetPosition(temp, paramPF);
234  } else {
235  positionCalc_->calculateAndSetPosition(temp, paramPF);
236  }
237  out.emplace_back(std::move(temp));
238  }
239  }
240 
241  event.emplace(legacyPfClustersToken_, std::move(out));
242 }
LegacyPFClusterProducer::positionCalc_
std::unique_ptr< PFCPositionCalculatorBase > positionCalc_
Definition: LegacyPFClusterProducer.cc:188
mps_fire.i
i
Definition: mps_fire.py:429
LegacyPFClusterProducer::pfClusterSoAToken_
const edm::EDGetTokenT< reco::PFClusterHostCollection > pfClusterSoAToken_
Definition: LegacyPFClusterProducer.cc:180
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:181
LegacyPFClusterProducer::InputPFRecHitSoA_Token_
const edm::EDGetTokenT< reco::PFRecHitHostCollection > InputPFRecHitSoA_Token_
Definition: LegacyPFClusterProducer.cc:182
reco::PFRecHitFraction
Fraction of a PFRecHit (rechits can be shared between several PFCluster&#39;s)
Definition: PFRecHitFraction.h:18
LegacyPFClusterProducer::legacyPfClustersToken_
const edm::EDPutTokenT< reco::PFClusterCollection > legacyPfClustersToken_
Definition: LegacyPFClusterProducer.cc:183
LegacyPFClusterProducer::allCellsPositionCalc_
std::unique_ptr< PFCPositionCalculatorBase > allCellsPositionCalc_
Definition: LegacyPFClusterProducer.cc:189
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
HLT_2024v11_cff.pfRecHits
pfRecHits
Definition: HLT_2024v11_cff.py:10422
LegacyPFClusterProducer::hcalCutsToken_
const edm::ESGetToken< HcalPFCuts, HcalPFCutsRcd > hcalCutsToken_
Definition: LegacyPFClusterProducer.cc:185
LegacyPFClusterProducer::recHitsLabel_
const edm::EDGetTokenT< reco::PFRecHitCollection > recHitsLabel_
Definition: LegacyPFClusterProducer.cc:184
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 189 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ cutsFromDB_

const bool LegacyPFClusterProducer::cutsFromDB_
private

Definition at line 186 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ hcalCutsToken_

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

Definition at line 185 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ InputPFRecHitSoA_Token_

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

Definition at line 182 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ legacyPfClustersToken_

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

Definition at line 183 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ pfClusterSoAToken_

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

Definition at line 180 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ pfRecHitFractionSoAToken_

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

Definition at line 181 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ positionCalc_

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

Definition at line 188 of file LegacyPFClusterProducer.cc.

Referenced by produce().

◆ recHitsLabel_

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

Definition at line 184 of file LegacyPFClusterProducer.cc.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.14