UgliStripHandle.cxx

Go to the documentation of this file.

// $Id: UgliStripHandle.cxx,v 1.10 2005/09/01 08:07:08 rhatcher Exp $
//
// UgliStripHandle
//
// UgliStripHandle is handle to a single strip
//
// Author:  R. Hatcher 2000.11.15
//

#include "GeoGeometry/GeoScintMdlNode.h"
#include "UgliGeometry/UgliStripHandle.h"
#include "UgliGeometry/UgliStripNode.h"
//#include "UgliGeometry/UgliGeometry.h"

#include "MessageService/MsgService.h"
CVSID("$Id: UgliStripHandle.cxx,v 1.10 2005/09/01 08:07:08 rhatcher Exp $");

#include <cassert>

ClassImp(UgliStripHandle)

//_____________________________________________________________________________
UgliStripHandle::UgliStripHandle()
  : fUgliStripNode(0),fGeoStripNode(0)
{
   // Default constructor

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

//_____________________________________________________________________________
UgliStripHandle::~UgliStripHandle()
{
   // handle deletion updates reference count
  if (fUgliStripNode) fUgliStripNode->DecrementRef();
  if (fGeoStripNode) fGeoStripNode->DecrementRef();
  
}

//_____________________________________________________________________________
UgliStripHandle::UgliStripHandle(UgliStripNode* node)
  : fUgliStripNode(node),fGeoStripNode(0)
{
   // ctor keeps reference count up-to-date

   if (fUgliStripNode) fUgliStripNode->IncrementRef();
} 

//_____________________________________________________________________________
UgliStripHandle::UgliStripHandle(GeoStripNode* node)
  : fUgliStripNode(0),fGeoStripNode(node)
{
   // ctor keeps reference count up-to-date

   if (fGeoStripNode) fGeoStripNode->IncrementRef();
}

//_____________________________________________________________________________
UgliStripHandle::UgliStripHandle(const UgliStripHandle &that)
   : UgliStripABC(),
     fUgliStripNode(that.fUgliStripNode),
     fGeoStripNode(that.fGeoStripNode)
{
   // copy constructor keeps reference counts up-to-date

   if (fUgliStripNode) fUgliStripNode->IncrementRef();
   if (fGeoStripNode) fGeoStripNode->IncrementRef();
}

//_____________________________________________________________________________
UgliStripHandle&  UgliStripHandle::operator=(const UgliStripHandle& that)
{
   // assignment keeps reference counts up-to-date

   
   if (fUgliStripNode) fUgliStripNode->DecrementRef();
   if (fGeoStripNode) fGeoStripNode->DecrementRef();
   
   fUgliStripNode = that.fUgliStripNode;
   fGeoStripNode = that.fGeoStripNode;
   
   if (fUgliStripNode) fUgliStripNode->IncrementRef();
   if (fGeoStripNode) fGeoStripNode->IncrementRef();

   return *this;
}

//_____________________________________________________________________________
// inlining these causes circular dependencies

//inline 
UgliScintMdlHandle UgliStripHandle::GetScintMdlHandle(void) const
{ if (!IsGeo()) return UgliScintMdlHandle(fUgliStripNode->GetScintMdlNode());
  else return UgliScintMdlHandle(fGeoStripNode->GetScintMdlNode()); 
}

//inline 
UgliScintPlnHandle UgliStripHandle::GetScintPlnHandle(void) const
{ if (!IsGeo()) return 
  UgliScintPlnHandle(fUgliStripNode->GetScintMdlNode()->GetScintPlnNode()); 
 else return 
  UgliScintPlnHandle(fGeoStripNode->GetScintMdlNode()->GetScintPlnNode()); 
}

//_____________________________________________________________________________
Bool_t UgliStripHandle::AttemptSet(const char* what) const
{
  if (!IsFrozen()) return true;

  MSG("Ugli",Msg::kError) 
    << "AttemptSet " << what << " on Frozen " 
    << GetSEId() << " twarted!" << endl;
  return false;
}

void UgliStripHandle::SetZRotRelMdlRad(Float_t radians)
{
  if (!AttemptSet("ZRotRelMdlRad")) return;
  if (!IsGeo())fUgliStripNode->SetZRotRelMdlRad(radians);
  else fGeoStripNode->SetZRotRelMdlRad(radians);
  
}

void UgliStripHandle::SetLTPosRelMdl(Float_t lpos, Float_t tpos)
{
  if (!AttemptSet("LTPos")) return;
  if (!IsGeo())fUgliStripNode->SetLTPosRelMdl(lpos,tpos);
  else fGeoStripNode->SetLTPosRelMdl(lpos/Munits::cm,tpos/Munits::cm);
  
}

//_____________________________________________________________________________
TVector3 UgliStripHandle::GlobalPos(const Float_t alongLength, 
                                    const Bool_t onWLS,
                                    const Bool_t globalInXYZ) const {
  // convert a distance along the strip into global position
  // onWLS takes into account the extra waveshifter length in bypass region
  // globalInXYZ determines whether result is in XYZ or UVZ
  TVector3 globalXYZ;
  if ( !IsGeo() ) 
      globalXYZ = fUgliStripNode->GlobalPos(alongLength,onWLS);
  else {
      TVector3 gpos = fGeoStripNode->GlobalPos(alongLength/Munits::cm,onWLS);
      globalXYZ = (gpos *= Munits::cm);
  }

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

//_____________________________________________________________________________
TVector3 UgliStripHandle::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(GetSEId().GetDetector(),global);

  if ( !IsGeo() ) 
      return fUgliStripNode->GlobalToLocal(globalXYZ);
  else {
      TVector3 scaledglobal(globalXYZ); scaledglobal *= 1./Munits::cm;
      return (fGeoStripNode->GlobalToLocal(scaledglobal) *= Munits::cm); 
  }
}

//_____________________________________________________________________________
TVector3 UgliStripHandle::LocalToGlobal(const TVector3& local,
                                        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;
  if ( !IsGeo() ) 
      globalXYZ = fUgliStripNode->LocalToGlobal(local); 
  else { 
      TVector3 scaledlocal(local); scaledlocal *= 1./Munits::cm;
      globalXYZ = (fGeoStripNode->LocalToGlobal(scaledlocal) *= Munits::cm); 
  }

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

//_____________________________________________________________________________

---