CMS 3D CMS Logo
PATLostTracks.cc
Go to the documentation of this file.
1 #include "DataFormats/RecoCandidate/interface/RecoChargedCandidate.h"
2 #include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
3 #include "DataFormats/ParticleFlowCandidate/interface/PFCandidateFwd.h"
4 #include "DataFormats/PatCandidates/interface/PackedCandidate.h"
5 #include "DataFormats/Common/interface/Association.h"
6 #include "DataFormats/MuonReco/interface/Muon.h"
7 #include "DataFormats/VertexReco/interface/Vertex.h"
8 #include "DataFormats/VertexReco/interface/VertexFwd.h"
9 #include "FWCore/Framework/interface/global/EDProducer.h"
10 #include "FWCore/Framework/interface/ESHandle.h"
11 #include "FWCore/Framework/interface/Event.h"
12 #include "FWCore/ParameterSet/interface/ParameterSet.h"
13 #include "CommonTools/Utils/interface/StringCutObjectSelector.h"
14 
15 namespace {
16  bool passesQuality(const reco::Track& trk, const std::vector<reco::TrackBase::TrackQuality>& allowedQuals) {
17  for (const auto& qual : allowedQuals) {
18  if (trk.quality(qual))
19  return true;
20  }
21  return false;
22  }
23 } // namespace
24 
25 namespace pat {
26  class PATLostTracks : public edm::global::EDProducer<> {
27  public:
28  explicit PATLostTracks(const edm::ParameterSet&);
29  ~PATLostTracks() override;
30 
31  void produce(edm::StreamID, edm::Event&, const edm::EventSetup&) const override;
32 
33  private:
34  enum class TrkStatus { NOTUSED = 0, PFCAND, PFCANDNOTRKPROPS, PFELECTRON, PFPOSITRON, VTX };
35  bool passTrkCuts(const reco::Track& tr) const;
36  void addPackedCandidate(std::vector<pat::PackedCandidate>& cands,
37  const reco::TrackRef& trk,
38  const reco::VertexRef& pvSlimmed,
39  const reco::VertexRefProd& pvSlimmedColl,
40  const TrkStatus& trkStatus,
41  const pat::PackedCandidate::PVAssociationQuality& pvAssocQuality,
42  edm::Handle<reco::MuonCollection> muons) const;
43  std::pair<int, pat::PackedCandidate::PVAssociationQuality> associateTrkToVtx(const reco::VertexCollection& vertices,
44  const reco::TrackRef& trk) const;
45 
46  private:
47  const edm::EDGetTokenT<reco::PFCandidateCollection> cands_;
48  const edm::EDGetTokenT<edm::Association<pat::PackedCandidateCollection>> map_;
49  const edm::EDGetTokenT<reco::TrackCollection> tracks_;
50  const edm::EDGetTokenT<reco::VertexCollection> vertices_;
51  const edm::EDGetTokenT<reco::VertexCompositeCandidateCollection> kshorts_;
52  const edm::EDGetTokenT<reco::VertexCompositeCandidateCollection> lambdas_;
53  const edm::EDGetTokenT<reco::VertexCollection> pv_;
54  const edm::EDGetTokenT<reco::VertexCollection> pvOrigs_;
55  const double minPt_;
56  const double minHits_;
57  const double minPixelHits_;
58  const double minPtToStoreProps_;
59  const double minPtToStoreLowQualityProps_;
60  const int covarianceVersion_;
61  const std::vector<int> covariancePackingSchemas_;
62  std::vector<reco::TrackBase::TrackQuality> qualsToAutoAccept_;
63  const edm::EDGetTokenT<reco::MuonCollection> muons_;
64  StringCutObjectSelector<reco::Track, false> passThroughCut_;
65  const double maxDzForPrimaryAssignment_;
66  const double maxDzSigForPrimaryAssignment_;
67  const double maxDzErrorForPrimaryAssignment_;
68  const double maxDxyForNotReconstructedPrimary_;
69  const double maxDxySigForNotReconstructedPrimary_;
70  const bool useLegacySetup_;
71  const bool xiSelection_;
72  const double xiMassCut_;
73  };
74 } // namespace pat
75 
76 pat::PATLostTracks::PATLostTracks(const edm::ParameterSet& iConfig)
77  : cands_(consumes<reco::PFCandidateCollection>(iConfig.getParameter<edm::InputTag>("inputCandidates"))),
78  map_(consumes<edm::Association<pat::PackedCandidateCollection>>(
79  iConfig.getParameter<edm::InputTag>("packedPFCandidates"))),
80  tracks_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("inputTracks"))),
81  vertices_(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("secondaryVertices"))),
82  kshorts_(consumes<reco::VertexCompositeCandidateCollection>(iConfig.getParameter<edm::InputTag>("kshorts"))),
83  lambdas_(consumes<reco::VertexCompositeCandidateCollection>(iConfig.getParameter<edm::InputTag>("lambdas"))),
84  pv_(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("primaryVertices"))),
85  pvOrigs_(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("originalVertices"))),
86  minPt_(iConfig.getParameter<double>("minPt")),
87  minHits_(iConfig.getParameter<uint32_t>("minHits")),
88  minPixelHits_(iConfig.getParameter<uint32_t>("minPixelHits")),
89  minPtToStoreProps_(iConfig.getParameter<double>("minPtToStoreProps")),
90  minPtToStoreLowQualityProps_(iConfig.getParameter<double>("minPtToStoreLowQualityProps")),
91  covarianceVersion_(iConfig.getParameter<int>("covarianceVersion")),
92  covariancePackingSchemas_(iConfig.getParameter<std::vector<int>>("covariancePackingSchemas")),
93  muons_(consumes<reco::MuonCollection>(iConfig.getParameter<edm::InputTag>("muons"))),
94  passThroughCut_(iConfig.getParameter<std::string>("passThroughCut")),
95  maxDzForPrimaryAssignment_(
96  iConfig.getParameter<edm::ParameterSet>("pvAssignment").getParameter<double>("maxDzForPrimaryAssignment")),
97  maxDzSigForPrimaryAssignment_(
98  iConfig.getParameter<edm::ParameterSet>("pvAssignment").getParameter<double>("maxDzSigForPrimaryAssignment")),
99  maxDzErrorForPrimaryAssignment_(iConfig.getParameter<edm::ParameterSet>("pvAssignment")
100  .getParameter<double>("maxDzErrorForPrimaryAssignment")),
101  maxDxyForNotReconstructedPrimary_(iConfig.getParameter<edm::ParameterSet>("pvAssignment")
102  .getParameter<double>("maxDxyForNotReconstructedPrimary")),
103  maxDxySigForNotReconstructedPrimary_(iConfig.getParameter<edm::ParameterSet>("pvAssignment")
104  .getParameter<double>("maxDxySigForNotReconstructedPrimary")),
105  useLegacySetup_(iConfig.getParameter<bool>("useLegacySetup")),
106  xiSelection_(iConfig.getParameter<bool>("xiSelection")),
107  xiMassCut_(iConfig.getParameter<double>("xiMassCut")) {
108  std::vector<std::string> trkQuals(iConfig.getParameter<std::vector<std::string>>("qualsToAutoAccept"));
109  std::transform(
110  trkQuals.begin(), trkQuals.end(), std::back_inserter(qualsToAutoAccept_), reco::TrackBase::qualityByName);
111 
112  if (std::find(qualsToAutoAccept_.begin(), qualsToAutoAccept_.end(), reco::TrackBase::undefQuality) !=
113  qualsToAutoAccept_.end()) {
114  std::ostringstream msg;
115  msg << " PATLostTracks has a quality requirement which resolves to undefQuality. This usually means a typo and is "
116  "therefore treated a config error\nquality requirements:\n ";
117  for (const auto& trkQual : trkQuals)
118  msg << trkQual << " ";
119  throw cms::Exception("Configuration") << msg.str();
120  }
121 
122  produces<std::vector<reco::Track>>();
123  produces<std::vector<pat::PackedCandidate>>();
124  produces<std::vector<pat::PackedCandidate>>("eleTracks");
125  produces<edm::Association<pat::PackedCandidateCollection>>();
126 }
127 
128 pat::PATLostTracks::~PATLostTracks() {}
129 
130 void pat::PATLostTracks::produce(edm::StreamID, edm::Event& iEvent, const edm::EventSetup& iSetup) const {
131  edm::Handle<reco::PFCandidateCollection> cands;
132  iEvent.getByToken(cands_, cands);
133 
134  edm::Handle<edm::Association<pat::PackedCandidateCollection>> pf2pc;
135  iEvent.getByToken(map_, pf2pc);
136 
137  edm::Handle<reco::TrackCollection> tracks;
138  iEvent.getByToken(tracks_, tracks);
139 
140  edm::Handle<reco::VertexCollection> vertices;
141  iEvent.getByToken(vertices_, vertices);
142 
143  edm::Handle<reco::MuonCollection> muons;
144  iEvent.getByToken(muons_, muons);
145 
146  edm::Handle<reco::VertexCompositeCandidateCollection> kshorts;
147  iEvent.getByToken(kshorts_, kshorts);
148  edm::Handle<reco::VertexCompositeCandidateCollection> lambdas;
149  iEvent.getByToken(lambdas_, lambdas);
150 
151  edm::Handle<reco::VertexCollection> pvs;
152  iEvent.getByToken(pv_, pvs);
153  reco::VertexRef pv(pvs.id());
154  reco::VertexRefProd pvRefProd(pvs);
155  edm::Handle<reco::VertexCollection> pvOrigs;
156  iEvent.getByToken(pvOrigs_, pvOrigs);
157 
158  auto outPtrTrks = std::make_unique<std::vector<reco::Track>>();
159  auto outPtrTrksAsCands = std::make_unique<std::vector<pat::PackedCandidate>>();
160  auto outPtrEleTrksAsCands = std::make_unique<std::vector<pat::PackedCandidate>>();
161 
162  std::vector<TrkStatus> trkStatus(tracks->size(), TrkStatus::NOTUSED);
163  //Mark all tracks used in candidates
164  //check if packed candidates are storing the tracks by seeing if number of hits >0
165  //currently we dont use that information though
166  //electrons will never store their track (they store the GSF track)
167  for (unsigned int ic = 0, nc = cands->size(); ic < nc; ++ic) {
168  edm::Ref<reco::PFCandidateCollection> r(cands, ic);
169  const reco::PFCandidate& cand = (*cands)[ic];
170  if (cand.charge() && cand.trackRef().isNonnull() && cand.trackRef().id() == tracks.id()) {
171  if (cand.pdgId() == 11)
172  trkStatus[cand.trackRef().key()] = TrkStatus::PFELECTRON;
173  else if (cand.pdgId() == -11)
174  trkStatus[cand.trackRef().key()] = TrkStatus::PFPOSITRON;
175  else if ((*pf2pc)[r]->numberOfHits() > 0)
176  trkStatus[cand.trackRef().key()] = TrkStatus::PFCAND;
177  else
178  trkStatus[cand.trackRef().key()] = TrkStatus::PFCANDNOTRKPROPS;
179  }
180  }
181 
182  //Mark all tracks used in secondary vertices
183  for (const auto& secVert : *vertices) {
184  for (auto trkIt = secVert.tracks_begin(); trkIt != secVert.tracks_end(); trkIt++) {
185  if (trkStatus[trkIt->key()] == TrkStatus::NOTUSED)
186  trkStatus[trkIt->key()] = TrkStatus::VTX;
187  }
188  }
189  for (const auto& v0 : *kshorts) {
190  for (size_t dIdx = 0; dIdx < v0.numberOfDaughters(); dIdx++) {
191  size_t key = (dynamic_cast<const reco::RecoChargedCandidate*>(v0.daughter(dIdx)))->track().key();
192  if (trkStatus[key] == TrkStatus::NOTUSED)
193  trkStatus[key] = TrkStatus::VTX;
194  }
195  }
196  for (const auto& v0 : *lambdas) {
197  double protonCharge = 0;
198  for (size_t dIdx = 0; dIdx < v0.numberOfDaughters(); dIdx++) {
199  size_t key = (dynamic_cast<const reco::RecoChargedCandidate*>(v0.daughter(dIdx)))->track().key();
200  if (trkStatus[key] == TrkStatus::NOTUSED)
201  trkStatus[key] = TrkStatus::VTX;
202  protonCharge += v0.daughter(dIdx)->charge() * v0.daughter(dIdx)->momentum().mag2();
203  }
204  if (xiSelection_) {
205  // selecting potential Xi- -> Lambda pi candidates
206  math::XYZTLorentzVector p4Lambda(v0.px(), v0.py(), v0.pz(), sqrt(v0.momentum().mag2() + v0.mass() * v0.mass()));
207 
208  for (unsigned int trkIndx = 0; trkIndx < tracks->size(); trkIndx++) {
209  reco::TrackRef trk(tracks, trkIndx);
210  if ((*trk).charge() * protonCharge < 0 || trkStatus[trkIndx] != TrkStatus::NOTUSED)
211  continue;
212 
213  math::XYZTLorentzVector p4pi(
214  trk->px(), trk->py(), trk->pz(), sqrt(trk->momentum().mag2() + 0.01947995518)); // pion mass ^2
215  if ((p4Lambda + p4pi).M() < xiMassCut_)
216  trkStatus[trkIndx] = TrkStatus::VTX; // selecting potential Xi- candidates
217  }
218  }
219  }
220  std::vector<int> mapping(tracks->size(), -1);
221  int lostTrkIndx = 0;
222  for (unsigned int trkIndx = 0; trkIndx < tracks->size(); trkIndx++) {
223  reco::TrackRef trk(tracks, trkIndx);
224  if (trkStatus[trkIndx] == TrkStatus::VTX || (trkStatus[trkIndx] == TrkStatus::NOTUSED && passTrkCuts(*trk))) {
225  outPtrTrks->emplace_back(*trk);
226  //association to PV
227  std::pair<int, pat::PackedCandidate::PVAssociationQuality> pvAsso = associateTrkToVtx(*pvOrigs, trk);
228  const reco::VertexRef& pvOrigRef = reco::VertexRef(pvOrigs, pvAsso.first);
229  if (pvOrigRef.isNonnull()) {
230  pv = reco::VertexRef(pvs, pvOrigRef.key()); // WARNING: assume the PV slimmer is keeping same order
231  } else if (!pvs->empty()) {
232  pv = reco::VertexRef(pvs, 0);
233  }
234  addPackedCandidate(*outPtrTrksAsCands, trk, pv, pvRefProd, trkStatus[trkIndx], pvAsso.second, muons);
235 
236  //for creating the reco::Track -> pat::PackedCandidate map
237  //not done for the lostTrack:eleTracks collection
238  mapping[trkIndx] = lostTrkIndx;
239  lostTrkIndx++;
240  } else if ((trkStatus[trkIndx] == TrkStatus::PFELECTRON || trkStatus[trkIndx] == TrkStatus::PFPOSITRON) &&
241  passTrkCuts(*trk)) {
242  //association to PV
243  std::pair<int, pat::PackedCandidate::PVAssociationQuality> pvAsso = associateTrkToVtx(*pvOrigs, trk);
244  const reco::VertexRef& pvOrigRef = reco::VertexRef(pvOrigs, pvAsso.first);
245  if (pvOrigRef.isNonnull()) {
246  pv = reco::VertexRef(pvs, pvOrigRef.key()); // WARNING: assume the PV slimmer is keeping same order
247  } else if (!pvs->empty()) {
248  pv = reco::VertexRef(pvs, 0);
249  }
250  addPackedCandidate(*outPtrEleTrksAsCands, trk, pv, pvRefProd, trkStatus[trkIndx], pvAsso.second, muons);
251  }
252  }
253 
254  iEvent.put(std::move(outPtrTrks));
255  iEvent.put(std::move(outPtrEleTrksAsCands), "eleTracks");
256  edm::OrphanHandle<pat::PackedCandidateCollection> oh = iEvent.put(std::move(outPtrTrksAsCands));
257  auto tk2pc = std::make_unique<edm::Association<pat::PackedCandidateCollection>>(oh);
258  edm::Association<pat::PackedCandidateCollection>::Filler tk2pcFiller(*tk2pc);
259  tk2pcFiller.insert(tracks, mapping.begin(), mapping.end());
260  tk2pcFiller.fill();
261  iEvent.put(std::move(tk2pc));
262 }
263 
264 bool pat::PATLostTracks::passTrkCuts(const reco::Track& tr) const {
265  const bool passTrkHits = tr.pt() > minPt_ && tr.numberOfValidHits() >= minHits_ &&
266  tr.hitPattern().numberOfValidPixelHits() >= minPixelHits_;
267  const bool passTrkQual = passesQuality(tr, qualsToAutoAccept_);
268 
269  return passTrkHits || passTrkQual || passThroughCut_(tr);
270 }
271 
272 void pat::PATLostTracks::addPackedCandidate(std::vector<pat::PackedCandidate>& cands,
273  const reco::TrackRef& trk,
274  const reco::VertexRef& pvSlimmed,
275  const reco::VertexRefProd& pvSlimmedColl,
276  const pat::PATLostTracks::TrkStatus& trkStatus,
277  const pat::PackedCandidate::PVAssociationQuality& pvAssocQuality,
278  edm::Handle<reco::MuonCollection> muons) const {
279  const float mass = 0.13957018;
280 
281  int id = 211 * trk->charge();
282  if (trkStatus == TrkStatus::PFELECTRON)
283  id = 11;
284  else if (trkStatus == TrkStatus::PFPOSITRON)
285  id = -11;
286 
287  // assign the proper pdgId for tracks that are reconstructed as a muon
288  for (auto& mu : *muons) {
289  if (reco::TrackRef(mu.innerTrack()) == trk) {
290  id = -13 * trk->charge();
291  break;
292  }
293  }
294 
295  pat::PackedCandidate::LostInnerHits lostHits = pat::PackedCandidate::noLostInnerHits;
296  int nlost = trk->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS);
297  if (nlost == 0) {
298  if (trk->hitPattern().hasValidHitInPixelLayer(PixelSubdetector::SubDetector::PixelBarrel, 1)) {
299  lostHits = pat::PackedCandidate::validHitInFirstPixelBarrelLayer;
300  }
301  } else {
302  lostHits = (nlost == 1 ? pat::PackedCandidate::oneLostInnerHit : pat::PackedCandidate::moreLostInnerHits);
303  }
304 
305  reco::Candidate::PolarLorentzVector p4(trk->pt(), trk->eta(), trk->phi(), mass);
306  cands.emplace_back(
307  pat::PackedCandidate(p4, trk->vertex(), trk->pt(), trk->eta(), trk->phi(), id, pvSlimmedColl, pvSlimmed.key()));
308 
309  cands.back().setTrackHighPurity(trk->quality(reco::TrackBase::highPurity));
310  cands.back().setCovarianceVersion(covarianceVersion_);
311  cands.back().setLostInnerHits(lostHits);
312  if (trk->pt() > minPtToStoreProps_ || trkStatus == TrkStatus::VTX) {
313  if (useLegacySetup_ || std::abs(id) == 11 || trkStatus == TrkStatus::VTX) {
314  cands.back().setTrackProperties(*trk, covariancePackingSchemas_[4], covarianceVersion_);
315  } else {
316  if (trk->hitPattern().numberOfValidPixelHits() > 0) {
317  cands.back().setTrackProperties(
318  *trk, covariancePackingSchemas_[0], covarianceVersion_); // high quality with pixels
319  } else {
320  cands.back().setTrackProperties(
321  *trk, covariancePackingSchemas_[1], covarianceVersion_); // high quality without pixels
322  }
323  }
324  } else if (!useLegacySetup_ && trk->pt() > minPtToStoreLowQualityProps_) {
325  if (trk->hitPattern().numberOfValidPixelHits() > 0) {
326  cands.back().setTrackProperties(
327  *trk, covariancePackingSchemas_[2], covarianceVersion_); // low quality with pixels
328  } else {
329  cands.back().setTrackProperties(
330  *trk, covariancePackingSchemas_[3], covarianceVersion_); // low quality without pixels
331  }
332  }
333  cands.back().setAssociationQuality(pvAssocQuality);
334 }
335 
336 std::pair<int, pat::PackedCandidate::PVAssociationQuality> pat::PATLostTracks::associateTrkToVtx(
337  const reco::VertexCollection& vertices, const reco::TrackRef& trk) const {
338  //For legacy setup check only if the track is used in fit of the PV, i.e. vertices[0],
339  //and associate quality if weight > 0.5. Otherwise return invalid vertex index (-1)
340  //and default quality flag (NotReconstructedPrimary = 0)
341  if (useLegacySetup_) {
342  float w = vertices[0].trackWeight(trk);
343  if (w > 0.5) {
344  return std::pair<int, pat::PackedCandidate::PVAssociationQuality>(0, pat::PackedCandidate::UsedInFitTight);
345  } else if (w > 0.) {
346  return std::pair<int, pat::PackedCandidate::PVAssociationQuality>(0,
347  pat::PackedCandidate::NotReconstructedPrimary);
348  } else {
349  return std::pair<int, pat::PackedCandidate::PVAssociationQuality>(-1,
350  pat::PackedCandidate::NotReconstructedPrimary);
351  }
352  }
353 
354  //Inspired by CommonTools/RecoAlgos/interface/PrimaryVertexAssignment.h
355  //but without specific association for secondaries in jets and option to use timing
356 
357  int iVtxMaxWeight = -1;
358  int iVtxMinDzDist = -1;
359  size_t idx = 0;
360  float maxWeight = 0;
361  double minDz = std::numeric_limits<double>::max();
362  double minDzSig = std::numeric_limits<double>::max();
363  for (auto const& vtx : vertices) {
364  float w = vtx.trackWeight(trk);
365  double dz = std::abs(trk->dz(vtx.position()));
366  double dzSig = dz / trk->dzError();
367  if (w > maxWeight) {
368  maxWeight = w;
369  iVtxMaxWeight = idx;
370  }
371  if (dzSig < minDzSig) {
372  minDzSig = dzSig;
373  minDz = dz;
374  iVtxMinDzDist = idx;
375  }
376  idx++;
377  }
378  // vertex in which fit the track was used
379  if (iVtxMaxWeight >= 0) {
380  if (maxWeight > 0.5) {
381  return std::pair<int, pat::PackedCandidate::PVAssociationQuality>(iVtxMaxWeight,
382  pat::PackedCandidate::UsedInFitTight);
383  } else {
384  return std::pair<int, pat::PackedCandidate::PVAssociationQuality>(iVtxMaxWeight,
385  pat::PackedCandidate::UsedInFitLoose);
386  }
387  }
388  // vertex "closest in Z" with tight cuts (targetting primary particles)
389  if (minDz < maxDzForPrimaryAssignment_) {
390  const double add_cov = vertices[iVtxMinDzDist].covariance(2, 2);
391  const double dzErr = sqrt(trk->dzError() * trk->dzError() + add_cov);
392  if (minDz / dzErr < maxDzSigForPrimaryAssignment_ && trk->dzError() < maxDzErrorForPrimaryAssignment_) {
393  return std::pair<int, pat::PackedCandidate::PVAssociationQuality>(iVtxMinDzDist,
394  pat::PackedCandidate::CompatibilityDz);
395  }
396  }
397  // if the track is not compatible with other PVs but is compatible with the BeamSpot, we may simply have not reco'ed the PV!
398  // we still point it to the closest in Z, but flag it as possible orphan-primary
399  if (!vertices.empty() && std::abs(trk->dxy(vertices[0].position())) < maxDxyForNotReconstructedPrimary_ &&
400  std::abs(trk->dxy(vertices[0].position()) / trk->dxyError()) < maxDxySigForNotReconstructedPrimary_)
401  return std::pair<int, pat::PackedCandidate::PVAssociationQuality>(iVtxMinDzDist,
402  pat::PackedCandidate::NotReconstructedPrimary);
403  // for tracks not associated to any PV return the closest in dz
404  return std::pair<int, pat::PackedCandidate::PVAssociationQuality>(iVtxMinDzDist, pat::PackedCandidate::OtherDeltaZ);
405 }
406 
407 using pat::PATLostTracks;
408 #include "FWCore/Framework/interface/MakerMacros.h"
409 DEFINE_FWK_MODULE(PATLostTracks);
edm::RefProd< VertexCollection >
electrons_cff.dzErr
dzErr
Definition: electrons_cff.py:354
edm::StreamID
Definition: StreamID.h:30
PDWG_BPHSkim_cff.muons
muons
Definition: PDWG_BPHSkim_cff.py:47
electrons_cff.bool
bool
Definition: electrons_cff.py:366
V0Monitor_cff.v0
v0
Definition: V0Monitor_cff.py:7
Muon.h
edm::Association::Filler
Definition: Association.h:78
pat::PATLostTracks::TrkStatus::PFCAND
pat::PATLostTracks::maxDzForPrimaryAssignment_
const double maxDzForPrimaryAssignment_
Definition: PATLostTracks.cc:65
pat::PATLostTracks::minPtToStoreProps_
const double minPtToStoreProps_
Definition: PATLostTracks.cc:58
pat::PATLostTracks::tracks_
const edm::EDGetTokenT< reco::TrackCollection > tracks_
Definition: PATLostTracks.cc:49
ESHandle.h
PFCandidate.h
amptDefaultParameters_cff.mu
mu
Definition: amptDefaultParameters_cff.py:16
pat::PATLostTracks::TrkStatus::PFPOSITRON
pat::PackedCandidate::oneLostInnerHit
Definition: PackedCandidate.h:803
edm::EDGetTokenT< reco::PFCandidateCollection >
edm
HLT enums.
Definition: AlignableModifier.h:19
reco::TrackBase::undefQuality
Definition: TrackBase.h:151
reco::Candidate::PolarLorentzVector
math::PtEtaPhiMLorentzVector PolarLorentzVector
Lorentz vector.
Definition: Candidate.h:38
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
lostTracks_cfi.lambdas
lambdas
Definition: lostTracks_cfi.py:11
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89301
pat::PATLostTracks::pv_
const edm::EDGetTokenT< reco::VertexCollection > pv_
Definition: PATLostTracks.cc:53
pat::PackedCandidate::moreLostInnerHits
Definition: PackedCandidate.h:804
pat::PATLostTracks::PATLostTracks
PATLostTracks(const edm::ParameterSet &)
Definition: PATLostTracks.cc:76
pat::PATLostTracks::passThroughCut_
StringCutObjectSelector< reco::Track, false > passThroughCut_
Definition: PATLostTracks.cc:64
pat::PATLostTracks::xiMassCut_
const double xiMassCut_
Definition: PATLostTracks.cc:72
mps_check.msg
tuple msg
Definition: mps_check.py:285
pat::PATLostTracks::addPackedCandidate
void addPackedCandidate(std::vector< pat::PackedCandidate > &cands, const reco::TrackRef &trk, const reco::VertexRef &pvSlimmed, const reco::VertexRefProd &pvSlimmedColl, const TrkStatus &trkStatus, const pat::PackedCandidate::PVAssociationQuality &pvAssocQuality, edm::Handle< reco::MuonCollection > muons) const
Definition: PATLostTracks.cc:272
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:46
pat::PATLostTracks::covariancePackingSchemas_
const std::vector< int > covariancePackingSchemas_
Definition: PATLostTracks.cc:61
Association.h
pat::PATLostTracks::covarianceVersion_
const int covarianceVersion_
Definition: PATLostTracks.cc:60
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
edm::Handle< reco::MuonCollection >
pat::PATLostTracks::maxDxyForNotReconstructedPrimary_
const double maxDxyForNotReconstructedPrimary_
Definition: PATLostTracks.cc:68
reco::TrackBase::numberOfValidHits
unsigned short numberOfValidHits() const
number of valid hits found
Definition: TrackBase.h:798
pat::PATLostTracks
Definition: PATLostTracks.cc:26
pat::PATLostTracks::minHits_
const double minHits_
Definition: PATLostTracks.cc:56
AlignmentTracksFromVertexSelector_cfi.vertices
vertices
Definition: AlignmentTracksFromVertexSelector_cfi.py:5
pat::PackedCandidate::UsedInFitLoose
Definition: PackedCandidate.h:727
pat::PATLostTracks::TrkStatus
TrkStatus
Definition: PATLostTracks.cc:34
edm::Ref< TrackCollection >
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
pat::PackedCandidate::validHitInFirstPixelBarrelLayer
Definition: PackedCandidate.h:800
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:637
GeomDetEnumerators::PixelBarrel
Definition: GeomDetEnumerators.h:11
MakerMacros.h
pat::PATLostTracks::maxDxySigForNotReconstructedPrimary_
const double maxDxySigForNotReconstructedPrimary_
Definition: PATLostTracks.cc:69
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
pat::PATLostTracks::maxDzSigForPrimaryAssignment_
const double maxDzSigForPrimaryAssignment_
Definition: PATLostTracks.cc:66
w
const double w
Definition: UKUtility.cc:23
pat::PATLostTracks::maxDzErrorForPrimaryAssignment_
const double maxDzErrorForPrimaryAssignment_
Definition: PATLostTracks.cc:67
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
pat::PATLostTracks::TrkStatus::VTX
reco::Track
Definition: Track.h:27
pat::PATLostTracks::minPixelHits_
const double minPixelHits_
Definition: PATLostTracks.cc:57
pat::PackedCandidate::PVAssociationQuality
PVAssociationQuality
Definition: PackedCandidate.h:722
pat::MuonCollection
std::vector< Muon > MuonCollection
Definition: Muon.h:35
pat::PATLostTracks::xiSelection_
const bool xiSelection_
Definition: PATLostTracks.cc:71
HcalDetIdTransform::transform
unsigned transform(const HcalDetId &id, unsigned transformCode)
Definition: HcalDetIdTransform.cc:7
pat::PATLostTracks::minPt_
const double minPt_
Definition: PATLostTracks.cc:55
edm::global::EDProducer
Definition: EDProducer.h:32
pat::PackedCandidate::noLostInnerHits
Definition: PackedCandidate.h:801
pat::PATLostTracks::produce
void produce(edm::StreamID, edm::Event &, const edm::EventSetup &) const override
Definition: PATLostTracks.cc:130
pat::PATLostTracks::muons_
const edm::EDGetTokenT< reco::MuonCollection > muons_
Definition: PATLostTracks.cc:63
HLT_FULL_cff.cands
cands
Definition: HLT_FULL_cff.py:15161
pat::PackedCandidate::CompatibilityDz
Definition: PackedCandidate.h:726
edm::ParameterSet
Definition: ParameterSet.h:47
pat::PackedCandidate::OtherDeltaZ
Definition: PackedCandidate.h:724
Event.h
tracks
const uint32_t *__restrict__ const HitContainer *__restrict__ TkSoA *__restrict__ tracks
Definition: CAHitNtupletGeneratorKernelsImpl.h:176
ParameterSet
Definition: Functions.h:16
pat::PackedCandidate::LostInnerHits
LostInnerHits
Enumerator specifying the.
Definition: PackedCandidate.h:799
SiStripPI::max
Definition: SiStripPayloadInspectorHelper.h:169
pat::PackedCandidate
Definition: PackedCandidate.h:22
PackedCandidate.h
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:238
reco::VertexRef
edm::Ref< VertexCollection > VertexRef
persistent reference to a Vertex
Definition: VertexFwd.h:13
lostTracks_cfi.kshorts
kshorts
Definition: lostTracks_cfi.py:10
cand
Definition: decayParser.h:32
pat::PATLostTracks::qualsToAutoAccept_
std::vector< reco::TrackBase::TrackQuality > qualsToAutoAccept_
Definition: PATLostTracks.cc:62
createfilelist.int
int
Definition: createfilelist.py:10
MetAnalyzer.pv
def pv(vc)
Definition: MetAnalyzer.py:7
iEvent
int iEvent
Definition: GenABIO.cc:224
pat::PATLostTracks::pvOrigs_
const edm::EDGetTokenT< reco::VertexCollection > pvOrigs_
Definition: PATLostTracks.cc:54
p4
double p4[4]
Definition: TauolaWrapper.h:92
trackerHitRTTI::vector
Definition: trackerHitRTTI.h:21
edm::EventSetup
Definition: EventSetup.h:58
pat
Definition: HeavyIon.h:7
FSQDQM_cfi.pvs
pvs
Definition: FSQDQM_cfi.py:12
l1t::PFCandidateCollection
std::vector< l1t::PFCandidate > PFCandidateCollection
Definition: PFCandidate.h:57
pat::PATLostTracks::lambdas_
const edm::EDGetTokenT< reco::VertexCompositeCandidateCollection > lambdas_
Definition: PATLostTracks.cc:52
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
pat::PackedCandidateCollection
std::vector< pat::PackedCandidate > PackedCandidateCollection
Definition: PackedCandidate.h:1131
alignCSCRings.r
r
Definition: alignCSCRings.py:93
reco::VertexCompositeCandidateCollection
std::vector< VertexCompositeCandidate > VertexCompositeCandidateCollection
collection of Candidate objects
Definition: VertexCompositeCandidateFwd.h:19
pat::PATLostTracks::TrkStatus::PFCANDNOTRKPROPS
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
VertexFwd.h
pat::PATLostTracks::TrkStatus::NOTUSED
reco::TrackBase::hitPattern
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:504
pat::PATLostTracks::vertices_
const edm::EDGetTokenT< reco::VertexCollection > vertices_
Definition: PATLostTracks.cc:50
eostools.move
def move(src, dest)
Definition: eostools.py:511
std
Definition: JetResolutionObject.h:76
edm::OrphanHandle
Definition: EDProductfwd.h:39
StringCutObjectSelector.h
extraflags_cff.vtx
vtx
Definition: extraflags_cff.py:19
RecoChargedCandidate.h
Vertex.h
PVValHelper::dz
Definition: PVValidationHelpers.h:51
pat::PATLostTracks::TrkStatus::PFELECTRON
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:29
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
electrons_cff.lostHits
lostHits
Definition: electrons_cff.py:391
pat::PATLostTracks::kshorts_
const edm::EDGetTokenT< reco::VertexCompositeCandidateCollection > kshorts_
Definition: PATLostTracks.cc:51
StringCutObjectSelector< reco::Track, false >
Exception
Definition: hltDiff.cc:245
l1t::VertexCollection
std::vector< Vertex > VertexCollection
Definition: Vertex.h:12
reco::HitPattern::MISSING_INNER_HITS
Definition: HitPattern.h:155
EgHLTOffHistBins_cfi.mass
mass
Definition: EgHLTOffHistBins_cfi.py:34
pat::PackedCandidate::NotReconstructedPrimary
Definition: PackedCandidate.h:723
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
reco::HitPattern::numberOfValidPixelHits
int numberOfValidPixelHits() const
Definition: HitPattern.h:818
reco::PFCandidate
Particle reconstructed by the particle flow algorithm.
Definition: PFCandidate.h:41
edm::Ref::key
key_type key() const
Accessor for product key.
Definition: Ref.h:250
pat::PATLostTracks::useLegacySetup_
const bool useLegacySetup_
Definition: PATLostTracks.cc:70
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
pat::PATLostTracks::associateTrkToVtx
std::pair< int, pat::PackedCandidate::PVAssociationQuality > associateTrkToVtx(const reco::VertexCollection &vertices, const reco::TrackRef &trk) const
Definition: PATLostTracks.cc:336
ParameterSet.h
EDProducer.h
pat::PATLostTracks::~PATLostTracks
~PATLostTracks() override
Definition: PATLostTracks.cc:128
pat::PATLostTracks::minPtToStoreLowQualityProps_
const double minPtToStoreLowQualityProps_
Definition: PATLostTracks.cc:59
pat::PATLostTracks::cands_
const edm::EDGetTokenT< reco::PFCandidateCollection > cands_
Definition: PATLostTracks.cc:47
edm::Event
Definition: Event.h:73
pat::PackedCandidate::UsedInFitTight
Definition: PackedCandidate.h:728
crabWrapper.key
key
Definition: crabWrapper.py:19
pat::PATLostTracks::passTrkCuts
bool passTrkCuts(const reco::Track &tr) const
Definition: PATLostTracks.cc:264
taus_updatedMVAIds_cff.mapping
mapping
Definition: taus_updatedMVAIds_cff.py:29
reco::TrackBase::quality
bool quality(const TrackQuality) const
Track quality.
Definition: TrackBase.h:552
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:14
pat::PATLostTracks::map_
const edm::EDGetTokenT< edm::Association< pat::PackedCandidateCollection > > map_
Definition: PATLostTracks.cc:48
PFCandidateFwd.h
pfMETCorrectionType0_cfi.minDz
minDz
Definition: pfMETCorrectionType0_cfi.py:52
reco::TrackBase::highPurity
Definition: TrackBase.h:154
Generated for CMSSW Reference Manual by doxygen 1.8.16