GeoScintPlnNode.cxx

Go to the documentation of this file.

//
// GeoScintPlnNode
//
// GeoScintPlnNode is a plane node for a scint plane
//
// Author:  S. Kasahara 07/04

#include "MessageService/MsgService.h"
#include "GeoGeometry/GeoScintPlnNode.h"
#include "GeoGeometry/GeoScintMdlNode.h"
#include "GeoGeometry/GeoStripNode.h"
#include <TMath.h>
#include <TGeoMatrix.h>

ClassImp(GeoScintPlnNode)
CVSID("$Id: GeoScintPlnNode.cxx,v 1.15 2007/12/11 06:04:16 rhatcher Exp $");

//_____________________________________________________________________________
GeoScintPlnNode::GeoScintPlnNode(GeoGeometry* geo, TGeoVolume* plnvol,
                           TGeoMatrix* plnmatrix, TGeoVolume* parVol,
                           std::string globalpath, std::string nodename, 
                           const PlexPlaneId& planeid)
  : GeoPlnNode(geo,plnvol,plnmatrix,parVol,globalpath,nodename,planeid) {
  // Normal constructor
  // Protected

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

//_____________________________________________________________________________
GeoStripNode* GeoScintPlnNode::GetStripNode(const PlexStripEndId& seid) const {
  // Returns strip node corresponding to specified seid

  MSG("Geo",Msg::kDebug) << " PlaneId " << GetPlexPlaneId().AsString() 
                         << " GetStripNode " << seid.AsString() << endl;

  UpdateGlobalManager();
  
  TGeoVolume* vol = this -> GetVolume();
  Int_t nmodule = vol -> GetNdaughters();
  for ( int imdl = 0; imdl < nmodule; imdl++ ) {
    GeoScintMdlNode* mdlNode = dynamic_cast<GeoScintMdlNode*>
                               (vol -> GetNode(imdl));
    if ( mdlNode ) {
      GeoStripNode* stpNode = mdlNode -> GetStripNode(seid);
      if ( stpNode ) return stpNode;
    }
  }

  MSG("Geo",Msg::kWarning) << "GetStripNode for plane " 
       << GetPlexPlaneId().AsString() << " failed to find strip " 
       << seid.AsString() << "." << endl;
  
  return 0;
  
}

//_____________________________________________________________________________
Float_t GeoScintPlnNode::GetZRotRelSteelRad() const {
  // Return zrot relative to steel

  UpdateGlobalManager();
  
  const Double_t* rotmatrix = GetMatrix() -> GetRotationMatrix();
  Float_t zrotrelsteelrad = TMath::ATan2(rotmatrix[3],rotmatrix[0]);
  
  return zrotrelsteelrad;
  
}

//_____________________________________________________________________________
Float_t GeoScintPlnNode::GetX0RelSteel() const {
  // Return the x pos relative to steel (x-coord from translation matrix)

  UpdateGlobalManager();
  
  const Double_t* trans = GetMatrix() -> GetTranslation();
  
  return trans[0];
  
}

//_____________________________________________________________________________
Float_t GeoScintPlnNode::GetY0RelSteel() const {
  // Return the y pos relative to steel (y-coord from translation matrix)

  UpdateGlobalManager();

  const Double_t* trans = GetMatrix() -> GetTranslation();
  
  return trans[1];
  
}

//_____________________________________________________________________________
Int_t GeoScintPlnNode::NumberOfStrips() const {
  // Return total number of strips in plane

  UpdateGlobalManager();
  
  // Should cache this data once calculated
  return (GetScintPlnVolume() -> NumberOfStrips());
  
}

//_____________________________________________________________________________
Int_t GeoScintPlnNode::NumberOfModules() const {
  // Return total number of strips in plane

  UpdateGlobalManager();
  
  // Should cache this data once calculated
  return (GetScintPlnVolume() -> NumberOfModules());
  
}

//_____________________________________________________________________________
vector<GeoStripNode*> GeoScintPlnNode::GetStripNodePtrVector() const {
  // Return vector of GeoStripNode's in this plane

  UpdateGlobalManager();
  
  return (GetScintPlnVolume()->GetStripNodePtrVector());
  
}

//_____________________________________________________________________________
vector<GeoScintMdlNode*> GeoScintPlnNode::GetScintMdlNodePtrVector() const {
  // Return vector of GeoScintMdlNode's in this plane

  UpdateGlobalManager();
  
  return (GetScintPlnVolume()->GetScintMdlNodePtrVector());
  
}

//_____________________________________________________________________________
void GeoScintPlnNode::SetZRotRelSteelRad(Float_t /* radians */) {
  // Dummy

  UpdateGlobalManager();
  
  MSG("Geo",Msg::kWarning) 
   << "GeoScintPlnNode::SetZRotRelSteelRad not implemented yet." << endl;
  return;
  
}

//_____________________________________________________________________________
void GeoScintPlnNode::SetXY0RelSteel(Float_t /* x0 */, Float_t /* y0 */) {
  // Dummy

  UpdateGlobalManager();
  
  MSG("Geo",Msg::kWarning) 
    << "GeoScintPlnNode::SetXY0RelSteel not implemented yet." << endl;
  return;
  
}


//_____________________________________________________________________________
GeoStripNode* GeoScintPlnNode::GetClosestStrip(Float_t tpos,
                                               Float_t orthCoord) const {
  // Return the node of the strip in the plane with a position
  // closest to "tpos", orthCoord is position in orthogonal coordinate

  UpdateGlobalManager();

  vector<GeoStripNode*> stripNodes = this -> GetStripNodePtrVector();
  GeoStripNode* closestNode = 0;
  Float_t closestDtpos = 999999.;
  
  for ( unsigned int i = 0; i < stripNodes.size(); ++i ) {
    Float_t strip_tpos = stripNodes[i] -> GetTPos(orthCoord);
    Float_t dtpos = TMath::Abs(strip_tpos-tpos);
    if ( dtpos < closestDtpos ) {
      closestNode = stripNodes[i];
      closestDtpos = dtpos;
    }
    else if ( closestNode != 0 ) {
      // Moving away from closet node
      return closestNode;
    }
  }
  
  return closestNode;
  
}

---