#include <TrkPriVxPurityTool.h>

Inheritance diagram for Trk::TrkPriVxPurityTool:

Collaboration diagram for Trk::TrkPriVxPurityTool:

Public Member Functions
StatusCode	initialize ()
	AlgTool method. More...

StatusCode	finalize ()
	AlgTool method. More...

	TrkPriVxPurityTool (const std::string &t, const std::string &n, const IInterface *p)
	AlgTool constructor. More...

virtual	~TrkPriVxPurityTool ()
	AlgTool destructor. More...

const TrkPriVxPurity *	purity (const Trk::VxCandidate *vertex) const
	Actual analysis method, returning a Trk::TrkPriVxPurity object for a given TrkVxCandidate. More...

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()
	AlgTool interface method. More...

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)

void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
double	m_r_tol
	Primary vertex definition for MC: tolerances taken in r (mm) around the simulated point of proton-proton interaction. More...

double	m_z_tol
	Primary vertex definition for MC: tolerances taken in z (mm) around the simulated point of proton-proton interaction. More...

std::string	m_trackParticleTruthCollName
	Name of the track true <-> rec map to be read from the storegate. More...

std::string	m_mc_collection_name
	Name of the MC event collection to read. More...

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

A tool estimating the purity of the reconstructed primary vertex. The Generator level tracks, corresponding to those fitted to the primary vertex candidate are examined. The fraction of outliers fitted to the vertex and tracks originating from the pileup is determined.

Author: Kirill Prokofiev, September 2006

Definition at line 30 of file TrkPriVxPurityTool.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ TrkPriVxPurityTool()

Trk::TrkPriVxPurityTool::TrkPriVxPurityTool	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

AlgTool constructor.

Definition at line 39 of file TrkPriVxPurityTool.cxx.

                                                                                                          : AthAlgTool ( t,n,p ) {
 //tolerances around the signal pp interaction  to create a primary vertex candidate:
 //defaul values are: sigma z = 1um (0.001 mm), sigma r = 1um (0.001 mm)
         declareProperty ( "VertexRTolerance", m_r_tol = 0.001 );
  
         declareProperty ( "VertexZTolerance", m_z_tol = 0.001 );
  
 //name of the rec <-> sim truth map to be read from the storegate
 //default one is the output of the  InDetTrackTruthMaker (see InDetTruthAlgs in InnerDetector packasge)
         declareProperty ( "SimTrackMapName", m_trackParticleTruthCollName = "TrackParticleTruthCollection" );
  
 //name of the Mc event collection. Required to get the signal event.
         declareProperty ( "MonteCarloCollection", m_mc_collection_name = "TruthEvent" );
  
 // declaring the interface fo the tool itself
         declareInterface<TrkPriVxPurityTool> ( this );
  
     }//end of constructor

◆ ~TrkPriVxPurityTool()

Trk::TrkPriVxPurityTool::~TrkPriVxPurityTool ( )

virtualdefault

AlgTool destructor.

Member Function Documentation

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyArrayType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
   }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	t,
		const SG::NotHandleType &
	)

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleBase &	hndl,
		const std::string &	doc,
		const SG::VarHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKey &	hndl,
		const std::string &	doc,
		const SG::VarHandleKeyType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKeyArray &	hndArr,
		const std::string &	doc,
		const SG::VarHandleKeyArrayType &
	)

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc,
		const SG::NotHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(name, property, doc);
   }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc = `"none"`
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

   {
     typedef typename SG::HandleClassifier<T>::type htype;
     return declareProperty (name, property, doc, htype());
   }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                  {
     typedef typename SG::HandleClassifier<T>::type htype;
     return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
   }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ finalize()

StatusCode Trk::TrkPriVxPurityTool::finalize ( )

AlgTool method.

Definition at line 34 of file TrkPriVxPurityTool.cxx.

                                             {
         msg(MSG::INFO)  << "Finalize successful" << endmsg;
         return StatusCode::SUCCESS;
     }

◆ initialize()

StatusCode Trk::TrkPriVxPurityTool::initialize ( )

AlgTool method.

Definition at line 30 of file TrkPriVxPurityTool.cxx.

                                               {
         return StatusCode::SUCCESS;
     }//end of initialize method

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

static const InterfaceID& Trk::TrkPriVxPurityTool::interfaceID ( )

inlinestatic

AlgTool interface method.

Definition at line 57 of file TrkPriVxPurityTool.h.

57 { return IID_TrkPriVxPurityTool; }

◆ msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ purity()

const TrkPriVxPurity * Trk::TrkPriVxPurityTool::purity ( const Trk::VxCandidate * vertex ) const

Actual analysis method, returning a Trk::TrkPriVxPurity object for a given TrkVxCandidate.

Definition at line 62 of file TrkPriVxPurityTool.cxx.

                                                                                             {
 //protection and related
         if ( vertex != nullptr ) {
             const std::vector<Trk::VxTrackAtVertex *> * tracks = vertex->vxTrackAtVertex();
             if ( tracks !=nullptr ) {
  
 // first getting the Mc Event collection
                 const McEventCollection * mcCollection ( nullptr );
                 StatusCode sc = evtStore()->retrieve ( mcCollection, m_mc_collection_name );
                 if ( sc.isFailure() ) {
                     if (msgLvl(MSG::DEBUG))                                      {
                         msg() << "Unable to retrieve MC collection: " << m_mc_collection_name << endmsg;
                         msg() << "Zero pointer returned." << endmsg;
                     }
                     return nullptr;
                 }
  
 //getting the signal event itself
                 McEventCollection::const_iterator it = mcCollection->begin();
                 const HepMC::GenEvent* genEvent= ( *it );
  
                 if( genEvent->vertices_empty() ) {
                   ATH_MSG_DEBUG( "No vertices found in first GenEvent" );
                   return nullptr;
                 }
 #ifdef HEPMC3
                auto pv = genEvent->vertices()[0];
 #else
                auto pv = *(genEvent->vertices_begin());
 #endif
  
 //analysing the MC event to create PV candidate
 //first finding the vertex of primary pp interaction
  
 //and storing its position
                 const auto& pv_pos = pv ->position();
                 double pv_r = pv_pos.perp();
                 double pv_z = pv_pos.z();
  
 // storing all the ids of vertices reasonably close to the primary one.
 // here the region of interest is selected.
                 std::map<int,HepMC::ConstGenVertexPtr> vertex_ids;
 #ifdef HEPMC3
                 for (const auto& vtx: genEvent->vertices()){
                     const auto& lv_pos = vtx->position();
                     if ( std::abs ( lv_pos.perp() - pv_r ) <m_r_tol  && std::abs ( lv_pos.z() - pv_z ) <m_z_tol ) {vertex_ids[vtx->id()] = vtx;}
                 }//end  of loop over all the vertices
 #else
                 for ( HepMC::GenEvent::vertex_const_iterator i = genEvent->vertices_begin(); i != genEvent->vertices_end()  ;++i ) {
                     auto vtx=*i;
                     const auto& lv_pos = vtx->position();
                     if ( std::abs ( lv_pos.perp() - pv_r ) <m_r_tol  && std::abs ( lv_pos.z() - pv_z ) <m_z_tol ) {vertex_ids[ HepMC::barcode(vtx) ]= vtx;}
                 }//end  of loop over all the vertices
 #endif
  
  
 //getting the track truth collection
                 const TrackParticleTruthCollection * trackParticleTruthCollection ( nullptr );
                 sc = evtStore()->retrieve ( trackParticleTruthCollection, m_trackParticleTruthCollName );
  
                 if ( sc.isFailure() ) {
                     if (msgLvl(MSG::DEBUG)) {
                         msg() << "Cannot retrieve " << m_trackParticleTruthCollName << endmsg;
                         msg() << "Zero pointer returned" << endmsg;
                     }
                     return nullptr;
                 }
  
 //looping over the tracks to find those matched to the GenParticle originating from signal PV
                 std::vector<Trk::VxTrackAtVertex *>::const_iterator vt = tracks->begin();
                 std::vector<Trk::VxTrackAtVertex *>::const_iterator ve = tracks->end();
                 unsigned int n_failed = 0;
                 std::vector<double> in_weights ( 0 );
                 std::vector<double> out_weights ( 0 );
                 std::vector<double> pu_weights ( 0 );
                 std::vector<double> no_correspondance ( 0 );
  
  
                 for ( ;vt!=ve;++vt ) {
 //original element link
 //ugly so far, we'll correct the EDM later
                     if ( ( *vt ) !=nullptr ) {
  
                         ITrackLink * origLink = ( **vt ).trackOrParticleLink();
  
                         if ( origLink !=nullptr ) {
                             // get to the original track particle
                             LinkToTrackParticleBase * tr_part = dynamic_cast< LinkToTrackParticleBase * > ( origLink );
                             if ( tr_part !=nullptr  && tr_part->isValid()) {
                 
  
                                 std::map< Rec::TrackParticleTruthKey, TrackParticleTruth>::const_iterator ttItr = trackParticleTruthCollection->end();
                   
                                 const Rec::TrackParticle* tp = dynamic_cast<const Rec::TrackParticle*>(**tr_part);
                                 if(tp) {
                                     const Rec::TrackParticleContainer* tpCont = dynamic_cast<const Rec::TrackParticleContainer*>(tr_part->getStorableObjectPointer());
                                     if (tpCont) {
                                         ElementLink<Rec::TrackParticleContainer> linkTruth;
                                         linkTruth.setElement(tp);
                                         linkTruth.setStorableObject(*tpCont);
                                         ttItr = trackParticleTruthCollection->find(Rec::TrackParticleTruthKey(linkTruth));
                                     }
                                 }
                 
  
                                 if (ttItr != trackParticleTruthCollection->end() ) {
                                     const HepMcParticleLink& particleLink = ttItr->second.particleLink();
 #ifdef HEPMC3
                                     HepMC::ConstGenParticlePtr genParticle = particleLink.scptr();
 #else
                                     HepMC::ConstGenParticlePtr genParticle = particleLink.cptr();
 #endif
  
                                     if(genParticle) {
                                         const auto *tpEvent = genParticle->parent_event();
                                         if(tpEvent==genEvent) { 
                                             HepMC::ConstGenVertexPtr pVertex{nullptr};
                                             if (genParticle) pVertex = genParticle->production_vertex();
                                             if ( pVertex) {
                                                 int link_pid = genParticle->pdg_id();
                                                 bool primary_track = false;
                                                 bool secondary_track = false;
                   
 //loop over the particles until decision is really taken
                                                 do {
 #ifdef HEPMC3
                                                     auto idf_res = vertex_ids.find ( pVertex->id() );
 #else
                                                     auto idf_res = vertex_ids.find ( HepMC::barcode(pVertex) );
 #endif
  
 //for the HepMcParticle Link, the signal event has an index 0.
 // tagging on it
                                                     if ( idf_res != vertex_ids.end() ) {
 //this is a primary track, storing its weight
                                                         primary_track = true;
                                                         in_weights.push_back ( ( *vt )->weight() );
                                                     }else {
 //this vertex is not from the central region.
 //checking whether it is a bremsstrahlung
 //if so, propagating track to its origin, otherwise rejecting it completely.
                                                         if ( pVertex->particles_in_size() == 1 ) {
 // one mother particle: is it a brem of some kind?
 #ifdef HEPMC3
                                                             auto inp = pVertex->particles_in()[0] ;
 #else
                                                             auto inp =*(pVertex->particles_in_const_begin()) ;
 #endif
                                                             auto tmpVertex_loc = inp ->production_vertex();
                                                             if ( inp ->pdg_id() == link_pid  && tmpVertex_loc) {
 // seems like a brem (this can be generator/simulation dependent unfortunately)
 // continue iterating
                                                                 pVertex = tmpVertex_loc;
                                                             }else {
                                                                 secondary_track = true;
                                                                 out_weights.push_back ( ( *vt )->weight() );
                                                             }//end of "no id change" check
                                                         }else {
  
 //has more than one mother particle; comes from reaction, not from central region
 // it is an outlier
                                                             secondary_track = true;
                                                             out_weights.push_back ( ( *vt )->weight() );
                                                         }//end of reaction check
                                                     }//end of central region check
                                                 }while ( !primary_track && !secondary_track );//end of loop over gen particles
                                             }else {
           
 // this track has no production vertex. Whatever it is, it does not
 // come from the primary interaction, consider as pileup
                                                 if (msgLvl(MSG::INFO)) {
                                                     msg() <<"A truth particle with no production vertex found."<<endmsg;
                                                     msg() <<"Since it does not come from PV, consider as PileUp"<<endmsg;
                                                 }
                                                 pu_weights.push_back ( ( *vt )->weight() );
                                             } //end of particle link without production vertex check.
                                         }else {
            
 //        std::cout<<"Not equal events "<< std::endl; 
                                             pu_weights.push_back ( ( *vt )->weight() );
                                         }//end of event pointers comparison
                                     }else{
          
 //the corresponding link to the track particles is zero. 
 //This seems to be a  fake (broken link in any case).           
                                         no_correspondance.push_back ( ( *vt )->weight() );
                                     }//end of genParticle check 
                                 }//end of search for correspondance in the trurth   collection
                             }else{
                                 if (msgLvl(MSG::DEBUG)) msg() <<"This track at vertex has no valid link to its original trackparticle "<<endmsg;
                                 ++ n_failed;
                             }//end of search for the truth link
                         }else{
                             if (msgLvl(MSG::DEBUG)) msg() <<"There is an empty pointer in the vector<VxTrackAtVertex *> for this vertex"<<endmsg;
                             ++n_failed;
                         }
                     }//end of valid link check
                 } //end of check of the empty pointer in the vector
  
 //debug
 //       std::cout<<"Total number of tracks fitted to the vertex: "<<tracks->size() <<std::endl;
 //       std::cout<<"Number of tracks w/o truth correspondance: "<<no_correspondance.size() <<std::endl;
 //       std::cout<<"Number of correct tracks: "<<total_size<<std::endl;
 //       std::cout<<"Number of tracks with broken origin links "<<n_failed<<std::endl;
 //       std::cout<<"Number of fitted outliers "<<out_weights.size() <<std::endl;
 //       std::cout<<"Number of fitted intliers "<<in_weights.size() <<std::endl;
       
                 return new Trk::TrkPriVxPurity ( tracks->size(), pu_weights, no_correspondance,
                                                  n_failed, in_weights, out_weights );    
             } // end of if tracks != 0
         }else{
             if (msgLvl(MSG::ERROR)) {
                 msg()<<"Empty vertex pointer received"<<endmsg;
                 msg()<<"Empty pointer returned."<<endmsg;
             }
             return nullptr;
         }//end of empty vertex check
         return nullptr;
     }//end of purity method

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

   {
     h.renounce();
     PBASE::renounce (h);
   }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                                {
     // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
     //         << "  size: " << m_vhka.size() << endmsg;
     for (auto &a : m_vhka) {
       std::vector<SG::VarHandleKey*> keys = a->keys();
       for (auto k : keys) {
         k->setOwner(this);
       }
     }
   }

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_mc_collection_name

std::string Trk::TrkPriVxPurityTool::m_mc_collection_name

private

Name of the MC event collection to read.

Definition at line 91 of file TrkPriVxPurityTool.h.

◆ m_r_tol

double Trk::TrkPriVxPurityTool::m_r_tol

private

Primary vertex definition for MC: tolerances taken in r (mm) around the simulated point of proton-proton interaction.

All the vertices in this region will be merged to form a single simulated PV with associated charged tracks.

Definition at line 73 of file TrkPriVxPurityTool.h.

◆ m_trackParticleTruthCollName

std::string Trk::TrkPriVxPurityTool::m_trackParticleTruthCollName

private

Name of the track true <-> rec map to be read from the storegate.

Definition at line 86 of file TrkPriVxPurityTool.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_z_tol

double Trk::TrkPriVxPurityTool::m_z_tol

private

Primary vertex definition for MC: tolerances taken in z (mm) around the simulated point of proton-proton interaction.

All the vertices in this region will be merged to form a single simulated PV with associated charged tracks.

Definition at line 81 of file TrkPriVxPurityTool.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ TrkPriVxPurityTool()

◆ ~TrkPriVxPurityTool()

Member Function Documentation

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ finalize()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ purity()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_detStore

◆ m_evtStore

◆ m_mc_collection_name

◆ m_r_tol

◆ m_trackParticleTruthCollName

◆ m_varHandleArraysDeclared

◆ m_vhka

◆ m_z_tol