-
Notifications
You must be signed in to change notification settings - Fork 6
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.
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
Updating Stub Format for VMSMERouter Changes #279
Changes from all commits
b1e6517
3fae98f
ee50448
db108e7
99ffbb9
f5a8118
f4e51ce
ca46e84
4758a36
c4a8c5d
3368d92
d46656e
8131c53
4f0e15a
a9d917f
8d510d2
c4fafbb
38c9f38
d84665a
874a688
41c3e31
e833919
02276dc
c52e649
File filter
Filter by extension
Conversations
Jump to
Diff view
Diff view
There are no files selected for viewing
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -1,4 +1,4 @@ | ||
To run the L1 tracking & create a TTree of tracking performance: | ||
To run the L1 tracking & create a TTree of tracking performance: | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. Remove extra space |
||
|
||
cmsRun L1TrackNtupleMaker_cfg.py | ||
|
||
|
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -88,6 +88,7 @@ namespace trklet { | |
void setTableTEDFile(std::string tableTEDFileName) { tableTEDFile_ = tableTEDFileName; } | ||
void setTableTREFile(std::string tableTREFileName) { tableTREFile_ = tableTREFileName; } | ||
|
||
unsigned int nndbitsstub(unsigned int layerdisk) const { return nndbitsstub_[layerdisk]; } | ||
unsigned int nzbitsstub(unsigned int layerdisk) const { return nzbitsstub_[layerdisk]; } | ||
unsigned int nphibitsstub(unsigned int layerdisk) const { return nphibitsstub_[layerdisk]; } | ||
unsigned int nrbitsstub(unsigned int layerdisk) const { return nrbitsstub_[layerdisk]; } | ||
|
@@ -170,12 +171,11 @@ namespace trklet { | |
double rphicut2S(unsigned int iSeed, unsigned int idisk) const { return rphicut2S_[idisk][iSeed]; } | ||
double rcut2S(unsigned int iSeed, unsigned int idisk) const { return rcut2S_[idisk][iSeed]; } | ||
|
||
|
||
unsigned int irmean(unsigned int iLayer) const { return irmean_[iLayer]; } | ||
double rmean(unsigned int iLayer) const { return irmean_[iLayer] * rmaxdisk_ / 4096; } | ||
double rmax(unsigned int iLayer) const { return rmean(iLayer) + drmax(); } | ||
double rmin(unsigned int iLayer) const { return rmean(iLayer) - drmax(); } | ||
unsigned int izmean(unsigned int iDisk) const { return izmean_[iDisk]; } | ||
unsigned int izmean(unsigned int iDisk) const { return izmean_[iDisk]; } | ||
double zmean(unsigned int iDisk) const { return izmean_[iDisk] * zlength_ / 2048; } | ||
double zmax(unsigned int iDisk) const { return zmean(iDisk) + dzmax(); } | ||
double zmin(unsigned int iDisk) const { return zmean(iDisk) - dzmax(); } | ||
|
@@ -344,7 +344,9 @@ namespace trklet { | |
|
||
double kz() const { return 2.0 * zlength_ / (1 << nzbitsstub_[0]); } | ||
double kz(unsigned int layerdisk) const { return 2.0 * zlength_ / (1 << nzbitsstub_[layerdisk]); } | ||
double kr() const { return rmaxdisk_ / (1 << nrbitsstub_[N_LAYER]); } | ||
double kr() const { | ||
return rmaxdisk_ / (1 << (nrbitsstub_[N_LAYER] + 1)); | ||
} // + 1 required to offset artificial decrease in # of diskps r bits from 12 -> 11 to make space for negDisk bit | ||
double krbarrel() const { return 2.0 * drmax() / (1 << nrbitsstub_[0]); } | ||
|
||
double maxrinv() const { return maxrinv_; } | ||
|
@@ -527,9 +529,10 @@ namespace trklet { | |
std::array<unsigned int, N_LAYER> irmean_{{851, 1269, 1784, 2347, 2936, 3697}}; | ||
std::array<unsigned int, N_DISK> izmean_{{2239, 2645, 3163, 3782, 4523}}; | ||
|
||
std::array<unsigned int, N_LAYER + N_DISK> nndbitsstub_{{0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1}}; | ||
std::array<unsigned int, N_LAYER + N_DISK> nzbitsstub_{{12, 12, 12, 8, 8, 8, 7, 7, 7, 7, 7}}; | ||
std::array<unsigned int, N_LAYER + N_DISK> nphibitsstub_{{14, 14, 14, 17, 17, 17, 14, 14, 14, 14, 14}}; | ||
std::array<unsigned int, N_LAYER + N_DISK> nrbitsstub_{{7, 7, 7, 7, 7, 7, 12, 12, 12, 12, 12}}; | ||
std::array<unsigned int, N_LAYER + N_DISK> nrbitsstub_{{7, 7, 7, 7, 7, 7, 11, 11, 11, 11, 11}}; | ||
|
||
unsigned int nrbitsprojderdisk_{9}; | ||
unsigned int nbitsphiprojderL123_{10}; | ||
|
@@ -1051,18 +1054,23 @@ namespace trklet { | |
double stripLength_PS_{0.1467}; | ||
double stripLength_2S_{5.0250}; | ||
|
||
// The DR binning below disabled, as doesn't match latest FW. | ||
|
||
//Following values are used for duplicate removal | ||
//Rinv bins were optimised to ensure a similar number of tracks in each bin prior to DR | ||
//Rinv bin edges for 6 bins. | ||
std::vector<double> rinvBins_{-rinvcut(), -0.004968, -0.003828, 0, 0.003828, 0.004968, rinvcut()}; | ||
//std::vector<double> rinvBins_{-rinvcut(), -0.004968, -0.003828, 0, 0.003828, 0.004968, rinvcut()}; | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. Are these changes related to this updated? |
||
std::vector<double> rinvBins_{-rinvcut(), rinvcut()}; | ||
//Phi bin edges for 2 bins. | ||
std::vector<double> phiBins_{0, dphisectorHG() / 2, dphisectorHG()}; | ||
//std::vector<double> phiBins_{0, dphisectorHG() / 2, dphisectorHG()}; | ||
std::vector<double> phiBins_{0, dphisectorHG()}; | ||
//Overlap size for the overlap rinv bins in DR | ||
double rinvOverlapSize_{0.0004}; | ||
//Overlap size for the overlap phi bins in DR | ||
double phiOverlapSize_{M_PI / 360}; | ||
//The maximum number of tracks that are compared to all the other tracks per rinv bin | ||
int numTracksComparedPerBin_{32}; | ||
//int numTracksComparedPerBin_{32}; | ||
int numTracksComparedPerBin_{9999}; | ||
|
||
double sensorSpacing_2S_{0.18}; | ||
}; | ||
|
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -221,6 +221,8 @@ namespace trklet { | |
|
||
unsigned int PSseed() const { return ((layer() == 1) || (layer() == 2) || (disk() != 0)) ? 1 : 0; } | ||
|
||
// Seed type: | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. Why is this added? |
||
// L1L2=0,L2L3=1,L3L4=2,L5L6=3,D1D2=4,D3D4=5,L1D1=6,L2D1=7 | ||
unsigned int seedIndex() const { return seedIndex_; } | ||
|
||
private: | ||
|
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is this related to these changes?