PRDHandle_TRT Class Reference

#include <PRDHandle_TRT.h>

Inheritance diagram for PRDHandle_TRT:

Collaboration diagram for PRDHandle_TRT:

Public Member Functions
	PRDHandle_TRT (PRDCollHandle_TRT *, const InDet::TRT_DriftCircle *)
virtual	~PRDHandle_TRT ()
void	buildShapes (SoNode *&shape_simple, SoNode *&shape_detailed)
int	regionIndex ()
bool	highLight () const
bool	transformUsesSurfacePositionOnly () const
const InDet::TRT_DriftCircle *	driftCircle () const
const Trk::PrepRawData *	getPRD () const
bool	isBarrel () const
bool	isPositiveZ () const
bool	highLevel () const
virtual QStringList	clicked () const
void	setVisible (bool)
bool	visible () const
SoSeparator *	sepSimple () const
SoSeparator *	sepDetailed () const
PRDSysCommonData *	common () const
virtual Amg::Vector3D	center () const
Amg::Vector3D	positionPRD () const
Amg::Vector3D	positionSecondPRD () const
Amg::Transform3D	getTransform_CLHEP () const
virtual SoTransform *	createTransform () const
void	update3DObjects ()
void	updateMaterial ()
virtual bool	inMuonChamber () const
virtual GeoPVConstLink	parentMuonChamberPV () const
PRDCollHandleBase *	collHandle () const
virtual const Trk::PrepRawData *	getSecondPRD () const
virtual bool	isSane () const
	Returns false if the PRD is not safe to draw.

Private Attributes
const InDet::TRT_DriftCircle *	m_driftcircle
Imp *	m_d
PRDSysCommonData *	m_common
bool	m_visible

Detailed Description

Definition at line 12 of file PRDHandle_TRT.h.

Constructor & Destructor Documentation

PRDHandle_TRT()

PRDHandle_TRT::PRDHandle_TRT	(	PRDCollHandle_TRT *	collhandle,
		const InDet::TRT_DriftCircle *	dc )

Definition at line 24 of file PRDHandle_TRT.cxx.

  : PRDHandleBase(static_cast<PRDCollHandleBase*>(collhandle)), m_driftcircle(dc)
{
}

~PRDHandle_TRT()

virtual PRDHandle_TRT::~PRDHandle_TRT ( )

inlinevirtual

Definition at line 16 of file PRDHandle_TRT.h.

{};

Member Function Documentation

buildShapes()

void PRDHandle_TRT::buildShapes ( SoNode *& shape_simple, SoNode *& shape_detailed )

virtual

Implements PRDHandleBase.

Definition at line 30 of file PRDHandle_TRT.cxx.

{
  const Trk::CylinderBounds* ccbo = dynamic_cast<const Trk::CylinderBounds*>(&(m_driftcircle->detectorElement()->surface(m_driftcircle->identify()).bounds()));
  if (not ccbo) return;
  const double radius = m_driftcircle->localPosition()[0];
  const double halflength = ccbo->halflengthZ();
  const double mintuberadius = 0.1;
 
  if ( m_driftcircle->rdoList().size() > 1)
    VP1Msg::message("Warning: TRT has additional elements in rdoList: NOT IMPLEMENTED");
 
  const PRDCollHandle_TRT* trtcollhandle(static_cast<PRDCollHandle_TRT*>(collHandle()));
 
  if (trtcollhandle->project()) {

    //  Projections:  //
 
    const InDetProjFlags::InDetProjPartsFlags projparts(trtcollhandle->appropriateProjections());
    bool inner,outer;
    if (isBarrel()) {
      inner = (projparts&InDetProjFlags::BarrelCentral);
      outer = (projparts&(isPositiveZ()?InDetProjFlags::BarrelPositive:InDetProjFlags::BarrelNegative));
    } else {
      const bool posz(isPositiveZ());
      if ( (projparts&InDetProjFlags::TRT_EndCapZToRCentral)
       || (projparts&(posz?InDetProjFlags::TRT_EndCapZToRPositive:InDetProjFlags::TRT_EndCapZToRNegative)) ) {
    //Special TRT EndCap Z->R projections.
 
    //Find center of tube:
    Amg::Transform3D prdtransform(getTransform_CLHEP());
    Amg::Vector3D center(prdtransform*Amg::Vector3D(0,0,0));
 
    // --> Get points for special projections
    const double eps = InDetProjParams::trt_data_disttosurface_epsilon();
    std::vector<Amg::Vector3D > points;
 
    if (posz) {
      if ( projparts & InDetProjFlags::TRT_EndCapZToRPositive ) {
        Amg::Vector3D p = center;
        InDetProjHelper::transformECPointToZPlane_specialZtoR(p,
                                  InDetProjParams::trt_barrel_posneg_z()-0.5* InDetProjParams::surfacethickness() - eps,
                                  InDetProjParams::trt_endcap_zasr_innerradius(),
                                  InDetProjParams::trt_endcap_zasr_endcapz_begin(),
                                  InDetProjParams::trt_endcap_zasr_squeezefact());
        points.push_back(p);
      }
      if ( projparts & InDetProjFlags::TRT_EndCapZToRCentral ) {
        Amg::Vector3D p = center;
        //Fixme: Use same parameters here as in InDetProjHelper!!
        InDetProjHelper::transformECPointToZPlane_specialZtoR(p,
                                  0.5* InDetProjParams::surfacethickness() + eps,
                                  InDetProjParams::trt_endcap_zasr_innerradius(),
                                  InDetProjParams::trt_endcap_zasr_endcapz_begin(),
                                  InDetProjParams::trt_endcap_zasr_squeezefact() );
        points.push_back(p);
      }
    } else {
      if ( projparts & InDetProjFlags::TRT_EndCapZToRNegative ) {
        Amg::Vector3D p = center;
        InDetProjHelper::transformECPointToZPlane_specialZtoR(p,
                                  - InDetProjParams::trt_barrel_posneg_z()+0.5* InDetProjParams::surfacethickness() + eps,
                                  InDetProjParams::trt_endcap_zasr_innerradius(),
                                  InDetProjParams::trt_endcap_zasr_endcapz_begin(),
                                  InDetProjParams::trt_endcap_zasr_squeezefact());
        points.push_back(p);
      }
      if ( projparts & InDetProjFlags::TRT_EndCapZToRCentral ) {
        Amg::Vector3D p = center;
        InDetProjHelper::transformECPointToZPlane_specialZtoR(p,
                                  -0.5* InDetProjParams::surfacethickness() - eps,
                                  InDetProjParams::trt_endcap_zasr_innerradius(),
                                  InDetProjParams::trt_endcap_zasr_endcapz_begin(),
                                  InDetProjParams::trt_endcap_zasr_squeezefact());
        points.push_back(p);
      }
    }
    // --> Build scene objects for special projections
 
    const unsigned npoints = points.size();
    if (npoints<1||npoints>2) {
      collHandle()->message("PRDHandle_TRT::buildShapes ERROR: Unexpected number of points in special endcap projections!!");
    } else {
      SoSeparator * sepsimple = new SoSeparator;
      SoTransform * transf1(VP1LinAlgUtils::toSoTransform((prdtransform.inverse())
                                  * Amg::Translation3D(points.at(0).x(),points.at(0).y(),points.at(0).z())));
      SoNode* point = common()->nodeManager()->getShapeNode_Point();
      sepsimple->addChild(transf1);
      sepsimple->addChild(point);
      SoTranslation * transf2(0);
      if (npoints==2) {
        transf2 = new SoTranslation;
        transf2->translation.setValue(points[1].x()-points[0].x(),points[1].y()-points[0].y(),points[1].z()-points[0].z());
        sepsimple->addChild(transf2);
        sepsimple->addChild(point);
      }
      shape_simple = sepsimple;
 
      if (radius<mintuberadius) {
        shape_detailed = shape_simple;
        return;
      }
      SoSeparator * sepdetailed = new SoSeparator;
      SoNode* disc = common()->nodeManager()->getShapeNode_DriftDisc(radius);
      sepdetailed->addChild(transf1);
      sepdetailed->addChild(disc);
      if (npoints==2) {
        sepdetailed->addChild(transf2);
        sepdetailed->addChild(disc);
      }
      shape_detailed = sepdetailed;
      return;
    }
    inner = false;outer = false;
    //end special projections
      } else {
    inner = (projparts&(posz?InDetProjFlags::EndCapInnerPositive:InDetProjFlags::EndCapInnerNegative));
    outer = (projparts&(posz?InDetProjFlags::EndCapOuterPositive:InDetProjFlags::EndCapOuterNegative));
      }
    }//end endcap
    if (!inner&&!outer) {
      collHandle()->message("PRDHandle_TRT::buildShapes ERROR: Inconsistent projection settings! Won't project!");
    } else {
      //Perform projections
      shape_simple = common()->nodeManager()->getShapeNode_ProjectedDriftTube( halflength, 0.0, inner, outer );
      if (radius<mintuberadius)
    shape_detailed = shape_simple;
      else
    shape_detailed = common()->nodeManager()->getShapeNode_ProjectedDriftTube( halflength, radius, inner, outer );
      return;
    }
  }//endif projections

  //  No projections:  //
 
  shape_simple = common()->nodeManager()->getShapeNode_DriftTube( halflength, 0.0 );
  if (radius<mintuberadius)
    shape_detailed = shape_simple;
  else
    shape_detailed = common()->nodeManager()->getShapeNode_DriftTube( halflength, radius );
 
}

center()

Amg::Vector3D PRDHandleBase::center ( ) const

virtualinherited

Reimplemented in PRDHandle_SpacePoint.

Definition at line 316 of file PRDHandleBase.cxx.

{
  return positionPRD();
}

clicked()

QStringList PRDHandle_TRT::clicked ( ) const

virtual

Reimplemented from PRDHandleBase.

Definition at line 188 of file PRDHandle_TRT.cxx.

{
  QStringList l;
  l << PRDHandleBase::clicked();
  const unsigned int data = m_driftcircle->getWord();
 
  QString s;
  int i;
  for (i = 26; i >=0; i--)
  {
    if (i!=26&&(i+1)%9==0)
      s+=" ";
    s += ( data & (0x1<<i) ? "1" : "0");
    if ((i+1)%9==0)
      s+=" ";
  }
  l << "Raw data bits: "+s;
 
 
  std::ostringstream os;
  os << std::hex<<data;
  l << QString(os.str().c_str());
 
 
//   l << "Drift radius: "+QString::number(m_driftcircle->localPosition()[0]/CLHEP::mm)+" CLHEP::mm";
//   const Trk::CylinderBounds* ccbo = dynamic_cast<const Trk::CylinderBounds*>(&(m_driftcircle->detectorElement()->surface(m_driftcircle->identify()).bounds()));
//   if (ccbo)
//     l << "Straw length: "+QString::number(ccbo->halflengthZ()*2/CLHEP::cm)+" CLHEP::cm";
 
  return l;
}

collHandle()

PRDCollHandleBase * PRDHandleBase::collHandle ( ) const

inherited

Definition at line 228 of file PRDHandleBase.cxx.

{
  return m_d->colhandle;
}

common()

PRDSysCommonData * PRDHandleBase::common ( ) const

inlineinherited

Definition at line 53 of file PRDHandleBase.h.

{ return m_common; }

createTransform()

SoTransform * PRDHandleBase::createTransform ( ) const

virtualinherited

Reimplemented in PRDHandle_sTGC, and PRDHandle_TGC.

Definition at line 371 of file PRDHandleBase.cxx.

{  
  //NB: Some code duplicated here and in the getTransform_CLHEP() method.
  const Trk::PrepRawData * prd = getPRD();
  if (!prd) {
    std::cerr<<"PRDHandleBase::createTransform() No prd!"<<std::endl;
    return new SoTransform;
  }
  const Trk::Surface& theSurface = prd->detectorElement()->surface(prd->identify());
  SoTransform * theHitTransform = VP1LinAlgUtils::toSoTransform(theSurface.transform());
 
  Amg::Vector3D theHitGPos;
  if ( transformUsesSurfacePositionOnly() ) {
    // for tubes, should use position of center of tube (if drawing full tube)
    theHitGPos= Amg::Vector3D (theSurface.center());
  } else {
    // for strips, clusters or short tubes, use position of hit.
    theHitGPos = theSurface.localToGlobal(prd->localPosition());
  }
  if ((theHitGPos)[0]!=(theHitGPos)[0] || (theHitGPos)[1]!=(theHitGPos)[1] || (theHitGPos)[2]!=(theHitGPos)[2]){
    std::cerr<<"PRDHandleBase::createTransform() NaN in globalposition"<<std::endl;
  }
  theHitTransform->translation.setValue((theHitGPos)[0], (theHitGPos)[1], (theHitGPos)[2]);
  return theHitTransform;
}

driftCircle()

const InDet::TRT_DriftCircle * PRDHandle_TRT::driftCircle ( ) const

inline

Definition at line 24 of file PRDHandle_TRT.h.

{ return m_driftcircle; }

getPRD()

const Trk::PrepRawData * PRDHandle_TRT::getPRD ( ) const

inlinevirtual

Implements PRDHandleBase.

Definition at line 25 of file PRDHandle_TRT.h.

{ return m_driftcircle; }

getSecondPRD()

virtual const Trk::PrepRawData * PRDHandleBase::getSecondPRD ( ) const

inlinevirtualinherited

Reimplemented in PRDHandle_SpacePoint.

Definition at line 73 of file PRDHandleBase.h.

{ return 0; }//For SCT spacepoints

getTransform_CLHEP()

Amg::Transform3D PRDHandleBase::getTransform_CLHEP ( ) const

inherited

Definition at line 322 of file PRDHandleBase.cxx.

{  
  if (m_d->sepDetailed&&m_d->sepDetailed->getNumChildren()>0) {
    //For efficiency/consistency we simply copy the transform from the
    //inventor transformation (only loose slight precision in the
    //double->float->double conversions):
    SoTransform * transform = static_cast<SoTransform*>(m_d->sepDetailed->getChild(0));
    float tx, ty, tz;
    transform->translation.getValue().getValue(tx,ty,tz);
    SbVec3f so_rotaxis;
    float so_rotangle;
    transform->rotation.getValue().getValue(so_rotaxis, so_rotangle);
    float rx, ry, rz;
    so_rotaxis.getValue(rx,ry,rz);
    Amg::Vector3D axis(rx,ry,rz);
//    Amg::AngleAxis3D angleAxis = Amg::getRotation3DfromAngleAxis( so_rotangle, axis );
//    Amg::Rotation3D rot;
//    rot = angleAxis;
    Amg::Rotation3D rot = Amg::getRotation3DfromAngleAxis( so_rotangle, axis );
    return Amg::Translation3D(tx,ty,tz) * rot;
  }
  std::cout<<"VP1 PRDHandleBase Warning: getTransform_CLHEP called at inopportune moment!!"<<std::endl;
 
  //NB: Some code duplicated here and in the createTransform() method.
 
  const Trk::PrepRawData * prd = getPRD();
  if (!prd) {
    return Amg::Transform3D::Identity();
  }
  const Trk::Surface& theSurface = prd->detectorElement()->surface(prd->identify());
  Amg::Vector3D theHitGPos;
  if ( transformUsesSurfacePositionOnly() ) {
    // for tubes, should use position of center of tube (if drawing full tube)
    theHitGPos= Amg::Vector3D(theSurface.center());
  } else {
    // for clusters or short tubes, use position of hit.
    theHitGPos = theSurface.localToGlobal(prd->localPosition()); //theSurface 'new s' a Vector3D and returns pointer
  }
 
//  CLHEP::Hep3Vector t(theSurface.transform().getTranslation());
//  return HepGeom::Transform3D( HepGeom::Translate3D(theHitGPos->x()-t.x(),theHitGPos->y()-t.y(),theHitGPos->z()-t.z()) * (theSurface.transform()) );
  Amg::Vector3D t;
  t = Amg::getTranslationVectorFromTransform(theSurface.transform());
  Amg::Translation3D transl = Amg::Translation3D(theHitGPos.x()-t.x(), theHitGPos.y()-t.y(), theHitGPos.z()-t.z());
  Amg::Transform3D transf = transl * (theSurface.transform());
  return transf;
}

highLevel()

bool PRDHandle_TRT::highLevel ( ) const

inline

Definition at line 29 of file PRDHandle_TRT.h.

{ return m_driftcircle->highLevel(); }

highLight()

bool PRDHandle_TRT::highLight ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 43 of file PRDHandle_TRT.h.

{
  return highLevel() && static_cast<PRDCollHandle_TRT*>(collHandle())->highLightHighThreshold();
}

inMuonChamber()

virtual bool PRDHandleBase::inMuonChamber ( ) const

inlinevirtualinherited

Reimplemented in PRDHandle_CSC, PRDHandle_CSC_Strip, PRDHandle_MDT, PRDHandle_MM, PRDHandle_RPC, PRDHandle_sTGC, and PRDHandle_TGC.

Definition at line 67 of file PRDHandleBase.h.

{ return false; }

isBarrel()

bool PRDHandle_TRT::isBarrel ( ) const

inline

Definition at line 27 of file PRDHandle_TRT.h.

{ return m_driftcircle->detectorElement()->type()==InDetDD::TRT_BaseElement::BARREL; }

isPositiveZ()

bool PRDHandle_TRT::isPositiveZ ( ) const

inline

Definition at line 28 of file PRDHandle_TRT.h.

{ return m_driftcircle->detectorElement()->center().z() > 0.0; }

isSane()

virtual bool PRDHandleBase::isSane ( ) const

inlinevirtualinherited

Returns false if the PRD is not safe to draw.

Reimplemented in PRDHandle_Pixel.

Definition at line 75 of file PRDHandleBase.h.

parentMuonChamberPV()

virtual GeoPVConstLink PRDHandleBase::parentMuonChamberPV ( ) const

inlinevirtualinherited

Reimplemented in PRDHandle_CSC, PRDHandle_CSC_Strip, PRDHandle_MDT, PRDHandle_MM, PRDHandle_RPC, PRDHandle_sTGC, and PRDHandle_TGC.

Definition at line 68 of file PRDHandleBase.h.

{ return GeoPVConstLink(); }

positionPRD()

Amg::Vector3D PRDHandleBase::positionPRD ( ) const

inherited

Definition at line 303 of file PRDHandleBase.cxx.

{
  return Imp::positionPRD(getPRD());
}

positionSecondPRD()

Amg::Vector3D PRDHandleBase::positionSecondPRD ( ) const

inherited

Definition at line 309 of file PRDHandleBase.cxx.

{
  return Imp::positionPRD(getSecondPRD());
}

regionIndex()

int PRDHandle_TRT::regionIndex ( )

virtual

Reimplemented from PRDHandleBase.

Definition at line 177 of file PRDHandle_TRT.cxx.

{
  const TRT_ID * id_helper = VP1DetInfo::trtIDHelper();
  if (!id_helper)
    return 0;
  Identifier id = m_driftcircle->identify();
  return id_helper->barrel_ec(id)*10000+id_helper->layer_or_wheel(id)*100+(id_helper->phi_module(id)%4);
}

sepDetailed()

SoSeparator * PRDHandleBase::sepDetailed ( ) const

inherited

Definition at line 60 of file PRDHandleBase.cxx.

{
  return m_d->sepDetailed;
}

sepSimple()

SoSeparator * PRDHandleBase::sepSimple ( ) const

inherited

Definition at line 54 of file PRDHandleBase.cxx.

{
  return m_d->sepSimple;
}

setVisible()

void PRDHandleBase::setVisible ( bool vis )

inherited

Definition at line 234 of file PRDHandleBase.cxx.

{
  if (vis==m_visible)
    return;
  //std::cout<<"Changing visible status from "<<m_visible<<" to "<<vis<<" for: "<<*getPRD()<<std::endl;
  m_visible=vis;
  if (vis) {
    m_d->colhandle->incrementNShownHandles();
    if (inMuonChamber())
      common()->touchedMuonChamberHelper()->incrementNumberOfObjectsForPV(parentMuonChamberPV());
    if (!m_d->sepSimple||!m_d->sepDetailed)
      m_d->rebuild3DObjects(this);//The call to rebuild also fixes attached state.
    else
      m_d->attach3DObjects(this);
  } else {
    m_d->colhandle->decrementNShownHandles();
    if (inMuonChamber())
      common()->touchedMuonChamberHelper()->decrementNumberOfObjectsForPV(parentMuonChamberPV());
    m_d->detach3DObjects();
  }
}

transformUsesSurfacePositionOnly()

bool PRDHandle_TRT::transformUsesSurfacePositionOnly ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 22 of file PRDHandle_TRT.h.

{ return true; }

update3DObjects()

void PRDHandleBase::update3DObjects ( )

inherited

Definition at line 257 of file PRDHandleBase.cxx.

                                    {
  //Fixme: If selected we really need to redo selection updates!!!
  if (m_visible) {
    m_d->rebuild3DObjects(this);
  } else {
    //Simply clear the present 3D objects. They will only be recreated if/when the prd becomes visible again.
    if (m_d->sepDetailed) {
      m_d->sepDetailed->unref();
      m_d->sepDetailed = 0;
    }
    if (m_d->sepSimple) {
      m_d->sepSimple->unref();
      m_d->sepSimple = 0;
    }
  }
 
}

updateMaterial()

void PRDHandleBase::updateMaterial ( )

inherited

Definition at line 160 of file PRDHandleBase.cxx.

{
  //We need to change the current material. This means we need to
  //clear the cache. Only if we are visible do we need to do a
  //detach->clear material->attach cycle to trigger attachment under new material.
 
  if (!m_d->currentmaterial)
    return;//We have no material cached and is thus not attached either.
  if (m_visible) {
    //See if the material changed. If it did, detach, update the material, attach again.
    SoMaterial * newmat = m_d->determineMaterial(this);
    newmat->ref();
    if (newmat!=m_d->currentmaterial) {
      m_d->detach3DObjects();
      m_d->currentmaterial->unref();
      m_d->currentmaterial = newmat;
      m_d->attach3DObjects(this);
    } else {
      newmat->unref();
    }
  } else {
    //Just clear material.
    m_d->currentmaterial->unref();
    m_d->currentmaterial = 0;
  }
}

visible()

bool PRDHandleBase::visible ( ) const

inlineinherited

Definition at line 42 of file PRDHandleBase.h.

{ return m_visible; }

Member Data Documentation

m_common

PRDSysCommonData* PRDHandleBase::m_common

privateinherited

Definition at line 84 of file PRDHandleBase.h.

m_d

Imp* PRDHandleBase::m_d

privateinherited

Definition at line 83 of file PRDHandleBase.h.

m_driftcircle

const InDet::TRT_DriftCircle* PRDHandle_TRT::m_driftcircle

private

Definition at line 35 of file PRDHandle_TRT.h.

m_visible

bool PRDHandleBase::m_visible

privateinherited

Definition at line 86 of file PRDHandleBase.h.

---

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

• PRDHandle_TRT.h
• PRDHandle_TRT.cxx