FitFacetRelationVarSet Class Reference

Class to compute floating point variables from a facet relation which can be recorded as a flat TNtuple or serve as input to a MVA method. More...

#include <FitFacetRelationVarSet.h>

Inheritance diagram for FitFacetRelationVarSet:

Public Types
using	Object = typename Super::Object
	Type from which variables should be extracted.

Public Member Functions
bool	extract (const Relation< const CDCFacet > *ptrFacetRelation) final
	Generate and assign the contained variables.
std::vector< Named< Float_t * > >	getNamedVariables (const std::string &prefix) override
	Getter for the named references to the individual variables Base implementation returns empty vector.
std::vector< Named< Float_t * > >	getNamedVariables ()
	Getter for the named references to the individual variables.
MayBePtr< Float_t >	find (const std::string &varName) override
	Pointer to the variable with the given name.
virtual bool	extract (const Object *obj)
	Main method that extracts the variable values from the complex object.
bool	extract (const Object &obj)
	Method for extraction from an object instead of a pointer.
std::map< std::string, Float_t >	getNamedValues (const std::string &prefix) const
	Getter for a map of names to float values.
std::map< std::string, Float_t >	getNamedValues () const
	Getter for a map of names to float values.
void	initialize () override
	Receive and dispatch signal before the start of the event processing.
void	beginRun () override
	Receive and dispatch signal for the beginning of a new run.
void	beginEvent () override
	Receive and dispatch signal for the start of a new event.
void	endRun () override
	Receive and dispatch signal for the end of the run.
void	terminate () override
	Receive and dispatch Signal for termination of the event processing.

Protected Member Functions
Float_t	get () const
	Getter for the value of the ith variable. Static version.
Float_t &	var ()
	Reference getter for the value of the ith variable. Static version.
AssignFinite< Float_t >	finitevar ()
	Reference getter for the value of the ith variable. Transforms non-finite values to finite value.
void	addProcessingSignalListener (ProcessingSignalListener *psl)
	Register a processing signal listener to be notified.
int	getNProcessingSignalListener ()
	Get the number of currently registered listeners.

Static Protected Member Functions
static constexpr int	named (const char *name)
	Getter for the index from the name.

Private Types
using	Super = BaseVarSet< typename AVarNames::Object >
	Type of the super class.

Private Attributes
FixedSizeNamedFloatTuple< FitFacetRelationVarNames >	m_variables
	Memory for nVars floating point values.
std::vector< ProcessingSignalListener * >	m_subordinaryProcessingSignalListeners
	References to subordinary signal processing listener contained in this findlet.
bool	m_initialized = false
	Flag to keep track whether initialization happened before.
bool	m_terminated = false
	Flag to keep track whether termination happened before.
std::string	m_initializedAs
	Name of the type during initialisation.

Static Private Attributes
static const size_t	nVars
	Number of floating point values represented by this class.

Detailed Description

Class to compute floating point variables from a facet relation which can be recorded as a flat TNtuple or serve as input to a MVA method.

Definition at line 93 of file FitFacetRelationVarSet.h.

Member Typedef Documentation

Object

using Object = typename Super::Object

inherited

Type from which variables should be extracted.

Definition at line 44 of file VarSet.h.

Super

using Super = BaseVarSet<typename AVarNames::Object>

privateinherited

Type of the super class.

Definition at line 40 of file VarSet.h.

Member Function Documentation

addProcessingSignalListener()

void addProcessingSignalListener ( ProcessingSignalListener *  psl )

protectedinherited

Register a processing signal listener to be notified.

Definition at line 55 of file CompositeProcessingSignalListener.cc.

{
  m_subordinaryProcessingSignalListeners.push_back(psl);
}

beginEvent()

void beginEvent ( )

overridevirtualinherited

Receive and dispatch signal for the start of a new event.

Reimplemented from ProcessingSignalListener.

Reimplemented in SpacePointTagger< Belle2::CKFToPXDResult, Belle2::PXDCluster >, SpacePointTagger< Belle2::CKFToSVDResult, Belle2::SVDCluster >, BaseEventTimeExtractor< RecoTrack * >, BaseEventTimeExtractor< TrackFindingCDC::CDCWireHit & >, SharingHitsMatcher< Belle2::TrackFindingCDC::CDCTrack, Belle2::TrackFindingCDC::CDCSegment2D >, MCSymmetric< BaseAxialSegmentPairFilter >, MCSymmetric< BaseFacetFilter >, MCSymmetric< BaseFacetRelationFilter >, MCSymmetric< BaseSegmentPairFilter >, MCSymmetric< BaseSegmentPairRelationFilter >, MCSymmetric< BaseSegmentRelationFilter >, MCSymmetric< BaseSegmentTripleFilter >, MCSymmetric< BaseSegmentTripleRelationFilter >, MCSymmetric< BaseTrackRelationFilter >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCFacet >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCWireHit, true >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCSegment2D >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCTrack >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCSegmentPair >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCSegmentTriple >, RecoTrackStorer, ROIFinder, and SVDHoughTracking.

Definition at line 31 of file CompositeProcessingSignalListener.cc.

{
  Super::beginEvent();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->beginEvent();
  }
}

beginRun()

void beginRun ( )

overridevirtualinherited

Receive and dispatch signal for the beginning of a new run.

Reimplemented from ProcessingSignalListener.

Reimplemented in LayerRelationFilter< AFilter >, FourHitFilter, QualityIndicatorFilter, ThreeHitFilter, TwoHitVirtualIPFilter, TwoHitVirtualIPQIFilter, RecoTrackStorer, ROIFinder, SpacePointLoaderAndPreparer, and TrackCandidateResultRefiner.

Definition at line 23 of file CompositeProcessingSignalListener.cc.

{
  Super::beginRun();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->beginRun();
  }
}

endRun()

void endRun ( )

overridevirtualinherited

Receive and dispatch signal for the end of the run.

Reimplemented from ProcessingSignalListener.

Definition at line 39 of file CompositeProcessingSignalListener.cc.

{
  for (ProcessingSignalListener* psl : reversedRange(m_subordinaryProcessingSignalListeners)) {
    psl->endRun();
  }
  Super::endRun();
}

extract() [1/3]

bool extract ( const Object &  obj )

inlineinherited

Method for extraction from an object instead of a pointer.

Definition at line 56 of file BaseVarSet.h.

      {
        return extract(&obj);
      }

extract() [2/3]

virtual bool extract ( const Object *  obj )

inlinevirtualinherited

Main method that extracts the variable values from the complex object.

Returns indication whether the extraction could be completed successfully. Base implementation returns always true.

Parameters

obj	dummy object not used

Returns

: always true for the base implementation

Definition at line 50 of file BaseVarSet.h.

      {
        return true;
      }

extract() [3/3]

bool extract ( const Relation< const CDCFacet > *  ptrFacetRelation )

final

Generate and assign the contained variables.

Definition at line 27 of file FitFacetRelationVarSet.cc.

{
  if (not ptrFacetRelation) return false;
 
  const CDCFacet* fromFacet = ptrFacetRelation->first;
  const CDCFacet* toFacet   = ptrFacetRelation->second;
 
  const UncertainParameterLine2D& fromFitLine = fromFacet->getFitLine();
  const UncertainParameterLine2D& toFitLine   = toFacet->getFitLine();
 
  LineCovariance fromCov = fromFitLine.lineCovariance();
  LineParameters fromPar = fromFitLine.lineParameters();
 
  LineCovariance toCov   = toFitLine.lineCovariance();
  LineParameters toPar   = toFitLine.lineParameters();
 
  Vector2D fromTangential = fromFacet->getStartToEndLine().tangential();
  Vector2D toTangential   = toFacet->getStartToEndLine().tangential();
 
  Vector2D tangential = Vector2D::average(fromTangential, toTangential);
 
  double fromMiddleCos = fromFacet->getStartToMiddleLine().tangential().cosWith(toTangential);
  double toMiddleCos = fromTangential.cosWith(toFacet->getMiddleToEndLine().tangential());
 
  var<named("cos_delta")>() = fromTangential.cosWith(toTangential);
 
  var<named("from_middle_cos_delta")>() = fromMiddleCos;
  var<named("to_middle_cos_delta")>() = toMiddleCos;
 
  Vector2D frontWirePos2D = fromFacet->getStartWireHit().getRefPos2D();
  Vector2D backWirePos2D  = toFacet->getEndWireHit().getRefPos2D();
  {
    int nSteps = 0;
    UncertainParameterLine2D fitLine = FacetFitter::fit(*fromFacet, *toFacet, nSteps);
    double s = fitLine->lengthOnCurve(frontWirePos2D, backWirePos2D);
    double alpha = fitLine->support().angleWith(fitLine->tangential());
    var<named("alpha_0")>() = alpha;
    var<named("chi2_0")>() = fitLine.chi2();
    var<named("chi2_0_per_s")>() = fitLine.chi2() / s;
    var<named("erf_0")>() = std::erf(fitLine.chi2() / 800);
    var<named("fit_0_phi0")>() = fitLine->tangential().phi();
    var<named("fit_0_cos_delta")>() = fitLine->tangential().cosWith(tangential);
  }
 
  {
    int nSteps = 1;
    UncertainParameterLine2D fitLine = FacetFitter::fit(*fromFacet, *toFacet, nSteps);
    double s = fitLine->lengthOnCurve(frontWirePos2D, backWirePos2D);
    var<named("chi2_1")>() = fitLine.chi2();
    var<named("chi2_1_per_s")>() = fitLine.chi2() / s;
    var<named("fit_1_phi0")>() = fitLine->tangential().phi();
    var<named("fit_1_cos_delta")>() = fitLine->tangential().cosWith(tangential);
  }
 
  {
    UncertainParameterLine2D fitLine = FacetFitter::fit(*fromFacet, *toFacet);
    double s = fitLine->lengthOnCurve(frontWirePos2D, backWirePos2D);
    var<named("chi2")>() = fitLine.chi2();
    var<named("chi2_per_s")>() = fitLine.chi2() / s;
    var<named("fit_phi0")>() = fitLine->tangential().phi();
    var<named("fit_cos_delta")>() = fitLine->tangential().cosWith(tangential);
  }
 
  // Combination fit
  {
    using namespace NLineParameterIndices;
    {
      double phi0_var = fromFacet->getFitLine().lineCovariance()(c_Phi0, c_Phi0);
      if (not std::isfinite(phi0_var)) {
        B2INFO("from addr " << *fromFacet);
        B2INFO("From cov " << std::endl << fromFacet->getFitLine().lineCovariance());
        B2INFO("From cov " << std::endl << fromFitLine.lineCovariance());
      }
    }
    {
      double phi0_var = toFacet->getFitLine().lineCovariance()(c_Phi0, c_Phi0);
      if (not std::isfinite(phi0_var)) {
        B2INFO("to addr " << *toFacet);
        B2INFO("To cov " << std::endl << toFacet->getFitLine().lineCovariance());
        B2INFO("To cov " << std::endl << toFitLine.lineCovariance());
      }
    }
 
    var<named("phi0_from_sigma")>() = std::sqrt(fromCov(c_Phi0, c_Phi0));
    var<named("phi0_to_sigma")>() = std::sqrt(toCov(c_Phi0, c_Phi0));
    var<named("phi0_ref_sigma")>() = std::sqrt(fromCov(c_Phi0, c_Phi0) + toCov(c_Phi0, c_Phi0));
    var<named("phi0_ref_diff")>() = AngleUtil::normalised(toPar(c_Phi0) - fromPar(c_Phi0));
    var<named("phi0_ref_pull")>() =
      std::fabs(AngleUtil::normalised(toPar(c_Phi0) - fromPar(c_Phi0)) /
                std::sqrt((toCov(c_Phi0, c_Phi0) + fromCov(c_Phi0, c_Phi0))));
 
    LineParameters avgPar;
    LineCovariance avgCov;
    double chi2 = LineUtil::average(fromPar, fromCov, toPar, toCov, avgPar, avgCov);
 
    LineParameters meanPar = (fromPar + toPar) / 2.0;
    meanPar(c_Phi0) = AngleUtil::average(fromPar(c_Phi0), toPar(c_Phi0));
 
    LineParameters relAvgPar = avgPar - meanPar;
    relAvgPar(c_Phi0) = AngleUtil::normalised(relAvgPar(c_Phi0));
 
    var<named("chi2_comb")>() = chi2;
    var<named("phi0_comb_pull")>() =
      std::fabs(relAvgPar(c_Phi0) / std::sqrt(avgCov(c_Phi0, c_Phi0)));
    var<named("phi0_comb_diff")>() = relAvgPar(c_Phi0);
    var<named("phi0_comb_sigma")>() = std::sqrt(avgCov(c_Phi0, c_Phi0));
  }
 
  // Fitter
  {
    CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_Unit);
    observations2D.append(*fromFacet);
    observations2D.append(*toFacet);
 
    CDCTrajectory2D fittedTrajectory = CDCKarimakiFitter::getFitter().fit(std::move(observations2D));
    var<named("chi2_kari_unit")>() = fittedTrajectory.getChi2();
    var<named("abs_curv_unit")>() = std::fabs(fittedTrajectory.getCurvature());
  }
 
  {
    CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_DriftLength);
    observations2D.append(*fromFacet);
    observations2D.append(*toFacet);
 
    CDCTrajectory2D fittedTrajectory = CDCKarimakiFitter::getFitter().fit(std::move(observations2D));
    var<named("chi2_kari_l")>() = fittedTrajectory.getChi2();
    var<named("abs_curv_l")>() = std::fabs(fittedTrajectory.getCurvature());
  }
 
  {
    CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_Pseudo);
    observations2D.append(*fromFacet);
    observations2D.append(*toFacet);
 
    CDCTrajectory2D fittedTrajectory = CDCKarimakiFitter::getFitter().fit(std::move(observations2D));
    var<named("chi2_kari_pseudo")>() = fittedTrajectory.getChi2();
    var<named("abs_curv_pseudo")>() = std::fabs(fittedTrajectory.getCurvature());
  }
 
  {
    CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_Proper);
    observations2D.append(*fromFacet);
    observations2D.append(*toFacet);
 
    CDCTrajectory2D fittedTrajectory = CDCKarimakiFitter::getFitter().fit(std::move(observations2D));
    var<named("chi2_kari_proper")>() = fittedTrajectory.getChi2();
    var<named("abs_curv_proper")>() = std::fabs(fittedTrajectory.getCurvature());
  }
 
  return true;
}

find()

MayBePtr< Float_t > find ( const std::string &  varName )

inlineoverridevirtualinherited

Pointer to the variable with the given name.

Returns nullptr if not found.

Reimplemented from BaseVarSet< AVarNames::Object >.

Definition at line 64 of file VarSet.h.

      {
        return m_variables.find(varName);
      }

finitevar()

AssignFinite< Float_t > finitevar ( )

inlineprotectedinherited

Reference getter for the value of the ith variable. Transforms non-finite values to finite value.

Definition at line 130 of file VarSet.h.

      {
        static_assert(I < nVars, "Requested variable index exceeds number of variables.");
        return AssignFinite<Float_t>(m_variables[I]);
      }

get()

Float_t get ( ) const

inlineprotectedinherited

Getter for the value of the ith variable. Static version.

Definition at line 85 of file VarSet.h.

      {
        static_assert(I < nVars, "Requested variable index exceeds number of variables.");
        return m_variables.get(I);
      }

getNamedValues() [1/2]

std::map< std::string, Float_t > getNamedValues ( ) const

inlineinherited

Getter for a map of names to float values.

Definition at line 95 of file BaseVarSet.h.

      {
        const std::string prefix = "";
        return this->getNamedValues(prefix);
      }

getNamedValues() [2/2]

std::map< std::string, Float_t > getNamedValues ( const std::string &  prefix ) const

inlineinherited

Getter for a map of names to float values.

Parameters

prefix

Name prefix to apply to all variable names.

Definition at line 83 of file BaseVarSet.h.

      {
        std::map<std::string, Float_t> result;
        std::vector<Named<Float_t*> > namedVariables = this->getNamedVariables(prefix);
        for (const Named<Float_t*>& namedVariable : namedVariables) {
          Float_t* variable = namedVariable;
          result[namedVariable.getName()] = *variable;
        }
        return result;
      }

getNamedVariables() [1/2]

std::vector< Named< Float_t * > > getNamedVariables ( )

inlineinherited

Getter for the named references to the individual variables.

Definition at line 73 of file BaseVarSet.h.

      {
        const std::string prefix = "";
        return this->getNamedVariables(prefix);
      }

getNamedVariables() [2/2]

std::vector< Named< Float_t * > > getNamedVariables ( const std::string &  prefix )

inlineoverridevirtualinherited

Getter for the named references to the individual variables Base implementation returns empty vector.

Reimplemented from BaseVarSet< AVarNames::Object >.

Definition at line 55 of file VarSet.h.

      {
        return m_variables.getNamedVariables(prefix);
      }

getNProcessingSignalListener()

int getNProcessingSignalListener ( )

protectedinherited

Get the number of currently registered listeners.

Definition at line 60 of file CompositeProcessingSignalListener.cc.

{
  return m_subordinaryProcessingSignalListeners.size();
}

initialize()

void initialize ( )

overridevirtualinherited

Receive and dispatch signal before the start of the event processing.

Reimplemented from ProcessingSignalListener.

Reimplemented in UnionVarSet< AObject >, UnionVarSet< Object >, VariadicUnionVarSet< AVarSets >, ResultStorer< Belle2::CKFToPXDResult >, ResultStorer< Belle2::CKFToSVDResult >, BaseEventTimeExtractor< RecoTrack * >, BaseEventTimeExtractor< TrackFindingCDC::CDCWireHit & >, StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::HyperHough >, StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::QuadraticLegendre >, StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::Z0TanLambdaLegendre >, OnVarSet< Filter< ATruthVarSet::Object > >, OnVarSet< Filter< AVarSet::Object > >, OnVarSet< BaseFacetFilter >, OnVarSet< BaseFacetRelationFilter >, OnVarSet< BaseAxialSegmentPairFilter >, OnVarSet< BaseSegmentRelationFilter >, OnVarSet< BaseTrackRelationFilter >, OnVarSet< BaseSegmentPairRelationFilter >, MCSymmetric< BaseAxialSegmentPairFilter >, MCSymmetric< BaseFacetFilter >, MCSymmetric< BaseFacetRelationFilter >, MCSymmetric< BaseSegmentPairFilter >, MCSymmetric< BaseSegmentPairRelationFilter >, MCSymmetric< BaseSegmentRelationFilter >, MCSymmetric< BaseSegmentTripleFilter >, MCSymmetric< BaseSegmentTripleRelationFilter >, MCSymmetric< BaseTrackRelationFilter >, StoreArrayLoader< const Belle2::SpacePoint >, StoreArrayLoader< DataStoreInputTypeRefType >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCFacet >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCWireHit, true >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCSegment2D >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCTrack >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCSegmentPair >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCSegmentTriple >, RelationVarSet< ABaseVarSet >, QualityIndicatorFilter, TwoHitVirtualIPQIFilter, MultiHoughSpaceFastInterceptFinder, RawTrackCandCleaner< AHit >, RawTrackCandCleaner< Belle2::vxdHoughTracking::VXDHoughState >, RecoTrackStorer, ROIFinder, SingleHoughSpaceFastInterceptFinder, SpacePointLoaderAndPreparer, TrackCandidateOverlapResolver, and TrackCandidateResultRefiner.

Definition at line 15 of file CompositeProcessingSignalListener.cc.

{
  Super::initialize();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->initialize();
  }
}

named()

static constexpr int named ( const char *  name )

inlinestaticconstexprprotectedinherited

Getter for the index from the name.

Looks through the associated names and returns the right index if found Returns nVars (one after the last element) if not found.

Parameters

name	Name of the sought variable

Returns

Index of the name, nVars if not found.

Definition at line 78 of file VarSet.h.

      {
        return index<nVars>(AVarNames::getName, name);
      }

terminate()

void terminate ( )

overridevirtualinherited

Receive and dispatch Signal for termination of the event processing.

Reimplemented from ProcessingSignalListener.

Reimplemented in StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::HyperHough >, StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::QuadraticLegendre >, and StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::Z0TanLambdaLegendre >.

Definition at line 47 of file CompositeProcessingSignalListener.cc.

{
  for (ProcessingSignalListener* psl : reversedRange(m_subordinaryProcessingSignalListeners)) {
    psl->terminate();
  }
  Super::terminate();
}

var()

Float_t & var ( )

inlineprotectedinherited

Reference getter for the value of the ith variable. Static version.

Definition at line 93 of file VarSet.h.

      {
        static_assert(I < nVars, "Requested variable index exceeds number of variables.");
        return m_variables[I];
      }

Member Data Documentation

m_initialized

bool m_initialized = false

privateinherited

Flag to keep track whether initialization happened before.

Definition at line 52 of file ProcessingSignalListener.h.

m_initializedAs

std::string m_initializedAs

privateinherited

Name of the type during initialisation.

Definition at line 58 of file ProcessingSignalListener.h.

m_subordinaryProcessingSignalListeners

std::vector<ProcessingSignalListener*> m_subordinaryProcessingSignalListeners

privateinherited

References to subordinary signal processing listener contained in this findlet.

Definition at line 52 of file CompositeProcessingSignalListener.h.

m_terminated

bool m_terminated = false

privateinherited

Flag to keep track whether termination happened before.

Definition at line 55 of file ProcessingSignalListener.h.

m_variables

FixedSizeNamedFloatTuple<FitFacetRelationVarNames > m_variables

privateinherited

Memory for nVars floating point values.

Definition at line 138 of file VarSet.h.

nVars

const size_t nVars

staticprivateinherited

Number of floating point values represented by this class.

Definition at line 48 of file VarSet.h.

---

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

• tracking/trackFindingCDC/filters/facetRelation/include/FitFacetRelationVarSet.h
• tracking/trackFindingCDC/filters/facetRelation/src/FitFacetRelationVarSet.cc