PRDHandle_TGC Class Reference

#include <PRDHandle_TGC.h>

Inheritance diagram for PRDHandle_TGC:

Collaboration diagram for PRDHandle_TGC:

Public Member Functions
	PRDHandle_TGC (PRDCollHandle_TGC *, const Muon::TgcPrepData *)
virtual	~PRDHandle_TGC ()
SoTransform *	createTransform () const
void	buildShapes (SoNode *&shape_simple, SoNode *&shape_detailed)
int	regionIndex ()
const Muon::TgcPrepData *	tgc () const
const Trk::PrepRawData *	getPRD () const
bool	inMuonChamber () const
GeoPVConstLink	parentMuonChamberPV () const
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
void	update3DObjects ()
void	updateMaterial ()
virtual QStringList	clicked () const
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::TgcPrepData *	m_tgc
Imp *	m_d
PRDSysCommonData *	m_common
bool	m_visible

Detailed Description

Definition at line 12 of file PRDHandle_TGC.h.

Constructor & Destructor Documentation

PRDHandle_TGC()

PRDHandle_TGC::PRDHandle_TGC	(	PRDCollHandle_TGC *	collhandle,
		const Muon::TgcPrepData *	tgc
	)

Definition at line 24 of file PRDHandle_TGC.cxx.

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

~PRDHandle_TGC()

virtual PRDHandle_TGC::~PRDHandle_TGC ( )

inlinevirtual

Definition at line 16 of file PRDHandle_TGC.h.

{};

Member Function Documentation

buildShapes()

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

virtual

Implements PRDHandleBase.

Definition at line 97 of file PRDHandle_TGC.cxx.

{
  const TgcIdHelper * idhelper = VP1DetInfo::tgcIDHelper();
  if (!idhelper)
    return;
 
  Identifier id = m_tgc->identify();
 
  int plane = idhelper->gasGap( id );
  int strip = idhelper->channel( id );
  int isStrip = idhelper->isStrip( id );
 
  double striplength =0.0, stripWidth = 0.0;
 
  if (isStrip){
    striplength = m_tgc->detectorElement()->stripLength();
    stripWidth = m_tgc->detectorElement()->stripWidth(plane, strip);
  } else {    
    striplength = m_tgc->detectorElement()->gangShortWidth(plane, strip);
    stripWidth = m_tgc->detectorElement()->gangRadialLength(plane, strip);
  }
 
  if (static_cast<PRDCollHandle_TGC*>(collHandle())->project())
    striplength += 300.0;//Fixme: Rough extension for now
 
  shape_simple = common()->nodeManager()->getShapeNode_Strip(striplength);
 
  const bool settingsShowRDOs = true; //FIXME: get from controller
  SoSeparator * errDetailed  = new SoSeparator;
  const std::vector<Identifier> rdolist = m_tgc->rdoList();
  if (rdolist.size() == 1 || !settingsShowRDOs)
  {
    errDetailed->addChild(common()->nodeManager()->getShapeNode_Strip(striplength,
      std::max(10.0,stripWidth), //strip width -> fixme: std::max hack for now since stripWidth returns 0.0
      3*0.8)); //strip thickness - hardcoded to be ~= the gas gap
    } else
    {
      VP1Msg::message("Warning: TGC has additional elements in rdoList: THIS HAS NEVER BEEN TESTED");
      SoSeparator * rdos = new SoSeparator;
 
      const Amg::Vector3D& globalposHIT = m_tgc->detectorElement()->channelPos( id );
    // get local position on wire plane, here we have to use a tolarance as the wire plane is located 2.5 CLHEP::mm
    // from the strip plane
      double tolerance = 3.;
      std::optional<Amg::Vector2D>localposHIT = m_tgc->detectorElement()->surface( id ).Trk::Surface::globalToLocal(globalposHIT,tolerance);
      if( !localposHIT )
      {
        localposHIT.emplace();
        localposHIT->setZero();
        VP1Msg::message("Warning: Local wire position is NULL");
      }
 
      rdos->addChild(common()->nodeManager()->getShapeNode_Strip(striplength,
        std::max(10.0,stripWidth), //strip width -> fixme: std::max hack for now since stripWidth returns 0.0
        3*0.8)); //strip thickness - hardcoded to be ~= the gas gap
 
      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 == id )
          continue;
        const Amg::Vector3D& globalposRDO = m_tgc->detectorElement()->channelPos( rdo_id );
        std::optional<Amg::Vector2D> localposRDO = m_tgc->detectorElement()->surface( rdo_id ).Trk::Surface::globalToLocal(globalposRDO,tolerance);
        if (!localposRDO)
        {
          VP1Msg::message("Warning: Local wire position is NULL");
          continue;
        }
 
        SoTranslation * localtrans = new SoTranslation;
        localtrans->translation.setValue((*localposRDO)[Trk::locX]-localposOLD[Trk::locX],(*localposRDO)[Trk::locY]-localposOLD[Trk::locY],0);
        rdos->addChild(localtrans);
 
        rdos->addChild(common()->nodeManager()->getShapeNode_Strip(striplength,
          std::max(10.0,stripWidth), //strip width -> fixme: std::max hack for now since stripWidth returns 0.0
          3*0.8)); //strip thickness - hardcoded to be ~= the gas gap
 
        localposOLD = *localposRDO;
      }
      errDetailed->addChild(rdos);
    }
    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");
}

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 * PRDHandle_TGC::createTransform ( ) const

virtual

Reimplemented from PRDHandleBase.

Definition at line 31 of file PRDHandle_TGC.cxx.

{
  const TgcIdHelper * idhelper = VP1DetInfo::tgcIDHelper();
  const Muon::TgcPrepData * prd = tgc();
  if (!prd || !idhelper) {
    return new SoTransform;
  }
 
  const Trk::Surface& theSurface = prd->detectorElement()->surface(prd->identify());
  SoTransform * theHitTransform =0;
   if (idhelper->isStrip( prd->identify() )){
    
    int stripNo = idhelper->channel(prd->identify()); 
    int gasGap  = idhelper->gasGap(prd->identify()); 
    
    const MuonGM::TgcReadoutElement* detEl = prd->detectorElement(); 
    
    // calculate two points along the tgc phi strip in the local tgc reference frame 
    const Amg::Vector3D lposTGC = detEl->stripCenterLocX(gasGap, stripNo, 0.) * Amg::Vector3D::UnitX();
    const double shift = detEl->getStationEta() > 0 ? 1.*Gaudi::Units::cm : - 1.*Gaudi::Units::cm;
    const  Amg::Vector3D lposTgcShifted = detEl->stripCenterLocX(gasGap, stripNo, shift) * Amg::Vector3D::UnitX()+
                                           shift * Amg::Vector3D::UnitY();
    const double angle = (lposTGC - lposTgcShifted).phi();
     // for phi strips, use sinstereo to get correct orientation
     Amg::RotationMatrix3D localRot{Amg::RotationMatrix3D::Identity()};
 
 
     // std::ostream os;
     // std::cout<<localRot.print(os)<<std::endl;
     //double angle  = M_PI/2.0 - asin(prd->detectorElement()->sinStereo(prd->identify())) ;
     // double angle  = asin(prd->detectorElement()->sinStereo(prd->identify())) ;
     // trying by trial and error to set right component of matrix!
     // setting theta=PI/4 rotated strips 45deg around z axis, but also a bit around x and y too! So clearly not working as expected.
      
    // localRot.setPhi (M_PI/4); 
    
     if (0 < detEl->getStationEta()) { 
       Amg::setPhi (localRot, angle);
     } else {
       Amg::setPhi (localRot, -angle);
     }
     // std::cout<<localRot.print(os)<<std::endl;
       // localRot.setPhi (angle); 
     
     // VP1Msg::message(QString::number(prd->identify().get_compact())+": angle="+QString::number(angle));
     localRot*=theSurface.transform().rotation ();
 
     Amg::Transform3D trans (localRot, theSurface.transform().translation());
     // Trk::GlobalDirection difPosGlobal = gpos_shift-gpos;
     // HepTransform3D toLocal = trans.inverse();
     // Trk::GlobalDirection difPosLocal = toLocal*difPosGlobal;
     // Trk::GlobalDirection difPosLocal2 = trans*difPosGlobal;
     
     // std::cout << " global strip direction " << difPosGlobal << " local " << difPosLocal << " local2 " << difPosLocal2 << std::endl;
     theHitTransform = VP1LinAlgUtils::toSoTransform(trans);
   } else {
     // eta strips just use surface orientiation
     theHitTransform = VP1LinAlgUtils::toSoTransform(theSurface.transform());
   }
 
   Amg::Vector3D theHitGPos= theSurface.localToGlobal(prd->localPosition());
   theHitTransform->translation.setValue((theHitGPos)[0], (theHitGPos)[1], (theHitGPos)[2]);
   return theHitTransform;
}

getPRD()

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

inlinevirtual

Implements PRDHandleBase.

Definition at line 24 of file PRDHandle_TGC.h.

{ return m_tgc; }

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

bool PRDHandle_TGC::inMuonChamber ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 26 of file PRDHandle_TGC.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.

parentMuonChamberPV()

GeoPVConstLink PRDHandle_TGC::parentMuonChamberPV ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 27 of file PRDHandle_TGC.h.

{ return m_tgc->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_TGC::regionIndex ( )

virtual

Reimplemented from PRDHandleBase.

Definition at line 185 of file PRDHandle_TGC.cxx.

  {
  //Unique for each station.
    const MuonGM::MuonStation* station = m_tgc->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();
  }
}

tgc()

const Muon::TgcPrepData* PRDHandle_TGC::tgc ( ) const

inline

Definition at line 23 of file PRDHandle_TGC.h.

{ return m_tgc; }

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_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_tgc

const Muon::TgcPrepData* PRDHandle_TGC::m_tgc

private

Definition at line 30 of file PRDHandle_TGC.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_TGC.h
• PRDHandle_TGC.cxx