UgliStripNode Class Reference

#include <UgliStripNode.h>

Inheritance diagram for UgliStripNode:

List of all members.

Public Member Functions
virtual	~UgliStripNode ()
virtual PlexStripEndId	GetSEId () const
virtual Float_t	GetHalfLength () const
virtual Float_t	GetHalfThickness () const
virtual Float_t	GetHalfWidth () const
virtual Float_t	GetTPos (Float_t orthCoord=9999.) const
virtual Float_t	PartialLength (const StripEnd::StripEnd_t end) const
virtual Bool_t	IsMirrored (const StripEnd::StripEnd_t end) const
virtual Float_t	WlsPigtail (const StripEnd::StripEnd_t end) const
virtual Float_t	WlsBypass () const
virtual Float_t	ClearFiber (const StripEnd::StripEnd_t end) const
virtual Float_t	TotalAttenuation (const StripEnd::StripEnd_t end, const Float_t alongLength) const
virtual TVector3	GlobalPos (const Float_t alongLength, const Bool_t onWLS=false) const
virtual TVector3	GlobalToLocal (const TVector3 &global) const
virtual TVector3	LocalToGlobal (const TVector3 &local) const
virtual Float_t	DistanceAlong (const PlexStripEndId &orthogonalStrip) const
virtual TVector3	Intersection (const PlexStripEndId &orthogonalStrip) const
virtual UgliGeometry *	GetUgliGeometry () const
UgliScintMdlNode *	GetScintMdlNode (void) const
Float_t	GetZRotRelMdlRad () const
Float_t	GetLPosRelMdl () const
Float_t	GetTPosRelMdl () const
void	SetZRotRelMdlRad (Float_t radians)
void	SetLTPosRelMdl (Float_t lpos, Float_t tpos)
	UgliStripNode ()
Protected Member Functions
virtual void	IncrementRef ()
virtual void	DecrementRef ()
	UgliStripNode (const PlexStripEndId seid, UgliScintMdlNode *mdl_node, const UgliDbiTables &ugliTables)
UgliStripShape *	GetUgliStripShape () const
Protected Attributes
UgliGeometry *	fUgliGeometry
PlexStripEndId	fSEId
UInt_t	fMirrorMask
Private Member Functions
VldContext	BackConstructVldContext () const
	UgliStripNode (const UgliStripNode &that)
UgliStripNode &	operator= (const UgliStripNode &ugh)
Friends
class	UgliGeometry
class	UgliScintPlnNode
class	UgliStripHandle

---

Detailed Description

Definition at line 27 of file UgliStripNode.h.

---

Constructor & Destructor Documentation

UgliStripNode::~UgliStripNode

(

)

[virtual]

Definition at line 98 of file UgliStripNode.cxx.

References UgliRefCnt::CountRef(), Msg::kWarning, and MSG.

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

UgliStripNode::UgliStripNode

(

)

Definition at line 91 of file UgliStripNode.cxx.

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

UgliStripNode::UgliStripNode	(	const PlexStripEndId	seid,
		UgliScintMdlNode *	mdl_node,
		const UgliDbiTables &	ugliTables
	)			[protected]

Definition at line 109 of file UgliStripNode.cxx.

References PlexStripEndId::AsString(), fMirrorMask, UgliDbiTables::fStripStructTbl, fUgliGeometry, UgliDbiTables::GetDbiScintMdlById(), UgliDbiTables::GetDbiStripById(), PlexPlaneId::GetDetector(), UgliDbiStripStruct::GetInModule(), UgliDbiStrip::GetLPosRelMdl(), UgliDbiScintMdl::GetModule(), Nav::GetName(), DbiResultPtr< T >::GetRowByIndex(), UgliDbiScintMdl::GetScintMdlId(), UgliGeometry::GetShape(), UgliDbiStrip::GetTPosRelMdl(), UgliDbiScintMdl::GetTPosRelPln(), UgliDbiStrip::GetZRotRelMdlRad(), UgliDbiScintMdl::GetZRotRelPlnDeg(), PlexStripEndId::IsValid(), StripEnd::kEast, Msg::kError, Msg::kInfo, Detector::kNear, Msg::kWarning, MSG, SetLTPosRelMdl(), TNodeX::SetParent(), and SetZRotRelMdlRad().

   : TNodeX(seid.AsString("c"),seid.AsString("c"),noshapeName), //note below
     fUgliGeometry(mdlnode->GetUgliGeometry()), fSEId(seid)     
{
   // DBI driven ctor under UgliScintMdlNodes

   SetParent(mdlnode);

   fMirrorMask = 0;
   if (seid.GetDetector() == Detector::kNear) 
      fMirrorMask = StripEnd::kEast;

   // NOTE: we set the shape name here because seid.AsString
   // uses a static space such that having both "s" and "c"
   // options on the initialization line means we only see one
   // of the two possibilities
   fShape = fUgliGeometry->GetShape(seid.AsString("s"));

   UgliStripShape* sshape = dynamic_cast<UgliStripShape*>(fShape);
   if (!sshape) {
      gGeometry->ls("s");
      const char* name = "(nil)";
      if (fShape) name = fShape->GetName();
      MSG("Ugli",Msg::kWarning)
         << "ctor for " << seid.AsString("c") 
         << "(" << GetName() << "," << GetTitle() << ")" 
         << endl
         << "   found shape wasn't the UgliStripShape " 
         << seid.AsString("s") 
         << " instead was named " << name << endl;
   }

   unsigned int stripstruct_indx = UgliDbiStructHash(seid).HashAsStrip();
   const UgliDbiStripStruct *stripStructRow = 
      ugliTables.fStripStructTbl.GetRowByIndex(stripstruct_indx);
   int inmdl = stripStructRow->GetInModule();

   // for now some module level info lives at this level
   // since we don't have modules to ask
   const UgliDbiScintMdl *mdlRow = 
      ugliTables.GetDbiScintMdlById(PlexScintMdlId(seid,inmdl));
   if (!mdlRow) {
      MSG("Ugli",Msg::kError)
         << "UgliStripNode ctor failed to get UgliDbiScintMdl" 
         << PlexScintMdlId(seid,inmdl)
         << endl;
      return;
   }

   const UgliDbiStrip *stripRow = ugliTables.GetDbiStripById(seid);
   if (!stripRow) {
      MSG("Ugli",Msg::kError)
         << "UgliStripNode ctor failed to get UgliDbiStrip" 
         << seid
         << endl;
      return;
   }
   // for now we're positioning ourselves in the plane
   // as there are no modules
   if (mdlRow->GetZRotRelPlnDeg() != 0) {
      MSG("Ugli",Msg::kWarning)
         << " module  " << mdlRow->GetScintMdlId()
         << ", " << mdlRow->GetModule()
         << " has rotation of " << mdlRow->GetZRotRelPlnDeg()
         << " degrees, but we're not placing strip in module" << endl
         << "     rotation is effectively ignored" << endl;
   }
   Float_t tpos = stripRow->GetTPosRelMdl();
   Float_t lpos = stripRow->GetLPosRelMdl();

   SetLTPosRelMdl(lpos,tpos);
   
   if (!seid.IsValid()) {
      static int nmsg = 20;
      if (nmsg) {
        MSG("Ugli",Msg::kInfo) 
          << "UgliStripNode::ctor: "
          << seid.AsString("c") << " @ "
          << " tpos " << tpos << " = " 
          << mdlRow->GetTPosRelPln() << " + " << stripRow->GetTPosRelMdl()
          << endl
          << "  not a valid strip: "
          << "  seid.IsValid() " << (int)seid.IsValid()
          // << ", encoded 0x" << hex << seid.GetEncoded() << dec 
          // << ", masked out plane 0x" << hex << (seid.GetEncoded()&maskPlexIdPlane) << dec 
          // << ", !=mask is " << hex << (int)((seid.GetEncoded()&maskPlexIdPlane)!=maskPlexIdPlane) << dec 
          << endl;
        nmsg--;
        if (nmsg==0) {
          MSG("Ugli",Msg::kInfo) << "  ... last message of this type." << endl;
        }
      }
   }

   Float_t zrot = stripRow->GetZRotRelMdlRad();
   SetZRotRelMdlRad(zrot);

   SetLineColor(18);  // a nice light gray
   SetVisibility(0);  // visible
//   if (seid.GetStrip()%8 != 0) SetVisibility(0);

}

UgliStripNode::UgliStripNode

(

const UgliStripNode &

that

)

[private]

---

Member Function Documentation

VldContext UgliStripNode::BackConstructVldContext

(

)

const [private]

Definition at line 806 of file UgliStripNode.cxx.

References det, fUgliGeometry, VldRange::GetDetectorMask(), VldRange::GetSimMask(), VldRange::GetTimeStart(), and UgliGeometry::GetVldRange().

Referenced by ClearFiber(), PartialLength(), WlsBypass(), and WlsPigtail().

{
   // build a fake VldContext from the VldRange
   // so we can access the DBI

   VldRange vrng = fUgliGeometry->GetVldRange();
   int det = vrng.GetDetectorMask();
   int sim = vrng.GetSimMask();
   // trim "det" and "sim" down to single bit
   for (int i=0; i<32; ++i) {
      int bit = 1<<i;
      if (det&bit) det=bit;
      if (sim&bit) sim=bit;
   }
   
   return VldContext((Detector::Detector_t)det,
                     (SimFlag::SimFlag_t)sim,
                     vrng.GetTimeStart());

}

Float_t UgliStripNode::ClearFiber

(

const StripEnd::StripEnd_t

end

)

const [virtual]

Implements UgliStripABC.

Definition at line 749 of file UgliStripNode.cxx.

References PlexStripEndId::AsString(), BackConstructVldContext(), REROOT_CellPos::ClearFiberA(), REROOT_CellPos::ClearFiberB(), Munits::cm, fSEId, fUgliGeometry, RerootExodus::GetCellPos(), UgliScintMdlNode::GetClearFiber(), UgliDbiScintMdl::GetClearLenEast(), UgliDbiScintMdl::GetClearLenWest(), UgliDbiStripStruct::GetInModule(), PlexPlaneId::GetPlaneView(), DbiResultPtr< T >::GetRowByIndex(), VldRange::GetSimMask(), UgliGeometry::GetVldRange(), UgliDbiScintMdl::HashToIndex(), StripEnd::kEast, SimFlag::kReroot, Msg::kWarning, StripEnd::kWest, and MSG.

Referenced by UgliStripHandle::ClearFiber().

{

   UgliScintMdlNode* mdl = dynamic_cast<UgliScintMdlNode*>(fParent);

   if (mdl) return mdl->GetClearFiber(end);
      
   MSG("Ugli",Msg::kWarning) << "ClearFiber() " << endl
                             << "  parent isn't a UgliScintMdlNode" << endl;

   // fake up a context from the range
   VldContext vldc = BackConstructVldContext();

   DbiResultPtr<UgliDbiStripStruct> stripStructTbl(vldc);
   unsigned int stripstruct_indx = UgliDbiStructHash(fSEId).HashAsStrip();
   const UgliDbiStripStruct *stripStructRow = 
      stripStructTbl.GetRowByIndex(stripstruct_indx);
   int inmdl = stripStructRow->GetInModule();

   DbiResultPtr<UgliDbiScintMdl> mdlTbl(vldc);
   int mdl_indx = UgliDbiScintMdl::HashToIndex(PlexScintMdlId(fSEId,inmdl));
   const UgliDbiScintMdl *mdlRow = mdlTbl.GetRowByIndex(mdl_indx);

   if (mdlRow) {
      switch (end) {
      case StripEnd::kEast: return mdlRow->GetClearLenEast(); break;
      case StripEnd::kWest: return mdlRow->GetClearLenWest(); break;
      default: return huge_length;
      }
   }
#ifdef UGLI_USE_EXODUS
   else if (fUgliGeometry->GetVldRange().GetSimMask() & SimFlag::kReroot) {

      REROOT_CellPos *cellpos = RerootExodus::GetCellPos(fSEId);
      if (!cellpos) {
         MSG("Ugli",Msg::kWarning) 
            << "No CellPos for " << fSEId.AsString() << endl;
         return huge_length;
      }

      Float_t fClearFiber[3];
      fClearFiber[0] = cellpos->ClearFiberA() * Munits::cm;
      fClearFiber[1] = cellpos->ClearFiberB() * Munits::cm;
      fClearFiber[2] = huge_length;
      PlaneView::PlaneView_t view = fSEId.GetPlaneView();

      return fClearFiber[isABside(view,end)];

   }
#endif
   MSG("Ugli",Msg::kWarning) 
     << " ClearFiber() missing info:  MdlRow " << mdlRow << endl;
   return huge_length;

}

void UgliStripNode::DecrementRef

(

)

[protected, virtual]

Reimplemented from UgliRefCnt.

Definition at line 260 of file UgliStripNode.cxx.

References UgliRefCnt::DecrementRef(), UgliRefCnt::fRef, and fUgliGeometry.

Referenced by UgliStripHandle::operator=(), and UgliStripHandle::~UgliStripHandle().

{ 
   fRef--; 
   fUgliGeometry->DecrementRef();
//   cout << fSEId.AsString("c") << " -- ref to " << fRef 
//        << ( (!fRef) ? " make invisible" : " " ) << endl;
   if (!fRef) SetVisibility(0);
}

Float_t UgliStripNode::DistanceAlong

(

const PlexStripEndId &

orthogonalStrip

)

const [virtual]

Implements UgliStripABC.

Definition at line 479 of file UgliStripNode.cxx.

References UgliGeometry::GetStripNode(), GetUgliGeometry(), GlobalPos(), Msg::kError, and MSG.

Referenced by UgliStripHandle::DistanceAlong().

{ 
   // Get the distance along a strip represented by orthogonal strip.
   // Note this is relative to the strip's local co-ordinate system
   // origin (which for NearDet is not along the U or V axis), so
   // one can use it for attenuation corrections but not directly
   // for TPos position information.
  
   // The closest approach ponts on the "s" and "t" lines are:
   //    xyz_t(i) = base_t(i) + t*dir_t(i) for i=0..2 {x,y,z} 
   //    xyz_o(i) = base_o(i) + o*dir_o(i) for i=0..2 {x,y,z}
   // where base_x is  LocalToGlobal(TVector3(0,0,0))
   // and   dir_x  is (LocalToGlobal(TVector3(1,0,0))-base_x)
   //
   // Finding the intersection or closest appraoch of two skew lines
   // algorithm was adapted from:
   //   Graphic Gems, ed. Andrew S. Glassner, pg 304
   //   ISBN 0-12-286165-5   T385.G697(1990)
   //
   const Double_t farAway = 1.0e25;

   UgliStripNode* otherNode = GetUgliGeometry()->GetStripNode(orthogonalStrip);

   if ( ! otherNode ) {
     MSG("Ugli",Msg::kError)
       << "UgliStripNode::DistanceAlong(), other strip " 
       << orthogonalStrip << " does not exist" << endl;
     return farAway;
   }
   
   // s=this, t=other
   TVector3 baseThis  = this->GlobalPos(0.,false);
   TVector3 dirThis   = this->GlobalPos(1.,false) - baseThis;
   
   TVector3 baseOther = otherNode->GlobalPos(0.,false);
   TVector3 dirOther  = otherNode->GlobalPos(1.,false) - baseOther;
   
   TVector3 dirCross  = dirThis.Cross(dirOther);
   Double_t vxvm2 = dirCross.Mag2();
   if (0.0 == vxvm2) {
     // exactly parallel strips!
     return farAway;
   }
   
   TVector3 dp   = baseOther - baseThis;
   TVector3 dpxv = dp.Cross(dirOther);
   Double_t s = dpxv.Dot(dirCross) / vxvm2;
   
   return s;   
}

Float_t UgliStripNode::GetHalfLength

(

)

const [virtual]

Implements UgliStripABC.

Definition at line 270 of file UgliStripNode.cxx.

Referenced by UgliStripHandle::GetHalfLength().

{ 
   // get the strip's half length
   TBRIK* brik = dynamic_cast<TBRIK*>(fShape);
   if (brik) return brik->GetDx(); else return 0;
}

Float_t UgliStripNode::GetHalfThickness

(

)

const [virtual]

Implements UgliStripABC.

Definition at line 278 of file UgliStripNode.cxx.

Referenced by UgliStripHandle::GetHalfThickness().

{ 
   // get the strip's half length
   TBRIK* brik = dynamic_cast<TBRIK*>(fShape);
   if (brik) return brik->GetDz(); else return 0;
}

Float_t UgliStripNode::GetHalfWidth

(

)

const [virtual]

Implements UgliStripABC.

Definition at line 286 of file UgliStripNode.cxx.

Referenced by UgliStripHandle::GetHalfWidth().

{ 
   // get the strip's half length
   TBRIK* brik = dynamic_cast<TBRIK*>(fShape);
   if (brik) return brik->GetDy(); else return 0;
}

Float_t UgliStripNode::GetLPosRelMdl

(

)

const [inline]

Definition at line 131 of file UgliStripNode.h.

Referenced by UgliStripHandle::GetLPosRelMdl().

{ return fX; }

UgliScintMdlNode * UgliStripNode::GetScintMdlNode

(

void

)

const

Definition at line 831 of file UgliStripNode.cxx.

Referenced by UgliStripHandle::GetScintMdlHandle(), UgliStripHandle::GetScintPlnHandle(), and GetTPos().

{  return dynamic_cast<UgliScintMdlNode*>(fParent); }

PlexStripEndId UgliStripNode::GetSEId

(

)

const [inline, virtual]

Implements UgliStripABC.

Definition at line 120 of file UgliStripNode.h.

References fSEId.

Referenced by UgliStripHandle::GetSEId().

{ return fSEId; }

Float_t UgliStripNode::GetTPos

(

Float_t

orthCoord = 9999.

)

const [virtual]

Implements UgliStripABC.

Definition at line 293 of file UgliStripNode.cxx.

References bfld::AsString(), fSEId, PlexPlaneId::GetDetector(), PlexPlaneId::GetPlaneView(), GetScintMdlNode(), UgliScintMdlNode::GetScintPlnNode(), UgliPlnNode::GetZ0(), GlobalPos(), PlaneView::kA, PlaneView::kB, Detector::kCalib, Msg::kError, PlaneView::kU, PlaneView::kV, PlaneView::kVSTopEastSlant, PlaneView::kVSTopFlat, PlaneView::kVSTopWestSlant, PlaneView::kVSWallEastSlant, PlaneView::kVSWallOnEdge, PlaneView::kVSWallWestSlant, Msg::kWarning, PlaneView::kX, PlaneView::kY, and MSG.

Referenced by UgliStripHandle::GetTPos().

{ 
   // get the strip's tpos offset
   // orthCoord is the position in the orthogonal view, this allows one to
   // include rotation corrections.
   // if orthCoord>1000 use the center of the strip

   PlaneView::PlaneView_t view = fSEId.GetPlaneView();
   bool viewIsUV = PlaneView::kU==view || PlaneView::kV==view;
   Detector::Detector_t detector = fSEId.GetDetector();

   const Double_t rsqrt2 = 1.0/TMath::Sqrt(2.);
   const TVector3 dir_x(1.,0.,0.);
   const TVector3 dir_y(0.,1.,0.);
   const TVector3 dir_u(rsqrt2,rsqrt2,0.);
   const TVector3 dir_v(-rsqrt2,rsqrt2,0.);
   const TVector3 dir_u_caldet(1.,0.,0.);
   const TVector3 dir_v_caldet(0.,1.,0.);
   const TVector3 dir_a(0.,0.,1.);
   const TVector3 dir_b(1.,0.,0.);

   TVector3 dir_tpos, dir_orthCoord;

   switch (view) {
   case PlaneView::kX:
     dir_tpos = dir_x; dir_orthCoord = dir_y;
     break;
   case PlaneView::kY:
     dir_tpos = dir_y; dir_orthCoord = dir_x;
     break;
   case PlaneView::kU:
     if (Detector::kCalib != detector) {
       dir_tpos = dir_u; dir_orthCoord = dir_v;
     } else {
       dir_tpos = dir_u_caldet; dir_orthCoord = dir_v_caldet;
     }
     break;
   case PlaneView::kV:
     if (Detector::kCalib != detector) {
       dir_tpos = dir_v; dir_orthCoord = dir_u;
     } else {
       dir_tpos = dir_v_caldet; dir_orthCoord = dir_u_caldet;
     }
     break;
   case PlaneView::kA:
     dir_tpos = dir_a;
     break;
   case PlaneView::kB:
     dir_tpos = dir_b;
     break;
   case PlaneView::kVSTopFlat:
   case PlaneView::kVSTopEastSlant:
   case PlaneView::kVSTopWestSlant:
   case PlaneView::kVSWallOnEdge:
   case PlaneView::kVSWallEastSlant:
   case PlaneView::kVSWallWestSlant:
     MSG("Ugli",Msg::kError)
       << "UgliStripNode::GetTPos() not meaningful for veto shield strips "
       << fSEId << endl;
     return huge_length;
     break;
   default:
      MSG("Ugli",Msg::kWarning)
         << "UgliStripNode::GetTPos unknown view " 
         << PlaneView::AsString(fSEId.GetPlaneView())
         << " in " << fSEId
         << endl;
      return huge_length;
   }

   //cout << " " << fSEId.AsString("c");

   Float_t alongLength = 0.0;  // for wild orthCoord choose *strip* center
   if (orthCoord < 1000. && viewIsUV) {
     // orthCoord is sensible and we're a U or V plane
     // so we need to translate that orthogonal position
     // into a position along the strip in order to do rotation corrections

     TVector3 gstrip_center  = GlobalPos(0.,kFALSE);
     TVector3 gstrip_plusend = GlobalPos(1.,kFALSE);
     TVector3 gstrip_dir = gstrip_plusend - gstrip_center;

     TVector3 tpos_base = orthCoord * dir_orthCoord;
     Double_t z = this->GetScintMdlNode()->GetScintPlnNode()->GetZ0();
     tpos_base.SetZ(z);

     // aliases, just 'cause it's easier to translate the formula
     TVector3& base_s = gstrip_center;
     TVector3& dir_s  = gstrip_dir;
     TVector3& base_t = tpos_base;
     TVector3& dir_t  = dir_tpos;
     
     TVector3 txs  = dir_t.Cross(dir_s);
     TVector3 dpst = base_s - base_t;
     TVector3 dpxs = dpst.Cross(dir_s);
     Float_t  t    = dpxs.Dot(txs) / txs.Mag2();

     return t;

     /* 
        // the unnecessarily hard way ...
        TVector3 sxt  = dir_s.Cross(dir_t);
        TVector3 dpts = base_t - base_s;
        TVector3 dpxt = dpts.Cross(dir_t);
        Float_t  s    = dpxt.Dot(sxt) / sxt.Mag2();
        
        alongLength  = s;
     */

   }
   TVector3 globalpos = GlobalPos(alongLength,kFALSE);
   //return   globalpos.Dot(dir_tpos);
   Float_t tpos = globalpos.Dot(dir_tpos);
   //cout << " return " << tpos << endl;
   return  tpos;
}

Float_t UgliStripNode::GetTPosRelMdl

(

)

const [inline]

Definition at line 132 of file UgliStripNode.h.

Referenced by UgliStripHandle::GetTPosRelMdl().

{ return fY; }

virtual UgliGeometry* UgliStripNode::GetUgliGeometry

(

)

const [inline, virtual]

Definition at line 65 of file UgliStripNode.h.

References fUgliGeometry.

Referenced by DistanceAlong(), Intersection(), and UgliStripHandle::IsFrozen().

{ return fUgliGeometry; }

UgliStripShape * UgliStripNode::GetUgliStripShape

(

)

const [inline, protected]

Definition at line 117 of file UgliStripNode.h.

Referenced by PartialLength(), WlsBypass(), and WlsPigtail().

{ return dynamic_cast<UgliStripShape*>(fShape); }

Float_t UgliStripNode::GetZRotRelMdlRad

(

)

const [inline]

Definition at line 128 of file UgliStripNode.h.

References Ugli::deg2rad.

Referenced by UgliStripHandle::GetZRotRelMdlRad().

{ return fMatrix->GetPhi()*Ugli::deg2rad; }

TVector3 UgliStripNode::GlobalPos	(	const Float_t	alongLength,
		const Bool_t	onWLS = false
	)			const [virtual]

Definition at line 421 of file UgliStripNode.cxx.

References Msg::kFatal, LocalToGlobal(), and MSG.

Referenced by DistanceAlong(), GetTPos(), UgliStripHandle::GlobalPos(), and Intersection().

{ 
   // get the position based on distance along centerline *or* WLS fiber
   // global position is assumed to be in XYZ (hall) coordinates

   if (onWLS) {
      MSG("Ugli",Msg::kFatal)
         << "UgliStripNode::GlobalPos not yet implemented with onWLS=kTRUE" << endl;
      return LocalToGlobal(TVector3(alongLength,0.,0.));
   } else {
      return LocalToGlobal(TVector3(alongLength,0.,0.));
   }
}

TVector3 UgliStripNode::GlobalToLocal

(

const TVector3 &

global

)

const [virtual]

Definition at line 437 of file UgliStripNode.cxx.

Referenced by UgliStripHandle::GlobalToLocal().

{ 
   // get the local position based on global position
   // global position is assumed to be in XYZ (hall) coordinates

   Double_t gxyz[3], lxyz[3];

   gxyz[0] = global.X();
   gxyz[1] = global.Y();
   gxyz[2] = global.Z();

   // deal with fact that TNode::Master2Local isn't declared const
   UgliStripNode* self = const_cast<UgliStripNode*>(this);
   self->cd();
   self->UpdateMatrix();
   self->Master2Local(gxyz,lxyz);

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

void UgliStripNode::IncrementRef

(

)

[protected, virtual]

Reimplemented from UgliRefCnt.

Definition at line 252 of file UgliStripNode.cxx.

References UgliRefCnt::fRef, fUgliGeometry, and UgliRefCnt::IncrementRef().

Referenced by UgliStripHandle::operator=(), and UgliStripHandle::UgliStripHandle().

{
   fRef++; 
   fUgliGeometry->IncrementRef();
//   cout << fSEId.AsString("c") << " ++ ref to " << fRef << endl;
   SetVisibility(1);
}

TVector3 UgliStripNode::Intersection

(

const PlexStripEndId &

orthogonalStrip

)

const [virtual]

Implements UgliStripABC.

Definition at line 531 of file UgliStripNode.cxx.

References UgliGeometry::GetStripNode(), GetUgliGeometry(), GlobalPos(), Msg::kError, and MSG.

Referenced by UgliStripHandle::Intersection().

{ 
   // get the 3-space "intersection" of two strips

   // The closest approach ponts on the "s" and "t" lines are:
   //    xyz_t(i) = base_t(i) + t*dir_t(i) for i=0..2 {x,y,z} 
   //    xyz_o(i) = base_o(i) + o*dir_o(i) for i=0..2 {x,y,z}
   // where base_x is  LocalToGlobal(TVector3(0,0,0))
   // and   dir_x  is (LocalToGlobal(TVector3(1,0,0))-base_x)
   //
   // Finding the intersection or closest appraoch of two skew lines
   // algorithm was adapted from:
   //   Graphic Gems, ed. Andrew S. Glassner, pg 304
   //   ISBN 0-12-286165-5   T385.G697(1990)
   //
   const Double_t farAway = 1.0e25;
   const TVector3 farAway3(farAway,farAway,farAway);

   UgliStripNode* otherNode = GetUgliGeometry()->GetStripNode(orthogonalStrip);
   
   if ( ! otherNode ) {
     MSG("Ugli",Msg::kError)
       << "UgliStripNode::Intersection() other strip " 
       << orthogonalStrip << " does not exist" << endl;
     return farAway3;
   }
   
   // s=this, t=other
   TVector3 baseThis  = this->GlobalPos(0.,false);
   TVector3 dirThis   = this->GlobalPos(1.,false) - baseThis;
   
   TVector3 baseOther = otherNode->GlobalPos(0.,false);
   TVector3 dirOther  = otherNode->GlobalPos(1.,false) - baseOther;
   
   TVector3 dirCross  = dirThis.Cross(dirOther);
   Double_t vxvm2 = dirCross.Mag2();
   if (0.0 == vxvm2) {
     // Exactly parallel strips!
     return 0.5*(baseThis+baseOther) + farAway*dirThis;
   }
   
   TVector3 dp   = baseOther - baseThis;
   TVector3 dpxvThis = dp.Cross(dirOther);
   Double_t s = dpxvThis.Dot(dirCross) / vxvm2;
   
   TVector3 thisPoint = baseThis + s*dirThis;

   dirCross *= -1;
   dp       *= -1;
   TVector3 dpxvOther = dp.Cross(dirThis);
   Double_t t = dpxvOther.Dot(dirCross) / vxvm2;

   TVector3 otherPoint = baseOther + t*dirOther;

   return 0.5*(thisPoint+otherPoint);
}

Bool_t UgliStripNode::IsMirrored

(

const StripEnd::StripEnd_t

end

)

const [inline, virtual]

Implements UgliStripABC.

Definition at line 122 of file UgliStripNode.h.

References fMirrorMask.

Referenced by UgliStripHandle::IsMirrored().

{ return ( fMirrorMask & end ); }

TVector3 UgliStripNode::LocalToGlobal

(

const TVector3 &

local

)

const [virtual]

Definition at line 458 of file UgliStripNode.cxx.

Referenced by GlobalPos(), and UgliStripHandle::LocalToGlobal().

{ 
   // get the global position based on local position
   // global position is assumed to be in XYZ (hall) coordinates

   Double_t gxyz[3], lxyz[3];

   lxyz[0] = local.X();
   lxyz[1] = local.Y();
   lxyz[2] = local.Z();

   // deal with fact that TNode::Local2Master isn't declared const
   UgliStripNode* self = const_cast<UgliStripNode*>(this);
   self->cd();
   self->UpdateMatrix();
   self->Local2Master(lxyz,gxyz);

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

UgliStripNode& UgliStripNode::operator=

(

const UgliStripNode &

ugh

)

[private]

Reimplemented from UgliStripABC.

Float_t UgliStripNode::PartialLength

(

const StripEnd::StripEnd_t

end

)

const [virtual]

Definition at line 589 of file UgliStripNode.cxx.

References BackConstructVldContext(), fSEId, UgliDbiStripStruct::GetLenEastPart(), UgliStripShape::GetLenPart(), UgliDbiStripStruct::GetLenWestPart(), DbiResultPtr< T >::GetRowByIndex(), GetUgliStripShape(), StripEnd::kEast, StripEnd::kNegative, StripEnd::kPositive, Msg::kWarning, StripEnd::kWest, and MSG.

Referenced by UgliStripHandle::PartialLength().

{
   // return the (total) length of one of the parts of a split strip
   // (if strip is not split return total length regardless of end)

   // this is, to be proper, info of the strip shape
   UgliStripShape* strip_shape = GetUgliStripShape();

   if (strip_shape) {

      switch (end) {
      case StripEnd::kNegative: 
      case StripEnd::kPositive: 
         return strip_shape->GetLenPart(end);
         break;
      default: return huge_length;
      }
   }
   
   MSG("Ugli",Msg::kWarning) << "PartialLength() " << endl
      << "  shape isn't a UgliStripShape" << endl;

   // fake up a context from the range
   VldContext vldc = BackConstructVldContext();

   DbiResultPtr<UgliDbiStripStruct> stripStructTbl(vldc);
   unsigned int stripstruct_indx = UgliDbiStructHash(fSEId).HashAsStrip();
   const UgliDbiStripStruct *stripStructRow = 
      stripStructTbl.GetRowByIndex(stripstruct_indx);

   if (stripStructRow) {
      switch (end) {
      case StripEnd::kEast: 
         return stripStructRow->GetLenEastPart();
         break;
      case StripEnd::kWest: 
         return stripStructRow->GetLenWestPart();
          break;
      default: return huge_length;
      }
   } 

   MSG("Ugli",Msg::kWarning) 
     << " PartialLength() missing UgliDbiStripStruct" << endl;
   return huge_length;

}

void UgliStripNode::SetLTPosRelMdl	(	Float_t	lpos,
		Float_t	tpos
	)			[inline]

Definition at line 125 of file UgliStripNode.h.

Referenced by UgliStripHandle::SetLTPosRelMdl(), and UgliStripNode().

{ SetPosition(lpos,tpos,0.0); }

void UgliStripNode::SetZRotRelMdlRad

(

Float_t

radians

)

Definition at line 215 of file UgliStripNode.cxx.

References fUgliGeometry, UgliGeometry::GetRotMatrix(), Msg::kDebug, MSG, and Ugli::rad2deg.

Referenced by UgliStripHandle::SetZRotRelMdlRad(), and UgliStripNode().

{
   // set the rotation matrix for this strip
   // limit ourselves to -xx.x mrad

   // build the rotation matrices as we need them so we can share
   // if we expect a range of +/- 4cm at the furthest extreme of the far
   // theta = asin(4cm/400cm) ~ .01 rad = 10mrad
   // if we limit ourself to +/- 0.1 mrad then the smallest non-zero
   // deviation is .04cm which should be sufficient
   const float rad2mrad = 1000.;
   const float mrad2rad = 0.001;

   const char* rname = "Identity";
   char rname_mrad[16];
   if (radians != 0) {
      sprintf(rname_mrad,"%+5.1fmrad",radians*rad2mrad);
      rname = rname_mrad;
   }
   TRotMatrix* rotm = fUgliGeometry->GetRotMatrix(rname);
   if (!rotm) {
      Double_t zrot_trunc;
      sscanf(rname,"%5lf",&zrot_trunc);
      zrot_trunc *= Ugli::rad2deg * mrad2rad;
      MSG("Ugli",Msg::kDebug) 
         << "UgliStripNode::SetZRotRelMdlRad no rotm '" 
         << rname << "' zrot = " 
         << zrot_trunc << " deg" << endl;
      Double_t deg90 = 90;
      rotm = new TRotMatrix(rname,rname,deg90,zrot_trunc,deg90,
                            zrot_trunc+deg90,0,0);
   }
   fMatrix = rotm;
}

Float_t UgliStripNode::TotalAttenuation	(	const StripEnd::StripEnd_t	end,
		const Float_t	alongLength
	)			const [virtual]

Implements UgliStripABC.

Definition at line 411 of file UgliStripNode.cxx.

References Msg::kFatal, and MSG.

Referenced by UgliStripHandle::TotalAttenuation().

{ 
   // get the attenuation from some point towards pixel at one end
   MSG("Ugli",Msg::kFatal)
      << "UgliStripNode::GetTotalAttenuation not yet implemented" << endl;
   return 0;
}

Float_t UgliStripNode::WlsBypass

(

)

const [virtual]

Definition at line 719 of file UgliStripNode.cxx.

References BackConstructVldContext(), fSEId, DbiResultPtr< T >::GetRowByIndex(), GetUgliStripShape(), UgliDbiStripStruct::GetWlsLenBypass(), UgliStripShape::GetWlsLenBypass(), Msg::kWarning, and MSG.

Referenced by UgliStripHandle::WlsBypass().

{
   // return the length of the WLS of the central bypass

   // this is, to be proper, info of the strip shape
   UgliStripShape* strip_shape = GetUgliStripShape();

   if (strip_shape) return strip_shape->GetWlsLenBypass();
   
   MSG("Ugli",Msg::kWarning) << "WlsBypass() " << endl
      << "  shape isn't a UgliStripShape" << endl;

   // fake up a context from the range
   VldContext vldc = BackConstructVldContext();

   DbiResultPtr<UgliDbiStripStruct> stripStructTbl(vldc);
   unsigned int stripstruct_indx = UgliDbiStructHash(fSEId).HashAsStrip();
   const UgliDbiStripStruct *stripStructRow = 
      stripStructTbl.GetRowByIndex(stripstruct_indx);

   if (stripStructRow) return stripStructRow->GetWlsLenBypass();

   MSG("Ugli",Msg::kWarning) 
     << " WlsBypass() missing UgliDbiStripStruct " << endl;

   return 0.0;

}

Float_t UgliStripNode::WlsPigtail

(

const StripEnd::StripEnd_t

end

)

const [virtual]

Implements UgliStripABC.

Definition at line 638 of file UgliStripNode.cxx.

References PlexStripEndId::AsString(), BackConstructVldContext(), Munits::cm, REROOT_CellPos::FiberTailA(), REROOT_CellPos::FiberTailB(), fSEId, fUgliGeometry, RerootExodus::GetCellPos(), UgliScintMdlNode::GetExtraWlsFiber(), UgliDbiStripStruct::GetInModule(), PlexPlaneId::GetPlaneView(), DbiResultPtr< T >::GetRowByIndex(), VldRange::GetSimMask(), GetUgliStripShape(), UgliGeometry::GetVldRange(), UgliStripShape::GetWls(), UgliDbiScintMdl::GetWlsLenEast(), UgliDbiStripStruct::GetWlsLenEast(), UgliDbiScintMdl::GetWlsLenWest(), UgliDbiStripStruct::GetWlsLenWest(), UgliDbiScintMdl::HashToIndex(), StripEnd::kEast, StripEnd::kNegative, StripEnd::kPositive, SimFlag::kReroot, Msg::kWarning, StripEnd::kWest, and MSG.

Referenced by UgliStripHandle::WlsPigtail().

{
   // return the length of the WLS pigtail

   // this is, to be proper, info of the strip shape
   UgliStripShape* strip_shape = GetUgliStripShape();
   UgliScintMdlNode* mdl = dynamic_cast<UgliScintMdlNode*>(fParent);

   if (strip_shape && mdl) {

      switch (end) {
      case StripEnd::kNegative: 
      case StripEnd::kPositive: 
         return 
            strip_shape->GetWls(end) + 
            mdl->GetExtraWlsFiber(end);
         break;
      default: return huge_length;
      }
   }
   
   MSG("Ugli",Msg::kWarning) << "WlsPigtail() " << endl
      << "  shape isn't a UgliStripShape or parent isn't a UgliScintMdlNode" << endl;

   // fake up a context from the range
   VldContext vldc = BackConstructVldContext();

   DbiResultPtr<UgliDbiStripStruct> stripStructTbl(vldc);
   unsigned int stripstruct_indx = UgliDbiStructHash(fSEId).HashAsStrip();
   const UgliDbiStripStruct *stripStructRow = 
      stripStructTbl.GetRowByIndex(stripstruct_indx);
   int inmdl = stripStructRow->GetInModule();

   DbiResultPtr<UgliDbiScintMdl> mdlTbl(vldc);
   int mdl_indx = UgliDbiScintMdl::HashToIndex(PlexScintMdlId(fSEId,inmdl));
   const UgliDbiScintMdl *mdlRow = mdlTbl.GetRowByIndex(mdl_indx);

   if (stripStructRow && mdlRow) {
      switch (end) {
      case StripEnd::kEast: 
         return 
            stripStructRow->GetWlsLenEast() +
            mdlRow->GetWlsLenEast();
         break;
      case StripEnd::kWest: 
         return 
            stripStructRow->GetWlsLenWest() +
            mdlRow->GetWlsLenWest();
         break;
      default: return huge_length;
      }
   } 
#ifdef UGLI_USE_EXODUS
   else if (fUgliGeometry->GetVldRange().GetSimMask() & SimFlag::kReroot) {

      REROOT_CellPos *cellpos = RerootExodus::GetCellPos(fSEId);
      if (!cellpos) {
         MSG("Ugli",Msg::kWarning) 
            << "No CellPos for " << fSEId.AsString() << endl;
         return huge_length;
      }

      Float_t fWlsPigtail[3];
      fWlsPigtail[0] = cellpos->FiberTailA() * Munits::cm;
      fWlsPigtail[1] = cellpos->FiberTailB() * Munits::cm;
      fWlsPigtail[2] = huge_length;
      PlaneView::PlaneView_t view = fSEId.GetPlaneView();

      return fWlsPigtail[isABside(view,end)];
      
   }
#endif
   MSG("Ugli",Msg::kWarning) 
     << " WlsPigtail() missing info: " << endl
     << "     stripStructRow " << stripStructRow
     << " MdlRow " << mdlRow << endl;
   return huge_length;

}

---

Friends And Related Function Documentation

friend class UgliGeometry [friend]

Definition at line 29 of file UgliStripNode.h.

friend class UgliScintPlnNode [friend]

Definition at line 30 of file UgliStripNode.h.

friend class UgliStripHandle [friend]

Definition at line 31 of file UgliStripNode.h.

---

Member Data Documentation

UInt_t UgliStripNode::fMirrorMask [protected]

Definition at line 98 of file UgliStripNode.h.

Referenced by IsMirrored(), and UgliStripNode().

PlexStripEndId UgliStripNode::fSEId [protected]

Definition at line 97 of file UgliStripNode.h.

Referenced by ClearFiber(), GetSEId(), GetTPos(), PartialLength(), WlsBypass(), and WlsPigtail().

UgliGeometry* UgliStripNode::fUgliGeometry [protected]

Definition at line 96 of file UgliStripNode.h.

Referenced by BackConstructVldContext(), ClearFiber(), DecrementRef(), GetUgliGeometry(), IncrementRef(), SetZRotRelMdlRad(), UgliStripNode(), and WlsPigtail().

---

The documentation for this class was generated from the following files:

• UgliStripNode.h
• UgliStripNode.cxx