slimmer of PAT Taus More...

Inheritance diagram for pat::PATPhotonSlimmer:

Public Member Functions
	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
	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

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_

edm::EDGetTokenT < pat::PackedCandidateCollection >	pc_

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

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 26 of file PATPhotonSlimmer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 48 of file PATPhotonSlimmer.cc.

References edm::ParameterSet::getParameter(), linkToPackedPF_, pc_, pf2pc_, and reco2pf_.

                                                                      :
     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")))
 {
     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 29 of file PATPhotonSlimmer.cc.

29 { }

Member Function Documentation

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

virtual

Implements edm::EDProducer.

Definition at line 71 of file PATPhotonSlimmer.cc.

References pat::PATObject< ObjectType >::addUserFloat(), pat::Photon::basicClusters_, EcalClusterLazyToolsT< EcalClusterToolsImpl >::e1x5(), EcalClusterLazyToolsT< EcalClusterToolsImpl >::e3x3(), EcalClusterLazyToolsT< EcalClusterToolsImpl >::e5x5(), pat::Photon::embeddedRecHits_, pat::Photon::embeddedSeedCluster_, pat::Photon::embeddedSuperCluster_, edm::Event::getByToken(), edm::isNotFinite(), EcalClusterLazyToolsT< EcalClusterToolsImpl >::localCovariances(), dbtoconf::out, pat::Photon::preshowerClusters_, edm::RefVector< C, T, F >::push_back(), edm::Event::put(), pat::Photon::recHits_, pat::PATObject< ObjectType >::refToOrig_, edm::refToPtr(), pat::Photon::seedCluster_, pat::Photon::setAssociatedPackedPFCandidates(), pat::Photon::setPackedPFCandidateCollection(), edm::RefVector< C, T, F >::size(), mathSSE::sqrt(), alcazmumu_cfi::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());
 
     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 (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_) {
             photon.setPackedPFCandidateCollection(edm::RefProd<pat::PackedCandidateCollection>(pc));
             //std::cout << " PAT  photon in  " << src.id() << " comes from " << photon.refToOrig_.id() << ", " << photon.refToOrig_.key() << std::endl;
             edm::RefVector<pat::PackedCandidateCollection> origs;
             for (const reco::PFCandidateRef & pf : (*reco2pf)[photon.refToOrig_]) {
                 if (pf2pc->contains(pf.id())) {
                     origs.push_back((*pf2pc)[pf]);
                 } //else std::cerr << " Photon linked to a PFCand in " << pf.id() << " while we expect them in " << pf2pc->ids().front().first << "\n";
             }
             //std::cout << "Photon with pt " << photon.pt() << " associated to " << origs.size() << " PF Candidates\n";
             photon.setAssociatedPackedPFCandidates(origs);
             //if there's just one PF Cand then it's me, otherwise I have no univoque parent so my ref will be null 
             if (origs.size() == 1) {
                 photon.refToOrig_ = refToPtr(origs[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);
 }