|
|
Definition at line 77 of file GsfElectronProducer.cc.
◆ GsfElectronProducer()
Definition at line 291 of file GsfElectronProducer.cc.
311 if (
cfg.getParameter<
bool>(
"fillConvVtxFitProb"))
318 cfg.getParameter<
bool>(
"ecalDrivenEcalEnergyFromClassBasedParameterization");
320 cfg.getParameter<
bool>(
"ecalDrivenEcalErrorFromClassBasedParameterization");
322 cfg.getParameter<
bool>(
"pureTrackerDrivenEcalErrorFromSimpleParameterization");
366 auto const& flagnamesbarrel =
cfg.getParameter<std::vector<std::string>>(
"recHitFlagsToBeExcludedBarrel");
367 recHitsCfg.recHitFlagsToBeExcludedBarrel = StringToEnumValue<EcalRecHit::Flags>(flagnamesbarrel);
368 auto const& flagnamesendcaps =
cfg.getParameter<std::vector<std::string>>(
"recHitFlagsToBeExcludedEndcaps");
369 recHitsCfg.recHitFlagsToBeExcludedEndcaps = StringToEnumValue<EcalRecHit::Flags>(flagnamesendcaps);
370 auto const& severitynamesbarrel =
cfg.getParameter<std::vector<std::string>>(
"recHitSeverityToBeExcludedBarrel");
371 recHitsCfg.recHitSeverityToBeExcludedBarrel =
372 StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesbarrel);
373 auto const& severitynamesendcaps =
cfg.getParameter<std::vector<std::string>>(
"recHitSeverityToBeExcludedEndcaps");
374 recHitsCfg.recHitSeverityToBeExcludedEndcaps =
375 StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesendcaps);
380 .etMinHcal =
cfg.getParameter<
double>(
"etMinHcal"),
381 .intRadiusEcalBarrel =
cfg.getParameter<
double>(
"intRadiusEcalBarrel"),
382 .intRadiusEcalEndcaps =
cfg.getParameter<
double>(
"intRadiusEcalEndcaps"),
383 .jurassicWidth =
cfg.getParameter<
double>(
"jurassicWidth"),
385 .eMinBarrel =
cfg.getParameter<
double>(
"eMinBarrel"),
386 .etMinEndcaps =
cfg.getParameter<
double>(
"etMinEndcaps"),
387 .eMinEndcaps =
cfg.getParameter<
double>(
"eMinEndcaps"),
389 .useNumCrystals =
cfg.getParameter<
bool>(
"useNumCrystals")};
393 .ecalWeightsFromDB =
cfg.getParameter<
bool>(
"ecalWeightsFromDB"),
394 .ecalRegressionWeightFiles =
cfg.getParameter<std::vector<std::string>>(
"ecalRefinedRegressionWeightFiles"),
396 cfg.getParameter<std::vector<std::string>>(
"combinationRegressionWeightLabels"),
397 .combinationWeightsFromDB =
cfg.getParameter<
bool>(
"combinationWeightsFromDB"),
398 .combinationRegressionWeightFiles =
399 cfg.getParameter<std::vector<std::string>>(
"combinationRegressionWeightFile")};
402 algo_ = std::make_unique<GsfElectronAlgo>(
416 consumesCollector());
References looper::cfg.
◆ checkEcalSeedingParameters()
void GsfElectronProducer::checkEcalSeedingParameters |
( |
edm::ParameterSet const & |
pset | ) |
|
|
private |
◆ fillDescriptions()
Definition at line 121 of file GsfElectronProducer.cc.
128 desc.add<
edm::InputTag>(
"barrelRecHitCollectionTag", {
"ecalRecHit",
"EcalRecHitsEB"});
129 desc.add<
edm::InputTag>(
"endcapRecHitCollectionTag", {
"ecalRecHit",
"EcalRecHitsEE"});
135 desc.add<
bool>(
"useDefaultEnergyCorrection",
true);
136 desc.add<
bool>(
"useCombinationRegression",
false);
137 desc.add<
bool>(
"ecalDrivenEcalEnergyFromClassBasedParameterization",
true);
138 desc.add<
bool>(
"ecalDrivenEcalErrorFromClassBasedParameterization",
true);
139 desc.add<
bool>(
"applyPreselection",
false);
140 desc.add<
bool>(
"useEcalRegression",
false);
141 desc.add<
bool>(
"applyAmbResolution",
false);
142 desc.add<
bool>(
"ignoreNotPreselected",
true);
143 desc.add<
bool>(
"useGsfPfRecTracks",
true);
144 desc.add<
bool>(
"pureTrackerDrivenEcalErrorFromSimpleParameterization",
true);
145 desc.add<
unsigned int>(
"ambSortingStrategy", 1);
146 desc.add<
unsigned int>(
"ambClustersOverlapStrategy", 1);
147 desc.add<
bool>(
"fillConvVtxFitProb",
true);
148 desc.add<
bool>(
"resetMvaValuesUsingPFCandidates",
false);
151 desc.add<std::vector<std::string>>(
"recHitFlagsToBeExcludedBarrel");
152 desc.add<std::vector<std::string>>(
"recHitFlagsToBeExcludedEndcaps");
153 desc.add<std::vector<std::string>>(
"recHitSeverityToBeExcludedBarrel");
154 desc.add<std::vector<std::string>>(
"recHitSeverityToBeExcludedEndcaps");
157 desc.add<
bool>(
"checkHcalStatus",
true);
159 desc.add<std::vector<double>>(
"recHitEThresholdHB", {0., 0., 0., 0.});
160 desc.add<std::vector<double>>(
"recHitEThresholdHE", {0., 0., 0., 0., 0., 0., 0.});
161 desc.add<
int>(
"maxHcalRecHitSeverity", 999999);
162 desc.add<
bool>(
"hcalRun2EffDepth",
false);
169 desc.add<
bool>(
"useNumCrystals",
true);
170 desc.add<
double>(
"etMinBarrel", 0.0);
171 desc.add<
double>(
"etMinEndcaps", 0.11);
172 desc.add<
double>(
"etMinHcal", 0.0);
173 desc.add<
double>(
"eMinBarrel", 0.095);
174 desc.add<
double>(
"eMinEndcaps", 0.0);
175 desc.add<
double>(
"intRadiusEcalBarrel", 3.0);
176 desc.add<
double>(
"intRadiusEcalEndcaps", 3.0);
177 desc.add<
double>(
"intRadiusHcal", 0.15);
178 desc.add<
double>(
"jurassicWidth", 1.5);
179 desc.add<
bool>(
"vetoClustered",
false);
182 desc.add<
bool>(
"ctfTracksCheck",
true);
185 desc.add<
double>(
"MaxElePtForOnlyMVA", 50.0);
186 desc.add<
double>(
"PreSelectMVA", -0.1);
190 psd0.
add<
double>(
"minSCEtBarrel", 4.0);
191 psd0.add<
double>(
"minSCEtEndcaps", 4.0);
192 psd0.add<
double>(
"minEOverPBarrel", 0.0);
193 psd0.add<
double>(
"minEOverPEndcaps", 0.0);
194 psd0.add<
double>(
"maxEOverPBarrel", 999999999.0);
195 psd0.add<
double>(
"maxEOverPEndcaps", 999999999.0);
196 psd0.add<
double>(
"maxDeltaEtaBarrel", 0.02);
197 psd0.add<
double>(
"maxDeltaEtaEndcaps", 0.02);
198 psd0.add<
double>(
"maxDeltaPhiBarrel", 0.15);
199 psd0.add<
double>(
"maxDeltaPhiEndcaps", 0.15);
200 psd0.add<
double>(
"hOverEConeSize", 0.15);
201 psd0.add<
double>(
"maxHOverEBarrelCone", 0.15);
202 psd0.add<
double>(
"maxHOverEEndcapsCone", 0.15);
203 psd0.add<
double>(
"maxHBarrelCone", 0.0);
204 psd0.add<
double>(
"maxHEndcapsCone", 0.0);
205 psd0.add<
double>(
"maxHOverEBarrelBc", 0.15);
206 psd0.add<
double>(
"maxHOverEEndcapsBc", 0.15);
207 psd0.add<
double>(
"maxHBarrelBc", 0.0);
208 psd0.add<
double>(
"maxHEndcapsBc", 0.0);
209 psd0.add<
double>(
"maxSigmaIetaIetaBarrel", 999999999.0);
210 psd0.add<
double>(
"maxSigmaIetaIetaEndcaps", 999999999.0);
211 psd0.add<
double>(
"maxFbremBarrel", 999999999.0);
212 psd0.add<
double>(
"maxFbremEndcaps", 999999999.0);
213 psd0.add<
bool>(
"isBarrel",
false);
214 psd0.add<
bool>(
"isEndcaps",
false);
215 psd0.add<
bool>(
"isFiducial",
false);
216 psd0.add<
bool>(
"seedFromTEC",
true);
217 psd0.add<
double>(
"maxTIP", 999999999.0);
225 desc.add<
std::string>(
"crackCorrectionFunction",
"EcalClusterCrackCorrection");
227 desc.add<
bool>(
"ecalWeightsFromDB",
true);
228 desc.add<std::vector<std::string>>(
"ecalRefinedRegressionWeightFiles", {})
229 ->setComment(
"if not from DB. Otherwise, keep empty");
230 desc.add<
bool>(
"combinationWeightsFromDB",
true);
231 desc.add<std::vector<std::string>>(
"combinationRegressionWeightFile", {})
232 ->setComment(
"if not from DB. Otherwise, keep empty");
235 desc.add<std::vector<std::string>>(
"ecalRefinedRegressionWeightLabels", {});
236 desc.add<std::vector<std::string>>(
"combinationRegressionWeightLabels", {});
238 desc.add<std::vector<std::string>>(
239 "ElecMVAFilesString",
241 "RecoEgamma/ElectronIdentification/data/TMVA_Category_BDTSimpleCat_10_17Feb2011.weights.xml",
242 "RecoEgamma/ElectronIdentification/data/TMVA_Category_BDTSimpleCat_12_17Feb2011.weights.xml",
243 "RecoEgamma/ElectronIdentification/data/TMVA_Category_BDTSimpleCat_20_17Feb2011.weights.xml",
244 "RecoEgamma/ElectronIdentification/data/TMVA_Category_BDTSimpleCat_22_17Feb2011.weights.xml",
246 desc.add<std::vector<std::string>>(
247 "SoftElecMVAFilesString",
249 "RecoEgamma/ElectronIdentification/data/TMVA_BDTSoftElectrons_7Feb2014.weights.xml",
252 descriptions.
add(
"gsfElectronProducerDefault",
desc);
References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, EgammaLocalCovParamDefaults::kMultThresEB, EgammaLocalCovParamDefaults::kMultThresEE, EleTkIsolFromCands::pSetDescript(), and AlCaHLTBitMon_QueryRunRegistry::string.
◆ globalEndJob()
◆ initializeGlobalCache()
Definition at line 83 of file GsfElectronProducer.cc.
84 return std::make_unique<GsfElectronAlgo::HeavyObjectCache>(conf);
◆ isPreselected()
◆ produce()
Definition at line 571 of file GsfElectronProducer.cc.
576 if (!
seeds.isValid()) {
578 <<
"Cannot check consistency of parameters with ecal seeding ones,"
579 <<
" because the original collection of seeds is not any more available.";
590 el.setMvaInput(gsfMVAInputMap.find(el.gsfTrack())->
second);
594 logElectrons(
electrons,
event,
"GsfElectronAlgo Info (before preselection)");
600 logElectrons(
electrons,
event,
"GsfElectronAlgo Info (after preselection)");
607 logElectrons(
electrons,
event,
"GsfElectronAlgo Info (after amb. solving)");
612 ele.hcalToRun2EffDepth();
References algo_, reco::GsfElectron::ambiguous(), GsfElectronAlgo::StrategyConfiguration::applyAmbResolution, GsfElectronAlgo::StrategyConfiguration::applyPreselection, checkEcalSeedingParameters(), ecalSeedingParametersChecked_, egmPFCandidateCollection_, electronPutToken_, pwdgSkimBPark_cfi::electrons, hcalRun2EffDepth_, GsfElectronAlgo::StrategyConfiguration::ignoreNotPreselected, inputCfg_, isPreselected(), eostools::move(), edm::parameterSet(), resetMvaValuesUsingPFCandidates_, edm::second(), DetachedQuadStep_cff::seeds, GsfElectronAlgo::Tokens::seedsTag, setAmbiguityData(), singleTopDQM_cfi::setup, strategyCfg_, and GsfElectronAlgo::Tokens::vtxCollectionTag.
◆ setAmbiguityData()
Definition at line 456 of file GsfElectronProducer.cc.
471 throw cms::Exception(
"GsfElectronAlgo|UnknownAmbiguitySortingStrategy")
486 if (gsfPfRecTrack.gsfTrackRef() ==
e1.gsfTrack()) {
488 edm::LogWarning(
"GsfElectronAlgo") <<
"associated gsfPfRecTrack already found";
491 for (
auto const& duplicate : gsfPfRecTrack.convBremGsfPFRecTrackRef()) {
492 e1.addAmbiguousGsfTrack(duplicate->gsfTrackRef());
509 LogDebug(
"GsfElectronAlgo") <<
"Blessing electron with E/P " <<
e1->eSuperClusterOverP() <<
", cluster "
510 << scRef1.
get() <<
" & track " <<
e1->gsfTrack().get();
512 for (
auto e2 =
e1 + 1; e2 !=
electrons.end(); ++e2) {
522 bool sameCluster =
false;
524 sameCluster = (scRef1 == scRef2);
534 throw cms::Exception(
"GsfElectronAlgo|UnknownAmbiguityClustersOverlapStrategy")
540 LogDebug(
"GsfElectronAlgo") <<
"Discarding electron with E/P " << e2->eSuperClusterOverP() <<
", cluster "
541 << scRef2.
get() <<
" and track " << e2->gsfTrack().get();
542 e1->addAmbiguousGsfTrack(e2->gsfTrack());
543 e2->setAmbiguous(
true);
544 }
else if (
e1->gsfTrack() == e2->gsfTrack()) {
545 edm::LogWarning(
"GsfElectronAlgo") <<
"Forgetting electron with E/P " << e2->eSuperClusterOverP()
546 <<
", cluster " << scRef2.
get() <<
" and track " << e2->gsfTrack().get();
547 e2->setAmbiguous(
true);
References GsfElectronAlgo::StrategyConfiguration::ambClustersOverlapStrategy, GsfElectronAlgo::StrategyConfiguration::ambSortingStrategy, GsfElectronAlgo::Tokens::barrelRecHitCollection, HLT_FULL_cff::barrelRecHits, GsfElectronAlgo::Tokens::beamSpotTag, StorageManager_cfg::e1, HPSPFTauProducerPuppi_cfi::electron, pwdgSkimBPark_cfi::electrons, photonAnalyzer_cfi::eMin, GsfElectronAlgo::Tokens::endcapRecHitCollection, HLT_FULL_cff::endcapRecHits, PVValHelper::eta, Exception, newFWLiteAna::found, edm::Ref< C, T, F >::get(), uncleanedOnlyElectronSequence_cff::gsfPfRecTracks, gsfPfRecTracksTag_, gsfElectrons_cfi::ignoreNotPreselected, inputCfg_, egamma::isBetterElectron(), egamma::isInnermostElectron(), isPreselected(), LogDebug, egamma::sharedEnergy(), jetUpdater_cfi::sort, strategyCfg_, remoteMonitoring_LED_IterMethod_cfg::threshold, and useGsfPfRecTracks_.
Referenced by produce().
◆ algo_
◆ cutsCfg_
◆ ecalSeedingParametersChecked_
bool GsfElectronProducer::ecalSeedingParametersChecked_ |
|
private |
◆ egmPFCandidateCollection_
◆ electronPutToken_
◆ gsfPfRecTracksTag_
◆ hcalCfg_
◆ hcalCfgBc_
◆ hcalRun2EffDepth_
bool GsfElectronProducer::hcalRun2EffDepth_ |
|
private |
◆ inputCfg_
◆ resetMvaValuesUsingPFCandidates_
const bool GsfElectronProducer::resetMvaValuesUsingPFCandidates_ |
|
private |
◆ strategyCfg_
◆ useGsfPfRecTracks_
const bool GsfElectronProducer::useGsfPfRecTracks_ |
|
private |
edm::EDGetTokenT< HBHERecHitCollection > hbheRecHits
double maxHOverEBarrelCone
ParameterDescriptionBase * add(U const &iLabel, T const &value)
std::vector< std::string > ecalRegressionWeightLabels
std::array< double, 4 > arrayHB
edm::EDGetTokenT< reco::GsfElectronCoreCollection > gsfElectronCores
bool useCombinationRegression
ElectronHcalHelper::Configuration hcalCfgBc_
T const * get() const
Returns C++ pointer to the item.
bool ignoreNotPreselected
GsfElectronAlgo::Tokens inputCfg_
U second(std::pair< T, U > const &p)
unsigned ambClustersOverlapStrategy
const GsfElectronAlgo::CutsConfiguration cutsCfg_
def create(alignables, pedeDump, additionalData, outputFile, config)
edm::EDGetTokenT< reco::ElectronSeedCollection > seedsTag
static constexpr float kMultThresEB
edm::EDGetTokenT< reco::BeamSpot > beamSpotTag
edm::EDGetTokenT< EcalRecHitCollection > barrelRecHitCollection
Log< level::Warning, false > LogWarning
bool isPreselected(reco::GsfElectron const &ele) const
bool ecalDrivenEcalErrorFromClassBasedParameterization
double maxHOverEEndcapsCone
const bool resetMvaValuesUsingPFCandidates_
unsigned ambSortingStrategy
ParameterSet const & parameterSet(StableProvenance const &provenance, ProcessHistory const &history)
bool pureTrackerDrivenEcalErrorFromSimpleParameterization
void add(std::string const &label, ParameterSetDescription const &psetDescription)
edm::EDGetTokenT< EcalRecHitCollection > endcapRecHitCollection
static constexpr float kMultThresEE
std::unique_ptr< GsfElectronAlgo > algo_
const bool useGsfPfRecTracks_
edm::EDGetTokenT< reco::ConversionCollection > conversions
edm::EDGetTokenT< reco::PFCandidateCollection > egmPFCandidateCollection_
static edm::ParameterSetDescription pSetDescript()
float sharedEnergy(reco::CaloCluster const &clu1, reco::CaloCluster const &clu2, EcalRecHitCollection const &barrelRecHits, EcalRecHitCollection const &endcapRecHits)
void checkEcalSeedingParameters(edm::ParameterSet const &)
const edm::EDPutTokenT< reco::GsfElectronCollection > electronPutToken_
EgammaHcalIsolation::arrayHB eThresHB
bool ecalSeedingParametersChecked_
edm::EDGetTokenT< reco::VertexCollection > vtxCollectionTag
ElectronHcalHelper::Configuration hcalCfg_
bool useDefaultEnergyCorrection
T getParameter(std::string const &) const
bool isBetterElectron(reco::GsfElectron const &, reco::GsfElectron const &)
bool isInnermostElectron(reco::GsfElectron const &, reco::GsfElectron const &)
GsfElectronAlgo::StrategyConfiguration strategyCfg_
edm::EDGetTokenT< reco::TrackCollection > ctfTracks
edm::EDGetTokenT< HBHERecHitCollection > hbheRecHitsTag
bool ecalDrivenEcalEnergyFromClassBasedParameterization
combinationRegressionWeightLabels
void setAmbiguityData(reco::GsfElectronCollection &electrons, edm::Event const &event, bool ignoreNotPreselected=true) const
EgammaHcalIsolation::arrayHE eThresHE
std::array< double, 7 > arrayHE
const edm::EDGetTokenT< reco::GsfPFRecTrackCollection > gsfPfRecTracksTag_