|
|
- Author
- Nancy Marinelli, U. of Notre Dame, US
Definition at line 49 of file ConversionTrackCandidateProducer.cc.
◆ ConversionTrackCandidateProducer()
ConversionTrackCandidateProducer::ConversionTrackCandidateProducer |
( |
const edm::ParameterSet & |
ps | ) |
|
Definition at line 132 of file ConversionTrackCandidateProducer.cc.
169 auto const& flagnamesEB =
config.getParameter<std::vector<std::string>>(
"RecHitFlagToBeExcludedEB");
170 auto const& flagnamesEE =
config.getParameter<std::vector<std::string>>(
"RecHitFlagToBeExcludedEE");
172 flagsexclEB_ = StringToEnumValue<EcalRecHit::Flags>(flagnamesEB);
173 flagsexclEE_ = StringToEnumValue<EcalRecHit::Flags>(flagnamesEE);
175 auto const& severitynamesEB =
config.getParameter<std::vector<std::string>>(
"RecHitSeverityToBeExcludedEB");
176 auto const& severitynamesEE =
config.getParameter<std::vector<std::string>>(
"RecHitSeverityToBeExcludedEE");
178 severitiesexclEB_ = StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEB);
179 severitiesexclEE_ = StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEE);
◆ beginRun()
◆ buildCollections()
Definition at line 329 of file ConversionTrackCandidateProducer.cc.
342 for (
auto const& aClus : scHandle->ptrs()) {
344 if (aClus->energy() / cosh(aClus->eta()) <=
minSCEt_)
346 if (aClus->eta() > 1.479 && aClus->eta() < 1.556)
351 double scEt = sc->
energy() / cosh(sc->
eta());
353 double HoE = towerIso.getTowerESum(sc) / sc->
energy();
393 for (
auto it = theOutInTracks.begin(); it != theOutInTracks.end(); ++it) {
394 vecRecOI.push_back(aClus);
398 for (
auto it = theInOutTracks.begin(); it != theInOutTracks.end(); ++it) {
399 vecRecIO.push_back(aClus);
References EgammaRecHitIsolation::doFlagChecks(), EgammaRecHitIsolation::doSeverityChecks(), DetId::Ecal, ecalIsoCut_offset_, ecalIsoCut_slope_, RecoTauCleanerPlugins::ecalIsolation, reco::CaloCluster::energy(), reco::CaloCluster::eta(), flagsexclEB_, flagsexclEE_, EgammaRecHitIsolation::getEtSum(), conversions_cfi::hcalTowers, hOverEConeSize_, inOutSeedFinder_, inOutTrackFinder_, PixelPluginsPhase0_cfi::isBarrel, isoConeR_, isoEMin_, isoEtaSlice_, isoEtMin_, isoInnerConeR_, OutInConversionSeedFinder::makeSeeds(), InOutConversionSeedFinder::makeSeeds(), maxHOverE_, minSCEt_, outInSeedFinder_, outInTrackFinder_, reco::CaloCluster::position(), ConversionSeedFinder::seeds(), ConversionSeedFinder::setCandidate(), InOutConversionSeedFinder::setTracks(), EgammaRecHitIsolation::setUseNumCrystals(), EgammaRecHitIsolation::setVetoClustered(), severitiesexclEB_, severitiesexclEE_, theCaloGeom_, InOutConversionTrackFinder::tracks(), OutInConversionTrackFinder::tracks(), useNumXtals_, and vetoClusteredHits_.
Referenced by produce().
◆ fillDescriptions()
Definition at line 405 of file ConversionTrackCandidateProducer.cc.
409 desc.add<
edm::InputTag>(
"bcBarrelCollection", {
"particleFlowSuperClusterECAL",
"particleFlowBasicClusterECALBarrel"});
410 desc.add<
edm::InputTag>(
"bcEndcapCollection", {
"particleFlowSuperClusterECAL",
"particleFlowBasicClusterECALEndcap"});
412 {
"particleFlowSuperClusterECAL",
"particleFlowSuperClusterECALBarrel"});
414 {
"particleFlowSuperClusterECAL",
"particleFlowSuperClusterECALEndcapWithPreshower"});
416 desc.add<
std::string>(
"outInTrackCandidateSCAssociationCollection",
"outInTrackCandidateSCAssociationCollection");
417 desc.add<
std::string>(
"inOutTrackCandidateSCAssociationCollection",
"inOutTrackCandidateSCAssociationCollection");
418 desc.add<
std::string>(
"outInTrackCandidateCollection",
"outInTracksFromConversions");
419 desc.add<
std::string>(
"inOutTrackCandidateCollection",
"inOutTracksFromConversions");
421 desc.add<
edm::InputTag>(
"barrelEcalRecHitCollection", {
"ecalRecHit",
"EcalRecHitsEB"});
422 desc.add<
edm::InputTag>(
"endcapEcalRecHitCollection", {
"ecalRecHit",
"EcalRecHitsEE"});
424 desc.add<
std::string>(
"OutInRedundantSeedCleaner",
"CachingSeedCleanerBySharedInput");
425 desc.add<
std::string>(
"InOutRedundantSeedCleaner",
"CachingSeedCleanerBySharedInput");
426 desc.add<
bool>(
"useHitsSplitting",
false);
427 desc.add<
int>(
"maxNumOfSeedsOutIn", 50);
428 desc.add<
int>(
"maxNumOfSeedsInOut", 50);
429 desc.add<
double>(
"bcEtCut", 1.5);
430 desc.add<
double>(
"bcECut", 1.5);
431 desc.add<
bool>(
"useEtCut",
true);
433 desc.add<
double>(
"minSCEt", 20.0);
434 desc.add<
double>(
"hOverEConeSize", 0.15);
435 desc.add<
double>(
"maxHOverE", 0.15);
436 desc.add<
double>(
"isoInnerConeR", 3.5);
437 desc.add<
double>(
"isoConeR", 0.4);
438 desc.add<
double>(
"isoEtaSlice", 2.5);
439 desc.add<
double>(
"isoEtMin", 0.0);
440 desc.add<
double>(
"isoEMin", 0.08);
441 desc.add<
bool>(
"vetoClusteredHits",
false);
442 desc.add<
bool>(
"useNumXstals",
true);
443 desc.add<
double>(
"ecalIsoCut_offset", 999999999);
444 desc.add<
double>(
"ecalIsoCut_slope", 0.0);
446 desc.add<std::vector<std::string>>(
"RecHitFlagToBeExcludedEB", {});
447 desc.add<std::vector<std::string>>(
"RecHitSeverityToBeExcludedEB", {});
448 desc.add<std::vector<std::string>>(
"RecHitFlagToBeExcludedEE", {});
449 desc.add<std::vector<std::string>>(
"RecHitSeverityToBeExcludedEE", {});
451 desc.add<
double>(
"fractionShared", 0.5);
452 desc.add<
std::string>(
"TrajectoryBuilder",
"TrajectoryBuilderForConversions");
460 psd0.
add<
std::string>(
"propagatorAlongTISE",
"alongMomElePropagator");
461 psd0.
add<
int>(
"numberMeasurementsForFit", 4);
462 psd0.
add<
std::string>(
"propagatorOppositeTISE",
"oppositeToMomElePropagator");
465 desc.add<
bool>(
"allowSharedFirstHit",
true);
466 desc.add<
double>(
"ValidHitBonus", 5.0);
467 desc.add<
double>(
"MissingHitPenalty", 20.0);
469 descriptions.
add(
"conversionTrackCandidates",
desc);
References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, edm::ParameterSetDescription::setUnknown(), and AlCaHLTBitMon_QueryRunRegistry::string.
◆ produce()
Definition at line 205 of file ConversionTrackCandidateProducer.cc.
219 auto outInTrackCandidate_p = std::make_unique<TrackCandidateCollection>();
221 auto inOutTrackCandidate_p = std::make_unique<TrackCandidateCollection>();
223 auto outInAssoc_p = std::make_unique<reco::TrackCandidateCaloClusterPtrAssociation>();
224 auto inOutAssoc_p = std::make_unique<reco::TrackCandidateCaloClusterPtrAssociation>();
227 bool validBarrelBCHandle =
true;
229 if (!bcBarrelHandle.isValid()) {
230 edm::LogError(
"ConversionTrackCandidateProducer") <<
"Error! Can't get the Barrel Basic Clusters!";
231 validBarrelBCHandle =
false;
235 bool validEndcapBCHandle =
true;
237 if (!bcEndcapHandle.isValid()) {
238 edm::LogError(
"CoonversionTrackCandidateProducer") <<
"Error! Can't get the Endcap Basic Clusters";
239 validEndcapBCHandle =
false;
243 bool validBarrelSCHandle =
true;
245 if (!scBarrelHandle.isValid()) {
246 edm::LogError(
"CoonversionTrackCandidateProducer") <<
"Error! Can't get the barrel superclusters!";
247 validBarrelSCHandle =
false;
251 bool validEndcapSCHandle =
true;
253 if (!scEndcapHandle.isValid()) {
254 edm::LogError(
"CoonversionTrackCandidateProducer") <<
"Error! Can't get the endcap superclusters!";
255 validEndcapSCHandle =
false;
274 if (validBarrelBCHandle && validBarrelSCHandle)
281 *outInTrackCandidate_p,
282 *inOutTrackCandidate_p,
286 if (validEndcapBCHandle && validEndcapSCHandle) {
294 *outInTrackCandidate_p,
295 *inOutTrackCandidate_p,
References barrelecalCollection_, bcBarrelCollection_, bcEndcapCollection_, buildCollections(), caloPtrVecInOut_, caloPtrVecOutIn_, ConversionSeedFinder::clear(), endcapecalCollection_, edm::helper::Filler< Map >::fill(), edm::EventSetup::get(), edm::Event::get(), get, edm::Event::getHandle(), conversions_cfi::hcalTowers, hcalTowers_, inOutSeedFinder_, InOutTrackCandidateCollection_, InOutTrackSCAssociationCollection_, edm::helper::Filler< Map >::insert(), PixelPluginsPhase0_cfi::isBarrel, measurementTrkEvtToken_, eostools::move(), outInSeedFinder_, OutInTrackCandidateCollection_, OutInTrackSCAssociationCollection_, edm::ESHandle< T >::product(), edm::Event::put(), scHybridBarrelProducer_, scIslandEndcapProducer_, ConversionSeedFinder::setEvent(), setEventSetup(), theCaloGeom_, and theTrajectoryBuilder_.
◆ setEventSetup()
void ConversionTrackCandidateProducer::setEventSetup |
( |
const edm::EventSetup & |
es | ) |
|
|
private |
◆ barrelecalCollection_
◆ bcBarrelCollection_
◆ bcEndcapCollection_
◆ caloPtrVecInOut_
◆ caloPtrVecOutIn_
◆ ecalIsoCut_offset_
double ConversionTrackCandidateProducer::ecalIsoCut_offset_ |
|
private |
◆ ecalIsoCut_slope_
double ConversionTrackCandidateProducer::ecalIsoCut_slope_ |
|
private |
◆ endcapecalCollection_
◆ flagsexclEB_
std::vector<int> ConversionTrackCandidateProducer::flagsexclEB_ |
|
private |
◆ flagsexclEE_
std::vector<int> ConversionTrackCandidateProducer::flagsexclEE_ |
|
private |
◆ hcalTowers_
◆ hOverEConeSize_
double ConversionTrackCandidateProducer::hOverEConeSize_ |
|
private |
◆ inOutSeedFinder_
◆ InOutTrackCandidateCollection_
std::string ConversionTrackCandidateProducer::InOutTrackCandidateCollection_ |
|
private |
◆ inOutTrackFinder_
◆ InOutTrackSCAssociationCollection_
std::string ConversionTrackCandidateProducer::InOutTrackSCAssociationCollection_ |
|
private |
◆ isoConeR_
double ConversionTrackCandidateProducer::isoConeR_ |
|
private |
◆ isoEMin_
double ConversionTrackCandidateProducer::isoEMin_ |
|
private |
◆ isoEtaSlice_
double ConversionTrackCandidateProducer::isoEtaSlice_ |
|
private |
◆ isoEtMin_
double ConversionTrackCandidateProducer::isoEtMin_ |
|
private |
◆ isoInnerConeR_
double ConversionTrackCandidateProducer::isoInnerConeR_ |
|
private |
◆ maxHOverE_
double ConversionTrackCandidateProducer::maxHOverE_ |
|
private |
◆ measurementTrkEvtToken_
◆ minSCEt_
double ConversionTrackCandidateProducer::minSCEt_ |
|
private |
◆ outInSeedFinder_
◆ OutInTrackCandidateCollection_
std::string ConversionTrackCandidateProducer::OutInTrackCandidateCollection_ |
|
private |
◆ outInTrackFinder_
◆ OutInTrackSCAssociationCollection_
std::string ConversionTrackCandidateProducer::OutInTrackSCAssociationCollection_ |
|
private |
◆ scHybridBarrelProducer_
◆ scIslandEndcapProducer_
◆ severitiesexclEB_
std::vector<int> ConversionTrackCandidateProducer::severitiesexclEB_ |
|
private |
◆ severitiesexclEE_
std::vector<int> ConversionTrackCandidateProducer::severitiesexclEE_ |
|
private |
◆ theCaloGeom_
◆ theTrajectoryBuilder_
◆ useNumXtals_
bool ConversionTrackCandidateProducer::useNumXtals_ |
|
private |
◆ vecOfSCRefForInOut
◆ vecOfSCRefForOutIn
◆ vetoClusteredHits_
bool ConversionTrackCandidateProducer::vetoClusteredHits_ |
|
private |
OutInConversionSeedFinder outInSeedFinder_
T const * product() const
void setEventSetup(const edm::EventSetup &es)
Initialize EventSetup objects at each event.
std::vector< Trajectory > tracks(const TrajectorySeedCollection &seeds, TrackCandidateCollection &candidate) const override
void setTracks(std::vector< Trajectory > const &in)
edm::EDGetTokenT< edm::View< reco::CaloCluster > > scIslandEndcapProducer_
InOutConversionTrackFinder inOutTrackFinder_
std::vector< int > severitiesexclEB_
edm::EDGetTokenT< edm::View< reco::CaloCluster > > bcBarrelCollection_
std::unique_ptr< BaseCkfTrajectoryBuilder > theTrajectoryBuilder_
ParameterDescriptionBase * add(U const &iLabel, T const &value)
InOutConversionSeedFinder inOutSeedFinder_
std::string OutInTrackSCAssociationCollection_
std::string OutInTrackCandidateCollection_
std::vector< int > severitiesexclEE_
edm::EDGetTokenT< EcalRecHitCollection > barrelecalCollection_
edm::EDGetTokenT< EcalRecHitCollection > endcapecalCollection_
std::vector< Trajectory > tracks(const TrajectorySeedCollection &seeds, TrackCandidateCollection &candidates) const override
std::vector< int > flagsexclEE_
void add(std::string const &label, ParameterSetDescription const &psetDescription)
edm::EDGetTokenT< CaloTowerCollection > hcalTowers_
config
parse the configuration file
std::vector< edm::Ptr< reco::CaloCluster > > caloPtrVecOutIn_
Global3DPoint GlobalPoint
void setEventSetup(const edm::EventSetup &es)
Initialize EventSetup objects at each event.
void setEventSetup(const edm::EventSetup &es)
Initialize EventSetup objects at each event.
std::string InOutTrackCandidateCollection_
std::vector< edm::Ptr< reco::CaloCluster > > caloPtrVecInOut_
void makeSeeds(const edm::Handle< edm::View< reco::CaloCluster > > &allBc) override
void setNavigationSchool(const NavigationSchool *navigation)
std::string InOutTrackSCAssociationCollection_
double ecalIsoCut_offset_
edm::EDGetTokenT< MeasurementTrackerEvent > measurementTrkEvtToken_
double eta() const
pseudorapidity of cluster centroid
OutInConversionTrackFinder outInTrackFinder_
TrajectorySeedCollection & seeds()
virtual void setCandidate(float e, GlobalPoint pos)
Log< level::Error, false > LogError
void makeSeeds(const edm::Handle< edm::View< reco::CaloCluster > > &allBc) override
const math::XYZPoint & position() const
cluster centroid position
void setEvent(const edm::Event &e)
edm::EDGetTokenT< edm::View< reco::CaloCluster > > scHybridBarrelProducer_
edm::ESHandle< CaloGeometry > theCaloGeom_
edm::EDGetTokenT< edm::View< reco::CaloCluster > > bcEndcapCollection_
double energy() const
cluster energy
void buildCollections(bool detector, const edm::Handle< edm::View< reco::CaloCluster >> &scHandle, const edm::Handle< edm::View< reco::CaloCluster >> &bcHandle, const EcalRecHitCollection &ecalRecHits, const EcalSeverityLevelAlgo *sevLev, CaloTowerCollection const &hcalTowers, TrackCandidateCollection &outInTracks, TrackCandidateCollection &inOutTracks, std::vector< edm::Ptr< reco::CaloCluster >> &vecRecOI, std::vector< edm::Ptr< reco::CaloCluster >> &vecRecIO)
std::vector< int > flagsexclEB_