IParticleHandle_TrackParticle Class Reference

#include <IParticleHandle_TrackParticle.h>

Inheritance diagram for IParticleHandle_TrackParticle:

Collaboration diagram for IParticleHandle_TrackParticle:

Classes
class	Imp

Public Member Functions
	IParticleHandle_TrackParticle (IParticleCollHandleBase *, const xAOD::TrackParticle *)
virtual	~IParticleHandle_TrackParticle ()
virtual bool	has3DObjects ()
	Returns true if the 3D objects have been created.
virtual void	clear3DObjects ()
	Delete objects.
virtual SoNode *	nodes ()
	Returns the 3Dobjects.
void	addParameterShapes ()
void	addLine_FromTrackParticle ()
	Make the line, but only using information on the track (a spline fit will be done from the parameters)
void	addLine_Extrapolated ()
	Extrapolate the initial track parameters to get the line (a spline fit is done between the extrapolated parameters)
void	fillLineFromSplineFit (const std::vector< Amg::Vector3D > &positions, const std::vector< Amg::Vector3D > &momenta)
	Fill d->line with a SoLineSet made from a spline fit between the information in positions and momenta.
void	linear (Amg::Vector3D &output, Amg::Vector3D &a, Amg::Vector3D &b, float t)
void	bezier (Amg::Vector3D &output, Amg::Vector3D &a, Amg::Vector3D &b, Amg::Vector3D &c, Amg::Vector3D &d, float t)
virtual QStringList	clicked () const
	Called when user selects the node (stringlist is displayed in messagebox).
virtual Amg::Vector3D	momentum () const
virtual Amg::Vector3D	position () const
const xAOD::IParticle &	iParticle () const
virtual double	charge () const
	Returns unknown() in case of trouble.
virtual unsigned	summaryValue (xAOD::SummaryType type) const
virtual unsigned	getNPixelHits () const
virtual unsigned	getNSCTHits () const
virtual unsigned	getNTRTHits () const
virtual unsigned	getNMuonPrecisionLayers () const
virtual unsigned	getNMuonPrecisionHoleLayers () const
virtual unsigned	getNMuonPhiLayers () const
virtual unsigned	getNMuonPhiHoleLayers () const
virtual QString	type () const
	return very short word with type (maybe link with collection type?)
QString	shortInfo () const
	This returns the information shown about the object in the object browser.
void	fillObjectBrowser (QList< QTreeWidgetItem * > &listOfItems)
	Create and fill the object browser QTreeWidgetItem.
const QList< std::pair< xAOD::ParameterPosition, Amg::Vector3D > > &	parametersAndPositions () const
virtual void	dumpToJSON (std::ofstream &) const
SoMaterial *	determineMaterial ()
	Should be implemented by children, in order to change the material depending on the interface etc.
virtual QStringList	baseInfo () const
virtual double	mass () const
	Returns unknown() in case of trouble.
bool	hasCharge () const
bool	hasMass () const
const AODSysCommonData *	common () const
void	setVisible (bool)
	use by the collection handle.
bool	visible () const
const AODCollHandleBase *	collHandle () const
virtual void	update3DObjects ()
	Called after some configuration changes, or when the object is first created. Must be overloaded by children.
virtual void	rebuild3DObjects ()
virtual void	attach3DObjects ()
	Connect to the scene graph.
virtual void	detach3DObjects ()
	Remove from the scene graph.
void	updateMaterial ()
	Called after some configuration related to material changes.
SoMaterial *	currentMaterial () const
	Returns the current material of the handle.
void	updateObjectBrowser ()
	Update object browser QTreeWidgetItem.
QTreeWidgetItem *	browserTreeItem () const
	Return the QTreeWidgetItem;.

Static Public Member Functions
static int	numberOfInstances ()
static double	unknown ()

Protected Member Functions
virtual void	visibleStateChanged ()
	Dummy implementation does nothing.
virtual void	currentMaterialChanged ()

Private Member Functions
void	ensureInitSubSysHitInfo () const

Private Attributes
Imp *	m_d
bool	m_visible
AODCollHandleBase *	m_collhandle
SoMaterial *	m_currentmaterial

Detailed Description

Definition at line 27 of file IParticleHandle_TrackParticle.h.

Constructor & Destructor Documentation

IParticleHandle_TrackParticle()

IParticleHandle_TrackParticle::IParticleHandle_TrackParticle	(	IParticleCollHandleBase *	ch,
		const xAOD::TrackParticle *	tp )

Definition at line 62 of file IParticleHandle_TrackParticle.cxx.

  : IParticleHandleBase(ch), m_d(new Imp)
{ 
  m_d->theclass = this;
  m_d->trackparticle = tp;
  m_d->collHandle = dynamic_cast<IParticleCollHandle_TrackParticle*>(ch);
  m_d->sep=0;
  m_d->line=0;
}

~IParticleHandle_TrackParticle()

IParticleHandle_TrackParticle::~IParticleHandle_TrackParticle ( )

virtual

Definition at line 73 of file IParticleHandle_TrackParticle.cxx.

{
  if (m_d->line) m_d->line->unref();
  if (m_d->sep) m_d->sep->unref();
  delete m_d;
}

Member Function Documentation

addLine_Extrapolated()

void IParticleHandle_TrackParticle::addLine_Extrapolated

(

)

Extrapolate the initial track parameters to get the line (a spline fit is done between the extrapolated parameters)

Definition at line 295 of file IParticleHandle_TrackParticle.cxx.

                                                        {
  VP1Msg::messageVerbose("IParticleHandle_TrackParticle::addLine_Extrapolated().");
 
  const Trk::Perigee& peri = m_d->trackparticle->perigeeParameters (); // FIXME - I'd quite like not to use anything which requires Athena ...
  
  Trk::CurvilinearParameters startParameters(peri.position(),peri.momentum(),peri.charge());
  Trk::ExtrapolationCell<Trk::TrackParameters> ecc(startParameters);
  // ecc.addConfigurationMode(Trk::ExtrapolationMode::StopAtBoundary);
  ecc.addConfigurationMode(Trk::ExtrapolationMode::CollectPassive);
  ecc.addConfigurationMode(Trk::ExtrapolationMode::CollectSensitive);
  
  // we want to extrapolate outwards only for muon tracks
  if (getNMuonPrecisionLayers()==0 &&  getNMuonPhiLayers()==0 )
    ecc.addConfigurationMode(Trk::ExtrapolationMode::StopAtBoundary);
  
  // call the extrapolation engine
  const Trk::IExtrapolationEngine * engine = collHandle()->common()->extrapolator();
  if (!engine) {
     VP1Msg::messageVerbose("No extrapolator engine available. Aborting.");
     return;
  }
  
  std::vector<Amg::Vector3D> positions; 
  std::vector<Amg::Vector3D> momenta;
 
  // Add start positions
  positions.push_back(Amg::Vector3D(peri.position().x(),peri.position().y(),peri.position().z()));
  momenta.  push_back(Amg::Vector3D(peri.momentum().x(),peri.momentum().y(),peri.momentum().z()));
 
  
  Trk::ExtrapolationCode eCode = engine->extrapolate(ecc);
  
  // end the parameters if there
  if (eCode.isSuccess()){
    // loop over the collected information
    for (auto& es : ecc.extrapolationSteps){
      
      // continue if we have parameters
      const Trk::TrackParameters* parameters = es.parameters;
      if (parameters){
        // take what you  need -> position etc
 
        // finally delete the parameters
        Amg::Vector3D pos = parameters->position();
        Amg::Vector3D mom = parameters->momentum();
        positions.push_back(pos);
        momenta.push_back(mom);
        
        delete parameters;
      }
    }//
    
    Amg::Vector3D pos = ecc.endParameters->position();
    Amg::Vector3D mom = ecc.endParameters->momentum();
    positions.push_back(pos);
    momenta.push_back(mom);
    // and delete the final parameter
    delete ecc.endParameters;
  } else {
    VP1Msg::messageVerbose("IParticleHandle_TrackParticle::addLine_Extrapolated() - extrapolation failed.");
  }
  fillLineFromSplineFit(positions, momenta);
  // std::cout<<"About to add "<<collHandleTrk->material()<<std::endl;
  m_d->sep->addChild(m_d->collHandle->material()); 
  // ^^ FIXME - should rearrange so we don't need to reset material
 
 
  m_d->sep->addChild(m_d->line);
}

addLine_FromTrackParticle()

void IParticleHandle_TrackParticle::addLine_FromTrackParticle

(

)

Make the line, but only using information on the track (a spline fit will be done from the parameters)

Definition at line 173 of file IParticleHandle_TrackParticle.cxx.

                                                             {
 
  std::vector<Amg::Vector3D> positions, momenta;
 
#if defined BUILDVP1LIGHT
  positions.push_back(position());
  momenta.  push_back(momentum());
#else
  const Trk::Perigee& peri = m_d->trackparticle->perigeeParameters (); // FIXME - I'd quite like not to use anything which requires Athena ...
 
  positions.push_back(Amg::Vector3D(peri.position().x(),peri.position().y(),peri.position().z()));
  momenta.  push_back(Amg::Vector3D(peri.momentum().x(),peri.momentum().y(),peri.momentum().z()));
  // std::cout<<"i:"<<0<<"/"<<m_d->trackparticle->numberOfParameters()+1<<": ("<<peri.position().x()<<","<<peri.position().y()<<","<<peri.position().z()<<")"<<std::endl;
#endif // BUILDVP1LIGHT
 
 
  VP1Msg::messageVerbose("IParticleHandle_TrackParticle::addLine_FromTrackParticle - has "+QString::number(m_d->trackparticle->numberOfParameters())+" extra parameters.");
   for (unsigned int i=0; i<m_d->trackparticle->numberOfParameters() ; ++i){
    // std::cout<<"i:"<<i+1<<"/"<<m_d->trackparticle->numberOfParameters()+1<<": ("<<m_d->trackparticle->parameterX(i)<<","<<m_d->trackparticle->parameterY(i)<<","<<m_d->trackparticle->parameterZ(i)<<")"<<std::endl;
 
    positions.push_back(Amg::Vector3D(m_d->trackparticle->parameterX(i),
                                      m_d->trackparticle->parameterY(i),
                                      m_d->trackparticle->parameterZ(i)));      
    momenta.push_back(Amg::Vector3D(  m_d->trackparticle->parameterPX(i),
                                      m_d->trackparticle->parameterPY(i),
                                      m_d->trackparticle->parameterPZ(i)));      
 
  } // end of loop.
 
  // if ( positions.size()<2 ) VP1Msg::messageVerbose("IParticleHandle_TrackParticle::addLine_FromTrackParticle - WARNING - not enough points to make a line.");
 
  fillLineFromSplineFit(positions, momenta);
 
  // std::cout<<"About to add "<<collHandleTrk->material()<<std::endl;
  m_d->sep->addChild(m_d->collHandle->material()); 
  // ^^ FIXME - should rearrange so we don't need to reset material
 
  m_d->sep->addChild(m_d->line);
}

addParameterShapes()

void IParticleHandle_TrackParticle::addParameterShapes

(

)

Definition at line 128 of file IParticleHandle_TrackParticle.cxx.

                                                      {
  // Points:
  SoPointSet* parameterPoints = 0;
  SoVertexProperty * vertices2 = 0;
 
  // std::cout<<"showing parameters... have "<<m_d->trackparticle->numberOfParameters()+1<<" to draw"<<std::endl;
  
      // Okay. so let's make the containers to hold them...
  parameterPoints = new SoPointSet;
  vertices2 = new SoVertexProperty;
 
  int i = 0;
  bool colourByParamType = m_d->collHandle->collSettingsButton().colourParametersByType();
  QList<std::pair<xAOD::ParameterPosition, Amg::Vector3D> > parAndPoses = parametersAndPositions();
  
  for (std::pair<xAOD::ParameterPosition, Amg::Vector3D> parpos : parAndPoses ) {
    Amg::Vector3D pos = parpos.second;
 
    if (colourByParamType) {
      // Okay, so we treat each point separately & add materials as we go...
      m_d->sep->addChild(m_d->collHandle->collSettingsButton().parameterMaterial(parpos.first));
      vertices2->vertex.set1Value(0,pos.x(),pos.y(),pos.z());
      parameterPoints->numPoints=1;
      parameterPoints->vertexProperty = vertices2;
      m_d->sep->addChild(parameterPoints);
      parameterPoints = new SoPointSet; // reset for next point
      vertices2 = new SoVertexProperty;
    } else { // Not colouring by type
      vertices2->vertex.set1Value(i++,pos.x(),pos.y(),pos.z());
    } 
  }
 
  if (colourByParamType) {
    // In this case, we don't need the last Nodes created in the loop.
    parameterPoints->unref();
    vertices2->unref();
  
  } else {  
    parameterPoints->numPoints=i;
    parameterPoints->vertexProperty = vertices2;
    m_d->sep->addChild(static_cast<const IParticleCollHandle_TrackParticle*>(collHandle())->collSettingsButton().defaultParameterMaterial());
    m_d->sep->addChild(parameterPoints);
  }
}

attach3DObjects()

void AODHandleBase::attach3DObjects ( )

virtualinherited

Connect to the scene graph.

Definition at line 157 of file AODHandleBase.cxx.

{
  VP1Msg::messageVerbose("AODHandleBase::attach3DObjects");
 
  if (!m_currentmaterial) {
    m_currentmaterial = determineMaterial();
    m_currentmaterial->ref();
    currentMaterialChanged();
  }
  if (has3DObjects() && m_collhandle->sepHelper()) {
    VP1Msg::messageVerbose("Adding nodes to sep helper...");
  
    m_collhandle->common()->registerHandle(this);
 
    // debug
    // std::cout<<"m_collhandle->sepHelper() = " << m_collhandle->sepHelper()<<std::endl;
    SoNode* tmp = nodes();
    // std::cout << "tmp nodes = "<< tmp << std::endl;
 
    //      m_collhandle->sepHelper()->addNodeUnderMaterial( nodes(), m_currentmaterial );
    m_collhandle->sepHelper()->addNodeUnderMaterial( tmp, m_currentmaterial );
 
  } else {
    if (!has3DObjects()) {
      VP1Msg::message("AODHandleBase::attach3DObjects() - No 3D objects to attach!");
    }
    if (!m_collhandle->sepHelper()) VP1Msg::message("AODHandleBase::attach3DObjects() - No sep helper!");
  }
  // if (label_sep && theclass->common()->textSep() && theclass->common()->controller()->doTrackLabels())
  //   theclass->common()->textSep()->addChild(label_sep);
  //
}

baseInfo()

QStringList IParticleHandleBase::baseInfo ( ) const

virtualinherited

Reimplemented from AODHandleBase.

Definition at line 85 of file IParticleHandleBase.cxx.

{
  QStringList l;
  Amg::Vector3D mom = momentum()/1000;
  if (mom.mag2()==0.0) {
    l << "Momentum : 0 (undefined)";
  } else {
    l << "Momentum [GeV]: "+VP1Msg::str(mom);
    l << "|Pt|/|P| [GeV]: "+VP1Msg::str(mom.perp())+" / " + VP1Msg::str(mom.mag());
    l << VP1Msg::str("(")+QChar(0x03B7)+","+QChar(0x03D5)+VP1Msg::str(")=(")
      +VP1Msg::str(mom.eta())+VP1Msg::str(",")+VP1Msg::str(VP1LinAlgUtils::phiFromXY(mom.x(), mom.y() ))+VP1Msg::str(")");
    
    l << "Eta: "+VP1Msg::str(mom.eta());
    l << "Phi: "+VP1Msg::str(VP1LinAlgUtils::phiFromXY(mom.x(), mom.y() ));
  }
  // l<<"Hit summary: Pix["+VP1Msg::str(getNPixelHits())+"], SCT["+VP1Msg::str(getNSCTHits())+"], TRT["+VP1Msg::str(getNTRTHits())
  //   +"], MDT["+QString::number(getNMDTHits())+"], RPC["+QString::number(getNRPCHits())+"], TGC["+QString::number(getNTGCHits())+"], CSC["+QString::number(getNCSCHits())+"].";
  // 
  // return l;
  // int pdg = pdgCode();
  // if (pdg) {
  //   bool ok;
  //   QString name = VP1ParticleData::particleName(pdg,ok);
  //   if (ok)
  //     l << "pdg: "+QString::number(pdg)+" ["+name+"]";
  //   else
  //     l << "pdg: "+QString::number(pdg);
  // }
 
  return l;
}

bezier()

void IParticleHandle_TrackParticle::bezier	(	Amg::Vector3D &	output,
		Amg::Vector3D &	a,
		Amg::Vector3D &	b,
		Amg::Vector3D &	c,
		Amg::Vector3D &	d,
		float	t )

Definition at line 220 of file IParticleHandle_TrackParticle.cxx.

                                                                                                                                      {
  // DeCasteljau subdivision
  Amg::Vector3D ab,bc,cd,abbc,bccd;
  linear(ab, a,b,t);           // point between a and b (green)
  linear(bc, b,c,t);           // point between b and c (green)
  linear(cd, c,d,t);           // point between c and d (green)
  linear(abbc, ab,bc,t);       // point between ab and bc (blue)
  linear(bccd, bc,cd,t);       // point between bc and cd (blue)
  linear(output, abbc,bccd,t);   // point on the bezier-curve (black)
}

browserTreeItem()

QTreeWidgetItem * AODHandleBase::browserTreeItem ( ) const

inherited

Return the QTreeWidgetItem;.

Definition at line 311 of file AODHandleBase.cxx.

                                                      {
  return m_d->m_objBrowseTree;
}

charge()

double IParticleHandle_TrackParticle::charge ( ) const

virtual

Returns unknown() in case of trouble.

Implements IParticleHandleBase.

Definition at line 448 of file IParticleHandle_TrackParticle.cxx.

{
  return m_d->trackparticle->charge();
}

clear3DObjects()

void IParticleHandle_TrackParticle::clear3DObjects ( )

virtual

Delete objects.

Implements AODHandleBase.

Definition at line 83 of file IParticleHandle_TrackParticle.cxx.

                                                  {
  VP1Msg::messageVerbose("IParticleHandle_TrackParticle::clear3DObjects()");
  
  if (m_d->line) {
    m_d->line->unref();
    m_d->line=0;
  }
  if (m_d->sep) {
    m_d->sep->unref();
    m_d->sep=0;
  }
  
  // Wipe caches
  m_d->positionsToWrite.clear();
  m_d->parametersAndPositions.clear();
}

clicked()

QStringList IParticleHandle_TrackParticle::clicked ( ) const

virtual

Called when user selects the node (stringlist is displayed in messagebox).

Reimplementations should add info momInfo() to get momentum, eta, phi, pid,charge, mass printed.

Implements AODHandleBase.

Definition at line 367 of file IParticleHandle_TrackParticle.cxx.

{
  QStringList l;
  l << "TrackParticle:";
  l << IParticleHandleBase::baseInfo();
  return l;
}

collHandle()

const AODCollHandleBase * AODHandleBase::collHandle ( ) const

inlineinherited

Definition at line 55 of file AODHandleBase.h.

{ return m_collhandle; }

common()

const AODSysCommonData * AODHandleBase::common ( ) const

inherited

Definition at line 93 of file AODHandleBase.cxx.

{
  return std::as_const(m_collhandle)->common();
}

currentMaterial()

SoMaterial * AODHandleBase::currentMaterial ( ) const

inlineinherited

Returns the current material of the handle.

Definition at line 76 of file AODHandleBase.h.

{ return m_currentmaterial; }

currentMaterialChanged()

virtual void AODHandleBase::currentMaterialChanged ( )

inlineprotectedvirtualinherited

Definition at line 111 of file AODHandleBase.h.

{};

detach3DObjects()

void AODHandleBase::detach3DObjects ( )

virtualinherited

Remove from the scene graph.

Definition at line 191 of file AODHandleBase.cxx.

{
  VP1Msg::messageVerbose("AODHandleBase::detach3DObjects");
    
  if (!m_currentmaterial)
    return;//Can never have been attached!
  if (has3DObjects() && m_collhandle->sepHelper()) {
    m_collhandle->sepHelper()->removeNodeUnderMaterial(nodes(),m_currentmaterial);
  }
  
  m_collhandle->common()->deregisterHandle(this);
  
  // if (label_sep && theclass->common()->textSep()) 
  //   theclass->common()->textSep()->removeChild(label_sep);
}

determineMaterial()

SoMaterial * IParticleHandleBase::determineMaterial ( )

virtualinherited

Should be implemented by children, in order to change the material depending on the interface etc.

This will (by default) return the collection material.

Implements AODHandleBase.

Definition at line 71 of file IParticleHandleBase.cxx.

                                                    {
  // By default we use the collection material.
  // std::cout<<"AODHandleBase::determineMaterial() - collHandle()->material()"<<collHandle()->material()<<std::endl;
  return collHandle()->material();
}

dumpToJSON()

void IParticleHandle_TrackParticle::dumpToJSON ( std::ofstream & str ) const

virtual

Reimplemented from AODHandleBase.

Definition at line 597 of file IParticleHandle_TrackParticle.cxx.

                                                                      {
  str << "\"chi2\":"<<m_d->trackparticle->chiSquared () <<", ";
  str << "\"dof\":" <<m_d->trackparticle->numberDoF () <<", ";
  str << "\"dparams\": [" 
      << m_d->trackparticle->d0() <<", "
      << m_d->trackparticle->z0() <<", "
      << m_d->trackparticle->phi0() <<", "
      << m_d->trackparticle->theta() <<", "
      << m_d->trackparticle->qOverP()
      << "], \"pos\": [";
  bool first = true;
  for (auto pos: m_d->positionsToWrite){
    if (first){
      first = false;
    } else {
      str <<",";
    }
    str << "["<<pos.x()<<","<<pos.y()<<","<<pos.z()<<"]";
  }
  str << "] ";
}

ensureInitSubSysHitInfo()

void IParticleHandle_TrackParticle::ensureInitSubSysHitInfo ( ) const

private

fillLineFromSplineFit()

void IParticleHandle_TrackParticle::fillLineFromSplineFit	(	const std::vector< Amg::Vector3D > &	positions,
		const std::vector< Amg::Vector3D > &	momenta )

Fill d->line with a SoLineSet made from a spline fit between the information in positions and momenta.

Definition at line 231 of file IParticleHandle_TrackParticle.cxx.

                                                                                                                                      {
  int iver(0);
  unsigned int npointsused(0);
  // Create new lineset etc
  m_d->line = new SoLineSet();
  m_d->line->ref();
  
  // For the moment, lets just draw the lines.
  SoVertexProperty * vertices = new SoVertexProperty();
  
  //  std::cout<<"About to loop over "<<positions.size()<<" points to make line"<<std::endl;
  // for (unsigned int i=0; i<positions.size(); ++i){
  //   std::cout<<"i: "<<i<<" at ("<<positions[i].x()<<","<<positions[i].y()<<","<<positions[i].z()<<")"<<" \tmom = ("<<momenta[i].x()<<","<<momenta[i].y()<<","<<momenta[i].z()<<")"<<std::endl;
  // }
 
  // std::cout<<"Now do interpolations."<<std::endl;
 
  bool useExistingParams = m_d->collHandle->collSettingsButton().useExistingParameters();
 
  for (unsigned int i=1; i<positions.size(); ++i){
    // std::cout<<"i: "<<i<<" at ("<<positions[i].x()<<","<<positions[i].y()<<","<<positions[i].z()<<")"<<" \tmom = ("<<momenta[i].x()<<","<<momenta[i].y()<<","<<momenta[i].z()<<")"<<std::endl;
    
    if (useExistingParams){ 
      float scale = ( (positions[i]-positions[i-1]).mag() )/3.0; // want to split curve into three parts
      Amg::Vector3D p0(positions[i-1]); // first position
      Amg::Vector3D p1(positions[i-1]+(scale * momenta[i-1].unit() ) ); // move 1/3 or the way towards next point, but following momentum direction
      Amg::Vector3D p2(positions[i]  -(scale * momenta[i].unit() ) ); // back 1/3 from second point, but following momentum direction
      Amg::Vector3D p3(positions[i]);
      // std::cout<<"p0: ("<<p0[0]<<", "<<p0[1]<<", "<<p0[2]<<") "<<std::endl;
      // std::cout<<"p1: ("<<p1[0]<<", "<<p1[1]<<", "<<p1[2]<<") "<<std::endl;
      // std::cout<<"p2: ("<<p2[0]<<", "<<p2[1]<<", "<<p2[2]<<") "<<std::endl;
      // std::cout<<"p3: ("<<p3[0]<<", "<<p3[1]<<", "<<p3[2]<<") "<<std::endl;
    
      unsigned int maxCount=m_d->collHandle->collSettingsButton().numOfStepsForInterpolation();
      for (unsigned int count=1; count<=maxCount;++count){
        float t=static_cast<float>(count)/static_cast<float>(maxCount);
        Amg::Vector3D pos;
        bezier(pos, p0, p1, p2, p3, t);
        // std::cout<<"j: "<<count<<" \tat ("<<pos.x()<<","<<pos.y()<<","<<pos.z()<<")"<<std::endl;
      
        // pos = pos + momdelta;
        vertices->vertex.set1Value(iver++,pos.x(),pos.y(),pos.z());
        m_d->positionsToWrite.append(pos);
      
        npointsused++;        
      }
    } else {
      vertices->vertex.set1Value(iver++,positions[i].x(),positions[i].y(),positions[i].z());
      m_d->positionsToWrite.append(positions[i]);
      npointsused++;        
    }
  }
  
  // Add final point.
  
  
    
  // Add to SoLine set
  m_d->line->numVertices.set1Value(0,npointsused);
  m_d->line->vertexProperty = vertices;
  
}

fillObjectBrowser()

void IParticleHandle_TrackParticle::fillObjectBrowser ( QList< QTreeWidgetItem * > & list )

virtual

Create and fill the object browser QTreeWidgetItem.

Reimplemented from AODHandleBase.

Definition at line 514 of file IParticleHandle_TrackParticle.cxx.

{
  IParticleHandleBase::fillObjectBrowser(listOfItems); // Obligatory!
 
  QTreeWidgetItem* TSOSitem = new QTreeWidgetItem(browserTreeItem());
  TSOSitem->setText(0, QString("Def. Parameters " ) );
  QString dParameters("(");
  dParameters+=QString::number(m_d->trackparticle->d0());
  dParameters+=", ";     
  dParameters+=QString::number(m_d->trackparticle->z0());
  dParameters+=", ";     
  dParameters+=QString::number(m_d->trackparticle->phi0());
  dParameters+=", ";     
  dParameters+=QString::number(m_d->trackparticle->theta());
  dParameters+=", ";     
  dParameters+=QString::number(m_d->trackparticle->qOverP());
  dParameters+=")";     
  TSOSitem->setText(1, dParameters );
  
  for (unsigned int i=0; i<m_d->trackparticle->numberOfParameters() ; ++i){
    QTreeWidgetItem* TSOSitem = new QTreeWidgetItem(browserTreeItem());
    TSOSitem->setText(0, QString("Parameter "+QString::number( i+1 ) ) );
    QString pos(", Position = (");  
    pos+=QString::number(m_d->trackparticle->parameterX(i));
    pos+=", ";    
    pos+=QString::number(m_d->trackparticle->parameterY(i));
    pos+=", ";
    pos+=QString::number(m_d->trackparticle->parameterZ(i));
    pos+=")";
    
    switch (m_d->trackparticle->parameterPosition(i)){
      case xAOD::BeamLine: 
        TSOSitem->setText(1, QString("BeamLine" )+pos );
        break;
      case xAOD::FirstMeasurement:
        TSOSitem->setText(1, QString("FirstMeasurement")+pos );
        break;
      case xAOD::LastMeasurement: 
        TSOSitem->setText(1, QString("LastMeasurement" )+pos );
        break;
      case xAOD::CalorimeterEntrance:
        TSOSitem->setText(1, QString("CalorimeterEntrance")+pos );
        break;
      case xAOD::CalorimeterExit: 
        TSOSitem->setText(1, QString("CalorimeterExit" )+pos );
        break;
      case xAOD::MuonSpectrometerEntrance:
        TSOSitem->setText(1, QString("MuonSpectrometerEntrance")+pos );
        break;
      default:
        TSOSitem->setText(1, QString("Undefined")+pos );
    }
  }
  
  // TODO - add more.
}

getNMuonPhiHoleLayers()

unsigned IParticleHandle_TrackParticle::getNMuonPhiHoleLayers ( ) const

virtual

Definition at line 497 of file IParticleHandle_TrackParticle.cxx.

{ 
  return summaryValue(xAOD::numberOfPhiHoleLayers);
}

getNMuonPhiLayers()

unsigned IParticleHandle_TrackParticle::getNMuonPhiLayers ( ) const

virtual

Definition at line 492 of file IParticleHandle_TrackParticle.cxx.

{ 
  return summaryValue(xAOD::numberOfPhiLayers);
}

getNMuonPrecisionHoleLayers()

unsigned IParticleHandle_TrackParticle::getNMuonPrecisionHoleLayers ( ) const

virtual

Definition at line 487 of file IParticleHandle_TrackParticle.cxx.

{ 
  return summaryValue(xAOD::numberOfPrecisionHoleLayers);
}

getNMuonPrecisionLayers()

unsigned IParticleHandle_TrackParticle::getNMuonPrecisionLayers ( ) const

virtual

Definition at line 482 of file IParticleHandle_TrackParticle.cxx.

{ 
  return summaryValue(xAOD::numberOfPrecisionLayers);
}

getNPixelHits()

unsigned IParticleHandle_TrackParticle::getNPixelHits ( ) const

virtual

Definition at line 466 of file IParticleHandle_TrackParticle.cxx.

{ 
  return summaryValue(xAOD::numberOfPixelHits);
}

getNSCTHits()

unsigned IParticleHandle_TrackParticle::getNSCTHits ( ) const

virtual

Definition at line 472 of file IParticleHandle_TrackParticle.cxx.

{ 
  return summaryValue(xAOD::numberOfSCTHits);
}

getNTRTHits()

unsigned IParticleHandle_TrackParticle::getNTRTHits ( ) const

virtual

Definition at line 477 of file IParticleHandle_TrackParticle.cxx.

{ 
  return summaryValue(xAOD::numberOfTRTHits);
}

has3DObjects()

bool IParticleHandle_TrackParticle::has3DObjects ( )

virtual

Returns true if the 3D objects have been created.

Reimplemented from AODHandleBase.

Definition at line 80 of file IParticleHandle_TrackParticle.cxx.

                                                {
  return 0!=m_d->sep;
}

hasCharge()

bool IParticleHandleBase::hasCharge ( ) const

inlineinherited

Definition at line 59 of file IParticleHandleBase.h.

{ return charge()!=unknown(); }

hasMass()

bool IParticleHandleBase::hasMass ( ) const

inlineinherited

Definition at line 60 of file IParticleHandleBase.h.

{ return mass()!=unknown(); }

iParticle()

const xAOD::IParticle & IParticleHandle_TrackParticle::iParticle ( ) const

virtual

Implements IParticleHandleBase.

Definition at line 442 of file IParticleHandle_TrackParticle.cxx.

{
  return *(m_d->trackparticle);
}

linear()

void IParticleHandle_TrackParticle::linear	(	Amg::Vector3D &	output,
		Amg::Vector3D &	a,
		Amg::Vector3D &	b,
		float	t )

Definition at line 213 of file IParticleHandle_TrackParticle.cxx.

                                                                                                      {
  float x = a.x() + (b.x()-a.x())*t;
  float y = a.y() + (b.y()-a.y())*t;
  float z = a.z() + (b.z()-a.z())*t;
  output = Amg::Vector3D(x,y,z);
}

mass()

double IParticleHandleBase::mass ( ) const

virtualinherited

Returns unknown() in case of trouble.

Definition at line 78 of file IParticleHandleBase.cxx.

{
  return 0.0;
}

momentum()

Amg::Vector3D IParticleHandle_TrackParticle::momentum ( ) const

virtual

Reimplemented from IParticleHandleBase.

Definition at line 407 of file IParticleHandle_TrackParticle.cxx.

  {
    const Trk::Perigee& p = m_d->trackparticle->perigeeParameters();
    return p.momentum();
  }

nodes()

SoNode * IParticleHandle_TrackParticle::nodes ( )

virtual

Returns the 3Dobjects.

Implements AODHandleBase.

Definition at line 100 of file IParticleHandle_TrackParticle.cxx.

                                            {
  VP1Msg::messageVerbose("IParticleHandle_TrackParticle::nodes with d->sep="+VP1Msg::str(m_d->sep));
 
  if (m_d->sep) {
    VP1Msg::messageVerbose("IParticleHandle_TrackParticle::nodes() - already have sep so not doing anything.");
    return m_d->sep; 
  }
  if (!m_d->sep) {
    m_d->sep = new SoSeparator();
    m_d->sep->ref();
  }
 
  bool fromTrack=m_d->collHandle->collSettingsButton().useExistingParameters();
  
  if (fromTrack) {
    addLine_FromTrackParticle();
  } else {
    #ifndef BUILDVP1LIGHT
      addLine_Extrapolated();
    #endif // BUILDVP1LIGHT
  }
 
  if (m_d->collHandle->collSettingsButton().showParameters()) 
      addParameterShapes();  
 
  return m_d->sep;
}

numberOfInstances()

int AODHandleBase::numberOfInstances ( )

staticinherited

Definition at line 87 of file AODHandleBase.cxx.

{
  return Imp::nobjhandles;
}

parametersAndPositions()

const QList< std::pair< xAOD::ParameterPosition, Amg::Vector3D > > & IParticleHandle_TrackParticle::parametersAndPositions

(

)

const

Definition at line 571 of file IParticleHandle_TrackParticle.cxx.

                                                                                                                       {
  // 1st position
  if (!m_d->parametersAndPositions.isEmpty()) 
    return m_d->parametersAndPositions;
  
  typedef std::pair<xAOD::ParameterPosition, Amg::Vector3D> paramAndPos;
  
#if defined BUILDVP1LIGHT
  m_d->parametersAndPositions.append(paramAndPos(xAOD::BeamLine, position() ) );
#else
  const Trk::Perigee& peri = m_d->trackparticle->perigeeParameters (); // FIXME - I'd quite like not to use anything which requires Athena ...
  m_d->parametersAndPositions.append(paramAndPos(xAOD::BeamLine, Amg::Vector3D(peri.position().x(),peri.position().y(),peri.position().z()) ) );
#endif // BUILDVP1LIGHT
 
  float x,y,z;
  for (unsigned int i=0; i<m_d->trackparticle->numberOfParameters() ; ++i){
    x = m_d->trackparticle->parameterX(i);
    y = m_d->trackparticle->parameterY(i);
    z = m_d->trackparticle->parameterZ(i);
  
    m_d->parametersAndPositions.append(paramAndPos(m_d->trackparticle->parameterPosition(i), Amg::Vector3D(x,y,z)));
  }
  VP1Msg::messageVerbose("TrackParticle parametersAndPositions() - returning this many positions : "+VP1Msg::str(m_d->parametersAndPositions.size()) );
  return m_d->parametersAndPositions;
}

position()

Amg::Vector3D IParticleHandle_TrackParticle::position ( ) const

virtual

Definition at line 433 of file IParticleHandle_TrackParticle.cxx.

  {
    const Trk::Perigee& p = m_d->trackparticle->perigeeParameters();
    return p.position(); // TODO: check what this returns
  }

rebuild3DObjects()

void AODHandleBase::rebuild3DObjects ( )

virtualinherited

Definition at line 137 of file AODHandleBase.cxx.

{  
  QString tmp2 = (m_visible==true)?"True":"False";
  VP1Msg::messageDebug("AODHandleBase::rebuild3DObjects() - with m_visible="+tmp2 );
  
  //Ensure we are always detached while updating.
  if (m_visible) {
    detach3DObjects();
  }
 
  clear3DObjects();
  nodes(); // rebuilds the 3D shapes
 
  //Attach if visible:
  if (m_visible) {
    attach3DObjects();
  }
}

setVisible()

void AODHandleBase::setVisible ( bool vis )

inherited

use by the collection handle.

Definition at line 99 of file AODHandleBase.cxx.

{
  QString tmp = (vis==true)?"True":"False";
  QString tmp2 = (m_visible==true)?"True":"False";
  VP1Msg::messageDebug(QString("AODHandleBase calling setVisible with vis=")+tmp+QString(", and m_visible=")+tmp2 );
  if (vis==m_visible)
    return;
 
  m_visible=vis;
  if (vis) {
    m_collhandle->incrementNShownHandles();
 
    if (!has3DObjects())
      rebuild3DObjects();//The call to rebuild also fixes attached state.
    else
      attach3DObjects();
  } else {
    m_collhandle->decrementNShownHandles();
    detach3DObjects();
  }
 
  visibleStateChanged();
}

shortInfo()

QString IParticleHandle_TrackParticle::shortInfo ( ) const

virtual

This returns the information shown about the object in the object browser.

Reimplemented from AODHandleBase.

Definition at line 502 of file IParticleHandle_TrackParticle.cxx.

{
  QString l("");
  //l+= "|P|="+VP1Msg::str(momentum().mag()/SYSTEM_OF_UNITS::GeV)+" [GeV], ";
  l+= "Pix["+QString::number(getNPixelHits())+"], ";
  l+= "SCT["+QString::number(getNSCTHits())+"], " ;
  l+= "TRT["+QString::number(getNTRTHits())+"], ";
  l+= "Muon prec. layers/holes [";
  l+= QString::number(getNMuonPrecisionLayers())+"/"+QString::number(getNMuonPrecisionHoleLayers())+"]";
  return l;
}

summaryValue()

unsigned IParticleHandle_TrackParticle::summaryValue ( xAOD::SummaryType type ) const

virtual

Definition at line 454 of file IParticleHandle_TrackParticle.cxx.

{
  uint8_t num = 0;
  if (m_d->trackparticle->summaryValue(num,type)){
    return num;
  }
  // else...
  VP1Msg::message("IParticleHandle_TrackParticle::getSummaryValue - unable to retrieve the requested enum: "+VP1Msg::str(type));
  return 999999;
}

type()

virtual QString IParticleHandle_TrackParticle::type ( ) const

inlinevirtual

return very short word with type (maybe link with collection type?)

Reimplemented from AODHandleBase.

Definition at line 65 of file IParticleHandle_TrackParticle.h.

unknown()

double AODHandleBase::unknown ( )

inlinestaticinherited

Definition at line 93 of file AODHandleBase.h.

{ return -999.0; }

update3DObjects()

void AODHandleBase::update3DObjects ( )

virtualinherited

Called after some configuration changes, or when the object is first created. Must be overloaded by children.

Definition at line 125 of file AODHandleBase.cxx.

{
  if (m_visible) {
    rebuild3DObjects();
  } else {    
    //Simply clear the present 3D objects. They will only be recreated if/when the track/jet becomes visible again.
    clear3DObjects();
  }
}

updateMaterial()

void AODHandleBase::updateMaterial ( )

inherited

Called after some configuration related to material changes.

Definition at line 208 of file AODHandleBase.cxx.

{
  VP1Msg::messageVerbose("AODHandleBase::updateMaterial()");
    
  //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_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 = determineMaterial();
    newmat->ref();
    if (newmat!=m_currentmaterial) {
      detach3DObjects();
      m_currentmaterial->unref();
      m_currentmaterial = newmat;
      attach3DObjects();
      currentMaterialChanged();
    } else {
      newmat->unref();
    }
  } else {
    //Just clear material.
    m_currentmaterial->unref();
    m_currentmaterial = 0;
    currentMaterialChanged();
  }
}

updateObjectBrowser()

void AODHandleBase::updateObjectBrowser ( )

inlineinherited

Update object browser QTreeWidgetItem.

Definition at line 102 of file AODHandleBase.h.

{};

visible()

bool AODHandleBase::visible ( ) const

inlineinherited

Definition at line 54 of file AODHandleBase.h.

{ return m_visible; }

visibleStateChanged()

void AODHandleBase::visibleStateChanged ( )

protectedvirtualinherited

Dummy implementation does nothing.

override if need to take action in this case.

Definition at line 316 of file AODHandleBase.cxx.

{  
  if ( !browserTreeItem()) {
    VP1Msg::messageVerbose("visibleStateChanged: No m_objBrowseTree!");
    return;
  }
    
  if (!visible()) {
    browserTreeItem()->setFlags(Qt::ItemFlags()); // not selectable, not enabled
  } else {
    browserTreeItem()->setFlags(Qt::ItemIsSelectable | Qt::ItemIsEnabled); //  selectable,  enabled
  }  
  QFont itemFont = browserTreeItem()->font(0);
  itemFont.setStrikeOut(!visible());
  browserTreeItem()->setFont(0, itemFont);
  browserTreeItem()->setFont(1, itemFont);
}

Member Data Documentation

m_collhandle

AODCollHandleBase* AODHandleBase::m_collhandle

privateinherited

Definition at line 126 of file AODHandleBase.h.

m_currentmaterial

SoMaterial* AODHandleBase::m_currentmaterial

privateinherited

Definition at line 127 of file AODHandleBase.h.

m_d

Imp* IParticleHandle_TrackParticle::m_d

private

Definition at line 80 of file IParticleHandle_TrackParticle.h.

m_visible

bool AODHandleBase::m_visible

privateinherited

Definition at line 125 of file AODHandleBase.h.

---

The documentation for this class was generated from the following files:

• IParticleHandle_TrackParticle.h
• IParticleHandle_TrackParticle.cxx