FCCAnalyses::VertexFinderLCFIPlus Namespace Reference

Primary and Seconday Vertex Finder interface using vertex fitter from VertexFitterSimple. More...

Functions
ROOT::VecOps::RVec< ROOT::VecOps::RVec< VertexingUtils::FCCAnalysesVertex > >	get_SV_jets (ROOT::VecOps::RVec< edm4hep::ReconstructedParticleData > recoparticles, ROOT::VecOps::RVec< edm4hep::TrackState > thetracks, VertexingUtils::FCCAnalysesVertex PV, ROOT::VecOps::RVec< bool > isInPrimary, ROOT::VecOps::RVec< fastjet::PseudoJet > jets, std::vector< std::vector< int > > jet_consti, bool V0_rej=true, double chi2_cut=9., double invM_cut=10., double chi2Tr_cut=5.)
	returns SVs reconstructed from non-primary tracks of jets non-primary separated from all tracks using isInPrimary (bool) vector currently not separating SVs by jet
ROOT::VecOps::RVec< VertexingUtils::FCCAnalysesVertex >	get_SV_event (ROOT::VecOps::RVec< edm4hep::ReconstructedParticleData > recoparticles, ROOT::VecOps::RVec< edm4hep::TrackState > thetracks, VertexingUtils::FCCAnalysesVertex PV, ROOT::VecOps::RVec< bool > isInPrimary, bool V0_rej=true, double chi2_cut=9., double invM_cut=10., double chi2Tr_cut=5.)
	returns SVs reconstructed from non-primary tracks of the event SV finding done before jet clustering non-primary separated from all tracks using isInPrimary (bool) vector
ROOT::VecOps::RVec< VertexingUtils::FCCAnalysesVertex >	get_SV_event (ROOT::VecOps::RVec< edm4hep::TrackState > np_tracks, ROOT::VecOps::RVec< edm4hep::TrackState > thetracks, VertexingUtils::FCCAnalysesVertex PV, bool V0_rej=true, double chi2_cut=9., double invM_cut=10., double chi2Tr_cut=5.)
	returns SVs reconstructed from non-primary tracks of the event SV finding done before jet clustering
ROOT::VecOps::RVec< int >	VertexSeed_best (ROOT::VecOps::RVec< edm4hep::TrackState > tracks, VertexingUtils::FCCAnalysesVertex PV, double chi2_cut=9., double invM_cut=10.)
	returns indices of the best pair of tracks from a vector of (non-primary) tracks default chi2 threshold is 9 and default invariant mass threshold is 10GeV
ROOT::VecOps::RVec< int >	addTrack_best (ROOT::VecOps::RVec< edm4hep::TrackState > tracks, ROOT::VecOps::RVec< int > vtx_tr, VertexingUtils::FCCAnalysesVertex PV, double chi2_cut=9., double invM_cut=10., double chi2Tr_cut=5.)
	adds index of the best track (from the remaining tracks) to the (seed) vtx default chi2 threshold is 9 and default invariant mass threshold is 10GeV default threshold for track's chi2 contribution is 5 (?)
ROOT::VecOps::RVec< edm4hep::TrackState >	V0rejection_tight (ROOT::VecOps::RVec< edm4hep::TrackState > tracks, VertexingUtils::FCCAnalysesVertex PV, bool V0_rej=true)
	V0 rejection (tight) takes all (non-primary tracks) & removes tracks coming from V0s if user chooses by default V0 rejection is done.
ROOT::VecOps::RVec< VertexingUtils::FCCAnalysesVertex >	findSVfromTracks (ROOT::VecOps::RVec< edm4hep::TrackState > tracks_fin, const ROOT::VecOps::RVec< edm4hep::TrackState > &alltracks, VertexingUtils::FCCAnalysesVertex PV, double chi2_cut=9., double invM_cut=10., double chi2Tr_cut=5.)
	find SVs from a set of tracks default values of thresholds for the constraints are set
bool	check_constraints (VertexingUtils::FCCAnalysesVertex vtx, ROOT::VecOps::RVec< edm4hep::TrackState > tracks, VertexingUtils::FCCAnalysesVertex PV, bool seed=true, double chi2_cut=9., double invM_cut=10., double chi2Tr_cut=5.)
	check constraints of vertex candidates default values of thresholds for the constraints are set default constraint check is that for finding vertex seed seed=true -> constraints for seed; seed=false -> constraints for adding tracks
ROOT::VecOps::RVec< bool >	isV0 (ROOT::VecOps::RVec< edm4hep::TrackState > np_tracks, VertexingUtils::FCCAnalysesVertex PV, bool tight=false)
	V0 rejection/identification takes all (non-primary) tracks & assigns "true" to pairs that form a V0 if(tight) -> tight constraints if(!tight) -> loose constraints by default loose constraints.
VertexingUtils::FCCAnalysesV0	get_V0s (ROOT::VecOps::RVec< edm4hep::TrackState > np_tracks, VertexingUtils::FCCAnalysesVertex PV, bool tight, double chi2_cut=9.)
	returns V0s reconstructed from a set of tracks (as an FCCAnalysesV0 object) constraint thresholds can be chosen out of two sets
VertexingUtils::FCCAnalysesV0	get_V0s (ROOT::VecOps::RVec< edm4hep::TrackState > np_tracks, VertexingUtils::FCCAnalysesVertex PV, double Ks_invM_low=0.493, double Ks_invM_high=0.503, double Ks_dis=0.5, double Ks_cosAng=0.999, double Lambda_invM_low=1.111, double Lambda_invM_high=1.121, double Lambda_dis=0.5, double Lambda_cosAng=0.99995, double Gamma_invM_low=0., double Gamma_invM_high=0.005, double Gamma_dis=9, double Gamma_cosAng=0.99995, double chi2_cut=9.)
	returns V0s reconstructed from a set of tracks (as an FCCAnalysesV0 object) constraint thresholds can be set manually
VertexingUtils::FCCAnalysesV0	get_V0s_jet (ROOT::VecOps::RVec< edm4hep::ReconstructedParticleData > recoparticles, ROOT::VecOps::RVec< edm4hep::TrackState > thetracks, ROOT::VecOps::RVec< bool > isInPrimary, ROOT::VecOps::RVec< fastjet::PseudoJet > jets, std::vector< std::vector< int > > jet_consti, VertexingUtils::FCCAnalysesVertex PV, bool tight=true, double chi2_cut=9.)
	returns V0s reconstructed in each jet of the event (as an FCCAnalysesV0 object) need to perform jet clustering before calling this function
ROOT::VecOps::RVec< double >	get_V0candidate (VertexingUtils::FCCAnalysesVertex &V0_vtx, ROOT::VecOps::RVec< edm4hep::TrackState > tr_pair, VertexingUtils::FCCAnalysesVertex PV, bool chi2, double chi2_cut=9.)
	returns invariant mass, distance from PV, and colliniarity variables for all V0 candidates [0] -> invM_Ks [GeV] [1] -> invM_Lambda1 [GeV] [2] -> invM_Lambda2 [GeV] [3] -> invM_Gamma [GeV] [4] -> r [mm] [5] -> r.p [unit vector] boolean check for if chi2 constraint needs to be checked skip the candidate with output size 0 - doesn't pass the chi2 cut
ROOT::VecOps::RVec< double >	constraints_Ks (bool tight)
	functions to fill constraint thresholds tight -> tight constraints !tight -> loose constraints also an option to choose constraint threshold
ROOT::VecOps::RVec< double >	constraints_Lambda0 (bool tight)
ROOT::VecOps::RVec< double >	constraints_Gamma (bool tight)
ROOT::VecOps::RVec< double >	constraints_Ks (double invM_low, double invM_high, double dis, double cosAng)
ROOT::VecOps::RVec< double >	constraints_Lambda0 (double invM_low, double invM_high, double dis, double cosAng)
ROOT::VecOps::RVec< double >	constraints_Gamma (double invM_low, double invM_high, double dis, double cosAng)

Variables
bool	debug_me = false
const double	m_pi = 0.13957039
const double	m_p = 0.93827208
const double	m_e = 0.00051099

Detailed Description

Primary and Seconday Vertex Finder interface using vertex fitter from VertexFitterSimple.

This represents a set functions and utilities to find vertices from a list of tracks following the algorithm from LCFIPlus framework.

Function Documentation

addTrack_best()

ROOT::VecOps::RVec< int > FCCAnalyses::VertexFinderLCFIPlus::addTrack_best	(	ROOT::VecOps::RVec< edm4hep::TrackState >	tracks,
		ROOT::VecOps::RVec< int >	vtx_tr,
		VertexingUtils::FCCAnalysesVertex	PV,
		double	chi2_cut = 9.,
		double	invM_cut = 10.,
		double	chi2Tr_cut = 5. )

adds index of the best track (from the remaining tracks) to the (seed) vtx default chi2 threshold is 9 and default invariant mass threshold is 10GeV default threshold for track's chi2 contribution is 5 (?)

check_constraints()

bool FCCAnalyses::VertexFinderLCFIPlus::check_constraints	(	VertexingUtils::FCCAnalysesVertex	vtx,
		ROOT::VecOps::RVec< edm4hep::TrackState >	tracks,
		VertexingUtils::FCCAnalysesVertex	PV,
		bool	seed = true,
		double	chi2_cut = 9.,
		double	invM_cut = 10.,
		double	chi2Tr_cut = 5. )

check constraints of vertex candidates default values of thresholds for the constraints are set default constraint check is that for finding vertex seed seed=true -> constraints for seed; seed=false -> constraints for adding tracks

constraints_Gamma() [1/2]

ROOT::VecOps::RVec< double > FCCAnalyses::VertexFinderLCFIPlus::constraints_Gamma

(

bool

tight

)

constraints_Gamma() [2/2]

ROOT::VecOps::RVec< double > FCCAnalyses::VertexFinderLCFIPlus::constraints_Gamma	(	double	invM_low,
		double	invM_high,
		double	dis,
		double	cosAng )

constraints_Ks() [1/2]

ROOT::VecOps::RVec< double > FCCAnalyses::VertexFinderLCFIPlus::constraints_Ks

(

bool

tight

)

functions to fill constraint thresholds tight -> tight constraints !tight -> loose constraints also an option to choose constraint threshold

[0] -> invariant mass lower limit [GeV] [1] -> invariant mass upper limit [GeV] [2] -> distance from PV [mm] [3] -> colinearity

constraints_Ks() [2/2]

ROOT::VecOps::RVec< double > FCCAnalyses::VertexFinderLCFIPlus::constraints_Ks	(	double	invM_low,
		double	invM_high,
		double	dis,
		double	cosAng )

constraints_Lambda0() [1/2]

ROOT::VecOps::RVec< double > FCCAnalyses::VertexFinderLCFIPlus::constraints_Lambda0

(

bool

tight

)

constraints_Lambda0() [2/2]

ROOT::VecOps::RVec< double > FCCAnalyses::VertexFinderLCFIPlus::constraints_Lambda0	(	double	invM_low,
		double	invM_high,
		double	dis,
		double	cosAng )

findSVfromTracks()

ROOT::VecOps::RVec< VertexingUtils::FCCAnalysesVertex > FCCAnalyses::VertexFinderLCFIPlus::findSVfromTracks	(	ROOT::VecOps::RVec< edm4hep::TrackState >	tracks_fin,
		const ROOT::VecOps::RVec< edm4hep::TrackState > &	alltracks,
		VertexingUtils::FCCAnalysesVertex	PV,
		double	chi2_cut = 9.,
		double	invM_cut = 10.,
		double	chi2Tr_cut = 5. )

find SVs from a set of tracks default values of thresholds for the constraints are set

get_SV_event() [1/2]

ROOT::VecOps::RVec< VertexingUtils::FCCAnalysesVertex > FCCAnalyses::VertexFinderLCFIPlus::get_SV_event	(	ROOT::VecOps::RVec< edm4hep::ReconstructedParticleData >	recoparticles,
		ROOT::VecOps::RVec< edm4hep::TrackState >	thetracks,
		VertexingUtils::FCCAnalysesVertex	PV,
		ROOT::VecOps::RVec< bool >	isInPrimary,
		bool	V0_rej = true,
		double	chi2_cut = 9.,
		double	invM_cut = 10.,
		double	chi2Tr_cut = 5. )

returns SVs reconstructed from non-primary tracks of the event SV finding done before jet clustering non-primary separated from all tracks using isInPrimary (bool) vector

get_SV_event() [2/2]

ROOT::VecOps::RVec< VertexingUtils::FCCAnalysesVertex > FCCAnalyses::VertexFinderLCFIPlus::get_SV_event	(	ROOT::VecOps::RVec< edm4hep::TrackState >	np_tracks,
		ROOT::VecOps::RVec< edm4hep::TrackState >	thetracks,
		VertexingUtils::FCCAnalysesVertex	PV,
		bool	V0_rej = true,
		double	chi2_cut = 9.,
		double	invM_cut = 10.,
		double	chi2Tr_cut = 5. )

returns SVs reconstructed from non-primary tracks of the event SV finding done before jet clustering

get_SV_jets()

ROOT::VecOps::RVec< ROOT::VecOps::RVec< VertexingUtils::FCCAnalysesVertex > > FCCAnalyses::VertexFinderLCFIPlus::get_SV_jets	(	ROOT::VecOps::RVec< edm4hep::ReconstructedParticleData >	recoparticles,
		ROOT::VecOps::RVec< edm4hep::TrackState >	thetracks,
		VertexingUtils::FCCAnalysesVertex	PV,
		ROOT::VecOps::RVec< bool >	isInPrimary,
		ROOT::VecOps::RVec< fastjet::PseudoJet >	jets,
		std::vector< std::vector< int > >	jet_consti,
		bool	V0_rej = true,
		double	chi2_cut = 9.,
		double	invM_cut = 10.,
		double	chi2Tr_cut = 5. )

returns SVs reconstructed from non-primary tracks of jets non-primary separated from all tracks using isInPrimary (bool) vector currently not separating SVs by jet

get_V0candidate()

ROOT::VecOps::RVec< double > FCCAnalyses::VertexFinderLCFIPlus::get_V0candidate	(	VertexingUtils::FCCAnalysesVertex &	V0_vtx,
		ROOT::VecOps::RVec< edm4hep::TrackState >	tr_pair,
		VertexingUtils::FCCAnalysesVertex	PV,
		bool	chi2,
		double	chi2_cut = 9. )

returns invariant mass, distance from PV, and colliniarity variables for all V0 candidates [0] -> invM_Ks [GeV] [1] -> invM_Lambda1 [GeV] [2] -> invM_Lambda2 [GeV] [3] -> invM_Gamma [GeV] [4] -> r [mm] [5] -> r.p [unit vector] boolean check for if chi2 constraint needs to be checked skip the candidate with output size 0 - doesn't pass the chi2 cut

get_V0s() [1/2]

VertexingUtils::FCCAnalysesV0 FCCAnalyses::VertexFinderLCFIPlus::get_V0s	(	ROOT::VecOps::RVec< edm4hep::TrackState >	np_tracks,
		VertexingUtils::FCCAnalysesVertex	PV,
		bool	tight,
		double	chi2_cut = 9. )

returns V0s reconstructed from a set of tracks (as an FCCAnalysesV0 object) constraint thresholds can be chosen out of two sets

get_V0s() [2/2]

VertexingUtils::FCCAnalysesV0 FCCAnalyses::VertexFinderLCFIPlus::get_V0s	(	ROOT::VecOps::RVec< edm4hep::TrackState >	np_tracks,
		VertexingUtils::FCCAnalysesVertex	PV,
		double	Ks_invM_low = 0.493,
		double	Ks_invM_high = 0.503,
		double	Ks_dis = 0.5,
		double	Ks_cosAng = 0.999,
		double	Lambda_invM_low = 1.111,
		double	Lambda_invM_high = 1.121,
		double	Lambda_dis = 0.5,
		double	Lambda_cosAng = 0.99995,
		double	Gamma_invM_low = 0.,
		double	Gamma_invM_high = 0.005,
		double	Gamma_dis = 9,
		double	Gamma_cosAng = 0.99995,
		double	chi2_cut = 9. )

returns V0s reconstructed from a set of tracks (as an FCCAnalysesV0 object) constraint thresholds can be set manually

get_V0s_jet()

VertexingUtils::FCCAnalysesV0 FCCAnalyses::VertexFinderLCFIPlus::get_V0s_jet	(	ROOT::VecOps::RVec< edm4hep::ReconstructedParticleData >	recoparticles,
		ROOT::VecOps::RVec< edm4hep::TrackState >	thetracks,
		ROOT::VecOps::RVec< bool >	isInPrimary,
		ROOT::VecOps::RVec< fastjet::PseudoJet >	jets,
		std::vector< std::vector< int > >	jet_consti,
		VertexingUtils::FCCAnalysesVertex	PV,
		bool	tight = true,
		double	chi2_cut = 9. )

returns V0s reconstructed in each jet of the event (as an FCCAnalysesV0 object) need to perform jet clustering before calling this function

isV0()

ROOT::VecOps::RVec< bool > FCCAnalyses::VertexFinderLCFIPlus::isV0	(	ROOT::VecOps::RVec< edm4hep::TrackState >	np_tracks,
		VertexingUtils::FCCAnalysesVertex	PV,
		bool	tight = false )

V0 rejection/identification takes all (non-primary) tracks & assigns "true" to pairs that form a V0 if(tight) -> tight constraints if(!tight) -> loose constraints by default loose constraints.

V0rejection_tight()

ROOT::VecOps::RVec< edm4hep::TrackState > FCCAnalyses::VertexFinderLCFIPlus::V0rejection_tight	(	ROOT::VecOps::RVec< edm4hep::TrackState >	tracks,
		VertexingUtils::FCCAnalysesVertex	PV,
		bool	V0_rej = true )

V0 rejection (tight) takes all (non-primary tracks) & removes tracks coming from V0s if user chooses by default V0 rejection is done.

VertexSeed_best()

ROOT::VecOps::RVec< int > FCCAnalyses::VertexFinderLCFIPlus::VertexSeed_best	(	ROOT::VecOps::RVec< edm4hep::TrackState >	tracks,
		VertexingUtils::FCCAnalysesVertex	PV,
		double	chi2_cut = 9.,
		double	invM_cut = 10. )

returns indices of the best pair of tracks from a vector of (non-primary) tracks default chi2 threshold is 9 and default invariant mass threshold is 10GeV

Variable Documentation

debug_me

bool FCCAnalyses::VertexFinderLCFIPlus::debug_me = false

m_e

const double FCCAnalyses::VertexFinderLCFIPlus::m_e = 0.00051099

m_p

const double FCCAnalyses::VertexFinderLCFIPlus::m_p = 0.93827208

m_pi

const double FCCAnalyses::VertexFinderLCFIPlus::m_pi = 0.13957039