CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Attributes
pat::PATPhotonSlimmer Class Reference

slimmer of PAT Taus More...

Inheritance diagram for pat::PATPhotonSlimmer:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

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< ConsumesInfoconsumesInfo () 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.

52  :
54  dropSuperClusters_(iConfig.getParameter<std::string>("dropSuperCluster")),
55  dropBasicClusters_(iConfig.getParameter<std::string>("dropBasicClusters")),
56  dropPreshowerClusters_(iConfig.getParameter<std::string>("dropPreshowerClusters")),
57  dropSeedCluster_(iConfig.getParameter<std::string>("dropSeedCluster")),
58  dropRecHits_(iConfig.getParameter<std::string>("dropRecHits")),
59  linkToPackedPF_(iConfig.getParameter<bool>("linkToPackedPFCandidates")),
60  saveNonZSClusterShapes_(iConfig.getParameter<std::string>("saveNonZSClusterShapes")),
61  reducedBarrelRecHitCollectionToken_(consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("reducedBarrelRecHitCollection"))),
62  reducedEndcapRecHitCollectionToken_(consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("reducedEndcapRecHitCollection"))),
63  modifyPhoton_(iConfig.getParameter<bool>("modifyPhotons"))
64 {
66  if( modifyPhoton_ ) {
67  const edm::ParameterSet& mod_config = iConfig.getParameter<edm::ParameterSet>("modifierConfig");
68  photonModifier_.reset(new pat::ObjectModifier<pat::Photon>(mod_config) );
69  photonModifier_->setConsumes(sumes);
70  } else {
71  photonModifier_.reset(nullptr);
72  }
73 
74  produces<std::vector<pat::Photon> >();
75  if (linkToPackedPF_) {
76  reco2pf_ = consumes<edm::ValueMap<std::vector<reco::PFCandidateRef>>>(iConfig.getParameter<edm::InputTag>("recoToPFMap"));
77  pf2pc_ = consumes<edm::Association<pat::PackedCandidateCollection>>(iConfig.getParameter<edm::InputTag>("packedPFCandidates"));
78  pc_ = consumes<pat::PackedCandidateCollection>(iConfig.getParameter<edm::InputTag>("packedPFCandidates"));
79  }
80  mayConsume<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEB"));
81  mayConsume<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEE"));
82 }
T getParameter(std::string const &) const
StringCutObjectSelector< pat::Photon > dropSeedCluster_
StringCutObjectSelector< pat::Photon > dropPreshowerClusters_
edm::EDGetTokenT< edm::View< pat::Photon > > src_
StringCutObjectSelector< pat::Photon > dropSuperClusters_
edm::EDGetTokenT< edm::ValueMap< std::vector< reco::PFCandidateRef > > > reco2pf_
StringCutObjectSelector< pat::Photon > dropRecHits_
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
edm::EDGetTokenT< EcalRecHitCollection > reducedEndcapRecHitCollectionToken_
edm::EDGetTokenT< EcalRecHitCollection > reducedBarrelRecHitCollectionToken_
edm::EDGetTokenT< edm::Association< pat::PackedCandidateCollection > > pf2pc_
edm::EDGetTokenT< pat::PackedCandidateCollection > pc_
StringCutObjectSelector< pat::Photon > saveNonZSClusterShapes_
std::unique_ptr< pat::ObjectModifier< pat::Photon > > photonModifier_
StringCutObjectSelector< pat::Photon > dropBasicClusters_
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.

References modifyPhoton_, and photonModifier_.

85  {
86  if( modifyPhoton_ ) photonModifier_->setEventContent(iSetup);
87 }
std::unique_ptr< pat::ObjectModifier< pat::Photon > > photonModifier_
void pat::PATPhotonSlimmer::produce ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
virtual

Implements edm::EDProducer.

Definition at line 90 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_.

90  {
91  using namespace edm;
92  using namespace std;
93 
95  iEvent.getByToken(src_, src);
96 
100  if (linkToPackedPF_) {
101  iEvent.getByToken(reco2pf_, reco2pf);
102  iEvent.getByToken(pf2pc_, pf2pc);
103  iEvent.getByToken(pc_, pc);
104  }
106 
107  auto_ptr<vector<pat::Photon> > out(new vector<pat::Photon>());
108  out->reserve(src->size());
109 
110  if( modifyPhoton_ ) { photonModifier_->setEvent(iEvent); }
111 
112  std::vector<unsigned int> keys;
113  for (View<pat::Photon>::const_iterator it = src->begin(), ed = src->end(); it != ed; ++it) {
114  out->push_back(*it);
115  pat::Photon & photon = out->back();
116 
117  if( modifyPhoton_ ) { photonModifier_->modify(photon); }
118 
119  if (dropSuperClusters_(photon)) { photon.superCluster_.clear(); photon.embeddedSuperCluster_ = false; }
120  if (dropBasicClusters_(photon)) { photon.basicClusters_.clear(); }
121  if (dropPreshowerClusters_(photon)) { photon.preshowerClusters_.clear(); }
122  if (dropSeedCluster_(photon)) { photon.seedCluster_.clear(); photon.embeddedSeedCluster_ = false; }
123  if (dropRecHits_(photon)) { photon.recHits_ = EcalRecHitCollection(); photon.embeddedRecHits_ = false; }
124 
125  if (linkToPackedPF_) {
126  //std::cout << " PAT photon in " << src.id() << " comes from " << photon.refToOrig_.id() << ", " << photon.refToOrig_.key() << std::endl;
127  keys.clear();
128  for(auto const& pf: (*reco2pf)[photon.refToOrig_]) {
129  if( pf2pc->contains(pf.id()) ) {
130  keys.push_back( (*pf2pc)[pf].key());
131  }
132  }
134  keys.begin(), keys.end());
135  //std::cout << "Photon with pt " << photon.pt() << " associated to " << photon.associatedPackedFCandidateIndices_.size() << " PF Candidates\n";
136  //if there's just one PF Cand then it's me, otherwise I have no univoque parent so my ref will be null
137  if (keys.size() == 1) {
138  photon.refToOrig_ = photon.sourceCandidatePtr(0);
139  } else {
140  photon.refToOrig_ = reco::CandidatePtr(pc.id());
141  }
142  }
143  if (saveNonZSClusterShapes_(photon)) {
144  std::vector<float> vCov = lazyToolsNoZS.localCovariances(*( photon.superCluster()->seed()));
145  float r9 = lazyToolsNoZS.e3x3( *( photon.superCluster()->seed())) / photon.superCluster()->rawEnergy() ;
146  float sigmaIetaIeta = ( !edm::isNotFinite(vCov[0]) ) ? sqrt(vCov[0]) : 0;
147  float sigmaIetaIphi = vCov[1];
148  float sigmaIphiIphi = ( !edm::isNotFinite(vCov[2]) ) ? sqrt(vCov[2]) : 0;
149  float e15o55 = lazyToolsNoZS.e1x5( *( photon.superCluster()->seed()) ) / lazyToolsNoZS.e5x5( *( photon.superCluster()->seed()) );
150  photon.addUserFloat("sigmaIetaIeta_NoZS", sigmaIetaIeta);
151  photon.addUserFloat("sigmaIetaIphi_NoZS", sigmaIetaIphi);
152  photon.addUserFloat("sigmaIphiIphi_NoZS", sigmaIphiIphi);
153  photon.addUserFloat("r9_NoZS", r9);
154  photon.addUserFloat("e1x5_over_e5x5_NoZS", e15o55);
155  }
156  }
157 
158  iEvent.put(out);
159 }
StringCutObjectSelector< pat::Photon > dropSeedCluster_
Analysis-level Photon class.
Definition: Photon.h:47
reco::SuperClusterRef superCluster() const
override the superCluster method from CaloJet, to access the internal storage of the supercluster ...
StringCutObjectSelector< pat::Photon > dropPreshowerClusters_
edm::EDGetTokenT< edm::View< pat::Photon > > src_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
StringCutObjectSelector< pat::Photon > dropSuperClusters_
void setAssociatedPackedPFCandidates(const edm::RefProd< pat::PackedCandidateCollection > &refprod, T beginIndexItr, T endIndexItr)
References to PFCandidates linked to this object (e.g. for isolation vetos or masking before jet recl...
Definition: Photon.h:308
std::vector< reco::SuperCluster > superCluster_
Definition: Photon.h:329
std::vector< reco::CaloCluster > preshowerClusters_
Place to store electron&#39;s preshower clusters internally.
Definition: Photon.h:335
edm::EDGetTokenT< edm::ValueMap< std::vector< reco::PFCandidateRef > > > reco2pf_
std::vector< reco::CaloCluster > seedCluster_
Place to store electron&#39;s seed cluster internally.
Definition: Photon.h:339
StringCutObjectSelector< pat::Photon > dropRecHits_
bool isNotFinite(T x)
Definition: isFinite.h:10
reco::CandidatePtr sourceCandidatePtr(size_type i) const
get the source candidate pointer with index i
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:113
void addUserFloat(const std::string &label, float data)
Set user-defined float.
Definition: PATObject.h:765
T sqrt(T t)
Definition: SSEVec.h:48
edm::EDGetTokenT< EcalRecHitCollection > reducedEndcapRecHitCollectionToken_
edm::EDGetTokenT< EcalRecHitCollection > reducedBarrelRecHitCollectionToken_
bool embeddedSuperCluster_
Definition: Photon.h:328
tuple out
Definition: dbtoconf.py:99
EcalRecHitCollection recHits_
Place to store electron&#39;s RecHits internally (5x5 around seed+ all RecHits)
Definition: Photon.h:343
bool embeddedRecHits_
True if RecHits stored internally.
Definition: Photon.h:341
edm::Ptr< Candidate > CandidatePtr
persistent reference to an object in a collection of Candidate objects
Definition: CandidateFwd.h:25
edm::EDGetTokenT< edm::Association< pat::PackedCandidateCollection > > pf2pc_
edm::EDGetTokenT< pat::PackedCandidateCollection > pc_
StringCutObjectSelector< pat::Photon > saveNonZSClusterShapes_
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:81
bool embeddedSeedCluster_
True if seed cluster is stored internally.
Definition: Photon.h:337
std::unique_ptr< pat::ObjectModifier< pat::Photon > > photonModifier_
StringCutObjectSelector< pat::Photon > dropBasicClusters_
edm::Ptr< reco::Candidate > refToOrig_
Definition: PATObject.h:407
std::vector< reco::CaloCluster > basicClusters_
Place to store electron&#39;s basic clusters internally.
Definition: Photon.h:333
edm::SortedCollection< EcalRecHit > EcalRecHitCollection

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 beginLuminosityBlock(), and 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 beginLuminosityBlock(), and 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().