UgliSteelPlnHandle.cxx

Go to the documentation of this file.

// $Id: UgliSteelPlnHandle.cxx,v 1.12 2006/06/27 20:08:23 rhatcher Exp $
//
// UgliSteelPlnHandle
//
// UgliSteelPlnHandle is handle to a single steel plane
//
// Author:  R. Hatcher 2000.11.15
//

#include "UgliGeometry/UgliSteelPlnHandle.h"
#include "UgliGeometry/Ugli.h"

#include "MessageService/MsgService.h"
CVSID("$Id: UgliSteelPlnHandle.cxx,v 1.12 2006/06/27 20:08:23 rhatcher Exp $");

#include <cassert>

ClassImp(UgliSteelPlnHandle)

//_____________________________________________________________________________
UgliSteelPlnHandle::UgliSteelPlnHandle()
  : UgliPlnHandle((UgliPlnNode*)0), fUgliSteelPlnNode(0), fGeoSteelPlnNode(0)
{
   // Default constructor

   MSG("Ugli",Msg::kDebug) 
     << "UgliSteelPlnHandle created by default ctor" << endl;
//   assert(0);
}

//_____________________________________________________________________________
UgliSteelPlnHandle::~UgliSteelPlnHandle()
{
   // handle deletion updates reference count
   // done by super-class UgliPlnHandle
   // don't do it here or we'll get two counts for every handle

   // if (fUgliSteelPlnNode) fUgliSteelPlnNode->DecrementRef();
   // if (fGeoSteelPlnNode) fGeoSteelPlnNode->DecrementRef();
}

//_____________________________________________________________________________
UgliSteelPlnHandle::UgliSteelPlnHandle(UgliSteelPlnNode* node)
  : UgliPlnHandle(node), fUgliSteelPlnNode(node), fGeoSteelPlnNode(0)
{
   // ctor keeps reference count up-to-date
   // done by super-class UgliPlnHandle
   // don't do it here or we'll get two counts for every handle

   // if (fUgliSteelPlnNode) fUgliSteelPlnNode->IncrementRef();
   // if (fGeoSteelPlnNode) fGeoSteelPlnNode->IncrementRef();
}

//_____________________________________________________________________________
UgliSteelPlnHandle::UgliSteelPlnHandle(GeoSteelPlnNode* node)
  : UgliPlnHandle(node), fUgliSteelPlnNode(0), fGeoSteelPlnNode(node)
{
   // ctor keeps reference count up-to-date
   // done by super-class UgliPlnHandle
   // don't do it here or we'll get two counts for every handle

   // if (fUgliSteelPlnNode) fUgliSteelPlnNode->IncrementRef();
   // if (fGeoSteelPlnNode) fGeoSteelPlnNode->IncrementRef();
}

//_____________________________________________________________________________
UgliSteelPlnHandle::UgliSteelPlnHandle(const UgliSteelPlnHandle &that)
   : UgliPlnHandle(that),
     UgliSteelPlnABC(),
     fUgliSteelPlnNode(that.fUgliSteelPlnNode),
     fGeoSteelPlnNode(that.fGeoSteelPlnNode)
{
   // copy constructor keeps reference counts up-to-date
   // done by super-class UgliPlnHandle
   // don't do it here or we'll get two counts for every handle

   // if (fUgliSteelPlnNode) fUgliSteelPlnNode->IncrementRef();
   // if (fGeoSteelPlnNode) fGeoSteelPlnNode->IncrementRef();
}

//_____________________________________________________________________________
UgliSteelPlnHandle&  UgliSteelPlnHandle::operator=(const UgliSteelPlnHandle& that)
{
   // assignment keeps reference counts up-to-date
   // done by super-class UgliPlnHandle
   // don't do it here or we'll get two counts for every handle
   
   // if (fUgliSteelPlnNode) fUgliSteelPlnNode->DecrementRef();
   // fUgliSteelPlnNode = that.fUgliSteelPlnNode;
   // if (fUgliSteelPlnNode) fUgliSteelPlnNode->IncrementRef();

   if ( this != &that ) {
     UgliPlnHandle::operator=(that);
     fUgliSteelPlnNode = that.fUgliSteelPlnNode;
     //fUgliPlnNode      = that.fUgliPlnNode; // redundant? s.k.3/16/05
     fGeoSteelPlnNode = that.fGeoSteelPlnNode;
   }
   return *this;
}

//_____________________________________________________________________________
TVector3 UgliSteelPlnHandle::GlobalToLocal(const TVector3& global,
                                           const Bool_t globalInXYZ) const
{
  // convert a global position into a local (volume) one
  // globalInXYZ determines whether result is in XYZ or UVZ
  TVector3 globalXYZ = (globalInXYZ) ? 
      global : Ugli::uvz2xyz(GetPlexPlaneId().GetDetector(),global);

  if (!IsGeo())
      return fUgliSteelPlnNode->GlobalToLocal(global);
  else {
      TVector3 scaledglobal(global); scaledglobal *= 1./Munits::cm;
      return (fGeoSteelPlnNode->GlobalToLocal(scaledglobal)*Munits::cm); 
  }
}

//_____________________________________________________________________________
TVector3 UgliSteelPlnHandle::LocalToGlobal(const TVector3& local,
                                           const Bool_t globalInXYZ) const
{
  // convert a local (volume) position into a global one
  // globalInXYZ determines whether result is in XYZ or UVZ
  TVector3 globalXYZ;
  if (!IsGeo())
      globalXYZ = fUgliSteelPlnNode->LocalToGlobal(local); 
  else {
      TVector3 scaledlocal(local); scaledlocal *= 1./Munits::cm;
      globalXYZ = (fGeoSteelPlnNode->LocalToGlobal(scaledlocal)*Munits::cm); 
  }

  if (globalInXYZ ) return globalXYZ;
  else              return Ugli::xyz2uvz(GetPlexPlaneId().GetDetector(),
                                         globalXYZ);
}

//_____________________________________________________________________________
TVector3 UgliSteelPlnHandle::GlobalToLocalVect(const TVector3& global,
                                               const Bool_t globalInXYZ) const
{
  // convert a global direction vector into a local (volume) one
  // globalInXYZ determines whether result is in XYZ or UVZ
  TVector3 globalXYZ = (globalInXYZ) ? 
      global : Ugli::uvz2xyz(GetPlexPlaneId().GetDetector(),global);

  if (!IsGeo())
      return fUgliSteelPlnNode->GlobalToLocalVect(global);
  else {
      TVector3 scaledglobal(global); scaledglobal *= 1./Munits::cm;
      return (fGeoSteelPlnNode->GlobalToLocalVect(scaledglobal)*Munits::cm); 
  }
}

//_____________________________________________________________________________
TVector3 UgliSteelPlnHandle::LocalToGlobalVect(const TVector3& local,
                                               const Bool_t globalInXYZ) const
{
  // convert a local (volume) direction into a global one
  // globalInXYZ determines whether result is in XYZ or UVZ
  TVector3 globalXYZ;
  if (!IsGeo())
      globalXYZ = fUgliSteelPlnNode->LocalToGlobalVect(local); 
  else {
      TVector3 scaledlocal(local); scaledlocal *= 1./Munits::cm;
      globalXYZ = (fGeoSteelPlnNode->LocalToGlobalVect(scaledlocal)*Munits::cm); 
  }

  if (globalInXYZ ) return globalXYZ;
  else              return Ugli::xyz2uvz(GetPlexPlaneId().GetDetector(),
                                         globalXYZ);
}

//_____________________________________________________________________________

---