reco::JetSignalVertexCompatibilityAlgo Class Reference

#include <JetSignalVertexCompatibilityAlgo.h>

Classes
struct	RefToBaseLess

Public Member Functions
std::vector< float >	compatibility (const reco::VertexCollection &vertices, const reco::TrackRefVector &tracks) const
	JetSignalVertexCompatibilityAlgo (double cut, double temperature)
void	resetEvent (const TransientTrackBuilder *trackBuilder)
	~JetSignalVertexCompatibilityAlgo ()

Private Types
typedef std::map< reco::TrackBaseRef, reco::TransientTrack, RefToBaseLess< reco::Track > >	TransientTrackMap

Private Member Functions
double	activation (double compat) const
const TransientTrack &	convert (const reco::TrackBaseRef &track) const

Static Private Member Functions
static double	trackVertexCompat (const reco::Vertex &vtx, const TransientTrack &track)

Private Attributes
const double	cut
const double	temperature
const TransientTrackBuilder *	trackBuilder
TransientTrackMap	trackMap

Detailed Description

Definition at line 15 of file JetSignalVertexCompatibilityAlgo.h.

Member Typedef Documentation

TransientTrackMap

typedef std::map<reco::TrackBaseRef, reco::TransientTrack, RefToBaseLess<reco::Track> > reco::JetSignalVertexCompatibilityAlgo::TransientTrackMap

private

Definition at line 30 of file JetSignalVertexCompatibilityAlgo.h.

Constructor & Destructor Documentation

JetSignalVertexCompatibilityAlgo()

JetSignalVertexCompatibilityAlgo::JetSignalVertexCompatibilityAlgo	(	double	cut,
		double	temperature
	)

Definition at line 29 of file JetSignalVertexCompatibilityAlgo.cc.

    : cut(cut), temperature(temperature) {}

~JetSignalVertexCompatibilityAlgo()

JetSignalVertexCompatibilityAlgo::~JetSignalVertexCompatibilityAlgo

(

)

Definition at line 32 of file JetSignalVertexCompatibilityAlgo.cc.

{}

Member Function Documentation

activation()

double JetSignalVertexCompatibilityAlgo::activation ( double  compat ) const

private

Definition at line 59 of file JetSignalVertexCompatibilityAlgo.cc.

References cut, JetChargeProducer_cfi::exp, and temperature.

Referenced by compatibility().

                                                                       {
  return 1. / (std::exp((compat - cut) / temperature) + 1.);
}

compatibility()

std::vector< float > JetSignalVertexCompatibilityAlgo::compatibility	(	const reco::VertexCollection &	vertices,
		const reco::TrackRefVector &	tracks
	)		const

Definition at line 63 of file JetSignalVertexCompatibilityAlgo.cc.

References activation(), convert(), MillePedeFileConverter_cfg::e, mps_fire::i, mps_fire::result, HLT_2023v12_cff::track, pwdgSkimBPark_cfi::tracks, trackVertexCompat(), and AlignmentTracksFromVertexSelector_cfi::vertices.

Referenced by JetSignalVertexCompatibility::produce().

                                                                                                           {
  std::vector<float> result(vertices.size(), 0.);
  float sum = 0.;

  for (TrackRefVector::const_iterator track = tracks.begin(); track != tracks.end(); ++track) {
    const TransientTrack &transientTrack = convert(TrackBaseRef(*track));

    for (unsigned int i = 0; i < vertices.size(); i++) {
      double compat = trackVertexCompat(vertices[i], transientTrack);
      double contribution = activation(compat) * (*track)->pt();

      result[i] += contribution;
      sum += contribution;
    }
  }

  if (sum < 1.0e-9) {
    for (unsigned int i = 0; i < result.size(); i++)
      result[i] = 1.0 / result.size();
  } else {
    for (unsigned int i = 0; i < result.size(); i++)
      result[i] /= sum;
  }

  return result;
}

convert()

const TransientTrack & JetSignalVertexCompatibilityAlgo::convert ( const reco::TrackBaseRef &  track ) const

private

Definition at line 49 of file JetSignalVertexCompatibilityAlgo.cc.

References TransientTrackBuilder::build(), HLT_2023v12_cff::track, trackBuilder, and trackMap.

Referenced by compatibility().

                                                                                               {
  TransientTrackMap::iterator pos = trackMap.lower_bound(track);
  if (pos != trackMap.end() && pos->first == track)
    return pos->second;

  // the castTo will only work with regular, i.e. no GsfTracks
  // the interface is not intrinsically polymorph...
  return trackMap.insert(pos, std::make_pair(track, trackBuilder->build(track.castTo<TrackRef>())))->second;
}

resetEvent()

void JetSignalVertexCompatibilityAlgo::resetEvent

(

const TransientTrackBuilder *

trackBuilder

)

Definition at line 91 of file JetSignalVertexCompatibilityAlgo.cc.

References trackBuilder, and trackMap.

Referenced by JetSignalVertexCompatibility::produce().

                                                                                      {
  trackMap.clear();
  trackBuilder = builder;
}

trackVertexCompat()

double JetSignalVertexCompatibilityAlgo::trackVertexCompat ( const reco::Vertex &  vtx, const TransientTrack &  track  )

staticprivate

Definition at line 34 of file JetSignalVertexCompatibilityAlgo.cc.

References RecoVertex::convertError(), RecoVertex::convertPos(), DeadROC_duringRun::dir, TrajectoryStateOnSurface::globalPosition(), mathSSE::sqrt(), HLT_2023v12_cff::track, and L1BJetProducer_cff::vtx.

Referenced by compatibility().

                                                                                                         {
  GlobalPoint point1 = RecoVertex::convertPos(vtx.position());
  AnalyticalImpactPointExtrapolator extrap(track.field());
  TrajectoryStateOnSurface tsos = extrap.extrapolate(track.impactPointState(), point1);

  if (!tsos.isValid())
    return 1.0e6;

  GlobalPoint point2 = tsos.globalPosition();
  ROOT::Math::SVector<double, 3> dir(point1.x() - point2.x(), point1.y() - point2.y(), point1.z() - point2.z());
  GlobalError cov = RecoVertex::convertError(vtx.covariance()) + tsos.cartesianError().position();

  return ROOT::Math::Mag2(dir) / std::sqrt(ROOT::Math::Similarity(cov.matrix(), dir));
}

Member Data Documentation

cut

const double reco::JetSignalVertexCompatibilityAlgo::cut

private

Definition at line 40 of file JetSignalVertexCompatibilityAlgo.h.

Referenced by activation().

temperature

const double reco::JetSignalVertexCompatibilityAlgo::temperature

private

Definition at line 41 of file JetSignalVertexCompatibilityAlgo.h.

Referenced by activation().

trackBuilder

const TransientTrackBuilder* reco::JetSignalVertexCompatibilityAlgo::trackBuilder

private

Definition at line 38 of file JetSignalVertexCompatibilityAlgo.h.

Referenced by convert(), and resetEvent().

trackMap

TransientTrackMap reco::JetSignalVertexCompatibilityAlgo::trackMap

mutableprivate

Definition at line 37 of file JetSignalVertexCompatibilityAlgo.h.

Referenced by convert(), and resetEvent().