trklet::Tracklet Class Reference

#include <Tracklet.h>

Public Member Functions
void	addMatch (unsigned int layerdisk, int ideltaphi, int ideltarz, double dphi, double drz, double dphiapprox, double drzapprox, int stubid, const trklet::Stub *stubptr)
std::string	addressstr ()
double	chiSqfit () const
double	d0 () const
double	d0approx () const
double	d0fit () const
double	d0fitexact () const
int	disk () const
const std::string	diskstubstr (const unsigned disk) const
bool	fit () const
const FPGAWord &	fpgad0 () const
const FPGAWord &	fpgaphi0 () const
const FPGAWord &	fpgarinv () const
const FPGAWord &	fpgat () const
const FPGAWord &	fpgaz0 () const
std::string	fullmatchdiskstr (int disk)
std::string	fullmatchstr (int layer)
int	getISeed () const
int	getITC () const
std::vector< const L1TStub * >	getL1Stubs ()
std::map< int, int >	getStubIDs ()
Track *	getTrack ()
FPGAWord	ichiSqfit () const
const FPGAWord &	id0fit () const
const Stub *	innerFPGAStub ()
const FPGAWord &	iphi0fit () const
const FPGAWord &	irinvfit () const
bool	isBarrel () const
int	isDisk () const
bool	isOverlap () const
const FPGAWord &	itfit () const
const FPGAWord &	iz0fit () const
int	layer () const
const std::string	layerstubstr (const unsigned layer) const
Track	makeTrack (const std::vector< const L1TStub * > &l1stubs)
bool	match (unsigned int layerdisk)
const Stub *	middleFPGAStub ()
const Stub *	outerFPGAStub ()
double	phi0 () const
double	phi0approx () const
double	phi0fit () const
double	phi0fitexact () const
Projection &	proj (int layerdisk)
unsigned int	PSseed () const
const Residual &	resid (unsigned int layerdisk)
double	rinv () const
double	rinvapprox () const
double	rinvfit () const
double	rinvfitexact () const
unsigned int	seedIndex () const
void	setFitPars (double rinvfit, double phi0fit, double d0fit, double tfit, double z0fit, double chisqrphifit, double chisqrzfit, double rinvfitexact, double phi0fitexact, double d0fitexact, double tfitexact, double z0fitexact, double chisqrphifitexact, double chisqrzfitexact, int irinvfit, int iphi0fit, int id0fit, int itfit, int iz0fit, int ichisqrphifit, int ichisqrzfit, int hitpattern, const std::vector< const L1TStub * > &l1stubs=std::vector< const L1TStub * >())
void	setTCIndex (int index)
void	setTrackIndex (int index)
void	setTrackletIndex (unsigned int index)
bool	stubtruthmatch (const L1TStub *stub)
double	t () const
double	tapprox () const
int	TCID () const
int	TCIndex () const
double	tfit () const
double	tfitexact () const
int	tpseed ()
std::string	trackfitstr () const
int	trackIndex () const
	Tracklet (Settings const &settings, unsigned int iSeed, const Stub *innerFPGAStub, const Stub *middleFPGAStub, const Stub *outerFPGAStub, double rinv, double phi0, double d0, double z0, double t, double rinvapprox, double phi0approx, double d0approx, double z0approx, double tapprox, int irinv, int iphi0, int id0, int iz0, int it, Projection projs[N_LAYER+N_DISK], bool disk, bool overlap=false)
int	trackletIndex () const
std::string	trackletparstr ()
std::string	trackletprojstr (int layer) const
std::string	trackletprojstrD (int disk) const
std::string	trackletprojstrdisk (int disk) const
std::string	trackletprojstrlayer (int layer) const
bool	validProj (int layerdisk) const
std::string	vmstrdisk (int disk, unsigned int allstubindex)
std::string	vmstrlayer (int layer, unsigned int allstubindex)
double	z0 () const
double	z0approx () const
double	z0fit () const
double	z0fitexact () const
	~Tracklet ()=default

Private Attributes
bool	barrel_
double	chisqrphifit_
double	chisqrphifitexact_
double	chisqrzfit_
double	chisqrzfitexact_
bool	disk_
TrackPars< double >	fitpars_
TrackPars< double >	fitparsexact_
TrackPars< FPGAWord >	fpgafitpars_
TrackPars< FPGAWord >	fpgapars_
std::unique_ptr< Track >	fpgatrack_
int	hitpattern_
FPGAWord	ichisqrphifit_
FPGAWord	ichisqrzfit_
const Stub *	innerFPGAStub_
const Stub *	middleFPGAStub_
const Stub *	outerFPGAStub_
bool	overlap_
Projection	proj_ [N_LAYER+N_DISK]
int	projdisk_ [N_DISK]
int	projlayer_ [N_LAYER-2]
Residual	resid_ [N_LAYER+N_DISK]
unsigned int	seedIndex_
Settings const &	settings_
int	TCIndex_
int	trackIndex_
int	trackletIndex_
TrackPars< double >	trackpars_
TrackPars< double >	trackparsapprox_
bool	triplet_

Detailed Description

Definition at line 27 of file Tracklet.h.

Constructor & Destructor Documentation

Tracklet::Tracklet	(	Settings const &	settings,
		unsigned int	iSeed,
		const Stub *	innerFPGAStub,
		const Stub *	middleFPGAStub,
		const Stub *	outerFPGAStub,
		double	rinv,
		double	phi0,
		double	d0,
		double	z0,
		double	t,
		double	rinvapprox,
		double	phi0approx,
		double	d0approx,
		double	z0approx,
		double	tapprox,
		int	irinv,
		int	iphi0,
		int	id0,
		int	iz0,
		int	it,
		Projection	projs[N_LAYER+N_DISK],
		bool	disk,
		bool	overlap = false
	)

Definition at line 16 of file Tracklet.cc.

References cms::cuda::assert(), barrel_, disk(), disk_, fpgapars_, fpgatrack_, mps_fire::i, ichisqrphifit_, ichisqrzfit_, trklet::TrackPars< T >::init(), innerFPGAStub(), innerFPGAStub_, trklet::Stub::l1tstub(), trklet::L1TStub::layer(), middleFPGAStub(), middleFPGAStub_, trklet::N_DISK, trklet::N_LAYER, trklet::Settings::nbitsd0(), trklet::Settings::nbitsphi0(), trklet::Settings::nbitsrinv(), trklet::Settings::nbitst(), trklet::Settings::nbitsz0(), outerFPGAStub(), outerFPGAStub_, muon::overlap(), overlap_, proj_, projdisk_, trklet::Settings::projdisks(), projlayer_, trklet::Settings::projlayers(), seedIndex_, trklet::FPGAWord::set(), settings_, TCIndex_, trackletIndex_, trackpars_, trackparsapprox_, and triplet_.

    : settings_(settings) {
  seedIndex_ = iSeed;

  overlap_ = overlap;
  disk_ = disk;
  assert(!(disk && overlap));
  barrel_ = (!disk) && (!overlap);
  triplet_ = false;

  trackletIndex_ = -1;
  TCIndex_ = -1;

  assert(disk_ || barrel_ || overlap_);

  if (barrel_ && middleFPGAStub == nullptr)
    assert(innerFPGAStub->l1tstub()->layer() < N_LAYER);

  innerFPGAStub_ = innerFPGAStub;
  middleFPGAStub_ = middleFPGAStub;
  outerFPGAStub_ = outerFPGAStub;

  trackpars_.init(rinv, phi0, d0, t, z0);

  trackparsapprox_.init(rinvapprox, phi0approx, d0approx, tapprox, z0approx);

  fpgapars_.rinv().set(irinv, settings_.nbitsrinv(), false, __LINE__, __FILE__);
  fpgapars_.phi0().set(iphi0, settings_.nbitsphi0(), false, __LINE__, __FILE__);
  fpgapars_.d0().set(id0, settings_.nbitsd0(), false, __LINE__, __FILE__);
  fpgapars_.z0().set(iz0, settings_.nbitsz0(), false, __LINE__, __FILE__);
  fpgapars_.t().set(it, settings_.nbitst(), false, __LINE__, __FILE__);

  fpgatrack_ = nullptr;

  triplet_ = (seedIndex_ >= 8);

  //fill projection layers
  for (unsigned int i = 0; i < N_LAYER - 2; i++) {
    projlayer_[i] = settings.projlayers(seedIndex_, i);
  }

  //fill projection disks
  for (unsigned int i = 0; i < N_DISK; i++) {
    projdisk_[i] = settings.projdisks(seedIndex_, i);
  }

  //Handle projections to the layers
  for (unsigned int i = 0; i < N_LAYER - 2; i++) {
    if (projlayer_[i] == 0)
      continue;
    if (!projs[projlayer_[i] - 1].valid())
      continue;

    proj_[projlayer_[i] - 1] = projs[projlayer_[i] - 1];
  }
  //Now handle projections to the disks
  for (unsigned int i = 0; i < N_DISK; i++) {
    if (projdisk_[i] == 0)
      continue;
    if (!projs[N_LAYER + projdisk_[i] - 1].valid())
      continue;

    proj_[N_LAYER + projdisk_[i] - 1] = projs[N_LAYER + projdisk_[i] - 1];
  }

  ichisqrphifit_.set(-1, 8, false);
  ichisqrzfit_.set(-1, 8, false);
}

trklet::Tracklet::~Tracklet ( )

default

Member Function Documentation

void Tracklet::addMatch	(	unsigned int	layerdisk,
		int	ideltaphi,
		int	ideltarz,
		double	dphi,
		double	drz,
		double	dphiapprox,
		double	drzapprox,
		int	stubid,
		const trklet::Stub *	stubptr
	)

Definition at line 295 of file Tracklet.cc.

References cms::cuda::assert(), trklet::Residual::init(), trklet::N_DISK, trklet::N_LAYER, resid_, and settings_.

Referenced by trklet::MatchCalculator::execute(), and trklet::MatchProcessor::matchCalculator().

                                                   {
  assert(layerdisk < N_LAYER + N_DISK);
  resid_[layerdisk].init(settings_, layerdisk, ideltaphi, ideltarz, stubid, dphi, drz, dphiapprox, drzapprox, stubptr);
}

std::string Tracklet::addressstr

(

)

Definition at line 177 of file Tracklet.cc.

References innerFPGAStub_, middleFPGAStub_, outerFPGAStub_, trklet::Stub::phiregionaddressstr(), str, and AlCaHLTBitMon_QueryRunRegistry::string.

                               {
  std::string str;
  str = innerFPGAStub_->phiregionaddressstr() + "|";
  if (middleFPGAStub_) {
    str += middleFPGAStub_->phiregionaddressstr() + "|";
  }
  str += outerFPGAStub_->phiregionaddressstr();

  return str;
}

double trklet::Tracklet::chiSqfit ( ) const

inline

Definition at line 143 of file Tracklet.h.

References chisqrphifit_, and chisqrzfit_.

{ return chisqrphifit_ + chisqrzfit_; }

double trklet::Tracklet::d0 ( ) const

inline

Definition at line 122 of file Tracklet.h.

References trklet::TrackPars< T >::d0(), and trackpars_.

Referenced by trklet::FitTrack::trackFitFake().

{ return trackpars_.d0(); }

double trklet::Tracklet::d0approx ( ) const

inline

Definition at line 128 of file Tracklet.h.

References trklet::TrackPars< T >::d0(), and trackparsapprox_.

Referenced by trklet::FitTrack::trackFitFake().

{ return trackparsapprox_.d0(); }

double trklet::Tracklet::d0fit ( ) const

inline

Definition at line 140 of file Tracklet.h.

References trklet::TrackPars< T >::d0(), and fitpars_.

{ return fitpars_.d0(); }

double trklet::Tracklet::d0fitexact ( ) const

inline

Definition at line 147 of file Tracklet.h.

References trklet::TrackPars< T >::d0(), and fitparsexact_.

{ return fitparsexact_.d0(); }

int Tracklet::disk

(

)

const

Definition at line 779 of file Tracklet.cc.

References funct::abs(), ztail::d, d1, trklet::L1TStub::disk(), innerFPGAStub_, trklet::Stub::l1tstub(), trklet::Stub::layerdisk(), middleFPGAStub_, min(), trklet::N_LAYER, and outerFPGAStub_.

Referenced by diskstubstr(), trklet::MatchCalculator::execute(), trklet::MatchProcessor::matchCalculator(), PSseed(), trklet::FitTrack::trackFitChisq(), and Tracklet().

                         {
  int d1 = (innerFPGAStub_ && (innerFPGAStub_->layerdisk() >= N_LAYER)) ? innerFPGAStub_->l1tstub()->disk() : 999,
      d2 = (middleFPGAStub_ && (middleFPGAStub_->layerdisk() >= N_LAYER)) ? middleFPGAStub_->l1tstub()->disk() : 999,
      d3 = (outerFPGAStub_ && (outerFPGAStub_->layerdisk() >= N_LAYER)) ? outerFPGAStub_->l1tstub()->disk() : 999,
      d = 999;
  if (abs(d1) < min(abs(d2), abs(d3)))
    d = d1;
  if (abs(d2) < min(abs(d1), abs(d3)))
    d = d2;
  if (abs(d3) < min(abs(d1), abs(d2)))
    d = d3;
  return (d < 999 ? d : 0);
}

const std::string Tracklet::diskstubstr

(

const unsigned

disk

)

const

Definition at line 595 of file Tracklet.cc.

References cms::cuda::assert(), gather_cfg::cout, disk(), trklet::Residual::fpgaphiresid(), trklet::Residual::fpgarzresid(), trklet::Residual::fpgastubid(), trklet::Stub::isPSmodule(), trklet::N_DISK, trklet::N_LAYER, trklet::Settings::nbitstrackletindex(), trklet::Settings::ntrackletmax(), trklet::Stub::r(), resid_, settings_, trklet::FPGAWord::str(), trklet::Residual::stubptr(), createJobs::tmp, and trackIndex_.

Referenced by trackfitstr().

                                                               {
  assert(disk < N_DISK);

  std::stringstream oss("");
  if (!resid_[N_LAYER + disk].valid())
    oss << "0|0000000|0000000000|000000000000|000000000000|0000000";
  else {
    if (trackIndex_ < 0 || trackIndex_ > (int)settings_.ntrackletmax()) {
      cout << "trackIndex_ = " << trackIndex_ << endl;
      assert(0);
    }
    const FPGAWord tmp(trackIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
    const FPGAWord& stubr = resid_[N_LAYER + disk].stubptr()->r();
    const bool isPS = resid_[N_LAYER + disk].stubptr()->isPSmodule();
    oss << "1|";  // valid bit
    oss << tmp.str() << "|";
    oss << resid_[N_LAYER + disk].fpgastubid().str() << "|";
    oss << (isPS ? stubr.str() : ("00000000" + stubr.str())) << "|";
    oss << resid_[N_LAYER + disk].fpgaphiresid().str() << "|";
    oss << resid_[N_LAYER + disk].fpgarzresid().str();
  }

  return oss.str();
}

bool trklet::Tracklet::fit ( ) const

inline

Definition at line 195 of file Tracklet.h.

References ichisqrphifit_, and trklet::FPGAWord::value().

Referenced by trklet::PurgeDuplicate::execute(), trklet::FitTrack::execute(), makeTrack(), and trackingPlots.Iteration::modules().

{ return ichisqrphifit_.value() != -1; }

const FPGAWord& trklet::Tracklet::fpgad0 ( ) const

inline

Definition at line 134 of file Tracklet.h.

References fpgapars_.

Referenced by trklet::FitTrack::trackFitFake().

{ return fpgapars_.d0(); }

const FPGAWord& trklet::Tracklet::fpgaphi0 ( ) const

inline

Definition at line 133 of file Tracklet.h.

References fpgapars_.

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return fpgapars_.phi0(); }

const FPGAWord& trklet::Tracklet::fpgarinv ( ) const

inline

Definition at line 132 of file Tracklet.h.

References fpgapars_.

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return fpgapars_.rinv(); }

const FPGAWord& trklet::Tracklet::fpgat ( ) const

inline

Definition at line 135 of file Tracklet.h.

References fpgapars_.

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return fpgapars_.t(); }

const FPGAWord& trklet::Tracklet::fpgaz0 ( ) const

inline

Definition at line 136 of file Tracklet.h.

References fpgapars_.

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return fpgapars_.z0(); }

std::string Tracklet::fullmatchdiskstr

(

int

disk

)

Definition at line 324 of file Tracklet.cc.

References cms::cuda::assert(), Exception, trklet::Residual::fpgaphiresid(), trklet::Residual::fpgarzresid(), trklet::Residual::fpgastubid(), trklet::Stub::isPSmodule(), trklet::N_DISK, trklet::N_LAYER, trklet::Settings::nbitstcindex(), trklet::Settings::nbitstrackletindex(), trklet::Settings::ntrackletmax(), trklet::Stub::r(), resid_, trklet::FPGAWord::set(), settings_, trklet::FPGAWord::str(), AlCaHLTBitMon_QueryRunRegistry::string, trklet::Residual::stubptr(), TCIndex_, createJobs::tmp, and trackletIndex_.

                                             {
  assert(disk > 0 && disk <= N_DISK);

  FPGAWord tmp;
  if (trackletIndex_ < 0 || trackletIndex_ > (int)settings_.ntrackletmax()) {
    throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__ << " trackletIndex_ = " << trackletIndex_;
  }
  tmp.set(trackletIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
  FPGAWord tcid;
  tcid.set(TCIndex_, settings_.nbitstcindex(), true, __LINE__, __FILE__);
  const FPGAWord& stubr = resid_[N_LAYER + disk - 1].stubptr()->r();
  const bool isPS = resid_[N_LAYER + disk - 1].stubptr()->isPSmodule();
  std::string oss = tcid.str() + "|" + tmp.str() + "|" + resid_[N_LAYER + disk - 1].fpgastubid().str() + "|" +
                    (isPS ? stubr.str() : ("00000000" + stubr.str())) + "|" +
                    resid_[N_LAYER + disk - 1].fpgaphiresid().str() + "|" +
                    resid_[N_LAYER + disk - 1].fpgarzresid().str();
  return oss;
}

std::string Tracklet::fullmatchstr

(

int

layer

)

Definition at line 308 of file Tracklet.cc.

References cms::cuda::assert(), Exception, trklet::Residual::fpgaphiresid(), trklet::Residual::fpgarzresid(), trklet::Residual::fpgastubid(), trklet::N_LAYER, trklet::Settings::nbitstcindex(), trklet::Settings::nbitstrackletindex(), trklet::Settings::ntrackletmax(), trklet::Stub::r(), resid_, trklet::FPGAWord::set(), settings_, trklet::FPGAWord::str(), AlCaHLTBitMon_QueryRunRegistry::string, trklet::Residual::stubptr(), TCIndex_, createJobs::tmp, and trackletIndex_.

                                          {
  assert(layer > 0 && layer <= N_LAYER);

  FPGAWord tmp;
  if (trackletIndex_ < 0 || trackletIndex_ > (int)settings_.ntrackletmax()) {
    throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__ << " trackletIndex_ = " << trackletIndex_;
  }
  tmp.set(trackletIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
  FPGAWord tcid;
  tcid.set(TCIndex_, settings_.nbitstcindex(), true, __LINE__, __FILE__);
  std::string oss = tcid.str() + "|" + tmp.str() + "|" + resid_[layer - 1].fpgastubid().str() + "|" +
                    resid_[layer - 1].stubptr()->r().str() + "|" + resid_[layer - 1].fpgaphiresid().str() + "|" +
                    resid_[layer - 1].fpgarzresid().str();
  return oss;
}

int Tracklet::getISeed

(

)

const

Definition at line 798 of file Tracklet.cc.

References cms::cuda::assert(), trklet::N_SEED, trklet::Settings::nbitsitc(), settings_, and TCIndex_.

Referenced by trklet::TrackletCalculatorBase::barrelSeeding(), trklet::TrackletCalculatorDisplaced::DDLSeeding(), trklet::TrackletCalculatorBase::diskSeeding(), trklet::MatchEngine::execute(), trklet::MatchCalculator::execute(), trklet::FitTrack::execute(), getITC(), trklet::TrackletCalculatorDisplaced::LLDSeeding(), trklet::TrackletCalculatorDisplaced::LLLSeeding(), makeTrack(), trklet::MatchProcessor::matchCalculator(), trklet::TrackletCalculatorBase::overlapSeeding(), and trackfitstr().

                             {
  const int iSeed = TCIndex_ >> settings_.nbitsitc();
  assert(iSeed >= 0 && iSeed <= (int)N_SEED);
  return iSeed;
}

int Tracklet::getITC

(

)

const

Definition at line 804 of file Tracklet.cc.

References cms::cuda::assert(), getISeed(), trklet::Settings::nbitsitc(), settings_, and TCIndex_.

                           {
  const int iSeed = getISeed(), iTC = TCIndex_ - (iSeed << settings_.nbitsitc());
  assert(iTC >= 0 && iTC <= 14);
  return iTC;
}

std::vector< const L1TStub * > Tracklet::getL1Stubs

(

)

Definition at line 343 of file Tracklet.cc.

References innerFPGAStub_, trklet::Stub::l1tstub(), middleFPGAStub_, outerFPGAStub_, resid_, and createJobs::tmp.

Referenced by makeTrack().

                                               {
  std::vector<const L1TStub*> tmp;

  if (innerFPGAStub_)
    tmp.push_back(innerFPGAStub_->l1tstub());
  if (middleFPGAStub_)
    tmp.push_back(middleFPGAStub_->l1tstub());
  if (outerFPGAStub_)
    tmp.push_back(outerFPGAStub_->l1tstub());

  for (const auto& iresid : resid_) {
    if (iresid.valid())
      tmp.push_back(iresid.stubptr()->l1tstub());
  }

  return tmp;
}

std::map< int, int > Tracklet::getStubIDs

(

)

Definition at line 361 of file Tracklet.cc.

References cms::cuda::assert(), barrel_, trklet::L1TStub::disk(), trklet::Stub::disk(), disk_, trklet::Residual::fpgastubid(), mps_fire::i, innerFPGAStub_, itfit(), trklet::Stub::l1tstub(), trklet::Stub::layer(), middleFPGAStub_, trklet::N_BITSMEMADDRESS, trklet::N_DISK, trklet::N_LAYER, trklet::FPGAWord::nbits(), outerFPGAStub_, overlap_, trklet::Stub::phiregionaddress(), resid_, trklet::Stub::stubindex(), trklet::Residual::valid(), trklet::FPGAWord::value(), and relativeConstraints::value.

Referenced by makeTrack().

                                      {
  std::map<int, int> stubIDs;

  // For future reference, *resid_[i] uses i as the absolute stub index. (0-5 for barrel, 0-4 for disk)
  // On the other hand, proj*_[i] uses i almost like *resid_[i], except the *seeding* layer indices are removed entirely.
  // E.g. An L3L4 track has 0=L1, 1=L2, 2=L4, 3=L5 for the barrels (for proj*_[i])

  if (innerFPGAStub_)
    assert(innerFPGAStub_->stubindex().nbits() == N_BITSMEMADDRESS);
  if (middleFPGAStub_)
    assert(middleFPGAStub_->stubindex().nbits() == N_BITSMEMADDRESS);
  if (outerFPGAStub_)
    assert(outerFPGAStub_->stubindex().nbits() == N_BITSMEMADDRESS);

  if (barrel_) {
    for (int i = 0; i < N_LAYER; i++) {
      //check barrel
      if (resid_[i].valid()) {
        // two extra bits to indicate if the matched stub is local or from neighbor
        int location = 1;  // local
        location <<= resid_[i].fpgastubid().nbits();

        stubIDs[1 + i] = resid_[i].fpgastubid().value() + location;
      }

      //check disk
      if (i >= N_DISK)
        continue;  //i=[0..4] for disks
      if (resid_[N_LAYER + i].valid()) {
        if (i == 3 && resid_[0].valid() && innerFPGAStub_->layer().value() == 1)
          continue;  // Don't add D4 if track has L1 stub
        // two extra bits to indicate if the matched stub is local or from neighbor
        int location = 1;  // local
        location <<= resid_[N_LAYER + i].fpgastubid().nbits();

        if (itfit().value() < 0) {
          stubIDs[-(N_LAYER + N_DISK) - i] = resid_[N_LAYER + i].fpgastubid().value() + location;
        } else {
          stubIDs[N_LAYER + N_DISK + i] = resid_[N_LAYER + i].fpgastubid().value() + location;
        }
      }
    }

    //get stubs making up tracklet
    if (innerFPGAStub_)
      stubIDs[innerFPGAStub_->layer().value() + 1] = innerFPGAStub_->phiregionaddress() + (1 << 10);
    if (middleFPGAStub_)
      stubIDs[middleFPGAStub_->layer().value() + 1] = middleFPGAStub_->phiregionaddress() + (1 << 10);
    if (outerFPGAStub_)
      stubIDs[outerFPGAStub_->layer().value() + 1] = outerFPGAStub_->phiregionaddress() + (1 << 10);

  } else if (disk_) {
    for (int i = 0; i < N_DISK; i++) {
      //check barrel
      if (resid_[i].valid()) {
        // two extra bits to indicate if the matched stub is local or from neighbor
        int location = 1;  // local
        location <<= resid_[i].fpgastubid().nbits();

        stubIDs[1 + i] = resid_[i].fpgastubid().value() + location;
      }

      //check disks
      if (i == 4 && resid_[1].valid())
        continue;  // Don't add D5 if track has L2 stub
      if (resid_[N_LAYER + i].valid()) {
        // two extra bits to indicate if the matched stub is local or from neighbor
        int location = 1;  // local
        location <<= resid_[N_LAYER + i].fpgastubid().nbits();

        if (innerFPGAStub_->l1tstub()->disk() < 0) {
          stubIDs[-11 - i] = resid_[N_LAYER + i].fpgastubid().value() + location;
        } else {
          stubIDs[11 + i] = resid_[N_LAYER + i].fpgastubid().value() + location;
        }
      }
    }

    //get stubs making up tracklet
    if (innerFPGAStub_->disk().value() < 0) {  //negative side runs 6-10
      if (innerFPGAStub_)
        stubIDs[innerFPGAStub_->disk().value() - 10] = innerFPGAStub_->phiregionaddress() + (1 << 10);
      if (middleFPGAStub_)
        stubIDs[middleFPGAStub_->disk().value() - 10] = middleFPGAStub_->phiregionaddress() + (1 << 10);
      if (outerFPGAStub_)
        stubIDs[outerFPGAStub_->disk().value() - 10] = outerFPGAStub_->phiregionaddress() + (1 << 10);
    } else {  // positive side runs 11-15]
      if (innerFPGAStub_)
        stubIDs[innerFPGAStub_->disk().value() + 10] = innerFPGAStub_->phiregionaddress() + (1 << 10);
      if (middleFPGAStub_)
        stubIDs[middleFPGAStub_->disk().value() + 10] = middleFPGAStub_->phiregionaddress() + (1 << 10);
      if (outerFPGAStub_)
        stubIDs[outerFPGAStub_->disk().value() + 10] = outerFPGAStub_->phiregionaddress() + (1 << 10);
    }

  } else if (overlap_) {
    for (int i = 0; i < N_DISK; i++) {
      //check barrel
      if (resid_[i].valid()) {
        // two extra bits to indicate if the matched stub is local or from neighbor
        int location = 1;  // local
        location <<= resid_[i].fpgastubid().nbits();

        stubIDs[1 + i] = resid_[i].fpgastubid().value() + location;
      }

      //check disks
      if (resid_[N_LAYER + i].valid()) {
        // two extra bits to indicate if the matched stub is local or from neighbor
        int location = 1;  // local
        location <<= resid_[N_LAYER + i].fpgastubid().nbits();

        if (innerFPGAStub_->l1tstub()->disk() < 0) {  // if negative overlap
          if (innerFPGAStub_->layer().value() != 2 || !resid_[0].valid() ||
              i != 3) {  // Don't add D4 if this is an L3L2 track with an L1 stub
            stubIDs[-11 - i] = resid_[N_LAYER + i].fpgastubid().value() + location;
          }
        } else {
          if (innerFPGAStub_->layer().value() != 2 || !resid_[0].valid() || i != 3) {
            stubIDs[11 + i] = resid_[N_LAYER + i].fpgastubid().value() + location;
          }
        }
      }
    }

    //get stubs making up tracklet

    if (innerFPGAStub_->layer().value() == 2) {  // L3L2 track
      if (innerFPGAStub_)
        stubIDs[innerFPGAStub_->layer().value() + 1] = innerFPGAStub_->phiregionaddress() + (1 << 10);
      if (middleFPGAStub_)
        stubIDs[middleFPGAStub_->layer().value() + 1] = middleFPGAStub_->phiregionaddress() + (1 << 10);
      if (outerFPGAStub_)
        stubIDs[outerFPGAStub_->layer().value() + 1] = outerFPGAStub_->phiregionaddress() + (1 << 10);
    } else if (innerFPGAStub_->disk().value() < 0) {  //negative side runs -11 - -15
      if (innerFPGAStub_)
        stubIDs[innerFPGAStub_->disk().value() - 10] = innerFPGAStub_->phiregionaddress() + (1 << 10);
      if (middleFPGAStub_)
        stubIDs[middleFPGAStub_->layer().value() + 1] = middleFPGAStub_->phiregionaddress() + (1 << 10);
      if (outerFPGAStub_)
        stubIDs[outerFPGAStub_->layer().value() + 1] = outerFPGAStub_->phiregionaddress() + (1 << 10);
    } else {  // positive side runs 11-15]
      if (innerFPGAStub_)
        stubIDs[innerFPGAStub_->disk().value() + 10] = innerFPGAStub_->phiregionaddress() + (1 << 10);
      if (middleFPGAStub_)
        stubIDs[middleFPGAStub_->layer().value() + 1] = middleFPGAStub_->phiregionaddress() + (1 << 10);
      if (outerFPGAStub_)
        stubIDs[outerFPGAStub_->layer().value() + 1] = outerFPGAStub_->phiregionaddress() + (1 << 10);
    }
  }

  return stubIDs;
}

Track* trklet::Tracklet::getTrack ( )

inline

Definition at line 190 of file Tracklet.h.

References cms::cuda::assert(), and fpgatrack_.

Referenced by trklet::PurgeDuplicate::execute().

                      {
      assert(fpgatrack_ != nullptr);
      return fpgatrack_.get();
    }

FPGAWord trklet::Tracklet::ichiSqfit ( ) const

inline

Definition at line 156 of file Tracklet.h.

References ichisqrphifit_, ichisqrzfit_, trklet::FPGAWord::nbits(), and trklet::FPGAWord::value().

                               {
      return FPGAWord(ichisqrphifit_.value() + ichisqrzfit_.value(), ichisqrphifit_.nbits());
    }

const FPGAWord& trklet::Tracklet::id0fit ( ) const

inline

Definition at line 153 of file Tracklet.h.

References fpgafitpars_.

{ return fpgafitpars_.d0(); }

const Stub* trklet::Tracklet::innerFPGAStub ( )

inline

Definition at line 61 of file Tracklet.h.

References innerFPGAStub_.

Referenced by trklet::FitTrack::trackFitChisq(), and Tracklet().

{ return innerFPGAStub_; }

const FPGAWord& trklet::Tracklet::iphi0fit ( ) const

inline

Definition at line 152 of file Tracklet.h.

References fpgafitpars_.

{ return fpgafitpars_.phi0(); }

const FPGAWord& trklet::Tracklet::irinvfit ( ) const

inline

Definition at line 151 of file Tracklet.h.

References fpgafitpars_.

{ return fpgafitpars_.rinv(); }

bool trklet::Tracklet::isBarrel ( ) const

inline

Definition at line 200 of file Tracklet.h.

References barrel_.

Referenced by trklet::FitTrack::trackFitChisq().

{ return barrel_; }

int trklet::Tracklet::isDisk ( ) const

inline

Definition at line 202 of file Tracklet.h.

References disk_.

Referenced by trklet::FitTrack::trackFitChisq().

{ return disk_; }

bool trklet::Tracklet::isOverlap ( ) const

inline

Definition at line 201 of file Tracklet.h.

References overlap_.

Referenced by trklet::FitTrack::trackFitChisq().

{ return overlap_; }

const FPGAWord& trklet::Tracklet::itfit ( ) const

inline

Definition at line 154 of file Tracklet.h.

References fpgafitpars_.

Referenced by getStubIDs().

{ return fpgafitpars_.t(); }

const FPGAWord& trklet::Tracklet::iz0fit ( ) const

inline

Definition at line 155 of file Tracklet.h.

References fpgafitpars_.

{ return fpgafitpars_.z0(); }

int Tracklet::layer

(

)

const

Definition at line 770 of file Tracklet.cc.

References innerFPGAStub_, cmsLHEtoEOSManager::l, trklet::Stub::l1tstub(), trklet::Stub::layerdisk(), trklet::L1TStub::layerdisk(), middleFPGAStub_, min(), trklet::N_LAYER, and outerFPGAStub_.

Referenced by trklet::MatchCalculator::execute(), geometryXMLparser.DTAlignable::index(), geometryXMLparser.CSCAlignable::index(), layerstubstr(), trklet::MatchProcessor::matchCalculator(), PSseed(), and trklet::FitTrack::trackFitChisq().

                          {
  int l1 = (innerFPGAStub_ && innerFPGAStub_->layerdisk() < N_LAYER) ? innerFPGAStub_->l1tstub()->layerdisk() + 1 : 999,
      l2 = (middleFPGAStub_ && middleFPGAStub_->layerdisk() < N_LAYER) ? middleFPGAStub_->l1tstub()->layerdisk() + 1
                                                                       : 999,
      l3 = (outerFPGAStub_ && outerFPGAStub_->layerdisk() < N_LAYER) ? outerFPGAStub_->l1tstub()->layerdisk() + 1 : 999,
      l = min(min(l1, l2), l3);
  return (l < 999 ? l : 0);
}

const std::string Tracklet::layerstubstr

(

const unsigned

layer

)

const

Definition at line 572 of file Tracklet.cc.

References cms::cuda::assert(), gather_cfg::cout, trklet::Residual::fpgaphiresid(), trklet::Residual::fpgarzresid(), trklet::Residual::fpgastubid(), layer(), trklet::N_LAYER, trklet::Settings::nbitstrackletindex(), trklet::Settings::ntrackletmax(), trklet::Stub::r(), resid_, settings_, trklet::FPGAWord::str(), trklet::Residual::stubptr(), createJobs::tmp, and trackIndex_.

Referenced by trackfitstr().

                                                                 {
  assert(layer < N_LAYER);

  std::stringstream oss("");
  if (!resid_[layer].valid())
    oss << "0|0000000|0000000000|0000000|000000000000|000000000";
  else {
    if (trackIndex_ < 0 || trackIndex_ > (int)settings_.ntrackletmax()) {
      cout << "trackIndex_ = " << trackIndex_ << endl;
      assert(0);
    }
    const FPGAWord tmp(trackIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
    oss << "1|";  // valid bit
    oss << tmp.str() << "|";
    oss << resid_[layer].fpgastubid().str() << "|";
    oss << resid_[layer].stubptr()->r().str() << "|";
    oss << resid_[layer].fpgaphiresid().str() << "|";
    oss << resid_[layer].fpgarzresid().str();
  }

  return oss.str();
}

Track Tracklet::makeTrack

(

const std::vector< const L1TStub * > &

l1stubs

)

Definition at line 738 of file Tracklet.cc.

References cms::cuda::assert(), chisqrphifit_, chisqrzfit_, fit(), fpgafitpars_, getISeed(), getL1Stubs(), getStubIDs(), hitpattern_, ichisqrphifit_, ichisqrzfit_, createJobs::tmp, and trklet::FPGAWord::value().

Referenced by setFitPars().

                                                               {
  assert(fit());

  TrackPars<int> ipars(fpgafitpars_.rinv().value(),
                       fpgafitpars_.phi0().value(),
                       fpgafitpars_.d0().value(),
                       fpgafitpars_.t().value(),
                       fpgafitpars_.z0().value());

  // If fitter produced no stub list, take it from original tracklet.
  vector<const L1TStub*> tmp = l1stubs.empty() ? getL1Stubs() : l1stubs;

  vector<L1TStub> tmp2;

  tmp2.reserve(tmp.size());
  for (auto stub : tmp) {
    tmp2.push_back(*stub);
  }

  Track tmpTrack(ipars,
                 ichisqrphifit_.value(),
                 ichisqrzfit_.value(),
                 chisqrphifit_,
                 chisqrzfit_,
                 hitpattern_,
                 getStubIDs(),
                 tmp2,
                 getISeed());

  return tmpTrack;
}

bool trklet::Tracklet::match ( unsigned int  layerdisk )

inline

Definition at line 105 of file Tracklet.h.

References cms::cuda::assert(), trklet::N_DISK, trklet::N_LAYER, resid_, and trklet::Residual::valid().

Referenced by trklet::FitTrack::trackFitChisq().

                                       {
      assert(layerdisk < N_LAYER + N_DISK);
      return resid_[layerdisk].valid();
    }

const Stub* trklet::Tracklet::middleFPGAStub ( )

inline

Definition at line 63 of file Tracklet.h.

References middleFPGAStub_.

Referenced by Tracklet().

{ return middleFPGAStub_; }

const Stub* trklet::Tracklet::outerFPGAStub ( )

inline

Definition at line 65 of file Tracklet.h.

References outerFPGAStub_.

Referenced by trklet::FitTrack::trackFitChisq(), and Tracklet().

{ return outerFPGAStub_; }

double trklet::Tracklet::phi0 ( ) const

inline

Definition at line 121 of file Tracklet.h.

References trklet::TrackPars< T >::phi0(), and trackpars_.

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return trackpars_.phi0(); }

double trklet::Tracklet::phi0approx ( ) const

inline

Definition at line 127 of file Tracklet.h.

References trklet::TrackPars< T >::phi0(), and trackparsapprox_.

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return trackparsapprox_.phi0(); }

double trklet::Tracklet::phi0fit ( ) const

inline

Definition at line 139 of file Tracklet.h.

References fitpars_, and trklet::TrackPars< T >::phi0().

{ return fitpars_.phi0(); }

double trklet::Tracklet::phi0fitexact ( ) const

inline

Definition at line 146 of file Tracklet.h.

References fitparsexact_, and trklet::TrackPars< T >::phi0().

{ return fitparsexact_.phi0(); }

Projection& trklet::Tracklet::proj ( int  layerdisk )

inline

Definition at line 87 of file Tracklet.h.

References cms::cuda::assert(), proj_, and validProj().

Referenced by trklet::TrackletCalculatorDisplaced::addDiskProj(), trklet::TrackletCalculatorBase::addDiskProj(), trklet::TrackletCalculatorDisplaced::addLayerProj(), trklet::TrackletCalculatorBase::addLayerProj(), trklet::ProjectionRouter::execute(), trklet::MatchEngine::execute(), trklet::MatchCalculator::execute(), trklet::MatchProcessor::execute(), trklet::PurgeDuplicate::getPhiRes(), and trklet::MatchProcessor::matchCalculator().

                                    {
      assert(validProj(layerdisk));
      return proj_[layerdisk];
    }

unsigned int trklet::Tracklet::PSseed ( ) const

inline

Definition at line 220 of file Tracklet.h.

References disk(), and layer().

Referenced by trklet::MatchEngine::execute(), trklet::MatchProcessor::execute(), and vmstrlayer().

{ return ((layer() == 1) || (layer() == 2) || (disk() != 0)) ? 1 : 0; }

const Residual& trklet::Tracklet::resid ( unsigned int  layerdisk )

inline

Definition at line 110 of file Tracklet.h.

References cms::cuda::assert(), trklet::N_DISK, trklet::N_LAYER, and resid_.

Referenced by trklet::FitTrack::trackFitChisq().

                                                  {
      assert(layerdisk < N_LAYER + N_DISK);
      assert(resid_[layerdisk].valid());
      return resid_[layerdisk];
    }

double trklet::Tracklet::rinv ( ) const

inline

Definition at line 120 of file Tracklet.h.

References trklet::TrackPars< T >::rinv(), and trackpars_.

Referenced by trklet::MatchCalculator::execute(), trklet::MatchProcessor::matchCalculator(), trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return trackpars_.rinv(); }

double trklet::Tracklet::rinvapprox ( ) const

inline

Definition at line 126 of file Tracklet.h.

References trklet::TrackPars< T >::rinv(), and trackparsapprox_.

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return trackparsapprox_.rinv(); }

double trklet::Tracklet::rinvfit ( ) const

inline

Definition at line 138 of file Tracklet.h.

References fitpars_, and trklet::TrackPars< T >::rinv().

{ return fitpars_.rinv(); }

double trklet::Tracklet::rinvfitexact ( ) const

inline

Definition at line 145 of file Tracklet.h.

References fitparsexact_, and trklet::TrackPars< T >::rinv().

{ return fitparsexact_.rinv(); }

unsigned int trklet::Tracklet::seedIndex ( ) const

inline

Definition at line 222 of file Tracklet.h.

References seedIndex_.

Referenced by trklet::PurgeDuplicate::execute(), trklet::PurgeDuplicate::getPhiRes(), and trackfitstr().

{ return seedIndex_; }

void Tracklet::setFitPars	(	double	rinvfit,
		double	phi0fit,
		double	d0fit,
		double	tfit,
		double	z0fit,
		double	chisqrphifit,
		double	chisqrzfit,
		double	rinvfitexact,
		double	phi0fitexact,
		double	d0fitexact,
		double	tfitexact,
		double	z0fitexact,
		double	chisqrphifitexact,
		double	chisqrzfitexact,
		int	irinvfit,
		int	iphi0fit,
		int	id0fit,
		int	itfit,
		int	iz0fit,
		int	ichisqrphifit,
		int	ichisqrzfit,
		int	hitpattern,
		const std::vector< const L1TStub * > &	l1stubs = std::vector<const L1TStub*>()
	)

Definition at line 515 of file Tracklet.cc.

References chisqrphifit_, chisqrphifitexact_, chisqrzfit_, chisqrzfitexact_, fitpars_, fitparsexact_, fpgafitpars_, fpgatrack_, hitpattern_, ichisqrphifit_, ichisqrzfit_, trklet::TrackPars< T >::init(), makeTrack(), trklet::Settings::nbitsz0(), trklet::FPGAWord::set(), and settings_.

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

                                                                 {
  fitpars_.init(rinvfit, phi0fit, d0fit, tfit, z0fit);
  chisqrphifit_ = chisqrphifit;
  chisqrzfit_ = chisqrzfit;

  fitparsexact_.init(rinvfitexact, phi0fitexact, d0fitexact, tfitexact, z0fitexact);
  chisqrphifitexact_ = chisqrphifitexact;
  chisqrzfitexact_ = chisqrzfitexact;

  if (irinvfit > (1 << 14))
    irinvfit = (1 << 14);
  if (irinvfit <= -(1 << 14))
    irinvfit = -(1 << 14) + 1;
  fpgafitpars_.rinv().set(irinvfit, 15, false, __LINE__, __FILE__);
  fpgafitpars_.phi0().set(iphi0fit, 19, false, __LINE__, __FILE__);
  fpgafitpars_.d0().set(id0fit, 19, false, __LINE__, __FILE__);
  fpgafitpars_.t().set(itfit, 14, false, __LINE__, __FILE__);

  if (iz0fit >= (1 << (settings_.nbitsz0() - 1))) {
    iz0fit = (1 << (settings_.nbitsz0() - 1)) - 1;
  }

  if (iz0fit <= -(1 << (settings_.nbitsz0() - 1))) {
    iz0fit = 1 - (1 << (settings_.nbitsz0() - 1));
  }

  fpgafitpars_.z0().set(iz0fit, settings_.nbitsz0(), false, __LINE__, __FILE__);
  ichisqrphifit_.set(ichisqrphifit, 8, true, __LINE__, __FILE__);
  ichisqrzfit_.set(ichisqrzfit, 8, true, __LINE__, __FILE__);

  hitpattern_ = hitpattern;

  fpgatrack_ = std::make_unique<Track>(makeTrack(l1stubs));
}

void trklet::Tracklet::setTCIndex ( int  index )

inline

Definition at line 208 of file Tracklet.h.

References TCIndex_.

Referenced by trklet::TrackletCalculatorBase::barrelSeeding(), trklet::TrackletCalculatorDisplaced::DDLSeeding(), trklet::TrackletCalculatorBase::diskSeeding(), trklet::TrackletCalculatorDisplaced::LLDSeeding(), trklet::TrackletCalculatorDisplaced::LLLSeeding(), and trklet::TrackletCalculatorBase::overlapSeeding().

{ TCIndex_ = index; }

void Tracklet::setTrackIndex

(

int

index

)

Definition at line 810 of file Tracklet.cc.

References cms::cuda::assert(), trklet::Settings::ntrackletmax(), settings_, and trackIndex_.

Referenced by trklet::FitTrack::execute().

                                      {
  trackIndex_ = index;
  assert(index <= (int)settings_.ntrackletmax());
}

void Tracklet::setTrackletIndex

(

unsigned int

index

)

Definition at line 793 of file Tracklet.cc.

References cms::cuda::assert(), trklet::Settings::ntrackletmax(), settings_, and trackletIndex_.

Referenced by trklet::TrackletCalculatorBase::barrelSeeding(), trklet::TrackletCalculatorDisplaced::DDLSeeding(), trklet::TrackletCalculatorBase::diskSeeding(), trklet::TrackletCalculatorDisplaced::LLDSeeding(), trklet::TrackletCalculatorDisplaced::LLLSeeding(), and trklet::TrackletCalculatorBase::overlapSeeding().

                                                  {
  trackletIndex_ = index;
  assert(index <= settings_.ntrackletmax());
}

bool Tracklet::stubtruthmatch

(

const L1TStub *

stub

)

Definition at line 138 of file Tracklet.cc.

References funct::abs(), innerFPGAStub_, trklet::Stub::l1tstub(), outerFPGAStub_, cmsswSequenceInfo::tp, and trklet::L1TStub::tps().

Referenced by trklet::MatchCalculator::execute(), and trklet::MatchProcessor::matchCalculator().

                                                 {
  set<int> tpset;
  set<int> tpsetstub;
  set<int> tpsetstubinner;
  set<int> tpsetstubouter;

  vector<int> tps = stub->tps();
  for (auto tp : tps) {
    if (tp != 0) {
      tpsetstub.insert(tp);
      tpset.insert(abs(tp));
    }
  }
  tps = innerFPGAStub_->l1tstub()->tps();
  for (auto tp : tps) {
    if (tp != 0) {
      tpsetstubinner.insert(tp);
      tpset.insert(abs(tp));
    }
  }
  tps = outerFPGAStub_->l1tstub()->tps();
  for (auto tp : tps) {
    if (tp != 0) {
      tpsetstubouter.insert(tp);
      tpset.insert(abs(tp));
    }
  }

  for (auto tp : tpset) {
    if (tpsetstub.find(tp) != tpsetstub.end() && tpsetstub.find(-tp) != tpsetstub.end() &&
        tpsetstubinner.find(tp) != tpsetstubinner.end() && tpsetstubinner.find(-tp) != tpsetstubinner.end() &&
        tpsetstubouter.find(tp) != tpsetstubouter.end() && tpsetstubouter.find(-tp) != tpsetstubouter.end()) {
      return true;
    }
  }

  return false;
}

double trklet::Tracklet::t ( ) const

inline

Definition at line 123 of file Tracklet.h.

References trklet::TrackPars< T >::t(), and trackpars_.

Referenced by trklet::MatchCalculator::execute(), trklet::FitTrack::execute(), trklet::MatchProcessor::matchCalculator(), trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return trackpars_.t(); }

double trklet::Tracklet::tapprox ( ) const

inline

Definition at line 129 of file Tracklet.h.

References trklet::TrackPars< T >::t(), and trackparsapprox_.

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return trackparsapprox_.t(); }

int trklet::Tracklet::TCID ( ) const

inline

Definition at line 212 of file Tracklet.h.

References trklet::Settings::nbitstrackletindex(), settings_, TCIndex_, and trackletIndex_.

Referenced by trklet::FullMatchMemory::addMatch(), trklet::VMProjectionsMemory::addTracklet(), trklet::ProjectionRouter::execute(), trklet::MatchCalculator::execute(), trklet::MatchProcessor::execute(), trklet::FitTrack::execute(), and trklet::MatchEngineUnit::TCID().

{ return TCIndex_ * (1 << settings_.nbitstrackletindex()) + trackletIndex_; }

int trklet::Tracklet::TCIndex ( ) const

inline

Definition at line 210 of file Tracklet.h.

References TCIndex_.

Referenced by trklet::MatchEngineUnit::step().

{ return TCIndex_; }

double trklet::Tracklet::tfit ( ) const

inline

Definition at line 141 of file Tracklet.h.

References fitpars_, and trklet::TrackPars< T >::t().

{ return fitpars_.t(); }

double trklet::Tracklet::tfitexact ( ) const

inline

Definition at line 148 of file Tracklet.h.

References fitparsexact_, and trklet::TrackPars< T >::t().

{ return fitparsexact_.t(); }

int Tracklet::tpseed

(

)

Definition at line 107 of file Tracklet.cc.

References funct::abs(), innerFPGAStub_, trklet::Stub::l1tstub(), outerFPGAStub_, cmsswSequenceInfo::tp, and trklet::L1TStub::tps().

Referenced by trklet::TrackletCalculatorBase::barrelSeeding().

                     {
  set<int> tpset;

  set<int> tpsetstubinner;
  set<int> tpsetstubouter;

  vector<int> tps = innerFPGAStub_->l1tstub()->tps();
  for (auto tp : tps) {
    if (tp != 0) {
      tpsetstubinner.insert(tp);
      tpset.insert(abs(tp));
    }
  }

  tps = outerFPGAStub_->l1tstub()->tps();
  for (auto tp : tps) {
    if (tp != 0) {
      tpsetstubouter.insert(tp);
      tpset.insert(abs(tp));
    }
  }

  for (auto& tp : tpset) {
    if (tpsetstubinner.find(tp) != tpsetstubinner.end() && tpsetstubinner.find(-tp) != tpsetstubinner.end() &&
        tpsetstubouter.find(tp) != tpsetstubouter.end() && tpsetstubouter.find(-tp) != tpsetstubouter.end()) {
      return tp;
    }
  }
  return 0;
}

std::string Tracklet::trackfitstr

(

)

const

Definition at line 620 of file Tracklet.cc.

References diskstubstr(), fpgapars_, getISeed(), mps_fire::i, layerstubstr(), trklet::Settings::nbitsseed(), seedIndex(), settings_, AlCaHLTBitMon_QueryRunRegistry::string, createJobs::tmp, and trklet::Settings::writeoutReal().

                                      {
  const unsigned maxNHits = 8;
  const unsigned nBitsPerHit = 3;
  vector<string> stub(maxNHits, "0");
  string hitmap(maxNHits * nBitsPerHit, '0');

  // Assign stub strings for each of the possible projections for each seed.
  // The specific layers/disks for a given seed are determined by the wiring.
  switch (seedIndex()) {
    case 0:                       // L1L2
      stub[0] = layerstubstr(2);  // L3
      stub[1] = layerstubstr(3);  // L4
      stub[2] = layerstubstr(4);  // L5
      stub[3] = layerstubstr(5);  // L6

      stub[4] = diskstubstr(0);  // D1
      stub[5] = diskstubstr(1);  // D2
      stub[6] = diskstubstr(2);  // D3
      stub[7] = diskstubstr(3);  // D4

      break;

    case 1:                       // L2L3
      stub[0] = layerstubstr(0);  // L1
      stub[1] = layerstubstr(3);  // L4
      stub[2] = layerstubstr(4);  // L5

      stub[3] = diskstubstr(0);  // D1
      stub[4] = diskstubstr(1);  // D2
      stub[5] = diskstubstr(2);  // D3
      stub[6] = diskstubstr(3);  // D4

      break;

    case 2:                       // L3L4
      stub[0] = layerstubstr(0);  // L1
      stub[1] = layerstubstr(1);  // L2
      stub[2] = layerstubstr(4);  // L5
      stub[3] = layerstubstr(5);  // L6

      stub[4] = diskstubstr(0);  // D1
      stub[5] = diskstubstr(1);  // D2

      break;

    case 3:                       // L5L6
      stub[0] = layerstubstr(0);  // L1
      stub[1] = layerstubstr(1);  // L2
      stub[2] = layerstubstr(2);  // L3
      stub[3] = layerstubstr(3);  // L4

      break;

    case 4:                       // D1D2
      stub[0] = layerstubstr(0);  // L1
      stub[1] = layerstubstr(1);  // L2

      stub[2] = diskstubstr(2);  // D3
      stub[3] = diskstubstr(3);  // D4
      stub[4] = diskstubstr(4);  // D5

      break;

    case 5:                       // D3D4
      stub[0] = layerstubstr(0);  // L1

      stub[1] = diskstubstr(0);  // D1
      stub[2] = diskstubstr(1);  // D2
      stub[3] = diskstubstr(4);  // D5

      break;

    case 6:                      // L1D1
      stub[0] = diskstubstr(1);  // D2
      stub[1] = diskstubstr(2);  // D3
      stub[2] = diskstubstr(3);  // D4
      stub[3] = diskstubstr(4);  // D5

      break;

    case 7:                       // L2D1
      stub[0] = layerstubstr(0);  // L1

      stub[1] = diskstubstr(1);  // D2
      stub[2] = diskstubstr(2);  // D3
      stub[3] = diskstubstr(3);  // D4

      break;
  }

  // Only one hit per layer/disk is allowed currently, so the hit map for a
  // given layer/disk is just equal to the valid bit of the corresponding stub
  // string, which is the first character.
  for (unsigned i = 0; i < maxNHits; i++)
    hitmap[i * nBitsPerHit + 2] = stub[i][0];

  std::string oss("");
  //Binary print out
  if (!settings_.writeoutReal()) {
    const FPGAWord tmp(getISeed(), settings_.nbitsseed(), true, __LINE__, __FILE__);

    oss += "1|";  // valid bit
    oss += tmp.str() + "|";
    oss += fpgapars_.rinv().str() + "|";
    oss += fpgapars_.phi0().str() + "|";
    oss += fpgapars_.z0().str() + "|";
    oss += fpgapars_.t().str() + "|";
    oss += hitmap;
    for (unsigned i = 0; i < maxNHits; i++)
      // If a valid stub string was never assigned, then that stub is not
      // included in the output.
      if (stub[i] != "0")
        oss += "|" + stub[i];
  }

  return oss;
}

int Tracklet::trackIndex

(

)

const

Definition at line 815 of file Tracklet.cc.

References trackIndex_.

{ return trackIndex_; }

int trklet::Tracklet::trackletIndex ( ) const

inline

Definition at line 206 of file Tracklet.h.

References trackletIndex_.

Referenced by trklet::FullMatchMemory::addMatch(), and trklet::MatchEngineUnit::step().

{ return trackletIndex_; }

std::string Tracklet::trackletparstr

(

)

Definition at line 188 of file Tracklet.cc.

References fpgapars_, innerFPGAStub_, trklet::Settings::kd0pars(), trklet::Settings::kphi0pars(), trklet::Settings::krinvpars(), trklet::Settings::ktpars(), trklet::Settings::kz(), middleFPGAStub_, outerFPGAStub_, settings_, trklet::FPGAWord::str(), str, AlCaHLTBitMon_QueryRunRegistry::string, trklet::Stub::stubindex(), and trklet::Settings::writeoutReal().

                                   {
  if (settings_.writeoutReal()) {
    std::string oss = std::to_string(fpgapars_.rinv().value() * settings_.krinvpars()) + " " +
                      std::to_string(fpgapars_.phi0().value() * settings_.kphi0pars()) + " " +
                      std::to_string(fpgapars_.d0().value() * settings_.kd0pars()) + " " +
                      std::to_string(fpgapars_.z0().value() * settings_.kz()) + " " +
                      std::to_string(fpgapars_.t().value() * settings_.ktpars());
    return oss;
  } else {
    std::string str = innerFPGAStub_->stubindex().str() + "|";
    if (middleFPGAStub_) {
      str += middleFPGAStub_->stubindex().str() + "|";
    }
    str += outerFPGAStub_->stubindex().str() + "|" + fpgapars_.rinv().str() + "|" + fpgapars_.phi0().str() + "|";
    if (middleFPGAStub_)
      str += fpgapars_.d0().str() + "|";
    str += fpgapars_.z0().str() + "|" + fpgapars_.t().str();
    return str;
  }
}

std::string Tracklet::trackletprojstr

(

int

layer

)

const

Definition at line 259 of file Tracklet.cc.

References cms::cuda::assert(), Exception, trklet::Projection::fpgaphiproj(), trklet::Projection::fpgaphiprojder(), trklet::Projection::fpgarzproj(), trklet::Projection::fpgarzprojder(), trklet::N_LAYER, trklet::Settings::nbitstcindex(), trklet::Settings::nbitstrackletindex(), trklet::Settings::ntrackletmax(), proj_, trklet::FPGAWord::set(), settings_, trklet::FPGAWord::str(), AlCaHLTBitMon_QueryRunRegistry::string, TCIndex_, createJobs::tmp, and trackletIndex_.

Referenced by trackletprojstrlayer().

                                                   {
  assert(layer > 0 && layer <= N_LAYER);
  FPGAWord tmp;
  if (trackletIndex_ < 0 || trackletIndex_ > (int)settings_.ntrackletmax()) {
    throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__ << " trackletIndex_ = " << trackletIndex_;
  }
  tmp.set(trackletIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
  FPGAWord tcid;
  tcid.set(TCIndex_, settings_.nbitstcindex(), true, __LINE__, __FILE__);

  std::string oss = tcid.str() + "|" + tmp.str() + "|" + proj_[layer - 1].fpgaphiproj().str() + "|" +
                    proj_[layer - 1].fpgarzproj().str() + "|" + proj_[layer - 1].fpgaphiprojder().str() + "|" +
                    proj_[layer - 1].fpgarzprojder().str();
  return oss;
}

std::string Tracklet::trackletprojstrD

(

int

disk

)

const

Definition at line 275 of file Tracklet.cc.

References funct::abs(), cms::cuda::assert(), Exception, trklet::Settings::extended(), trklet::Projection::fpgaphiproj(), trklet::Projection::fpgaphiprojder(), trklet::Projection::fpgarzproj(), trklet::Projection::fpgarzprojder(), trklet::N_DISK, trklet::N_LAYER, trklet::Settings::nbitstrackletindex(), trklet::Settings::ntrackletmax(), proj_, trklet::FPGAWord::set(), settings_, trklet::FPGAWord::str(), AlCaHLTBitMon_QueryRunRegistry::string, TCIndex_, createJobs::tmp, and trackletIndex_.

Referenced by trackletprojstrdisk().

                                                   {
  assert(abs(disk) <= N_DISK);
  FPGAWord tmp;
  if (trackletIndex_ < 0 || trackletIndex_ > (int)settings_.ntrackletmax()) {
    throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__ << " trackletIndex_ = " << trackletIndex_;
  }
  tmp.set(trackletIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
  FPGAWord tcid;
  if (settings_.extended()) {
    tcid.set(TCIndex_, 8, true, __LINE__, __FILE__);
  } else {
    tcid.set(TCIndex_, 7, true, __LINE__, __FILE__);
  }
  std::string oss = tcid.str() + "|" + tmp.str() + "|" + proj_[N_LAYER + abs(disk) - 1].fpgaphiproj().str() + "|" +
                    proj_[N_LAYER + abs(disk) - 1].fpgarzproj().str() + "|" +
                    proj_[N_LAYER + abs(disk) - 1].fpgaphiprojder().str() + "|" +
                    proj_[N_LAYER + abs(disk) - 1].fpgarzprojder().str();
  return oss;
}

std::string trklet::Tracklet::trackletprojstrdisk ( int  disk ) const

inline

Definition at line 80 of file Tracklet.h.

References trackletprojstrD().

{ return trackletprojstrD(disk); }

std::string trklet::Tracklet::trackletprojstrlayer ( int  layer ) const

inline

Definition at line 79 of file Tracklet.h.

References trackletprojstr().

{ return trackletprojstr(layer); }

bool trklet::Tracklet::validProj ( int  layerdisk ) const

inline

Definition at line 82 of file Tracklet.h.

References cms::cuda::assert(), trklet::N_DISK, trklet::N_LAYER, proj_, and trklet::Projection::valid().

Referenced by trklet::TrackletProjectionsMemory::addProj(), trklet::TrackletCalculatorBase::barrelSeeding(), trklet::TrackletCalculatorDisplaced::DDLSeeding(), trklet::TrackletCalculatorBase::diskSeeding(), trklet::TrackletCalculatorDisplaced::LLDSeeding(), trklet::TrackletCalculatorDisplaced::LLLSeeding(), trklet::TrackletCalculatorBase::overlapSeeding(), and proj().

                                        {
      assert(layerdisk >= 0 && layerdisk < N_LAYER + N_DISK);
      return proj_[layerdisk].valid();
    }

std::string Tracklet::vmstrdisk	(	int	disk,
		unsigned int	allstubindex
	)

Definition at line 244 of file Tracklet.cc.

References trklet::Projection::fpgafinephivm(), trklet::Projection::fpgafinerzvm(), trklet::Projection::fpgarzbin1projvm(), trklet::Projection::fpgarzbin2projvm(), trklet::Projection::getBendIndex(), trklet::N_LAYER, proj_, trklet::FPGAWord::set(), trklet::FPGAWord::str(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                 {
  FPGAWord index;
  if (allstubindex >= (1 << 7)) {
    edm::LogPrint("Tracklet") << "Warning projection number too large!";
    index.set((1 << 7) - 1, 7, true, __LINE__, __FILE__);
  } else {
    index.set(allstubindex, 7, true, __LINE__, __FILE__);
  }
  std::string oss =
      index.str() + "|" + proj_[N_LAYER + disk - 1].fpgarzbin1projvm().str() + "|" +
      proj_[N_LAYER + disk - 1].fpgarzbin2projvm().str() + "|" + proj_[N_LAYER + disk - 1].fpgafinerzvm().str() + "|" +
      proj_[N_LAYER + disk - 1].fpgafinephivm().str() + "|" + proj_[N_LAYER + disk - 1].getBendIndex().str();
  return oss;
}

std::string Tracklet::vmstrlayer	(	int	layer,
		unsigned int	allstubindex
	)

Definition at line 209 of file Tracklet.cc.

References cms::cuda::assert(), trklet::Settings::extended(), trklet::Projection::fpgafinephivm(), trklet::Projection::fpgafinerzvm(), trklet::Projection::fpgaphiprojder(), trklet::Projection::fpgarzbin1projvm(), trklet::Projection::fpgarzbin2projvm(), trklet::FPGAWord::nbits(), proj_, PSseed(), trklet::FPGAWord::set(), settings_, trklet::FPGAWord::str(), AlCaHLTBitMon_QueryRunRegistry::string, createJobs::tmp, and trklet::FPGAWord::value().

                                                                   {
  FPGAWord index;
  if (allstubindex >= (1 << 7)) {
    edm::LogPrint("Tracklet") << "Warning projection number too large!";
    index.set((1 << 7) - 1, 7, true, __LINE__, __FILE__);
  } else {
    index.set(allstubindex, 7, true, __LINE__, __FILE__);
  }

  // This is a shortcut.
  //int irinvvm=16+(fpgarinv().value()>>(fpgarinv().nbits()-5));
  // rinv is not directly available in the TrackletProjection.
  // can be inferred from phi derivative: rinv = - phider * 2
  int tmp_irinv = proj_[layer - 1].fpgaphiprojder().value() * (-2);
  int nbits_irinv = proj_[layer - 1].fpgaphiprojder().nbits() + 1;

  // irinv in VMProjection:
  // top 5 bits of rinv and shifted to be positive
  int irinvvm = 16 + (tmp_irinv >> (nbits_irinv - 5));

  if (settings_.extended() && (irinvvm > 31)) {  //TODO - displaced tracking should protect against this
    edm::LogPrint("Tracklet") << "Warning irinvvm too large:" << irinvvm;
    irinvvm = 31;
  }

  assert(irinvvm >= 0);
  assert(irinvvm < 32);
  FPGAWord tmp;
  tmp.set(irinvvm, 5, true, __LINE__, __FILE__);
  std::string oss = index.str() + "|" + proj_[layer - 1].fpgarzbin1projvm().str() + "|" +
                    proj_[layer - 1].fpgarzbin2projvm().str() + "|" + proj_[layer - 1].fpgafinerzvm().str() + "|" +
                    proj_[layer - 1].fpgafinephivm().str() + "|" + tmp.str() + "|" + std::to_string(PSseed());
  return oss;
}

double trklet::Tracklet::z0 ( ) const

inline

Definition at line 124 of file Tracklet.h.

References trackpars_, and trklet::TrackPars< T >::z0().

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return trackpars_.z0(); }

double trklet::Tracklet::z0approx ( ) const

inline

Definition at line 130 of file Tracklet.h.

References trackparsapprox_, and trklet::TrackPars< T >::z0().

Referenced by trklet::FitTrack::trackFitChisq(), and trklet::FitTrack::trackFitFake().

{ return trackparsapprox_.z0(); }

double trklet::Tracklet::z0fit ( ) const

inline

Definition at line 142 of file Tracklet.h.

References fitpars_, and trklet::TrackPars< T >::z0().

{ return fitpars_.z0(); }

double trklet::Tracklet::z0fitexact ( ) const

inline

Definition at line 149 of file Tracklet.h.

References fitparsexact_, and trklet::TrackPars< T >::z0().

{ return fitparsexact_.z0(); }

Member Data Documentation

bool trklet::Tracklet::barrel_

private

Definition at line 228 of file Tracklet.h.

Referenced by getStubIDs(), isBarrel(), and Tracklet().

double trklet::Tracklet::chisqrphifit_

private

Definition at line 257 of file Tracklet.h.

Referenced by chiSqfit(), makeTrack(), and setFitPars().

double trklet::Tracklet::chisqrphifitexact_

private

Definition at line 261 of file Tracklet.h.

Referenced by setFitPars().

double trklet::Tracklet::chisqrzfit_

private

Definition at line 258 of file Tracklet.h.

Referenced by chiSqfit(), makeTrack(), and setFitPars().

double trklet::Tracklet::chisqrzfitexact_

private

Definition at line 262 of file Tracklet.h.

Referenced by setFitPars().

bool trklet::Tracklet::disk_

private

Definition at line 229 of file Tracklet.h.

Referenced by getStubIDs(), isDisk(), and Tracklet().

TrackPars<double> trklet::Tracklet::fitpars_

private

Definition at line 256 of file Tracklet.h.

Referenced by d0fit(), phi0fit(), rinvfit(), setFitPars(), tfit(), and z0fit().

TrackPars<double> trklet::Tracklet::fitparsexact_

private

Definition at line 260 of file Tracklet.h.

Referenced by d0fitexact(), phi0fitexact(), rinvfitexact(), setFitPars(), tfitexact(), and z0fitexact().

TrackPars<FPGAWord> trklet::Tracklet::fpgafitpars_

private

Definition at line 252 of file Tracklet.h.

Referenced by id0fit(), iphi0fit(), irinvfit(), itfit(), iz0fit(), makeTrack(), and setFitPars().

TrackPars<FPGAWord> trklet::Tracklet::fpgapars_

private

Definition at line 242 of file Tracklet.h.

Referenced by fpgad0(), fpgaphi0(), fpgarinv(), fpgat(), fpgaz0(), trackfitstr(), Tracklet(), and trackletparstr().

std::unique_ptr<Track> trklet::Tracklet::fpgatrack_

private

Definition at line 266 of file Tracklet.h.

Referenced by getTrack(), setFitPars(), and Tracklet().

int trklet::Tracklet::hitpattern_

private

Definition at line 264 of file Tracklet.h.

Referenced by makeTrack(), and setFitPars().

FPGAWord trklet::Tracklet::ichisqrphifit_

private

Definition at line 253 of file Tracklet.h.

Referenced by fit(), ichiSqfit(), makeTrack(), setFitPars(), and Tracklet().

FPGAWord trklet::Tracklet::ichisqrzfit_

private

Definition at line 254 of file Tracklet.h.

Referenced by ichiSqfit(), makeTrack(), setFitPars(), and Tracklet().

const Stub* trklet::Tracklet::innerFPGAStub_

private

Definition at line 233 of file Tracklet.h.

Referenced by addressstr(), disk(), getL1Stubs(), getStubIDs(), innerFPGAStub(), layer(), stubtruthmatch(), tpseed(), Tracklet(), and trackletparstr().

const Stub* trklet::Tracklet::middleFPGAStub_

private

Definition at line 234 of file Tracklet.h.

Referenced by addressstr(), disk(), getL1Stubs(), getStubIDs(), layer(), middleFPGAStub(), Tracklet(), and trackletparstr().

const Stub* trklet::Tracklet::outerFPGAStub_

private

Definition at line 235 of file Tracklet.h.

Referenced by addressstr(), disk(), getL1Stubs(), getStubIDs(), layer(), outerFPGAStub(), stubtruthmatch(), tpseed(), Tracklet(), and trackletparstr().

bool trklet::Tracklet::overlap_

private

Definition at line 230 of file Tracklet.h.

Referenced by getStubIDs(), isOverlap(), and Tracklet().

Projection trklet::Tracklet::proj_[N_LAYER+N_DISK]

private

Definition at line 268 of file Tracklet.h.

Referenced by proj(), Tracklet(), trackletprojstr(), trackletprojstrD(), validProj(), vmstrdisk(), and vmstrlayer().

int trklet::Tracklet::projdisk_[N_DISK]

private

Definition at line 249 of file Tracklet.h.

Referenced by Tracklet().

int trklet::Tracklet::projlayer_[N_LAYER-2]

private

Definition at line 248 of file Tracklet.h.

Referenced by Tracklet().

Residual trklet::Tracklet::resid_[N_LAYER+N_DISK]

private

Definition at line 270 of file Tracklet.h.

Referenced by addMatch(), diskstubstr(), fullmatchdiskstr(), fullmatchstr(), getL1Stubs(), getStubIDs(), layerstubstr(), match(), and resid().

unsigned int trklet::Tracklet::seedIndex_

private

Definition at line 225 of file Tracklet.h.

Referenced by seedIndex(), and Tracklet().

Settings const& trklet::Tracklet::settings_

private

Definition at line 272 of file Tracklet.h.

Referenced by addMatch(), diskstubstr(), fullmatchdiskstr(), fullmatchstr(), getISeed(), getITC(), layerstubstr(), setFitPars(), setTrackIndex(), setTrackletIndex(), TCID(), trackfitstr(), Tracklet(), trackletparstr(), trackletprojstr(), trackletprojstrD(), and vmstrlayer().

int trklet::Tracklet::TCIndex_

private

Definition at line 238 of file Tracklet.h.

Referenced by fullmatchdiskstr(), fullmatchstr(), getISeed(), getITC(), setTCIndex(), TCID(), TCIndex(), Tracklet(), trackletprojstr(), and trackletprojstrD().

int trklet::Tracklet::trackIndex_

private

Definition at line 239 of file Tracklet.h.

Referenced by diskstubstr(), layerstubstr(), setTrackIndex(), and trackIndex().

int trklet::Tracklet::trackletIndex_

private

Definition at line 237 of file Tracklet.h.

Referenced by fullmatchdiskstr(), fullmatchstr(), setTrackletIndex(), TCID(), Tracklet(), trackletIndex(), trackletprojstr(), and trackletprojstrD().

TrackPars<double> trklet::Tracklet::trackpars_

private

Definition at line 244 of file Tracklet.h.

Referenced by d0(), phi0(), rinv(), t(), Tracklet(), and z0().

TrackPars<double> trklet::Tracklet::trackparsapprox_

private

Definition at line 245 of file Tracklet.h.

Referenced by d0approx(), phi0approx(), rinvapprox(), tapprox(), Tracklet(), and z0approx().

bool trklet::Tracklet::triplet_

private

Definition at line 231 of file Tracklet.h.

Referenced by Tracklet().