#include <IParticleHandle_Electron.h>

Inheritance diagram for IParticleHandle_Electron:

Collaboration diagram for IParticleHandle_Electron:

Classes
class	Imp

Public Member Functions
	IParticleHandle_Electron (IParticleCollHandleBase , const xAOD::Electron )

virtual	~IParticleHandle_Electron ()

virtual bool	has3DObjects ()
	Returns true if the 3D objects have been created. More...

virtual void	clear3DObjects ()
	Delete objects. More...

virtual SoNode *	nodes ()
	Returns the 3Dobjects. More...

void	addParameterShapes ()

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

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

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. More...

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). More...

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. More...

const xAOD::Electron &	electron () const

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?) More...

QString	shortInfo () const
	This returns the information shown about the object in the object browser. More...

void	fillObjectBrowser (QList< QTreeWidgetItem * > &listOfItems)
	Create and fill the object browser QTreeWidgetItem. More...

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. More...

virtual QStringList	baseInfo () const

virtual double	mass () const
	Returns unknown() in case of trouble. More...

bool	hasCharge () const

bool	hasMass () const

const AODSysCommonData *	common () const

void	setVisible (bool)
	use by the collection handle. More...

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. More...

virtual void	rebuild3DObjects ()

virtual void	attach3DObjects ()
	Connect to the scene graph. More...

virtual void	detach3DObjects ()
	Remove from the scene graph. More...

void	updateMaterial ()
	Called after some configuration related to material changes. More...

SoMaterial *	currentMaterial () const
	Returns the current material of the handle. More...

void	updateObjectBrowser ()
	Update object browser QTreeWidgetItem. More...

QTreeWidgetItem *	browserTreeItem () const
	Return the QTreeWidgetItem;. More...

Static Public Member Functions
static int	numberOfInstances ()

static double	unknown ()

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

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_Electron.h.

Constructor & Destructor Documentation

◆ IParticleHandle_Electron()

IParticleHandle_Electron::IParticleHandle_Electron	(	IParticleCollHandleBase *	ch,
		const xAOD::Electron *	tp
	)

Definition at line 57 of file IParticleHandle_Electron.cxx.

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

◆ ~IParticleHandle_Electron()

IParticleHandle_Electron::~IParticleHandle_Electron ( )

virtual

Definition at line 68 of file IParticleHandle_Electron.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_Electron::addLine_Extrapolated ( )

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

Definition at line 299 of file IParticleHandle_Electron.cxx.

                                                    {
   VP1Msg::messageVerbose("IParticleHandle_Electron::addLine_Extrapolated().");
  
   const Trk::Perigee& peri = m_d->electron->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())); // TODO: move to position()?
   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_Electron::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_FromElectron()

void IParticleHandle_Electron::addLine_FromElectron ( )

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_Electron.cxx.

                                                    {
  
   std::vector<Amg::Vector3D> positions, momenta;
  
   const xAOD::TrackParticle* trackparticle = m_d->electron->trackParticle();
   if (!trackparticle) {
     VP1Msg::message("ERROR : no primary track particle associated with this Electron. Can't visualise it.");
     return;
   }
  
 #if defined BUILDVP1LIGHT
   positions.push_back(position());
   momenta.  push_back(momentum());
 #else
   const Trk::Perigee& peri = 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())); // TODO: move to position()?
   momenta.  push_back(Amg::Vector3D(peri.momentum().x(),peri.momentum().y(),peri.momentum().z()));
   // std::cout<<"i:"<<0<<"/"<<m_d->electron->numberOfParameters()+1<<": ("<<peri.position().x()<<","<<peri.position().y()<<","<<peri.position().z()<<")"<<std::endl;
 #endif // BUILDVP1LIGHT
  
   VP1Msg::messageVerbose("IParticleHandle_Electron::addLine_FromElectron - has "+QString::number(trackparticle->numberOfParameters())+" extra parameters.");
    for (unsigned int i=0; i<trackparticle->numberOfParameters() ; ++i){
     // std::cout<<"i:"<<i+1<<"/"<<m_d->electron->numberOfParameters()+1<<": ("<<m_d->electron->parameterX(i)<<","<<m_d->electron->parameterY(i)<<","<<m_d->electron->parameterZ(i)<<")"<<std::endl;
  
     positions.push_back(Amg::Vector3D(trackparticle->parameterX(i),
                                       trackparticle->parameterY(i),
                                       trackparticle->parameterZ(i)));
     momenta.push_back(Amg::Vector3D(  trackparticle->parameterPX(i),
                                       trackparticle->parameterPY(i),
                                       trackparticle->parameterPZ(i)));
   } // end of loop.
  
   // if ( positions.size()<2 ) VP1Msg::messageVerbose("IParticleHandle_Electron::addLine_FromElectron - 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_Electron::addParameterShapes ( )

Definition at line 128 of file IParticleHandle_Electron.cxx.

                                                  {
   // Points:
   SoPointSet* parameterPoints = 0;
   SoVertexProperty * vertices2 = 0;
  
   // std::cout<<"showing parameters... have "<<m_d->electron->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_Electron*>(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_Electron::bezier	(	Amg::Vector3D &	output,
		Amg::Vector3D &	a,
		Amg::Vector3D &	b,
		Amg::Vector3D &	c,
		Amg::Vector3D &	d,
		float	t
	)

Definition at line 224 of file IParticleHandle_Electron.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_Electron::charge ( ) const

virtual

Returns unknown() in case of trouble.

Implements IParticleHandleBase.

Definition at line 452 of file IParticleHandle_Electron.cxx.

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

◆ clear3DObjects()

void IParticleHandle_Electron::clear3DObjects ( )

virtual

Delete objects.

Implements AODHandleBase.

Definition at line 82 of file IParticleHandle_Electron.cxx.

                                              {
   VP1Msg::messageVerbose("IParticleHandle_Electron::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_Electron::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 371 of file IParticleHandle_Electron.cxx.

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

◆ collHandle()

const AODCollHandleBase* AODHandleBase::collHandle ( ) const

inlineinherited

Definition at line 55 of file AODHandleBase.h.

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

76 { return m_currentmaterial; }

◆ currentMaterialChanged()

virtual void AODHandleBase::currentMaterialChanged ( )

inlineprotectedvirtualinherited

Definition at line 111 of file AODHandleBase.h.

111 {};

◆ 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_Electron::dumpToJSON ( std::ofstream & str ) const

virtual

Reimplemented from AODHandleBase.

Definition at line 597 of file IParticleHandle_Electron.cxx.

                                                                  {
   str << "\"chi2\":"<<m_d->electron->trackParticle()->chiSquared () <<", ";
   str << "\"dof\":" <<m_d->electron->trackParticle()->numberDoF () <<", ";
   str << "\"dparams\": ["
       << m_d->electron->trackParticle()->d0() <<", "
       << m_d->electron->trackParticle()->z0() <<", "
       << m_d->electron->trackParticle()->phi0() <<", "
       << m_d->electron->trackParticle()->theta() <<", "
       << m_d->electron->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 << "] ";
 }

◆ electron()

const xAOD::Electron & IParticleHandle_Electron::electron ( ) const

Definition at line 75 of file IParticleHandle_Electron.cxx.

                                                              {
     return *(m_d->electron);
 }

◆ ensureInitSubSysHitInfo()

void IParticleHandle_Electron::ensureInitSubSysHitInfo ( ) const

private

◆ fillLineFromSplineFit()

void IParticleHandle_Electron::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 235 of file IParticleHandle_Electron.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_Electron::fillObjectBrowser ( QList< QTreeWidgetItem * > & list )

virtual

Create and fill the object browser QTreeWidgetItem.

Reimplemented from AODHandleBase.

Definition at line 514 of file IParticleHandle_Electron.cxx.

 {
   IParticleHandleBase::fillObjectBrowser(listOfItems); // Obligatory!
  
   QTreeWidgetItem* TSOSitem = new QTreeWidgetItem(browserTreeItem());
   TSOSitem->setText(0, QString("Def. Parameters " ) );
   QString dParameters("(");
   dParameters+=QString::number(m_d->electron->trackParticle()->d0());
   dParameters+=", ";
   dParameters+=QString::number(m_d->electron->trackParticle()->z0());
   dParameters+=", ";
   dParameters+=QString::number(m_d->electron->trackParticle()->phi0());
   dParameters+=", ";
   dParameters+=QString::number(m_d->electron->trackParticle()->theta());
   dParameters+=", ";
   dParameters+=QString::number(m_d->electron->trackParticle()->qOverP());
   dParameters+=")";
   TSOSitem->setText(1, dParameters );
  
   for (unsigned int i=0; i<m_d->electron->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->electron->trackParticle()->parameterX(i));
     pos+=", ";
     pos+=QString::number(m_d->electron->trackParticle()->parameterY(i));
     pos+=", ";
     pos+=QString::number(m_d->electron->trackParticle()->parameterZ(i));
     pos+=")";
  
     switch (m_d->electron->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_Electron::getNMuonPhiHoleLayers ( ) const

virtual

Definition at line 501 of file IParticleHandle_Electron.cxx.

 {
   return summaryValue(xAOD::numberOfPhiHoleLayers);
 }

◆ getNMuonPhiLayers()

unsigned IParticleHandle_Electron::getNMuonPhiLayers ( ) const

virtual

Definition at line 496 of file IParticleHandle_Electron.cxx.

 {
   return summaryValue(xAOD::numberOfPhiLayers);
  
 }//____________________________________________________________________

◆ getNMuonPrecisionHoleLayers()

unsigned IParticleHandle_Electron::getNMuonPrecisionHoleLayers ( ) const

virtual

Definition at line 491 of file IParticleHandle_Electron.cxx.

 {
   return summaryValue(xAOD::numberOfPrecisionHoleLayers);
 }

◆ getNMuonPrecisionLayers()

unsigned IParticleHandle_Electron::getNMuonPrecisionLayers ( ) const

virtual

Definition at line 486 of file IParticleHandle_Electron.cxx.

 {
   return summaryValue(xAOD::numberOfPrecisionLayers);
 }

◆ getNPixelHits()

unsigned IParticleHandle_Electron::getNPixelHits ( ) const

virtual

Definition at line 470 of file IParticleHandle_Electron.cxx.

 {
   return summaryValue(xAOD::numberOfPixelHits);
 }

◆ getNSCTHits()

unsigned IParticleHandle_Electron::getNSCTHits ( ) const

virtual

Definition at line 476 of file IParticleHandle_Electron.cxx.

 {
   return summaryValue(xAOD::numberOfSCTHits);
 }

◆ getNTRTHits()

unsigned IParticleHandle_Electron::getNTRTHits ( ) const

virtual

Definition at line 481 of file IParticleHandle_Electron.cxx.

 {
   return summaryValue(xAOD::numberOfTRTHits);
 }

◆ has3DObjects()

bool IParticleHandle_Electron::has3DObjects ( )

virtual

Returns true if the 3D objects have been created.

Reimplemented from AODHandleBase.

Definition at line 79 of file IParticleHandle_Electron.cxx.

                                            {
   return 0!=m_d->sep;
 }

◆ hasCharge()

bool IParticleHandleBase::hasCharge ( ) const

inlineinherited

Definition at line 59 of file IParticleHandleBase.h.

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

◆ hasMass()

bool IParticleHandleBase::hasMass ( ) const

inlineinherited

Definition at line 60 of file IParticleHandleBase.h.

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

◆ iParticle()

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

virtual

Implements IParticleHandleBase.

Definition at line 446 of file IParticleHandle_Electron.cxx.

 {
   return *(m_d->electron);
 }

◆ linear()

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

Definition at line 217 of file IParticleHandle_Electron.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_Electron::momentum ( ) const

virtual

Reimplemented from IParticleHandleBase.

Definition at line 411 of file IParticleHandle_Electron.cxx.

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

◆ nodes()

SoNode * IParticleHandle_Electron::nodes ( )

virtual

Returns the 3Dobjects.

Implements AODHandleBase.

Definition at line 99 of file IParticleHandle_Electron.cxx.

                                        {
   VP1Msg::messageVerbose("IParticleHandle_Electron::nodes with m_d->sep="+VP1Msg::str(m_d->sep));
  
   if (m_d->sep) {
     VP1Msg::messageVerbose("IParticleHandle_Electron::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_FromElectron();
   } else {
     #ifndef BUILDVP1LIGHT
       VP1Msg::messageVerbose( "addLine_Extrapolated()" );
       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_Electron::parametersAndPositions ( ) const

Definition at line 571 of file IParticleHandle_Electron.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->electron->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()) ) ); // TODO: move to position()?
 #endif // BUILDVP1LIGHT
  
   float x,y,z;
   for (unsigned int i=0; i<m_d->electron->trackParticle()->numberOfParameters() ; ++i){
     x = m_d->electron->trackParticle()->parameterX(i);
     y = m_d->electron->trackParticle()->parameterY(i);
     z = m_d->electron->trackParticle()->parameterZ(i);
  
     m_d->parametersAndPositions.append(paramAndPos(m_d->electron->trackParticle()->parameterPosition(i), Amg::Vector3D(x,y,z)));
   }
   VP1Msg::messageVerbose("Electron parametersAndPositions() - returning this many positions : "+VP1Msg::str(m_d->parametersAndPositions.size()) );
   return m_d->parametersAndPositions;
 }

◆ position()

Amg::Vector3D IParticleHandle_Electron::position ( ) const

virtual

Definition at line 436 of file IParticleHandle_Electron.cxx.

   {
     const Trk::Perigee& p = m_d->electron->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_Electron::shortInfo ( ) const

virtual

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

Reimplemented from AODHandleBase.

Definition at line 506 of file IParticleHandle_Electron.cxx.

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

◆ summaryValue()

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

virtual

Definition at line 458 of file IParticleHandle_Electron.cxx.

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

◆ type()

virtual QString IParticleHandle_Electron::type ( ) const

inlinevirtual

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

Reimplemented from AODHandleBase.

Definition at line 66 of file IParticleHandle_Electron.h.

◆ unknown()

static double AODHandleBase::unknown ( )

inlinestaticinherited

Definition at line 93 of file AODHandleBase.h.

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

102 {};

◆ visible()

bool AODHandleBase::visible ( ) const

inlineinherited

Definition at line 54 of file AODHandleBase.h.

54 { 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_Electron::m_d

private

Definition at line 81 of file IParticleHandle_Electron.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:

Classes

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ IParticleHandle_Electron()

◆ ~IParticleHandle_Electron()

Member Function Documentation

◆ addLine_Extrapolated()

◆ addLine_FromElectron()

◆ addParameterShapes()

◆ attach3DObjects()

◆ baseInfo()

◆ bezier()

◆ browserTreeItem()

◆ charge()

◆ clear3DObjects()

◆ clicked()

◆ collHandle()

◆ common()

◆ currentMaterial()

◆ currentMaterialChanged()

◆ detach3DObjects()

◆ determineMaterial()

◆ dumpToJSON()

◆ electron()

◆ ensureInitSubSysHitInfo()

◆ fillLineFromSplineFit()

◆ fillObjectBrowser()

◆ getNMuonPhiHoleLayers()

◆ getNMuonPhiLayers()

◆ getNMuonPrecisionHoleLayers()

◆ getNMuonPrecisionLayers()

◆ getNPixelHits()

◆ getNSCTHits()

◆ getNTRTHits()

◆ has3DObjects()

◆ hasCharge()

◆ hasMass()

◆ iParticle()

◆ linear()

◆ mass()

◆ momentum()

◆ nodes()

◆ numberOfInstances()

◆ parametersAndPositions()

◆ position()

◆ rebuild3DObjects()

◆ setVisible()

◆ shortInfo()

◆ summaryValue()

◆ type()

◆ unknown()

◆ update3DObjects()

◆ updateMaterial()

◆ updateObjectBrowser()

◆ visible()

◆ visibleStateChanged()

Member Data Documentation

◆ m_collhandle

◆ m_currentmaterial

◆ m_d

◆ m_visible