slimmer of PAT Taus More...

Inheritance diagram for pat::PATPhotonSlimmer:

Public Member Functions
virtual void	beginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) overridefinal

	PATPhotonSlimmer (const edm::ParameterSet &iConfig)

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

virtual	~PATPhotonSlimmer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
StringCutObjectSelector < pat::Photon >	dropBasicClusters_

StringCutObjectSelector < pat::Photon >	dropPreshowerClusters_

StringCutObjectSelector < pat::Photon >	dropRecHits_

StringCutObjectSelector < pat::Photon >	dropSeedCluster_

StringCutObjectSelector < pat::Photon >	dropSuperClusters_

bool	linkToPackedPF_

bool	modifyPhoton_

edm::EDGetTokenT < pat::PackedCandidateCollection >	pc_

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

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

edm::EDGetTokenT < edm::ValueMap< std::vector < reco::PFCandidateRef > > >	reco2pf_

edm::EDGetTokenT < EcalRecHitCollection >	reducedBarrelRecHitCollectionToken_

edm::EDGetTokenT < EcalRecHitCollection >	reducedEndcapRecHitCollectionToken_

StringCutObjectSelector < pat::Photon >	saveNonZSClusterShapes_

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

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

slimmer of PAT Taus

Definition at line 27 of file PATPhotonSlimmer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 52 of file PATPhotonSlimmer.cc.

References edm::EDConsumerBase::consumesCollector(), edm::ParameterSet::getParameter(), and HLT_25ns14e33_v1_cff::InputTag.

                                                                      :
     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")),
     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"))
 {
     edm::ConsumesCollector sumes(consumesCollector());
     if( modifyPhoton_ ) {
       const edm::ParameterSet& mod_config = iConfig.getParameter<edm::ParameterSet>("modifierConfig");
       photonModifier_.reset(new pat::ObjectModifier<pat::Photon>(mod_config) );
       photonModifier_->setConsumes(sumes);
     } else {
       photonModifier_.reset(nullptr);
     }
     
     produces<std::vector<pat::Photon> >();
     if (linkToPackedPF_) {
         reco2pf_ = consumes<edm::ValueMap<std::vector<reco::PFCandidateRef>>>(iConfig.getParameter<edm::InputTag>("recoToPFMap"));
         pf2pc_   = consumes<edm::Association<pat::PackedCandidateCollection>>(iConfig.getParameter<edm::InputTag>("packedPFCandidates"));
         pc_   = consumes<pat::PackedCandidateCollection>(iConfig.getParameter<edm::InputTag>("packedPFCandidates"));
     }
     mayConsume<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEB"));
     mayConsume<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEE"));
 }

virtual pat::PATPhotonSlimmer::~PATPhotonSlimmer ( )

inlinevirtual

Definition at line 30 of file PATPhotonSlimmer.cc.

30 { }

Member Function Documentation

void pat::PATPhotonSlimmer::beginLuminosityBlock	(	const edm::LuminosityBlock &	,
		const edm::EventSetup &	iSetup
	)

finaloverridevirtual

Reimplemented from edm::EDProducer.

Definition at line 85 of file PATPhotonSlimmer.cc.

85 {

86 }

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

virtual

Implements edm::EDProducer.

Definition at line 89 of file PATPhotonSlimmer.cc.

References pat::PATObject< ObjectType >::addUserFloat(), pat::Photon::basicClusters_, dropBasicClusters_, dropPreshowerClusters_, dropRecHits_, dropSeedCluster_, dropSuperClusters_, EcalClusterLazyToolsT< EcalClusterToolsImpl >::e1x5(), EcalClusterLazyToolsT< EcalClusterToolsImpl >::e3x3(), EcalClusterLazyToolsT< EcalClusterToolsImpl >::e5x5(), pat::Photon::embeddedRecHits_, pat::Photon::embeddedSeedCluster_, pat::Photon::embeddedSuperCluster_, edm::Event::getByToken(), edm::isNotFinite(), relval_steps::key, relativeConstraints::keys, linkToPackedPF_, EcalClusterLazyToolsT< EcalClusterToolsImpl >::localCovariances(), modifyPhoton_, dbtoconf::out, pc_, pf2pc_, photonModifier_, pat::Photon::preshowerClusters_, edm::Event::put(), pat::Photon::recHits_, reco2pf_, reducedBarrelRecHitCollectionToken_, reducedEndcapRecHitCollectionToken_, pat::PATObject< ObjectType >::refToOrig_, saveNonZSClusterShapes_, pat::Photon::seedCluster_, pat::Photon::setAssociatedPackedPFCandidates(), pat::Photon::sourceCandidatePtr(), mathSSE::sqrt(), alcazmumu_cfi::src, src_, pat::Photon::superCluster(), and pat::Photon::superCluster_.

                                                                             {
     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, iSetup, reducedBarrelRecHitCollectionToken_, reducedEndcapRecHitCollectionToken_);
 
     auto_ptr<vector<pat::Photon> >  out(new 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 (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)) {
             std::vector<float> 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(out);
 }

Member Data Documentation

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

private

Definition at line 38 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 38 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 38 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 38 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 38 of file PATPhotonSlimmer.cc.

Referenced by produce().

bool pat::PATPhotonSlimmer::linkToPackedPF_

private

Definition at line 43 of file PATPhotonSlimmer.cc.

Referenced by produce().

bool pat::PATPhotonSlimmer::modifyPhoton_

private

Definition at line 46 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 42 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 41 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 47 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 40 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 45 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 45 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 44 of file PATPhotonSlimmer.cc.

Referenced by produce().

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

private

Definition at line 36 of file PATPhotonSlimmer.cc.

Referenced by produce().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation