CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Types | Private Member Functions | Private Attributes
pat::PATPhotonProducer Class Reference

Produces the pat::Photon. More...

Inheritance diagram for pat::PATPhotonProducer:
edm::stream::EDProducer<>

Public Member Functions

 PATPhotonProducer (const edm::ParameterSet &iConfig)
 
void produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
 
 ~PATPhotonProducer () override
 
- 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 EDProduceroperator= (const EDProducer &)=delete
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 

Private Types

typedef std::vector< edm::Handle< edm::ValueMap< IsoDeposit > > > IsoDepositMaps
 
typedef std::pair< pat::IsolationKeys, edm::InputTagIsolationLabel
 
typedef std::vector< IsolationLabelIsolationLabels
 
typedef std::vector< edm::Handle< edm::ValueMap< double > > > IsolationValueMaps
 
typedef std::pair< std::string, edm::InputTagNameTag
 

Private Member Functions

template<typename T >
void readIsolationLabels (const edm::ParameterSet &iConfig, const char *psetName, IsolationLabels &labels, std::vector< edm::EDGetTokenT< edm::ValueMap< T >>> &tokens)
 

Private Attributes

bool addEfficiencies_
 
bool addGenMatch_
 
bool addPFClusterIso_
 
bool addPhotonID_
 
bool addPuppiIsolation_
 
bool addResolutions_
 
edm::EDGetTokenT< reco::BeamSpotbeamLineToken_
 
const EcalClusterLazyTools::ESGetTokens ecalClusterToolsESGetTokens_
 
const CaloGeometryecalGeometry_
 
const edm::ESGetToken< CaloGeometry, CaloGeometryRecordecalGeometryToken_
 
edm::EDGetTokenT< edm::ValueMap< float > > ecalPFClusterIsoT_
 
const CaloTopologyecalTopology_
 
const edm::ESGetToken< CaloTopology, CaloTopologyRecordecalTopologyToken_
 
pat::helper::EfficiencyLoader efficiencyLoader_
 
edm::EDGetTokenT< reco::GsfElectronCollectionelectronToken_
 
bool embedBasicClusters_
 
bool embedGenMatch_
 
bool embedPreshowerClusters_
 
bool embedRecHits_
 
bool embedSeedCluster_
 
bool embedSuperCluster_
 
GreaterByEt< PhotoneTComparator_
 
std::vector< edm::EDGetTokenT< edm::Association< reco::GenParticleCollection > > > genMatchTokens_
 
edm::EDGetTokenT< edm::ValueMap< float > > hcalPFClusterIsoT_
 
edm::EDGetTokenT< reco::ConversionCollectionhConversionsToken_
 
IsolationLabels isoDepositLabels_
 
std::vector< edm::EDGetTokenT< edm::ValueMap< IsoDeposit > > > isoDepositTokens_
 
IsolationLabels isolationValueLabels_
 
std::vector< edm::EDGetTokenT< edm::ValueMap< double > > > isolationValueTokens_
 
pat::helper::MultiIsolator isolator_
 
pat::helper::MultiIsolator::IsolationValuePairs isolatorTmpStorage_
 
std::vector< NameTagphotIDSrcs_
 
std::vector< edm::EDGetTokenT< edm::ValueMap< Bool_t > > > photIDTokens_
 
edm::EDGetTokenT< edm::View< reco::Photon > > photonToken_
 
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_charged_hadrons_
 
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_neutral_hadrons_
 
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_photons_
 
edm::InputTag reducedBarrelRecHitCollection_
 
edm::EDGetTokenT< EcalRecHitCollectionreducedBarrelRecHitCollectionToken_
 
edm::InputTag reducedEndcapRecHitCollection_
 
edm::EDGetTokenT< EcalRecHitCollectionreducedEndcapRecHitCollectionToken_
 
pat::helper::KinResolutionsLoader resolutionLoader_
 
bool saveRegressionData_
 
pat::PATUserDataHelper< pat::PhotonuserDataHelper_
 
bool useUserData_
 

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

Produces the pat::Photon.

The PATPhotonProducer produces the analysis-level pat::Photon starting from a collection of objects of PhotonType.

Author
Steven Lowette
Version
Id
PATPhotonProducer.h,v 1.19 2009/06/25 23:49:35 gpetrucc Exp

Definition at line 49 of file PATPhotonProducer.cc.

Member Typedef Documentation

◆ IsoDepositMaps

Definition at line 90 of file PATPhotonProducer.cc.

◆ IsolationLabel

Definition at line 92 of file PATPhotonProducer.cc.

◆ IsolationLabels

Definition at line 93 of file PATPhotonProducer.cc.

◆ IsolationValueMaps

typedef std::vector<edm::Handle<edm::ValueMap<double> > > pat::PATPhotonProducer::IsolationValueMaps
private

Definition at line 91 of file PATPhotonProducer.cc.

◆ NameTag

typedef std::pair<std::string, edm::InputTag> pat::PATPhotonProducer::NameTag
private

Definition at line 118 of file PATPhotonProducer.cc.

Constructor & Destructor Documentation

◆ PATPhotonProducer()

PATPhotonProducer::PATPhotonProducer ( const edm::ParameterSet iConfig)
explicit

Definition at line 195 of file PATPhotonProducer.cc.

196  :
197 
198  ecalClusterToolsESGetTokens_{consumesCollector()},
199  isolator_(iConfig.exists("userIsolation") ? iConfig.getParameter<edm::ParameterSet>("userIsolation")
200  : edm::ParameterSet(),
201  consumesCollector(),
202  false),
203  useUserData_(iConfig.exists("userData")),
206  // initialize the configurables
207  photonToken_ = consumes<edm::View<reco::Photon>>(iConfig.getParameter<edm::InputTag>("photonSource"));
208  electronToken_ = consumes<reco::GsfElectronCollection>(iConfig.getParameter<edm::InputTag>("electronSource"));
209  hConversionsToken_ = consumes<reco::ConversionCollection>(iConfig.getParameter<edm::InputTag>("conversionSource"));
210  beamLineToken_ = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamLineSrc"));
211  embedSuperCluster_ = iConfig.getParameter<bool>("embedSuperCluster");
212  embedSeedCluster_ = iConfig.getParameter<bool>("embedSeedCluster");
213  embedBasicClusters_ = iConfig.getParameter<bool>("embedBasicClusters");
214  embedPreshowerClusters_ = iConfig.getParameter<bool>("embedPreshowerClusters");
215  embedRecHits_ = iConfig.getParameter<bool>("embedRecHits");
216  reducedBarrelRecHitCollection_ = iConfig.getParameter<edm::InputTag>("reducedBarrelRecHitCollection");
217  reducedBarrelRecHitCollectionToken_ = mayConsume<EcalRecHitCollection>(reducedBarrelRecHitCollection_);
218  reducedEndcapRecHitCollection_ = iConfig.getParameter<edm::InputTag>("reducedEndcapRecHitCollection");
219  reducedEndcapRecHitCollectionToken_ = mayConsume<EcalRecHitCollection>(reducedEndcapRecHitCollection_);
220  // MC matching configurables
221  addGenMatch_ = iConfig.getParameter<bool>("addGenMatch");
222  if (addGenMatch_) {
223  embedGenMatch_ = iConfig.getParameter<bool>("embedGenMatch");
224  if (iConfig.existsAs<edm::InputTag>("genParticleMatch")) {
226  iConfig.getParameter<edm::InputTag>("genParticleMatch")));
227  } else {
229  iConfig.getParameter<std::vector<edm::InputTag>>("genParticleMatch"),
230  [this](edm::InputTag const& tag) { return consumes<edm::Association<reco::GenParticleCollection>>(tag); });
231  }
232  }
233  // Efficiency configurables
234  addEfficiencies_ = iConfig.getParameter<bool>("addEfficiencies");
235  if (addEfficiencies_) {
237  pat::helper::EfficiencyLoader(iConfig.getParameter<edm::ParameterSet>("efficiencies"), consumesCollector());
238  }
239  // PFCluster Isolation maps
240  addPuppiIsolation_ = iConfig.getParameter<bool>("addPuppiIsolation");
241  if (addPuppiIsolation_) {
243  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiIsolationChargedHadrons"));
245  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiIsolationNeutralHadrons"));
247  consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("puppiIsolationPhotons"));
248  }
249  addPFClusterIso_ = iConfig.getParameter<bool>("addPFClusterIso");
250  if (addPFClusterIso_) {
251  ecalPFClusterIsoT_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("ecalPFClusterIsoMap"));
252  auto hcPFC = iConfig.getParameter<edm::InputTag>("hcalPFClusterIsoMap");
253  if (not hcPFC.label().empty())
254  hcalPFClusterIsoT_ = consumes<edm::ValueMap<float>>(hcPFC);
255  }
256 
257  // photon ID configurables
258  addPhotonID_ = iConfig.getParameter<bool>("addPhotonID");
259  if (addPhotonID_) {
260  // it might be a single photon ID
261  if (iConfig.existsAs<edm::InputTag>("photonIDSource")) {
262  photIDSrcs_.push_back(NameTag("", iConfig.getParameter<edm::InputTag>("photonIDSource")));
263  }
264  // or there might be many of them
265  if (iConfig.existsAs<edm::ParameterSet>("photonIDSources")) {
266  // please don't configure me twice
267  if (!photIDSrcs_.empty()) {
268  throw cms::Exception("Configuration")
269  << "PATPhotonProducer: you can't specify both 'photonIDSource' and 'photonIDSources'\n";
270  }
271  // read the different photon ID names
272  edm::ParameterSet idps = iConfig.getParameter<edm::ParameterSet>("photonIDSources");
273  std::vector<std::string> names = idps.getParameterNamesForType<edm::InputTag>();
274  for (std::vector<std::string>::const_iterator it = names.begin(), ed = names.end(); it != ed; ++it) {
275  photIDSrcs_.push_back(NameTag(*it, idps.getParameter<edm::InputTag>(*it)));
276  }
277  }
278  // but in any case at least once
279  if (photIDSrcs_.empty())
280  throw cms::Exception("Configuration") << "PATPhotonProducer: id addPhotonID is true, you must specify either:\n"
281  << "\tInputTag photonIDSource = <someTag>\n"
282  << "or\n"
283  << "\tPSet photonIDSources = { \n"
284  << "\t\tInputTag <someName> = <someTag> // as many as you want \n "
285  << "\t}\n";
286  }
288  photIDSrcs_, [this](NameTag const& tag) { return mayConsume<edm::ValueMap<Bool_t>>(tag.second); });
289  // Resolution configurables
290  addResolutions_ = iConfig.getParameter<bool>("addResolutions");
291  if (addResolutions_) {
293  pat::helper::KinResolutionsLoader(iConfig.getParameter<edm::ParameterSet>("resolutions"), consumesCollector());
294  }
295  // Check to see if the user wants to add user data
296  if (useUserData_) {
298  PATUserDataHelper<Photon>(iConfig.getParameter<edm::ParameterSet>("userData"), consumesCollector());
299  }
300  // produces vector of photons
301  produces<std::vector<Photon>>();
302 
303  // read isoDeposit labels, for direct embedding
304  readIsolationLabels(iConfig, "isoDeposits", isoDepositLabels_, isoDepositTokens_);
305  // read isolation value labels, for direct embedding
307 
308  saveRegressionData_ = iConfig.getParameter<bool>("saveRegressionData");
309 }
Assists in assimilating all pat::UserData into pat objects.
ESGetTokenH3DDVariant esConsumes(std::string const &Record, edm::ConsumesCollector &)
Definition: DeDxTools.cc:283
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
pat::PATUserDataHelper< pat::Photon > userDataHelper_
edm::EDGetTokenT< EcalRecHitCollection > reducedEndcapRecHitCollectionToken_
std::pair< std::string, edm::InputTag > NameTag
void readIsolationLabels(const edm::ParameterSet &iConfig, const char *psetName, IsolationLabels &labels, std::vector< edm::EDGetTokenT< edm::ValueMap< T >>> &tokens)
std::vector< edm::EDGetTokenT< edm::ValueMap< Bool_t > > > photIDTokens_
bool exists(std::string const &parameterName) const
checks if a parameter exists
edm::EDGetTokenT< reco::BeamSpot > beamLineToken_
edm::EDGetTokenT< edm::View< reco::Photon > > photonToken_
edm::EDGetTokenT< edm::ValueMap< float > > ecalPFClusterIsoT_
std::vector< edm::EDGetTokenT< edm::Association< reco::GenParticleCollection > > > genMatchTokens_
pat::helper::KinResolutionsLoader resolutionLoader_
edm::EDGetTokenT< edm::ValueMap< float > > hcalPFClusterIsoT_
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:171
auto vector_transform(std::vector< InputType > const &input, Function predicate) -> std::vector< typename std::remove_cv< typename std::remove_reference< decltype(predicate(input.front()))>::type >::type >
Definition: transform.h:11
const std::string names[nVars_]
edm::InputTag reducedEndcapRecHitCollection_
IsolationLabels isoDepositLabels_
IsolationLabels isolationValueLabels_
std::vector< edm::EDGetTokenT< edm::ValueMap< IsoDeposit > > > isoDepositTokens_
pat::helper::EfficiencyLoader efficiencyLoader_
std::vector< std::string > getParameterNamesForType(bool trackiness=true) const
Definition: ParameterSet.h:179
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_photons_
std::vector< NameTag > photIDSrcs_
pat::helper::MultiIsolator isolator_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_charged_hadrons_
const EcalClusterLazyTools::ESGetTokens ecalClusterToolsESGetTokens_
const edm::ESGetToken< CaloTopology, CaloTopologyRecord > ecalTopologyToken_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_neutral_hadrons_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord > ecalGeometryToken_
edm::EDGetTokenT< EcalRecHitCollection > reducedBarrelRecHitCollectionToken_
std::vector< edm::EDGetTokenT< edm::ValueMap< double > > > isolationValueTokens_
edm::InputTag reducedBarrelRecHitCollection_
edm::EDGetTokenT< reco::ConversionCollection > hConversionsToken_
edm::EDGetTokenT< reco::GsfElectronCollection > electronToken_

◆ ~PATPhotonProducer()

PATPhotonProducer::~PATPhotonProducer ( )
override

Definition at line 311 of file PATPhotonProducer.cc.

311 {}

Member Function Documentation

◆ fillDescriptions()

void PATPhotonProducer::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 642 of file PATPhotonProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addNode(), edm::ParameterSetDescription::addOptional(), pat::helper::KinResolutionsLoader::fillDescription(), pat::PATUserDataHelper< ObjectType >::fillDescription(), edm::ParameterSetDescription::ifValue(), HLT_2022v12_cff::InputTag, or, edm::ParameterSetDescription::setAllowAnything(), edm::ParameterSetDescription::setComment(), and edm::ParameterDescriptionNode::setComment().

642  {
644  iDesc.setComment("PAT photon producer module");
645 
646  // input source
647  iDesc.add<edm::InputTag>("photonSource", edm::InputTag("no default"))->setComment("input collection");
648  iDesc.add<edm::InputTag>("electronSource", edm::InputTag("no default"))->setComment("input collection");
649  iDesc.add<edm::InputTag>("conversionSource", edm::InputTag("allConversions"))->setComment("input collection");
650 
651  iDesc.add<edm::InputTag>("reducedBarrelRecHitCollection", edm::InputTag("reducedEcalRecHitsEB"));
652  iDesc.add<edm::InputTag>("reducedEndcapRecHitCollection", edm::InputTag("reducedEcalRecHitsEE"));
653 
654  iDesc.ifValue(
655  edm::ParameterDescription<bool>("addPFClusterIso", false, true),
657  "ecalPFClusterIsoMap", edm::InputTag("photonEcalPFClusterIsolationProducer"), true) and
659  "hcalPFClusterIsoMap", edm::InputTag("photonHcalPFClusterIsolationProducer"), true)) or
660  false >> (edm::ParameterDescription<edm::InputTag>("ecalPFClusterIsoMap", edm::InputTag(""), true) and
661  edm::ParameterDescription<edm::InputTag>("hcalPFClusterIsoMap", edm::InputTag(""), true)));
662 
663  iDesc.ifValue(
664  edm::ParameterDescription<bool>("addPuppiIsolation", false, true),
666  "puppiIsolationChargedHadrons", edm::InputTag("egmPhotonPUPPIIsolation", "h+-DR030-"), true) and
668  "puppiIsolationNeutralHadrons", edm::InputTag("egmPhotonPUPPIIsolation", "h0-DR030-"), true) and
670  "puppiIsolationPhotons", edm::InputTag("egmPhotonPUPPIIsolation", "gamma-DR030-"), true)) or
671  false >> edm::EmptyGroupDescription());
672 
673  iDesc.add<bool>("embedSuperCluster", true)->setComment("embed external super cluster");
674  iDesc.add<bool>("embedSeedCluster", true)->setComment("embed external seed cluster");
675  iDesc.add<bool>("embedBasicClusters", true)->setComment("embed external basic clusters");
676  iDesc.add<bool>("embedPreshowerClusters", true)->setComment("embed external preshower clusters");
677  iDesc.add<bool>("embedRecHits", true)->setComment("embed external RecHits");
678 
679  // MC matching configurables
680  iDesc.add<bool>("addGenMatch", true)->setComment("add MC matching");
681  iDesc.add<bool>("embedGenMatch", false)->setComment("embed MC matched MC information");
682  std::vector<edm::InputTag> emptySourceVector;
683  iDesc
684  .addNode(edm::ParameterDescription<edm::InputTag>("genParticleMatch", edm::InputTag(), true) xor
685  edm::ParameterDescription<std::vector<edm::InputTag>>("genParticleMatch", emptySourceVector, true))
686  ->setComment("input with MC match information");
687 
689 
690  // photon ID configurables
691  iDesc.add<bool>("addPhotonID", true)->setComment("add photon ID variables");
692  edm::ParameterSetDescription photonIDSourcesPSet;
693  photonIDSourcesPSet.setAllowAnything();
694  iDesc
695  .addNode(edm::ParameterDescription<edm::InputTag>("photonIDSource", edm::InputTag(), true) xor
696  edm::ParameterDescription<edm::ParameterSetDescription>("photonIDSources", photonIDSourcesPSet, true))
697  ->setComment("input with photon ID variables");
698 
699  // IsoDeposit configurables
700  edm::ParameterSetDescription isoDepositsPSet;
701  isoDepositsPSet.addOptional<edm::InputTag>("tracker");
702  isoDepositsPSet.addOptional<edm::InputTag>("ecal");
703  isoDepositsPSet.addOptional<edm::InputTag>("hcal");
704  isoDepositsPSet.addOptional<edm::InputTag>("pfAllParticles");
705  isoDepositsPSet.addOptional<edm::InputTag>("pfChargedHadrons");
706  isoDepositsPSet.addOptional<edm::InputTag>("pfChargedAll");
707  isoDepositsPSet.addOptional<edm::InputTag>("pfPUChargedHadrons");
708  isoDepositsPSet.addOptional<edm::InputTag>("pfNeutralHadrons");
709  isoDepositsPSet.addOptional<edm::InputTag>("pfPhotons");
710  isoDepositsPSet.addOptional<std::vector<edm::InputTag>>("user");
711  iDesc.addOptional("isoDeposits", isoDepositsPSet);
712 
713  // isolation values configurables
714  edm::ParameterSetDescription isolationValuesPSet;
715  isolationValuesPSet.addOptional<edm::InputTag>("tracker");
716  isolationValuesPSet.addOptional<edm::InputTag>("ecal");
717  isolationValuesPSet.addOptional<edm::InputTag>("hcal");
718  isolationValuesPSet.addOptional<edm::InputTag>("pfAllParticles");
719  isolationValuesPSet.addOptional<edm::InputTag>("pfChargedHadrons");
720  isolationValuesPSet.addOptional<edm::InputTag>("pfChargedAll");
721  isolationValuesPSet.addOptional<edm::InputTag>("pfPUChargedHadrons");
722  isolationValuesPSet.addOptional<edm::InputTag>("pfNeutralHadrons");
723  isolationValuesPSet.addOptional<edm::InputTag>("pfPhotons");
724  isolationValuesPSet.addOptional<std::vector<edm::InputTag>>("user");
725  iDesc.addOptional("isolationValues", isolationValuesPSet);
726 
727  // Efficiency configurables
728  edm::ParameterSetDescription efficienciesPSet;
729  efficienciesPSet.setAllowAnything(); // TODO: the pat helper needs to implement a description.
730  iDesc.add("efficiencies", efficienciesPSet);
731  iDesc.add<bool>("addEfficiencies", false);
732 
733  // Check to see if the user wants to add user data
734  edm::ParameterSetDescription userDataPSet;
736  iDesc.addOptional("userData", userDataPSet);
737 
738  edm::ParameterSetDescription isolationPSet;
739  isolationPSet.setAllowAnything(); // TODO: the pat helper needs to implement a description.
740  iDesc.add("userIsolation", isolationPSet);
741 
742  iDesc.addNode(edm::ParameterDescription<edm::InputTag>("beamLineSrc", edm::InputTag(), true))
743  ->setComment("input with high level selection");
744 
745  iDesc.add<bool>("saveRegressionData", true)->setComment("save regression input variables");
746 
747  descriptions.add("PATPhotonProducer", iDesc);
748 }
void setComment(std::string const &value)
ParameterDescriptionNode * ifValue(ParameterDescription< T > const &switchParameter, std::unique_ptr< ParameterDescriptionCases< T >> cases)
ParameterDescriptionBase * addOptional(U const &iLabel, T const &value)
void setAllowAnything()
allow any parameter label/value pairs
ParameterDescriptionNode * addNode(ParameterDescriptionNode const &node)
static void fillDescription(edm::ParameterSetDescription &iDesc)
void setComment(std::string const &value)
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
ParameterDescriptionBase * add(U const &iLabel, T const &value)
static void fillDescription(edm::ParameterSetDescription &iDesc)
Method for documentation and validation of PSet.
void add(std::string const &label, ParameterSetDescription const &psetDescription)

◆ produce()

void PATPhotonProducer::produce ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
override

Definition at line 313 of file PATPhotonProducer.cc.

References pat::PATUserDataHelper< ObjectType >::add(), addGenMatch_, pat::PATObject< ObjectType >::addGenParticleRef(), addPFClusterIso_, addPhotonID_, addPuppiIsolation_, Reference_intrackfit_cff::barrel, beamLineToken_, pat::helper::MultiIsolator::beginEvent(), EcalRegressionData::C1, EcalRegressionData::C2, EcalRegressionData::C3, CandIsolatorFromDeposits_cfi::deposits, EcalRegressionData::e2nd(), EcalRegressionData::e3x3(), EcalRegressionData::eBottom(), DetId::Ecal, EcalBarrel, ecalClusterToolsESGetTokens_, EcalEndcap, ecalGeometry_, ecalGeometryToken_, ecalPFClusterIsoT_, ecalTopology_, ecalTopologyToken_, efficiencyLoader_, electronToken_, EcalRegressionData::eLeft(), EcalRegressionData::eMax(), pat::Photon::embedBasicClusters(), embedBasicClusters_, embedGenMatch_, pat::PATObject< ObjectType >::embedGenParticle(), pat::Photon::embedPreshowerClusters(), embedPreshowerClusters_, pat::Photon::embedRecHits(), embedRecHits_, pat::Photon::embedSeedCluster(), embedSeedCluster_, pat::Photon::embedSuperCluster(), embedSuperCluster_, pat::helper::EfficiencyLoader::enabled(), pat::helper::KinResolutionsLoader::enabled(), pat::helper::MultiIsolator::enabled(), pat::helper::MultiIsolator::endEvent(), EcalRegressionData::eRight(), eTComparator_, EcalRegressionData::eTop(), pat::helper::MultiIsolator::fill(), EcalRegressionData::fill(), first, genMatchTokens_, EcalClusterLazyToolsBase::ESGetTokens::get(), edm::EventSetup::getData(), reco::Photon::getPflowIsolationVariables(), CaloTopology::getSubdetectorTopology(), CaloSubdetectorTopology::getWindow(), ConversionTools::hasMatchedPromptElectron(), reco::Photon::hasPixelSeed(), hcalPFClusterIsoT_, hConversionsToken_, mps_fire::i, heavyIonCSV_trainingSettings::idx, iEvent, isoDepositLabels_, isoDepositTokens_, isolationValueLabels_, displacedMuonProducer_cff::isolationValues, isolationValueTokens_, isolator_, isolatorTmpStorage_, edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), dqmiolumiharvest::j, EcalRegressionData::maxSubClusDR(), EcalRegressionData::maxSubClusDRDEta(), EcalRegressionData::maxSubClusDRDPhi(), EcalRegressionData::maxSubClusDRRawEnergy(), eostools::move(), dqmiodumpmetadata::n, pat::helper::EfficiencyLoader::newEvent(), pat::helper::KinResolutionsLoader::newEvent(), photIDSrcs_, photIDTokens_, BPHMonitor_cfi::photons, photonToken_, reco::BeamSpot::position(), edm::Handle< T >::product(), PUPPIIsolation_charged_hadrons_, PUPPIIsolation_neutral_hadrons_, PUPPIIsolation_photons_, edm::SortedCollection< T, SORT >::push_back(), FastTrackerRecHitMaskProducer_cfi::recHits, reducedBarrelRecHitCollectionToken_, reducedEndcapRecHitCollectionToken_, resolutionLoader_, saveRegressionData_, fileCollector::seed, EcalRegressionData::seedCrysEtaOrX(), EcalRegressionData::seedCrysIEtaOrIX(), EcalRegressionData::seedCrysIPhiOrIY(), EcalRegressionData::seedCrysPhiOrY(), pat::Photon::setCryEta(), pat::Photon::setCryPhi(), pat::Photon::setE2nd(), pat::Photon::setE3x3(), pat::Photon::setEBottom(), pat::helper::EfficiencyLoader::setEfficiencies(), pat::Photon::setELeft(), pat::Photon::setEMax(), pat::Photon::setERight(), pat::Photon::setETop(), pat::Photon::setHasPixelSeed(), pat::Photon::setIEta(), pat::Photon::setIPhi(), pat::Photon::setIsoDeposit(), pat::Photon::setIsolation(), pat::Photon::setIsolationPUPPI(), pat::Photon::setMaxDR(), pat::Photon::setMaxDRDEta(), pat::Photon::setMaxDRDPhi(), pat::Photon::setMaxDRRawEnergy(), pat::Photon::setPassElectronVeto(), reco::Photon::setPflowIsolationVariables(), pat::Photon::setPhotonIDs(), pat::helper::KinResolutionsLoader::setResolutions(), pat::Photon::setSee(), pat::Photon::setSeedEnergy(), pat::Photon::setSep(), pat::Photon::setSpp(), pat::Photon::setSubClusDEta1(), pat::Photon::setSubClusDEta2(), pat::Photon::setSubClusDEta3(), pat::Photon::setSubClusDPhi1(), pat::Photon::setSubClusDPhi2(), pat::Photon::setSubClusDPhi3(), pat::Photon::setSubClusRawE1(), pat::Photon::setSubClusRawE2(), pat::Photon::setSubClusRawE3(), EcalRegressionData::sigmaIEtaIEta(), EcalRegressionData::sigmaIEtaIPhi(), EcalRegressionData::sigmaIPhiIPhi(), edm::SortedCollection< T, SORT >::size(), jetUpdater_cfi::sort, edm::SortedCollection< T, SORT >::sort(), EcalRegressionData::subClusDEta(), EcalRegressionData::subClusDPhi(), EcalRegressionData::subClusRawEnergy(), reco::Photon::PflowIsolationVariables::sumEcalClusterEt, reco::Photon::PflowIsolationVariables::sumHcalClusterEt, reco::Photon::superCluster(), tier0::unique(), userDataHelper_, and useUserData_.

313  {
314  // switch off embedding (in unschedules mode)
315  if (iEvent.isRealData()) {
316  addGenMatch_ = false;
317  embedGenMatch_ = false;
318  }
319 
322 
323  // Get the vector of Photon's from the event
325  iEvent.getByToken(photonToken_, photons);
326 
327  // for conversion veto selection
329  iEvent.getByToken(hConversionsToken_, hConversions);
330 
331  // Get the collection of electrons from the event
333  iEvent.getByToken(electronToken_, hElectrons);
334 
335  // Get the beamspot
336  edm::Handle<reco::BeamSpot> beamSpotHandle;
337  iEvent.getByToken(beamLineToken_, beamSpotHandle);
338 
339  EcalClusterLazyTools lazyTools(iEvent,
343 
344  // prepare the MC matching
345  std::vector<edm::Handle<edm::Association<reco::GenParticleCollection>>> genMatches(genMatchTokens_.size());
346  if (addGenMatch_) {
347  for (size_t j = 0, nd = genMatchTokens_.size(); j < nd; ++j) {
348  iEvent.getByToken(genMatchTokens_[j], genMatches[j]);
349  }
350  }
351 
352  if (isolator_.enabled())
353  isolator_.beginEvent(iEvent, iSetup);
354 
359 
361  for (size_t j = 0, nd = isoDepositTokens_.size(); j < nd; ++j) {
362  iEvent.getByToken(isoDepositTokens_[j], deposits[j]);
363  }
364 
366  for (size_t j = 0; j < isolationValueTokens_.size(); ++j) {
368  }
369 
370  // prepare ID extraction
371  std::vector<edm::Handle<edm::ValueMap<Bool_t>>> idhandles;
372  std::vector<pat::Photon::IdPair> ids;
373  if (addPhotonID_) {
374  idhandles.resize(photIDSrcs_.size());
375  ids.resize(photIDSrcs_.size());
376  for (size_t i = 0; i < photIDSrcs_.size(); ++i) {
377  iEvent.getByToken(photIDTokens_[i], idhandles[i]);
378  ids[i].first = photIDSrcs_[i].first;
379  }
380  }
381 
382  //value maps for puppi isolation
383  edm::Handle<edm::ValueMap<float>> PUPPIIsolation_charged_hadrons;
384  edm::Handle<edm::ValueMap<float>> PUPPIIsolation_neutral_hadrons;
385  edm::Handle<edm::ValueMap<float>> PUPPIIsolation_photons;
386  if (addPuppiIsolation_) {
387  iEvent.getByToken(PUPPIIsolation_charged_hadrons_, PUPPIIsolation_charged_hadrons);
388  iEvent.getByToken(PUPPIIsolation_neutral_hadrons_, PUPPIIsolation_neutral_hadrons);
389  iEvent.getByToken(PUPPIIsolation_photons_, PUPPIIsolation_photons);
390  }
391 
392  // loop over photons
393  std::vector<Photon>* PATPhotons = new std::vector<Photon>();
394  for (edm::View<reco::Photon>::const_iterator itPhoton = photons->begin(); itPhoton != photons->end(); itPhoton++) {
395  // construct the Photon from the ref -> save ref to original object
396  unsigned int idx = itPhoton - photons->begin();
397  edm::RefToBase<reco::Photon> photonRef = photons->refAt(idx);
398  edm::Ptr<reco::Photon> photonPtr = photons->ptrAt(idx);
399  Photon aPhoton(photonRef);
400  auto phoPtr = photons->ptrAt(idx);
401  if (embedSuperCluster_)
402  aPhoton.embedSuperCluster();
403  if (embedSeedCluster_)
404  aPhoton.embedSeedCluster();
406  aPhoton.embedBasicClusters();
408  aPhoton.embedPreshowerClusters();
409 
410  std::vector<DetId> selectedCells;
411  bool barrel = itPhoton->isEB();
412  //loop over sub clusters
413  if (embedBasicClusters_) {
414  for (reco::CaloCluster_iterator clusIt = itPhoton->superCluster()->clustersBegin();
415  clusIt != itPhoton->superCluster()->clustersEnd();
416  ++clusIt) {
417  //get seed (max energy xtal)
418  DetId seed = lazyTools.getMaximum(**clusIt).first;
419  //get all xtals in 5x5 window around the seed
420  std::vector<DetId> dets5x5 =
423  selectedCells.insert(selectedCells.end(), dets5x5.begin(), dets5x5.end());
424 
425  //get all xtals belonging to cluster
426  for (const std::pair<DetId, float>& hit : (*clusIt)->hitsAndFractions()) {
427  selectedCells.push_back(hit.first);
428  }
429  }
430  }
431 
432  //remove duplicates
433  std::sort(selectedCells.begin(), selectedCells.end());
434  std::unique(selectedCells.begin(), selectedCells.end());
435 
436  // Retrieve the corresponding RecHits
437 
438  edm::Handle<EcalRecHitCollection> recHitsEBHandle;
439  iEvent.getByToken(reducedBarrelRecHitCollectionToken_, recHitsEBHandle);
440  edm::Handle<EcalRecHitCollection> recHitsEEHandle;
441  iEvent.getByToken(reducedEndcapRecHitCollectionToken_, recHitsEEHandle);
442 
443  //orginal code would throw an exception via the handle not being valid but now it'll just have a null pointer error
444  //should have little effect, if its not barrel or endcap, something very bad has happened elsewhere anyways
445  const EcalRecHitCollection* recHits = nullptr;
446  if (photonRef->superCluster()->seed()->hitsAndFractions().at(0).first.subdetId() == EcalBarrel)
447  recHits = recHitsEBHandle.product();
448  else if (photonRef->superCluster()->seed()->hitsAndFractions().at(0).first.subdetId() == EcalEndcap)
449  recHits = recHitsEEHandle.product();
450 
451  EcalRecHitCollection selectedRecHits;
452 
453  unsigned nSelectedCells = selectedCells.size();
454  for (unsigned icell = 0; icell < nSelectedCells; ++icell) {
455  EcalRecHitCollection::const_iterator it = recHits->find(selectedCells[icell]);
456  if (it != recHits->end()) {
457  selectedRecHits.push_back(*it);
458  }
459  }
460  selectedRecHits.sort();
461  if (embedRecHits_)
462  aPhoton.embedRecHits(&selectedRecHits);
463 
464  // store the match to the generated final state muons
465  if (addGenMatch_) {
466  for (size_t i = 0, n = genMatches.size(); i < n; ++i) {
467  reco::GenParticleRef genPhoton = (*genMatches[i])[photonRef];
468  aPhoton.addGenParticleRef(genPhoton);
469  }
470  if (embedGenMatch_)
471  aPhoton.embedGenParticle();
472  }
473 
474  if (efficiencyLoader_.enabled()) {
475  efficiencyLoader_.setEfficiencies(aPhoton, photonRef);
476  }
477 
478  if (resolutionLoader_.enabled()) {
480  }
481 
482  // here comes the extra functionality
483  if (isolator_.enabled()) {
485  typedef pat::helper::MultiIsolator::IsolationValuePairs IsolationValuePairs;
486  // better to loop backwards, so the vector is resized less times
487  for (IsolationValuePairs::const_reverse_iterator it = isolatorTmpStorage_.rbegin(),
488  ed = isolatorTmpStorage_.rend();
489  it != ed;
490  ++it) {
491  aPhoton.setIsolation(it->first, it->second);
492  }
493  }
494 
495  for (size_t j = 0, nd = deposits.size(); j < nd; ++j) {
496  aPhoton.setIsoDeposit(isoDepositLabels_[j].first, (*deposits[j])[photonRef]);
497  }
498 
499  for (size_t j = 0; j < isolationValues.size(); ++j) {
500  aPhoton.setIsolation(isolationValueLabels_[j].first, (*isolationValues[j])[photonRef]);
501  }
502 
503  // add photon ID info
504  if (addPhotonID_) {
505  for (size_t i = 0; i < photIDSrcs_.size(); ++i) {
506  ids[i].second = (*idhandles[i])[photonRef];
507  }
508  aPhoton.setPhotonIDs(ids);
509  }
510 
511  if (useUserData_) {
512  userDataHelper_.add(aPhoton, iEvent, iSetup);
513  }
514 
515  // set conversion veto selection
516  bool passelectronveto = false;
517  if (hConversions.isValid()) {
518  // this is recommended method
519  passelectronveto = !ConversionTools::hasMatchedPromptElectron(
520  photonRef->superCluster(), *hElectrons, *hConversions, beamSpotHandle->position());
521  }
522  aPhoton.setPassElectronVeto(passelectronveto);
523 
524  // set electron veto using pixel seed (not recommended but many analysis groups are still using since it is powerful method to remove electrons)
525  aPhoton.setHasPixelSeed(photonRef->hasPixelSeed());
526 
527  // set seed energy
528  aPhoton.setSeedEnergy(photonRef->superCluster()->seed()->energy());
529 
530  // set input variables for regression energy correction
531  if (saveRegressionData_) {
532  EcalRegressionData ecalRegData;
533  ecalRegData.fill(*(photonRef->superCluster()),
534  recHitsEBHandle.product(),
535  recHitsEEHandle.product(),
538  -1);
539 
540  aPhoton.setEMax(ecalRegData.eMax());
541  aPhoton.setE2nd(ecalRegData.e2nd());
542  aPhoton.setE3x3(ecalRegData.e3x3());
543  aPhoton.setETop(ecalRegData.eTop());
544  aPhoton.setEBottom(ecalRegData.eBottom());
545  aPhoton.setELeft(ecalRegData.eLeft());
546  aPhoton.setERight(ecalRegData.eRight());
547  aPhoton.setSee(ecalRegData.sigmaIEtaIEta());
548  aPhoton.setSep(
549  ecalRegData.sigmaIEtaIPhi() * ecalRegData.sigmaIEtaIEta() *
550  ecalRegData
551  .sigmaIPhiIPhi()); //there is a conflict on what sigmaIEtaIPhi actually is, regression and ID have it differently, this may change in later releases
552  aPhoton.setSpp(ecalRegData.sigmaIPhiIPhi());
553 
554  aPhoton.setMaxDR(ecalRegData.maxSubClusDR());
555  aPhoton.setMaxDRDPhi(ecalRegData.maxSubClusDRDPhi());
556  aPhoton.setMaxDRDEta(ecalRegData.maxSubClusDRDEta());
557  aPhoton.setMaxDRRawEnergy(ecalRegData.maxSubClusDRRawEnergy());
558  aPhoton.setSubClusRawE1(ecalRegData.subClusRawEnergy(EcalRegressionData::SubClusNr::C1));
559  aPhoton.setSubClusRawE2(ecalRegData.subClusRawEnergy(EcalRegressionData::SubClusNr::C2));
560  aPhoton.setSubClusRawE3(ecalRegData.subClusRawEnergy(EcalRegressionData::SubClusNr::C3));
561  aPhoton.setSubClusDPhi1(ecalRegData.subClusDPhi(EcalRegressionData::SubClusNr::C1));
562  aPhoton.setSubClusDPhi2(ecalRegData.subClusDPhi(EcalRegressionData::SubClusNr::C2));
563  aPhoton.setSubClusDPhi3(ecalRegData.subClusDPhi(EcalRegressionData::SubClusNr::C3));
564  aPhoton.setSubClusDEta1(ecalRegData.subClusDEta(EcalRegressionData::SubClusNr::C1));
565  aPhoton.setSubClusDEta2(ecalRegData.subClusDEta(EcalRegressionData::SubClusNr::C2));
566  aPhoton.setSubClusDEta3(ecalRegData.subClusDEta(EcalRegressionData::SubClusNr::C3));
567 
568  aPhoton.setCryPhi(ecalRegData.seedCrysPhiOrY());
569  aPhoton.setCryEta(ecalRegData.seedCrysEtaOrX());
570  aPhoton.setIEta(ecalRegData.seedCrysIEtaOrIX());
571  aPhoton.setIPhi(ecalRegData.seedCrysIPhiOrIY());
572  } else {
573  aPhoton.setEMax(0);
574  aPhoton.setE2nd(0);
575  aPhoton.setE3x3(0);
576  aPhoton.setETop(0);
577  aPhoton.setEBottom(0);
578  aPhoton.setELeft(0);
579  aPhoton.setERight(0);
580  aPhoton.setSee(0);
581  aPhoton.setSep(0);
582  aPhoton.setSpp(0);
583 
584  aPhoton.setMaxDR(0);
585  aPhoton.setMaxDRDPhi(0);
586  aPhoton.setMaxDRDEta(0);
587  aPhoton.setMaxDRRawEnergy(0);
588  aPhoton.setSubClusRawE1(0);
589  aPhoton.setSubClusRawE2(0);
590  aPhoton.setSubClusRawE3(0);
591  aPhoton.setSubClusDPhi1(0);
592  aPhoton.setSubClusDPhi2(0);
593  aPhoton.setSubClusDPhi3(0);
594  aPhoton.setSubClusDEta1(0);
595  aPhoton.setSubClusDEta2(0);
596  aPhoton.setSubClusDEta3(0);
597 
598  aPhoton.setCryPhi(0);
599  aPhoton.setCryEta(0);
600  aPhoton.setIEta(0);
601  aPhoton.setIPhi(0);
602  }
603 
604  if (addPuppiIsolation_)
605  aPhoton.setIsolationPUPPI((*PUPPIIsolation_charged_hadrons)[phoPtr],
606  (*PUPPIIsolation_neutral_hadrons)[phoPtr],
607  (*PUPPIIsolation_photons)[phoPtr]);
608  else
609  aPhoton.setIsolationPUPPI(-999., -999., -999.);
610 
611  // Get PFCluster Isolation
612  if (addPFClusterIso_) {
613  reco::Photon::PflowIsolationVariables newPFIsol = aPhoton.getPflowIsolationVariables();
614  edm::Handle<edm::ValueMap<float>> ecalPFClusterIsoMapH;
615  iEvent.getByToken(ecalPFClusterIsoT_, ecalPFClusterIsoMapH);
616  newPFIsol.sumEcalClusterEt = (*ecalPFClusterIsoMapH)[photonRef];
617  edm::Handle<edm::ValueMap<float>> hcalPFClusterIsoMapH;
619  iEvent.getByToken(hcalPFClusterIsoT_, hcalPFClusterIsoMapH);
620  newPFIsol.sumHcalClusterEt = (*hcalPFClusterIsoMapH)[photonRef];
621  } else {
622  newPFIsol.sumHcalClusterEt = -999.;
623  }
624  aPhoton.setPflowIsolationVariables(newPFIsol);
625  }
626 
627  // add the Photon to the vector of Photons
628  PATPhotons->push_back(aPhoton);
629  }
630 
631  // sort Photons in ET
632  std::sort(PATPhotons->begin(), PATPhotons->end(), eTComparator_);
633 
634  // put genEvt object in Event
635  std::unique_ptr<std::vector<Photon>> myPhotons(PATPhotons);
636  iEvent.put(std::move(myPhotons));
637  if (isolator_.enabled())
639 }
float maxSubClusDRDEta() const
bool hasPixelSeed() const
Bool flagging photons having a non-zero size vector of Ref to electornPixel seeds.
Definition: Photon.h:79
bool enabled() const
&#39;true&#39; if this there is at least one efficiency configured
float sigmaIPhiIPhi() const
void newEvent(const edm::Event &event)
To be called for each new event, reads in the ValueMaps for efficiencies.
float seedCrysPhiOrY() const
pat::PATUserDataHelper< pat::Photon > userDataHelper_
Definition: Photon.py:1
edm::EDGetTokenT< EcalRecHitCollection > reducedEndcapRecHitCollectionToken_
float maxSubClusDRDPhi() const
bool enabled() const
True if it has a non null configuration.
Definition: MultiIsolator.h:55
std::vector< edm::Handle< edm::ValueMap< IsoDeposit > > > IsoDepositMaps
const Point & position() const
position
Definition: BeamSpot.h:59
float seedCrysIPhiOrIY() const
void setEfficiencies(pat::PATObject< T > &obj, const R &originalRef) const
Sets the efficiencies for this object, using the reference to the original objects.
std::vector< edm::EDGetTokenT< edm::ValueMap< Bool_t > > > photIDTokens_
void fill(const reco::SuperCluster &superClus, const EcalRecHitCollection *ebRecHits, const EcalRecHitCollection *eeRecHits, const CaloGeometry *geom, const CaloTopology *topology, const reco::VertexCollection *vertices)
size_type size() const
float maxSubClusDRRawEnergy() const
static bool hasMatchedPromptElectron(const reco::SuperClusterRef &sc, const reco::GsfElectronCollection &eleCol, const reco::ConversionCollection &convCol, const math::XYZPoint &beamspot, bool allowCkfMatch=true, float lxyMin=2.0, float probMin=1e-6, unsigned int nHitsBeforeVtxMax=0)
edm::EDGetTokenT< reco::BeamSpot > beamLineToken_
T const * product() const
Definition: Handle.h:70
constexpr bool isUninitialized() const noexcept
Definition: EDGetToken.h:99
edm::EDGetTokenT< edm::View< reco::Photon > > photonToken_
std::vector< EcalRecHit >::const_iterator const_iterator
edm::EDGetTokenT< edm::ValueMap< float > > ecalPFClusterIsoT_
void push_back(T const &t)
std::vector< edm::EDGetTokenT< edm::Association< reco::GenParticleCollection > > > genMatchTokens_
float seedCrysEtaOrX() const
pat::helper::KinResolutionsLoader resolutionLoader_
edm::EDGetTokenT< edm::ValueMap< float > > hcalPFClusterIsoT_
GreaterByEt< Photon > eTComparator_
void setResolutions(pat::PATObject< T > &obj) const
Sets the efficiencies for this object, using the reference to the original objects.
const CaloGeometry * ecalGeometry_
float seedCrysIEtaOrIX() const
void fill(const edm::View< T > &coll, int idx, IsolationValuePairs &isolations) const
Definition: MultiIsolator.h:84
IsolationLabels isoDepositLabels_
const std::vector< float > & subClusDEta() const
ESData get(edm::EventSetup const &eventSetup) const
pat::helper::MultiIsolator::IsolationValuePairs isolatorTmpStorage_
int iEvent
Definition: GenABIO.cc:224
const std::vector< float > & subClusDPhi() const
void beginEvent(const edm::Event &event, const edm::EventSetup &eventSetup)
const CaloTopology * ecalTopology_
reco::SuperClusterRef superCluster() const override
Ref to SuperCluster.
IsolationLabels isolationValueLabels_
std::vector< edm::EDGetTokenT< edm::ValueMap< IsoDeposit > > > isoDepositTokens_
pat::helper::EfficiencyLoader efficiencyLoader_
void newEvent(const edm::Event &event, const edm::EventSetup &setup)
To be called for each new event, reads in the EventSetup object.
def unique(seq, keepstr=True)
Definition: tier0.py:24
bool enabled() const
&#39;true&#39; if this there is at least one efficiency configured
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_photons_
bool getData(T &iHolder) const
Definition: EventSetup.h:122
virtual std::vector< DetId > getWindow(const DetId &id, const int &northSouthSize, const int &eastWestSize) const
std::vector< NameTag > photIDSrcs_
pat::helper::MultiIsolator isolator_
void add(ObjectType &patObject, edm::Event const &iEvent, edm::EventSetup const &iSetup)
float sigmaIEtaIEta() const
Definition: DetId.h:17
std::vector< std::pair< pat::IsolationKeys, float > > IsolationValuePairs
Definition: MultiIsolator.h:17
float sigmaIEtaIPhi() const
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_charged_hadrons_
const CaloSubdetectorTopology * getSubdetectorTopology(const DetId &id) const
access the subdetector Topology for the given subdetector directly
Definition: CaloTopology.cc:17
float maxSubClusDR() const
bool isValid() const
Definition: HandleBase.h:70
const EcalClusterLazyTools::ESGetTokens ecalClusterToolsESGetTokens_
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
const edm::ESGetToken< CaloTopology, CaloTopologyRecord > ecalTopologyToken_
edm::EDGetTokenT< edm::ValueMap< float > > PUPPIIsolation_neutral_hadrons_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord > ecalGeometryToken_
edm::EDGetTokenT< EcalRecHitCollection > reducedBarrelRecHitCollectionToken_
std::vector< edm::EDGetTokenT< edm::ValueMap< double > > > isolationValueTokens_
std::vector< edm::Handle< edm::ValueMap< double > > > IsolationValueMaps
const std::vector< float > & subClusRawEnergy() const
def move(src, dest)
Definition: eostools.py:511
edm::EDGetTokenT< reco::ConversionCollection > hConversionsToken_
edm::EDGetTokenT< reco::GsfElectronCollection > electronToken_

◆ readIsolationLabels()

template<typename T >
void pat::PATPhotonProducer::readIsolationLabels ( const edm::ParameterSet iConfig,
const char *  psetName,
IsolationLabels labels,
std::vector< edm::EDGetTokenT< edm::ValueMap< T >>> &  tokens 
)
private

fill the labels vector from the contents of the parameter set, for the isodeposit or isolation values embedding

Definition at line 141 of file PATPhotonProducer.cc.

References pat::EcalIso, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), pat::HcalIso, crabWrapper::key, label, SummaryClient_cfi::labels, pat::PfAllParticleIso, pat::PfChargedAllIso, pat::PfChargedHadronIso, pat::PfGammaIso, pat::PfNeutralHadronIso, pat::PfPUChargedHadronIso, crab_NZS2018C_1::psetName, pat::TrackIso, pat::UserBaseIso, and edm::vector_transform().

144  {
145  labels.clear();
146 
147  if (iConfig.exists(psetName)) {
149 
150  if (depconf.exists("tracker"))
151  labels.push_back(std::make_pair(pat::TrackIso, depconf.getParameter<edm::InputTag>("tracker")));
152  if (depconf.exists("ecal"))
153  labels.push_back(std::make_pair(pat::EcalIso, depconf.getParameter<edm::InputTag>("ecal")));
154  if (depconf.exists("hcal"))
155  labels.push_back(std::make_pair(pat::HcalIso, depconf.getParameter<edm::InputTag>("hcal")));
156  if (depconf.exists("pfAllParticles")) {
157  labels.push_back(std::make_pair(pat::PfAllParticleIso, depconf.getParameter<edm::InputTag>("pfAllParticles")));
158  }
159  if (depconf.exists("pfChargedHadrons")) {
160  labels.push_back(
161  std::make_pair(pat::PfChargedHadronIso, depconf.getParameter<edm::InputTag>("pfChargedHadrons")));
162  }
163  if (depconf.exists("pfChargedAll")) {
164  labels.push_back(std::make_pair(pat::PfChargedAllIso, depconf.getParameter<edm::InputTag>("pfChargedAll")));
165  }
166  if (depconf.exists("pfPUChargedHadrons")) {
167  labels.push_back(
168  std::make_pair(pat::PfPUChargedHadronIso, depconf.getParameter<edm::InputTag>("pfPUChargedHadrons")));
169  }
170  if (depconf.exists("pfNeutralHadrons")) {
171  labels.push_back(
172  std::make_pair(pat::PfNeutralHadronIso, depconf.getParameter<edm::InputTag>("pfNeutralHadrons")));
173  }
174  if (depconf.exists("pfPhotons")) {
175  labels.push_back(std::make_pair(pat::PfGammaIso, depconf.getParameter<edm::InputTag>("pfPhotons")));
176  }
177  if (depconf.exists("user")) {
178  std::vector<edm::InputTag> userdeps = depconf.getParameter<std::vector<edm::InputTag>>("user");
179  std::vector<edm::InputTag>::const_iterator it = userdeps.begin(), ed = userdeps.end();
181  for (; it != ed; ++it, ++key) {
182  labels.push_back(std::make_pair(pat::IsolationKeys(key), *it));
183  }
184  }
185 
186  tokens = edm::vector_transform(
187  labels, [this](IsolationLabel const& label) { return consumes<edm::ValueMap<T>>(label.second); });
188  }
189  tokens = edm::vector_transform(
190  labels, [this](IsolationLabel const& label) { return consumes<edm::ValueMap<T>>(label.second); });
191 }
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
bool exists(std::string const &parameterName) const
checks if a parameter exists
IsolationKeys
Enum defining isolation keys.
Definition: Isolation.h:9
auto vector_transform(std::vector< InputType > const &input, Function predicate) -> std::vector< typename std::remove_cv< typename std::remove_reference< decltype(predicate(input.front()))>::type >::type >
Definition: transform.h:11
char const * label
std::pair< pat::IsolationKeys, edm::InputTag > IsolationLabel

Member Data Documentation

◆ addEfficiencies_

bool pat::PATPhotonProducer::addEfficiencies_
private

Definition at line 111 of file PATPhotonProducer.cc.

◆ addGenMatch_

bool pat::PATPhotonProducer::addGenMatch_
private

Definition at line 83 of file PATPhotonProducer.cc.

Referenced by produce().

◆ addPFClusterIso_

bool pat::PATPhotonProducer::addPFClusterIso_
private

Definition at line 78 of file PATPhotonProducer.cc.

Referenced by produce().

◆ addPhotonID_

bool pat::PATPhotonProducer::addPhotonID_
private

Definition at line 117 of file PATPhotonProducer.cc.

Referenced by produce().

◆ addPuppiIsolation_

bool pat::PATPhotonProducer::addPuppiIsolation_
private

Definition at line 79 of file PATPhotonProducer.cc.

Referenced by produce().

◆ addResolutions_

bool pat::PATPhotonProducer::addResolutions_
private

Definition at line 114 of file PATPhotonProducer.cc.

◆ beamLineToken_

edm::EDGetTokenT<reco::BeamSpot> pat::PATPhotonProducer::beamLineToken_
private

Definition at line 63 of file PATPhotonProducer.cc.

Referenced by produce().

◆ ecalClusterToolsESGetTokens_

const EcalClusterLazyTools::ESGetTokens pat::PATPhotonProducer::ecalClusterToolsESGetTokens_
private

Definition at line 76 of file PATPhotonProducer.cc.

Referenced by produce().

◆ ecalGeometry_

const CaloGeometry* pat::PATPhotonProducer::ecalGeometry_
private

Definition at line 130 of file PATPhotonProducer.cc.

Referenced by produce().

◆ ecalGeometryToken_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> pat::PATPhotonProducer::ecalGeometryToken_
private

Definition at line 135 of file PATPhotonProducer.cc.

Referenced by produce().

◆ ecalPFClusterIsoT_

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATPhotonProducer::ecalPFClusterIsoT_
private

Definition at line 80 of file PATPhotonProducer.cc.

Referenced by produce().

◆ ecalTopology_

const CaloTopology* pat::PATPhotonProducer::ecalTopology_
private

Definition at line 129 of file PATPhotonProducer.cc.

Referenced by produce().

◆ ecalTopologyToken_

const edm::ESGetToken<CaloTopology, CaloTopologyRecord> pat::PATPhotonProducer::ecalTopologyToken_
private

Definition at line 134 of file PATPhotonProducer.cc.

Referenced by produce().

◆ efficiencyLoader_

pat::helper::EfficiencyLoader pat::PATPhotonProducer::efficiencyLoader_
private

Definition at line 112 of file PATPhotonProducer.cc.

Referenced by produce().

◆ electronToken_

edm::EDGetTokenT<reco::GsfElectronCollection> pat::PATPhotonProducer::electronToken_
private

Definition at line 61 of file PATPhotonProducer.cc.

Referenced by produce().

◆ embedBasicClusters_

bool pat::PATPhotonProducer::embedBasicClusters_
private

Definition at line 67 of file PATPhotonProducer.cc.

Referenced by produce().

◆ embedGenMatch_

bool pat::PATPhotonProducer::embedGenMatch_
private

Definition at line 84 of file PATPhotonProducer.cc.

Referenced by produce().

◆ embedPreshowerClusters_

bool pat::PATPhotonProducer::embedPreshowerClusters_
private

Definition at line 68 of file PATPhotonProducer.cc.

Referenced by produce().

◆ embedRecHits_

bool pat::PATPhotonProducer::embedRecHits_
private

Definition at line 69 of file PATPhotonProducer.cc.

Referenced by produce().

◆ embedSeedCluster_

bool pat::PATPhotonProducer::embedSeedCluster_
private

Definition at line 66 of file PATPhotonProducer.cc.

Referenced by produce().

◆ embedSuperCluster_

bool pat::PATPhotonProducer::embedSuperCluster_
private

Definition at line 65 of file PATPhotonProducer.cc.

Referenced by produce().

◆ eTComparator_

GreaterByEt<Photon> pat::PATPhotonProducer::eTComparator_
private

Definition at line 88 of file PATPhotonProducer.cc.

Referenced by produce().

◆ genMatchTokens_

std::vector<edm::EDGetTokenT<edm::Association<reco::GenParticleCollection> > > pat::PATPhotonProducer::genMatchTokens_
private

Definition at line 85 of file PATPhotonProducer.cc.

Referenced by produce().

◆ hcalPFClusterIsoT_

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATPhotonProducer::hcalPFClusterIsoT_
private

Definition at line 81 of file PATPhotonProducer.cc.

Referenced by produce().

◆ hConversionsToken_

edm::EDGetTokenT<reco::ConversionCollection> pat::PATPhotonProducer::hConversionsToken_
private

Definition at line 62 of file PATPhotonProducer.cc.

Referenced by produce().

◆ isoDepositLabels_

IsolationLabels pat::PATPhotonProducer::isoDepositLabels_
private

Definition at line 100 of file PATPhotonProducer.cc.

Referenced by produce().

◆ isoDepositTokens_

std::vector<edm::EDGetTokenT<edm::ValueMap<IsoDeposit> > > pat::PATPhotonProducer::isoDepositTokens_
private

Definition at line 97 of file PATPhotonProducer.cc.

Referenced by produce().

◆ isolationValueLabels_

IsolationLabels pat::PATPhotonProducer::isolationValueLabels_
private

Definition at line 101 of file PATPhotonProducer.cc.

Referenced by produce().

◆ isolationValueTokens_

std::vector<edm::EDGetTokenT<edm::ValueMap<double> > > pat::PATPhotonProducer::isolationValueTokens_
private

Definition at line 98 of file PATPhotonProducer.cc.

Referenced by produce().

◆ isolator_

pat::helper::MultiIsolator pat::PATPhotonProducer::isolator_
private

Definition at line 95 of file PATPhotonProducer.cc.

Referenced by produce().

◆ isolatorTmpStorage_

pat::helper::MultiIsolator::IsolationValuePairs pat::PATPhotonProducer::isolatorTmpStorage_
private

Definition at line 96 of file PATPhotonProducer.cc.

Referenced by produce().

◆ photIDSrcs_

std::vector<NameTag> pat::PATPhotonProducer::photIDSrcs_
private

Definition at line 119 of file PATPhotonProducer.cc.

Referenced by produce().

◆ photIDTokens_

std::vector<edm::EDGetTokenT<edm::ValueMap<Bool_t> > > pat::PATPhotonProducer::photIDTokens_
private

Definition at line 120 of file PATPhotonProducer.cc.

Referenced by produce().

◆ photonToken_

edm::EDGetTokenT<edm::View<reco::Photon> > pat::PATPhotonProducer::photonToken_
private

Definition at line 60 of file PATPhotonProducer.cc.

Referenced by produce().

◆ PUPPIIsolation_charged_hadrons_

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATPhotonProducer::PUPPIIsolation_charged_hadrons_
private

Definition at line 124 of file PATPhotonProducer.cc.

Referenced by produce().

◆ PUPPIIsolation_neutral_hadrons_

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATPhotonProducer::PUPPIIsolation_neutral_hadrons_
private

Definition at line 125 of file PATPhotonProducer.cc.

Referenced by produce().

◆ PUPPIIsolation_photons_

edm::EDGetTokenT<edm::ValueMap<float> > pat::PATPhotonProducer::PUPPIIsolation_photons_
private

Definition at line 126 of file PATPhotonProducer.cc.

Referenced by produce().

◆ reducedBarrelRecHitCollection_

edm::InputTag pat::PATPhotonProducer::reducedBarrelRecHitCollection_
private

Definition at line 71 of file PATPhotonProducer.cc.

◆ reducedBarrelRecHitCollectionToken_

edm::EDGetTokenT<EcalRecHitCollection> pat::PATPhotonProducer::reducedBarrelRecHitCollectionToken_
private

Definition at line 72 of file PATPhotonProducer.cc.

Referenced by produce().

◆ reducedEndcapRecHitCollection_

edm::InputTag pat::PATPhotonProducer::reducedEndcapRecHitCollection_
private

Definition at line 73 of file PATPhotonProducer.cc.

◆ reducedEndcapRecHitCollectionToken_

edm::EDGetTokenT<EcalRecHitCollection> pat::PATPhotonProducer::reducedEndcapRecHitCollectionToken_
private

Definition at line 74 of file PATPhotonProducer.cc.

Referenced by produce().

◆ resolutionLoader_

pat::helper::KinResolutionsLoader pat::PATPhotonProducer::resolutionLoader_
private

Definition at line 115 of file PATPhotonProducer.cc.

Referenced by produce().

◆ saveRegressionData_

bool pat::PATPhotonProducer::saveRegressionData_
private

Definition at line 132 of file PATPhotonProducer.cc.

Referenced by produce().

◆ userDataHelper_

pat::PATUserDataHelper<pat::Photon> pat::PATPhotonProducer::userDataHelper_
private

Definition at line 127 of file PATPhotonProducer.cc.

Referenced by produce().

◆ useUserData_

bool pat::PATPhotonProducer::useUserData_
private

Definition at line 122 of file PATPhotonProducer.cc.

Referenced by produce().