Functions
timer	Start ()

TFile	filerun ("testrun.root")

treepnt	SetBranchAddress ("STTPoint",&pnt)

TFile	filedigi ("testdigi.root")

treedigi	SetBranchAddress ("STTHit",&digi)

TFile	filehelix ("testreco.root")

treereco	SetBranchAddress ("SttHelixHit",&hh)

treereco	SetBranchAddress ("STTTrack",&track)

c	Divide (1, 2)

treepnt	GetEntry (evt)

c	cd (1)

hxy	Draw ()

c	cd (2)

	if (!stttrack) continue

cout<< "evt"<< evt<< " track " << k<< " num hits "<< stttrack->	GetNofHelixHits ()

Variables
TStopwatch	timer

TTree *	treepnt = (TTree*) filerun.Get("pndsim")

TClonesArray *	pnt = new TClonesArray("PndSttPoint")

TTree *	treedigi = (TTree*) filedigi.Get("pndsim")

TClonesArray *	digi = new TClonesArray("PndSttHit")

TTree *	treereco = (TTree*) filehelix.Get("pndsim")

TClonesArray *	hh = new TClonesArray("PndSttHelixHit")

TClonesArray *	track = new TClonesArray("PndSttTrack")

TCanvas *	c = new TCanvas("c", "c", 0, 0, 600, 600)

TH2F *	hxy = new TH2F("hxy","hxy",100,-42,42, 100,-42,42)

TH2F *	hyz = new TH2F("hyz","hyz",100,0,50, 100,-40,110)

int	evt = 0

cout<< "blue = Monte Carlo " << endl;cout<< "red = Helix Hit "<< endl;cout<< "green = Center Of Tubes "<< endl;for(Int_t k=0;k< track->GetEntriesFast();k++){PndSttTrack stttrack=(PndSttTrack ) track->	At (k)

Double_t	d0 = stttrack->GetDist()

Double_t	phi0 = stttrack->GetPhi()

Double_t	R = stttrack->GetRad()

Double_t	z0 = stttrack->GetZ()

Double_t	tanl = stttrack->GetTanL()

Double_t	h = -(Int_t) stttrack->GetCharge()

Double_t	ptran = 0.003 * 2 * R

Double_t	plong = ptran * tanl

Double_t	ptot = sqrt(plongplong + ptranptran)

Double_t	x0 = (d0 + R) * cos(phi0)

Double_t	y0 = (d0 + R) * sin(phi0)

Function Documentation

c cd ( 1 )

c cd ( 2 )

c Divide	(	1	,
		2
	)

hxy Draw ( )

Definition at line 64 of file NHitsPerEvent.C.

            {
 
         std::cout << "mvdHitArray: " << mvdHitArray->GetEntries() << " projectedHitArray " << mvdProjectedArray->GetEntries() << std::endl;
 
 //      gHitPoints->Clear();
         ExtractData(mvdHitArray, gHitPoints);
         ExtractData(mvdProjectedArray, gProjections);
         //if (gHitPoints->GetN() != x.size()) gHitPoints->Expand(x.size());
 
         std::cout << "NHits: " << gHitPoints->GetN() << std::endl;
         gHitPoints->Draw("apl");
         gProjections->Draw("pl");
 }

TFile filedigi ( "testdigi.root" )

TFile filehelix ( "testreco.root" )

TFile filerun ( "testrun.root" )

treepnt GetEntry ( evt )

cout<< "evt" << evt << " track " << k << " num hits " << stttrack-> GetNofHelixHits ( )

TMarker *mrkpnt = new TMarker(mcpoint->GetYtot(), mcpoint->GetZtot(), 2); mrkpnt->SetMarkerColor(4); mrkpnt->Draw("SAME"); TMarker *mrkdigi = new TMarker(hit->GetY(), hit->GetZ(), 6); mrkdigi->SetMarkerColor(3); mrkdigi->Draw("SAME"); TMarker *mrkhh = new TMarker(helixhit->GetY(), helixhit->GetZ(), 5); mrkhh->SetMarkerColor(2); mrkhh->Draw("SAME");

Definition at line 73 of file checkhelixhit.C.

                                                              {
       Int_t iHHit = stttrack->GetHelixHitIndex(j);
       PndSttHelixHit *helixhit = (PndSttHelixHit*) hh->At(iHHit);
 
       Int_t hitindex = helixhit->GetHitIndex();
       PndSttHit* hit = (PndSttHit*) digi->At(hitindex);
       PndSttPoint *mcpoint = (PndSttPoint*) pnt->At(hit->GetRefIndex());
 
       c->cd(1);
 
       // points/hits
       TMarker *mrkpnt = new TMarker(mcpoint->GetXtot(), mcpoint->GetYtot(), 2);
       mrkpnt->SetMarkerColor(4);
       mrkpnt->Draw("SAME");
       TMarker *mrkdigi = new TMarker(hit->GetX(), hit->GetY(), 6);
       mrkdigi->SetMarkerColor(3);
       mrkdigi->Draw("SAME");
       TMarker *mrkhh = new TMarker(helixhit->GetX(), helixhit->GetY(), 5);
       mrkhh->SetMarkerColor(2);
       mrkhh->Draw("SAME");
 
       c->cd(2);
 
       //    cout << "MC " <<  mcpoint->GetXtot()
       //           << " " <<             mcpoint->GetYtot()
       //           << " " <<             mcpoint->GetZtot() << endl;
 
       //   cout << "helixhit " <<  helixhit->GetX()
       //        << " " <<             helixhit->GetY()
       //        << " " <<             helixhit->GetZ() << endl;
 
 
 
       Double_t scoslMC = stttrack->CalculateScosl(mcpoint->GetXtot(), mcpoint->GetYtot());
       TMarker *mrkpnt = new TMarker(scoslMC, mcpoint->GetZtot(), 2);
       mrkpnt->SetMarkerColor(4);
       mrkpnt->Draw("SAME");
 
       Double_t scoslH  = stttrack->CalculateScosl(hit->GetX(), hit->GetY());
       TMarker *mrkdigi = new TMarker(scoslH, hit->GetZ(), 6);
       mrkdigi->SetMarkerColor(3);
       //  mrkdigi->Draw("SAME");
 
       Double_t scoslHH = stttrack->CalculateScosl(helixhit->GetX(), helixhit->GetY());
       TMarker *mrkhh = new TMarker(scoslHH, helixhit->GetZ(),  5);
       mrkhh->SetMarkerColor(2);
       mrkhh->Draw("SAME");
 
       //      cout <<  scoslMC<< " " << mcpoint->GetZtot() << endl;
     }

if ( ! stttrack )

Definition at line 4 of file checktubes.C.

                     {
   file->Get("FairBaseParSet"); 
   TGeoManager *geoMan = gGeoManager; 
   if(!geoMan) { 
     file->Get("FairGeoParSet"); 
     geoMan = gGeoManager; 
     if(!geoMan) { 
       cout << "Could not find valid GeoManager. Abort now!" << endl; 
       exit(1); 
     } 
   }
 
   PndGeoSttPar *parameters = file.Get("PndGeoSttPar");
   PndSttMapCreator *mapper = new PndSttMapCreator(parameters);
   TClonesArray *fTubeArray = mapper->FillTubeArray();
 
   TCanvas *c = new TCanvas("c", "", 400, 0, 900, 900);
   TH2F *h2 = new TH2F("h2", "", 100, -42, 42, 100, -42, 42);
   h2->SetStats(kFALSE);
   h2->Draw();
 
   TArc *arc = NULL;
   TMarker *mrk = NULL;
   PndSttTube *tube = NULL;
   int itube = 0;
   while(itube != -1) {
     c->Clear();
     c->Draw();
     h2->Draw();
 
 
     // draw all the tubes
     for(itube = 1; itube < fTubeArray->GetEntriesFast(); itube++) {
       tube = (PndSttTube*) fTubeArray->At(itube);
       if(tube->IsParallel()) {
         arc = new TArc(tube->GetPosition().X(), tube->GetPosition().Y(), 0.5);
         arc->Draw("SAME");
       }
       else {
         mrk = new TMarker(tube->GetPosition().X(), tube->GetPosition().Y(), 6);
         mrk->Draw("SAME");
       }
     }
     // ............................
 
     
     
     cout << "WHICH tube do you want to test? [to exit type -1]" << endl;
     cin >> itube;
     if(itube == -1) exit(0);
     cout << "TESTING tube " << itube << endl;
     tube = (PndSttTube*) fTubeArray.At(itube);
     if(!tube) {
       cout << "this tube does not exist" << endl; 
       continue; 
     }
 
     TArrayI neighborings = tube->GetNeighborings();
     cout << "nof neighboring tubes: " << neighborings->GetSize() << endl;
     arc = new TArc(tube->GetPosition().X(), tube->GetPosition().Y(), 0.5);
     arc->SetFillColor(kYellow);
     arc->Draw("SAME");
 
     for(int inei = 0; inei < neighborings.GetSize(); inei++) { 
       tube = (PndSttTube*) fTubeArray.At(neighborings.At(inei));
       TVector3 position = tube->GetPosition(); 
       cout << " " << neighborings.At(inei);
       if(tube->IsParallel()) {
         arc = new TArc(position.X(), position.Y(), 0.5); 
         arc->SetFillColor(3);
         arc->Draw("SAME");
       }
       else {
         mrk = new TMarker(position.X(), position.Y(), 6);
 
         TVector3 extr1 = position - tube->GetWireDirection() * tube->GetHalfLength();
         TVector3 extr2 = position + tube->GetWireDirection() * tube->GetHalfLength();
         //      mrk = new TMarker(extr1.X(), extr1.Y(), 6);
         //      mrk = new TMarker(extr2.X(), extr2.Y(), 6);
 
         mrk->SetMarkerColor(3);
         mrk->Draw("SAME");
       }
     }
 
     cout << endl;
     c->Update();
     c->Modified();
   }
  
 }

treepnt SetBranchAddress	(	"STTPoint"	,
		&	pnt
	)

treedigi SetBranchAddress	(	"STTHit"	,
		&	digi
	)

treereco SetBranchAddress	(	"SttHelixHit"	,
		&	hh
	)

treereco SetBranchAddress	(	"STTTrack"	,
		&	track
	)

timer Start ( )

Variable Documentation

cout<< "blue = Monte Carlo " << endl; cout << "red = Helix Hit " << endl; cout << "green = Center Of Tubes " << endl; for (Int_t k = 0; k < track->GetEntriesFast(); k++) { PndSttTrack *stttrack = (PndSttTrack*) track-> At(k)

Definition at line 56 of file checkhelixhit.C.

TCanvas* c = new TCanvas("c", "c", 0, 0, 600, 600)

Definition at line 31 of file checkhelixhit.C.