PRDHandle_Pixel Class Reference

#include <PRDHandle_Pixel.h>

Inheritance diagram for PRDHandle_Pixel:

Collaboration diagram for PRDHandle_Pixel:

Public Member Functions
	PRDHandle_Pixel (PRDCollHandle_Pixel *, const InDet::PixelCluster *)
virtual	~PRDHandle_Pixel ()
void	buildShapes (SoNode *&shape_simple, SoNode *&shape_detailed)
int	regionIndex ()
const InDet::PixelCluster *	cluster () const
const Trk::PrepRawData *	getPRD () const
bool	isBarrel () const
bool	isPositiveZ () const
bool	isSane () const
	returns false if anything seems weird about the Pixelcluster More...
void	setVisible (bool)
bool	visible () const
virtual bool	highLight () const
virtual bool	transformUsesSurfacePositionOnly () 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 QStringList	clicked () const
virtual bool	inMuonChamber () const
virtual GeoPVConstLink	parentMuonChamberPV () const
PRDCollHandleBase *	collHandle () const
virtual const Trk::PrepRawData *	getSecondPRD () const

Private Member Functions
SoGroup *	createErrorAtPixelCluster (const double &sigmaScale, const int &numPoints) const

Private Attributes
const InDet::PixelCluster *	m_cluster
Imp *	m_d
PRDSysCommonData *	m_common
bool	m_visible

Detailed Description

Definition at line 14 of file PRDHandle_Pixel.h.

Constructor & Destructor Documentation

PRDHandle_Pixel()

PRDHandle_Pixel::PRDHandle_Pixel	(	PRDCollHandle_Pixel *	collhandle,
		const InDet::PixelCluster *	clus
	)

Definition at line 23 of file PRDHandle_Pixel.cxx.

  : PRDHandleBase(static_cast<PRDCollHandleBase*>(collhandle)), m_cluster(clus)
{
    SoTransparency::initClass();
}

~PRDHandle_Pixel()

virtual PRDHandle_Pixel::~PRDHandle_Pixel ( )

inlinevirtual

Definition at line 18 of file PRDHandle_Pixel.h.

{};

Member Function Documentation

buildShapes()

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

virtual

Implements PRDHandleBase.

Definition at line 30 of file PRDHandle_Pixel.cxx.

{
  SoSeparator * errSimple  = new SoSeparator;
  SoSeparator * errDetailed  = new SoSeparator;
  
  if (!isSane()) {
    std::cerr<<"isSane check failed for PRD: "<<*m_cluster<<std::endl;
    return;
  }  
  errSimple->addChild(common()->nodeManager()->getShapeNode_Point());
 
  Amg::Vector2D localpos = m_cluster->localPosition();
  Amg::Vector2D localposHIT = m_cluster->detectorElement()->rawLocalPositionOfCell( m_cluster->identify() ); // Lorentz shift is not corrected because this class cannot use SiLorentzAngleTool
  SoTranslation * localtrans0 = new SoTranslation;
 
  double xdiff = localposHIT[Trk::locX]-localpos[Trk::locX];
  double ydiff = localposHIT[Trk::locY]-localpos[Trk::locY];
  
  if (xdiff!=xdiff || ydiff!=ydiff) {
    std::cerr<<"NaN is in local pos calc"<<std::endl;
    if (xdiff!=xdiff) std::cerr<<"X diff"<<std::endl;
    if (ydiff!=ydiff) std::cerr<<"Y diff"<<std::endl;
    std::cerr<<"localposHIT:"<<localposHIT<<std::endl;
    std::cerr<<"localpos:"<<localpos<<std::endl;
    std::cerr<<"localposHIT[Trk::locX]-localpos[Trk::locX]:"<<xdiff<<std::endl;
    std::cerr<<"localposHIT[Trk::locY]-localpos[Trk::locY]:"<<ydiff<<std::endl;
    localtrans0->unref();
    return;
  } else {
    localtrans0->translation.setValue(localposHIT[Trk::locX]-localpos[Trk::locX],localposHIT[Trk::locY]-localpos[Trk::locY],0);
  }
 
  const std::vector<Identifier> rdolist = m_cluster->rdoList();
  if (rdolist.size() == 1 || !collHandle()->drawRDOs())
  {
    errDetailed->addChild(localtrans0);
    //FIXME: dont hardcode thickness
    errDetailed->addChild(common()->nodeManager()->getShapeNode_Strip(m_cluster->detectorElement()->etaPitch(),m_cluster->detectorElement()->etaPitch(),m_cluster->detectorElement()->thickness()/10.0));
 
    //Translate back so errDetailed is left sane (eg. when drawing errors later)
    SoTranslation * localtransBack = new SoTranslation;
    localtransBack->translation.setValue(-(localtrans0->translation.getValue()));
    errDetailed->addChild(localtransBack);
  } else {
    SoSeparator * rdos = new SoSeparator;
    rdos->addChild(localtrans0);
 
    rdos->addChild(common()->nodeManager()->getShapeNode_Strip(m_cluster->detectorElement()->etaPitch(),m_cluster->detectorElement()->phiPitch( m_cluster->localPosition() ),m_cluster->detectorElement()->thickness()/10.0));
 
    SoTransparency * transparent = new SoTransparency;
    transparent->transparency.setValue(0.5);
    rdos->addChild( transparent );
    Amg::Vector2D localposOLD = localposHIT;
    for (const Identifier& rdo_id : rdolist)
    {
      if (rdo_id == m_cluster->identify() )
        continue;
 
      const Amg::Vector2D localposRDO = m_cluster->detectorElement()->rawLocalPositionOfCell(rdo_id); // Lorentz shift is not corrected because this class cannot use SiLorentzAngleTool
 
      SoTranslation * localtrans = new SoTranslation;
      localtrans->translation.setValue(localposRDO[Trk::locX]-localposOLD[Trk::locX],
                                       localposRDO[Trk::locY]-localposOLD[Trk::locY],
                                       0);
      rdos->addChild(localtrans);
 
      //FIXME: dont hardcode thickness
      rdos->addChild(common()->nodeManager()->getShapeNode_Strip(m_cluster->detectorElement()->etaPitch(),
                     m_cluster->detectorElement()->phiPitch( m_cluster->localPosition() ),
                     m_cluster->detectorElement()->thickness()/10.0));
 
      localposOLD = localposRDO;
    }
    errDetailed->addChild(rdos);
  }
 
  const double settingsSIGMASCALE=3.0;
  if (collHandle()->drawErrors())
  {
    SoGroup * gr = createErrorAtPixelCluster( settingsSIGMASCALE, 36);
    if (gr->getNumChildren()==0) {
      gr->unref();
      std::cerr<<"Something seems to have gone wrong: no error added. Dumping PRD: "<<*m_cluster<<std::endl;
    } else {
      errSimple->addChild( gr );
      errDetailed->addChild( gr );
    }
  }
 
  //FIXME: hardcoded colour  
  SoMaterial * mat = new SoMaterial;
  mat->diffuseColor.setValue(1.0,0,0);
  errDetailed->addChild(mat);
  errDetailed->addChild(common()->nodeManager()->getShapeNode_Cross(0.1));
 
  shape_simple = errSimple;
  shape_detailed = errDetailed;
}

center()

Amg::Vector3D PRDHandleBase::center ( ) const

virtualinherited

Reimplemented in PRDHandle_SpacePoint.

Definition at line 316 of file PRDHandleBase.cxx.

{
  return positionPRD();
}

clicked()

QStringList PRDHandleBase::clicked ( ) const

virtualinherited

Reimplemented in PRDHandle_TRT, PRDHandle_CSC, PRDHandle_CSC_Strip, PRDHandle_MDT, PRDHandle_MM, and PRDHandle_RPC.

Definition at line 276 of file PRDHandleBase.cxx.

{
  const Trk::PrepRawData * prd = getPRD();
  if (!prd)
    return QStringList("Null PRD");
  std::ostringstream os;
  os << *(prd);
  return QString(os.str().c_str()).split("\n");
}

cluster()

const InDet::PixelCluster* PRDHandle_Pixel::cluster ( ) const

inline

Definition at line 23 of file PRDHandle_Pixel.h.

{ return m_cluster; }

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; }

createErrorAtPixelCluster()

SoGroup * PRDHandle_Pixel::createErrorAtPixelCluster ( const double &  sigmaScale, const int &  numPoints  ) const

private

Definition at line 140 of file PRDHandle_Pixel.cxx.

{
  SoGroup * ellipseGrp = new SoSeparator;
 
  //FIXME: Try to automatically use the pixel height as in PRDHandle_Pixel.cxx (this is just copy paste)
  // change the last true to false above to get rid of one of the error ellipses
//  VP1ErrorUtils::addRotatedErrorEllipse(ellipseGrp, m_cluster->localErrorMatrix().covariance(), sigmaScale, numPoints,m_cluster->detectorElement()->thickness()/20.0 + 0.001, true);
  VP1ErrorUtils::addRotatedErrorEllipse(ellipseGrp, m_cluster->localCovariance(), sigmaScale, numPoints,m_cluster->detectorElement()->thickness()/20.0 + 0.001, true);
 
  return ellipseGrp;
}

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;
}

getPRD()

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

inlinevirtual

Implements PRDHandleBase.

Definition at line 24 of file PRDHandle_Pixel.h.

{ return m_cluster; }

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()

virtual bool PRDHandleBase::highLight ( ) const

inlinevirtualinherited

Reimplemented in PRDHandle_MDT, and PRDHandle_TRT.

Definition at line 45 of file PRDHandleBase.h.

{ return false; }//Return true to get hit highlighted (for e.g. a TRT hits with HT).

inMuonChamber()

virtual bool PRDHandleBase::inMuonChamber ( ) const

inlinevirtualinherited

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

Definition at line 67 of file PRDHandleBase.h.

{ return false; }

isBarrel()

bool PRDHandle_Pixel::isBarrel ( ) const

inline

Definition at line 25 of file PRDHandle_Pixel.h.

{ return m_cluster->detectorElement()->isBarrel(); }

isPositiveZ()

bool PRDHandle_Pixel::isPositiveZ ( ) const

inline

Definition at line 26 of file PRDHandle_Pixel.h.

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

isSane()

bool PRDHandle_Pixel::isSane ( ) const

virtual

returns false if anything seems weird about the Pixelcluster

Reimplemented from PRDHandleBase.

Definition at line 153 of file PRDHandle_Pixel.cxx.

                                   {
  // check for NaNs
  bool isSane=true;
  const InDet::PixelCluster * prd =cluster();
  if (!prd) return false;
  if (prd->localPosition()[0]!=prd->localPosition()[0] || prd->localPosition()[1]!=prd->localPosition()[1]){
    isSane=false;
    std::cerr<<"For pixel cluster with Id="<<prd->identify()<<", localposition is not sane:"<<(prd->localPosition())<<std::endl;
  }
  bool etaNan = std::isnan(m_cluster->detectorElement()->etaPitch());
  bool phiNan = std::isnan(m_cluster->detectorElement()->phiPitch( m_cluster->localPosition()));
  bool thickNan = std::isnan(m_cluster->detectorElement()->thickness());
 
  if ( etaNan || phiNan || thickNan){
    std::cerr<<"For pixel cluster with Id="<<prd->identify()<<", NaN is in detectorElement pitches "<<std::endl;
    if (etaNan) std::cerr<<" - eta (which is returning "<<m_cluster->detectorElement()->etaPitch()<<")"<<std::endl;
    if (phiNan) std::cerr<<" - phi (which is returning "<<m_cluster->detectorElement()->phiPitch(m_cluster->localPosition())<<")"<<std::endl;
    if (thickNan) std::cerr<<" - thickness (which is returning "<<m_cluster->detectorElement()->thickness()<<")"<<std::endl;
  }
  
  return isSane;
}

parentMuonChamberPV()

virtual GeoPVConstLink PRDHandleBase::parentMuonChamberPV ( ) const

inlinevirtualinherited

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

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

virtual

Reimplemented from PRDHandleBase.

Definition at line 130 of file PRDHandle_Pixel.cxx.

                                 { //instead of identifier juggling, we simply discretize by center coordinate.
  static const double l=200.0;//CLHEP::mm
  return       static_cast<int>(m_cluster->detectorElement()->center().z()/l)
         +1000*static_cast<int>(m_cluster->detectorElement()->center().y()/l)
      +1000000*static_cast<int>(m_cluster->detectorElement()->center().x()/l);
  //Fixme: Use identifiers instead for more intuitive regions.
}

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()

virtual bool PRDHandleBase::transformUsesSurfacePositionOnly ( ) const

inlinevirtualinherited

Reimplemented in PRDHandle_MDT, and PRDHandle_TRT.

Definition at line 47 of file PRDHandleBase.h.

{ return false; }//Override and return true if should use just the surface position.

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_cluster

const InDet::PixelCluster* PRDHandle_Pixel::m_cluster

private

Definition at line 30 of file PRDHandle_Pixel.h.

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_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_Pixel.h
• PRDHandle_Pixel.cxx