PRDHandle_MDT Class Reference

#include <PRDHandle_MDT.h>

Inheritance diagram for PRDHandle_MDT:

Collaboration diagram for PRDHandle_MDT:

Public Member Functions
	PRDHandle_MDT (PRDCollHandle_MDT *, const Muon::MdtPrepData *)
virtual	~PRDHandle_MDT ()
void	buildShapes (SoNode *&shape_simple, SoNode *&shape_detailed)
int	regionIndex ()
QStringList	clicked () const
bool	highLight () const
bool	transformUsesSurfacePositionOnly () const
const Muon::MdtPrepData *	driftCircle () const
const Trk::PrepRawData *	getPRD () const
bool	masked () const
QString	driftCircleStatus () const
int	ADC () const
bool	inMuonChamber () const
GeoPVConstLink	parentMuonChamberPV () 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 ()
PRDCollHandleBase *	collHandle () const
virtual const Trk::PrepRawData *	getSecondPRD () const
virtual bool	isSane () const
	Returns false if the PRD is not safe to draw. More...

Private Attributes
const Muon::MdtPrepData *	m_driftcircle
Imp *	m_d
PRDSysCommonData *	m_common
bool	m_visible

Detailed Description

Definition at line 12 of file PRDHandle_MDT.h.

Constructor & Destructor Documentation

PRDHandle_MDT()

PRDHandle_MDT::PRDHandle_MDT	(	PRDCollHandle_MDT *	collhandle,
		const Muon::MdtPrepData *	dc
	)

Definition at line 19 of file PRDHandle_MDT.cxx.

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

~PRDHandle_MDT()

virtual PRDHandle_MDT::~PRDHandle_MDT ( )

inlinevirtual

Definition at line 16 of file PRDHandle_MDT.h.

{};

Member Function Documentation

ADC()

int PRDHandle_MDT::ADC ( ) const

inline

Definition at line 49 of file PRDHandle_MDT.h.

{
  return m_driftcircle->adc();
}

buildShapes()

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

virtual

Implements PRDHandleBase.

Definition at line 25 of file PRDHandle_MDT.cxx.

{
  const Trk::CylinderBounds* ccbo = dynamic_cast<const Trk::CylinderBounds*>(&(m_driftcircle->detectorElement()->surface(m_driftcircle->identify()).bounds()));
  if (not ccbo) return;
 
  double radius = m_driftcircle->localPosition()[0];
  if (radius<0.15)
    radius = 0.0;//radius is so small it is better to collapse to line
  double halflength(ccbo->halflengthZ());
  PRDCollHandle_MDT::PROJECTION projection(static_cast<PRDCollHandle_MDT*>(collHandle())->projection());
  if (projection!=PRDCollHandle_MDT::NONE) {
    const double epsilon = 0.15;
    if (projection==PRDCollHandle_MDT::TOCHAMBERS) {
      double distanceToFirstEndPlane, distanceToSecondEndPlane;
      Amg::Transform3D transform(getTransform_CLHEP());
      Amg::Vector3D tubedir(0.0,0.0,1.0);
      MuonChamberProjectionHelper::applyTransformToVector( transform, tubedir);
      if (common()->muonChamberProjectionHelper()
        ->getDistancesToMDTChamberWallsAlongLine( parentMuonChamberPV(),
        transform * Amg::Vector3D(0.0,0.0,0.0), tubedir,
      distanceToFirstEndPlane, distanceToSecondEndPlane, radius )) {
        if (fabs(distanceToFirstEndPlane-distanceToSecondEndPlane)>0.1*epsilon)
          std::cout <<"PRDHandle_MDT: Warning tube is not centered in chamber! We should translate the prd object appropriately (todo)."<<std::endl;
        halflength = std::max(distanceToFirstEndPlane,distanceToSecondEndPlane)+epsilon;
      } else {
        std::cout <<"PRDHandle_MDT: Couldn't retrieve data for projections to chamber end planes."<<std::endl;
        halflength += 500.0;
      }
    } else {
      assert(projection==PRDCollHandle_MDT::TOTUBES);
      //To end of tubes.
      halflength += 56.0 + epsilon;
    }
  }
 
  shape_simple = common()->nodeManager()->getShapeNode_DriftTube( halflength, 0.0 );
  //Detailed shape uses tube (unless negligible radius):
  if (radius==0.0)
    shape_detailed = shape_simple;
  else
    shape_detailed = common()->nodeManager()->getShapeNode_DriftTube( halflength, radius );
 
  if ( m_driftcircle->rdoList().size() > 1)
    VP1Msg::message("Warning: MDT has additional elements in rdoList: NOT IMPLEMENTED");
}

center()

Amg::Vector3D PRDHandleBase::center ( ) const

virtualinherited

Reimplemented in PRDHandle_SpacePoint.

Definition at line 316 of file PRDHandleBase.cxx.

{
  return positionPRD();
}

clicked()

QStringList PRDHandle_MDT::clicked ( ) const

virtual

Reimplemented from PRDHandleBase.

Definition at line 80 of file PRDHandle_MDT.cxx.

{
    
    if (!m_driftcircle)
      return QStringList("Null PRD");
    const MdtIdHelper * idhelper = VP1DetInfo::mdtIDHelper();
    if (!idhelper)
      return QStringList("Null PRD");
 
    std::ostringstream os;
    Identifier id = m_driftcircle->identify();
 
    os << "MdtPrepData with Identifier ["<<id.get_compact()  ;
    os << "] = [" << (idhelper->print_to_string(id))  ;
    os << "]\n at global position = [" << m_driftcircle->globalPosition()<<"], local position = ["<<m_driftcircle->localPosition()<<"].";
    os <<"\nDrift radius = "<<m_driftcircle->localPosition()[0]<<" CLHEP::mm. ";
    os <<"ADC: "<<m_driftcircle->adc();
    os <<", TDC: "<<m_driftcircle->tdc();
    std::string status("Status: ");
    os<<status;
        
    QStringList l(QString(os.str().c_str())+driftCircleStatus());
 
    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 Muon::MdtPrepData* PRDHandle_MDT::driftCircle ( ) const

inline

Definition at line 25 of file PRDHandle_MDT.h.

{ return m_driftcircle; }

driftCircleStatus()

QString PRDHandle_MDT::driftCircleStatus

(

)

const

Definition at line 107 of file PRDHandle_MDT.cxx.

{
  QString status;
  switch (m_driftcircle->status()){
    case  Muon::MdtStatusMasked:
      status+="Masked"; 
      break;
    case Muon::MdtStatusDriftTime:
      status+="Drift Time";
      break;
    case Muon::MdtStatusBeforeSpectrum:
      status+="Before Spectrum";
      break;    
    case Muon::MdtStatusAfterSpectrum:
      status+="After Spectrum";
      break;
    default:
      status+=QString((" Unknown with value = "+std::to_string(m_driftcircle->status())).c_str());
  }
  return status;
}

getPRD()

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

inlinevirtual

Implements PRDHandleBase.

Definition at line 26 of file PRDHandle_MDT.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;
}

highLight()

bool PRDHandle_MDT::highLight ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 54 of file PRDHandle_MDT.h.

{
  PRDCollHandle_MDT* ch = static_cast<PRDCollHandle_MDT*>(collHandle());
  return ( masked() && ch->highLightMasked() ) || ( ADC() <= ch->highLightADCBelow() );
}

inMuonChamber()

bool PRDHandle_MDT::inMuonChamber ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 32 of file PRDHandle_MDT.h.

{ return true; }

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.

masked()

bool PRDHandle_MDT::masked ( ) const

inline

Definition at line 44 of file PRDHandle_MDT.h.

{
  return m_driftcircle->status() == Muon::MdtStatusMasked;
}

parentMuonChamberPV()

GeoPVConstLink PRDHandle_MDT::parentMuonChamberPV ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 33 of file PRDHandle_MDT.h.

{ return m_driftcircle->detectorElement()->parentStationPV(); }

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_MDT::regionIndex ( )

virtual

Reimplemented from PRDHandleBase.

Definition at line 72 of file PRDHandle_MDT.cxx.

{
  //Unique for each station.
  const MuonGM::MuonStation* station = m_driftcircle->detectorElement()->parentMuonStation();
  return station->getPhiIndex()-99999*station->getEtaIndex();//hopefully unique.
}

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_MDT::transformUsesSurfacePositionOnly ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 24 of file PRDHandle_MDT.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 Muon::MdtPrepData* PRDHandle_MDT::m_driftcircle

private

Definition at line 40 of file PRDHandle_MDT.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_MDT.h
• PRDHandle_MDT.cxx