Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fix KF Track Parameter Digitization #127

Merged
merged 2 commits into from
Feb 3, 2022
Merged

Fix KF Track Parameter Digitization #127

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
2697f94
Fix KF Track Parameter Digitization
SClarkPhysics Feb 1, 2022
4c44924
Add comments to Track.h
SClarkPhysics Feb 2, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
13 changes: 7 additions & 6 deletions L1Trigger/TrackFindingTracklet/interface/Track.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -58,11 +58,12 @@ namespace trklet {
// Get floating point helix params by undigitized digi helix params
double phi0(Settings const& settings) const;

double eta(Settings const& settings) const { return asinh(ipars_.t() * settings.ktpars()); }
double tanL(Settings const& settings) const { return ipars_.t() * settings.ktpars(); }
double z0(Settings const& settings) const { return ipars_.z0() * settings.kz0pars(); }
double rinv(Settings const& settings) const { return ipars_.rinv() * settings.krinvpars(); }
double d0(Settings const& settings) const { return ipars_.d0() * settings.kd0pars(); }
//The following return the floating point track parameters by undigitizing the digitized parameters
double eta(Settings const& settings) const { return (asinh(ipars_.t() + 0.5) * settings.ktpars()); }
double tanL(Settings const& settings) const { return (ipars_.t() + 0.5) * settings.ktpars(); }
double z0(Settings const& settings) const { return (ipars_.z0() + 0.5) * settings.kz0pars(); }
double rinv(Settings const& settings) const { return (ipars_.rinv() + 0.5) * settings.krinvpars(); }
double d0(Settings const& settings) const { return (ipars_.d0() + 0.5) * settings.kd0pars(); }
double chisq() const { return chisqrphi_ + chisqrz_; }

double chisqrphi() const { return chisqrphi_; }
Expand All @@ -78,7 +79,7 @@ namespace trklet {
}

private:
TrackPars<int> ipars_;
TrackPars<int> ipars_; //digitized track parameters
int ichisqrphi_;
int ichisqrz_;

Expand Down
10 changes: 5 additions & 5 deletions L1Trigger/TrackFindingTracklet/src/HybridFit.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -212,11 +212,11 @@ void HybridFit::Fit(Tracklet* tracklet, std::vector<const Stub*>& trackstublist)
double phi0fit = reco::reduceRange(phi0 - iSector_ * 2 * M_PI / N_SECTOR + 0.5 * settings_.dphisectorHG());
double rinvfit = 0.01 * settings_.c() * settings_.bfield() * qoverpt;

int irinvfit = rinvfit / settings_.krinvpars();
int iphi0fit = phi0fit / settings_.kphi0pars();
int itanlfit = trk.tanLambda() / settings_.ktpars();
int iz0fit = trk.z0() / settings_.kz0pars();
int id0fit = d0 / settings_.kd0pars();
int irinvfit = floor(rinvfit / settings_.krinvpars());
int iphi0fit = floor(phi0fit / settings_.kphi0pars());
int itanlfit = floor(trk.tanLambda() / settings_.ktpars());
int iz0fit = floor(trk.z0() / settings_.kz0pars());
int id0fit = floor(d0 / settings_.kd0pars());
int ichi2rphifit = chi2rphi / 16;
int ichi2rzfit = trk.chi2rz() / 16;

Expand Down