CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
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 FPGAWordfpgad0 () const
 
const FPGAWordfpgaphi0 () const
 
const FPGAWordfpgarinv () const
 
const FPGAWordfpgat () const
 
const FPGAWordfpgaz0 () 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 ()
 
TrackgetTrack ()
 
FPGAWord ichiSqfit () const
 
const FPGAWordid0fit () const
 
const StubinnerFPGAStub ()
 
const FPGAWordiphi0fit () const
 
const FPGAWordirinvfit () const
 
bool isBarrel () const
 
int isDisk () const
 
bool isOverlap () const
 
const FPGAWorditfit () const
 
const FPGAWordiz0fit () 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 StubmiddleFPGAStub ()
 
const StubouterFPGAStub ()
 
double phi0 () const
 
double phi0approx () const
 
double phi0fit () const
 
double phi0fitexact () const
 
Projectionproj (int layerdisk)
 
unsigned int PSseed () const
 
const Residualresid (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< FPGAWordfpgafitpars_
 
TrackPars< FPGAWordfpgapars_
 
std::unique_ptr< Trackfpgatrack_
 
int hitpattern_
 
FPGAWord ichisqrphifit_
 
FPGAWord ichisqrzfit_
 
const StubinnerFPGAStub_
 
const StubmiddleFPGAStub_
 
const StubouterFPGAStub_
 
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::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_, d0(), d0approx(), 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_, hltL1SingleMuFiltered5_cfi::overlap, overlap_, phi0(), phi0approx(), proj_, projdisk_, trklet::Settings::projdisks(), projlayer_, trklet::Settings::projlayers(), rinv(), rinvapprox(), seedIndex_, trklet::FPGAWord::set(), settings_, t(), tapprox(), TCIndex_, trackletIndex_, trackpars_, trackparsapprox_, triplet_, validateGeometry_cfg::valid, z0(), and z0approx().

39  : settings_(settings) {
40  seedIndex_ = iSeed;
41 
42  overlap_ = overlap;
43  disk_ = disk;
44  assert(!(disk && overlap));
45  barrel_ = (!disk) && (!overlap);
46  triplet_ = false;
47 
48  trackletIndex_ = -1;
49  TCIndex_ = -1;
50 
52 
53  if (barrel_ && middleFPGAStub == nullptr)
55 
59 
60  trackpars_.init(rinv, phi0, d0, t, z0);
61 
63 
64  fpgapars_.rinv().set(irinv, settings_.nbitsrinv(), false, __LINE__, __FILE__);
65  fpgapars_.phi0().set(iphi0, settings_.nbitsphi0(), false, __LINE__, __FILE__);
66  fpgapars_.d0().set(id0, settings_.nbitsd0(), false, __LINE__, __FILE__);
67  fpgapars_.z0().set(iz0, settings_.nbitsz0(), false, __LINE__, __FILE__);
68  fpgapars_.t().set(it, settings_.nbitst(), false, __LINE__, __FILE__);
69 
70  fpgatrack_ = nullptr;
71 
72  triplet_ = (seedIndex_ >= 8);
73 
74  //fill projection layers
75  for (unsigned int i = 0; i < N_LAYER - 2; i++) {
76  projlayer_[i] = settings.projlayers(seedIndex_, i);
77  }
78 
79  //fill projection disks
80  for (unsigned int i = 0; i < N_DISK; i++) {
81  projdisk_[i] = settings.projdisks(seedIndex_, i);
82  }
83 
84  //Handle projections to the layers
85  for (unsigned int i = 0; i < N_LAYER - 2; i++) {
86  if (projlayer_[i] == 0)
87  continue;
88  if (!projs[projlayer_[i] - 1].valid())
89  continue;
90 
91  proj_[projlayer_[i] - 1] = projs[projlayer_[i] - 1];
92  }
93  //Now handle projections to the disks
94  for (unsigned int i = 0; i < N_DISK; i++) {
95  if (projdisk_[i] == 0)
96  continue;
97  if (!projs[N_LAYER + projdisk_[i] - 1].valid())
98  continue;
99 
100  proj_[N_LAYER + projdisk_[i] - 1] = projs[N_LAYER + projdisk_[i] - 1];
101  }
102 
103  ichisqrphifit_.set(-1, 8, false);
104  ichisqrzfit_.set(-1, 8, false);
105 }
TrackPars< FPGAWord > fpgapars_
Definition: Tracklet.h:244
double t() const
Definition: Tracklet.h:123
const Stub * innerFPGAStub()
Definition: Tracklet.h:61
Projection proj_[N_LAYER+N_DISK]
Definition: Tracklet.h:270
unsigned int layer() const
Definition: L1TStub.h:45
constexpr int N_DISK
Definition: Settings.h:22
const Stub * outerFPGAStub()
Definition: Tracklet.h:65
int disk() const
Definition: Tracklet.cc:782
double phi0() const
Definition: Tracklet.h:121
TrackPars< double > trackparsapprox_
Definition: Tracklet.h:247
assert(be >=bs)
unsigned int nbitsd0() const
Definition: Settings.h:321
Settings const & settings_
Definition: Tracklet.h:274
FPGAWord ichisqrzfit_
Definition: Tracklet.h:256
FPGAWord ichisqrphifit_
Definition: Tracklet.h:255
int nbitst() const
Definition: Settings.h:348
const Stub * middleFPGAStub()
Definition: Tracklet.h:63
double phi0approx() const
Definition: Tracklet.h:127
double rinv() const
Definition: Tracklet.h:120
int projdisk_[N_DISK]
Definition: Tracklet.h:251
int projlayer_[N_LAYER - 2]
Definition: Tracklet.h:250
int nbitsrinv() const
Definition: Settings.h:346
unsigned int seedIndex_
Definition: Tracklet.h:227
std::unique_ptr< Track > fpgatrack_
Definition: Tracklet.h:268
double tapprox() const
Definition: Tracklet.h:129
L1TStub * l1tstub()
Definition: Stub.h:77
double rinvapprox() const
Definition: Tracklet.h:126
double d0() const
Definition: Tracklet.h:122
double z0approx() const
Definition: Tracklet.h:130
void set(int value, int nbits, bool positive=true, int line=-1, const char *file=nullptr)
Definition: FPGAWord.cc:14
double d0approx() const
Definition: Tracklet.h:128
double z0() const
Definition: Tracklet.h:124
const Stub * outerFPGAStub_
Definition: Tracklet.h:237
const Stub * innerFPGAStub_
Definition: Tracklet.h:235
const Stub * middleFPGAStub_
Definition: Tracklet.h:236
int nbitsphi0() const
Definition: Settings.h:347
TrackPars< double > trackpars_
Definition: Tracklet.h:246
int nbitsz0() const
Definition: Settings.h:349
void init(T rinv, T phi0, T d0, T t, T z0)
Definition: TrackPars.h:21
constexpr int N_LAYER
Definition: Settings.h:21

◆ ~Tracklet()

trklet::Tracklet::~Tracklet ( )
default

Member Function Documentation

◆ addMatch()

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().

303  {
304  assert(layerdisk < N_LAYER + N_DISK);
305  resid_[layerdisk].init(settings_, layerdisk, ideltaphi, ideltarz, stubid, dphi, drz, dphiapprox, drzapprox, stubptr);
306 }
constexpr int N_DISK
Definition: Settings.h:22
void init(Settings const &settings, unsigned int layerdisk, int iphiresid, int irzresid, int istubid, double phiresid, double rzresid, double phiresidapprox, double rzresidapprox, const Stub *stubptr)
Definition: Residual.cc:7
assert(be >=bs)
Settings const & settings_
Definition: Tracklet.h:274
Residual resid_[N_LAYER+N_DISK]
Definition: Tracklet.h:272
constexpr int N_LAYER
Definition: Settings.h:21

◆ addressstr()

std::string Tracklet::addressstr ( )

Definition at line 177 of file Tracklet.cc.

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

177  {
180  if (middleFPGAStub_) {
182  }
184 
185  return str;
186 }
std::string phiregionaddressstr() const
Definition: Stub.cc:117
const Stub * outerFPGAStub_
Definition: Tracklet.h:237
const Stub * innerFPGAStub_
Definition: Tracklet.h:235
const Stub * middleFPGAStub_
Definition: Tracklet.h:236
#define str(s)

◆ chiSqfit()

double trklet::Tracklet::chiSqfit ( ) const
inline

Definition at line 143 of file Tracklet.h.

References chisqrphifit_, and chisqrzfit_.

143 { return chisqrphifit_ + chisqrzfit_; }
double chisqrzfit_
Definition: Tracklet.h:260
double chisqrphifit_
Definition: Tracklet.h:259

◆ d0()

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(), and Tracklet().

122 { return trackpars_.d0(); }
const T & d0() const
Definition: TrackPars.h:31
TrackPars< double > trackpars_
Definition: Tracklet.h:246

◆ d0approx()

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(), and Tracklet().

128 { return trackparsapprox_.d0(); }
TrackPars< double > trackparsapprox_
Definition: Tracklet.h:247
const T & d0() const
Definition: TrackPars.h:31

◆ d0fit()

double trklet::Tracklet::d0fit ( ) const
inline

Definition at line 140 of file Tracklet.h.

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

Referenced by setFitPars().

140 { return fitpars_.d0(); }
TrackPars< double > fitpars_
Definition: Tracklet.h:258
const T & d0() const
Definition: TrackPars.h:31

◆ d0fitexact()

double trklet::Tracklet::d0fitexact ( ) const
inline

Definition at line 147 of file Tracklet.h.

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

Referenced by setFitPars().

147 { return fitparsexact_.d0(); }
const T & d0() const
Definition: TrackPars.h:31
TrackPars< double > fitparsexact_
Definition: Tracklet.h:262

◆ disk()

int Tracklet::disk ( ) const

Definition at line 782 of file Tracklet.cc.

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

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

782  {
783  int d1 = (innerFPGAStub_ && (innerFPGAStub_->layerdisk() >= N_LAYER)) ? innerFPGAStub_->l1tstub()->disk() : 999,
785  d3 = (outerFPGAStub_ && (outerFPGAStub_->layerdisk() >= N_LAYER)) ? outerFPGAStub_->l1tstub()->disk() : 999,
786  d = 999;
787  if (abs(d1) < min(abs(d2), abs(d3)))
788  d = d1;
789  if (abs(d2) < min(abs(d1), abs(d3)))
790  d = d2;
791  if (abs(d3) < min(abs(d1), abs(d2)))
792  d = d3;
793  return (d < 999 ? d : 0);
794 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
L1TStub * l1tstub()
Definition: Stub.h:77
int disk() const
Definition: L1TStub.h:46
unsigned int layerdisk() const
Definition: Stub.cc:185
d
Definition: ztail.py:151
const Stub * outerFPGAStub_
Definition: Tracklet.h:237
const Stub * innerFPGAStub_
Definition: Tracklet.h:235
const Stub * middleFPGAStub_
Definition: Tracklet.h:236
static constexpr float d1
constexpr int N_LAYER
Definition: Settings.h:21

◆ diskstubstr()

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, trackIndex_, and validateGeometry_cfg::valid.

Referenced by trackfitstr().

595  {
596  assert(disk < N_DISK);
597 
598  std::stringstream oss("");
599  if (!resid_[N_LAYER + disk].valid())
600  oss << "0|0000000|0000000000|000000000000|000000000000|0000000";
601  else {
602  if (trackIndex_ < 0 || trackIndex_ > (int)settings_.ntrackletmax()) {
603  cout << "trackIndex_ = " << trackIndex_ << endl;
604  assert(0);
605  }
606  const FPGAWord tmp(trackIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
607  const FPGAWord& stubr = resid_[N_LAYER + disk].stubptr()->r();
608  const bool isPS = resid_[N_LAYER + disk].stubptr()->isPSmodule();
609  oss << "1|"; // valid bit
610  oss << tmp.str() << "|";
611  oss << resid_[N_LAYER + disk].fpgastubid().str() << "|";
612  oss << (isPS ? stubr.str() : ("00000000" + stubr.str())) << "|";
613  oss << resid_[N_LAYER + disk].fpgaphiresid().str() << "|";
614  oss << resid_[N_LAYER + disk].fpgarzresid().str();
615  }
616 
617  return oss.str();
618 }
constexpr int N_DISK
Definition: Settings.h:22
const FPGAWord & r() const
Definition: Stub.h:60
int disk() const
Definition: Tracklet.cc:782
unsigned int nbitstrackletindex() const
Definition: Settings.h:281
const FPGAWord & fpgastubid() const
Definition: Residual.h:42
const FPGAWord & fpgaphiresid() const
Definition: Residual.h:32
assert(be >=bs)
Settings const & settings_
Definition: Tracklet.h:274
const FPGAWord & fpgarzresid() const
Definition: Residual.h:37
Residual resid_[N_LAYER+N_DISK]
Definition: Tracklet.h:272
bool isPSmodule() const
Definition: Stub.h:71
const Stub * stubptr() const
Definition: Residual.h:67
unsigned int ntrackletmax() const
Definition: Settings.h:343
tmp
align.sh
Definition: createJobs.py:716
std::string str() const
Definition: FPGAWord.cc:54
constexpr int N_LAYER
Definition: Settings.h:21

◆ fit()

bool trklet::Tracklet::fit ( ) const
inline

◆ fpgad0()

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

Definition at line 134 of file Tracklet.h.

References fpgapars_.

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

134 { return fpgapars_.d0(); }
TrackPars< FPGAWord > fpgapars_
Definition: Tracklet.h:244

◆ fpgaphi0()

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

Definition at line 133 of file Tracklet.h.

References fpgapars_.

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

133 { return fpgapars_.phi0(); }
TrackPars< FPGAWord > fpgapars_
Definition: Tracklet.h:244

◆ fpgarinv()

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

Definition at line 132 of file Tracklet.h.

References fpgapars_.

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

132 { return fpgapars_.rinv(); }
TrackPars< FPGAWord > fpgapars_
Definition: Tracklet.h:244

◆ fpgat()

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

Definition at line 135 of file Tracklet.h.

References fpgapars_.

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

135 { return fpgapars_.t(); }
TrackPars< FPGAWord > fpgapars_
Definition: Tracklet.h:244

◆ fpgaz0()

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

Definition at line 136 of file Tracklet.h.

References fpgapars_.

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

136 { return fpgapars_.z0(); }
TrackPars< FPGAWord > fpgapars_
Definition: Tracklet.h:244

◆ fullmatchdiskstr()

std::string Tracklet::fullmatchdiskstr ( int  disk)

Definition at line 324 of file Tracklet.cc.

References cms::cuda::assert(), disk(), 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_.

324  {
325  assert(disk > 0 && disk <= N_DISK);
326 
327  FPGAWord tmp;
328  if (trackletIndex_ < 0 || trackletIndex_ > (int)settings_.ntrackletmax()) {
329  throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__ << " trackletIndex_ = " << trackletIndex_;
330  }
331  tmp.set(trackletIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
332  FPGAWord tcid;
333  tcid.set(TCIndex_, settings_.nbitstcindex(), true, __LINE__, __FILE__);
334  const FPGAWord& stubr = resid_[N_LAYER + disk - 1].stubptr()->r();
335  const bool isPS = resid_[N_LAYER + disk - 1].stubptr()->isPSmodule();
336  std::string oss = tcid.str() + "|" + tmp.str() + "|" + resid_[N_LAYER + disk - 1].fpgastubid().str() + "|" +
337  (isPS ? stubr.str() : ("00000000" + stubr.str())) + "|" +
338  resid_[N_LAYER + disk - 1].fpgaphiresid().str() + "|" +
339  resid_[N_LAYER + disk - 1].fpgarzresid().str();
340  return oss;
341 }
constexpr int N_DISK
Definition: Settings.h:22
const FPGAWord & r() const
Definition: Stub.h:60
int disk() const
Definition: Tracklet.cc:782
unsigned int nbitstrackletindex() const
Definition: Settings.h:281
const FPGAWord & fpgastubid() const
Definition: Residual.h:42
const FPGAWord & fpgaphiresid() const
Definition: Residual.h:32
assert(be >=bs)
Settings const & settings_
Definition: Tracklet.h:274
const FPGAWord & fpgarzresid() const
Definition: Residual.h:37
unsigned int nbitstcindex() const
Definition: Settings.h:286
Residual resid_[N_LAYER+N_DISK]
Definition: Tracklet.h:272
void set(int value, int nbits, bool positive=true, int line=-1, const char *file=nullptr)
Definition: FPGAWord.cc:14
bool isPSmodule() const
Definition: Stub.h:71
const Stub * stubptr() const
Definition: Residual.h:67
unsigned int ntrackletmax() const
Definition: Settings.h:343
tmp
align.sh
Definition: createJobs.py:716
std::string str() const
Definition: FPGAWord.cc:54
constexpr int N_LAYER
Definition: Settings.h:21

◆ fullmatchstr()

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(), layer(), 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_.

308  {
309  assert(layer > 0 && layer <= N_LAYER);
310 
311  FPGAWord tmp;
312  if (trackletIndex_ < 0 || trackletIndex_ > (int)settings_.ntrackletmax()) {
313  throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__ << " trackletIndex_ = " << trackletIndex_;
314  }
315  tmp.set(trackletIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
316  FPGAWord tcid;
317  tcid.set(TCIndex_, settings_.nbitstcindex(), true, __LINE__, __FILE__);
318  std::string oss = tcid.str() + "|" + tmp.str() + "|" + resid_[layer - 1].fpgastubid().str() + "|" +
319  resid_[layer - 1].stubptr()->r().str() + "|" + resid_[layer - 1].fpgaphiresid().str() + "|" +
320  resid_[layer - 1].fpgarzresid().str();
321  return oss;
322 }
const FPGAWord & r() const
Definition: Stub.h:60
unsigned int nbitstrackletindex() const
Definition: Settings.h:281
const FPGAWord & fpgastubid() const
Definition: Residual.h:42
const FPGAWord & fpgaphiresid() const
Definition: Residual.h:32
assert(be >=bs)
Settings const & settings_
Definition: Tracklet.h:274
const FPGAWord & fpgarzresid() const
Definition: Residual.h:37
unsigned int nbitstcindex() const
Definition: Settings.h:286
Residual resid_[N_LAYER+N_DISK]
Definition: Tracklet.h:272
void set(int value, int nbits, bool positive=true, int line=-1, const char *file=nullptr)
Definition: FPGAWord.cc:14
const Stub * stubptr() const
Definition: Residual.h:67
unsigned int ntrackletmax() const
Definition: Settings.h:343
int layer() const
Definition: Tracklet.cc:773
tmp
align.sh
Definition: createJobs.py:716
std::string str() const
Definition: FPGAWord.cc:54
constexpr int N_LAYER
Definition: Settings.h:21

◆ getISeed()

int Tracklet::getISeed ( ) const

◆ getITC()

int Tracklet::getITC ( ) const

Definition at line 807 of file Tracklet.cc.

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

807  {
808  const int iSeed = getISeed(), iTC = TCIndex_ - (iSeed << settings_.nbitsitc());
809  assert(iTC >= 0 && iTC <= 14);
810  return iTC;
811 }
assert(be >=bs)
Settings const & settings_
Definition: Tracklet.h:274
unsigned int nbitsitc() const
Definition: Settings.h:284
int getISeed() const
Definition: Tracklet.cc:801

◆ getL1Stubs()

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().

343  {
344  std::vector<const L1TStub*> tmp;
345 
346  if (innerFPGAStub_)
347  tmp.push_back(innerFPGAStub_->l1tstub());
348  if (middleFPGAStub_)
349  tmp.push_back(middleFPGAStub_->l1tstub());
350  if (outerFPGAStub_)
351  tmp.push_back(outerFPGAStub_->l1tstub());
352 
353  for (const auto& iresid : resid_) {
354  if (iresid.valid())
355  tmp.push_back(iresid.stubptr()->l1tstub());
356  }
357 
358  return tmp;
359 }
L1TStub * l1tstub()
Definition: Stub.h:77
Residual resid_[N_LAYER+N_DISK]
Definition: Tracklet.h:272
const Stub * outerFPGAStub_
Definition: Tracklet.h:237
const Stub * innerFPGAStub_
Definition: Tracklet.h:235
const Stub * middleFPGAStub_
Definition: Tracklet.h:236
tmp
align.sh
Definition: createJobs.py:716

◆ getStubIDs()

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(), EcalCondDBWriter_cfi::location, middleFPGAStub_, trklet::N_BITSMEMADDRESS, trklet::N_DISK, trklet::N_LAYER, trklet::FPGAWord::nbits(), outerFPGAStub_, overlap_, trklet::Stub::phiregionaddress(), resid_, trklet::Stub::stubindex(), validateGeometry_cfg::valid, trklet::Residual::valid(), trklet::FPGAWord::value(), and relativeConstraints::value.

Referenced by makeTrack().

361  {
362  std::map<int, int> stubIDs;
363 
364  // For future reference, *resid_[i] uses i as the absolute stub index. (0-5 for barrel, 0-4 for disk)
365  // On the other hand, proj*_[i] uses i almost like *resid_[i], except the *seeding* layer indices are removed entirely.
366  // E.g. An L3L4 track has 0=L1, 1=L2, 2=L4, 3=L5 for the barrels (for proj*_[i])
367 
368  if (innerFPGAStub_)
370  if (middleFPGAStub_)
372  if (outerFPGAStub_)
374 
375  if (barrel_) {
376  for (int i = 0; i < N_LAYER; i++) {
377  //check barrel
378  if (resid_[i].valid()) {
379  // two extra bits to indicate if the matched stub is local or from neighbor
380  int location = 1; // local
381  location <<= resid_[i].fpgastubid().nbits();
382 
383  stubIDs[1 + i] = resid_[i].fpgastubid().value() + location;
384  }
385 
386  //check disk
387  if (i >= N_DISK)
388  continue; //i=[0..4] for disks
389  if (resid_[N_LAYER + i].valid()) {
390  if (i == 3 && resid_[0].valid() && innerFPGAStub_->layer().value() == 1)
391  continue; // Don't add D4 if track has L1 stub
392  // two extra bits to indicate if the matched stub is local or from neighbor
393  int location = 1; // local
395 
396  if (itfit().value() < 0) {
397  stubIDs[-(N_LAYER + N_DISK) - i] = resid_[N_LAYER + i].fpgastubid().value() + location;
398  } else {
399  stubIDs[N_LAYER + N_DISK + i] = resid_[N_LAYER + i].fpgastubid().value() + location;
400  }
401  }
402  }
403 
404  //get stubs making up tracklet
405  if (innerFPGAStub_)
406  stubIDs[innerFPGAStub_->layer().value() + 1] = innerFPGAStub_->phiregionaddress() + (1 << 10);
407  if (middleFPGAStub_)
408  stubIDs[middleFPGAStub_->layer().value() + 1] = middleFPGAStub_->phiregionaddress() + (1 << 10);
409  if (outerFPGAStub_)
410  stubIDs[outerFPGAStub_->layer().value() + 1] = outerFPGAStub_->phiregionaddress() + (1 << 10);
411 
412  } else if (disk_) {
413  for (int i = 0; i < N_DISK; i++) {
414  //check barrel
415  if (resid_[i].valid()) {
416  // two extra bits to indicate if the matched stub is local or from neighbor
417  int location = 1; // local
418  location <<= resid_[i].fpgastubid().nbits();
419 
420  stubIDs[1 + i] = resid_[i].fpgastubid().value() + location;
421  }
422 
423  //check disks
424  if (i == 4 && resid_[1].valid())
425  continue; // Don't add D5 if track has L2 stub
426  if (resid_[N_LAYER + i].valid()) {
427  // two extra bits to indicate if the matched stub is local or from neighbor
428  int location = 1; // local
430 
431  if (innerFPGAStub_->l1tstub()->disk() < 0) {
432  stubIDs[-11 - i] = resid_[N_LAYER + i].fpgastubid().value() + location;
433  } else {
434  stubIDs[11 + i] = resid_[N_LAYER + i].fpgastubid().value() + location;
435  }
436  }
437  }
438 
439  //get stubs making up tracklet
440  if (innerFPGAStub_->disk().value() < 0) { //negative side runs 6-10
441  if (innerFPGAStub_)
442  stubIDs[innerFPGAStub_->disk().value() - 10] = innerFPGAStub_->phiregionaddress() + (1 << 10);
443  if (middleFPGAStub_)
444  stubIDs[middleFPGAStub_->disk().value() - 10] = middleFPGAStub_->phiregionaddress() + (1 << 10);
445  if (outerFPGAStub_)
446  stubIDs[outerFPGAStub_->disk().value() - 10] = outerFPGAStub_->phiregionaddress() + (1 << 10);
447  } else { // positive side runs 11-15]
448  if (innerFPGAStub_)
449  stubIDs[innerFPGAStub_->disk().value() + 10] = innerFPGAStub_->phiregionaddress() + (1 << 10);
450  if (middleFPGAStub_)
451  stubIDs[middleFPGAStub_->disk().value() + 10] = middleFPGAStub_->phiregionaddress() + (1 << 10);
452  if (outerFPGAStub_)
453  stubIDs[outerFPGAStub_->disk().value() + 10] = outerFPGAStub_->phiregionaddress() + (1 << 10);
454  }
455 
456  } else if (overlap_) {
457  for (int i = 0; i < N_DISK; i++) {
458  //check barrel
459  if (resid_[i].valid()) {
460  // two extra bits to indicate if the matched stub is local or from neighbor
461  int location = 1; // local
462  location <<= resid_[i].fpgastubid().nbits();
463 
464  stubIDs[1 + i] = resid_[i].fpgastubid().value() + location;
465  }
466 
467  //check disks
468  if (resid_[N_LAYER + i].valid()) {
469  // two extra bits to indicate if the matched stub is local or from neighbor
470  int location = 1; // local
472 
473  if (innerFPGAStub_->l1tstub()->disk() < 0) { // if negative overlap
474  if (innerFPGAStub_->layer().value() != 2 || !resid_[0].valid() ||
475  i != 3) { // Don't add D4 if this is an L3L2 track with an L1 stub
476  stubIDs[-11 - i] = resid_[N_LAYER + i].fpgastubid().value() + location;
477  }
478  } else {
479  if (innerFPGAStub_->layer().value() != 2 || !resid_[0].valid() || i != 3) {
480  stubIDs[11 + i] = resid_[N_LAYER + i].fpgastubid().value() + location;
481  }
482  }
483  }
484  }
485 
486  //get stubs making up tracklet
487 
488  if (innerFPGAStub_->layer().value() == 2) { // L3L2 track
489  if (innerFPGAStub_)
490  stubIDs[innerFPGAStub_->layer().value() + 1] = innerFPGAStub_->phiregionaddress() + (1 << 10);
491  if (middleFPGAStub_)
492  stubIDs[middleFPGAStub_->layer().value() + 1] = middleFPGAStub_->phiregionaddress() + (1 << 10);
493  if (outerFPGAStub_)
494  stubIDs[outerFPGAStub_->layer().value() + 1] = outerFPGAStub_->phiregionaddress() + (1 << 10);
495  } else if (innerFPGAStub_->disk().value() < 0) { //negative side runs -11 - -15
496  if (innerFPGAStub_)
497  stubIDs[innerFPGAStub_->disk().value() - 10] = innerFPGAStub_->phiregionaddress() + (1 << 10);
498  if (middleFPGAStub_)
499  stubIDs[middleFPGAStub_->layer().value() + 1] = middleFPGAStub_->phiregionaddress() + (1 << 10);
500  if (outerFPGAStub_)
501  stubIDs[outerFPGAStub_->layer().value() + 1] = outerFPGAStub_->phiregionaddress() + (1 << 10);
502  } else { // positive side runs 11-15]
503  if (innerFPGAStub_)
504  stubIDs[innerFPGAStub_->disk().value() + 10] = innerFPGAStub_->phiregionaddress() + (1 << 10);
505  if (middleFPGAStub_)
506  stubIDs[middleFPGAStub_->layer().value() + 1] = middleFPGAStub_->phiregionaddress() + (1 << 10);
507  if (outerFPGAStub_)
508  stubIDs[outerFPGAStub_->layer().value() + 1] = outerFPGAStub_->phiregionaddress() + (1 << 10);
509  }
510  }
511 
512  return stubIDs;
513 }
constexpr int N_DISK
Definition: Settings.h:22
const FPGAWord & fpgastubid() const
Definition: Residual.h:42
assert(be >=bs)
const FPGAWord & disk() const
Definition: Stub.h:68
unsigned int phiregionaddress() const
Definition: Stub.cc:112
int value() const
Definition: FPGAWord.h:24
L1TStub * l1tstub()
Definition: Stub.h:77
int disk() const
Definition: L1TStub.h:46
const FPGAWord & stubindex() const
Definition: Stub.h:66
Residual resid_[N_LAYER+N_DISK]
Definition: Tracklet.h:272
int nbits() const
Definition: FPGAWord.h:25
bool valid() const
Definition: Residual.h:30
const FPGAWord & itfit() const
Definition: Tracklet.h:154
const Stub * outerFPGAStub_
Definition: Tracklet.h:237
const Stub * innerFPGAStub_
Definition: Tracklet.h:235
const FPGAWord & layer() const
Definition: Stub.h:67
const Stub * middleFPGAStub_
Definition: Tracklet.h:236
constexpr unsigned int N_BITSMEMADDRESS
Definition: Settings.h:39
constexpr int N_LAYER
Definition: Settings.h:21

◆ getTrack()

Track* trklet::Tracklet::getTrack ( )
inline

Definition at line 192 of file Tracklet.h.

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

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

192  {
193  assert(fpgatrack_ != nullptr);
194  return fpgatrack_.get();
195  }
assert(be >=bs)
std::unique_ptr< Track > fpgatrack_
Definition: Tracklet.h:268

◆ ichiSqfit()

FPGAWord trklet::Tracklet::ichiSqfit ( ) const
inline

Definition at line 156 of file Tracklet.h.

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

156  {
157  return FPGAWord(ichisqrphifit_.value() + ichisqrzfit_.value(), ichisqrphifit_.nbits());
158  }
FPGAWord ichisqrzfit_
Definition: Tracklet.h:256
FPGAWord ichisqrphifit_
Definition: Tracklet.h:255
int value() const
Definition: FPGAWord.h:24
int nbits() const
Definition: FPGAWord.h:25

◆ id0fit()

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

Definition at line 153 of file Tracklet.h.

References fpgafitpars_.

Referenced by setFitPars().

153 { return fpgafitpars_.d0(); }
TrackPars< FPGAWord > fpgafitpars_
Definition: Tracklet.h:254

◆ innerFPGAStub()

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

Definition at line 61 of file Tracklet.h.

References innerFPGAStub_.

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

61 { return innerFPGAStub_; }
const Stub * innerFPGAStub_
Definition: Tracklet.h:235

◆ iphi0fit()

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

Definition at line 152 of file Tracklet.h.

References fpgafitpars_.

Referenced by setFitPars().

152 { return fpgafitpars_.phi0(); }
TrackPars< FPGAWord > fpgafitpars_
Definition: Tracklet.h:254

◆ irinvfit()

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

Definition at line 151 of file Tracklet.h.

References fpgafitpars_.

Referenced by setFitPars().

151 { return fpgafitpars_.rinv(); }
TrackPars< FPGAWord > fpgafitpars_
Definition: Tracklet.h:254

◆ isBarrel()

bool trklet::Tracklet::isBarrel ( ) const
inline

Definition at line 202 of file Tracklet.h.

References barrel_.

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

202 { return barrel_; }

◆ isDisk()

int trklet::Tracklet::isDisk ( ) const
inline

Definition at line 204 of file Tracklet.h.

References disk_.

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

204 { return disk_; }

◆ isOverlap()

bool trklet::Tracklet::isOverlap ( ) const
inline

Definition at line 203 of file Tracklet.h.

References overlap_.

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

203 { return overlap_; }

◆ itfit()

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

Definition at line 154 of file Tracklet.h.

References fpgafitpars_.

Referenced by getStubIDs(), and setFitPars().

154 { return fpgafitpars_.t(); }
TrackPars< FPGAWord > fpgafitpars_
Definition: Tracklet.h:254

◆ iz0fit()

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

Definition at line 155 of file Tracklet.h.

References fpgafitpars_.

Referenced by setFitPars().

155 { return fpgafitpars_.z0(); }
TrackPars< FPGAWord > fpgafitpars_
Definition: Tracklet.h:254

◆ layer()

int Tracklet::layer ( ) const

Definition at line 773 of file Tracklet.cc.

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

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

773  {
774  int l1 = (innerFPGAStub_ && innerFPGAStub_->layerdisk() < N_LAYER) ? innerFPGAStub_->l1tstub()->layerdisk() + 1 : 999,
776  : 999,
778  l = min(min(l1, l2), l3);
779  return (l < 999 ? l : 0);
780 }
int layerdisk() const
Definition: L1TStub.h:112
L1TStub * l1tstub()
Definition: Stub.h:77
unsigned int layerdisk() const
Definition: Stub.cc:185
const Stub * outerFPGAStub_
Definition: Tracklet.h:237
const Stub * innerFPGAStub_
Definition: Tracklet.h:235
const Stub * middleFPGAStub_
Definition: Tracklet.h:236
constexpr int N_LAYER
Definition: Settings.h:21

◆ layerstubstr()

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, trackIndex_, and validateGeometry_cfg::valid.

Referenced by trackfitstr().

572  {
573  assert(layer < N_LAYER);
574 
575  std::stringstream oss("");
576  if (!resid_[layer].valid())
577  oss << "0|0000000|0000000000|0000000|000000000000|000000000";
578  else {
579  if (trackIndex_ < 0 || trackIndex_ > (int)settings_.ntrackletmax()) {
580  cout << "trackIndex_ = " << trackIndex_ << endl;
581  assert(0);
582  }
583  const FPGAWord tmp(trackIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
584  oss << "1|"; // valid bit
585  oss << tmp.str() << "|";
586  oss << resid_[layer].fpgastubid().str() << "|";
587  oss << resid_[layer].stubptr()->r().str() << "|";
588  oss << resid_[layer].fpgaphiresid().str() << "|";
589  oss << resid_[layer].fpgarzresid().str();
590  }
591 
592  return oss.str();
593 }
const FPGAWord & r() const
Definition: Stub.h:60
unsigned int nbitstrackletindex() const
Definition: Settings.h:281
const FPGAWord & fpgastubid() const
Definition: Residual.h:42
const FPGAWord & fpgaphiresid() const
Definition: Residual.h:32
assert(be >=bs)
Settings const & settings_
Definition: Tracklet.h:274
const FPGAWord & fpgarzresid() const
Definition: Residual.h:37
Residual resid_[N_LAYER+N_DISK]
Definition: Tracklet.h:272
const Stub * stubptr() const
Definition: Residual.h:67
unsigned int ntrackletmax() const
Definition: Settings.h:343
int layer() const
Definition: Tracklet.cc:773
tmp
align.sh
Definition: createJobs.py:716
std::string str() const
Definition: FPGAWord.cc:54
constexpr int N_LAYER
Definition: Settings.h:21

◆ makeTrack()

Track Tracklet::makeTrack ( const std::vector< const L1TStub *> &  l1stubs)

Definition at line 740 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().

740  {
741  assert(fit());
742 
743  // Digitized track helix params
744  TrackPars<int> ipars(fpgafitpars_.rinv().value(),
745  fpgafitpars_.phi0().value(),
746  fpgafitpars_.d0().value(),
747  fpgafitpars_.t().value(),
748  fpgafitpars_.z0().value());
749 
750  // If fitter produced no stub list, take it from original tracklet.
751  vector<const L1TStub*> tmp = l1stubs.empty() ? getL1Stubs() : l1stubs;
752 
753  vector<L1TStub> tmp2;
754 
755  tmp2.reserve(tmp.size());
756  for (auto stub : tmp) {
757  tmp2.push_back(*stub);
758  }
759 
760  Track tmpTrack(ipars,
764  chisqrzfit_,
765  hitpattern_,
766  getStubIDs(),
767  tmp2,
768  getISeed());
769 
770  return tmpTrack;
771 }
double chisqrzfit_
Definition: Tracklet.h:260
assert(be >=bs)
std::vector< const L1TStub * > getL1Stubs()
Definition: Tracklet.cc:343
FPGAWord ichisqrzfit_
Definition: Tracklet.h:256
bool fit() const
Definition: Tracklet.h:197
FPGAWord ichisqrphifit_
Definition: Tracklet.h:255
TrackPars< FPGAWord > fpgafitpars_
Definition: Tracklet.h:254
int value() const
Definition: FPGAWord.h:24
int getISeed() const
Definition: Tracklet.cc:801
tmp
align.sh
Definition: createJobs.py:716
double chisqrphifit_
Definition: Tracklet.h:259
std::map< int, int > getStubIDs()
Definition: Tracklet.cc:361

◆ match()

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::MatchCalculator::execute(), trklet::FitTrack::execute(), trklet::MatchProcessor::matchCalculator(), and trklet::FitTrack::trackFitChisq().

105  {
106  assert(layerdisk < N_LAYER + N_DISK);
107  return resid_[layerdisk].valid();
108  }
constexpr int N_DISK
Definition: Settings.h:22
assert(be >=bs)
Residual resid_[N_LAYER+N_DISK]
Definition: Tracklet.h:272
bool valid() const
Definition: Residual.h:30
constexpr int N_LAYER
Definition: Settings.h:21

◆ middleFPGAStub()

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

Definition at line 63 of file Tracklet.h.

References middleFPGAStub_.

Referenced by Tracklet().

63 { return middleFPGAStub_; }
const Stub * middleFPGAStub_
Definition: Tracklet.h:236

◆ outerFPGAStub()

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

Definition at line 65 of file Tracklet.h.

References outerFPGAStub_.

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

65 { return outerFPGAStub_; }
const Stub * outerFPGAStub_
Definition: Tracklet.h:237

◆ phi0()

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(), trklet::FitTrack::trackFitFake(), and Tracklet().

121 { return trackpars_.phi0(); }
const T & phi0() const
Definition: TrackPars.h:30
TrackPars< double > trackpars_
Definition: Tracklet.h:246

◆ phi0approx()

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(), trklet::FitTrack::trackFitFake(), and Tracklet().

127 { return trackparsapprox_.phi0(); }
const T & phi0() const
Definition: TrackPars.h:30
TrackPars< double > trackparsapprox_
Definition: Tracklet.h:247

◆ phi0fit()

double trklet::Tracklet::phi0fit ( ) const
inline

Definition at line 139 of file Tracklet.h.

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

Referenced by setFitPars().

139 { return fitpars_.phi0(); }
const T & phi0() const
Definition: TrackPars.h:30
TrackPars< double > fitpars_
Definition: Tracklet.h:258

◆ phi0fitexact()

double trklet::Tracklet::phi0fitexact ( ) const
inline

Definition at line 146 of file Tracklet.h.

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

Referenced by setFitPars().

146 { return fitparsexact_.phi0(); }
const T & phi0() const
Definition: TrackPars.h:30
TrackPars< double > fitparsexact_
Definition: Tracklet.h:262

◆ proj()

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

◆ PSseed()

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

Definition at line 222 of file Tracklet.h.

References disk(), and layer().

Referenced by vmstrlayer().

222 { return ((layer() == 1) || (layer() == 2) || (disk() != 0)) ? 1 : 0; }
int disk() const
Definition: Tracklet.cc:782
int layer() const
Definition: Tracklet.cc:773

◆ resid()

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, resid_, and validateGeometry_cfg::valid.

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

110  {
111  assert(layerdisk < N_LAYER + N_DISK);
112  assert(resid_[layerdisk].valid());
113  return resid_[layerdisk];
114  }
constexpr int N_DISK
Definition: Settings.h:22
assert(be >=bs)
Residual resid_[N_LAYER+N_DISK]
Definition: Tracklet.h:272
constexpr int N_LAYER
Definition: Settings.h:21

◆ rinv()

double trklet::Tracklet::rinv ( ) const
inline

◆ rinvapprox()

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(), trklet::FitTrack::trackFitFake(), and Tracklet().

126 { return trackparsapprox_.rinv(); }
TrackPars< double > trackparsapprox_
Definition: Tracklet.h:247
const T & rinv() const
Definition: TrackPars.h:29

◆ rinvfit()

double trklet::Tracklet::rinvfit ( ) const
inline

Definition at line 138 of file Tracklet.h.

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

Referenced by setFitPars().

138 { return fitpars_.rinv(); }
TrackPars< double > fitpars_
Definition: Tracklet.h:258
const T & rinv() const
Definition: TrackPars.h:29

◆ rinvfitexact()

double trklet::Tracklet::rinvfitexact ( ) const
inline

Definition at line 145 of file Tracklet.h.

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

Referenced by setFitPars().

145 { return fitparsexact_.rinv(); }
const T & rinv() const
Definition: TrackPars.h:29
TrackPars< double > fitparsexact_
Definition: Tracklet.h:262

◆ seedIndex()

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

Definition at line 224 of file Tracklet.h.

References seedIndex_.

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

224 { return seedIndex_; }
unsigned int seedIndex_
Definition: Tracklet.h:227

◆ setFitPars()

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_, d0fit(), d0fitexact(), fitpars_, fitparsexact_, fpgafitpars_, fpgatrack_, hitpattern_, ichisqrphifit_, ichisqrzfit_, id0fit(), trklet::TrackPars< T >::init(), iphi0fit(), irinvfit(), itfit(), iz0fit(), makeTrack(), trklet::Settings::nbitsz0(), phi0fit(), phi0fitexact(), rinvfit(), rinvfitexact(), trklet::FPGAWord::set(), settings_, tfit(), tfitexact(), z0fit(), and z0fitexact().

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

537  {
539  chisqrphifit_ = chisqrphifit;
540  chisqrzfit_ = chisqrzfit;
541 
543  chisqrphifitexact_ = chisqrphifitexact;
544  chisqrzfitexact_ = chisqrzfitexact;
545 
546  if (irinvfit > (1 << 14))
547  irinvfit = (1 << 14);
548  if (irinvfit <= -(1 << 14))
549  irinvfit = -(1 << 14) + 1;
550  fpgafitpars_.rinv().set(irinvfit, 15, false, __LINE__, __FILE__);
551  fpgafitpars_.phi0().set(iphi0fit, 19, false, __LINE__, __FILE__);
552  fpgafitpars_.d0().set(id0fit, 19, false, __LINE__, __FILE__);
553  fpgafitpars_.t().set(itfit, 14, false, __LINE__, __FILE__);
554 
555  if (iz0fit >= (1 << (settings_.nbitsz0() - 1))) {
556  iz0fit = (1 << (settings_.nbitsz0() - 1)) - 1;
557  }
558 
559  if (iz0fit <= -(1 << (settings_.nbitsz0() - 1))) {
560  iz0fit = 1 - (1 << (settings_.nbitsz0() - 1));
561  }
562 
563  fpgafitpars_.z0().set(iz0fit, settings_.nbitsz0(), false, __LINE__, __FILE__);
564  ichisqrphifit_.set(ichisqrphifit, 8, true, __LINE__, __FILE__);
565  ichisqrzfit_.set(ichisqrzfit, 8, true, __LINE__, __FILE__);
566 
567  hitpattern_ = hitpattern;
568 
569  fpgatrack_ = std::make_unique<Track>(makeTrack(l1stubs));
570 }
double d0fitexact() const
Definition: Tracklet.h:147
const FPGAWord & iz0fit() const
Definition: Tracklet.h:155
double rinvfitexact() const
Definition: Tracklet.h:145
double d0fit() const
Definition: Tracklet.h:140
double tfit() const
Definition: Tracklet.h:141
double chisqrzfit_
Definition: Tracklet.h:260
double chisqrzfitexact_
Definition: Tracklet.h:264
const FPGAWord & iphi0fit() const
Definition: Tracklet.h:152
Settings const & settings_
Definition: Tracklet.h:274
FPGAWord ichisqrzfit_
Definition: Tracklet.h:256
FPGAWord ichisqrphifit_
Definition: Tracklet.h:255
TrackPars< double > fitpars_
Definition: Tracklet.h:258
double z0fitexact() const
Definition: Tracklet.h:149
const FPGAWord & id0fit() const
Definition: Tracklet.h:153
TrackPars< FPGAWord > fpgafitpars_
Definition: Tracklet.h:254
std::unique_ptr< Track > fpgatrack_
Definition: Tracklet.h:268
double tfitexact() const
Definition: Tracklet.h:148
double phi0fitexact() const
Definition: Tracklet.h:146
double z0fit() const
Definition: Tracklet.h:142
void set(int value, int nbits, bool positive=true, int line=-1, const char *file=nullptr)
Definition: FPGAWord.cc:14
TrackPars< double > fitparsexact_
Definition: Tracklet.h:262
const FPGAWord & itfit() const
Definition: Tracklet.h:154
Track makeTrack(const std::vector< const L1TStub *> &l1stubs)
Definition: Tracklet.cc:740
int nbitsz0() const
Definition: Settings.h:349
double chisqrphifitexact_
Definition: Tracklet.h:263
const FPGAWord & irinvfit() const
Definition: Tracklet.h:151
double phi0fit() const
Definition: Tracklet.h:139
double rinvfit() const
Definition: Tracklet.h:138
double chisqrphifit_
Definition: Tracklet.h:259
void init(T rinv, T phi0, T d0, T t, T z0)
Definition: TrackPars.h:21

◆ setTCIndex()

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

◆ setTrackIndex()

void Tracklet::setTrackIndex ( int  index)

Definition at line 813 of file Tracklet.cc.

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

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

813  {
814  trackIndex_ = index;
815  assert(index <= (int)settings_.ntrackletmax());
816 }
assert(be >=bs)
Settings const & settings_
Definition: Tracklet.h:274
unsigned int ntrackletmax() const
Definition: Settings.h:343

◆ setTrackletIndex()

void Tracklet::setTrackletIndex ( unsigned int  index)

◆ stubtruthmatch()

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().

138  {
139  set<int> tpset;
140  set<int> tpsetstub;
141  set<int> tpsetstubinner;
142  set<int> tpsetstubouter;
143 
144  vector<int> tps = stub->tps();
145  for (auto tp : tps) {
146  if (tp != 0) {
147  tpsetstub.insert(tp);
148  tpset.insert(abs(tp));
149  }
150  }
151  tps = innerFPGAStub_->l1tstub()->tps();
152  for (auto tp : tps) {
153  if (tp != 0) {
154  tpsetstubinner.insert(tp);
155  tpset.insert(abs(tp));
156  }
157  }
158  tps = outerFPGAStub_->l1tstub()->tps();
159  for (auto tp : tps) {
160  if (tp != 0) {
161  tpsetstubouter.insert(tp);
162  tpset.insert(abs(tp));
163  }
164  }
165 
166  for (auto tp : tpset) {
167  if (tpsetstub.find(tp) != tpsetstub.end() && tpsetstub.find(-tp) != tpsetstub.end() &&
168  tpsetstubinner.find(tp) != tpsetstubinner.end() && tpsetstubinner.find(-tp) != tpsetstubinner.end() &&
169  tpsetstubouter.find(tp) != tpsetstubouter.end() && tpsetstubouter.find(-tp) != tpsetstubouter.end()) {
170  return true;
171  }
172  }
173 
174  return false;
175 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
L1TStub * l1tstub()
Definition: Stub.h:77
const Stub * outerFPGAStub_
Definition: Tracklet.h:237
const Stub * innerFPGAStub_
Definition: Tracklet.h:235
std::vector< int > tps() const
Definition: L1TStub.h:81

◆ t()

double trklet::Tracklet::t ( ) const
inline

◆ tapprox()

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(), trklet::FitTrack::trackFitFake(), and Tracklet().

129 { return trackparsapprox_.t(); }
TrackPars< double > trackparsapprox_
Definition: Tracklet.h:247
const T & t() const
Definition: TrackPars.h:32

◆ TCID()

int trklet::Tracklet::TCID ( ) const
inline

◆ TCIndex()

int trklet::Tracklet::TCIndex ( ) const
inline

Definition at line 212 of file Tracklet.h.

References TCIndex_.

212 { return TCIndex_; }

◆ tfit()

double trklet::Tracklet::tfit ( ) const
inline

Definition at line 141 of file Tracklet.h.

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

Referenced by setFitPars().

141 { return fitpars_.t(); }
TrackPars< double > fitpars_
Definition: Tracklet.h:258
const T & t() const
Definition: TrackPars.h:32

◆ tfitexact()

double trklet::Tracklet::tfitexact ( ) const
inline

Definition at line 148 of file Tracklet.h.

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

Referenced by setFitPars().

148 { return fitparsexact_.t(); }
TrackPars< double > fitparsexact_
Definition: Tracklet.h:262
const T & t() const
Definition: TrackPars.h:32

◆ tpseed()

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().

107  {
108  set<int> tpset;
109 
110  set<int> tpsetstubinner;
111  set<int> tpsetstubouter;
112 
113  vector<int> tps = innerFPGAStub_->l1tstub()->tps();
114  for (auto tp : tps) {
115  if (tp != 0) {
116  tpsetstubinner.insert(tp);
117  tpset.insert(abs(tp));
118  }
119  }
120 
121  tps = outerFPGAStub_->l1tstub()->tps();
122  for (auto tp : tps) {
123  if (tp != 0) {
124  tpsetstubouter.insert(tp);
125  tpset.insert(abs(tp));
126  }
127  }
128 
129  for (auto& tp : tpset) {
130  if (tpsetstubinner.find(tp) != tpsetstubinner.end() && tpsetstubinner.find(-tp) != tpsetstubinner.end() &&
131  tpsetstubouter.find(tp) != tpsetstubouter.end() && tpsetstubouter.find(-tp) != tpsetstubouter.end()) {
132  return tp;
133  }
134  }
135  return 0;
136 }
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
L1TStub * l1tstub()
Definition: Stub.h:77
const Stub * outerFPGAStub_
Definition: Tracklet.h:237
const Stub * innerFPGAStub_
Definition: Tracklet.h:235
std::vector< int > tps() const
Definition: L1TStub.h:81

◆ trackfitstr()

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().

620  {
621  const unsigned maxNHits = 8;
622  const unsigned nBitsPerHit = 3;
623  vector<string> stub(maxNHits, "0");
624  string hitmap(maxNHits * nBitsPerHit, '0');
625 
626  // Assign stub strings for each of the possible projections for each seed.
627  // The specific layers/disks for a given seed are determined by the wiring.
628  switch (seedIndex()) {
629  case 0: // L1L2
630  stub[0] = layerstubstr(2); // L3
631  stub[1] = layerstubstr(3); // L4
632  stub[2] = layerstubstr(4); // L5
633  stub[3] = layerstubstr(5); // L6
634 
635  stub[4] = diskstubstr(0); // D1
636  stub[5] = diskstubstr(1); // D2
637  stub[6] = diskstubstr(2); // D3
638  stub[7] = diskstubstr(3); // D4
639 
640  break;
641 
642  case 1: // L2L3
643  stub[0] = layerstubstr(0); // L1
644  stub[1] = layerstubstr(3); // L4
645  stub[2] = layerstubstr(4); // L5
646 
647  stub[3] = diskstubstr(0); // D1
648  stub[4] = diskstubstr(1); // D2
649  stub[5] = diskstubstr(2); // D3
650  stub[6] = diskstubstr(3); // D4
651 
652  break;
653 
654  case 2: // L3L4
655  stub[0] = layerstubstr(0); // L1
656  stub[1] = layerstubstr(1); // L2
657  stub[2] = layerstubstr(4); // L5
658  stub[3] = layerstubstr(5); // L6
659 
660  stub[4] = diskstubstr(0); // D1
661  stub[5] = diskstubstr(1); // D2
662 
663  break;
664 
665  case 3: // L5L6
666  stub[0] = layerstubstr(0); // L1
667  stub[1] = layerstubstr(1); // L2
668  stub[2] = layerstubstr(2); // L3
669  stub[3] = layerstubstr(3); // L4
670 
671  break;
672 
673  case 4: // D1D2
674  stub[0] = layerstubstr(0); // L1
675  stub[1] = layerstubstr(1); // L2
676 
677  stub[2] = diskstubstr(2); // D3
678  stub[3] = diskstubstr(3); // D4
679  stub[4] = diskstubstr(4); // D5
680 
681  break;
682 
683  case 5: // D3D4
684  stub[0] = layerstubstr(0); // L1
685 
686  stub[1] = diskstubstr(0); // D1
687  stub[2] = diskstubstr(1); // D2
688  stub[3] = diskstubstr(4); // D5
689 
690  break;
691 
692  case 6: // L1D1
693  stub[0] = diskstubstr(1); // D2
694  stub[1] = diskstubstr(2); // D3
695  stub[2] = diskstubstr(3); // D4
696  stub[3] = diskstubstr(4); // D5
697 
698  break;
699 
700  case 7: // L2D1
701  stub[0] = layerstubstr(0); // L1
702 
703  stub[1] = diskstubstr(1); // D2
704  stub[2] = diskstubstr(2); // D3
705  stub[3] = diskstubstr(3); // D4
706 
707  break;
708  }
709 
710  // Only one hit per layer/disk is allowed currently, so the hit map for a
711  // given layer/disk is just equal to the valid bit of the corresponding stub
712  // string, which is the first character.
713  for (unsigned i = 0; i < maxNHits; i++)
714  hitmap[i * nBitsPerHit + 2] = stub[i][0];
715 
716  std::string oss("");
717  //Binary print out
718  if (!settings_.writeoutReal()) {
719  const FPGAWord tmp(getISeed(), settings_.nbitsseed(), true, __LINE__, __FILE__);
720 
721  oss += "1|"; // valid bit
722  oss += tmp.str() + "|";
723  oss += fpgapars_.rinv().str() + "|";
724  oss += fpgapars_.phi0().str() + "|";
725  oss += fpgapars_.z0().str() + "|";
726  oss += fpgapars_.t().str() + "|";
727  oss += hitmap;
728  for (unsigned i = 0; i < maxNHits; i++)
729  // If a valid stub string was never assigned, then that stub is not
730  // included in the output.
731  if (stub[i] != "0")
732  oss += "|" + stub[i];
733  }
734 
735  return oss;
736 }
TrackPars< FPGAWord > fpgapars_
Definition: Tracklet.h:244
unsigned int seedIndex() const
Definition: Tracklet.h:224
const std::string layerstubstr(const unsigned layer) const
Definition: Tracklet.cc:572
unsigned int nbitsseed() const
Definition: Settings.h:285
Settings const & settings_
Definition: Tracklet.h:274
const std::string diskstubstr(const unsigned disk) const
Definition: Tracklet.cc:595
int getISeed() const
Definition: Tracklet.cc:801
bool writeoutReal() const
Definition: Settings.h:205
tmp
align.sh
Definition: createJobs.py:716

◆ trackIndex()

int Tracklet::trackIndex ( ) const

Definition at line 818 of file Tracklet.cc.

References trackIndex_.

818 { return trackIndex_; }

◆ trackletIndex()

int trklet::Tracklet::trackletIndex ( ) const
inline

Definition at line 208 of file Tracklet.h.

References trackletIndex_.

Referenced by trklet::FullMatchMemory::addMatch().

208 { return trackletIndex_; }

◆ trackletparstr()

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(), cond::impl::to_string(), and trklet::Settings::writeoutReal().

188  {
189  if (settings_.writeoutReal()) {
190  std::string oss = std::to_string(fpgapars_.rinv().value() * settings_.krinvpars()) + " " +
191  std::to_string(fpgapars_.phi0().value() * settings_.kphi0pars()) + " " +
192  std::to_string(fpgapars_.d0().value() * settings_.kd0pars()) + " " +
193  std::to_string(fpgapars_.z0().value() * settings_.kz()) + " " +
194  std::to_string(fpgapars_.t().value() * settings_.ktpars());
195  return oss;
196  } else {
198  if (middleFPGAStub_) {
199  str += middleFPGAStub_->stubindex().str() + "|";
200  }
201  str += outerFPGAStub_->stubindex().str() + "|" + fpgapars_.rinv().str() + "|" + fpgapars_.phi0().str() + "|";
202  if (middleFPGAStub_)
203  str += fpgapars_.d0().str() + "|";
204  str += fpgapars_.z0().str() + "|" + fpgapars_.t().str();
205  return str;
206  }
207 }
TrackPars< FPGAWord > fpgapars_
Definition: Tracklet.h:244
double kz() const
Definition: Settings.h:314
std::string to_string(const V &value)
Definition: OMSAccess.h:77
double kd0pars() const
Definition: Settings.h:404
double ktpars() const
Definition: Settings.h:402
double krinvpars() const
Definition: Settings.h:397
Settings const & settings_
Definition: Tracklet.h:274
double kphi0pars() const
Definition: Settings.h:401
const FPGAWord & stubindex() const
Definition: Stub.h:66
bool writeoutReal() const
Definition: Settings.h:205
const Stub * outerFPGAStub_
Definition: Tracklet.h:237
const Stub * innerFPGAStub_
Definition: Tracklet.h:235
const Stub * middleFPGAStub_
Definition: Tracklet.h:236
#define str(s)
std::string str() const
Definition: FPGAWord.cc:54

◆ trackletprojstr()

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(), layer(), 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().

259  {
260  assert(layer > 0 && layer <= N_LAYER);
261  FPGAWord tmp;
262  if (trackletIndex_ < 0 || trackletIndex_ > (int)settings_.ntrackletmax()) {
263  throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__ << " trackletIndex_ = " << trackletIndex_;
264  }
265  tmp.set(trackletIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
266  FPGAWord tcid;
267  tcid.set(TCIndex_, settings_.nbitstcindex(), true, __LINE__, __FILE__);
268 
269  std::string oss = tcid.str() + "|" + tmp.str() + "|" + proj_[layer - 1].fpgaphiproj().str() + "|" +
270  proj_[layer - 1].fpgarzproj().str() + "|" + proj_[layer - 1].fpgaphiprojder().str() + "|" +
271  proj_[layer - 1].fpgarzprojder().str();
272  return oss;
273 }
Projection proj_[N_LAYER+N_DISK]
Definition: Tracklet.h:270
unsigned int nbitstrackletindex() const
Definition: Settings.h:281
const FPGAWord & fpgarzprojder() const
Definition: Projection.h:54
assert(be >=bs)
const FPGAWord & fpgaphiproj() const
Definition: Projection.h:39
Settings const & settings_
Definition: Tracklet.h:274
const FPGAWord & fpgaphiprojder() const
Definition: Projection.h:49
unsigned int nbitstcindex() const
Definition: Settings.h:286
void set(int value, int nbits, bool positive=true, int line=-1, const char *file=nullptr)
Definition: FPGAWord.cc:14
const FPGAWord & fpgarzproj() const
Definition: Projection.h:44
unsigned int ntrackletmax() const
Definition: Settings.h:343
int layer() const
Definition: Tracklet.cc:773
tmp
align.sh
Definition: createJobs.py:716
std::string str() const
Definition: FPGAWord.cc:54
constexpr int N_LAYER
Definition: Settings.h:21

◆ trackletprojstrD()

std::string Tracklet::trackletprojstrD ( int  disk) const

Definition at line 275 of file Tracklet.cc.

References funct::abs(), cms::cuda::assert(), disk(), 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().

275  {
276  assert(abs(disk) <= N_DISK);
277  FPGAWord tmp;
278  if (trackletIndex_ < 0 || trackletIndex_ > (int)settings_.ntrackletmax()) {
279  throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__ << " trackletIndex_ = " << trackletIndex_;
280  }
281  tmp.set(trackletIndex_, settings_.nbitstrackletindex(), true, __LINE__, __FILE__);
282  FPGAWord tcid;
283  if (settings_.extended()) {
284  tcid.set(TCIndex_, 8, true, __LINE__, __FILE__);
285  } else {
286  tcid.set(TCIndex_, 7, true, __LINE__, __FILE__);
287  }
288  std::string oss = tcid.str() + "|" + tmp.str() + "|" + proj_[N_LAYER + abs(disk) - 1].fpgaphiproj().str() + "|" +
289  proj_[N_LAYER + abs(disk) - 1].fpgarzproj().str() + "|" +
290  proj_[N_LAYER + abs(disk) - 1].fpgaphiprojder().str() + "|" +
291  proj_[N_LAYER + abs(disk) - 1].fpgarzprojder().str();
292  return oss;
293 }
Projection proj_[N_LAYER+N_DISK]
Definition: Tracklet.h:270
constexpr int N_DISK
Definition: Settings.h:22
int disk() const
Definition: Tracklet.cc:782
unsigned int nbitstrackletindex() const
Definition: Settings.h:281
const FPGAWord & fpgarzprojder() const
Definition: Projection.h:54
assert(be >=bs)
const FPGAWord & fpgaphiproj() const
Definition: Projection.h:39
Settings const & settings_
Definition: Tracklet.h:274
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
const FPGAWord & fpgaphiprojder() const
Definition: Projection.h:49
void set(int value, int nbits, bool positive=true, int line=-1, const char *file=nullptr)
Definition: FPGAWord.cc:14
const FPGAWord & fpgarzproj() const
Definition: Projection.h:44
bool extended() const
Definition: Settings.h:256
unsigned int ntrackletmax() const
Definition: Settings.h:343
tmp
align.sh
Definition: createJobs.py:716
std::string str() const
Definition: FPGAWord.cc:54
constexpr int N_LAYER
Definition: Settings.h:21

◆ trackletprojstrdisk()

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

Definition at line 80 of file Tracklet.h.

References disk(), and trackletprojstrD().

80 { return trackletprojstrD(disk); }
int disk() const
Definition: Tracklet.cc:782
std::string trackletprojstrD(int disk) const
Definition: Tracklet.cc:275

◆ trackletprojstrlayer()

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

Definition at line 79 of file Tracklet.h.

References layer(), and trackletprojstr().

79 { return trackletprojstr(layer); }
std::string trackletprojstr(int layer) const
Definition: Tracklet.cc:259
int layer() const
Definition: Tracklet.cc:773

◆ validProj()

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

◆ vmstrdisk()

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

Definition at line 244 of file Tracklet.cc.

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

244  {
245  FPGAWord index;
246  if (allstubindex >= (1 << 7)) {
247  edm::LogPrint("Tracklet") << "Warning projection number too large!";
248  index.set((1 << 7) - 1, 7, true, __LINE__, __FILE__);
249  } else {
250  index.set(allstubindex, 7, true, __LINE__, __FILE__);
251  }
252  std::string oss =
253  index.str() + "|" + proj_[N_LAYER + disk - 1].fpgarzbin1projvm().str() + "|" +
254  proj_[N_LAYER + disk - 1].fpgarzbin2projvm().str() + "|" + proj_[N_LAYER + disk - 1].fpgafinerzvm().str() + "|" +
255  proj_[N_LAYER + disk - 1].fpgafinephivm().str() + "|" + proj_[N_LAYER + disk - 1].getBendIndex().str();
256  return oss;
257 }
const FPGAWord & fpgafinerzvm() const
Definition: Projection.h:69
Projection proj_[N_LAYER+N_DISK]
Definition: Tracklet.h:270
int disk() const
Definition: Tracklet.cc:782
const FPGAWord & fpgarzbin2projvm() const
Definition: Projection.h:64
const FPGAWord & getBendIndex() const
Definition: Projection.h:121
Log< level::Warning, true > LogPrint
const FPGAWord & fpgarzbin1projvm() const
Definition: Projection.h:59
const FPGAWord & fpgafinephivm() const
Definition: Projection.h:74
std::string str() const
Definition: FPGAWord.cc:54
constexpr int N_LAYER
Definition: Settings.h:21

◆ vmstrlayer()

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(), layer(), trklet::FPGAWord::nbits(), proj_, PSseed(), settings_, trklet::FPGAWord::str(), AlCaHLTBitMon_QueryRunRegistry::string, createJobs::tmp, cond::impl::to_string(), and trklet::FPGAWord::value().

209  {
210  FPGAWord index;
211  if (allstubindex >= (1 << 7)) {
212  edm::LogPrint("Tracklet") << "Warning projection number too large!";
213  index.set((1 << 7) - 1, 7, true, __LINE__, __FILE__);
214  } else {
215  index.set(allstubindex, 7, true, __LINE__, __FILE__);
216  }
217 
218  // This is a shortcut.
219  //int irinvvm=16+(fpgarinv().value()>>(fpgarinv().nbits()-5));
220  // rinv is not directly available in the TrackletProjection.
221  // can be inferred from phi derivative: rinv = - phider * 2
222  int tmp_irinv = proj_[layer - 1].fpgaphiprojder().value() * (-2);
223  int nbits_irinv = proj_[layer - 1].fpgaphiprojder().nbits() + 1;
224 
225  // irinv in VMProjection:
226  // top 5 bits of rinv and shifted to be positive
227  int irinvvm = 16 + (tmp_irinv >> (nbits_irinv - 5));
228 
229  if (settings_.extended() && (irinvvm > 31)) { //TODO - displaced tracking should protect against this
230  edm::LogPrint("Tracklet") << "Warning irinvvm too large:" << irinvvm;
231  irinvvm = 31;
232  }
233 
234  assert(irinvvm >= 0);
235  assert(irinvvm < 32);
236  FPGAWord tmp;
237  tmp.set(irinvvm, 5, true, __LINE__, __FILE__);
238  std::string oss = index.str() + "|" + proj_[layer - 1].fpgarzbin1projvm().str() + "|" +
239  proj_[layer - 1].fpgarzbin2projvm().str() + "|" + proj_[layer - 1].fpgafinerzvm().str() + "|" +
240  proj_[layer - 1].fpgafinephivm().str() + "|" + tmp.str() + "|" + std::to_string(PSseed());
241  return oss;
242 }
const FPGAWord & fpgafinerzvm() const
Definition: Projection.h:69
Projection proj_[N_LAYER+N_DISK]
Definition: Tracklet.h:270
const FPGAWord & fpgarzbin2projvm() const
Definition: Projection.h:64
std::string to_string(const V &value)
Definition: OMSAccess.h:77
assert(be >=bs)
Settings const & settings_
Definition: Tracklet.h:274
int value() const
Definition: FPGAWord.h:24
Log< level::Warning, true > LogPrint
const FPGAWord & fpgarzbin1projvm() const
Definition: Projection.h:59
const FPGAWord & fpgaphiprojder() const
Definition: Projection.h:49
int nbits() const
Definition: FPGAWord.h:25
bool extended() const
Definition: Settings.h:256
const FPGAWord & fpgafinephivm() const
Definition: Projection.h:74
int layer() const
Definition: Tracklet.cc:773
unsigned int PSseed() const
Definition: Tracklet.h:222
tmp
align.sh
Definition: createJobs.py:716
std::string str() const
Definition: FPGAWord.cc:54

◆ z0()

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(), trklet::FitTrack::trackFitFake(), and Tracklet().

124 { return trackpars_.z0(); }
const T & z0() const
Definition: TrackPars.h:33
TrackPars< double > trackpars_
Definition: Tracklet.h:246

◆ z0approx()

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(), trklet::FitTrack::trackFitFake(), and Tracklet().

130 { return trackparsapprox_.z0(); }
const T & z0() const
Definition: TrackPars.h:33
TrackPars< double > trackparsapprox_
Definition: Tracklet.h:247

◆ z0fit()

double trklet::Tracklet::z0fit ( ) const
inline

Definition at line 142 of file Tracklet.h.

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

Referenced by setFitPars().

142 { return fitpars_.z0(); }
const T & z0() const
Definition: TrackPars.h:33
TrackPars< double > fitpars_
Definition: Tracklet.h:258

◆ z0fitexact()

double trklet::Tracklet::z0fitexact ( ) const
inline

Definition at line 149 of file Tracklet.h.

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

Referenced by setFitPars().

149 { return fitparsexact_.z0(); }
const T & z0() const
Definition: TrackPars.h:33
TrackPars< double > fitparsexact_
Definition: Tracklet.h:262

Member Data Documentation

◆ barrel_

bool trklet::Tracklet::barrel_
private

Definition at line 230 of file Tracklet.h.

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

◆ chisqrphifit_

double trklet::Tracklet::chisqrphifit_
private

Definition at line 259 of file Tracklet.h.

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

◆ chisqrphifitexact_

double trklet::Tracklet::chisqrphifitexact_
private

Definition at line 263 of file Tracklet.h.

Referenced by setFitPars().

◆ chisqrzfit_

double trklet::Tracklet::chisqrzfit_
private

Definition at line 260 of file Tracklet.h.

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

◆ chisqrzfitexact_

double trklet::Tracklet::chisqrzfitexact_
private

Definition at line 264 of file Tracklet.h.

Referenced by setFitPars().

◆ disk_

bool trklet::Tracklet::disk_
private

Definition at line 231 of file Tracklet.h.

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

◆ fitpars_

TrackPars<double> trklet::Tracklet::fitpars_
private

Definition at line 258 of file Tracklet.h.

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

◆ fitparsexact_

TrackPars<double> trklet::Tracklet::fitparsexact_
private

Definition at line 262 of file Tracklet.h.

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

◆ fpgafitpars_

TrackPars<FPGAWord> trklet::Tracklet::fpgafitpars_
private

Definition at line 254 of file Tracklet.h.

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

◆ fpgapars_

TrackPars<FPGAWord> trklet::Tracklet::fpgapars_
private

Definition at line 244 of file Tracklet.h.

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

◆ fpgatrack_

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

Definition at line 268 of file Tracklet.h.

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

◆ hitpattern_

int trklet::Tracklet::hitpattern_
private

Definition at line 266 of file Tracklet.h.

Referenced by makeTrack(), and setFitPars().

◆ ichisqrphifit_

FPGAWord trklet::Tracklet::ichisqrphifit_
private

Definition at line 255 of file Tracklet.h.

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

◆ ichisqrzfit_

FPGAWord trklet::Tracklet::ichisqrzfit_
private

Definition at line 256 of file Tracklet.h.

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

◆ innerFPGAStub_

const Stub* trklet::Tracklet::innerFPGAStub_
private

◆ middleFPGAStub_

const Stub* trklet::Tracklet::middleFPGAStub_
private

◆ outerFPGAStub_

const Stub* trklet::Tracklet::outerFPGAStub_
private

◆ overlap_

bool trklet::Tracklet::overlap_
private

Definition at line 232 of file Tracklet.h.

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

◆ proj_

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

◆ projdisk_

int trklet::Tracklet::projdisk_[N_DISK]
private

Definition at line 251 of file Tracklet.h.

Referenced by Tracklet().

◆ projlayer_

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

Definition at line 250 of file Tracklet.h.

Referenced by Tracklet().

◆ resid_

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

◆ seedIndex_

unsigned int trklet::Tracklet::seedIndex_
private

Definition at line 227 of file Tracklet.h.

Referenced by seedIndex(), and Tracklet().

◆ settings_

Settings const& trklet::Tracklet::settings_
private

◆ TCIndex_

int trklet::Tracklet::TCIndex_
private

◆ trackIndex_

int trklet::Tracklet::trackIndex_
private

Definition at line 241 of file Tracklet.h.

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

◆ trackletIndex_

int trklet::Tracklet::trackletIndex_
private

◆ trackpars_

TrackPars<double> trklet::Tracklet::trackpars_
private

Definition at line 246 of file Tracklet.h.

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

◆ trackparsapprox_

TrackPars<double> trklet::Tracklet::trackparsapprox_
private

Definition at line 247 of file Tracklet.h.

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

◆ triplet_

bool trklet::Tracklet::triplet_
private

Definition at line 233 of file Tracklet.h.

Referenced by Tracklet().