DisplacedSolid.cxx

Go to the documentation of this file.

// $Id$
 
// -----------------------------------------------------------------------
// The RootGM package of the Virtual Geometry Model
// Copyright (C) 2007, Ivana Hrivnacova
// All rights reserved.
//
// For the licensing terms see vgm/LICENSE.
// Contact: ivana@ipno.in2p3.fr
// -----------------------------------------------------------------------
 
//
// Class DisplacedSolid
// ---------------------
// VGM implementation for Root displaced solid
//
// Author: Ivana Hrivnacova; IPN Orsay
 
#include "BaseVGM/common/utilities.h"
 
#include "RootGM/common/transform.h"
#include "RootGM/solids/Box.h"
#include "RootGM/solids/DisplacedSolid.h"
#include "RootGM/solids/SolidMap.h"
 
#include "TGeoBBox.h"
#include "TGeoBoolNode.h"
#include "TGeoCompositeShape.h"
#include "TGeoHalfSpace.h"
#include "TGeoMatrix.h"
#include "TString.h"
 
#include <cstdlib>
 
const std::string RootGM::DisplacedSolid::fgkNameExtension = "_displ";
 
//_____________________________________________________________________________
RootGM::DisplacedSolid::DisplacedSolid(
  const std::string& name, VGM::ISolid* solid, TGeoMatrix* displacement)
  : VGM::ISolid(),
    VGM::IDisplacedSolid(),
    BaseVGM::VDisplacedSolid(),
    fCompositeShape(0),
    fConstituentSolid(solid)
{
 
  // Get solid from the solid map
  TGeoShape* rootSolid = RootGM::SolidMap::Instance()->GetSolid(solid);
 
  // Define new name
  // std::string newName = name + fgkNameExtension;  // ??
  std::string newName = name;
 
  // Register TGeo matrix
  displacement->SetName(newName.data());
  displacement->RegisterYourself();
 
  TGeoBoolNode* boolNode = new TGeoUnion(rootSolid, rootSolid, displacement, 0);
 
  fCompositeShape = new TGeoCompositeShape(newName.data(), boolNode);
 
  RootGM::SolidMap::Instance()->AddSolid(this, fCompositeShape);
}
 
//_____________________________________________________________________________
RootGM::DisplacedSolid::DisplacedSolid(TGeoBBox* box)
  : VGM::ISolid(),
    VGM::IDisplacedSolid(),
    BaseVGM::VDisplacedSolid(),
    fCompositeShape(0),
    fConstituentSolid(0)
{
 
  // Give exception if box does not include offset
  const Double_t* origin = box->GetOrigin();
  if (!origin) {
    std::cerr << "    RootGM::DisplacedSolid::DisplacedSolid: " << std::endl;
    std::cerr << "    Cannot create displaced solid from a box without offset."
              << std::endl;
    std::cerr << "    (TGeoBBox name = " << box->GetName() << ")" << std::endl;
    std::cerr << "*** Error: Aborting execution  ***" << std::endl;
    exit(1);
  }
 
  // Define new name
  std::string newName = std::string(box->GetName()) + fgkNameExtension;
 
  // Create VGM box for constituent solid not registered in solid map
  fConstituentSolid = new RootGM::Box(box, false);
  // RootGM::SolidMap::Instance()->AddSolid(fConstituentSolid, box);
 
  // Create and register TGeo matrix
  TGeoMatrix* displacement =
    new TGeoTranslation(origin[0], origin[1], origin[2]);
  displacement->SetName(newName.data());
  displacement->RegisterYourself();
 
  TGeoBoolNode* boolNode = new TGeoUnion(box, box, displacement, 0);
 
  fCompositeShape = new TGeoCompositeShape(newName.data(), boolNode);
 
  RootGM::SolidMap::Instance()->AddSolidInRootMapOnly(this, fCompositeShape);
  RootGM::SolidMap::Instance()->AddSolidInVGMMapOnly(this, box);
}
 
//_____________________________________________________________________________
RootGM::DisplacedSolid::DisplacedSolid(TGeoHalfSpace* halfSpace)
  : VGM::ISolid(),
    VGM::IDisplacedSolid(),
    BaseVGM::VDisplacedSolid(),
    fCompositeShape(0),
    fConstituentSolid(0)
{
 
  // Define new name
  std::string newName = std::string(halfSpace->GetName()) + fgkNameExtension;
 
  // Create a big box for constituent solid not registered in solid map
  double size = 1.0e03; // 100 m
                        // when set 1.0e05 the test fails
  TGeoBBox* box = new TGeoBBox("bigBox", size, size, size);
  fConstituentSolid = new RootGM::Box(box, false);
  RootGM::SolidMap::Instance()->AddSolid(fConstituentSolid, halfSpace);
 
  // Get half space parameters
  Double_t* point = halfSpace->GetPoint();
  Double_t* norm = halfSpace->GetNorm();
 
  // Rotate box to get the norm of the z+ side (0,0,1)
  // in the norm direction
  Double_t mtx[3][3];
  Double_t from[3];
  from[0] = 0.;
  from[1] = 0.;
  from[2] = 1.0;
  RootGM::fromToRotation(from, norm, mtx);
 
  Double_t rot[9];
  rot[0] = mtx[0][0];
  rot[1] = mtx[0][1];
  rot[2] = mtx[0][2];
  rot[3] = mtx[1][0];
  rot[4] = mtx[1][1];
  rot[5] = mtx[1][2];
  rot[6] = mtx[2][0];
  rot[7] = mtx[2][1];
  rot[8] = mtx[2][2];
  TGeoRotation* rotation = new TGeoRotation();
  rotation->SetMatrix(rot);
 
  // Transform the box origin
  Double_t origin[3];
  origin[0] = 0;
  origin[1] = 0;
  origin[2] = -size;
  Double_t newOrigin[3];
  rotation->LocalToMaster(origin, newOrigin);
  for (Int_t i = 0; i < 3; i++) newOrigin[i] = newOrigin[i] + point[i];
 
  // Create and register TGeo matrix
  TGeoMatrix* displacement =
    new TGeoCombiTrans(newOrigin[0], newOrigin[1], newOrigin[2], rotation);
  displacement->SetName(newName.data());
  displacement->RegisterYourself();
 
  TGeoBoolNode* boolNode = new TGeoUnion(box, box, displacement, 0);
 
  fCompositeShape = new TGeoCompositeShape(newName.data(), boolNode);
 
  RootGM::SolidMap::Instance()->AddSolidInRootMapOnly(this, fCompositeShape);
  RootGM::SolidMap::Instance()->AddSolidInVGMMapOnly(this, halfSpace);
}
 
//_____________________________________________________________________________
RootGM::DisplacedSolid::DisplacedSolid()
  : VGM::ISolid(), VGM::IDisplacedSolid(), BaseVGM::VDisplacedSolid()
{
}
 
//_____________________________________________________________________________
RootGM::DisplacedSolid::DisplacedSolid(const DisplacedSolid& rhs)
  : VGM::ISolid(rhs), VGM::IDisplacedSolid(rhs), BaseVGM::VDisplacedSolid(rhs)
{
}
 
//_____________________________________________________________________________
RootGM::DisplacedSolid::~DisplacedSolid()
{
  //
}
 
//
// public methods
//
 
//_____________________________________________________________________________
std::string RootGM::DisplacedSolid::Name() const
{
  // Returns the Displaced solid name
  // ---
 
  return fCompositeShape->GetName();
}
 
//_____________________________________________________________________________
VGM::ISolid* RootGM::DisplacedSolid::ConstituentSolid() const
{
  // Returns the second constituent solid.
  // ---
 
  return fConstituentSolid;
}
 
//_____________________________________________________________________________
VGM::Transform RootGM::DisplacedSolid::Displacement() const
{
  // Returns the solid displacemnt transformation
  // in the frame of the first (left) solid.
  // ---
 
  TGeoBoolNode* boolNode = fCompositeShape->GetBoolNode();
 
  TGeoMatrix* matrixB = boolNode->GetLeftMatrix();
  TGeoHMatrix transformB(*matrixB);
 
  return Transform(transformB);
}