parallel_world.C

Go to the documentation of this file.

/// \file
/// \ingroup tutorial_geom
/// Misaligning geometry generate in many cases overlaps, due to the idealization
/// of the design and the fact that in real life movements of the geometry volumes
/// have constraints and are correlated.
///
/// This typically generates inconsistent
/// response of the navigation methods, leading to inefficiencies during tracking,
/// errors in the material budget calculations, and so on. Among those, there are
/// dangerous cases when the hidden volumes are sensitive.
/// This macro demonstrates how to use the "parallel world" feature to assign
/// highest navigation priority to some physical paths in geometry.
///
/// \macro_image
/// \macro_code
///
/// \author Andrei Gheata
 
void align();
 
//______________________________________________________________________________
void parallel_world(Bool_t usepw=kTRUE, Bool_t useovlp=kTRUE)
{
   // web geometry display does not support "parallel world" feature
   gROOT->SetWebDisplay("off");
 
   TGeoManager *geom = new TGeoManager("parallel_world", "Showcase for prioritized physical paths");
   TGeoMaterial *matV = new TGeoMaterial("Vac", 0,0,0);
   TGeoMedium *medV = new TGeoMedium("MEDVAC",1,matV);
   TGeoMaterial *matAl = new TGeoMaterial("Al", 26.98,13,2.7);
   TGeoMedium *medAl = new TGeoMedium("MEDAL",2,matAl);
   TGeoMaterial *matSi = new TGeoMaterial("Si", 28.085,14,2.329);
   TGeoMedium *medSi = new TGeoMedium("MEDSI",3,matSi);
   TGeoVolume *top = gGeoManager->MakeBox("TOP",medV,100,400,1000);
   gGeoManager->SetTopVolume(top);
 
   // Shape for the support block
   TGeoBBox *sblock = new TGeoBBox("sblock", 20,10,2);
   // The volume for the support
   TGeoVolume *support = new TGeoVolume("block",sblock, medAl);
   support->SetLineColor(kGreen);
 
   // Shape for the sensor to be prioritized in case of overlap
   TGeoBBox *ssensor = new TGeoBBox("sensor", 19,9,0.2);
   // The volume for the sensor
   TGeoVolume *sensor = new TGeoVolume("sensor",ssensor, medSi);
   sensor->SetLineColor(kRed);
 
   // Chip assembly of support+sensor
   TGeoVolumeAssembly *chip = new TGeoVolumeAssembly("chip");
   chip->AddNode(support, 1);
   chip->AddNode(sensor,1, new TGeoTranslation(0,0,-2.1));
 
   // A ladder that normally sags
   TGeoBBox *sladder = new TGeoBBox("sladder", 20,300,5);
   // The volume for the ladder
   TGeoVolume *ladder = new TGeoVolume("ladder",sladder, medAl);
   ladder->SetLineColor(kBlue);
 
   // Add nodes
   top->AddNode(ladder,1);
   for (Int_t i=0; i<10; i++)
      top->AddNode(chip, i+1, new TGeoTranslation(0, -225.+50.*i, 10));
 
   gGeoManager->CloseGeometry();
   TGeoParallelWorld *pw = nullptr;
   if (usepw) pw = gGeoManager->CreateParallelWorld("priority_sensors");
// Align chips
   align();
   if (usepw) {
      if (useovlp) pw->AddOverlap(ladder);
      pw->CloseGeometry();
      gGeoManager->SetUseParallelWorldNav(kTRUE);
   }
   TString cname;
   cname = usepw ? "cpw" : "cnopw";
   TCanvas *c = (TCanvas*)gROOT->GetListOfCanvases()->FindObject(cname);
   if (c) c->cd();
   else   c = new TCanvas(cname, "",800,600);
   top->Draw();
//   top->RandomRays(0,0,0,0,sensor->GetName());
   // Track random "particles" coming from the block side and draw only the tracklets
   // actually crossing one of the sensors. Note that some of the tracks coming
   // from the outer side may see the full sensor, while the others only part of it.
   TStopwatch timer;
   timer.Start();
   top->RandomRays(100000,0,0,-30,sensor->GetName());
   timer.Stop();
   timer.Print();
   TView3D *view = (TView3D*)gPad->GetView();
   if (view) {
      view->SetParallel();
      view->Side();
   }
   if (usepw) pw->PrintDetectedOverlaps();
}
 
//______________________________________________________________________________
void align()
{
// Aligning 2 sensors so they will overlap with the support. One sensor is positioned
// normally while the other using the shared matrix
   TGeoPhysicalNode *node;
   TGeoParallelWorld *pw = gGeoManager->GetParallelWorld();
   Double_t sag;
   for (Int_t i=0; i<10; i++) {
      node = gGeoManager->MakePhysicalNode(TString::Format("/TOP_1/chip_%d",i+1));
      sag = 8.-0.494*(i-4.5)*(i-4.5);
      TGeoTranslation *tr = new TGeoTranslation(0., -225.+50.*i, 10-sag);
      node->Align(tr);
      if (pw) pw->AddNode(TString::Format("/TOP_1/chip_%d",i+1));
   }
}