GeoNode.cxx

Go to the documentation of this file.

// $Id: GeoNode.cxx,v 1.8 2009/02/28 21:46:13 gmieg Exp $
//
// GeoNode
//
// GeoNode is a specialized TGeoNode to implement Minos specific methods
// that are generic to all Minos node types (plane,module,strip)
//
// Author:  S. Kasahara 06/04
#include <cassert>

#include <TGeoManager.h>
#include <TGeoVolume.h>
#include <TGeoMatrix.h>

#include "GeoGeometry/GeoNode.h"
#include "MessageService/MsgService.h"

ClassImp(GeoNode)

CVSID("$Id: GeoNode.cxx,v 1.8 2009/02/28 21:46:13 gmieg Exp $");

//_____________________________________________________________________________
GeoNode::GeoNode(GeoGeometry* geom, TGeoVolume* vol, TGeoMatrix* matrix, 
                 TGeoVolume* parVol,std::string globalpath,
                 std::string nodename) : 
                 TGeoNodeMatrix(vol,matrix), fGeoGeometry(geom),
                 fGlobalPath(globalpath),fGlobalMatrix(0) {
  
  // Normal constructor - Protected.

  UpdateGlobalManager();
  this -> SetName(nodename.c_str());
  this -> SetParent(parVol);
  
  matrix -> RegisterYourself();
  
}


//_____________________________________________________________________________
GeoNode::~GeoNode() {
  // Destructor should delete any owned objects

  UpdateGlobalManager();

  if ( CountRef() ) {
    MSG("Geo",Msg::kWarning)
      << "GeoNode destructor " << CountRef()
      << " outstanding references " << endl;
  }

  if ( fGlobalMatrix ) delete fGlobalMatrix; fGlobalMatrix = 0;

}

//_____________________________________________________________________________
void GeoNode::IncrementRef() {
  // Increment reference count to this node

  UpdateGlobalManager();

  fRef++; // increment own reference count
  fGeoGeometry->IncrementRef(); // and that of GeoGeometry
  SetVisibility(1);
  
}

//_____________________________________________________________________________
void GeoNode::DecrementRef() {
  // Decrement reference count to this node

  UpdateGlobalManager();

  fRef--; // decrement own reference count
  fGeoGeometry->DecrementRef(); // and that of GeoGeometry
  if (!CountRef() )SetVisibility(0);
}

//_____________________________________________________________________________
void GeoNode::SetParent(TGeoVolume* parentvol) {
  // Set the node in its parent volume and update the parent list of
  // child nodes
  // Private method used during geometry build.

  UpdateGlobalManager();

  if ( !parentvol ) {
    MSG("Geo",Msg::kWarning) << "Null parent volume arg" 
                             << "\n No action taken" << endl;
    return;
  }

  this -> SetMotherVolume(parentvol);
  TObjArray* listofnodes = parentvol->GetNodes();
   
  if (!listofnodes) { 
    listofnodes = new TObjArray();
    parentvol -> SetNodes(listofnodes); // SetBit(kVolumeImportNodes) called 
    ((TObject*)parentvol) -> 
              ResetBit(TGeoVolume::kVolumeImportNodes); // mark to be deleted
  }
  
  listofnodes -> Add(this);
  this -> SetNumber(1);
  
  return;

}

//_____________________________________________________________________________
void GeoNode::UpdateGlobalMatrix() {
  // Stores a copy of the matrix to perform local-to-global(MARS) 
  // and global-to-local transformations with this node.  Protected method
  // called by GeoGeometry after geometry is closed.

  UpdateGlobalManager();

  if ( fGlobalMatrix ) {
    // In principal this should only be called once at geometry closure.
    MSG("Geo",Msg::kWarning) << "UpdateGlobalMatrix deleting previous matrix\n"
                             << " for node " << this->GetName() << "." << endl;
    delete fGlobalMatrix;
    fGlobalMatrix = 0;
  }
  
  if ( fGlobalPath.empty() ) {
    MSG("Geo",Msg::kFatal) 
            << "UpdateGlobalMatrix called with empty path. Abort." << endl;
    abort();
  }

  if ( !gGeoManager->cd(fGlobalPath.c_str()) ) {
    MSG("Geo",Msg::kFatal)
      << "UpdateGlobalMatrix unable to cd to global path\n" 
      << fGlobalPath.c_str() << ". Abort." << endl;
    abort();
  }
  
  TGeoHMatrix* currentmatrix = gGeoManager->GetCurrentMatrix();
  if ( currentmatrix ) {
    // Take a copy to store with node. Store as TGeoCombiTrans to
    // avoid storing additional bytes associated with TGeoHMatrix::fScale 
    // since this is not used.
    fGlobalMatrix = new TGeoCombiTrans(*currentmatrix);
  }
  else {
    // should never happen
    MSG("Geo",Msg::kFatal) 
    << "UpdateGlobalMatrix retrieved null matrix ptr from TGeoManager." 
    << " Abort." << endl;
    abort();
  }

  return;
  
}

//_____________________________________________________________________________
TVector3 GeoNode::LocalToGlobal(const TVector3& local) const {
  // Convert from local to global (MARS) coordinates

  UpdateGlobalManager();
  
  Double_t lxyz[3] = {local.X(),local.Y(),local.Z()};
  Double_t gxyz[3] = {0};
  
  LocalToGlobal(lxyz,gxyz);
  
  return TVector3(gxyz[0],gxyz[1],gxyz[2]);
  
}

//_____________________________________________________________________________
TVector3 GeoNode::GlobalToLocal(const TVector3& global) const {
  // Convert from local to global (MARS) coordinates

  UpdateGlobalManager();
  
  Double_t gxyz[3] = {global.X(),global.Y(),global.Z()};
  Double_t lxyz[3] = {0};
  
  GlobalToLocal(gxyz,lxyz);
  
  return TVector3(lxyz[0],lxyz[1],lxyz[2]);
  
}

//_____________________________________________________________________________
void GeoNode::LocalToGlobal(Double_t* lxyz, Double_t* gxyz) const {
  
  UpdateGlobalManager();

  if ( !fGlobalMatrix ) {
    MSG("Geo",Msg::kFatal) << "LocalToGlobal for node "
                          << this->GetName() << " has null fGlobalMatrix ptr." 
                          << endl;
    abort();
  }
  
  fGlobalMatrix -> LocalToMaster(lxyz,gxyz);
  
  
}

//_____________________________________________________________________________
void GeoNode::GlobalToLocal(Double_t* gxyz, Double_t* lxyz) const {
  
  UpdateGlobalManager();

  if ( !fGlobalMatrix ) {
    MSG("Geo",Msg::kFatal) << "GlobalToLocal for node "
                          << this->GetName() << " has null fGlobalMatrix ptr." 
                          << endl;
    abort();
  }
  
  fGlobalMatrix -> MasterToLocal(gxyz,lxyz);
   
}

//_____________________________________________________________________________
TVector3 GeoNode::LocalToGlobalVect(const TVector3& local) const {
  // Convert from local to global (MARS) direction vector

  UpdateGlobalManager();
  
  Double_t lxyz[3] = {local.X(),local.Y(),local.Z()};
  Double_t gxyz[3] = {0};
  
  LocalToGlobalVect(lxyz,gxyz);
  
  return TVector3(gxyz[0],gxyz[1],gxyz[2]);
  
}

//_____________________________________________________________________________
TVector3 GeoNode::GlobalToLocalVect(const TVector3& global) const {
  // Convert from local to global (MARS) direction vector

  UpdateGlobalManager();
  
  Double_t gxyz[3] = {global.X(),global.Y(),global.Z()};
  Double_t lxyz[3] = {0};
  
  GlobalToLocalVect(gxyz,lxyz);
  
  return TVector3(lxyz[0],lxyz[1],lxyz[2]);
  
}

//_____________________________________________________________________________
void GeoNode::LocalToGlobalVect(Double_t* lxyz, Double_t* gxyz) const {
  
  UpdateGlobalManager();

  if ( !fGlobalMatrix ) {
    MSG("Geo",Msg::kFatal) << "LocalToGlobalVect for node "
                          << this->GetName() << " has null fGlobalMatrix ptr." 
                          << endl;
    abort();
  }
  
  fGlobalMatrix -> LocalToMasterVect(lxyz,gxyz);
  
}

//_____________________________________________________________________________
void GeoNode::GlobalToLocalVect(Double_t* gxyz, Double_t* lxyz) const {
  
  UpdateGlobalManager();

  if ( !fGlobalMatrix ) {
    MSG("Geo",Msg::kFatal) << "GlobalToLocalVect for node "
                          << this->GetName() << " has null fGlobalMatrix ptr." 
                          << endl;
    abort();
  }
  
  fGlobalMatrix -> MasterToLocalVect(gxyz,lxyz);
   
}

//_____________________________________________________________________________

---