Creator for VXDMeasurements with momentum estimation based on the dEdX information. More...

#include <VXDMomentumEstimationMeasurementCreator.h>

Inheritance diagram for VXDMomentumEstimationMeasurementCreator< HitType, detector >:

Public Member Functions
	VXDMomentumEstimationMeasurementCreator (const genfit::MeasurementFactory< genfit::AbsMeasurement > &measurementFactory)
	Constructor.

virtual	~VXDMomentumEstimationMeasurementCreator ()
	Destructor.

void	setParameter (const std::string &parameterName, const std::string &parameterValue) override
	Set the parameters of the fit functions and whether to use the thickness or not or the tracking seeds or not.

std::vector< genfit::TrackPoint * >	createMeasurementPoints (HitType *hit, RecoTrack &recoTrack, const RecoHitInformation &recoHitInformation) const override
	Create measurements based on coordinate measurements.

Protected Member Functions
virtual std::vector< genfit::AbsMeasurement * >	createMeasurementFromCoordinateMeasurement (HitType hit, const RecoTrack &recoTrack, const RecoHitInformation &, genfit::AbsMeasurement coordinateMeasurement) const override
	Create a measurement based on the momentum estimation given by the VXDMomentumEstimation class.

genfit::AbsMeasurement *	createCoordinateMeasurement (HitType *hit, const RecoHitInformation &recoHitInformation) const
	You probably need a coordinate measurement on which you can base your measurements in createMeasurements.

genfit::TrackPoint *	createTrackPointWithRecoHitInformation (genfit::AbsMeasurement *coordinateMeasurement, RecoTrack &recoTrack, const RecoHitInformation &recoHitInformation) const
	Helper: Create a TrackPoint from a measurement with a given RecoHitInformation.

Private Attributes
VXDMomentumEstimation< HitType >::FitParameters	m_fitParameters
	Parameters for the main function.

VXDMomentumEstimation< HitType >::CorrectionFitParameters	m_correctionFitParameters
	Parameters for the correction function.

bool	m_useMCInformation = false
	Use the seeds of the track finder or the seeds of the MC particles.

bool	m_useThickness = false
	Use the thickness of the clusters of the path length for estimating dX.

double	m_minimumMomentum = 0.1
	Minimal value for the momentum below the estimation is used.

double	m_sigma = 0.03
	Sigma of the measurement.

bool	m_useTrackingSeeds = false
	Use the tracking seeds in the origin for calculating the path length rather than the current state.

const genfit::MeasurementFactory< genfit::AbsMeasurement > &	m_measurementFactory
	A reference to the prefilled measurement factory.

Detailed Description

template<class HitType, Const::EDetector detector>
class Belle2::VXDMomentumEstimationMeasurementCreator< HitType, detector >

Creator for VXDMeasurements with momentum estimation based on the dEdX information.

Can use different possibilities to calculate the path length needed for the dEdX calculation. Needs the transformation fit parameters fitParameters and correctionFitParameters which can be calculated using the IPython notebook Analyse_Hitwise.ipynb.

Definition at line 27 of file VXDMomentumEstimationMeasurementCreator.h.

Constructor & Destructor Documentation

◆ VXDMomentumEstimationMeasurementCreator()

template<class HitType, Const::EDetector detector>

VXDMomentumEstimationMeasurementCreator ( const genfit::MeasurementFactory< genfit::AbsMeasurement > & measurementFactory )

inlineexplicit

Constructor.

Definition at line 30 of file VXDMomentumEstimationMeasurementCreator.h.

30 :

31 BaseMeasurementCreatorFromCoordinateMeasurement<HitType, detector>(measurementFactory) {}

◆ ~VXDMomentumEstimationMeasurementCreator()

template<class HitType, Const::EDetector detector>

virtual ~VXDMomentumEstimationMeasurementCreator ( )

inlinevirtual

Destructor.

Definition at line 34 of file VXDMomentumEstimationMeasurementCreator.h.

34{ }

Member Function Documentation

◆ createCoordinateMeasurement()

template<class HitType, Const::EDetector detector>

genfit::AbsMeasurement * createCoordinateMeasurement	(	HitType *	hit,
		const RecoHitInformation &	recoHitInformation ) const

protectedinherited

You probably need a coordinate measurement on which you can base your measurements in createMeasurements.

This function uses the measurementFactory to create one. Please be aware that this creates two new objects on the heap: the measurement and the track point. If you do not plan to use those two in a track, please delete them! The track point can be deleted by accessing the measurements GetTrackPoints function.

Parameters

hit
recoHitInformation

Returns: a coordinate AbsMeasurement as pointer.

Definition at line 21 of file BaseMeasurementCreatorFromHit.cc.

{
  // No one is using the trackCandHit (neither in SVDRecoHit nor in PXDRecoHit). So we do not provide any here!
  genfit::AbsMeasurement* coordinateMeasurement = m_measurementFactory.createOne(detector, hit->getArrayIndex(),
                                                  nullptr);
 
  return coordinateMeasurement;
}

◆ createMeasurementFromCoordinateMeasurement()

template<class HitType, Const::EDetector detector>

virtual std::vector< genfit::AbsMeasurement * > createMeasurementFromCoordinateMeasurement	(	HitType *	hit,
		const RecoTrack &	recoTrack,
		const RecoHitInformation &	,
		genfit::AbsMeasurement *	coordinateMeasurement ) const

inlineoverrideprotectedvirtual

Create a measurement based on the momentum estimation given by the VXDMomentumEstimation class.

Implements BaseMeasurementCreatorFromCoordinateMeasurement< HitType, detector >.

Definition at line 74 of file VXDMomentumEstimationMeasurementCreator.h.

    {
      genfit::PlanarMeasurement* planarMeasurement = dynamic_cast<genfit::PlanarMeasurement*>(coordinateMeasurement);
      if (planarMeasurement == nullptr) {
        B2FATAL("Can only add VXD hits which are based on PlanarMeasurements with momentum estimation!");
        // Make CPP check happy
        return {};
      }
 
      const ROOT::Math::XYZVector& momentum = recoTrack.getMomentumSeed();
 
      if (momentum.R() > m_minimumMomentum) {
        return {};
      }
 
      // cppcheck-suppress nullPointerRedundantCheck
      PlanarVXDMomentumMeasurement<HitType>* momentumMeasurement = new PlanarVXDMomentumMeasurement<HitType>(*planarMeasurement, hit,
          &recoTrack);
      momentumMeasurement->setCorrectionFitParameters(m_correctionFitParameters);
      momentumMeasurement->setFitParameters(m_fitParameters);
      momentumMeasurement->setSigma(m_sigma);
      momentumMeasurement->setUseMCInformation(m_useMCInformation);
      momentumMeasurement->setUseThickness(m_useThickness);
      momentumMeasurement->setUseTrackingSeeds(m_useTrackingSeeds);
      return {momentumMeasurement};
    }

◆ createMeasurementPoints()

template<class HitType, Const::EDetector detector>

std::vector< genfit::TrackPoint * > createMeasurementPoints	(	HitType *	hit,
		RecoTrack &	recoTrack,
		const RecoHitInformation &	recoHitInformation ) const

inlineoverridevirtualinherited

Create measurements based on coordinate measurements.

Implements BaseMeasurementCreatorFromHit< HitType, detector >.

Definition at line 33 of file BaseMeasurementCreatorFromCoordinateMeasurement.h.

    {
      genfit::AbsMeasurement* coordinateMeasurement = this->createCoordinateMeasurement(hit,
                                                      recoHitInformation);
 
      const std::vector<genfit::AbsMeasurement*>& measurements =
        this->createMeasurementFromCoordinateMeasurement(hit, recoTrack, recoHitInformation, coordinateMeasurement);
 
      // TODO: Do we want to create one track point for each or one track point for all?
      // TODO: This is especially important in the question to which the track point should be related!
      std::vector<genfit::TrackPoint*> trackPoints;
      trackPoints.reserve(measurements.size());
      for (genfit::AbsMeasurement* measurement : measurements) {
        trackPoints.push_back(this->createTrackPointWithRecoHitInformation(measurement, recoTrack, recoHitInformation));
      }
 
      return trackPoints;
    }

◆ createTrackPointWithRecoHitInformation()

template<class HitType, Const::EDetector detector>

genfit::TrackPoint * createTrackPointWithRecoHitInformation	(	genfit::AbsMeasurement *	coordinateMeasurement,
		RecoTrack &	recoTrack,
		const RecoHitInformation &	recoHitInformation ) const

protectedinherited

Helper: Create a TrackPoint from a measurement with a given RecoHitInformation.

Definition at line 32 of file BaseMeasurementCreatorFromHit.cc.

{
  genfit::TrackPoint* coordinateTrackPoint = new genfit::TrackPoint(coordinateMeasurement,
      &RecoTrackGenfitAccess::getGenfitTrack(recoTrack));
  coordinateTrackPoint->setSortingParameter(recoHitInformation.getSortingParameter());
 
  return coordinateTrackPoint;
}

◆ setParameter()

template<class HitType, Const::EDetector detector>

void setParameter	(	const std::string &	parameterName,
		const std::string &	parameterValue )

inlineoverridevirtual

Set the parameters of the fit functions and whether to use the thickness or not or the tracking seeds or not.

Reimplemented from BaseMeasurementCreator.

Definition at line 37 of file VXDMomentumEstimationMeasurementCreator.h.

    {
      if (parameterName == "minimumMomentum") {
        m_minimumMomentum = std::stod(parameterValue);
      } else if (parameterName == "aE") {
        m_fitParameters.aE = std::stod(parameterValue);
      } else if (parameterName == "bE") {
        m_fitParameters.bE = std::stod(parameterValue);
      } else if (parameterName == "cE") {
        m_fitParameters.cE = std::stod(parameterValue);
      } else if (parameterName == "dE") {
        m_fitParameters.dE = std::stod(parameterValue);
      } else if (parameterName == "aM") {
        m_correctionFitParameters.aM = std::stod(parameterValue);
      } else if (parameterName == "bM") {
        m_correctionFitParameters.bM = std::stod(parameterValue);
      } else if (parameterName == "cM") {
        m_correctionFitParameters.cM = std::stod(parameterValue);
      } else if (parameterName == "dM") {
        m_correctionFitParameters.dM = std::stod(parameterValue);
      } else if (parameterName == "useMCInformation") {
        m_useMCInformation = std::stoi(parameterValue);
      } else if (parameterName == "useThickness") {
        m_useThickness = std::stoi(parameterValue);
      } else if (parameterName == "sigma") {
        m_sigma = std::stod(parameterValue);
      } else if (parameterName == "useTrackingSeeds") {
        m_useTrackingSeeds = std::stoi(parameterValue);
      } else {
        B2FATAL("A parameter with the name " << parameterName << " and the value " << parameterValue << " could not be set.");
      }
    }

Member Data Documentation

◆ m_correctionFitParameters

template<class HitType, Const::EDetector detector>

VXDMomentumEstimation<HitType>::CorrectionFitParameters m_correctionFitParameters

private

Parameters for the correction function.

Set them to zero to not use a correction function.

Definition at line 107 of file VXDMomentumEstimationMeasurementCreator.h.

◆ m_fitParameters

template<class HitType, Const::EDetector detector>

VXDMomentumEstimation<HitType>::FitParameters m_fitParameters

private

Parameters for the main function.

Definition at line 105 of file VXDMomentumEstimationMeasurementCreator.h.

◆ m_measurementFactory

template<class HitType, Const::EDetector detector>

const genfit::MeasurementFactory<genfit::AbsMeasurement>& m_measurementFactory

privateinherited

A reference to the prefilled measurement factory.

Definition at line 72 of file BaseMeasurementCreatorFromHit.h.

◆ m_minimumMomentum

template<class HitType, Const::EDetector detector>

double m_minimumMomentum = 0.1

private

Minimal value for the momentum below the estimation is used.

Definition at line 113 of file VXDMomentumEstimationMeasurementCreator.h.

◆ m_sigma

template<class HitType, Const::EDetector detector>

double m_sigma = 0.03

private

Sigma of the measurement.

Definition at line 115 of file VXDMomentumEstimationMeasurementCreator.h.

◆ m_useMCInformation

template<class HitType, Const::EDetector detector>

bool m_useMCInformation = false

private

Use the seeds of the track finder or the seeds of the MC particles.

Definition at line 109 of file VXDMomentumEstimationMeasurementCreator.h.

◆ m_useThickness

template<class HitType, Const::EDetector detector>

bool m_useThickness = false

private

Use the thickness of the clusters of the path length for estimating dX.

Definition at line 111 of file VXDMomentumEstimationMeasurementCreator.h.

◆ m_useTrackingSeeds

template<class HitType, Const::EDetector detector>

bool m_useTrackingSeeds = false

private

Use the tracking seeds in the origin for calculating the path length rather than the current state.

Definition at line 117 of file VXDMomentumEstimationMeasurementCreator.h.

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

tracking/trackFitting/measurementCreator/creators/include/VXDMomentumEstimationMeasurementCreator.h

Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ VXDMomentumEstimationMeasurementCreator()

◆ ~VXDMomentumEstimationMeasurementCreator()

Member Function Documentation

◆ createCoordinateMeasurement()

◆ createMeasurementFromCoordinateMeasurement()

◆ createMeasurementPoints()

◆ createTrackPointWithRecoHitInformation()

◆ setParameter()

Member Data Documentation

◆ m_correctionFitParameters

◆ m_fitParameters

◆ m_measurementFactory

◆ m_minimumMomentum

◆ m_sigma

◆ m_useMCInformation

◆ m_useThickness

◆ m_useTrackingSeeds