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

Enable backward search for iterations seeded by pixel quadruplets #387

Merged
merged 4 commits into from
Nov 25, 2021
Merged

Enable backward search for iterations seeded by pixel quadruplets #387

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
5fe555f
Enable backward search for iterations seeded by pixel quadruplets, in…
mmasciov Nov 20, 2021
1a683ce
Improved seed hit merging
mmasciov Nov 20, 2021
cf77c34
Rebase and cleanup of PR #384
mmasciov Nov 24, 2021
0be6710
Fix Event.cc seed cleaning function (unused) to match iteration-speci…
mmasciov Nov 24, 2021
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
9 changes: 1 addition & 8 deletions Config.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -342,14 +342,7 @@ namespace Config
// config on seed cleaning
constexpr float track1GeVradius = 87.6; // = 1/(c*B)
constexpr float c_etamax_brl = 0.9;
constexpr float c_dpt_brl_0 = 0.025;
constexpr float c_dpt_ec_0 = 0.0125;
constexpr float c_ptmax_0 = 2.0;
constexpr float c_dpt_1 = 0.10;
constexpr float c_ptmax_1 = 5.0;
constexpr float c_dpt_2 = 0.20;
constexpr float c_ptmax_2 = 10.0;
constexpr float c_dpt_3 = 0.25;
constexpr float c_dpt_common = 0.25;
constexpr float c_dzmax_brl = 0.005;
constexpr float c_drmax_brl = 0.010;
constexpr float c_ptmin_hpt = 2.0;
Expand Down
30 changes: 2 additions & 28 deletions Event.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -574,14 +574,7 @@ void Event::print_tracks(const TrackVec& tracks, bool print_hits) const
int Event::clean_cms_seedtracks(TrackVec *seed_ptr)
{
const float etamax_brl = Config::c_etamax_brl;
const float dpt_brl_0 = Config::c_dpt_brl_0;
const float dpt_ec_0 = Config::c_dpt_ec_0;
const float ptmax_0 = Config::c_ptmax_0;
const float dpt_1 = Config::c_dpt_1;
const float ptmax_1 = Config::c_ptmax_1;
const float dpt_2 = Config::c_dpt_2;
const float ptmax_2 = Config::c_ptmax_2;
const float dpt_3 = Config::c_dpt_3;
const float dpt_common = Config::c_dpt_common;
const float dzmax_brl = Config::c_dzmax_brl;
const float drmax_brl = Config::c_drmax_brl;
const float ptmin_hpt = Config::c_ptmin_hpt;
Expand Down Expand Up @@ -654,28 +647,9 @@ int Event::clean_cms_seedtracks(TrackVec *seed_ptr)

const float thisDPt = std::abs(Pt2-Pt1);
////// Require pT consistency between seeds. If dpT is large, do not remove seed-track.
////// Adaptive thresholds, based on pT of reference seed-track (choice is a compromise between efficiency and duplicate rate):
////// - 2.5% if track is barrel and w/ pT<2 GeV
////// - 1.25% if track is non-barrel and w/ pT<2 GeV
////// - 10% if track w/ 2<pT<5 GeV
////// - 20% if track w/ 5<pT<10 GeV
////// - 25% if track w/ pT>10 GeV
if(thisDPt>dpt_brl_0*(Pt1) && Pt1<ptmax_0 && std::abs(Eta1)<etamax_brl)
if( thisDPt > dpt_common*(Pt1) )
continue;

else if(thisDPt>dpt_ec_0*(Pt1) && Pt1<ptmax_0 && std::abs(Eta1)>etamax_brl)
continue;

else if(thisDPt>dpt_1*(Pt1) && Pt1>ptmax_0 && Pt1<ptmax_1)
continue;

else if(thisDPt>dpt_2*(Pt1) && Pt1>ptmax_1 && Pt1<ptmax_2)
continue;

else if(thisDPt>dpt_3*(Pt1) && Pt1>ptmax_2)
continue;


const float Eta2 = eta[tss];
const float deta2 = std::pow(Eta1-Eta2, 2);

Expand Down
Loading