test
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 Member Functions | Private Attributes
pat::PATPackedCandidateProducer Class Reference
Inheritance diagram for pat::PATPackedCandidateProducer:
edm::global::EDProducer<> edm::global::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

bool candsOrdering (pat::PackedCandidate i, pat::PackedCandidate j) const
 
 PATPackedCandidateProducer (const edm::ParameterSet &)
 
virtual void produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override
 
template<typename T >
std::vector< size_t > sort_indexes (const std::vector< T > &v) const
 
 ~PATPackedCandidateProducer ()
 
- Public Member Functions inherited from edm::global::EDProducer<>
 EDProducer ()=default
 
- Public Member Functions inherited from edm::global::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducerBase ()
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector
< edm::ProductResolverIndex >
const & 
indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription
const &)> 
registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
 
virtual ~ProducerBase () noexcept(false)
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector
< ProductResolverIndexAndSkipBit >
const & 
itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

float calcDxy (float dx, float dy, float phi) const
 
float calcDz (reco::Candidate::Point p, reco::Candidate::Point v, const reco::Candidate &c) const
 

Private Attributes

const edm::EDGetTokenT
< reco::PFCandidateCollection
Cands_
 
const double minPtForTrackProperties_
 
const edm::EDGetTokenT
< std::vector
< reco::PFCandidate > > 
PuppiCands_
 
const edm::EDGetTokenT
< edm::ValueMap
< reco::CandidatePtr > > 
PuppiCandsMap_
 
const edm::EDGetTokenT
< std::vector
< reco::PFCandidate > > 
PuppiCandsNoLep_
 
const edm::EDGetTokenT
< edm::ValueMap< float > > 
PuppiWeight_
 
const edm::EDGetTokenT
< edm::ValueMap< float > > 
PuppiWeightNoLep_
 
const edm::EDGetTokenT
< edm::Association
< reco::VertexCollection > > 
PVAsso_
 
const edm::EDGetTokenT
< edm::ValueMap< int > > 
PVAssoQuality_
 
const edm::EDGetTokenT
< reco::VertexCollection
PVOrigs_
 
const edm::EDGetTokenT
< reco::VertexCollection
PVs_
 
std::vector< edm::EDGetTokenT
< edm::View
< reco::CompositePtrCandidate > > > 
SVWhiteLists_
 
const edm::EDGetTokenT
< reco::TrackCollection
TKOrigs_
 
const bool usePuppi_
 

Additional Inherited Members

- Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
typedef
ProductRegistryHelper::TypeLabelList 
TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- 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

Definition at line 36 of file PATPackedCandidateProducer.cc.

Constructor & Destructor Documentation

pat::PATPackedCandidateProducer::PATPackedCandidateProducer ( const edm::ParameterSet iConfig)
explicit

Definition at line 96 of file PATPackedCandidateProducer.cc.

References edm::EDConsumerBase::consumes(), edm::ParameterSet::getParameter(), and SVWhiteLists_.

96  :
97  usePuppi_(!iConfig.getParameter<edm::InputTag>("PuppiSrc").encode().empty() ||
98  !iConfig.getParameter<edm::InputTag>("PuppiNoLepSrc").encode().empty()),
99  Cands_(consumes<reco::PFCandidateCollection>(iConfig.getParameter<edm::InputTag>("inputCollection"))),
100  PVs_(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("inputVertices"))),
102  PVAssoQuality_(consumes<edm::ValueMap<int> >(iConfig.getParameter<edm::InputTag>("vertexAssociator"))),
103  PVOrigs_(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("originalVertices"))),
104  TKOrigs_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("originalTracks"))),
105  PuppiWeight_(usePuppi_ ? consumes<edm::ValueMap<float> >(iConfig.getParameter<edm::InputTag>("PuppiSrc")) : edm::EDGetTokenT< edm::ValueMap<float> >()),
106  PuppiWeightNoLep_(usePuppi_ ? consumes<edm::ValueMap<float> >(iConfig.getParameter<edm::InputTag>("PuppiNoLepSrc")) : edm::EDGetTokenT< edm::ValueMap<float> >()),
107  PuppiCandsMap_(usePuppi_ ? consumes<edm::ValueMap<reco::CandidatePtr> >(iConfig.getParameter<edm::InputTag>("PuppiSrc")) : edm::EDGetTokenT<edm::ValueMap<reco::CandidatePtr> >() ),
108  PuppiCands_(usePuppi_ ? consumes<std::vector< reco::PFCandidate > >(iConfig.getParameter<edm::InputTag>("PuppiSrc")) : edm::EDGetTokenT<std::vector< reco::PFCandidate > >() ),
109  PuppiCandsNoLep_(usePuppi_ ? consumes<std::vector< reco::PFCandidate > >(iConfig.getParameter<edm::InputTag>("PuppiNoLepSrc")) : edm::EDGetTokenT<std::vector< reco::PFCandidate > >()),
110  minPtForTrackProperties_(iConfig.getParameter<double>("minPtForTrackProperties"))
111 {
112  std::vector<edm::InputTag> sv_tags = iConfig.getParameter<std::vector<edm::InputTag> >("secondaryVerticesForWhiteList");
113  for(auto itag : sv_tags){
114  SVWhiteLists_.push_back(
116  );
117  }
118 
119  produces< std::vector<pat::PackedCandidate> > ();
120  produces< edm::Association<pat::PackedCandidateCollection> > ();
121  produces< edm::Association<reco::PFCandidateCollection> > ();
122 }
T getParameter(std::string const &) const
const edm::EDGetTokenT< edm::ValueMap< float > > PuppiWeight_
const edm::EDGetTokenT< reco::PFCandidateCollection > Cands_
std::string encode() const
Definition: InputTag.cc:164
const edm::EDGetTokenT< edm::ValueMap< int > > PVAssoQuality_
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
const edm::EDGetTokenT< reco::TrackCollection > TKOrigs_
const edm::EDGetTokenT< std::vector< reco::PFCandidate > > PuppiCands_
const edm::EDGetTokenT< edm::ValueMap< reco::CandidatePtr > > PuppiCandsMap_
const edm::EDGetTokenT< edm::ValueMap< float > > PuppiWeightNoLep_
std::vector< edm::EDGetTokenT< edm::View< reco::CompositePtrCandidate > > > SVWhiteLists_
edm::Ptr< Candidate > CandidatePtr
persistent reference to an object in a collection of Candidate objects
Definition: CandidateFwd.h:25
const edm::EDGetTokenT< std::vector< reco::PFCandidate > > PuppiCandsNoLep_
const edm::EDGetTokenT< reco::VertexCollection > PVOrigs_
const edm::EDGetTokenT< edm::Association< reco::VertexCollection > > PVAsso_
const edm::EDGetTokenT< reco::VertexCollection > PVs_
pat::PATPackedCandidateProducer::~PATPackedCandidateProducer ( )

Definition at line 124 of file PATPackedCandidateProducer.cc.

124 {}

Member Function Documentation

float pat::PATPackedCandidateProducer::calcDxy ( float  dx,
float  dy,
float  phi 
) const
inlineprivate

Definition at line 87 of file PATPackedCandidateProducer.cc.

References funct::cos(), and funct::sin().

87  {
88  return - dx * std::sin(phi) + dy * std::cos(phi);
89  }
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
float pat::PATPackedCandidateProducer::calcDz ( reco::Candidate::Point  p,
reco::Candidate::Point  v,
const reco::Candidate c 
) const
inlineprivate

Definition at line 90 of file PATPackedCandidateProducer.cc.

References reco::Candidate::pt(), reco::Candidate::px(), reco::Candidate::py(), and reco::Candidate::pz().

90  {
91  return p.Z()-v.Z() - ((p.X()-v.X()) * c.px() + (p.Y()-v.Y())*c.py()) * c.pz()/(c.pt()*c.pt());
92  }
virtual double pt() const =0
transverse momentum
virtual double pz() const =0
z coordinate of momentum vector
virtual double py() const =0
y coordinate of momentum vector
virtual double px() const =0
x coordinate of momentum vector
bool pat::PATPackedCandidateProducer::candsOrdering ( pat::PackedCandidate  i,
pat::PackedCandidate  j 
) const
inline

Definition at line 44 of file PATPackedCandidateProducer.cc.

References funct::abs(), pat::PackedCandidate::charge(), pat::PackedCandidate::eta(), edm::Ref< C, T, F >::key(), minPtForTrackProperties_, pat::PackedCandidate::pt(), and pat::PackedCandidate::vertexRef().

Referenced by sort_indexes().

44  {
45  if (std::abs(i.charge()) == std::abs(j.charge())) {
46  if(i.charge()!=0){
47  if(i.pt() > minPtForTrackProperties_ and j.pt() <= minPtForTrackProperties_ ) return true;
48  if(i.pt() <= minPtForTrackProperties_ and j.pt() > minPtForTrackProperties_ ) return false;
49  }
50  if(i.vertexRef() == j.vertexRef())
51  return i.eta() > j.eta();
52  else
53  return i.vertexRef().key() < j.vertexRef().key();
54  }
55  return std::abs(i.charge()) > std::abs(j.charge());
56  }
virtual double pt() const
transverse momentum
key_type key() const
Accessor for product key.
Definition: Ref.h:264
virtual double eta() const
momentum pseudorapidity
const reco::VertexRef vertexRef() const
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
virtual int charge() const
electric charge
void pat::PATPackedCandidateProducer::produce ( edm::StreamID  ,
edm::Event iEvent,
const edm::EventSetup iSetup 
) const
overridevirtual

Implements edm::global::EDProducerBase.

Definition at line 128 of file PATPackedCandidateProducer.cc.

References funct::abs(), EnergyCorrector::c, reco::TrackBase::dz(), reco::PFCandidate::ecalEnergy(), edm::helper::Filler< Map >::fill(), edm::Event::getByToken(), reco::PFCandidate::gsfTrackRef(), reco::HitPattern::hasValidHitInPixelLayer(), reco::PFCandidate::hcalEnergy(), reco::TrackBase::highPurity, reco::TrackBase::hitPattern(), i, edm::HandleBase::id(), edm::Ptr< T >::id(), edm::Ref< C, T, F >::id(), cuy::ii, edm::helper::Filler< Map >::insert(), edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), j, edm::Ptr< T >::key(), edm::Ref< C, T, F >::key(), relval_steps::key, reco::HitPattern::MISSING_INNER_HITS, pat::PackedCandidate::moreLostInnerHits, eostools::move(), reco::PFCandidate::muonRef(), pat::PackedCandidate::noLostInnerHits, reco::HitPattern::numberOfLostHits(), pat::PackedCandidate::oneLostInnerHit, reco::LeafCandidate::pdgId(), reco::LeafCandidate::phi(), reco::TrackBase::phi(), GeomDetEnumerators::PixelBarrel, reco::LeafCandidate::polarP4(), edm::Handle< T >::product(), reco::LeafCandidate::pt(), edm::Event::put(), HLT_25ns10e33_v2_cff::quality, pat::qualityMap, reco::TrackBase::referencePoint(), edm::View< T >::size(), reco::PFCandidate::trackRef(), pat::PackedCandidate::UsedInFitTight, pat::PackedCandidate::validHitInFirstPixelBarrelLayer, and reco::PFCandidate::vertex().

128  {
129 
131  iEvent.getByToken( Cands_, cands );
132  std::vector<reco::Candidate>::const_iterator cand;
133 
134  edm::Handle<edm::ValueMap<float> > puppiWeight;
137  edm::Handle<edm::ValueMap<float> > puppiWeightNoLep;
139  std::vector<reco::CandidatePtr> puppiCandsNoLepPtrs;
140  if(usePuppi_){
141  iEvent.getByToken( PuppiWeight_, puppiWeight );
142  iEvent.getByToken( PuppiCandsMap_, puppiCandsMap );
143  iEvent.getByToken( PuppiCands_, puppiCands );
144  iEvent.getByToken( PuppiWeightNoLep_, puppiWeightNoLep );
145  iEvent.getByToken( PuppiCandsNoLep_, puppiCandsNoLep );
146  for (auto pup : *puppiCandsNoLep){
147  puppiCandsNoLepPtrs.push_back(pup.sourceCandidatePtr(0));
148  }
149  }
150  std::vector<int> mappingPuppi(usePuppi_ ? puppiCands->size() : 0);
151 
153  iEvent.getByToken( PVOrigs_, PVOrigs );
154 
156  iEvent.getByToken(PVAsso_,assoHandle);
157  edm::Handle<edm::ValueMap<int> > assoQualityHandle;
158  iEvent.getByToken(PVAssoQuality_,assoQualityHandle);
159  const edm::Association<reco::VertexCollection> & associatedPV=*(assoHandle.product());
160  const edm::ValueMap<int> & associationQuality=*(assoQualityHandle.product());
161 
162 
163  std::set<unsigned int> whiteList;
164  for(auto itoken : SVWhiteLists_) {
166  iEvent.getByToken(itoken, svWhiteListHandle);
167  const edm::View<reco::CompositePtrCandidate > & svWhiteList=*(svWhiteListHandle.product());
168  for(unsigned int i=0; i<svWhiteList.size();i++) {
169  for(unsigned int j=0; j< svWhiteList[i].numberOfSourceCandidatePtrs(); j++) {
170  const edm::Ptr<reco::Candidate> & c = svWhiteList[i].sourceCandidatePtr(j);
171  if(c.id() == cands.id()) whiteList.insert(c.key());
172  }
173  }
174  }
175 
176 
178  iEvent.getByToken( PVs_, PVs );
179  reco::VertexRef PV(PVs.id());
180  reco::VertexRefProd PVRefProd(PVs);
181  math::XYZPoint PVpos;
182 
183 
185  iEvent.getByToken( TKOrigs_, TKOrigs );
186  auto outPtrP = std::make_unique<std::vector<pat::PackedCandidate>>();
187  std::vector<int> mapping(cands->size());
188  std::vector<int> mappingReverse(cands->size());
189  std::vector<int> mappingTk(TKOrigs->size(), -1);
190 
191  for(unsigned int ic=0, nc = cands->size(); ic < nc; ++ic) {
192  const reco::PFCandidate &cand=(*cands)[ic];
193  const reco::Track *ctrack = 0;
194  if ((abs(cand.pdgId()) == 11 || cand.pdgId() == 22) && cand.gsfTrackRef().isNonnull()) {
195  ctrack = &*cand.gsfTrackRef();
196  } else if (cand.trackRef().isNonnull()) {
197  ctrack = &*cand.trackRef();
198  }
199  if (ctrack) {
200  float dist=1e99;
201  int pvi=-1;
202  for(size_t ii=0;ii<PVs->size();ii++){
203  float dz=std::abs(ctrack->dz( ((*PVs)[ii]).position()));
204  if(dz<dist) {pvi=ii;dist=dz; }
205  }
206  PV = reco::VertexRef(PVs, pvi);
207  math::XYZPoint vtx = cand.vertex();
209  const reco::VertexRef & PVOrig = associatedPV[reco::CandidatePtr(cands,ic)];
210  if(PVOrig.isNonnull()) PV = reco::VertexRef(PVs, PVOrig.key()); // WARNING: assume the PV slimmer is keeping same order
211  int quality=associationQuality[reco::CandidatePtr(cands,ic)];
212 // if ((size_t)pvi!=PVOrig.key()) std::cout << "not closest in Z" << pvi << " " << PVOrig.key() << " " << cand.pt() << " " << quality << std::endl;
213  // TrajectoryStateOnSurface tsos = extrapolator.extrapolate(trajectoryStateTransform::initialFreeState(*ctrack,&*magneticField), RecoVertex::convertPos(PV->position()));
214  // vtx = tsos.globalPosition();
215  // phiAtVtx = tsos.globalDirection().phi();
216  vtx = ctrack->referencePoint();
217  float phiAtVtx = ctrack->phi();
219  if (nlost == 0) {
222  }
223  } else {
225  }
226 
227 
228  outPtrP->push_back( pat::PackedCandidate(cand.polarP4(), vtx, phiAtVtx, cand.pdgId(), PVRefProd, PV.key()));
229  outPtrP->back().setAssociationQuality(pat::PackedCandidate::PVAssociationQuality(qualityMap[quality]));
230  if(cand.trackRef().isNonnull() && PVOrig->trackWeight(cand.trackRef()) > 0.5 && quality == 7) {
231  outPtrP->back().setAssociationQuality(pat::PackedCandidate::UsedInFitTight);
232  }
233  // properties of the best track
234  outPtrP->back().setLostInnerHits( lostHits );
235  if(outPtrP->back().pt() > minPtForTrackProperties_ || whiteList.find(ic)!=whiteList.end()) {
236  outPtrP->back().setTrackProperties(*ctrack);
237  //outPtrP->back().setTrackProperties(*ctrack,tsos.curvilinearError());
238  }
239 
240  // these things are always for the CKF track
241  outPtrP->back().setTrackHighPurity( cand.trackRef().isNonnull() && cand.trackRef()->quality(reco::Track::highPurity) );
242  if (cand.muonRef().isNonnull()) {
243  outPtrP->back().setMuonID(cand.muonRef()->isStandAloneMuon(), cand.muonRef()->isGlobalMuon());
244  }
245  } else {
246 
247  if (!PVs->empty()) {
248  PV = reco::VertexRef(PVs, 0);
249  PVpos = PV->position();
250  }
251 
252  outPtrP->push_back( pat::PackedCandidate(cand.polarP4(), PVpos, cand.phi(), cand.pdgId(), PVRefProd, PV.key()));
254  }
255 
256  // neutrals
257 
258  if(abs(cand.pdgId()) == 1 || abs(cand.pdgId()) == 130) {
259  outPtrP->back().setHcalFraction(cand.hcalEnergy()/(cand.ecalEnergy()+cand.hcalEnergy()));
260  } else {
261  outPtrP->back().setHcalFraction(0);
262  }
263 
264 
265  if (usePuppi_){
266  reco::PFCandidateRef pkref( cands, ic );
267  // outPtrP->back().setPuppiWeight( (*puppiWeight)[pkref]);
268 
269  float puppiWeightVal = (*puppiWeight)[pkref];
270  float puppiWeightNoLepVal = 0.0;
271  // Check the "no lepton" puppi weights.
272  // If present, then it is not a lepton, use stored weight
273  // If absent, it is a lepton, so set the weight to 1.0
274  if ( puppiWeightNoLep.isValid() ) {
275  // Look for the pointer inside the "no lepton" candidate collection.
276  auto pkrefPtr = pkref->sourceCandidatePtr(0);
277 
278  bool foundNoLep = false;
279  for ( size_t ipcnl = 0; ipcnl < puppiCandsNoLepPtrs.size(); ipcnl++){
280  if (puppiCandsNoLepPtrs[ipcnl] == pkrefPtr){
281  foundNoLep = true;
282  puppiWeightNoLepVal = puppiCandsNoLep->at(ipcnl).pt()/cand.pt(); // a hack for now, should use the value map
283  break;
284  }
285  }
286  if ( !foundNoLep || puppiWeightNoLepVal > 1 ) {
287  puppiWeightNoLepVal = 1.0;
288  }
289  }
290  outPtrP->back().setPuppiWeight( puppiWeightVal, puppiWeightNoLepVal );
291 
292  mappingPuppi[((*puppiCandsMap)[pkref]).key()]=ic;
293  }
294 
295  mapping[ic] = ic; // trivial at the moment!
296  if (cand.trackRef().isNonnull() && cand.trackRef().id() == TKOrigs.id()) {
297  mappingTk[cand.trackRef().key()] = ic;
298  }
299 
300  }
301 
302  auto outPtrPSorted = std::make_unique<std::vector<pat::PackedCandidate>>();
303  std::vector<size_t> order=sort_indexes(*outPtrP);
304  std::vector<size_t> reverseOrder(order.size());
305  for(size_t i=0,nc=cands->size();i<nc;i++) {
306  outPtrPSorted->push_back((*outPtrP)[order[i]]);
307  reverseOrder[order[i]] = i;
308  mappingReverse[order[i]]=i;
309  }
310 
311  // Fix track association for sorted candidates
312  for(size_t i=0,ntk=mappingTk.size();i<ntk;i++){
313  if(mappingTk[i] >= 0)
314  mappingTk[i]=reverseOrder[mappingTk[i]];
315  }
316 
317  for(size_t i=0,ntk=mappingPuppi.size();i<ntk;i++){
318  mappingPuppi[i]=reverseOrder[mappingPuppi[i]];
319  }
320 
322 
323  // now build the two maps
324  auto pf2pc = std::make_unique<edm::Association<pat::PackedCandidateCollection>>(oh);
325  auto pc2pf = std::make_unique<edm::Association<reco::PFCandidateCollection>>(cands);
328  pf2pcFiller.insert(cands, mappingReverse.begin(), mappingReverse.end());
329  pc2pfFiller.insert(oh , order.begin(), order.end());
330  // include also the mapping track -> packed PFCand
331  pf2pcFiller.insert(TKOrigs, mappingTk.begin(), mappingTk.end());
332  if(usePuppi_) pf2pcFiller.insert(puppiCands, mappingPuppi.begin(), mappingPuppi.end());
333 
334  pf2pcFiller.fill();
335  pc2pfFiller.fill();
336  iEvent.put(std::move(pf2pc));
337  iEvent.put(std::move(pc2pf));
338 
339 }
double ecalEnergy() const
return corrected Ecal energy
Definition: PFCandidate.h:221
const Point & referencePoint() const
Reference point on the track.
Definition: TrackBase.h:676
int i
Definition: DBlmapReader.cc:9
const edm::EDGetTokenT< edm::ValueMap< float > > PuppiWeight_
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
bool hasValidHitInPixelLayer(enum PixelSubdetector::SubDetector, uint16_t layer) const
Definition: HitPattern.cc:327
key_type key() const
Definition: Ptr.h:186
ProductID id() const
Definition: HandleBase.cc:15
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
virtual double phi() const final
momentum azimuthal angle
size_type size() const
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:640
static const int qualityMap[8]
conversion map from quality flags used in PV association and miniAOD one
const edm::EDGetTokenT< reco::PFCandidateCollection > Cands_
key_type key() const
Accessor for product key.
Definition: Ref.h:264
int ii
Definition: cuy.py:588
ProductID id() const
Accessor for product ID.
Definition: Ref.h:258
const edm::EDGetTokenT< edm::ValueMap< int > > PVAssoQuality_
reco::TrackRef trackRef() const
Definition: PFCandidate.cc:438
def move
Definition: eostools.py:510
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
int j
Definition: DBlmapReader.cc:9
virtual const Point & vertex() const
vertex position (overwritten by PF...)
Definition: PFCandidate.cc:652
const edm::EDGetTokenT< reco::TrackCollection > TKOrigs_
const edm::EDGetTokenT< std::vector< reco::PFCandidate > > PuppiCands_
std::vector< size_t > sort_indexes(const std::vector< T > &v) const
const edm::EDGetTokenT< edm::ValueMap< reco::CandidatePtr > > PuppiCandsMap_
bool isValid() const
Definition: HandleBase.h:75
const edm::EDGetTokenT< edm::ValueMap< float > > PuppiWeightNoLep_
edm::Ref< VertexCollection > VertexRef
persistent reference to a Vertex
Definition: VertexFwd.h:13
string key
FastSim: produces sample of signal events, overlayed with premixed minbias events.
reco::MuonRef muonRef() const
Definition: PFCandidate.cc:455
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:604
LostInnerHits
Enumerator specifying the.
std::vector< edm::EDGetTokenT< edm::View< reco::CompositePtrCandidate > > > SVWhiteLists_
edm::Ptr< Candidate > CandidatePtr
persistent reference to an object in a collection of Candidate objects
Definition: CandidateFwd.h:25
T const * product() const
Definition: Handle.h:81
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:445
const edm::EDGetTokenT< std::vector< reco::PFCandidate > > PuppiCandsNoLep_
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
ProductID id() const
Accessor for product ID.
Definition: Ptr.h:181
const edm::EDGetTokenT< reco::VertexCollection > PVOrigs_
virtual const PolarLorentzVector & polarP4() const final
four-momentum Lorentz vector
int numberOfLostHits(HitCategory category) const
Definition: HitPattern.h:902
virtual int pdgId() const final
PDG identifier.
Particle reconstructed by the particle flow algorithm.
Definition: PFCandidate.h:39
reco::GsfTrackRef gsfTrackRef() const
Definition: PFCandidate.cc:476
double hcalEnergy() const
return corrected Hcal energy
Definition: PFCandidate.h:231
const edm::EDGetTokenT< edm::Association< reco::VertexCollection > > PVAsso_
const edm::EDGetTokenT< reco::VertexCollection > PVs_
virtual double pt() const final
transverse momentum
template<typename T >
std::vector<size_t> pat::PATPackedCandidateProducer::sort_indexes ( const std::vector< T > &  v) const
inline

Definition at line 58 of file PATPackedCandidateProducer.cc.

References candsOrdering(), i, and findQualityFiles::v.

58  {
59  std::vector<size_t> idx(v.size());
60  for (size_t i = 0; i != idx.size(); ++i) idx[i] = i;
61  std::sort(idx.begin(), idx.end(),[&v,this](size_t i1, size_t i2) { return candsOrdering(v[i1],v[i2]);});
62  return idx;
63  }
int i
Definition: DBlmapReader.cc:9
bool candsOrdering(pat::PackedCandidate i, pat::PackedCandidate j) const

Member Data Documentation

const edm::EDGetTokenT<reco::PFCandidateCollection> pat::PATPackedCandidateProducer::Cands_
private

Definition at line 71 of file PATPackedCandidateProducer.cc.

const double pat::PATPackedCandidateProducer::minPtForTrackProperties_
private

Definition at line 84 of file PATPackedCandidateProducer.cc.

Referenced by candsOrdering().

const edm::EDGetTokenT<std::vector< reco::PFCandidate > > pat::PATPackedCandidateProducer::PuppiCands_
private

Definition at line 80 of file PATPackedCandidateProducer.cc.

const edm::EDGetTokenT<edm::ValueMap<reco::CandidatePtr> > pat::PATPackedCandidateProducer::PuppiCandsMap_
private

Definition at line 79 of file PATPackedCandidateProducer.cc.

const edm::EDGetTokenT<std::vector< reco::PFCandidate > > pat::PATPackedCandidateProducer::PuppiCandsNoLep_
private

Definition at line 81 of file PATPackedCandidateProducer.cc.

const edm::EDGetTokenT< edm::ValueMap<float> > pat::PATPackedCandidateProducer::PuppiWeight_
private

Definition at line 77 of file PATPackedCandidateProducer.cc.

const edm::EDGetTokenT< edm::ValueMap<float> > pat::PATPackedCandidateProducer::PuppiWeightNoLep_
private

Definition at line 78 of file PATPackedCandidateProducer.cc.

const edm::EDGetTokenT<edm::Association<reco::VertexCollection> > pat::PATPackedCandidateProducer::PVAsso_
private

Definition at line 73 of file PATPackedCandidateProducer.cc.

const edm::EDGetTokenT<edm::ValueMap<int> > pat::PATPackedCandidateProducer::PVAssoQuality_
private

Definition at line 74 of file PATPackedCandidateProducer.cc.

const edm::EDGetTokenT<reco::VertexCollection> pat::PATPackedCandidateProducer::PVOrigs_
private

Definition at line 75 of file PATPackedCandidateProducer.cc.

const edm::EDGetTokenT<reco::VertexCollection> pat::PATPackedCandidateProducer::PVs_
private

Definition at line 72 of file PATPackedCandidateProducer.cc.

std::vector< edm::EDGetTokenT<edm::View<reco::CompositePtrCandidate> > > pat::PATPackedCandidateProducer::SVWhiteLists_
private

Definition at line 82 of file PATPackedCandidateProducer.cc.

Referenced by PATPackedCandidateProducer().

const edm::EDGetTokenT<reco::TrackCollection> pat::PATPackedCandidateProducer::TKOrigs_
private

Definition at line 76 of file PATPackedCandidateProducer.cc.

const bool pat::PATPackedCandidateProducer::usePuppi_
private

Definition at line 69 of file PATPackedCandidateProducer.cc.