#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.

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

◆ ~JetSignalVertexCompatibilityAlgo()

JetSignalVertexCompatibilityAlgo::~JetSignalVertexCompatibilityAlgo ( )

Definition at line 32 of file JetSignalVertexCompatibilityAlgo.cc.

32 {}

Member Function Documentation

◆ activation()

double JetSignalVertexCompatibilityAlgo::activation ( double compat ) const

private

Definition at line 59 of file JetSignalVertexCompatibilityAlgo.cc.

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

References cut, JetChargeProducer_cfi::exp, and temperature.

Referenced by compatibility().

◆ compatibility()

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

Definition at line 63 of file JetSignalVertexCompatibilityAlgo.cc.

                                                                                                            {
   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;
 }

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

Referenced by JetSignalVertexCompatibility::produce().

◆ convert()

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

private

Definition at line 49 of file JetSignalVertexCompatibilityAlgo.cc.

                                                                                                {
   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;
 }

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

Referenced by compatibility().

◆ resetEvent()

void JetSignalVertexCompatibilityAlgo::resetEvent ( const TransientTrackBuilder * trackBuilder )

Definition at line 91 of file JetSignalVertexCompatibilityAlgo.cc.

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

References trackBuilder, and trackMap.

Referenced by JetSignalVertexCompatibility::produce().

◆ trackVertexCompat()

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

staticprivate

Definition at line 34 of file JetSignalVertexCompatibilityAlgo.cc.

                                                                                                          {
   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));
 }