Hough.C

Go to the documentation of this file.

int Hough(int i)
{
gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
        gSystem->Load("libriemann.C");

        TString MCFile = "Mvd_Test.root";

        PndFileNameCreator creator(MCFile.Data());
        TString DigiFile = creator.GetDigiFileName(false).c_str();
        TString RecoFile = creator.GetRecoFileName(false).c_str();
        TString TrackFFile = creator.GetTrackFindingFileName(false).c_str();

        TH2D h2("h2","h2", 200,0,TMath::Pi(), 200,-TMath::Pi(), TMath::Pi());

        TFile* fReco = new TFile(RecoFile.Data());


        TTree* t = (TTree*)(fReco->Get("pndsim"));

        t->StartViewer();

        TClonesArray* PixReco = new TClonesArray("PndSdsHit");
        TClonesArray* StripReco = new TClonesArray("PndSdsHit");

        t->SetBranchAddress("MVDHitsPixel", &PixReco);
        t->SetBranchAddress("MVDHitsStrip", &StripReco);

        t->GetEntry(i);

        for (int l = 0; l < PixReco->GetEntriesFast(); l++)
        {
                PndSdsHit* myPixHit = (PndSdsHit*)PixReco->At(l);
                TVector3 hit;
                myPixHit->Position(hit);
                h2.Fill(hit.Theta(), hit.Phi());
                std::cout << "Pixel " << l << " : ";
                //myPixHit->Print();
                std::cout << "Theta: " << hit.Theta() << " Phi: " << hit.Phi() << std::endl;

        }

        for (int l = 0; l < StripReco->GetEntriesFast(); l++)
        {
                PndSdsHit* myStripHit = (PndSdsHit*)StripReco->At(l);
                TVector3 striphit;
                myStripHit->Position(striphit);
                h2.Fill(striphit.Theta(), striphit.Phi());
                std::cout << "Strip " << l << " : ";
                //myStripHit->Print();
                std::cout << "Theta: " << striphit.Theta() << " Phi: " << striphit.Phi() << std::endl;
        }
        gStyle->SetPalette(1);
        h2.DrawCopy("colz");
  return 0;
}