GeoScintMdlNode.cxx

Go to the documentation of this file.

// $Id: GeoScintMdlNode.cxx,v 1.7 2007/12/11 06:04:16 rhatcher Exp $
//
// GeoScintMdlNode
//
// GeoScintMdlNode is a module node
//
// Author:  S. Kasahara 03/05

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

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

//_____________________________________________________________________________
GeoScintMdlNode::GeoScintMdlNode(GeoGeometry* geom, TGeoVolume* mdlvol,
                          TGeoMatrix* mdlmatrix, GeoScintPlnNode* plnNode,
                          std::string nodename, const PlexScintMdlId& mdlId,
                          Float_t clearleneast, Float_t clearlenwest,
                          Float_t wlsleneast, Float_t wlslenwest) 
  : GeoNode(geom,mdlvol,mdlmatrix,plnNode->GetVolume(),
            plnNode->GetGlobalPath()+"/"+nodename,nodename),
    fScintMdlId(mdlId),fScintPlnNode(plnNode),fClearLenEast(clearleneast),
    fClearLenWest(clearlenwest),fWlsLenEast(wlsleneast),
    fWlsLenWest(wlslenwest) {
    // Normal  constructor
    // Protected
    
  UpdateGlobalManager();
  
}
  
//_____________________________________________________________________________
GeoScintMdlNode::~GeoScintMdlNode() {
  // Destructor should delete any owned objects

  UpdateGlobalManager();
  
  if ( CountRef() ) {
    MSG("Geo",Msg::kWarning)
      << "GeoScintMdlNode destructor " 
      << " still had " << CountRef()
      << " outstanding references " << endl;
  }

}

//_____________________________________________________________________________
Float_t GeoScintMdlNode::GetClearFiber(StripEnd::StripEnd_t stripend) const {
  // Return clear fiber length for specified module end

   Float_t clearfiber = 0;
   if      (stripend == StripEnd::kEast) clearfiber = fClearLenEast;
   else if (stripend == StripEnd::kWest) clearfiber = fClearLenWest;

   MSG("Geo",Msg::kDebug) << "GetClearFiber " << fScintMdlId.AsString()
                          << " end " << StripEnd::AsString(stripend)
                          << " clearfiber " << clearfiber << endl;

   return clearfiber;
}

//_____________________________________________________________________________
Float_t GeoScintMdlNode::GetExtraWlsFiber(StripEnd::StripEnd_t stripend)
                                                                   const {
  // Return extra wls fiber length (beyond strip specific wls fiber length)
  // for specified module end

  Float_t extrawlsfiber = 0;
  if       (stripend == StripEnd::kEast) extrawlsfiber = fWlsLenEast;
  else if  (stripend == StripEnd::kWest) extrawlsfiber = fWlsLenWest;
  
  MSG("Geo",Msg::kDebug) << "GetExtraWlsFiber " << fScintMdlId.AsString()
                         << " end " << StripEnd::AsString(stripend)
                         << " extrawlsfiber " << extrawlsfiber << endl;

  return extrawlsfiber;
  
}

//_____________________________________________________________________________
Float_t GeoScintMdlNode::GetLPosRelPln() const {
  // Return the lpos relative to plane (x-coord from translation matrix)

  UpdateGlobalManager();

  const Double_t* trans = GetMatrix() -> GetTranslation();
  
  MSG("Geo",Msg::kDebug) << "GetLPosRelPln " << fScintMdlId.AsString()
                         << " lposrelpln " << trans[0] << endl;
  return trans[0]; 
  
}

//_____________________________________________________________________________
Float_t GeoScintMdlNode::GetTPosRelPln() const {
  // Return the tpos relative to plane (y-coord from translation matrix)

  UpdateGlobalManager();

  const Double_t* trans = GetMatrix() -> GetTranslation();

  MSG("Geo",Msg::kDebug) << "GetTPosRelPln " << fScintMdlId.AsString()
                         << " tposrelpln " << trans[1] << endl;

  return trans[1];
  
}


//_____________________________________________________________________________
Float_t GeoScintMdlNode::GetZRotRelPlnRad() const {
  // Return the zrot relative to plane in radians

  UpdateGlobalManager();

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

  MSG("Geo",Msg::kDebug) << "GetZRotRelPlnRad " << fScintMdlId.AsString()
                         << " zrotrelplnrad " << zrotrelplnrad << endl;

  return zrotrelplnrad;
   
}

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

  UpdateGlobalManager();

  GeoStripNode* stripnode = GetScintMdlVolume()->GetStripNode(seid);

  MSG("Geo",Msg::kDebug) << "Module node " << fScintMdlId.AsString() 
                         << " GetStripNode "  << seid.AsString() 
                         << " node ptr " << stripnode << endl;
 
  return stripnode;
   
}

//_____________________________________________________________________________
std::vector<GeoStripNode*> GeoScintMdlNode::GetStripNodePtrVector() const {
  // Returns strip corresponding to the specified seid

  UpdateGlobalManager();

  return (GetScintMdlVolume()->GetStripNodePtrVector());
   
}

//_____________________________________________________________________________
void GeoScintMdlNode::SetZRotRelPlnRad(Float_t /*radians*/) {
  // Dummy

  UpdateGlobalManager();

  MSG("Geo",Msg::kWarning)
    << "GeoScintMdlNode::SetZRotRelPlnRad not implemented yet." << endl;
  return;

}

//_____________________________________________________________________________
void GeoScintMdlNode::SetLTPosRelPln(Float_t /*lpos*/, Float_t /*tpos*/) {
  // Dummy

  UpdateGlobalManager();

  MSG("Geo",Msg::kWarning)
    << "GeoScintMdlNode::SetLTPosRelPln not implemented yet." << endl;
  return;

}

---