PRDHandle_CSC_Strip Class Reference

#include <PRDHandle_CSC_Strip.h>

Inheritance diagram for PRDHandle_CSC_Strip:

Collaboration diagram for PRDHandle_CSC_Strip:

Public Member Functions
	PRDHandle_CSC_Strip (PRDCollHandle_CSC_Strip *, const Muon::CscStripPrepData *)
virtual	~PRDHandle_CSC_Strip ()
void	buildShapes (SoNode *&shape_simple, SoNode *&shape_detailed)
int	regionIndex ()
const Muon::CscStripPrepData *	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::CscStripPrepData *	m_csc
Imp *	m_d
PRDSysCommonData *	m_common
bool	m_visible

Detailed Description

Definition at line 13 of file PRDHandle_CSC_Strip.h.

Constructor & Destructor Documentation

PRDHandle_CSC_Strip()

PRDHandle_CSC_Strip::PRDHandle_CSC_Strip	(	PRDCollHandle_CSC_Strip *	collhandle,
		const Muon::CscStripPrepData *	csc )

Definition at line 25 of file PRDHandle_CSC_Strip.cxx.

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

~PRDHandle_CSC_Strip()

virtual PRDHandle_CSC_Strip::~PRDHandle_CSC_Strip ( )

inlinevirtual

Definition at line 17 of file PRDHandle_CSC_Strip.h.

{};

Member Function Documentation

buildShapes()

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

virtual

Implements PRDHandleBase.

Definition at line 32 of file PRDHandle_CSC_Strip.cxx.

{
  const CscIdHelper * idhelper = VP1DetInfo::cscIDHelper();
  if (!idhelper)
    return;
 
  Identifier id = m_csc->identify();
 
  SoSeparator * errDetailed  = new SoSeparator;
  SoSeparator * errSimple  = new SoSeparator;
  shape_simple  = errSimple;
 
  double striplength = m_csc->detectorElement()->stripLength( id );
  if (static_cast<PRDCollHandle_CSC_Strip*>(collHandle())->project())
    striplength += 300.0;//Fixme: Rough extension for now
  SoNode* line = common()->nodeManager()->getShapeNode_Strip(striplength);
  errSimple->addChild(line);
  errSimple->addChild(common()->nodeManager()->getShapeNode_Cross(10));
  const double maxCharge = 1e5;
  //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*CLHEP::mm; // [1,5]
 
  const std::vector< float > charges = m_csc->sampleCharges();
  if (charges.size() == 0 )
  {
        // Add arbitrary strip height & draw entire length of strip.
    errDetailed->addChild(common()->nodeManager()->getShapeNode_Strip(striplength,
      m_csc->detectorElement()->cathodeReadoutPitch(idhelper->chamberLayer( id ),
      idhelper->measuresPhi( id )), stripHeight));
  } else {
    SoSeparator * rdos = new SoSeparator;
 
    Amg::Vector2D localposCharge;
 
        // SoTransparency * transparent = new SoTransparency;
        // transparent->transparency.setValue(0.5);
        // rdos->addChild( transparent );
    double chargeLength=striplength/static_cast<double>(charges.size()); // size of one charge cell
        // std::cout<<"Number of charges="<<charges.size()<<", chargeLength="<<chargeLength<<" striplength="<<striplength<<std::endl;
    unsigned int i=0;
    double ypos=(chargeLength-striplength)/2.0;
    std::vector<float>::const_iterator st = charges.begin();
    std::vector<float>::const_iterator en = charges.end();
    double stripPitch=m_csc->detectorElement()->cathodeReadoutPitch(idhelper->chamberLayer( id ),idhelper->measuresPhi( id ));
    SoTranslation * localtrans = new SoTranslation;
    localtrans->translation.setValue(0.0,ypos,0.0);
    rdos->addChild(localtrans);
 
    for( std::vector<float>::const_iterator it = st; it!=en; ++it, ++i){
            // std::cout<<"Charge: "<<i<<" \ty="<<ypos<<" has charge=:"<<*it<<std::endl;  
      float charge = std::max(1.0f,(*it)); // keep charges positive until I understand if it is okay to be -ve
      const double stripHeightCharge = (1.0 + 4*sqrt(charge/maxCharge))*CLHEP::mm;
            // std::cout<<"stripHeightCharge: "<<stripHeightCharge<<std::endl;
      rdos->addChild(common()->nodeManager()->getShapeNode_Strip(chargeLength, stripPitch, stripHeightCharge));
            // ypos+=chargeLength;
      SoTranslation * localtrans2 = new SoTranslation;
      localtrans2->translation.setValue(0.0,chargeLength,0.0);
      rdos->addChild(localtrans2);
 
    }
    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_Strip::clicked ( ) const

virtual

Reimplemented from PRDHandleBase.

Definition at line 118 of file PRDHandle_CSC_Strip.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 << "CscStripPrepData with Identifier ["<<id<<"] = ["<<idhelper->print_to_string(id)
    <<"] at global position = "<<m_csc->globalPosition()<<", local position = "<<m_csc->localPosition();
  os << ". Samples = [";
  for( unsigned int i=0;i<m_csc->sampleCharges().size();++i) os << m_csc->sampleCharges()[i] << " , ";
 
  os << "]. Time of first sample= "<<m_csc->timeOfFirstSample() << ", sampling time (CLHEP::ns) " 
    << m_csc->samplingTime() << ", sampling phase " << m_csc->samplingPhase();
  os << ". RDO list= [";
  for( unsigned int i=0;i<m_csc->rdoList().size();++i) os << m_csc->rdoList()[i] << " , ";
  os << "].";
      // os <<"\n Dump follows:";
      // os <<"\n--------------";
    //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.

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

csc()

const Muon::CscStripPrepData * PRDHandle_CSC_Strip::csc ( ) const

inline

Definition at line 22 of file PRDHandle_CSC_Strip.h.

{ return m_csc; }

getPRD()

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

inlinevirtual

Implements PRDHandleBase.

Definition at line 23 of file PRDHandle_CSC_Strip.h.

{ return m_csc; }

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_CSC_Strip::inMuonChamber ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 26 of file PRDHandle_CSC_Strip.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_CSC_Strip::parentMuonChamberPV ( ) const

inlinevirtual

Reimplemented from PRDHandleBase.

Definition at line 27 of file PRDHandle_CSC_Strip.h.

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

virtual

Reimplemented from PRDHandleBase.

Definition at line 110 of file PRDHandle_CSC_Strip.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.

{ 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_csc

const Muon::CscStripPrepData* PRDHandle_CSC_Strip::m_csc

private

Definition at line 33 of file PRDHandle_CSC_Strip.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_CSC_Strip.h
• PRDHandle_CSC_Strip.cxx