Commit d92420fb authored by Johnston's avatar Johnston
Browse files

modified: hough_transform2.cxx

        purpose to isolate linefinding from helix finding in conformal
                            space
parent cf17443c
......@@ -85,7 +85,8 @@ void hough_transform2(
return;
}
for(auto ahit : all_chits) {
//for(auto ahit : all_chits) {
for(auto ahit : hits) {
//std::cout << ahit.X() << " , " << ahit.Y() << std::endl;
huv->Fill( ahit.X(), ahit.Y() );
hphi->Fill(TMath::ATan(ahit.Y()/ahit.X()));
......@@ -94,9 +95,13 @@ void hough_transform2(
hrphi->Fill(theta/degree, ahit.X()*TMath::Cos(theta) + ahit.Y()*TMath::Sin(theta) );
}
}
auto mg = ht->FillHoughTransforms(all_chits);
//auto mg = ht->FillHoughTransforms(all_chits);
std::cout << " n peaks : " << ht->FindPeaks(all_chits) << std::endl;
auto mg = ht->FillHoughTransforms(hits);
std::cout << " n peaks : " << ht->FindPeaks(hits) << std::endl;
//std::cout << " n peaks : " << ht->FindPeaks(all_chits) << std::endl;
auto pretrack_hits = ht->GetPreTrackHits(master_hits,4.0);
......@@ -166,6 +171,16 @@ void hough_transform2(
hphi->Draw();
//hrphi->Draw("lego2");
c = new TCanvas();
c->Divide(2,2);
int counter = 1;
for( auto ahist : alltrack ){
c->cd(counter);
ahist->Draw("box");
counter++;
if( counter > 4) {break;}
}
c = new TCanvas();
c->Divide(2,2);
c->cd(1);
......
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