CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Classes | Public Member Functions | Private Member Functions | Private Attributes
CombinedSVComputer Class Reference

#include <CombinedSVComputer.h>

Classes

struct  IterationRange
 

Public Member Functions

 CombinedSVComputer (const edm::ParameterSet &params)
 
reco::TaggingVariableList operator() (const reco::TrackIPTagInfo &ipInfo, const reco::SecondaryVertexTagInfo &svInfo) const
 

Private Member Functions

IterationRange flipIterate (int size, bool vertex) const
 
double flipValue (double value, bool vertex) const
 
const
reco::TrackIPTagInfo::TrackIPData & 		threshTrack (const reco::TrackIPTagInfo &trackIPTagInfo, const reco::TrackIPTagInfo::SortCriteria sort, const reco::Jet &jet, const GlobalPoint &pv) const
 

Private Attributes

double charmCut
 
double minTrackWeight
 
unsigned int pseudoMultiplicityMin
 
reco::V0Filter pseudoVertexV0Filter
 
reco::TrackIPTagInfo::SortCriteria sortCriterium
 
bool trackFlip
 
unsigned int trackMultiplicityMin
 
reco::TrackSelector trackNoDeltaRSelector
 
reco::V0Filter trackPairV0Filter
 
reco::TrackSelector trackPseudoSelector
 
reco::TrackSelector trackSelector
 
bool useTrackWeights
 
bool vertexFlip
 
bool vertexMassCorrection
 

Detailed Description

Definition at line 14 of file CombinedSVComputer.h.

Constructor & Destructor Documentation

CombinedSVComputer::CombinedSVComputer ( const edm::ParameterSet &  params)

Definition at line 51 of file CombinedSVComputer.cc.

51  :
52  trackFlip(params.getParameter<bool>("trackFlip")),
53  vertexFlip(params.getParameter<bool>("vertexFlip")),
54  charmCut(params.getParameter<double>("charmCut")),
55  sortCriterium(TrackSorting::getCriterium(params.getParameter<std::string>("trackSort"))),
56  trackSelector(params.getParameter<edm::ParameterSet>("trackSelection")),
57  trackNoDeltaRSelector(dropDeltaR(params.getParameter<edm::ParameterSet>("trackSelection"))),
58  trackPseudoSelector(params.getParameter<edm::ParameterSet>("trackPseudoSelection")),
59  pseudoMultiplicityMin(params.getParameter<unsigned int>("pseudoMultiplicityMin")),
60  trackMultiplicityMin(params.getParameter<unsigned int>("trackMultiplicityMin")),
61  minTrackWeight(params.getParameter<double>("minimumTrackWeight")),
62  useTrackWeights(params.getParameter<bool>("useTrackWeights")),
63  vertexMassCorrection(params.getParameter<bool>("correctVertexMass")),
64  pseudoVertexV0Filter(params.getParameter<edm::ParameterSet>("pseudoVertexV0Filter")),
65  trackPairV0Filter(params.getParameter<edm::ParameterSet>("trackPairV0Filter"))
66 {
67 }
CombinedSVComputer::trackSelector
reco::TrackSelector trackSelector
Definition: CombinedSVComputer.h:38
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
dropDeltaR
static edm::ParameterSet dropDeltaR(const edm::ParameterSet &pset)
Definition: CombinedSVComputer.cc:44
TrackSorting::getCriterium
reco::TrackIPTagInfo::SortCriteria getCriterium(const std::string &name)
Definition: TrackSorting.cc:11
CombinedSVComputer::pseudoMultiplicityMin
unsigned int pseudoMultiplicityMin
Definition: CombinedSVComputer.h:41
CombinedSVComputer::trackNoDeltaRSelector
reco::TrackSelector trackNoDeltaRSelector
Definition: CombinedSVComputer.h:39
CombinedSVComputer::pseudoVertexV0Filter
reco::V0Filter pseudoVertexV0Filter
Definition: CombinedSVComputer.h:46
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
CombinedSVComputer::trackPseudoSelector
reco::TrackSelector trackPseudoSelector
Definition: CombinedSVComputer.h:40
CombinedSVComputer::trackPairV0Filter
reco::V0Filter trackPairV0Filter
Definition: CombinedSVComputer.h:47
CombinedSVComputer::sortCriterium
reco::TrackIPTagInfo::SortCriteria sortCriterium
Definition: CombinedSVComputer.h:37
CombinedSVComputer::trackMultiplicityMin
unsigned int trackMultiplicityMin
Definition: CombinedSVComputer.h:42

Member Function Documentation

CombinedSVComputer::IterationRange CombinedSVComputer::flipIterate ( int  size,
bool  vertex 
) const
inlineprivate

Definition at line 74 of file CombinedSVComputer.cc.

References CombinedSVComputer::IterationRange::begin, CombinedSVComputer::IterationRange::end, CombinedSVComputer::IterationRange::increment, findQualityFiles::size, trackFlip, and vertexFlip.

Referenced by operator()(), and threshTrack().

76 {
77  IterationRange range;
78  if (vertex ? vertexFlip : trackFlip) {
79  range.begin = size - 1;
80  range.end = -1;
81  range.increment = -1;
82  } else {
83  range.begin = 0;
84  range.end = size;
85  range.increment = +1;
86  }
87 
88  return range;
89 }
findQualityFiles.size
tuple size
Write out results.
Definition: findQualityFiles.py:442
double CombinedSVComputer::flipValue ( double  value,
bool  vertex 
) const
inlineprivate

Definition at line 69 of file CombinedSVComputer.cc.

References trackFlip, relativeConstraints::value, and vertexFlip.

Referenced by operator()().

70 {
71  return (vertex ? vertexFlip : trackFlip) ? -value : value;
72 }
TaggingVariableList CombinedSVComputer::operator() ( const reco::TrackIPTagInfo &  ipInfo,
const reco::SecondaryVertexTagInfo &  svInfo 
) const

Definition at line 140 of file CombinedSVComputer.cc.

References edm::RefVector< C, T, F >::begin(), reco::TrackIPTagInfo::TrackIPData::closestToJetAxis, data, Geom::deltaR(), dir, reco::TrackIPTagInfo::TrackIPData::distanceToJetAxis, edm::RefVector< C, T, F >::end(), reco::LeafCandidate::et(), etaRel(), reco::SecondaryVertexTagInfo::flightDirection(), reco::SecondaryVertexTagInfo::flightDistance(), reco::btau::flightDistance2dSig, reco::btau::flightDistance2dVal, reco::btau::flightDistance3dSig, reco::btau::flightDistance3dVal, flipIterate(), flipValue(), reco::Vertex::hasRefittedTracks(), i, customizeTrackingMonitorSeedNumber::idx, reco::TrackIPTagInfo::impactParameterData(), reco::TaggingVariableList::insert(), reco::TrackIPTagInfo::TrackIPData::ip2d, reco::TrackIPTagInfo::TrackIPData::ip3d, edm::Ref< C, T, F >::isNonnull(), j, reco::JTATagInfo::jet(), metsig::jet, reco::btau::jetEta, reco::btau::jetNSecondaryVertices, reco::btau::jetPt, PV3DBase< T, PVType, FrameType >::mag2(), minTrackWeight, reco::TrackBase::momentum(), reco::btag::Vertices::NoVertex, reco::SecondaryVertexTagInfo::nVertexTracks(), reco::SecondaryVertexTagInfo::nVertices(), convertSQLiteXML::ok, reco::TrackIPTagInfo::primaryVertex(), pseudoMultiplicityMin, reco::btag::Vertices::PseudoVertex, pseudoVertexV0Filter, range_for, reco::btag::Vertices::RecoVertex, reco::Vertex::refittedTrack(), reco::SecondaryVertexTagInfo::secondaryVertex(), reco::TrackIPTagInfo::selectedTracks(), Measurement1D::significance(), edm::RefVector< C, T, F >::size(), sortCriterium, reco::TrackIPTagInfo::sortedIndexes(), mathSSE::sqrt(), threshTrack(), reco::btau::trackDecayLenVal, reco::btau::trackDeltaR, reco::btau::trackEta, reco::btau::trackEtaRel, reco::btau::trackJetDistVal, reco::btau::trackMomentum, trackMultiplicityMin, trackPairV0Filter, reco::btau::trackPPar, reco::btau::trackPParRatio, trackPseudoSelector, reco::btau::trackPtRatio, reco::btau::trackPtRel, reco::JTATagInfo::tracks(), testEve_cfg::tracks, trackSelector, reco::btau::trackSip2dSig, reco::btau::trackSip2dSigAboveCharm, reco::btau::trackSip2dVal, reco::btau::trackSip3dSig, reco::btau::trackSip3dSigAboveCharm, reco::btau::trackSip3dVal, reco::btau::trackSumJetDeltaR, reco::btau::trackSumJetEtRatio, reco::SecondaryVertexTagInfo::trackWeight(), useTrackWeights, Measurement1D::value(), reco::btau::vertexCategory, reco::btau::vertexEnergyRatio, reco::btau::vertexJetDeltaR, reco::btau::vertexMass, vertexMassCorrection, reco::btau::vertexNTracks, reco::SecondaryVertexTagInfo::vertexTracks(), w(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

142 {
143  using namespace ROOT::Math;
144 
145  edm::RefToBase<Jet> jet = ipInfo.jet();
146  math::XYZVector jetDir = jet->momentum().Unit();
147  bool havePv = ipInfo.primaryVertex().isNonnull();
148  GlobalPoint pv;
149  if (havePv)
150  pv = GlobalPoint(ipInfo.primaryVertex()->x(),
151  ipInfo.primaryVertex()->y(),
152  ipInfo.primaryVertex()->z());
153 
154  btag::Vertices::VertexType vtxType = btag::Vertices::NoVertex;
155 
156  TaggingVariableList vars; // = ipInfo.taggingVariables();
157 
158  vars.insert(btau::jetPt, jet->pt(), true);
159  vars.insert(btau::jetEta, jet->eta(), true);
160 
161  if (ipInfo.tracks().size() < trackMultiplicityMin)
162  return vars;
163 
164  TrackKinematics allKinematics;
165  TrackKinematics vertexKinematics;
166 
167  int vtx = -1;
168  IterationRange range = flipIterate(svInfo.nVertices(), true);
169  range_for(i, range) {
170  if (vtx < 0)
171  vtx = i;
172 
173  const Vertex &vertex = svInfo.secondaryVertex(i);
174  bool hasRefittedTracks = vertex.hasRefittedTracks();
175  TrackRefVector tracks = svInfo.vertexTracks(i);
176  for(TrackRefVector::const_iterator track = tracks.begin();
177  track != tracks.end(); track++) {
178  double w = svInfo.trackWeight(i, *track);
179  if (w < minTrackWeight)
180  continue;
181  if (hasRefittedTracks) {
182  Track actualTrack =
183  vertex.refittedTrack(*track);
184  vertexKinematics.add(actualTrack, w);
185  vars.insert(btau::trackEtaRel, etaRel(jetDir,
186  actualTrack.momentum()), true);
187  } else {
188  vertexKinematics.add(**track, w);
189  vars.insert(btau::trackEtaRel, etaRel(jetDir,
190  (*track)->momentum()), true);
191  }
192  }
193  }
194 
195  if (vtx >= 0) {
196  vtxType = btag::Vertices::RecoVertex;
197 
198  vars.insert(btau::flightDistance2dVal,
199  flipValue(
200  svInfo.flightDistance(vtx, true).value(),
201  true),
202  true);
203  vars.insert(btau::flightDistance2dSig,
204  flipValue(
205  svInfo.flightDistance(vtx, true).significance(),
206  true),
207  true);
208  vars.insert(btau::flightDistance3dVal,
209  flipValue(
210  svInfo.flightDistance(vtx, false).value(),
211  true),
212  true);
213  vars.insert(btau::flightDistance3dSig,
214  flipValue(
215  svInfo.flightDistance(vtx, false).significance(),
216  true),
217  true);
218  vars.insert(btau::vertexJetDeltaR,
219  Geom::deltaR(svInfo.flightDirection(vtx), jetDir),true);
220  vars.insert(btau::jetNSecondaryVertices, svInfo.nVertices(), true);
221  vars.insert(btau::vertexNTracks, svInfo.nVertexTracks(), true);
222  }
223 
224  std::vector<std::size_t> indices = ipInfo.sortedIndexes(sortCriterium);
225  const std::vector<TrackIPTagInfo::TrackIPData> &ipData =
226  ipInfo.impactParameterData();
227  const edm::RefVector<TrackCollection> &tracks =
228  ipInfo.selectedTracks();
229  std::vector<TrackRef> pseudoVertexTracks;
230 
231  TrackRef trackPairV0Test[2];
232  range = flipIterate(indices.size(), false);
233  range_for(i, range) {
234  std::size_t idx = indices[i];
235  const TrackIPTagInfo::TrackIPData &data = ipData[idx];
236  const TrackRef &trackRef = tracks[idx];
237  const Track &track = *trackRef;
238 
239  // filter track
240 
241  if (!trackSelector(track, data, *jet, pv))
242  continue;
243 
244  // add track to kinematics for all tracks in jet
245 
246  allKinematics.add(track);
247 
248  // if no vertex was reconstructed, attempt pseudo vertex
249 
250  if (vtxType == btag::Vertices::NoVertex &&
251  trackPseudoSelector(track, data, *jet, pv)) {
252  pseudoVertexTracks.push_back(trackRef);
253  vertexKinematics.add(track);
254  }
255 
256  // check against all other tracks for V0 track pairs
257 
258  trackPairV0Test[0] = tracks[idx];
259  bool ok = true;
260  range_for(j, range) {
261  if (i == j)
262  continue;
263 
264  std::size_t pairIdx = indices[j];
265  const TrackIPTagInfo::TrackIPData &pairTrackData =
266  ipData[pairIdx];
267  const TrackRef &pairTrackRef = tracks[pairIdx];
268  const Track &pairTrack = *pairTrackRef;
269 
270  if (!trackSelector(pairTrack, pairTrackData, *jet, pv))
271  continue;
272 
273  trackPairV0Test[1] = pairTrackRef;
274  if (!trackPairV0Filter(trackPairV0Test, 2)) {
275  ok = false;
276  break;
277  }
278  }
279  if (!ok)
280  continue;
281 
282  // add track variables
283 
284  math::XYZVector trackMom = track.momentum();
285  double trackMag = std::sqrt(trackMom.Mag2());
286 
287  vars.insert(btau::trackSip3dVal,
288  flipValue(data.ip3d.value(), false), true);
289  vars.insert(btau::trackSip3dSig,
290  flipValue(data.ip3d.significance(), false), true);
291  vars.insert(btau::trackSip2dVal,
292  flipValue(data.ip2d.value(), false), true);
293  vars.insert(btau::trackSip2dSig,
294  flipValue(data.ip2d.significance(), false), true);
295  vars.insert(btau::trackJetDistVal, data.distanceToJetAxis.value(), true);
296 // vars.insert(btau::trackJetDistSig, data.distanceToJetAxis.significance(), true);
297 // vars.insert(btau::trackFirstTrackDist, data.distanceToFirstTrack, true);
298 // vars.insert(btau::trackGhostTrackVal, data.distanceToGhostTrack.value(), true);
299 // vars.insert(btau::trackGhostTrackSig, data.distanceToGhostTrack.significance(), true);
300  vars.insert(btau::trackDecayLenVal, havePv ? (data.closestToJetAxis - pv).mag() : -1.0, true);
301 
302  vars.insert(btau::trackMomentum, trackMag, true);
303  vars.insert(btau::trackEta, trackMom.Eta(), true);
304 
305  vars.insert(btau::trackPtRel, VectorUtil::Perp(trackMom, jetDir), true);
306  vars.insert(btau::trackPPar, jetDir.Dot(trackMom), true);
307  vars.insert(btau::trackDeltaR, VectorUtil::DeltaR(trackMom, jetDir), true);
308  vars.insert(btau::trackPtRatio, VectorUtil::Perp(trackMom, jetDir) / trackMag, true);
309  vars.insert(btau::trackPParRatio, jetDir.Dot(trackMom) / trackMag, true);
310  }
311 
312  if (vtxType == btag::Vertices::NoVertex &&
313  vertexKinematics.numberOfTracks() >= pseudoMultiplicityMin &&
314  pseudoVertexV0Filter(pseudoVertexTracks)) {
315  vtxType = btag::Vertices::PseudoVertex;
316  for(std::vector<TrackRef>::const_iterator track =
317  pseudoVertexTracks.begin();
318  track != pseudoVertexTracks.end(); ++track)
319  vars.insert(btau::trackEtaRel, etaRel(jetDir,
320  (*track)->momentum()), true);
321  }
322 
323  vars.insert(btau::vertexCategory, vtxType, true);
324  vars.insert(btau::trackSumJetDeltaR,
325  VectorUtil::DeltaR(allKinematics.vectorSum(), jetDir), true);
326  vars.insert(btau::trackSumJetEtRatio,
327  allKinematics.vectorSum().Et() / ipInfo.jet()->et(), true);
328  vars.insert(btau::trackSip3dSigAboveCharm,
329  flipValue(
330  threshTrack(ipInfo, TrackIPTagInfo::IP3DSig, *jet, pv)
331  .ip3d.significance(),
332  false),
333  true);
334  vars.insert(btau::trackSip2dSigAboveCharm,
335  flipValue(
336  threshTrack(ipInfo, TrackIPTagInfo::IP2DSig, *jet, pv)
337  .ip2d.significance(),
338  false),
339  true);
340 
341  if (vtxType != btag::Vertices::NoVertex) {
342  math::XYZTLorentzVector allSum = useTrackWeights
343  ? allKinematics.weightedVectorSum()
344  : allKinematics.vectorSum();
345  math::XYZTLorentzVector vertexSum = useTrackWeights
346  ? vertexKinematics.weightedVectorSum()
347  : vertexKinematics.vectorSum();
348 
349  if (vtxType != btag::Vertices::RecoVertex) {
350  vars.insert(btau::vertexNTracks,
351  vertexKinematics.numberOfTracks(), true);
352  vars.insert(btau::vertexJetDeltaR,
353  VectorUtil::DeltaR(vertexSum, jetDir), true);
354  }
355 
356  double vertexMass = vertexSum.M();
357  if (vtxType == btag::Vertices::RecoVertex &&
358  vertexMassCorrection) {
359  GlobalVector dir = svInfo.flightDirection(vtx);
360  double vertexPt2 =
361  math::XYZVector(dir.x(), dir.y(), dir.z()).
362  Cross(vertexSum).Mag2() / dir.mag2();
363  vertexMass = std::sqrt(vertexMass * vertexMass +
364  vertexPt2) + std::sqrt(vertexPt2);
365  }
366  vars.insert(btau::vertexMass, vertexMass, true);
367  if (allKinematics.numberOfTracks())
368  vars.insert(btau::vertexEnergyRatio,
369  vertexSum.E() / allSum.E(), true);
370  else
371  vars.insert(btau::vertexEnergyRatio, 1, true);
372  }
373 
374  vars.finalize();
375 
376  return vars;
377 }
CombinedSVComputer::trackSelector
reco::TrackSelector trackSelector
Definition: CombinedSVComputer.h:38
reco::LeafCandidate::et
virtual double et() const GCC11_FINAL
transverse energy
Definition: LeafCandidate.h:136
i
int i
Definition: DBlmapReader.cc:9
PV3DBase::mag2
T mag2() const
Definition: PV3DBase.h:66
reco::TrackBase::momentum
const Vector & momentum() const
track momentum vector
Definition: TrackBase.h:150
reco::SecondaryVertexTagInfo::flightDistance
Measurement1D flightDistance(unsigned int index, bool in2d=false) const
Definition: SecondaryVertexTagInfo.h:98
CombinedSVComputer::pseudoMultiplicityMin
unsigned int pseudoMultiplicityMin
Definition: CombinedSVComputer.h:41
reco::TrackIPTagInfo::primaryVertex
const edm::Ref< VertexCollection > & primaryVertex() const
Definition: TrackIPTagInfo.h:107
reco::Vertex::refittedTrack
Track refittedTrack(const TrackBaseRef &track) const
reco::Vertex::hasRefittedTracks
bool hasRefittedTracks() const
Checks whether refitted tracks are stored.
Definition: Vertex.h:121
CombinedSVComputer::pseudoVertexV0Filter
reco::V0Filter pseudoVertexV0Filter
Definition: CombinedSVComputer.h:46
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
range_for
#define range_for(i, x)
Definition: CombinedSVComputer.cc:41
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
edm::RefVector::end
const_iterator end() const
Termination of iteration.
Definition: RefVector.h:249
edm::RefVector::begin
const_iterator begin() const
Initialize an iterator over the RefVector.
Definition: RefVector.h:244
reco::JTATagInfo::tracks
virtual TrackRefVector tracks(void) const
returns a list of tracks associated to the jet
Definition: JTATagInfo.h:21
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:30
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:250
reco::SecondaryVertexTagInfo::flightDirection
const GlobalVector & flightDirection(unsigned int index) const
Definition: SecondaryVertexTagInfo.h:100
reco::SecondaryVertexTagInfo::secondaryVertex
const Vertex & secondaryVertex(unsigned int index) const
Definition: SecondaryVertexTagInfo.h:72
CombinedSVComputer::trackPseudoSelector
reco::TrackSelector trackPseudoSelector
Definition: CombinedSVComputer.h:40
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
j
int j
Definition: DBlmapReader.cc:9
CombinedSVComputer::threshTrack
const reco::TrackIPTagInfo::TrackIPData & threshTrack(const reco::TrackIPTagInfo &trackIPTagInfo, const reco::TrackIPTagInfo::SortCriteria sort, const reco::Jet &jet, const GlobalPoint &pv) const
Definition: CombinedSVComputer.cc:92
reco::SecondaryVertexTagInfo::trackWeight
float trackWeight(unsigned int svIndex, unsigned int trackindex) const
Definition: SecondaryVertexTagInfo.cc:161
CombinedSVComputer::flipIterate
IterationRange flipIterate(int size, bool vertex) const
Definition: CombinedSVComputer.cc:74
reco::JTATagInfo::jet
virtual edm::RefToBase< Jet > jet(void) const
returns a polymorphic reference to the tagged jet
Definition: JTATagInfo.h:20
reco::SecondaryVertexTagInfo::vertexTracks
TrackRefVector vertexTracks() const
Definition: SecondaryVertexTagInfo.cc:80
CombinedSVComputer::trackPairV0Filter
reco::V0Filter trackPairV0Filter
Definition: CombinedSVComputer.h:47
reco::TrackIPTagInfo::sortedIndexes
std::vector< size_t > sortedIndexes(SortCriteria mode=IP3DSig) const
Definition: TrackIPTagInfo.cc:78
etaRel
static double etaRel(const math::XYZVector &dir, const math::XYZVector &track)
Definition: CombinedSVComputer.cc:130
Measurement1D::significance
double significance() const
Definition: Measurement1D.h:32
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:31
customizeTrackingMonitorSeedNumber.idx
tuple idx
DEBUGGING if hasattr(process,&quot;trackMonIterativeTracking2012&quot;): print &quot;trackMonIterativeTracking2012 D...
Definition: customizeTrackingMonitorSeedNumber.py:15
CombinedSVComputer::flipValue
double flipValue(double value, bool vertex) const
Definition: CombinedSVComputer.cc:69
Measurement1D::value
double value() const
Definition: Measurement1D.h:28
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
Geom::deltaR
double deltaR(const Vector1 &v1, const Vector2 &v2)
Definition: VectorUtil.h:84
reco::TrackIPTagInfo::impactParameterData
const std::vector< TrackIPData > & impactParameterData() const
Definition: TrackIPTagInfo.h:68
edm::RefVector::size
size_type size() const
Size of the RefVector.
Definition: RefVector.h:89
w
T w() const
Definition: extBasic3DVector.h:245
dir
dbl *** dir
Definition: mlp_gen.cc:35
reco::TrackIPTagInfo::selectedTracks
const edm::RefVector< TrackCollection > & selectedTracks() const
Definition: TrackIPTagInfo.h:75
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
CombinedSVComputer::sortCriterium
reco::TrackIPTagInfo::SortCriteria sortCriterium
Definition: CombinedSVComputer.h:37
reco::TaggingVariableList::insert
void insert(const TaggingVariable &variable, bool delayed=false)
Definition: TaggingVariable.cc:178
CombinedSVComputer::trackMultiplicityMin
unsigned int trackMultiplicityMin
Definition: CombinedSVComputer.h:42
const TrackIPTagInfo::TrackIPData & CombinedSVComputer::threshTrack ( const reco::TrackIPTagInfo &  trackIPTagInfo,
const reco::TrackIPTagInfo::SortCriteria  sort,
const reco::Jet &  jet,
const GlobalPoint &  pv 
) const
private

Definition at line 92 of file CombinedSVComputer.cc.

References charmCut, data, flipIterate(), i, customizeTrackingMonitorSeedNumber::idx, reco::TrackIPTagInfo::impactParameterData(), range_for, reco::TrackIPTagInfo::selectedTracks(), reco::TrackIPTagInfo::sortedIndexes(), trackNoDeltaRSelector, and testEve_cfg::tracks.

Referenced by operator()().

96 {
97  const edm::RefVector<TrackCollection> &tracks =
98  trackIPTagInfo.selectedTracks();
99  const std::vector<TrackIPTagInfo::TrackIPData> &ipData =
100  trackIPTagInfo.impactParameterData();
101  std::vector<std::size_t> indices = trackIPTagInfo.sortedIndexes(sort);
102 
103  IterationRange range = flipIterate(indices.size(), false);
104  TrackKinematics kin;
105  range_for(i, range) {
106  std::size_t idx = indices[i];
107  const TrackIPTagInfo::TrackIPData &data = ipData[idx];
108  const Track &track = *tracks[idx];
109 
110  if (!trackNoDeltaRSelector(track, data, jet, pv))
111  continue;
112 
113  kin.add(track);
114  if (kin.vectorSum().M() > charmCut)
115  return data;
116  }
117 
118  static const TrackIPTagInfo::TrackIPData dummy = {
119  GlobalPoint(),
120  GlobalPoint(),
121  Measurement1D(-1.0, 1.0),
122  Measurement1D(-1.0, 1.0),
123  Measurement1D(-1.0, 1.0),
124  Measurement1D(-1.0, 1.0),
125  0.
126  };
127  return dummy;
128 }
i
int i
Definition: DBlmapReader.cc:9
CombinedSVComputer::trackNoDeltaRSelector
reco::TrackSelector trackNoDeltaRSelector
Definition: CombinedSVComputer.h:39
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
range_for
#define range_for(i, x)
Definition: CombinedSVComputer.cc:41
CombinedSVComputer::flipIterate
IterationRange flipIterate(int size, bool vertex) const
Definition: CombinedSVComputer.cc:74
reco::TrackIPTagInfo::sortedIndexes
std::vector< size_t > sortedIndexes(SortCriteria mode=IP3DSig) const
Definition: TrackIPTagInfo.cc:78
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
customizeTrackingMonitorSeedNumber.idx
tuple idx
DEBUGGING if hasattr(process,&quot;trackMonIterativeTracking2012&quot;): print &quot;trackMonIterativeTracking2012 D...
Definition: customizeTrackingMonitorSeedNumber.py:15
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
reco::TrackIPTagInfo::impactParameterData
const std::vector< TrackIPData > & impactParameterData() const
Definition: TrackIPTagInfo.h:68
reco::TrackIPTagInfo::selectedTracks
const edm::RefVector< TrackCollection > & selectedTracks() const
Definition: TrackIPTagInfo.h:75

Member Data Documentation

double CombinedSVComputer::charmCut
private

Definition at line 36 of file CombinedSVComputer.h.

Referenced by threshTrack().

double CombinedSVComputer::minTrackWeight
private

Definition at line 43 of file CombinedSVComputer.h.

Referenced by operator()().

unsigned int CombinedSVComputer::pseudoMultiplicityMin
private

Definition at line 41 of file CombinedSVComputer.h.

Referenced by operator()().

reco::V0Filter CombinedSVComputer::pseudoVertexV0Filter
private

Definition at line 46 of file CombinedSVComputer.h.

Referenced by operator()().

reco::TrackIPTagInfo::SortCriteria CombinedSVComputer::sortCriterium
private

Definition at line 37 of file CombinedSVComputer.h.

Referenced by operator()().

bool CombinedSVComputer::trackFlip
private

Definition at line 34 of file CombinedSVComputer.h.

Referenced by flipIterate(), and flipValue().

unsigned int CombinedSVComputer::trackMultiplicityMin
private

Definition at line 42 of file CombinedSVComputer.h.

Referenced by operator()().

reco::TrackSelector CombinedSVComputer::trackNoDeltaRSelector
private

Definition at line 39 of file CombinedSVComputer.h.

Referenced by threshTrack().

reco::V0Filter CombinedSVComputer::trackPairV0Filter
private

Definition at line 47 of file CombinedSVComputer.h.

Referenced by operator()().

reco::TrackSelector CombinedSVComputer::trackPseudoSelector
private

Definition at line 40 of file CombinedSVComputer.h.

Referenced by operator()().

reco::TrackSelector CombinedSVComputer::trackSelector
private

Definition at line 38 of file CombinedSVComputer.h.

Referenced by operator()().

bool CombinedSVComputer::useTrackWeights
private

Definition at line 44 of file CombinedSVComputer.h.

Referenced by operator()().

bool CombinedSVComputer::vertexFlip
private

Definition at line 35 of file CombinedSVComputer.h.

Referenced by flipIterate(), and flipValue().

bool CombinedSVComputer::vertexMassCorrection
private

Definition at line 45 of file CombinedSVComputer.h.

Referenced by operator()().

Generated for CMSSW Reference Manual by doxygen 1.8.5