#include <PRDHandle_CSC.h>

Inheritance diagram for PRDHandle_CSC:

Collaboration diagram for PRDHandle_CSC:

Public Member Functions
	PRDHandle_CSC (PRDCollHandle_CSC , const Muon::CscPrepData )
virtual	~PRDHandle_CSC ()
void	buildShapes (SoNode &shape_simple, SoNode &shape_detailed)
int	regionIndex ()
const Muon::CscPrepData *	csc () const
const Trk::PrepRawData *	getPRD () const
bool	inMuonChamber () const
GeoPVConstLink	parentMuonChamberPV () const
virtual QStringList	clicked () 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
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.

Private Attributes
const Muon::CscPrepData *	m_csc
Imp *	m_d
PRDSysCommonData *	m_common
bool	m_visible

Detailed Description

Definition at line 13 of file PRDHandle_CSC.h.

Constructor & Destructor Documentation

◆ PRDHandle_CSC()

PRDHandle_CSC::PRDHandle_CSC	(	PRDCollHandle_CSC *	collhandle,
		const Muon::CscPrepData *	csc )

Definition at line 24 of file PRDHandle_CSC.cxx.

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

◆ ~PRDHandle_CSC()

virtual PRDHandle_CSC::~PRDHandle_CSC ( )

inlinevirtual

Definition at line 17 of file PRDHandle_CSC.h.

17{};

Member Function Documentation

◆ buildShapes()

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

virtual

Implements PRDHandleBase.

Definition at line 31 of file PRDHandle_CSC.cxx.

{
  const CscIdHelper * idhelper = VP1DetInfo::cscIDHelper();
  if (!idhelper)
    return;
 
  Identifier id = m_csc->identify();
 
  SoSeparator * errDetailed  = new SoSeparator;
 
  double striplength = m_csc->detectorElement()->stripLength( id );
  if (static_cast<PRDCollHandle_CSC*>(collHandle())->project())
    striplength += 300.0;//Fixme: Rough extension for now
  shape_simple = common()->nodeManager()->getShapeNode_Strip(striplength);
 
  // Translation from objects local position to the strip position
  const Amg::Vector2D& localpos = m_csc->localPosition();
  const Amg::Vector3D& globalposHIT = m_csc->detectorElement()->stripPos( 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_csc->detectorElement()->surface( id ).Trk::Surface::globalToLocal(globalposHIT,tolerance); 
 
  if( !localposHIT )
  {
    localposHIT.emplace();
    localposHIT->setZero();
    VP1Msg::message("Warning: Local wire position is NULL");
  }
  SoTranslation * localtrans0 = new SoTranslation;
  localtrans0->translation.setValue((*localposHIT)[Trk::locX]-localpos[Trk::locX],(*localposHIT)[Trk::locY]-localpos[Trk::locY],0);
 
  const double maxCharge = 1e7;
  //Thickness of highlighed strip
  //FIXME: the maxCharge/3 is just a quick fix to get a nice presentation
  //       maybe this should go like sqrt(charge/maxCharge) so that it
  //       is clear for small charges.
  //const double stripHeight = (1 + 4*sqrt(m_csc->charge()/maxCharge))*CLHEP::mm; // [1,5]
  //
  //       The minimum is chosen equal to the neighbouring strips, such
  //       that a zero charge will give the same stirp height
  const double stripHeight = (1 + 4*m_csc->charge()/(maxCharge/3))*CLHEP::mm; // [1,5]
  //Thickness of additional strips (from rdoList())
  const double stripHeightRDOs = 1*CLHEP::mm;
 
  const std::vector<Identifier> rdolist = m_csc->rdoList();
  if (rdolist.size() == 1 || !collHandle()->drawRDOs())
  {
    errDetailed->addChild(localtrans0);
    errDetailed->addChild(common()->nodeManager()->getShapeNode_Strip(striplength,
                          m_csc->detectorElement()->cathodeReadoutPitch(idhelper->chamberLayer( id ),
                            idhelper->measuresPhi( id )), stripHeight));
 
    //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(striplength,
                          m_csc->detectorElement()->cathodeReadoutPitch(idhelper->chamberLayer( id ),
                            idhelper->measuresPhi( id )), stripHeight));
 
    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_csc->detectorElement()->stripPos( rdo_id );
 
      //      const Amg::Vector2D * localposRDO = m_csc->detectorElement()->surface( rdo_id ).globalToLocal(globalposRDO,tolerance);
      std::optional<Amg::Vector2D>  localposRDO = m_csc->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);
      double striplengthRDO = m_csc->detectorElement()->stripLength( rdo_id );
 
      rdos->addChild(common()->nodeManager()->getShapeNode_Strip(striplengthRDO,
                     m_csc->detectorElement()->cathodeReadoutPitch(idhelper->chamberLayer( rdo_id ),
                     idhelper->measuresPhi( rdo_id )), stripHeightRDOs));
 
      localposOLD = *localposRDO;
    }
    errDetailed->addChild(rdos);
  }
  
  SoMaterial * mat = new SoMaterial;
  mat->diffuseColor.setValue(1.0,0,0);
  errDetailed->addChild(mat);
  errDetailed->addChild(common()->nodeManager()->getShapeNode_Cross(10));
 
  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 PRDHandle_CSC::clicked ( ) const

virtual

Reimplemented from PRDHandleBase.

Definition at line 147 of file PRDHandle_CSC.cxx.

{
 
  if (!m_csc)
    return QStringList("Null PRD");
  const CscIdHelper * idhelper = VP1DetInfo::cscIDHelper();
  if (!idhelper)
    return QStringList("Null PRD");
 
  Identifier id = m_csc->identify();
  std::ostringstream os;
  os << "CscPrepData with Identifier ["<<id.get_compact()  ;
  os << "] = [" << (idhelper->print_to_string(id))  ;
  os << "]\n at global position = [" << m_csc->globalPosition()<<"], local position = ["<<m_csc->localPosition()<<"].";
  os << ". Charge="<<m_csc->charge()<<", time="<< m_csc->time()<<", status="<<Muon::toString(m_csc->status());
  os << ". RDO list= [";
  for( unsigned int i=0;i<m_csc->rdoList().size();++i) os << m_csc->rdoList()[i] << " , ";
  os << "].";
  // os << *(m_csc);
  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.

53{ return m_common; }

PRDHandleBase::m_common

PRDSysCommonData * m_common

Definition PRDHandleBase.h:84

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

◆ csc()

const Muon::CscPrepData * PRDHandle_CSC::csc ( ) const

inline

Definition at line 22 of file PRDHandle_CSC.h.

22{ return m_csc; }

◆ getPRD()

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

inlinevirtual

Implements PRDHandleBase.

Definition at line 23 of file PRDHandle_CSC.h.

23{ return m_csc; }

◆ getSecondPRD()

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

inlinevirtualinherited

Reimplemented in PRDHandle_SpacePoint.

Definition at line 73 of file PRDHandleBase.h.

73{ 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.

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

◆ inMuonChamber()

bool PRDHandle_CSC::inMuonChamber ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 26 of file PRDHandle_CSC.h.

26{ 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_CSC::parentMuonChamberPV ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 27 of file PRDHandle_CSC.h.

27{ return m_csc->detectorElement()->getMaterialGeom()->getParent(); }

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

virtual

Reimplemented from PRDHandleBase.

Definition at line 139 of file PRDHandle_CSC.cxx.

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

virtual bool PRDHandleBase::transformUsesSurfacePositionOnly ( ) const

inlinevirtualinherited

Reimplemented in PRDHandle_MDT, and PRDHandle_TRT.

Definition at line 47 of file PRDHandleBase.h.

47{ 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.

42{ return m_visible; }

Member Data Documentation

◆ m_common

PRDSysCommonData* PRDHandleBase::m_common

privateinherited

Definition at line 84 of file PRDHandleBase.h.

◆ m_csc

const Muon::CscPrepData* PRDHandle_CSC::m_csc

private

Definition at line 33 of file PRDHandle_CSC.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:

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PRDHandle_CSC()

◆ ~PRDHandle_CSC()

Member Function Documentation

◆ buildShapes()

◆ center()

◆ clicked()

◆ collHandle()

◆ common()

◆ createTransform()

◆ csc()

◆ getPRD()

◆ getSecondPRD()

◆ getTransform_CLHEP()

◆ highLight()

◆ inMuonChamber()

◆ isSane()

◆ parentMuonChamberPV()

◆ positionPRD()

◆ positionSecondPRD()

◆ regionIndex()

◆ sepDetailed()

◆ sepSimple()

◆ setVisible()

◆ transformUsesSurfacePositionOnly()

◆ update3DObjects()

◆ updateMaterial()

◆ visible()

Member Data Documentation

◆ m_common

◆ m_csc

◆ m_d

◆ m_visible