pfalgo_common_ref.h File Reference

#include "../firmware/data.h"
#include "../firmware/pfalgo_common.h"
#include <algorithm>

Go to the source code of this file.

Classes
struct	pfalgo_config

Functions
template<typename CO_t >
int	best_match_with_pt_ref (int nCAL, int dR2MAX, const CO_t calo[], const TkObj &track)
void	pfalgo_mu_ref (const pfalgo_config &cfg, const TkObj track[], const MuObj mu[], bool isMu[], PFChargedObj outmu[], bool debug)
template<typename T >
void	ptsort_ref (int nIn, int nOut, const T in[], T out[])
template<typename T , typename TV >
void	ptsort_ref (int nIn, int nOut, const TV &in, T out[])
template<typename T >
int	sqr (const T &t)

Function Documentation

template<typename CO_t >

int best_match_with_pt_ref	(	int	nCAL,
		int	dR2MAX,
		const CO_t	calo[],
		const TkObj &	track
	)

Definition at line 56 of file pfalgo_common_ref.h.

References runTauDisplay::dr, dr2_int(), TkObj::hwEta, TkObj::hwPhi, TkObj::hwPt, TkObj::hwPtErr, and sqr().

                                                                                         {
  pt_t caloPtMin = track.hwPt - 2 * (track.hwPtErr);
  if (caloPtMin < 0)
    caloPtMin = 0;
  int dptscale = (dR2MAX << 8) / std::max<int>(1, sqr(track.hwPtErr));
  int drmin = 0, ibest = -1;
  for (int ic = 0; ic < nCAL; ++ic) {
    if (calo[ic].hwPt <= caloPtMin)
      continue;
    int dr = dr2_int(track.hwEta, track.hwPhi, calo[ic].hwEta, calo[ic].hwPhi);
    if (dr >= dR2MAX)
      continue;
    dr += ((sqr(std::max<int>(track.hwPt - calo[ic].hwPt, 0)) * dptscale) >> 8);
    if (ibest == -1 || dr < drmin) {
      drmin = dr;
      ibest = ic;
    }
  }
  return ibest;
}

void pfalgo_mu_ref	(	const pfalgo_config &	cfg,
		const TkObj	track[],
		const MuObj	mu[],
		bool	isMu[],
		PFChargedObj	outmu[],
		bool	debug
	)

template<typename T >

void ptsort_ref	(	int	nIn,
		int	nOut,
		const T	in[],
		T	out[]
	)

template<typename T , typename TV >

void ptsort_ref	(	int	nIn,
		int	nOut,
		const TV &	in,
		T	out[]
	)

Definition at line 78 of file pfalgo_common_ref.h.

                                                             {
  for (int iout = 0; iout < nOut; ++iout) {
    out[iout].hwPt = 0;
  }
  for (int it = 0; it < nIn; ++it) {
    for (int iout = 0; iout < nOut; ++iout) {
      if (in[it].hwPt >= out[iout].hwPt) {
        for (int i2 = nOut - 1; i2 > iout; --i2) {
          out[i2] = out[i2 - 1];
        }
        out[iout] = in[it];
        break;
      }
    }
  }
}

template<typename T >

int sqr ( const T &  t )

inline

Definition at line 9 of file pfalgo_common_ref.h.

References submitPVValidationJobs::t.

Referenced by ThirdHitPredictionFromInvLine::add(), MonopoleTransportation::AlongStepGetPhysicalInteractionLength(), ThirdHitPredictionFromCircle::angle(), best_match_with_pt_ref(), PixelTrackBuilder::build(), ThirdHitPrediction::calculateRangesBarrel(), ThirdHitPredictionFromInvLine::check(), GlobalTrackingRegion::checkRZ(), HitPairGeneratorFromLayerPairForPhotonConversion::checkRZCompatibilityWithSeedTrack(), RectangularEtaPhiTrackingRegion::checkRZOld(), ParabolaFit::chi2(), MuonSeedFromRecHits::createSeed(), ThirdHitPredictionFromInvLine::crossing(), ThirdHitPredictionFromCircle::curvature(), MSLayer::distance2(), RectangularEtaPhiTrackingRegion::estimator(), HitQuadrupletGeneratorFromLayerPairForPhotonConversion::failCheckRZCompatibility(), ThirdHitPrediction::findMinimalCircles(), ThirdHitPrediction::findRectangle(), ThirdHitPrediction::findTouchingCircles(), ThirdHitPrediction::fitParabola(), FullModelReactionDynamics::GenerateXandPt(), L1MuonPixelTrackFitter::getBendingError(), HitQuadrupletGeneratorFromLayerPairForPhotonConversion::getEffectiveErrorOnZ(), ClusterShapeTrackFilter::getGlobalDirs(), HitPairGeneratorFromLayerPairForPhotonConversion::getPhiRange(), HitQuadrupletGeneratorFromLayerPairForPhotonConversion::getSqrEffectiveErrorOnZ(), SeedForPhotonConversionFromQuadruplets::getSqrEffectiveErrorOnZ(), CAHitTripletGenerator::hitNtuplets(), CAHitQuadrupletGenerator::hitNtuplets(), HitPairGeneratorFromLayerPairForPhotonConversion::hitPairs(), PixelTripletNoTipGenerator::hitTriplets(), ThirdHitCorrection::init(), SeedFromConsecutiveHitsCreator::initialError(), SeedForPhotonConversion1Leg::initialError(), SeedForPhotonConversionFromQuadruplets::initialError(), ThirdHitPrediction::invertCircle(), ThirdHitPrediction::isCompatibleWithMultipleScattering(), JacobianLocalToCartesian::JacobianLocalToCartesian(), SeedFromConsecutiveHitsCreator::makeSeed(), MatchJet::matchCollections(), RKCurvilinearDistance< T, N >::momentum(), KFBasedPixelFitter::MyBeamSpotHit::MyBeamSpotHit(), PixelTrackFilterByKinematics::operator()(), RKCurvilinearDistance< T, N >::operator()(), ThirdHitZPrediction::operator()(), ThirdHitPredictionFromCircle::phi(), InnerDeltaPhi::phiRange(), SeedGeneratorFromProtoTracksEDProducer::produce(), ThirdHitPredictionFromCircle::HelixRZ::rAtZ(), reco::GhostTrackVertexFinder::reassignTracks(), ThirdHitPredictionFromInvLine::remove(), PrimaryVertexResolution::Resolution::Resolution(), ParabolaFit::result(), PixelFitterByConformalMappingAndLine::run(), KFBasedPixelFitter::run(), PhotonConversionTrajectorySeedProducerFromSingleLegAlgo::selectPriVtxCompatibleWithTrack(), MSLayersKeeperX0AtEta::setX0(), SeedForPhotonConversionFromQuadruplets::simpleGetSlope(), reco::SequentialGhostTrackFitter::stable(), MSLayersAtAngle::sumX0D(), ThirdHitPredictionFromCircle::transverseIP(), reco::GhostTrackVertexFinder::vertexCompat(), reco::TrackGhostTrackState::vertexStateOnGhostTrack(), reco::TrackGhostTrackState::vertexStateOnMeasurement(), and ThirdHitPredictionFromCircle::HelixRZ::zAtR().

                           {
  return t * t;
}