UgliPlnNode.cxx

Go to the documentation of this file.

// $Id: UgliPlnNode.cxx,v 1.4 2003/10/29 19:28:02 rhatcher Exp $
//
// UgliPlnNode
//
// UgliPlnNode is a plane node
//
// Author:  R. Hatcher 2000.11.15
//

#include "UgliGeometry/UgliPlnNode.h"

#include "UgliGeometry/UgliGeometry.h"
#include "UgliGeometry/TGeometryX.h"

#include "TBRIK.h"
#include "TXTRU.h"

#include "MessageService/MsgService.h"
CVSID("$Id: UgliPlnNode.cxx,v 1.4 2003/10/29 19:28:02 rhatcher Exp $");

#include <cassert>

ClassImp(UgliPlnNode)

//_____________________________________________________________________________
UgliPlnNode::UgliPlnNode()
  : fUgliGeometry(0)
{
   // Default constructor (for i/o)
}

//_____________________________________________________________________________
UgliPlnNode::UgliPlnNode(UgliGeometry *ugliGeometry,
                         const PlexPlaneId& planeid)
   : TNodeX(planeid.AsString("p"),planeid.AsString("p"),"noshape"),
     fUgliGeometry(ugliGeometry), fPlaneId(planeid)
{

}

//_____________________________________________________________________________
UgliPlnNode::~UgliPlnNode()
{
   // destructor should delete any owned objects
   if (CountRef()) {
      MSG("Ugli",Msg::kWarning)
        << "~UgliPlnNode " << GetPlexPlaneId()
        << " still had " << CountRef()
        << " outstanding references " << endl;
   }
}

//_____________________________________________________________________________

void UgliPlnNode::IncrementRef()
{
   fRef++; 
   fUgliGeometry->IncrementRef();
   SetVisibility(1);
}

void UgliPlnNode::DecrementRef()
{ 
   fRef--; 
   fUgliGeometry->DecrementRef();
   if (!fRef) SetVisibility(0);
}

//_____________________________________________________________________________
Float_t UgliPlnNode::GetHalfThickness() const
{

   TBRIK *brik = dynamic_cast<TBRIK*>(fShape);
   if (brik) return brik->GetDz();

   TXTRU *xtru = dynamic_cast<TXTRU*>(fShape);
   if (xtru) {
      Int_t    nz = xtru->GetNz();
      Float_t* fz = xtru->GetZ();
      return 0.5*TMath::Abs(fz[0]-fz[nz-1]);
   }

   MSG("Ugli",Msg::kError) <<
      "UgliPlnNode::GetHalfThickness: not a BRIK or TXTRU" << endl;

   return 0;
}

//_____________________________________________________________________________
void UgliPlnNode::TransformLocal2Global(Double_t *lxyz, Double_t *gxyz) const
{
   // deal with fact that TNode::Master2Local isn't declared const
   UgliPlnNode* self = const_cast<UgliPlnNode*>(this);
   self->cd();
   self->UpdateMatrix();
   self->Local2Master(lxyz,gxyz);
}

//_____________________________________________________________________________
TVector3 UgliPlnNode::GetCenter() const
{
   Double_t lxyz[3] = {0., 0., 0.};
   Double_t gxyz[3];

   TransformLocal2Global(lxyz,gxyz);
   return TVector3(gxyz[0],gxyz[1],gxyz[2]);
}

//_____________________________________________________________________________
Float_t UgliPlnNode::GetX0() const
{
   Double_t lxyz[3] = {0., 0., 0.};
   Double_t gxyz[3];

   TransformLocal2Global(lxyz,gxyz);
   return gxyz[0];
}

//_____________________________________________________________________________
Float_t UgliPlnNode::GetY0() const
{
   Double_t lxyz[3] = {0., 0., 0.};
   Double_t gxyz[3];

   TransformLocal2Global(lxyz,gxyz);
   return gxyz[1];
}

//_____________________________________________________________________________
Float_t UgliPlnNode::GetZ0() const
{
   Double_t lxyz[3] = {0., 0., 0.};
   Double_t gxyz[3];

   TransformLocal2Global(lxyz,gxyz);
   return gxyz[2];
}

//_____________________________________________________________________________

---