Collaboration diagram for VP1TruthVertexCollection::Imp::VertexHandle:

Public Member Functions
	VertexHandle (HepMC::ConstGenVertexPtr v, VP1TruthVertexCollection::Imp *dd)

	~VertexHandle ()

void	recheckCutStatus ()

void	updateShape ()

HepMC::ConstGenVertexPtr	vertex () const

SoLineSet *	line () const

double	quantityVal (const VertexCommonFlags::QUANTITY &q)

Private Member Functions
void	ensureAttach (SoSeparator *collsep)

void	ensureDetach (SoSeparator *collsep)

bool	cut ()

Private Attributes
bool	m_attached

HepMC::ConstGenVertexPtr	m_vertex

SoLineSet *	m_line

VP1TruthVertexCollection::Imp *	m_d

Detailed Description

Definition at line 67 of file VP1TruthVertexCollection.cxx.

Constructor & Destructor Documentation

◆ VertexHandle()

VP1TruthVertexCollection::Imp::VertexHandle::VertexHandle	(	HepMC::ConstGenVertexPtr	v,
		VP1TruthVertexCollection::Imp *	dd
	)

inline

Definition at line 69 of file VP1TruthVertexCollection.cxx.

69 : m_attached(false), m_vertex(v), m_line(0),m_d(dd) {}

◆ ~VertexHandle()

VP1TruthVertexCollection::Imp::VertexHandle::~VertexHandle ( )

inline

Definition at line 70 of file VP1TruthVertexCollection.cxx.

70 { if (m_line) m_line->unref(); }

Member Function Documentation

◆ cut()

bool VP1TruthVertexCollection::Imp::VertexHandle::cut ( )

inlineprivate

Definition at line 157 of file VP1TruthVertexCollection.cxx.

                {
  
       if (m_d->controller->truthCutPrimaryVertexOnly()&&this!=m_d->vertices.at(0))
     return false;
  
       if (!m_d->controller->truthCutAllowedEta().contains(m_vertex->position().eta()))
     return false;
  
       QPair<VertexCommonFlags::QUANTITY,VP1Interval> p = m_d->controller->truthCutQuantity();
       if (!p.second.isAllR()) {
     if (!p.second.contains(quantityVal(p.first)))
       return false;
       }
  
       //We handle phi cut last:
       double phi(m_vertex->position().phi());
       for (const VP1Interval& i : m_d->controller->truthCutAllowedPhi()) {
     if (i.contains(phi)||i.contains(phi+2*M_PI)||i.contains(phi-2*M_PI))
       return true;
       }
       return false;
     }

◆ ensureAttach()

void VP1TruthVertexCollection::Imp::VertexHandle::ensureAttach ( SoSeparator * collsep )

inlineprivate

Definition at line 140 of file VP1TruthVertexCollection.cxx.

                                              {
       if (m_attached)
     return;
       m_attached = true;
       if (!m_line) {
     m_line = new SoLineSet;
     m_line->ref();
       }
       updateShape();
       collsep->addChild(m_line);
     }

◆ ensureDetach()

void VP1TruthVertexCollection::Imp::VertexHandle::ensureDetach ( SoSeparator * collsep )

inlineprivate

Definition at line 151 of file VP1TruthVertexCollection.cxx.

                                              {
       if (!m_attached)
     return;
       m_attached = false;
       collsep->removeChild(m_line);
     }

◆ line()

SoLineSet* VP1TruthVertexCollection::Imp::VertexHandle::line ( ) const

inline

Definition at line 113 of file VP1TruthVertexCollection.cxx.

113 { return m_line; }

◆ quantityVal()

double VP1TruthVertexCollection::Imp::VertexHandle::quantityVal ( const VertexCommonFlags::QUANTITY & q )

inline

Definition at line 115 of file VP1TruthVertexCollection.cxx.

                                                            {
       //Find total 4-momentum of incoming particles
       double px(0), py(0), pz(0), e(0);//HepMC::FourVector does NOT support operators + or += !!!!!!
 #ifdef HEPMC3
       for (const auto& PartIn:  m_vertex->particles_in()){
 #else
       HepMC::GenVertex::particles_in_const_iterator itPartIn,itPartInE(m_vertex->particles_in_const_end());
       for ( itPartIn = m_vertex->particles_in_const_begin();itPartIn!=itPartInE;++itPartIn) {
       auto PartIn=*itPartIn;
 #endif
     px +=PartIn->momentum().px();
     py +=PartIn->momentum().py();
     pz +=PartIn->momentum().pz();
     e += PartIn->momentum().e();
       }
       const HepMC::FourVector mom(px,py,pz,e);
       switch(q) {
       case VertexCommonFlags::ENERGY: return e;
       case VertexCommonFlags::MOMENTUM: return m_d->mag(mom);
       case VertexCommonFlags::TRANSVERSE_MOM: return mom.perp();
       default: // VertexCommonFlags::MASS:
     return mom.m();
       }
     }

◆ recheckCutStatus()

void VP1TruthVertexCollection::Imp::VertexHandle::recheckCutStatus ( )

inline

Definition at line 72 of file VP1TruthVertexCollection.cxx.

                             {
       bool cutval(cut());
       if (cutval!=m_attached) {
     if (cutval)
       ensureAttach(m_d->theclass->collSep());
     else
       ensureDetach(m_d->theclass->collSep());
       }
     }

◆ updateShape()

void VP1TruthVertexCollection::Imp::VertexHandle::updateShape ( )

inline

Definition at line 81 of file VP1TruthVertexCollection.cxx.

                        {
       if (!m_line)
     return;
       if (!m_attached) {
     //invalidate...
       }
  
       SoVertexProperty *vertices = static_cast<SoVertexProperty *>(m_line->vertexProperty.getValue());
       if (!vertices) {
     vertices = new SoVertexProperty();
     m_line->numVertices.set1Value(0,2);
     m_line->numVertices.set1Value(1,2);
     m_line->numVertices.set1Value(2,2);
     m_line->vertexProperty = vertices;
       }
  
       QPair<VertexCommonFlags::QUANTITY,double> p = m_d->controller->truthCrossLength();
       double extent;
       if (p.second<0)
     extent = -p.second;
       else
     extent = p.second * quantityVal(p.first);
  
       double x(m_vertex->position().x()), y(m_vertex->position().y()), z(m_vertex->position().z());
       vertices->vertex.set1Value ( 0,x-extent, y, z );
       vertices->vertex.set1Value ( 1,x+extent, y,z );
       vertices->vertex.set1Value ( 2, x,y-extent, z);
       vertices->vertex.set1Value ( 3, x,y+extent, z );
       vertices->vertex.set1Value ( 4, x, y,z-extent );
       vertices->vertex.set1Value ( 5, x, y,z+extent );
     }