Commit a91020b0 authored by Johnston's avatar Johnston
Browse files

Changes to be committed:

	modified:   HoughTransform.cxx
        modified:   ../tests/HT_test.cxx
parent dbace9d3
......@@ -176,7 +176,7 @@ namespace genfind {
double x = ahit.X();
double y = ahit.Y();
double R = x*x + y*y;
res.push_back(ROOT::Math::XYZTVector{ 2.0*x/R, -2.0*y/R, ahit.Z(), ahit.T() });
res.push_back(ROOT::Math::XYZTVector{ 2000.0*x/R, -2000.0*y/R, ahit.Z(), ahit.T() });
});
return res;
}
......@@ -277,7 +277,7 @@ namespace genfind {
//For better effect, make sure that there is a triple coincidence
//fThresh = 3; // this works well for regular, not conformal hits
fThresh = 3;
fThresh = 1;
for(int k_func = 0; k_func < n_funcs; k_func++) {
if (k_func==i_func || k_func==j_func){continue;}
auto f3 = funcs[k_func];
......@@ -401,7 +401,7 @@ namespace genfind {
// find all points within a fairly relaxed threshold around max bin
// Accumulate list of all the hits that make up this group
//fThresh = 15; // this works for regular, not conformal hits
fThresh = 15;
fThresh = 1;
std::vector< int > trackhitlist;
for(int i_func = 0; i_func < n_funcs; i_func++) {
auto f1 = funcs[i_func];
......
......@@ -33,7 +33,7 @@ void HT_test(
TH2F * hxz = new TH2F("hxz", "hxz", 100, -1000, 1000, 100, -1000, 1000);
TH2F * hyz = new TH2F("hyz", "hyz", 100, -1000, 1000, 100, -1000, 1000);
TH2F * hxy2 = new TH2F("hxy2", "hxy2", 100, -1010 ,1010, 100, -1010, 1010);
TH2F * huv = new TH2F("huv", "huv", 100, -.02 , .02, 100, -.02, .02);
TH2F * huv = new TH2F("huv", "huv", 100, -120 , 120, 100, -120, 120);
......@@ -95,7 +95,7 @@ void HT_test(
hxy2->Fill(x_extra, y_extra);*/
auto mg = ht->FillHoughTransforms(hits);
// auto mg = ht->FillHoughTransforms(hits);
// use track truth and the tracknumber in the findtracks result to test performance.
......@@ -109,6 +109,7 @@ void HT_test(
}
auto mg = ht->FillHoughTransforms(conf_hits);
std::vector<std::vector<std::tuple<int,double,double>>> FoundTracks = ht->FindTracks(conf_hits);
//std::vector<std::vector<std::tuple<int,double,double>>> FoundTracks = ht->FindTracks(hits);
std::cout << "FindTracks found " << FoundTracks.size() << " tracks " << endl;
......@@ -175,6 +176,8 @@ void HT_test(
mg->Draw("a");
c->cd(3);
huv->Draw();
c->cd(4);
ht->fIntersect->Draw();
c = new TCanvas();
c->Divide(2,2);
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment