pat::PATPhotonSlimmer Class Reference

slimmer of PAT Taus More...

Inheritance diagram for pat::PATPhotonSlimmer:

Public Member Functions
	PATPhotonSlimmer (const edm::ParameterSet &iConfig)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) 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 EDProducer &	operator= (const EDProducer &)=delete

Private Attributes
const StringCutObjectSelector< pat::Photon >	dropBasicClusters_
const StringCutObjectSelector< pat::Photon >	dropPreshowerClusters_
const StringCutObjectSelector< pat::Photon >	dropRecHits_
const StringCutObjectSelector< pat::Photon >	dropRegressionData_
const StringCutObjectSelector< pat::Photon >	dropSaturation_
const StringCutObjectSelector< pat::Photon >	dropSeedCluster_
const StringCutObjectSelector< pat::Photon >	dropSuperClusters_
const EcalClusterLazyTools::ESGetTokens	ecalClusterToolsESGetTokens_
const bool	linkToPackedPF_
const bool	modifyPhoton_
const edm::EDGetTokenT< pat::PackedCandidateCollection >	pc_
const edm::EDGetTokenT< edm::Association< pat::PackedCandidateCollection > >	pf2pc_
std::unique_ptr< pat::ObjectModifier< pat::Photon > >	photonModifier_
const edm::EDGetTokenT< edm::ValueMap< std::vector< reco::PFCandidateRef > > >	reco2pf_
const edm::EDGetTokenT< EcalRecHitCollection >	reducedBarrelRecHitCollectionToken_
const edm::EDGetTokenT< EcalRecHitCollection >	reducedEndcapRecHitCollectionToken_
const StringCutObjectSelector< pat::Photon >	saveNonZSClusterShapes_
const edm::EDGetTokenT< edm::View< pat::Photon > >	src_

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

slimmer of PAT Taus

Definition at line 28 of file PATPhotonSlimmer.cc.

Constructor & Destructor Documentation

PATPhotonSlimmer()

pat::PATPhotonSlimmer::PATPhotonSlimmer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 54 of file PATPhotonSlimmer.cc.

    : src_(consumes<edm::View<pat::Photon>>(iConfig.getParameter<edm::InputTag>("src"))),
      dropSuperClusters_(iConfig.getParameter<std::string>("dropSuperCluster")),
      dropBasicClusters_(iConfig.getParameter<std::string>("dropBasicClusters")),
      dropPreshowerClusters_(iConfig.getParameter<std::string>("dropPreshowerClusters")),
      dropSeedCluster_(iConfig.getParameter<std::string>("dropSeedCluster")),
      dropRecHits_(iConfig.getParameter<std::string>("dropRecHits")),
      dropSaturation_(iConfig.getParameter<std::string>("dropSaturation")),
      dropRegressionData_(iConfig.getParameter<std::string>("dropRegressionData")),
      reco2pf_(mayConsume<edm::ValueMap<std::vector<reco::PFCandidateRef>>>(
          iConfig.getParameter<edm::InputTag>("recoToPFMap"))),
      pf2pc_(mayConsume<edm::Association<pat::PackedCandidateCollection>>(
          iConfig.getParameter<edm::InputTag>("packedPFCandidates"))),
      pc_(mayConsume<pat::PackedCandidateCollection>(iConfig.getParameter<edm::InputTag>("packedPFCandidates"))),
      linkToPackedPF_(iConfig.getParameter<bool>("linkToPackedPFCandidates")),
      saveNonZSClusterShapes_(iConfig.getParameter<std::string>("saveNonZSClusterShapes")),
      reducedBarrelRecHitCollectionToken_(
          consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("reducedBarrelRecHitCollection"))),
      reducedEndcapRecHitCollectionToken_(
          consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("reducedEndcapRecHitCollection"))),
      modifyPhoton_(iConfig.getParameter<bool>("modifyPhotons")),
      ecalClusterToolsESGetTokens_{consumesCollector()} {
  if (modifyPhoton_) {
    const edm::ParameterSet& mod_config = iConfig.getParameter<edm::ParameterSet>("modifierConfig");
    photonModifier_ = std::make_unique<pat::ObjectModifier<pat::Photon>>(mod_config, consumesCollector());
  }

  mayConsume<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEB"));
  mayConsume<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEE"));

  produces<std::vector<pat::Photon>>();
}

Member Function Documentation

produce()

void pat::PATPhotonSlimmer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 87 of file PATPhotonSlimmer.cc.

References edm::HandleBase::id(), iEvent, edm::isNotFinite(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, submitPVResolutionJobs::key, relativeConstraints::keys, eostools::move(), MillePedeFileConverter_cfg::out, packedPFCandidateRefMixer_cfi::pf, displacedMuons_cfi::photon, electrons_cff::r9, run3scouting_cff::sigmaIetaIeta, mathSSE::sqrt(), and TrackRefitter_38T_cff::src.

                                                                               {
  using namespace edm;
  using namespace std;

  Handle<View<pat::Photon>> src;
  iEvent.getByToken(src_, src);

  Handle<edm::ValueMap<std::vector<reco::PFCandidateRef>>> reco2pf;
  Handle<edm::Association<pat::PackedCandidateCollection>> pf2pc;
  Handle<pat::PackedCandidateCollection> pc;
  if (linkToPackedPF_) {
    iEvent.getByToken(reco2pf_, reco2pf);
    iEvent.getByToken(pf2pc_, pf2pc);
    iEvent.getByToken(pc_, pc);
  }
  noZS::EcalClusterLazyTools lazyToolsNoZS(iEvent,
                                           ecalClusterToolsESGetTokens_.get(iSetup),
                                           reducedBarrelRecHitCollectionToken_,
                                           reducedEndcapRecHitCollectionToken_);

  auto out = std::make_unique<std::vector<pat::Photon>>();
  out->reserve(src->size());

  if (modifyPhoton_) {
    photonModifier_->setEvent(iEvent);
  }
  if (modifyPhoton_)
    photonModifier_->setEventContent(iSetup);

  std::vector<unsigned int> keys;
  for (View<pat::Photon>::const_iterator it = src->begin(), ed = src->end(); it != ed; ++it) {
    out->push_back(*it);
    pat::Photon& photon = out->back();

    if (modifyPhoton_) {
      photonModifier_->modify(photon);
    }

    if (dropSuperClusters_(photon)) {
      photon.superCluster_.clear();
      photon.embeddedSuperCluster_ = false;
    }
    if (dropBasicClusters_(photon)) {
      photon.basicClusters_.clear();
    }
    if (dropPreshowerClusters_(photon)) {
      photon.preshowerClusters_.clear();
    }
    if (dropSeedCluster_(photon)) {
      photon.seedCluster_.clear();
      photon.embeddedSeedCluster_ = false;
    }
    if (dropRecHits_(photon)) {
      photon.recHits_ = EcalRecHitCollection();
      photon.embeddedRecHits_ = false;
    }
    if (dropSaturation_(photon)) {
      photon.setSaturationInfo(reco::Photon::SaturationInfo());
    }
    if (dropRegressionData_(photon)) {
      photon.setEMax(0);
      photon.setE2nd(0);
      photon.setE3x3(0);
      photon.setETop(0);
      photon.setEBottom(0);
      photon.setELeft(0);
      photon.setERight(0);
      photon.setSee(0);
      photon.setSep(0);
      photon.setSpp(0);
      photon.setMaxDR(0);
      photon.setMaxDRDPhi(0);
      photon.setMaxDRDEta(0);
      photon.setMaxDRRawEnergy(0);
      photon.setSubClusRawE1(0);
      photon.setSubClusRawE2(0);
      photon.setSubClusRawE3(0);
      photon.setSubClusDPhi1(0);
      photon.setSubClusDPhi2(0);
      photon.setSubClusDPhi3(0);
      photon.setSubClusDEta1(0);
      photon.setSubClusDEta2(0);
      photon.setSubClusDEta3(0);
      photon.setCryPhi(0);
      photon.setCryEta(0);
      photon.setIEta(0);
      photon.setIPhi(0);
    }

    if (linkToPackedPF_) {
      //std::cout << " PAT  photon in  " << src.id() << " comes from " << photon.refToOrig_.id() << ", " << photon.refToOrig_.key() << std::endl;
      keys.clear();
      for (auto const& pf : (*reco2pf)[photon.refToOrig_]) {
        if (pf2pc->contains(pf.id())) {
          keys.push_back((*pf2pc)[pf].key());
        }
      }
      photon.setAssociatedPackedPFCandidates(
          edm::RefProd<pat::PackedCandidateCollection>(pc), keys.begin(), keys.end());
      //std::cout << "Photon with pt " << photon.pt() << " associated to " << photon.associatedPackedFCandidateIndices_.size() << " PF Candidates\n";
      //if there's just one PF Cand then it's me, otherwise I have no univoque parent so my ref will be null
      if (keys.size() == 1) {
        photon.refToOrig_ = photon.sourceCandidatePtr(0);
      } else {
        photon.refToOrig_ = reco::CandidatePtr(pc.id());
      }
    }
    if (saveNonZSClusterShapes_(photon)) {
      const auto& vCov = lazyToolsNoZS.localCovariances(*(photon.superCluster()->seed()));
      float r9 = lazyToolsNoZS.e3x3(*(photon.superCluster()->seed())) / photon.superCluster()->rawEnergy();
      float sigmaIetaIeta = (!edm::isNotFinite(vCov[0])) ? sqrt(vCov[0]) : 0;
      float sigmaIetaIphi = vCov[1];
      float sigmaIphiIphi = (!edm::isNotFinite(vCov[2])) ? sqrt(vCov[2]) : 0;
      float e15o55 =
          lazyToolsNoZS.e1x5(*(photon.superCluster()->seed())) / lazyToolsNoZS.e5x5(*(photon.superCluster()->seed()));
      photon.addUserFloat("sigmaIetaIeta_NoZS", sigmaIetaIeta);
      photon.addUserFloat("sigmaIetaIphi_NoZS", sigmaIetaIphi);
      photon.addUserFloat("sigmaIphiIphi_NoZS", sigmaIphiIphi);
      photon.addUserFloat("r9_NoZS", r9);
      photon.addUserFloat("e1x5_over_e5x5_NoZS", e15o55);
    }
  }

  iEvent.put(std::move(out));
}

Member Data Documentation

dropBasicClusters_

const StringCutObjectSelector<pat::Photon> pat::PATPhotonSlimmer::dropBasicClusters_

private

Definition at line 37 of file PATPhotonSlimmer.cc.

dropPreshowerClusters_

const StringCutObjectSelector<pat::Photon> pat::PATPhotonSlimmer::dropPreshowerClusters_

private

Definition at line 37 of file PATPhotonSlimmer.cc.

dropRecHits_

const StringCutObjectSelector<pat::Photon> pat::PATPhotonSlimmer::dropRecHits_

private

Definition at line 37 of file PATPhotonSlimmer.cc.

dropRegressionData_

const StringCutObjectSelector<pat::Photon> pat::PATPhotonSlimmer::dropRegressionData_

private

Definition at line 37 of file PATPhotonSlimmer.cc.

dropSaturation_

const StringCutObjectSelector<pat::Photon> pat::PATPhotonSlimmer::dropSaturation_

private

Definition at line 37 of file PATPhotonSlimmer.cc.

dropSeedCluster_

const StringCutObjectSelector<pat::Photon> pat::PATPhotonSlimmer::dropSeedCluster_

private

Definition at line 37 of file PATPhotonSlimmer.cc.

dropSuperClusters_

const StringCutObjectSelector<pat::Photon> pat::PATPhotonSlimmer::dropSuperClusters_

private

Definition at line 37 of file PATPhotonSlimmer.cc.

ecalClusterToolsESGetTokens_

const EcalClusterLazyTools::ESGetTokens pat::PATPhotonSlimmer::ecalClusterToolsESGetTokens_

private

Definition at line 49 of file PATPhotonSlimmer.cc.

linkToPackedPF_

const bool pat::PATPhotonSlimmer::linkToPackedPF_

private

Definition at line 43 of file PATPhotonSlimmer.cc.

modifyPhoton_

const bool pat::PATPhotonSlimmer::modifyPhoton_

private

Definition at line 47 of file PATPhotonSlimmer.cc.

pc_

const edm::EDGetTokenT<pat::PackedCandidateCollection> pat::PATPhotonSlimmer::pc_

private

Definition at line 42 of file PATPhotonSlimmer.cc.

pf2pc_

const edm::EDGetTokenT<edm::Association<pat::PackedCandidateCollection> > pat::PATPhotonSlimmer::pf2pc_

private

Definition at line 41 of file PATPhotonSlimmer.cc.

photonModifier_

std::unique_ptr<pat::ObjectModifier<pat::Photon> > pat::PATPhotonSlimmer::photonModifier_

private

Definition at line 48 of file PATPhotonSlimmer.cc.

reco2pf_

const edm::EDGetTokenT<edm::ValueMap<std::vector<reco::PFCandidateRef> > > pat::PATPhotonSlimmer::reco2pf_

private

Definition at line 40 of file PATPhotonSlimmer.cc.

reducedBarrelRecHitCollectionToken_

const edm::EDGetTokenT<EcalRecHitCollection> pat::PATPhotonSlimmer::reducedBarrelRecHitCollectionToken_

private

Definition at line 45 of file PATPhotonSlimmer.cc.

reducedEndcapRecHitCollectionToken_

const edm::EDGetTokenT<EcalRecHitCollection> pat::PATPhotonSlimmer::reducedEndcapRecHitCollectionToken_

private

Definition at line 45 of file PATPhotonSlimmer.cc.

saveNonZSClusterShapes_

const StringCutObjectSelector<pat::Photon> pat::PATPhotonSlimmer::saveNonZSClusterShapes_

private

Definition at line 44 of file PATPhotonSlimmer.cc.

src_

const edm::EDGetTokenT<edm::View<pat::Photon> > pat::PATPhotonSlimmer::src_

private

Definition at line 35 of file PATPhotonSlimmer.cc.