CMS 3D CMS Logo
CombinedSVComputer.h
Go to the documentation of this file.
1 #ifndef RecoBTag_SecondaryVertex_CombinedSVComputer_h
2 #define RecoBTag_SecondaryVertex_CombinedSVComputer_h
3 
4 #include <iostream>
5 #include <cstddef>
6 #include <string>
7 #include <cmath>
8 #include <vector>
9 
10 #include <Math/VectorUtil.h>
11 
12 #include "FWCore/ParameterSet/interface/ParameterSet.h"
13 #include "FWCore/Utilities/interface/Exception.h"
14 #include "FWCore/ParameterSet/interface/ParameterSet.h"
15 
16 #include "DataFormats/Math/interface/Vector3D.h"
17 #include "DataFormats/Math/interface/LorentzVector.h"
18 #include "DataFormats/GeometryCommonDetAlgo/interface/Measurement1D.h"
19 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
20 #include "DataFormats/GeometryVector/interface/GlobalVector.h"
21 #include "DataFormats/GeometryVector/interface/VectorUtil.h"
22 #include "DataFormats/TrackReco/interface/Track.h"
23 #include "DataFormats/TrackReco/interface/TrackFwd.h"
24 #include "DataFormats/BTauReco/interface/SecondaryVertexTagInfo.h"
25 #include "DataFormats/BTauReco/interface/CandSecondaryVertexTagInfo.h"
26 #include "DataFormats/BTauReco/interface/TaggingVariable.h"
27 #include "DataFormats/BTauReco/interface/VertexTypes.h"
28 #include "DataFormats/BTauReco/interface/ParticleMasses.h"
29 #include "DataFormats/JetReco/interface/PFJet.h"
30 #include "DataFormats/PatCandidates/interface/Jet.h"
31 
32 #include "RecoBTag/SecondaryVertex/interface/TrackSelector.h"
33 #include "RecoBTag/SecondaryVertex/interface/V0Filter.h"
34 #include "RecoBTag/SecondaryVertex/interface/TrackSorting.h"
35 #include "RecoBTag/SecondaryVertex/interface/TrackKinematics.h"
36 
37 #define range_for(i, x) for (int i = (x).begin; i != (x).end; i += (x).increment)
38 
39 class CombinedSVComputer {
40 public:
41  explicit CombinedSVComputer(const edm::ParameterSet &params);
42  virtual ~CombinedSVComputer() = default;
43  virtual reco::TaggingVariableList operator()(const reco::TrackIPTagInfo &ipInfo,
44  const reco::SecondaryVertexTagInfo &svInfo) const;
45  virtual reco::TaggingVariableList operator()(const reco::CandIPTagInfo &ipInfo,
46  const reco::CandSecondaryVertexTagInfo &svInfo) const;
47 
48  struct IterationRange {
49  int begin, end, increment;
50  };
51  double flipValue(double value, bool vertex) const;
52  IterationRange flipIterate(int size, bool vertex) const;
53 
54  edm::ParameterSet dropDeltaR(const edm::ParameterSet &pset) const;
55 
56  const reco::btag::TrackIPData &threshTrack(const reco::CandIPTagInfo &trackIPTagInfo,
57  const reco::btag::SortCriteria sort,
58  const reco::Jet &jet,
59  const GlobalPoint &pv) const;
60  const reco::btag::TrackIPData &threshTrack(const reco::TrackIPTagInfo &trackIPTagInfo,
61  const reco::btag::SortCriteria sort,
62  const reco::Jet &jet,
63  const GlobalPoint &pv) const;
64  template <class SVTI, class IPTI>
65  void fillCommonVariables(reco::TaggingVariableList &vars,
66  reco::TrackKinematics &vertexKinematics,
67  const IPTI &ipInfo,
68  const SVTI &svInfo,
69  double &vtx_track_ptSum,
70  double &vtx_track_ESum) const;
71 
72 private:
73  bool trackFlip;
74  bool vertexFlip;
75  double charmCut;
76  reco::btag::SortCriteria sortCriterium;
77  reco::TrackSelector trackSelector;
78  reco::TrackSelector trackNoDeltaRSelector;
79  reco::TrackSelector trackPseudoSelector;
80  unsigned int pseudoMultiplicityMin;
81  unsigned int trackMultiplicityMin;
82  double minTrackWeight;
83  bool useTrackWeights;
84  bool vertexMassCorrection;
85  reco::V0Filter pseudoVertexV0Filter;
86  reco::V0Filter trackPairV0Filter;
87  std::vector<reco::btau::TaggingVariableName> taggingVariables;
88 };
89 
90 template <class SVTI, class IPTI>
91 void CombinedSVComputer::fillCommonVariables(reco::TaggingVariableList &vars,
92  reco::TrackKinematics &vertexKinematics,
93  const IPTI &ipInfo,
94  const SVTI &svInfo,
95  double &vtx_track_ptSum,
96  double &vtx_track_ESum) const {
97  using namespace ROOT::Math;
98  using namespace reco;
99 
100  typedef typename IPTI::input_container Container;
101  typedef typename Container::value_type TrackRef;
102 
103  edm::RefToBase<Jet> jet = ipInfo.jet();
104  math::XYZVector jetDir = jet->momentum().Unit();
105  bool havePv = ipInfo.primaryVertex().isNonnull();
106  GlobalPoint pv;
107  if (havePv)
108  pv = GlobalPoint(ipInfo.primaryVertex()->x(), ipInfo.primaryVertex()->y(), ipInfo.primaryVertex()->z());
109 
110  btag::Vertices::VertexType vtxType = btag::Vertices::NoVertex;
111 
112  vars.insert(btau::jetPt, jet->pt(), true);
113  vars.insert(btau::jetEta, jet->eta(), true);
114  vars.insert(btau::jetAbsEta, fabs(jet->eta()), true);
115 
116  if (ipInfo.selectedTracks().size() < trackMultiplicityMin)
117  return;
118 
119  vars.insert(btau::jetNTracks, ipInfo.selectedTracks().size(), true);
120 
121  TrackKinematics allKinematics;
122  TrackKinematics trackJetKinematics;
123 
124  double jet_track_ESum = 0.;
125 
126  int vtx = -1;
127 
128  IterationRange range = flipIterate(svInfo.nVertices(), true);
129  range_for(i, range) {
130  if (vtx < 0)
131  vtx = i;
132  }
133 
134  if (vtx >= 0) {
135  vtxType = btag::Vertices::RecoVertex;
136  vars.insert(btau::flightDistance1dVal, flipValue(svInfo.flightDistance(vtx, 1).value(), true), true);
137  vars.insert(btau::flightDistance1dSig, flipValue(svInfo.flightDistance(vtx, 1).significance(), true), true);
138  vars.insert(btau::flightDistance2dVal, flipValue(svInfo.flightDistance(vtx, 2).value(), true), true);
139  vars.insert(btau::flightDistance2dSig, flipValue(svInfo.flightDistance(vtx, 2).significance(), true), true);
140  vars.insert(btau::flightDistance3dVal, flipValue(svInfo.flightDistance(vtx, 3).value(), true), true);
141  vars.insert(btau::flightDistance3dSig, flipValue(svInfo.flightDistance(vtx, 3).significance(), true), true);
142  vars.insert(btau::vertexJetDeltaR, Geom::deltaR(svInfo.flightDirection(vtx), vertexFlip ? -jetDir : jetDir), true);
143  vars.insert(btau::jetNSecondaryVertices, svInfo.nVertices(), true);
144  }
145 
146  std::vector<std::size_t> indices = ipInfo.sortedIndexes(sortCriterium);
147  const std::vector<reco::btag::TrackIPData> &ipData = ipInfo.impactParameterData();
148 
149  const Container &tracks = ipInfo.selectedTracks();
150  std::vector<TrackRef> pseudoVertexTracks;
151 
152  std::vector<TrackRef> trackPairV0Test(2);
153  range = flipIterate(indices.size(), false);
154  range_for(i, range) {
155  std::size_t idx = indices[i];
156  const reco::btag::TrackIPData &data = ipData[idx];
157  const TrackRef &track = tracks[idx];
158 
159  jet_track_ESum += std::sqrt(track->momentum().Mag2() + ROOT::Math::Square(ParticleMasses::piPlus));
160 
161  // add track to kinematics for all tracks in jet
162  //allKinematics.add(track); // would make more sense for some variables, e.g. vertexEnergyRatio nicely between 0 and 1, but not necessarily the best option for the discriminating power...
163 
164  // filter track -> this track selection can be tighter than the vertex track selection (used to fill the track related variables...)
165  if (!trackSelector(track, data, *jet, pv))
166  continue;
167 
168  // add track to kinematics for all tracks in jet
169  allKinematics.add(track);
170 
171  // if no vertex was reconstructed, attempt pseudo vertex
172  if (vtxType == btag::Vertices::NoVertex && trackPseudoSelector(track, data, *jet, pv)) {
173  pseudoVertexTracks.push_back(track);
174  vertexKinematics.add(track);
175  }
176 
177  // check against all other tracks for V0 track pairs
178  trackPairV0Test[0] = track;
179  bool ok = true;
180  range_for(j, range) {
181  if (i == j)
182  continue;
183 
184  std::size_t pairIdx = indices[j];
185  const reco::btag::TrackIPData &pairTrackData = ipData[pairIdx];
186  const TrackRef &pairTrack = tracks[pairIdx];
187 
188  if (!trackSelector(pairTrack, pairTrackData, *jet, pv))
189  continue;
190 
191  trackPairV0Test[1] = pairTrack;
192  if (!trackPairV0Filter(trackPairV0Test)) {
193  ok = false;
194  break;
195  }
196  }
197  if (!ok)
198  continue;
199 
200  trackJetKinematics.add(track);
201 
202  // add track variables
203  math::XYZVector trackMom = track->momentum();
204  double trackMag = std::sqrt(trackMom.Mag2());
205 
206  vars.insert(btau::trackSip3dVal, flipValue(data.ip3d.value(), false), true);
207  vars.insert(btau::trackSip3dSig, flipValue(data.ip3d.significance(), false), true);
208  vars.insert(btau::trackSip2dVal, flipValue(data.ip2d.value(), false), true);
209  vars.insert(btau::trackSip2dSig, flipValue(data.ip2d.significance(), false), true);
210  vars.insert(btau::trackJetDistVal, data.distanceToJetAxis.value(), true);
211  // vars.insert(btau::trackJetDistSig, data.distanceToJetAxis.significance(), true);
212  // vars.insert(btau::trackFirstTrackDist, data.distanceToFirstTrack, true);
213  // vars.insert(btau::trackGhostTrackVal, data.distanceToGhostTrack.value(), true);
214  // vars.insert(btau::trackGhostTrackSig, data.distanceToGhostTrack.significance(), true);
215  vars.insert(btau::trackDecayLenVal, havePv ? (data.closestToJetAxis - pv).mag() : -1.0, true);
216 
217  vars.insert(btau::trackMomentum, trackMag, true);
218  vars.insert(btau::trackEta, trackMom.Eta(), true);
219 
220  // check for erroneous Perp^2 values before evaluating Perp (fix for DeepCSV NaN outputs)
221  double perp_trackMom_jetDir = VectorUtil::Perp2(trackMom, jetDir);
222  perp_trackMom_jetDir = (perp_trackMom_jetDir > 0. ? std::sqrt(perp_trackMom_jetDir) : 0.);
223 
224  vars.insert(btau::trackPtRel, perp_trackMom_jetDir, true);
225  vars.insert(btau::trackPPar, jetDir.Dot(trackMom), true);
226  vars.insert(btau::trackDeltaR, VectorUtil::DeltaR(trackMom, jetDir), true);
227  vars.insert(btau::trackPtRatio, perp_trackMom_jetDir / trackMag, true);
228  vars.insert(btau::trackPParRatio, jetDir.Dot(trackMom) / trackMag, true);
229  }
230 
231  if (vtxType == btag::Vertices::NoVertex && vertexKinematics.numberOfTracks() >= pseudoMultiplicityMin &&
232  pseudoVertexV0Filter(pseudoVertexTracks)) {
233  vtxType = btag::Vertices::PseudoVertex;
234  for (typename std::vector<TrackRef>::const_iterator trkIt = pseudoVertexTracks.begin();
235  trkIt != pseudoVertexTracks.end();
236  ++trkIt) {
237  vars.insert(btau::trackEtaRel, reco::btau::etaRel(jetDir, (*trkIt)->momentum()), true);
238  vtx_track_ptSum += std::sqrt((*trkIt)->momentum().Perp2());
239  vtx_track_ESum += std::sqrt((*trkIt)->momentum().Mag2() + ROOT::Math::Square(ParticleMasses::piPlus));
240  }
241  }
242 
243  vars.insert(btau::vertexCategory, vtxType, true);
244 
245  vars.insert(btau::trackJetPt, trackJetKinematics.vectorSum().Pt(), true);
246  vars.insert(btau::trackSumJetDeltaR, VectorUtil::DeltaR(allKinematics.vectorSum(), jetDir), true);
247  vars.insert(btau::trackSumJetEtRatio, allKinematics.vectorSum().Et() / ipInfo.jet()->et(), true);
248 
249  vars.insert(btau::trackSip3dSigAboveCharm,
250  flipValue(threshTrack(ipInfo, reco::btag::IP3DSig, *jet, pv).ip3d.significance(), false),
251  true);
252  vars.insert(btau::trackSip3dValAboveCharm,
253  flipValue(threshTrack(ipInfo, reco::btag::IP3DSig, *jet, pv).ip3d.value(), false),
254  true);
255  vars.insert(btau::trackSip2dSigAboveCharm,
256  flipValue(threshTrack(ipInfo, reco::btag::IP2DSig, *jet, pv).ip2d.significance(), false),
257  true);
258  vars.insert(btau::trackSip2dValAboveCharm,
259  flipValue(threshTrack(ipInfo, reco::btag::IP2DSig, *jet, pv).ip2d.value(), false),
260  true);
261 
262  if (vtxType != btag::Vertices::NoVertex) {
263  math::XYZTLorentzVector allSum = useTrackWeights ? allKinematics.weightedVectorSum() : allKinematics.vectorSum();
264  math::XYZTLorentzVector vertexSum =
265  useTrackWeights ? vertexKinematics.weightedVectorSum() : vertexKinematics.vectorSum();
266 
267  if (vtxType != btag::Vertices::RecoVertex) {
268  vars.insert(btau::vertexNTracks, vertexKinematics.numberOfTracks(), true);
269  vars.insert(btau::vertexJetDeltaR, VectorUtil::DeltaR(vertexSum, jetDir), true);
270  }
271 
272  double vertexMass = vertexSum.M();
273  if (vtxType == btag::Vertices::RecoVertex && vertexMassCorrection) {
274  const GlobalVector &dir = svInfo.flightDirection(vtx);
275  double vertexPt2 = math::XYZVector(dir.x(), dir.y(), dir.z()).Cross(vertexSum).Mag2() / dir.mag2();
276  vertexMass = std::sqrt(vertexMass * vertexMass + vertexPt2) + std::sqrt(vertexPt2);
277  }
278  vars.insert(btau::vertexMass, vertexMass, true);
279 
280  double varPi = (vertexMass / 5.2794) *
281  (vtx_track_ESum / jet_track_ESum); // 5.2794 should be the average B meson mass of PDG! CHECK!!!
282  vars.insert(btau::massVertexEnergyFraction, varPi / (varPi + 0.04), true);
283  double varB = (std::sqrt(5.2794) * vtx_track_ptSum) / (vertexMass * std::sqrt(jet->pt()));
284  vars.insert(btau::vertexBoostOverSqrtJetPt, varB * varB / (varB * varB + 10.), true);
285 
286  if (allKinematics.numberOfTracks()) {
287  vars.insert(btau::vertexEnergyRatio, vertexSum.E() / allSum.E(), true);
288  } else {
289  vars.insert(btau::vertexEnergyRatio, 1, true);
290  }
291  }
292 
293  reco::PFJet const *pfJet = dynamic_cast<reco::PFJet const *>(&*jet);
294  pat::Jet const *patJet = dynamic_cast<pat::Jet const *>(&*jet);
295  if (pfJet != nullptr) {
296  vars.insert(btau::chargedHadronEnergyFraction, pfJet->chargedHadronEnergyFraction(), true);
297  vars.insert(btau::neutralHadronEnergyFraction, pfJet->neutralHadronEnergyFraction(), true);
298  vars.insert(btau::photonEnergyFraction, pfJet->photonEnergyFraction(), true);
299  vars.insert(btau::electronEnergyFraction, pfJet->electronEnergyFraction(), true);
300  vars.insert(btau::muonEnergyFraction, pfJet->muonEnergyFraction(), true);
301  vars.insert(btau::chargedHadronMultiplicity, pfJet->chargedHadronMultiplicity(), true);
302  vars.insert(btau::neutralHadronMultiplicity, pfJet->neutralHadronMultiplicity(), true);
303  vars.insert(btau::photonMultiplicity, pfJet->photonMultiplicity(), true);
304  vars.insert(btau::electronMultiplicity, pfJet->electronMultiplicity(), true);
305  vars.insert(btau::muonMultiplicity, pfJet->muonMultiplicity(), true);
306  vars.insert(
307  btau::hadronMultiplicity, pfJet->chargedHadronMultiplicity() + pfJet->neutralHadronMultiplicity(), true);
308  vars.insert(btau::hadronPhotonMultiplicity,
309  pfJet->chargedHadronMultiplicity() + pfJet->neutralHadronMultiplicity() + pfJet->photonMultiplicity(),
310  true);
311  vars.insert(btau::totalMultiplicity,
312  pfJet->chargedHadronMultiplicity() + pfJet->neutralHadronMultiplicity() + pfJet->photonMultiplicity() +
313  pfJet->electronMultiplicity() + pfJet->muonMultiplicity(),
314  true);
315 
316  } else if (patJet != nullptr && patJet->isPFJet()) {
317  vars.insert(btau::chargedHadronEnergyFraction, patJet->chargedHadronEnergyFraction(), true);
318  vars.insert(btau::neutralHadronEnergyFraction, patJet->neutralHadronEnergyFraction(), true);
319  vars.insert(btau::photonEnergyFraction, patJet->photonEnergyFraction(), true);
320  vars.insert(btau::electronEnergyFraction, patJet->electronEnergyFraction(), true);
321  vars.insert(btau::muonEnergyFraction, patJet->muonEnergyFraction(), true);
322  vars.insert(btau::chargedHadronMultiplicity, patJet->chargedHadronMultiplicity(), true);
323  vars.insert(btau::neutralHadronMultiplicity, patJet->neutralHadronMultiplicity(), true);
324  vars.insert(btau::photonMultiplicity, patJet->photonMultiplicity(), true);
325  vars.insert(btau::electronMultiplicity, patJet->electronMultiplicity(), true);
326  vars.insert(btau::muonMultiplicity, patJet->muonMultiplicity(), true);
327  vars.insert(
328  btau::hadronMultiplicity, patJet->chargedHadronMultiplicity() + patJet->neutralHadronMultiplicity(), true);
329  vars.insert(
330  btau::hadronPhotonMultiplicity,
331  patJet->chargedHadronMultiplicity() + patJet->neutralHadronMultiplicity() + patJet->photonMultiplicity(),
332  true);
333  vars.insert(btau::totalMultiplicity,
334  patJet->chargedHadronMultiplicity() + patJet->neutralHadronMultiplicity() +
335  patJet->photonMultiplicity() + patJet->electronMultiplicity() + patJet->muonMultiplicity(),
336  true);
337 
338  } else {
339  vars.insert(btau::chargedHadronEnergyFraction, 0., true);
340  vars.insert(btau::neutralHadronEnergyFraction, 0., true);
341  vars.insert(btau::photonEnergyFraction, 0., true);
342  vars.insert(btau::electronEnergyFraction, 0., true);
343  vars.insert(btau::muonEnergyFraction, 0., true);
344  vars.insert(btau::chargedHadronMultiplicity, 0, true);
345  vars.insert(btau::neutralHadronMultiplicity, 0, true);
346  vars.insert(btau::photonMultiplicity, 0, true);
347  vars.insert(btau::electronMultiplicity, 0, true);
348  vars.insert(btau::muonMultiplicity, 0, true);
349  vars.insert(btau::hadronMultiplicity, 0, true);
350  vars.insert(btau::hadronPhotonMultiplicity, 0, true);
351  vars.insert(btau::totalMultiplicity, 0, true);
352  }
353 }
354 
355 #endif // RecoBTag_SecondaryVertex_CombinedSVComputer_h
Vector3DBase
Definition: Vector3DBase.h:8
pat::Jet::electronMultiplicity
int electronMultiplicity() const
electronMultiplicity
Definition: Jet.h:451
reco::btau::jetPt
Definition: TaggingVariable.h:32
reco::btau::vertexBoostOverSqrtJetPt
Definition: TaggingVariable.h:118
bTagCombinedSVVariables_cff.indices
indices
Definition: bTagCombinedSVVariables_cff.py:67
reco::btau::trackPtRatio
Definition: TaggingVariable.h:51
CombinedSVComputer
Definition: CombinedSVComputer.h:39
FastTimerService_cff.range
range
Definition: FastTimerService_cff.py:34
reco::btau::chargedHadronEnergyFraction
Definition: TaggingVariable.h:103
reco::TrackKinematics
Definition: TrackKinematics.h:16
CombinedSVComputer::charmCut
double charmCut
Definition: CombinedSVComputer.h:75
PDWG_EXOHSCP_cff.tracks
tracks
Definition: PDWG_EXOHSCP_cff.py:28
reco::IPTagInfo
Definition: IPTagInfo.h:48
reco::ParticleMasses::piPlus
const double piPlus
Definition: ParticleMasses.h:9
mps_fire.i
i
Definition: mps_fire.py:428
reco::btau::neutralHadronEnergyFraction
Definition: TaggingVariable.h:104
reco::Jet
Base class for all types of Jets.
Definition: Jet.h:20
HLT_FULL_cff.track
track
Definition: HLT_FULL_cff.py:11713
reco::PFJet::neutralHadronMultiplicity
int neutralHadronMultiplicity() const
neutralHadronMultiplicity
Definition: PFJet.h:126
CombinedSVComputer::useTrackWeights
bool useTrackWeights
Definition: CombinedSVComputer.h:83
reco::btau::trackSumJetEtRatio
Definition: TaggingVariable.h:66
ParticleMasses.h
reco::btau::flightDistance1dSig
Definition: TaggingVariable.h:83
reco::btau::neutralHadronMultiplicity
Definition: TaggingVariable.h:109
CombinedSVComputer::trackNoDeltaRSelector
reco::TrackSelector trackNoDeltaRSelector
Definition: CombinedSVComputer.h:78
reco::btau::vertexCategory
Definition: TaggingVariable.h:69
reco::PFJet::photonEnergyFraction
float photonEnergyFraction() const
photonEnergyFraction
Definition: PFJet.h:105
CalibrationSummaryClient_cfi.params
params
Definition: CalibrationSummaryClient_cfi.py:14
VectorUtil.h
reco::btau::totalMultiplicity
Definition: TaggingVariable.h:115
reco::btau::muonEnergyFraction
Definition: TaggingVariable.h:107
reco::btau::trackSip2dSigAboveCharm
Definition: TaggingVariable.h:90
reco::TemplatedSecondaryVertexTagInfo
Definition: TemplatedSecondaryVertexTagInfo.h:47
PFJet.h
V0Filter.h
reco::btau::photonEnergyFraction
Definition: TaggingVariable.h:105
reco::TaggingVariableList
Definition: TaggingVariable.h:194
pat::Jet::photonEnergyFraction
float photonEnergyFraction() const
photonEnergyFraction (relative to corrected jet energy)
Definition: Jet.h:420
pat::Jet::neutralHadronEnergyFraction
float neutralHadronEnergyFraction() const
neutralHadronEnergyFraction (relative to uncorrected jet energy)
Definition: Jet.h:404
CombinedSVComputer::IterationRange::end
int end
Definition: CombinedSVComputer.h:49
Measurement1D.h
reco::PFJet::muonEnergyFraction
float muonEnergyFraction() const
muonEnergyFraction
Definition: PFJet.h:113
reco::btag::Vertices::PseudoVertex
Definition: VertexTypes.h:18
reco::btau::trackJetPt
Definition: TaggingVariable.h:33
CombinedSVComputer::IterationRange::begin
int begin
Definition: CombinedSVComputer.h:49
reco::btau::trackSip2dValAboveCharm
Definition: TaggingVariable.h:89
reco::btau::flightDistance3dVal
Definition: TaggingVariable.h:86
reco::btau::trackEta
Definition: TaggingVariable.h:42
reco::btau::trackSip3dValAboveCharm
Definition: TaggingVariable.h:91
reco::TrackKinematics::numberOfTracks
unsigned int numberOfTracks() const
Definition: TrackKinematics.h:44
reco::PFJet::chargedHadronEnergyFraction
float chargedHadronEnergyFraction() const
chargedHadronEnergyFraction
Definition: PFJet.h:97
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
reco::btau::trackEtaRel
Definition: TaggingVariable.h:49
reco::PFJet::electronEnergyFraction
float electronEnergyFraction() const
electronEnergyFraction
Definition: PFJet.h:109
reco::btag::SortCriteria
SortCriteria
Definition: IPTagInfo.h:43
convertSQLiteXML.ok
bool ok
Definition: convertSQLiteXML.py:98
reco::btau::massVertexEnergyFraction
Definition: TaggingVariable.h:117
CombinedSVComputer::vertexMassCorrection
bool vertexMassCorrection
Definition: CombinedSVComputer.h:84
reco::btau::trackSip2dVal
Definition: TaggingVariable.h:54
reco::JetExtendedAssociation::Container
edm::AssociationVector< reco::JetRefBaseProd, Values > Container
Definition: JetExtendedAssociation.h:29
edm::Ref< TrackCollection >
VertexTypes.h
CombinedSVComputer::flipValue
double flipValue(double value, bool vertex) const
Definition: CombinedSVComputer.cc:27
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
pat::Jet::neutralHadronMultiplicity
int neutralHadronMultiplicity() const
neutralHadronMultiplicity
Definition: Jet.h:447
CombinedSVComputer::IterationRange
Definition: CombinedSVComputer.h:48
reco::btag::Vertices::NoVertex
Definition: VertexTypes.h:18
pat::Jet
Analysis-level calorimeter jet class.
Definition: Jet.h:77
range_for
#define range_for(i, x)
Definition: CombinedSVComputer.h:37
Track.h
CombinedSVComputer::minTrackWeight
double minTrackWeight
Definition: CombinedSVComputer.h:82
TrackFwd.h
reco::btau::flightDistance1dVal
Definition: TaggingVariable.h:82
vars
vars
Definition: DeepTauId.cc:163
reco::btag::Vertices::RecoVertex
Definition: VertexTypes.h:18
SecondaryVertexTagInfo.h
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
CombinedSVComputer::trackMultiplicityMin
unsigned int trackMultiplicityMin
Definition: CombinedSVComputer.h:81
reco::btau::trackDeltaR
Definition: TaggingVariable.h:50
reco::btau::trackPtRel
Definition: TaggingVariable.h:47
reco::btau::trackMomentum
Definition: TaggingVariable.h:41
reco::btau::trackSumJetDeltaR
Definition: TaggingVariable.h:67
reco::PFJet::chargedHadronMultiplicity
int chargedHadronMultiplicity() const
chargedHadronMultiplicity
Definition: PFJet.h:124
HGCalGeometryMode::Square
Definition: HGCalGeometryMode.h:26
reco::btau::muonMultiplicity
Definition: TaggingVariable.h:112
CombinedSVComputer::trackSelector
reco::TrackSelector trackSelector
Definition: CombinedSVComputer.h:77
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
reco::TrackKinematics::vectorSum
const math::XYZTLorentzVector & vectorSum() const
Definition: TrackKinematics.h:47
Point3DBase< float, GlobalTag >
CombinedSVComputer::fillCommonVariables
void fillCommonVariables(reco::TaggingVariableList &vars, reco::TrackKinematics &vertexKinematics, const IPTI &ipInfo, const SVTI &svInfo, double &vtx_track_ptSum, double &vtx_track_ESum) const
Definition: CombinedSVComputer.h:91
ROOT::Math
Definition: Transform3DPJ.h:41
reco::V0Filter
Definition: V0Filter.h:14
PbPb_ZMuSkimMuonDPG_cff.deltaR
deltaR
Definition: PbPb_ZMuSkimMuonDPG_cff.py:63
reco::TrackRef
edm::Ref< TrackCollection > TrackRef
persistent reference to a Track
Definition: TrackFwd.h:20
reco::btau::hadronPhotonMultiplicity
Definition: TaggingVariable.h:114
CombinedSVComputer::flipIterate
IterationRange flipIterate(int size, bool vertex) const
Definition: CombinedSVComputer.cc:31
reco::btau::trackJetDistVal
Definition: TaggingVariable.h:60
reco::btau::jetEta
Definition: TaggingVariable.h:34
bphysicsOniaDQM_cfi.vertex
vertex
Definition: bphysicsOniaDQM_cfi.py:7
edm::ParameterSet
Definition: ParameterSet.h:47
reco::btag::Vertices::VertexType
VertexType
Definition: VertexTypes.h:18
reco::PFJet::neutralHadronEnergyFraction
float neutralHadronEnergyFraction() const
neutralHadronEnergyFraction
Definition: PFJet.h:101
reco::btau::jetAbsEta
Definition: TaggingVariable.h:35
reco::btau::trackSip2dSig
Definition: TaggingVariable.h:55
CombinedSVComputer::taggingVariables
std::vector< reco::btau::TaggingVariableName > taggingVariables
Definition: CombinedSVComputer.h:87
reco::btau::trackSip3dSig
Definition: TaggingVariable.h:57
pat::Jet::electronEnergyFraction
float electronEnergyFraction() const
electronEnergyFraction (relative to corrected jet energy)
Definition: Jet.h:426
reco::btau::trackPParRatio
Definition: TaggingVariable.h:52
reco::btau::jetNSecondaryVertices
Definition: TaggingVariable.h:72
reco::PFJet::electronMultiplicity
int electronMultiplicity() const
electronMultiplicity
Definition: PFJet.h:130
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:31
electrons_cff.ip3d
ip3d
Definition: electrons_cff.py:384
LorentzVector.h
MetAnalyzer.pv
def pv(vc)
Definition: MetAnalyzer.py:7
reco::PFJet::photonMultiplicity
int photonMultiplicity() const
photonMultiplicity
Definition: PFJet.h:128
reco::btau::etaRel
double etaRel(const math::XYZVector &dir, const math::XYZVector &track)
Definition: TaggingVariable.h:22
value
Definition: value.py:1
TrackKinematics.h
TrackSelector.h
reco::btau::electronEnergyFraction
Definition: TaggingVariable.h:106
reco::JetExtendedAssociation::value_type
Container::value_type value_type
Definition: JetExtendedAssociation.h:30
reco::btau::vertexJetDeltaR
Definition: TaggingVariable.h:80
electronAnalyzer_cfi.DeltaR
DeltaR
Definition: electronAnalyzer_cfi.py:33
CombinedSVComputer::IterationRange::increment
int increment
Definition: CombinedSVComputer.h:49
reco::btau::trackPPar
Definition: TaggingVariable.h:48
reco::btau::chargedHadronMultiplicity
Definition: TaggingVariable.h:108
Jet.h
CombinedSVComputer::~CombinedSVComputer
virtual ~CombinedSVComputer()=default
CombinedSVComputer::pseudoMultiplicityMin
unsigned int pseudoMultiplicityMin
Definition: CombinedSVComputer.h:80
CandSecondaryVertexTagInfo.h
reco::btau::trackSip3dVal
Definition: TaggingVariable.h:56
pat::Jet::isPFJet
bool isPFJet() const
check to see if the jet is a reco::PFJet
Definition: Jet.h:275
pat::Jet::chargedHadronMultiplicity
int chargedHadronMultiplicity() const
chargedHadronMultiplicity
Definition: Jet.h:445
CombinedSVComputer::trackPairV0Filter
reco::V0Filter trackPairV0Filter
Definition: CombinedSVComputer.h:86
reco::PFJet::muonMultiplicity
int muonMultiplicity() const
muonMultiplicity
Definition: PFJet.h:132
reco::TrackSelector
Definition: TrackSelector.h:14
pat::Jet::muonMultiplicity
int muonMultiplicity() const
muonMultiplicity
Definition: Jet.h:736
reco::TrackKinematics::weightedVectorSum
const math::XYZTLorentzVector & weightedVectorSum() const
Definition: TrackKinematics.h:48
pat::Jet::chargedHadronEnergyFraction
float chargedHadronEnergyFraction() const
chargedHadronEnergyFraction (relative to uncorrected jet energy)
Definition: Jet.h:400
CombinedSVComputer::vertexFlip
bool vertexFlip
Definition: CombinedSVComputer.h:74
CombinedSVComputer::pseudoVertexV0Filter
reco::V0Filter pseudoVertexV0Filter
Definition: CombinedSVComputer.h:85
extraflags_cff.vtx
vtx
Definition: extraflags_cff.py:18
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
metsig::jet
Definition: SignAlgoResolutions.h:47
pat::Jet::photonMultiplicity
int photonMultiplicity() const
photonMultiplicity
Definition: Jet.h:449
TrackSorting.h
reco::btau::hadronMultiplicity
Definition: TaggingVariable.h:113
pat::Jet::muonEnergyFraction
float muonEnergyFraction() const
muonEnergyFraction (relative to corrected jet energy)
Definition: Jet.h:432
CombinedSVComputer::CombinedSVComputer
CombinedSVComputer(const edm::ParameterSet &params)
Definition: CombinedSVComputer.cc:11
reco::PFJet
Jets made from PFObjects.
Definition: PFJet.h:20
GlobalVector.h
reco::btau::jetNTracks
Definition: TaggingVariable.h:37
CombinedSVComputer::sortCriterium
reco::btag::SortCriteria sortCriterium
Definition: CombinedSVComputer.h:76
edm::RefToBase< Jet >
Exception.h
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:79
reco::btau::trackDecayLenVal
Definition: TaggingVariable.h:58
reco::btau::flightDistance3dSig
Definition: TaggingVariable.h:87
reco::btau::flightDistance2dVal
Definition: TaggingVariable.h:84
reco::btau::electronMultiplicity
Definition: TaggingVariable.h:111
reco::btag::IP3DSig
Definition: IPTagInfo.h:43
reco::btau::photonMultiplicity
Definition: TaggingVariable.h:110
CombinedSVComputer::trackPseudoSelector
reco::TrackSelector trackPseudoSelector
Definition: CombinedSVComputer.h:79
ParameterSet.h
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
TaggingVariable.h
CombinedSVComputer::threshTrack
const reco::btag::TrackIPData & threshTrack(const reco::CandIPTagInfo &trackIPTagInfo, const reco::btag::SortCriteria sort, const reco::Jet &jet, const GlobalPoint &pv) const
Vector3D.h
reco::btau::flightDistance2dSig
Definition: TaggingVariable.h:85
GlobalPoint.h
reco::btau::vertexMass
Definition: TaggingVariable.h:75
reco::btau::vertexNTracks
Definition: TaggingVariable.h:76
CombinedSVComputer::operator()
virtual reco::TaggingVariableList operator()(const reco::TrackIPTagInfo &ipInfo, const reco::SecondaryVertexTagInfo &svInfo) const
Definition: CombinedSVComputer.cc:133
reco::btau::vertexEnergyRatio
Definition: TaggingVariable.h:79
reco::btau::trackSip3dSigAboveCharm
Definition: TaggingVariable.h:92
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
reco::btag::IP2DSig
Definition: IPTagInfo.h:43
DeadROC_duringRun.dir
dir
Definition: DeadROC_duringRun.py:23
reco::btag::TrackIPData
Definition: IPTagInfo.h:28
CombinedSVComputer::trackFlip
bool trackFlip
Definition: CombinedSVComputer.h:73
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
CombinedSVComputer::dropDeltaR
edm::ParameterSet dropDeltaR(const edm::ParameterSet &pset) const
Definition: CombinedSVComputer.cc:5
reco::TrackKinematics::add
void add(const reco::Track &track, double weight=1.0)
Definition: TrackKinematics.cc:57
Generated for CMSSW Reference Manual by doxygen 1.8.16