checktubes.C

Go to the documentation of this file.

{
  
  TFile *file = new TFile("simparams.root");
  if (!gGeoManager) {
  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();
  }
 
}