Matrix inversion tools. More...

Namespaces
namespace	RKTools
	Array matrix multiplications used in RKTrackRep.

namespace	tools

Classes
class	AbsBField
	Abstract Interface to magnetic fields in GENFIT. More...

class	AbsFinitePlane
	Abstract base class for finite detector planes. More...

class	AbsFitter
	Abstract base class for fitters. More...

class	AbsFitterInfo
	This class collects all information needed and produced by a specific AbsFitter and is specific to one AbsTrackRep of the Track. More...

class	AbsHMatrix
	HMatrix for projecting from AbsTrackRep parameters to measured parameters in a DetPlane. More...

class	AbsKalmanFitter
	Abstract base class for Kalman fitter and derived fitting algorithms. More...

class	AbsMaterialInterface
	Abstract base class for geometry interfacing. More...

class	AbsMeasurement
	Contains the measurement and covariance in raw detector coordinates. More...

class	AbsMeasurementProducer
	Abstract interface class for MeasurementProducer. More...

class	AbsTrackRep
	Abstract base class for a track representation. More...

class	BellField
	Bell Field for SHiP. More...

class	ConstField
	Constant Magnetic field. More...

class	DAF
	Determinstic Annealing Filter (DAF) implementation. More...

class	DetPlane
	Detector plane. More...

class	EventDisplay
	Event display designed to run with Genfit. More...

class	Exception
	Exception class for error handling in GENFIT (provides storage for diagnostic information) More...

struct	ExtrapStep
	Helper for RKTrackRep. More...

class	FairShipFields
	Field for SHiP. More...

struct	fieldCache
	Cache B field at a position. Used by FieldManager. More...

class	FieldManager
	Singleton which provides access to magnetic field maps. More...

class	FitStatus
	Class where important numbers and properties of a fit can be stored. More...

class	FullMeasurement
	Measurement class implementing a measurement of all track parameters. More...

class	GFGbl
	Generic GBL implementation. More...

class	GFRaveMagneticField
	GFRaveMagneticField class Uses the FieldManager to provide a magnetic field to rave. More...

class	GFRavePropagator
	GFRavePropagator class. More...

class	GFRaveTrackParameters
	GFRaveTrackParameters class Contains a pointer to the original genfit::Track, the weight of the track in the vertex, and smoothed (with the vertex information) state and covariance of the track. More...

class	GFRaveVertex
	GFRaveVertex class. More...

class	GFRaveVertexFactory
	Vertex factory for producing GFRaveVertex objects from Track objects. More...

class	GoliathField
	Constant Magnetic field. More...

class	HelixTrackModel
	Helix track model for testing purposes. More...

class	HMatrixPhi
	AbsHMatrix implementation for one-dimensional MeasurementOnPlane and RKTrackRep parameterization. More...

class	HMatrixU
	AbsHMatrix implementation for one-dimensional MeasurementOnPlane and RKTrackRep parameterization. More...

class	HMatrixUnit
	AbsHMatrix implementation for 5-dimensional MeasurementOnPlane and RKTrackRep parameterization. More...

class	HMatrixUV
	AbsHMatrix implementation for two-dimensional MeasurementOnPlane and RKTrackRep parameterization. More...

class	HMatrixV
	AbsHMatrix implementation for one-dimensional MeasurementOnPlane and RKTrackRep parameterization. More...

class	KalmanFitStatus
	FitStatus for use with AbsKalmanFitter implementations. More...

class	KalmanFittedStateOnPlane
	MeasuredStateOnPlane with additional info produced by a Kalman filter or DAF. More...

class	KalmanFitter
	Simple Kalman filter implementation. More...

class	KalmanFitterInfo
	Collects information needed and produced by a AbsKalmanFitter implementations and is specific to one AbsTrackRep of the Track. More...

class	KalmanFitterRefTrack
	Kalman filter implementation with linearization around a reference track. More...

class	MaterialEffects
	Stepper and energy loss/noise matrix calculation. More...

class	MaterialProperties
	Material properties needed e.g. for material effects calculation. More...

struct	MatStep
	Simple struct containing MaterialProperties and stepsize in the material. More...

class	MeasuredStateOnPlane
	StateOnPlane with additional covariance matrix. More...

class	MeasurementCreator
	Create different measurement types along a HelixTrackModel for testing purposes. More...

class	MeasurementFactory
	Factory object to create AbsMeasurement objects from digitized and clustered data. More...

class	MeasurementOnPlane
	Measured coordinates on a plane. More...

class	MeasurementProducer
	Template class for a measurement producer module. More...

class	mySpacepointDetectorHit
	Example class for a spacepoint detector hit. More...

class	mySpacepointMeasurement
	Example class for a spacepoint measurement which can be created from mySpacepointDetectorHit via the MeasurementFactory. More...

class	PlanarMeasurement
	Measurement class implementing a planar hit geometry (1 or 2D). More...

class	ProlateSpacepointMeasurement
	Class for measurements implementing a space point hit geometry with a very prolate form of the covariance matrix. More...

class	RectangularFinitePlane
	Rectangular finite plane. More...

class	ReferenceStateOnPlane
	StateOnPlane with linearized transport to that ReferenceStateOnPlane from previous and next ReferenceStateOnPlane More...

struct	RKStep
	Helper for RKTrackRep. More...

class	RKTrackRep
	AbsTrackRep with 5D track parameterization in plane coordinates: (q/p, u', v', u, v) More...

class	SpacepointMeasurement
	Class for measurements implementing a space point hit geometry. More...

class	StateOnPlane
	A state with arbitrary dimension defined in a DetPlane. More...

class	StepLimits
	Helper to store different limits on the stepsize for the RKTRackRep. More...

class	TGeoMaterialInterface
	AbsMaterialInterface implementation for use with ROOT's TGeoManager. More...

class	ThinScatterer
	Thin or thick scatterer. More...

class	Track
	Collection of TrackPoint objects, AbsTrackRep objects and FitStatus objects. More...

struct	trackAndState
	Simple struct containing a Track pointer and a MeasuredStateOnPlane. Used in GFRave. More...

class	TrackCand
	Track candidate – seed values and indices. More...

class	TrackCandHit
	Hit object for use in TrackCand. Provides IDs and sorting parameters. More...

class	TrackPoint
	Object containing AbsMeasurement and AbsFitterInfo objects. More...

class	TrackPointComparator
	Helper class for TrackPoint sorting, used in Track::sort(). More...

class	WireMeasurement
	Class for measurements in wire detectors (Straw tubes and drift chambers) which do not measure the coordinate along the wire. More...

class	WirePointMeasurement
	Class for measurements in wire detectors (Straw tubes and drift chambers) which can measure the coordinate along the wire. More...

class	WireTrackCandHit
	Hit object for use in TrackCand. Provides additional left/right parameter. More...

Typedefs
typedef boost::shared_ptr< const genfit::MaterialProperties >	SharedMaterialPropertiesPtr

typedef boost::shared_ptr< genfit::DetPlane >	SharedPlanePtr
	Shared Pointer to a DetPlane.

typedef double	M1x3[1 *3]

typedef double	M1x4[1 *4]

typedef double	M1x6[1 *6]

typedef double	M1x7[1 *7]

typedef double	M5x5[5 *5]

typedef double	M6x6[6 *6]

typedef double	M7x7[7 *7]

typedef double	M8x7[8 *7]

typedef double	M6x5[6 *5]

typedef double	M7x5[7 *5]

typedef double	M5x6[5 *6]

typedef double	M5x7[5 *7]

Enumerations
enum	eFitterType { SimpleKalman , RefKalman , DafSimple , DafRef }

enum	eMultipleMeasurementHandling { weightedAverage , unweightedAverage , weightedClosestToReference , unweightedClosestToReference , weightedClosestToPrediction , unweightedClosestToPrediction , weightedClosestToReferenceWire , unweightedClosestToReferenceWire , weightedClosestToPredictionWire , unweightedClosestToPredictionWire }

enum	StepLimitType { stp_noLimit = 0 , stp_fieldCurv , stp_momLoss , stp_sMax , stp_sMaxArg , stp_boundary , stp_plane , ENUM_NR_ITEMS }

enum	eMeasurementType { Pixel = 0 , Spacepoint , ProlateSpacepoint , StripU , StripV , StripUV , Wire , WirePoint }

Functions
AbsFitterInfo *	new_clone (const AbsFitterInfo &a)
	Needed for boost cloneability:

MeasuredStateOnPlane	calcAverageState (const MeasuredStateOnPlane &forwardState, const MeasuredStateOnPlane &backwardState)
	Calculate weighted average between two MeasuredStateOnPlanes.

bool	operator== (const DetPlane &lhs, const DetPlane &rhs)

bool	operator!= (const DetPlane &lhs, const DetPlane &rhs)

bool	operator== (const MaterialProperties &lhs, const MaterialProperties &rhs)

bool	operator!= (const MaterialProperties &lhs, const MaterialProperties &rhs)

bool	operator== (const TrackCand &lhs, const TrackCand &rhs)

bool	operator== (const TrackCandHit &lhs, const TrackCandHit &rhs)

std::vector< rave::Track >	GFTracksToTracks (const std::vector< genfit::Track * > &GFTracks, std::vector< genfit::MeasuredStateOnPlane * > *GFStates, std::map< int, genfit::trackAndState > &IdGFTrackStateMap, int startID=0)

rave::Track	GFTrackToTrack (trackAndState, int id=-1, std::string tag="")

void	setData (const rave::Track &orig, MeasuredStateOnPlane *state)

GFRaveVertex *	RaveToGFVertex (const rave::Vertex &, const std::map< int, genfit::trackAndState > &IdGFTrackStateMap)

void	RaveToGFVertices (std::vector< GFRaveVertex * > *, const std::vector< rave::Vertex > &, const std::map< int, genfit::trackAndState > &IdGFTrackStateMap)

SharedPlanePtr	PlaneToGFDetPlane (const ravesurf::Plane &rplane)

TVector3	Point3DToTVector3 (const rave::Point3D &)

TVector3	Vector3DToTVector3 (const rave::Vector3D &)

TMatrixDSym	Covariance3DToTMatrixDSym (const rave::Covariance3D &)

TVectorD	Vector6DToTVectorD (const rave::Vector6D &)

TMatrixDSym	Covariance6DToTMatrixDSym (const rave::Covariance6D &)

rave::Point3D	TVector3ToPoint3D (const TVector3 &)

rave::Covariance3D	TMatrixDSymToCovariance3D (const TMatrixDSym &)

double	MeanExcEnergy_get (int Z)

double	MeanExcEnergy_get (TGeoMaterial *)

Variables
const int	MeanExcEnergy_NELEMENTS = 93

const double	MeanExcEnergy_vals [] = {1.e15, 19.2, 41.8, 40.0, 63.7, 76.0, 78., 82.0, 95.0, 115.0, 137.0, 149.0, 156.0, 166.0, 173.0, 173.0, 180.0, 174.0, 188.0, 190.0, 191.0, 216.0, 233.0, 245.0, 257.0, 272.0, 286.0, 297.0, 311.0, 322.0, 330.0, 334.0, 350.0, 347.0, 348.0, 343.0, 352.0, 363.0, 366.0, 379.0, 393.0, 417.0, 424.0, 428.0, 441.0, 449.0, 470.0, 470.0, 469.0, 488.0, 488.0, 487.0, 485.0, 491.0, 482.0, 488.0, 491.0, 501.0, 523.0, 535.0, 546.0, 560.0, 574.0, 580.0, 591.0, 614.0, 628.0, 650.0, 658.0, 674.0, 684.0, 694.0, 705.0, 718.0, 727.0, 736.0, 746.0, 757.0, 790.0, 790.0, 800.0, 810.0, 823.0, 823.0, 830.0, 825.0, 794.0, 827.0, 826.0, 841.0, 847.0, 878.0, 890.0}

Detailed Description

Matrix inversion tools.

Collection of converter functions.

Typedef Documentation

◆ M1x3

typedef double genfit::M1x3[1 *3]

Array Matrix typedefs. They are needed for SSE optimization: gcc can vectorize loops only if the array sizes are known.

Definition at line 33 of file RKTools.h.

◆ M1x4

typedef double genfit::M1x4[1 *4]

Definition at line 34 of file RKTools.h.

◆ M1x6

typedef double genfit::M1x6[1 *6]

Definition at line 35 of file RKTools.h.

◆ M1x7

typedef double genfit::M1x7[1 *7]

Definition at line 36 of file RKTools.h.

◆ M5x5

typedef double genfit::M5x5[5 *5]

Definition at line 37 of file RKTools.h.

◆ M5x6

typedef double genfit::M5x6[5 *6]

Definition at line 43 of file RKTools.h.

◆ M5x7

typedef double genfit::M5x7[5 *7]

Definition at line 44 of file RKTools.h.

◆ M6x5

typedef double genfit::M6x5[6 *5]

Definition at line 41 of file RKTools.h.

◆ M6x6

typedef double genfit::M6x6[6 *6]

Definition at line 38 of file RKTools.h.

◆ M7x5

typedef double genfit::M7x5[7 *5]

Definition at line 42 of file RKTools.h.

◆ M7x7

typedef double genfit::M7x7[7 *7]

Definition at line 39 of file RKTools.h.

◆ M8x7

typedef double genfit::M8x7[8 *7]

Definition at line 40 of file RKTools.h.

◆ SharedMaterialPropertiesPtr

typedef boost::shared_ptr< const genfit::MaterialProperties > genfit::SharedMaterialPropertiesPtr

Definition at line 38 of file SharedMaterialPropertiesPtr.h.

◆ SharedPlanePtr

typedef boost::shared_ptr< genfit::DetPlane > genfit::SharedPlanePtr

Shared Pointer to a DetPlane.

Ownership can be shared, e.g between multiple StateOnPlane objects. The DetPlane will automatically be deleted, if no owner remains.

Definition at line 43 of file SharedPlanePtr.h.

Enumeration Type Documentation

◆ eFitterType

enum genfit::eFitterType

Enumerator
SimpleKalman
RefKalman
DafSimple
DafRef

Definition at line 41 of file EventDisplay.h.

                 {
  SimpleKalman,
  RefKalman,
  DafSimple,
  DafRef
};

◆ eMeasurementType

enum genfit::eMeasurementType

Enumerator
Pixel
Spacepoint
ProlateSpacepoint
StripU
StripV
StripUV
Wire
WirePoint

Definition at line 42 of file MeasurementCreator.h.

                      { Pixel = 0,
        Spacepoint,
        ProlateSpacepoint,
        StripU,
        StripV,
        StripUV,
        Wire,
        WirePoint };

◆ eMultipleMeasurementHandling

enum genfit::eMultipleMeasurementHandling

Enumerator
weightedAverage	weighted average between measurements; used by DAF
unweightedAverage	average between measurements, all weighted with 1
weightedClosestToReference	closest to reference, weighted with its weight_
unweightedClosestToReference	closest to reference, weighted with 1
weightedClosestToPrediction	closest to prediction, weighted with its weight_
unweightedClosestToPrediction	closest to prediction, weighted with 1
weightedClosestToReferenceWire	if corresponding TrackPoint has one WireMeasurement, select closest to reference, weighted with its weight_. Otherwise use weightedAverage.
unweightedClosestToReferenceWire	if corresponding TrackPoint has one WireMeasurement, select closest to reference, weighted with 1. Otherwise use unweightedAverage.
weightedClosestToPredictionWire	if corresponding TrackPoint has one WireMeasurement, select closest to prediction, weighted with its weight_. Otherwise use weightedAverage.
unweightedClosestToPredictionWire	if corresponding TrackPoint has one WireMeasurement, select closest to prediction, weighted with 1. Otherwise use unweightedAverage. Recommended for KalmanFitter to 'resolve' l/r ambiguities

Definition at line 35 of file AbsKalmanFitter.h.

                                  {
  weightedAverage, 
  unweightedAverage, 
  weightedClosestToReference, 
  unweightedClosestToReference, 
  weightedClosestToPrediction, 
  unweightedClosestToPrediction, 
  weightedClosestToReferenceWire, 
  unweightedClosestToReferenceWire, 
  weightedClosestToPredictionWire, 
  unweightedClosestToPredictionWire 
};

◆ StepLimitType

enum genfit::StepLimitType

Enumerator
stp_noLimit
stp_fieldCurv
stp_momLoss
stp_sMax
stp_sMaxArg
stp_boundary
stp_plane
ENUM_NR_ITEMS

Definition at line 33 of file StepLimits.h.

                   {
  // soft limits (only rough estimation, can go beyond safely)
  stp_noLimit = 0,    // only for internal use
 
  // medium limits (can go a bit further if e.g. plane or boundary will be reached)
  stp_fieldCurv,  // stepsize limited by curvature and magnetic field inhomogenities
  stp_momLoss,    // stepsize limited by stepper because maximum momLoss is reached
  stp_sMax,       // stepsize limited by SMax defined in #estimateStep()
 
  // hard limits (must stop there at any case!)
  stp_sMaxArg,    // stepsize limited by argument maxStepArg passed to #estimateStep()
  stp_boundary,   // stepsize limited by stepper because material boundary is encountered
  stp_plane,      // stepsize limited because destination plane is reached
 
  ENUM_NR_ITEMS   // only for internal use
};

Function Documentation

◆ calcAverageState()

MeasuredStateOnPlane genfit::calcAverageState	(	const MeasuredStateOnPlane &	forwardState,
		const MeasuredStateOnPlane &	backwardState
	)

Calculate weighted average between two MeasuredStateOnPlanes.

Definition at line 65 of file MeasuredStateOnPlane.cc.

                                                                                                                           {
  // check if both states are defined in the same plane
  if (forwardState.getPlane() != backwardState.getPlane()) {
    Exception e("KalmanFitterInfo::calcAverageState: forwardState and backwardState are not defined in the same plane.", __LINE__,__FILE__);
    throw e;
  }
 
  // This code is called a lot, so some effort has gone into reducing
  // the number of temporary objects being constructed.
 
#if 0
  // For ease of understanding, here's a version of the code that is a
  // few percent slower, depending on the exact use-case, but which
  // makes the math more explicit:
  TMatrixDSym fCovInv, bCovInv, smoothed_cov;
  tools::invertMatrix(forwardState.getCov(), fCovInv);
  tools::invertMatrix(backwardState.getCov(), bCovInv);
 
  tools::invertMatrix(fCovInv + bCovInv, smoothed_cov);  // one temporary TMatrixDSym
 
  MeasuredStateOnPlane retVal(forwardState);
  retVal.setState(smoothed_cov*(fCovInv*forwardState.getState() + bCovInv*backwardState.getState())); // four temporary TVectorD's
  retVal.setCov(smoothed_cov);
  return retVal;
#endif
 
  static TMatrixDSym fCovInv, bCovInv;  // Static to avoid re-constructing for every call
  tools::invertMatrix(forwardState.getCov(), fCovInv);
  tools::invertMatrix(backwardState.getCov(), bCovInv);
 
  // Using a StateOnPlane here turned out at least as fast as
  // constructing a MeasuredStateOnPlane here, and then resetting its
  // covariance matrix below, even though it means another copy in the
  // return statement.
  StateOnPlane sop(forwardState); // copies auxInfo, plane, rep in the process
                                  // Using 'static' + subsequent assignment is measurably slower.
                                  // Surprisingly, using only TVectorD
                                  // sop(forwardState.getState()) with according
                                  // changes below measured slower.
 
  sop.getState() *= fCovInv;
  fCovInv += bCovInv;
  tools::invertMatrix(fCovInv);  // This is now the covariance of the average.
  sop.getState() += bCovInv*backwardState.getState();  // one temporary TVectorD
  sop.getState() *= fCovInv;
 
  return MeasuredStateOnPlane(sop, fCovInv);
}

◆ Covariance3DToTMatrixDSym()

TMatrixDSym genfit::Covariance3DToTMatrixDSym ( const rave::Covariance3D & ravecov )

Definition at line 243 of file GFRaveConverters.cc.

                                                          {
  TMatrixDSym cov(3);
 
  cov(0,0) = ravecov.dxx();
  cov(0,1) = ravecov.dxy();
  cov(0,2) = ravecov.dxz();
 
  cov(1,0) = ravecov.dxy();
  cov(1,1) = ravecov.dyy();
  cov(1,2) = ravecov.dyz();
 
  cov(2,0) = ravecov.dxz();
  cov(2,1) = ravecov.dyz();
  cov(2,2) = ravecov.dzz();
 
  return cov;
}

◆ Covariance6DToTMatrixDSym()

TMatrixDSym genfit::Covariance6DToTMatrixDSym ( const rave::Covariance6D & ravecov )

Definition at line 279 of file GFRaveConverters.cc.

                                                          {
  TMatrixDSym cov(6);
 
  cov(0,0) = ravecov.dxx();
  cov(0,1) = ravecov.dxy();
  cov(0,2) = ravecov.dxz();
  cov(0,3) = ravecov.dxpx();
  cov(0,4) = ravecov.dxpy();
  cov(0,5) = ravecov.dxpz();
 
  cov(1,0) = ravecov.dxy();
  cov(1,1) = ravecov.dyy();
  cov(1,2) = ravecov.dyz();
  cov(1,3) = ravecov.dypx();
  cov(1,4) = ravecov.dypy();
  cov(1,5) = ravecov.dypz();
 
  cov(2,0) = ravecov.dxz();
  cov(2,1) = ravecov.dyz();
  cov(2,2) = ravecov.dzz();
  cov(2,3) = ravecov.dzpx();
  cov(2,4) = ravecov.dzpy();
  cov(2,5) = ravecov.dzpz();
 
  cov(3,0) = ravecov.dxpx();
  cov(3,1) = ravecov.dypx();
  cov(3,2) = ravecov.dzpx();
  cov(3,3) = ravecov.dpxpx();
  cov(3,4) = ravecov.dpxpy();
  cov(3,5) = ravecov.dpxpz();
 
  cov(4,0) = ravecov.dxpy();
  cov(4,1) = ravecov.dypy();
  cov(4,2) = ravecov.dzpy();
  cov(4,3) = ravecov.dpxpy();
  cov(4,4) = ravecov.dpypy();
  cov(4,5) = ravecov.dpypz();
 
  cov(5,0) = ravecov.dxpz();
  cov(5,1) = ravecov.dypz();
  cov(5,2) = ravecov.dzpz();
  cov(5,3) = ravecov.dpxpz();
  cov(5,4) = ravecov.dpypz();
  cov(5,5) = ravecov.dpzpz();
 
  return cov;
}

◆ GFTracksToTracks()

std::vector< rave::Track > genfit::GFTracksToTracks	(	const std::vector< genfit::Track * > &	GFTracks,
		std::vector< genfit::MeasuredStateOnPlane * > *	GFStates,
		std::map< int, genfit::trackAndState > &	IdGFTrackStateMap,
		int	startID = `0`
	)

Convert a vector of genfit::Tracks to rave::Tracks also builds a map of unique ids to genfit::Tracks; These ids are stored in the rave::tracks. They are counted from startID. The map has to be passed to the GFRavePropagator, so that it can access the trackreps of the genfit::Tracks corresponding to the rave::tracks. If a vector of MeasuredStateOnPlane is provided, they will be used (instead of the tracks fitted states) to calculate the rave::Track parameters. Ownership over MeasuredStateOnPlane will be taken.

Definition at line 36 of file GFRaveConverters.cc.

                {
 
  unsigned int ntracks(GFTracks.size());
 
  if (GFStates != NULL)
    if (GFTracks.size() != GFStates->size()) {
      Exception exc("GFTracksToTracks ==> GFStates has not the same size as GFTracks!",__LINE__,__FILE__);
      throw exc;
    }
 
  std::vector < rave::Track > ravetracks;
  ravetracks.reserve(ntracks);
 
  for (unsigned int i=0; i<ntracks; ++i){
 
    if (GFTracks[i] == NULL) {
      Exception exc("GFTracksToTracks ==> genfit::Track is NULL",__LINE__,__FILE__);
      throw exc;
    }
 
    // only convert successfully fitted tracks!
    if (!GFTracks[i]->getFitStatus(GFTracks[i]->getCardinalRep())->isFitConverged()) continue;
 
    if (IdGFTrackStateMap.count(startID) > 0){
      Exception exc("GFTracksToTracks ==> IdGFTrackStateMap has already an entry for this id",__LINE__,__FILE__);
      throw exc;
    }
    IdGFTrackStateMap[startID].track_ = GFTracks[i];
    if (GFStates == NULL)
      IdGFTrackStateMap[startID].state_ = new MeasuredStateOnPlane(GFTracks[i]->getFittedState()); // here clones are made so that the state of the original GFTracks and their TrackReps will not be altered by the vertexing process
    else
      IdGFTrackStateMap[startID].state_ = (*GFStates)[i];
 
    ravetracks.push_back(GFTrackToTrack(IdGFTrackStateMap[startID], startID) );
 
    ++startID;
  }
 
  //std::cout << "IdGFTrackStateMap size " << IdGFTrackStateMap.size() << std::endl;
  return ravetracks;
}

◆ GFTrackToTrack()

rave::Track genfit::GFTrackToTrack	(	trackAndState	trackAndState,
		int	id = `-1`,
		std::string	tag = `""`
	)

Definition at line 83 of file GFRaveConverters.cc.

                                                                  {
 
  if (trackAndState.track_ == NULL) {
    Exception exc("GFTrackToTrack ==> originaltrack is NULL",__LINE__,__FILE__);
    throw exc;
  }
 
  if (! trackAndState.track_->getFitStatus()->isFitConverged()) {
    Exception exc("GFTrackToTrack ==> Trackfit is not converged",__LINE__,__FILE__);
    throw exc;
  }
 
  TVector3 pos, mom;
  TMatrixDSym cov;
 
  trackAndState.track_->getFittedState().getPosMomCov(pos, mom, cov);
 
  // state
  rave::Vector6D ravestate(pos.X(), pos.Y(), pos.Z(),
                           mom.X(), mom.Y(), mom.Z());
 
  // covariance
  rave::Covariance6D ravecov(cov(0,0), cov(1,0), cov(2,0),
                             cov(1,1), cov(2,1), cov(2,2),
                             cov(3,0), cov(4,0), cov(5,0),
                             cov(3,1), cov(4,1), cov(5,1),
                             cov(3,2), cov(4,2), cov(5,2),
                             cov(3,3), cov(4,3), cov(5,3),
                             cov(4,4), cov(5,4), cov(5,5));
 
  //std::cerr<<"create rave track with id " << id << std::endl;
  //std::cerr<<"  pos: "; Point3DToTVector3(ravestate.position()).Print();
  //std::cerr<<"  mom: "; Vector3DToTVector3(ravestate.momentum()).Print();
 
  rave::Track ret(id, ravestate, ravecov,
      trackAndState.track_->getFitStatus()->getCharge(),
      trackAndState.track_->getFitStatus()->getChi2(),
      trackAndState.track_->getFitStatus()->getNdf(),
      static_cast<void*>(const_cast<Track*>(trackAndState.track_)), tag);
 
  //std::cout << "ret.originalObject() " << ret.originalObject() << "\n";
  //std::cout << "ret.id() " << ret.id() << "\n";
 
  return ret;
}

◆ MeanExcEnergy_get() [1/2]

double genfit::MeanExcEnergy_get ( int Z )

Definition at line 223 of file TGeoMaterialInterface.cc.

                         {
  assert(Z >= 0 && Z < MeanExcEnergy_NELEMENTS);
  return MeanExcEnergy_vals[Z];
}

◆ MeanExcEnergy_get() [2/2]

double genfit::MeanExcEnergy_get ( TGeoMaterial * mat )

Definition at line 230 of file TGeoMaterialInterface.cc.

                                     {
  if (mat->IsMixture()) {
    double logMEE = 0.;
    double denom  = 0.;
    TGeoMixture* mix = (TGeoMixture*)mat;
    for (int i = 0; i < mix->GetNelements(); ++i) {
      int index = int(floor((mix->GetZmixt())[i]));
      logMEE += 1. / (mix->GetAmixt())[i] * (mix->GetWmixt())[i] * (mix->GetZmixt())[i] * log(MeanExcEnergy_get(index));
      denom  += (mix->GetWmixt())[i] * (mix->GetZmixt())[i] * 1. / (mix->GetAmixt())[i];
    }
    logMEE /= denom;
    return exp(logMEE);
  } else { // not a mixture
    int index = int(floor(mat->GetZ()));
    return MeanExcEnergy_get(index);
  }
}

◆ new_clone()

AbsFitterInfo * genfit::new_clone ( const AbsFitterInfo & a )

inline

Needed for boost cloneability:

Definition at line 107 of file AbsFitterInfo.h.

{
  return a.clone();
}

◆ operator!=() [1/2]

bool genfit::operator!=	(	const DetPlane &	lhs,
		const DetPlane &	rhs
	)

Definition at line 263 of file DetPlane.cc.

                                                          {
  return !(lhs==rhs);
}

◆ operator!=() [2/2]

bool genfit::operator!=	(	const MaterialProperties &	lhs,
		const MaterialProperties &	rhs
	)

Definition at line 39 of file MaterialProperties.cc.

                                                                               {
  return !(lhs==rhs);
}

◆ operator==() [1/4]

bool genfit::operator==	(	const DetPlane &	lhs,
		const DetPlane &	rhs
	)

Definition at line 241 of file DetPlane.cc.

                                                          {
  if (&lhs == &rhs)
    return true;
  static const double detplaneEpsilon = 1.E-5;
  if(
     fabs( (lhs.o_.X()-rhs.o_.X()) ) > detplaneEpsilon  ||
     fabs( (lhs.o_.Y()-rhs.o_.Y()) ) > detplaneEpsilon  ||
     fabs( (lhs.o_.Z()-rhs.o_.Z()) ) > detplaneEpsilon
     ) return false;
  else if(
    fabs( (lhs.u_.X()-rhs.u_.X()) ) > detplaneEpsilon  ||
    fabs( (lhs.u_.Y()-rhs.u_.Y()) ) > detplaneEpsilon  ||
    fabs( (lhs.u_.Z()-rhs.u_.Z()) ) > detplaneEpsilon
    ) return false;
  else if(
    fabs( (lhs.v_.X()-rhs.v_.X()) ) > detplaneEpsilon  ||
    fabs( (lhs.v_.Y()-rhs.v_.Y()) ) > detplaneEpsilon  ||
    fabs( (lhs.v_.Z()-rhs.v_.Z()) ) > detplaneEpsilon
    ) return false;
  return true;
}

◆ operator==() [2/4]

bool genfit::operator==	(	const MaterialProperties &	lhs,
		const MaterialProperties &	rhs
	)

Definition at line 26 of file MaterialProperties.cc.

                                                                              {
  if (&lhs == &rhs)
    return true;
  if (lhs.density_ != rhs.density_ ||
      lhs.Z_ != rhs.Z_ ||
      lhs.A_ != rhs.A_ ||
      lhs.radiationLength_ != rhs.radiationLength_ ||
      lhs.mEE_ != rhs.mEE_)
    return false;
 
  return true;
}

◆ operator==() [3/4]

bool genfit::operator==	(	const TrackCand &	lhs,
		const TrackCand &	rhs
	)

Definition at line 188 of file TrackCand.cc.

                                                            {
  if(lhs.getNHits() != rhs.getNHits()) return false;
  for (unsigned int i = 0; i < lhs.getNHits(); ++i){
    if (lhs.getHit(i) != rhs.getHit(i)) return false;
  }
  return true;
}

◆ operator==() [4/4]

bool genfit::operator==	(	const TrackCandHit &	lhs,
		const TrackCandHit &	rhs
	)

Definition at line 47 of file TrackCandHit.cc.

                                                                  {
  if(lhs.detId_ == rhs.detId_ &&
     lhs.hitId_ == rhs.hitId_ &&
     lhs.planeId_ == rhs.planeId_)
    return true;
  return false;
}

◆ PlaneToGFDetPlane()

SharedPlanePtr genfit::PlaneToGFDetPlane ( const ravesurf::Plane & rplane )

Definition at line 225 of file GFRaveConverters.cc.

                                               {
  return SharedPlanePtr(new DetPlane(Point3DToTVector3(rplane.position()),
                    Vector3DToTVector3(rplane.normalVector()) ));
}

◆ Point3DToTVector3()

TVector3 genfit::Point3DToTVector3 ( const rave::Point3D & v )

Definition at line 232 of file GFRaveConverters.cc.

                                        {
  return TVector3(v.x(), v.y(), v.z());
}

◆ RaveToGFVertex()

GFRaveVertex * genfit::RaveToGFVertex	(	const rave::Vertex &	raveVertex,
		const std::map< int, genfit::trackAndState > &	IdGFTrackStateMap
	)

Definition at line 141 of file GFRaveConverters.cc.

                                                              {
 
  if (!(raveVertex.isValid())) {
    Exception exc("RaveToGFVertex ==> rave Vertex is not valid!",__LINE__,__FILE__);
    throw exc;
  }
 
  std::vector < std::pair < float, rave::Track > > raveWeightedTracks(raveVertex.weightedTracks());
  std::vector < std::pair < float, rave::Track > > raveSmoothedTracks(raveVertex.weightedRefittedTracks());
 
  int id;
  unsigned int nTrks(raveWeightedTracks.size());
 
  // check if rave vertex has  refitted tracks
  bool smoothing(true);
  if (! (raveVertex.hasRefittedTracks()) ) {
    smoothing = false;
  }
 
  // check numbers of tracks and smoothed tracks
  if (smoothing && nTrks != raveSmoothedTracks.size()){
    Exception exc("RaveToGFVertex ==> number of smoothed tracks != number of tracks",__LINE__,__FILE__);
    throw exc;
  }
 
  // (smoothed) track parameters
  std::vector < GFRaveTrackParameters* > trackParameters;
  trackParameters.reserve(nTrks);
 
  // convert tracks
  for (unsigned int i=0; i<nTrks; ++i){
    id = raveWeightedTracks[i].second.id();
 
    if (IdGFTrackStateMap.count(id) == 0){
      Exception exc("RaveToGFVertex ==> rave track id is not present in IdGFTrackStateMap",__LINE__,__FILE__);
      throw exc;
    }
 
    GFRaveTrackParameters* trackparams;
 
    if(smoothing) {
      // convert smoothed track parameters
      trackparams = new GFRaveTrackParameters(IdGFTrackStateMap.at(id).track_, //track
          IdGFTrackStateMap.at(id).state_, //state
          raveWeightedTracks[i].first, //weight
          Vector6DToTVectorD(raveSmoothedTracks[i].second.state()), //smoothed state
          Covariance6DToTMatrixDSym(raveSmoothedTracks[i].second.error()), //smoothed cov
          true);
    }
    else {
      // convert track parameters, no smoothed tracks available
      trackparams = new GFRaveTrackParameters(IdGFTrackStateMap.at(id).track_, //track
          IdGFTrackStateMap.at(id).state_, //state
          raveWeightedTracks[i].first, //weight
          Vector6DToTVectorD(raveWeightedTracks[i].second.state()), //state
          Covariance6DToTMatrixDSym(raveWeightedTracks[i].second.error()), //cov
          false);
    }
    trackParameters.push_back(trackparams);
  }
 
  return new GFRaveVertex(Point3DToTVector3(raveVertex.position()),
                          Covariance3DToTMatrixDSym(raveVertex.error()),
                          trackParameters,
                          raveVertex.ndf(), raveVertex.chiSquared(), raveVertex.id());
}

◆ RaveToGFVertices()

void genfit::RaveToGFVertices	(	std::vector< GFRaveVertex * > *	GFVertices,
		const std::vector< rave::Vertex > &	raveVertices,
		const std::map< int, genfit::trackAndState > &	IdGFTrackStateMap
	)

Definition at line 210 of file GFRaveConverters.cc.

                                                              {
 
  unsigned int nVert(raveVertices.size());
 
  GFVertices->reserve(nVert);
 
  for (unsigned int i=0; i<nVert; ++i){
    GFVertices->push_back(RaveToGFVertex(raveVertices[i], IdGFTrackStateMap));
  }
}

◆ setData()

void genfit::setData	(	const rave::Track &	orig,
		MeasuredStateOnPlane *	state
	)

set state and cov of a MeasuredStateOnPlane according to rave track

Definition at line 131 of file GFRaveConverters.cc.

                                                           {
 
  state->setPosMomCov(TVector3(orig.state().x(), orig.state().y(), orig.state().z()),
                      TVector3(orig.state().px(), orig.state().py(), orig.state().pz()),
                      Covariance6DToTMatrixDSym(orig.error()));
 
}

◆ TMatrixDSymToCovariance3D()

rave::Covariance3D genfit::TMatrixDSymToCovariance3D ( const TMatrixDSym & matrix )

Definition at line 335 of file GFRaveConverters.cc.

                                                     {
  if (matrix.GetNrows()!=3) {
    Exception exc("TMatrixDSymToCovariance3D ==> TMatrixDSym is not 3x3!",__LINE__,__FILE__);
    throw exc;
  }
 
  return rave::Covariance3D(matrix(0,0), matrix(0,1), matrix(0,2),
                            matrix(1,1), matrix(1,2), matrix(2,2));
 
}

◆ TVector3ToPoint3D()

rave::Point3D genfit::TVector3ToPoint3D ( const TVector3 & vec )

Definition at line 329 of file GFRaveConverters.cc.

                                       {
  return rave::Point3D(vec.X(), vec.Y(), vec.Z());
}

◆ Vector3DToTVector3()

TVector3 genfit::Vector3DToTVector3 ( const rave::Vector3D & v )

Definition at line 237 of file GFRaveConverters.cc.

                                          {
  return TVector3(v.x(), v.y(), v.z());
}

◆ Vector6DToTVectorD()

TVectorD genfit::Vector6DToTVectorD ( const rave::Vector6D & ravevec )

Definition at line 263 of file GFRaveConverters.cc.

                                               {
  TVectorD vec(6);
 
  vec[0] = ravevec.x();
  vec[1] = ravevec.y();
  vec[2] = ravevec.z();
 
  vec[3] = ravevec.px();
  vec[4] = ravevec.py();
  vec[5] = ravevec.pz();
 
  return vec;
}

Variable Documentation

◆ MeanExcEnergy_NELEMENTS

const int genfit::MeanExcEnergy_NELEMENTS = 93

Definition at line 218 of file TGeoMaterialInterface.cc.

◆ MeanExcEnergy_vals

const double genfit::MeanExcEnergy_vals[] = {1.e15, 19.2, 41.8, 40.0, 63.7, 76.0, 78., 82.0, 95.0, 115.0, 137.0, 149.0, 156.0, 166.0, 173.0, 173.0, 180.0, 174.0, 188.0, 190.0, 191.0, 216.0, 233.0, 245.0, 257.0, 272.0, 286.0, 297.0, 311.0, 322.0, 330.0, 334.0, 350.0, 347.0, 348.0, 343.0, 352.0, 363.0, 366.0, 379.0, 393.0, 417.0, 424.0, 428.0, 441.0, 449.0, 470.0, 470.0, 469.0, 488.0, 488.0, 487.0, 485.0, 491.0, 482.0, 488.0, 491.0, 501.0, 523.0, 535.0, 546.0, 560.0, 574.0, 580.0, 591.0, 614.0, 628.0, 650.0, 658.0, 674.0, 684.0, 694.0, 705.0, 718.0, 727.0, 736.0, 746.0, 757.0, 790.0, 790.0, 800.0, 810.0, 823.0, 823.0, 830.0, 825.0, 794.0, 827.0, 826.0, 841.0, 847.0, 878.0, 890.0}

Definition at line 219 of file TGeoMaterialInterface.cc.

219{1.e15, 19.2, 41.8, 40.0, 63.7, 76.0, 78., 82.0, 95.0, 115.0, 137.0, 149.0, 156.0, 166.0, 173.0, 173.0, 180.0, 174.0, 188.0, 190.0, 191.0, 216.0, 233.0, 245.0, 257.0, 272.0, 286.0, 297.0, 311.0, 322.0, 330.0, 334.0, 350.0, 347.0, 348.0, 343.0, 352.0, 363.0, 366.0, 379.0, 393.0, 417.0, 424.0, 428.0, 441.0, 449.0, 470.0, 470.0, 469.0, 488.0, 488.0, 487.0, 485.0, 491.0, 482.0, 488.0, 491.0, 501.0, 523.0, 535.0, 546.0, 560.0, 574.0, 580.0, 591.0, 614.0, 628.0, 650.0, 658.0, 674.0, 684.0, 694.0, 705.0, 718.0, 727.0, 736.0, 746.0, 757.0, 790.0, 790.0, 800.0, 810.0, 823.0, 823.0, 830.0, 825.0, 794.0, 827.0, 826.0, 841.0, 847.0, 878.0, 890.0};

Namespaces

Classes

Typedefs

Enumerations

Functions

Variables