Skip to content

Optimizations in the CM tracker #1642

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

Merged
sawenzel merged 9 commits into from
Jan 30, 2019
Merged

Optimizations in the CM tracker #1642

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
5d9f568
Printing the number of processed clusters
iouribelikov Jan 18, 2019
adee05e
Adding the estimation of fake-track rate
iouribelikov Jan 16, 2019
f77b2a4
Optimizing the cluster selection, and fixing the discontinuity at phi…
iouribelikov Jan 13, 2019
16b06af
Optimizing the search in Z, and making it faster
iouribelikov Jan 12, 2019
e08d65f
An std::map-based track-label cooking
iouribelikov Jan 25, 2019
a5cbc7d
A command-line argument for setting the number of threads
iouribelikov Jan 28, 2019
3ae84ae
Clang-formatting
iouribelikov Jan 28, 2019
5ad2003
Making the tracking parameters float. Using a meaningful name for the…
iouribelikov Jan 30, 2019
0aa2b8d
Little fixes suggested by Ruben
iouribelikov Jan 30, 2019
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
41 changes: 32 additions & 9 deletions Detectors/ITSMFT/ITS/macros/test/CheckTracks.C
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -71,11 +71,22 @@ void CheckTracks(std::string tracfile = "o2trac_its.root", std::string clusfile
Int_t tf = 0, nrec = 0;
Int_t lastEventID = -1;
Int_t nev = mcTree->GetEntries();
TH1D* den = new TH1D("den", ";#it{p}_{T} (GeV/#it{c});Den", 100, 0.01, 10);
TH1D* num = new TH1D("num", ";#it{p}_{T} (GeV/#it{c});Num", 100, 0.01, 10);

Int_t nb = 100;
Double_t xbins[nb + 1], ptcutl = 0.01, ptcuth = 10.;
Double_t a = TMath::Log(ptcuth / ptcutl) / nb;
for (Int_t i = 0; i <= nb; i++)
xbins[i] = ptcutl * TMath::Exp(i * a);
TH1D* num = new TH1D("num", ";#it{p}_{T} (GeV/#it{c});Efficiency (fake-track rate)", nb, xbins);
num->Sumw2();
TH1D* fak = new TH1D("fak", ";#it{p}_{T} (GeV/#it{c});Fak", nb, xbins);
fak->Sumw2();
TH1D* den = new TH1D("den", ";#it{p}_{T} (GeV/#it{c});Den", nb, xbins);
den->Sumw2();

for (Int_t n = 0; n < nev; n++) {
std::cout << "\nMC event " << n << '/' << nev << std::endl;
Int_t nGen = 0, nGoo = 0;
Int_t nGen = 0, nGoo = 0, nFak = 0;
mcTree->GetEvent(n);
Int_t nmc = mcArr->size();
Int_t nmcrefs = mcTrackRefs->size();
Expand Down Expand Up @@ -106,7 +117,7 @@ void CheckTracks(std::string tracfile = "o2trac_its.root", std::string clusfile

int ok = 0;
// Check the availability of clusters
for (int i = 0; i < clusArr->size(); i++) {
for (uint i = 0; i < clusArr->size(); i++) {
const Cluster& c = (*clusArr)[i];
auto lab = (clusLabArr->getLabels(i))[0];
if (lab.getEventID() != n)
Expand Down Expand Up @@ -190,15 +201,21 @@ void CheckTracks(std::string tracfile = "o2trac_its.root", std::string clusfile
if (label > 0) {
nGoo++; // Good found tracks for the efficiency calculation
num->Fill(mcPt);
} else {
nFak++; // Fake-track rate calculation
fak->Fill(mcPt);
}
}

nt->Fill( // mcYOut,recYOut,
mcZOut, recZOut, mcPhiOut, recPhiOut, mcThetaOut, recThetaOut, mcPhi, recPhi, mcLam, recLam, mcPt, recPt, ip[0],
ip[1], label);
}
Float_t eff = (nGen > 0) ? nGoo / Float_t(nGen) : -1.;
std::cout << "Good found tracks: " << nGoo << ", efficiency: " << eff << std::endl;
if (nGen > 0) {
Float_t eff = nGoo / Float_t(nGen);
Float_t rat = nFak / Float_t(nGen);
std::cout << "Good found tracks: " << nGoo << ", efficiency: " << eff << ", fake-track rate: " << rat << std::endl;
}
}

// "recPt>0" means "found tracks only"
Expand All @@ -215,9 +232,15 @@ void CheckTracks(std::string tracfile = "o2trac_its.root", std::string clusfile
nt->Draw("mcPhiOut-recPhiOut", "recPt>0 && label>0");
new TCanvas;
nt->Draw("mcThetaOut-recThetaOut", "recPt>0 && label>0");
new TCanvas;
num->Divide(den);
num->Draw();
TCanvas* c1 = new TCanvas;
c1->SetLogx();
c1->SetGridx();
c1->SetGridy();
num->Divide(num, den, 1, 1, "b");
num->Draw("histe");
fak->Divide(fak, den, 1, 1, "b");
fak->SetLineColor(2);
fak->Draw("histesame");
f->Write();
f->Close();
}
4 changes: 2 additions & 2 deletions Detectors/ITSMFT/ITS/reconstruction/include/ITSReconstruction/CookedTracker.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -145,7 +145,7 @@ class CookedTracker::Layer
void setR(Double_t r) { mR = r; }
void unloadClusters();
void selectClusters(std::vector<Int_t>& s, Float_t phi, Float_t dy, Float_t z, Float_t dz);
Int_t findClusterIndex(Double_t z) const;
Int_t findClusterIndex(Float_t z) const;
Float_t getR() const { return mR; }
const Cluster* getCluster(Int_t i) const { return mClusters[i]; }
Float_t getAlphaRef(Int_t i) const { return mAlphaRef[i]; }
Expand All @@ -161,7 +161,7 @@ class CookedTracker::Layer
std::vector<const Cluster*> mClusters; ///< All clusters
std::vector<Float_t> mAlphaRef; ///< alpha of the reference plane
std::vector<Float_t> mPhi; ///< cluster phi
std::vector<Int_t> mSectors[kNSectors]; ///< Cluster indices sector-by-sector
std::vector<std::pair<int, float>> mSectors[kNSectors]; ///< Cluster indices sector-by-sector
};
} // namespace ITS
} // namespace o2
Expand Down
Loading