org.hps.recon.tracking

Class TrackUtils

java.lang.Object

org.hps.recon.tracking.TrackUtils

public class TrackUtils
extends Object

Assorted helper functions for the track and helix objects in lcsim. Re-use as much of HelixUtils as possible.

Author:

Omar Moreno

Nested Class Summary

Nested Classes

Modifier and Type	Class and Description
private static class	TrackUtils.LayerComparator

Field Summary

Fields

Modifier and Type	Field and Description
private static org.apache.commons.math3.util.Pair<org.lcsim.event.EventHeader,org.lcsim.event.RelationalTable>	hitToRotatedCache
private static org.apache.commons.math3.util.Pair<org.lcsim.event.EventHeader,org.lcsim.event.RelationalTable>	hitToStripsCache

Constructor Summary

Constructors

Modifier	Constructor and Description
private	TrackUtils() Private constructor to make class static only

Method Summary

Modifier and Type	Method and Description
static double	ArcTan(double x, double y)
static double	calculateCurvature(double p, double q, double B)
static double	calculatePhi(double px, double py)
static double	calculatePhi(double x, double y, double xc, double yc, double sign)
static double	calculateTanLambda(double pz, double p)
static double	Cos(double arg)
static void	cut(int[] cuts, EventQuality.Cut bit)
static Hep3Vector	extrapolateHelixToXPlane(org.lcsim.fit.helicaltrack.HelicalTrackFit htf, double x) Extrapolate helix to a position along the x-axis.
static Hep3Vector	extrapolateHelixToXPlane(org.lcsim.event.Track track, double x) Extrapolate track to a position along the x-axis.
static Hep3Vector	extrapolateHelixToXPlane(org.lcsim.event.TrackState track, double x)
static Hep3Vector	extrapolateTrack(org.lcsim.fit.helicaltrack.HelicalTrackFit helix, double z) Deprecated.
static Hep3Vector	extrapolateTrack(org.lcsim.event.Track track, double z) Deprecated.
static Hep3Vector	extrapolateTrack(org.lcsim.event.Track track, double z, org.lcsim.geometry.Detector detector) Extrapolate track to given position, using dipole position from geometry.
static Hep3Vector	extrapolateTrack(org.lcsim.event.TrackState track, double z) Deprecated.
static Hep3Vector	extrapolateTrackPositionToSensor(org.lcsim.event.Track track, HpsSiSensor sensor, List<HpsSiSensor> sensors, double bfield)
static org.lcsim.event.base.BaseTrackState	extrapolateTrackUsingFieldMap(org.lcsim.event.Track track, double startPositionX, double endPositionX, double stepSize, org.lcsim.geometry.FieldMap fieldMap) Backward compatibility function for extrapolateTrackUsingFieldMap .
static org.lcsim.event.base.BaseTrackState	extrapolateTrackUsingFieldMap(org.lcsim.event.TrackState track, double startPositionX, double endPosition, double stepSize, double epsilon, org.lcsim.geometry.FieldMap fieldMap)
static org.lcsim.event.base.BaseTrackState	extrapolateTrackUsingFieldMap(org.lcsim.event.TrackState track, double startPositionX, double endPosition, double stepSize, org.lcsim.geometry.FieldMap fieldMap) Iteratively extrapolates a track to a specified value of x (z in detector frame) using the full 3D field map.
static Pair<Hep3Vector,Hep3Vector>	extrapolateTrackUsingFieldMapRK(org.lcsim.event.TrackState ts, Hep3Vector startPosition, double distanceZ, double stepSize, org.lcsim.geometry.FieldMap fM)
static double	findTriangleArea(double x0, double y0, double x1, double y1, double x2, double y2)
static Hep3Vector	getBField(org.lcsim.geometry.Detector detector)
static int	getCharge(org.lcsim.event.Track track)
static double	getCosTheta(org.lcsim.event.Track track)
static double	getCosTheta(org.lcsim.event.TrackState track)
static SymmetricMatrix	getCovarianceAtNewRefPoint(double[] newRefPoint, double[] __refPoint, double[] par, SymmetricMatrix helixcov)
static double	getDoca(org.lcsim.event.Track track)
static double	getDoca(org.lcsim.event.TrackState track)
static Hep3Vector	getHelixPlaneIntercept(org.lcsim.fit.helicaltrack.HelicalTrackFit helfit, HelicalTrackStripGbl strip, double bfield) Calculate the point of interception between the helix and a plane in space.
static Hep3Vector	getHelixPlaneIntercept(org.lcsim.fit.helicaltrack.HelicalTrackFit helfit, Hep3Vector unit_vec_normal_to_plane, Hep3Vector point_on_plane, double bfield) Calculate the point of interception between the helix and a plane in space.
static Hep3Vector	getHelixPlaneIntercept(org.lcsim.fit.helicaltrack.HelicalTrackFit helfit, Hep3Vector unit_vec_normal_to_plane, Hep3Vector point_on_plane, double bfield, double initial_s)
static Hep3Vector	getHelixPlanePositionIter(org.lcsim.fit.helicaltrack.HelicalTrackFit helix, Hep3Vector origin, Hep3Vector normal, double eps) Deprecated.
private static Hep3Vector	getHelixPosAtX(org.lcsim.fit.helicaltrack.HelicalTrackFit helix, double x) Deprecated.
static int[]	getHitsInTopBottom(org.lcsim.event.Track track)
static org.lcsim.event.RelationalTable	getHitToRotatedTable(org.lcsim.event.EventHeader event)
static org.lcsim.event.RelationalTable	getHitToStripsTable(org.lcsim.event.EventHeader event)
static org.lcsim.fit.helicaltrack.HelicalTrackFit	getHTF(double[] par)
static org.lcsim.fit.helicaltrack.HelicalTrackFit	getHTF(org.lcsim.event.MCParticle mcp, double Bz) Transform MCParticle into a Helix object.
static org.lcsim.fit.helicaltrack.HelicalTrackFit	getHTF(org.lcsim.event.MCParticle mcp, Hep3Vector origin, double Bz) Transform MCParticle into a HelicalTrackFit object.
static org.lcsim.fit.helicaltrack.HelicalTrackFit	getHTF(org.lcsim.event.Track track)
static org.lcsim.fit.helicaltrack.HelicalTrackFit	getHTF(org.lcsim.event.TrackState state)
static double	getIsolation(org.lcsim.event.TrackerHit strip, org.lcsim.event.TrackerHit otherStrip, org.lcsim.event.RelationalTable hitToStrips, org.lcsim.event.RelationalTable hitToRotated) Compute strip isolation, defined as the minimum distance to another strip in the same sensor.
static Double[]	getIsolations(org.lcsim.event.Track trk, org.lcsim.event.RelationalTable hitToStrips, org.lcsim.event.RelationalTable hitToRotated)
static Double[]	getIsolations(org.lcsim.event.Track trk, org.lcsim.event.RelationalTable hitToStrips, org.lcsim.event.RelationalTable hitToRotated, int layers) Make an array with isolations for all 12 strip layers.
static int	getLayer(org.lcsim.event.TrackerHit strip)
static org.lcsim.event.MCParticle	getMatchedTruthParticle(org.lcsim.event.Track track)
static Hep3Vector	getMomentum(double omega, double phi0, double tanL, double magneticField)
static double[]	getParametersAtNewRefPoint(double[] newRefPoint, double[] __refPoint, double[] parameters) Change reference point of helix (following L3 Internal Note 1666.)
static double[]	getParametersAtNewRefPoint(double[] newRefPoint, HpsHelicalTrackFit helicalTrackFit) Change reference point of helix (following L3 Internal Note 1666.)
static double[]	getParametersAtNewRefPoint(double[] newRefPoint, org.lcsim.event.TrackState trkState) Change reference point of helix (following L3 Internal Note 1666.)
static double[]	getParametersFromPointAndMomentum(Hep3Vector point, Hep3Vector momentum, int charge, double BField)
static double[]	getParametersFromPointAndMomentum(Hep3Vector point, Hep3Vector momentum, int charge, double BField, boolean writeIt) Deprecated.
static double	getPhi(org.lcsim.event.Track track, Hep3Vector position)
static double	getPhi(org.lcsim.event.TrackState track, Hep3Vector position)
static double	getPhi0(org.lcsim.event.Track track)
static double	getPhi0(org.lcsim.event.TrackState track)
static Hep3Vector	getPoint(double doca, double phi0, double z0)
static double	getR(org.lcsim.event.Track track)
static double	getR(org.lcsim.event.TrackState track)
static double	getSinTheta(org.lcsim.event.Track track)
static double	getSinTheta(org.lcsim.event.TrackState track)
static List<org.lcsim.event.TrackerHit>	getStripHits(org.lcsim.event.Track track, org.lcsim.event.RelationalTable hitToStrips, org.lcsim.event.RelationalTable hitToRotated)
static double	getTanLambda(org.lcsim.event.Track track)
static double	getTanLambda(org.lcsim.event.TrackState track)
static org.lcsim.event.base.BaseTrackState	getTrackExtrapAtEcal(org.lcsim.event.Track track, org.lcsim.geometry.FieldMap fieldMap) Get position of a track extrapolated to the ECAL face in the HPS test run 2012
static org.lcsim.event.base.BaseTrackState	getTrackExtrapAtEcal(org.lcsim.event.TrackState track, org.lcsim.geometry.FieldMap fieldMap)
static org.lcsim.event.base.BaseTrackState	getTrackExtrapAtEcalRK(org.lcsim.event.Track trk, org.lcsim.geometry.FieldMap fM)
static org.lcsim.event.base.BaseTrackState	getTrackExtrapAtEcalRK(org.lcsim.event.Track trk, org.lcsim.geometry.FieldMap fM, double stepSize)
static org.lcsim.event.base.BaseTrackState	getTrackExtrapAtEcalRK(org.lcsim.event.TrackState ts, org.lcsim.geometry.FieldMap fM)
static org.lcsim.event.base.BaseTrackState	getTrackExtrapAtEcalRK(org.lcsim.event.TrackState ts, org.lcsim.geometry.FieldMap fM, double stepSize)
static Hep3Vector	getTrackPositionAtEcal(org.lcsim.event.Track track) Deprecated.
static Hep3Vector	getTrackPositionAtEcal(org.lcsim.event.TrackState track) Deprecated.
static Hep3Vector	getTrackPositionAtHarp(org.lcsim.event.Track track) Deprecated.
static org.lcsim.event.TrackState	getTrackStateAtECal(org.lcsim.event.Track trk)
static org.lcsim.event.TrackState	getTrackStateAtLocation(org.lcsim.event.Track trk, int location) Port of Track.getTrackState(int location) from the C++ LCIO API.
static double	getTrackTime(org.lcsim.event.Track track, org.lcsim.event.RelationalTable hitToStrips, org.lcsim.event.RelationalTable hitToRotated)
static double	getTrackTimeSD(org.lcsim.event.Track track, org.lcsim.event.RelationalTable hitToStrips, org.lcsim.event.RelationalTable hitToRotated)
static org.lcsim.util.swim.Trajectory	getTrajectory(Hep3Vector p0, org.lcsim.spacegeom.SpacePoint r0, double q, double B)
static double	getX0(double doca, double phi0)
static double	getX0(org.lcsim.event.Track track)
static double	getX0(org.lcsim.event.TrackState track)
static double	getY0(double doca, double phi0)
static double	getY0(org.lcsim.event.Track track)
static double	getY0(org.lcsim.event.TrackState track)
static double	getZ0(org.lcsim.event.Track track)
static double	getZ0(org.lcsim.event.TrackState track)
static Hep3Vector	globalToSensor(Hep3Vector trkpos, HpsSiSensor sensor)
static boolean	hasSharedHits(org.lcsim.event.Track track, List<org.lcsim.event.Track> tracklist)
static boolean	hasSharedStrips(org.lcsim.event.Track track1, org.lcsim.event.Track track2, org.lcsim.event.RelationalTable hitToStrips, org.lcsim.event.RelationalTable hitToRotated) Tells if there are shared 2D hits between tracks
static boolean	hasTopBotHit(org.lcsim.event.Track track)
static boolean	isBottomTrack(org.lcsim.fit.helicaltrack.HelicalTrackFit htf)
static boolean	isBottomTrack(org.lcsim.event.Track track, int minhits)
static boolean	isGoodTrack(org.lcsim.event.Track track, List<org.lcsim.event.Track> tracklist, EventQuality.Quality trk_quality)
static boolean	isSharedHit(org.lcsim.event.TrackerHit hit, List<org.lcsim.event.Track> othertracks)
static boolean	isSharedHit(org.lcsim.event.TrackerHit hit, org.lcsim.event.Track track)
static int	isTopOrBottomTrack(org.lcsim.event.Track track, int minhits)
static boolean	isTopTrack(org.lcsim.fit.helicaltrack.HelicalTrackFit htf)
static boolean	isTopTrack(org.lcsim.event.Track track, int minhits)
static org.lcsim.fit.helicaltrack.HelicalTrackHit	makeHelicalTrackHitFromTrackerHit(org.lcsim.event.TrackerHit hit)
static org.lcsim.recon.tracking.seedtracker.SeedTrack	makeSeedTrackFromBaseTrack(org.lcsim.event.Track track)
static org.lcsim.event.Track	mostSharedHitTrack(org.lcsim.event.Track track, List<org.lcsim.event.Track> tracklist) Number of shared 3D hits between tracks
static int	numberOfSharedHits(org.lcsim.event.Track track, List<org.lcsim.event.Track> tracklist) Number of shared 3D hits between tracks
static int	numberOfSharedHits(org.lcsim.event.Track track1, org.lcsim.event.Track track2) Number of shared 3D hits between two tracks
static int	numberOfSharedStrips(org.lcsim.event.Track track1, org.lcsim.event.Track track2, org.lcsim.event.RelationalTable hitToStrips, org.lcsim.event.RelationalTable hitToRotated) Number of 2D hits shared between tracks
static int	passTrackSelections(org.lcsim.event.Track track, List<org.lcsim.event.Track> tracklist, EventQuality.Quality trk_quality)
static double	Power(double arg, double pow)
static double	Sec(double arg)
static boolean	sensorContainsTrack(Hep3Vector trackPosition, org.lcsim.detector.tracker.silicon.SiSensor sensor)
static double	Sin(double arg)
static List<org.lcsim.event.TrackerHit>	sortHits(Collection<org.lcsim.event.TrackerHit> hits)
static double	Tan(double arg)
static void	writeTrajectory(Hep3Vector momentum, Hep3Vector position, int q, double bField, String outFile)

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Field Detail

hitToStripsCache

private static org.apache.commons.math3.util.Pair<org.lcsim.event.EventHeader,org.lcsim.event.RelationalTable> hitToStripsCache

hitToRotatedCache

private static org.apache.commons.math3.util.Pair<org.lcsim.event.EventHeader,org.lcsim.event.RelationalTable> hitToRotatedCache

Constructor Detail

TrackUtils

private TrackUtils()

Private constructor to make class static only

Method Detail

extrapolateTrackPositionToSensor

public static Hep3Vector extrapolateTrackPositionToSensor(org.lcsim.event.Track track,
                                                          HpsSiSensor sensor,
                                                          List<HpsSiSensor> sensors,
                                                          double bfield)

extrapolateHelixToXPlane

public static Hep3Vector extrapolateHelixToXPlane(org.lcsim.event.Track track,
                                                  double x)

Extrapolate track to a position along the x-axis. Turn the track into a helix object in order to use HelixUtils.

Parameters:

track -

x -

Returns:

the position along the x-axis

extrapolateHelixToXPlane

public static Hep3Vector extrapolateHelixToXPlane(org.lcsim.event.TrackState track,
                                                  double x)

getParametersFromPointAndMomentum

public static double[] getParametersFromPointAndMomentum(Hep3Vector point,
                                                         Hep3Vector momentum,
                                                         int charge,
                                                         double BField)

getParametersFromPointAndMomentum

@Deprecated
public static double[] getParametersFromPointAndMomentum(Hep3Vector point,
                                                                     Hep3Vector momentum,
                                                                     int charge,
                                                                     double BField,
                                                                     boolean writeIt)
                                                              throws FileNotFoundException

Deprecated.

Throws:

FileNotFoundException

getParametersAtNewRefPoint

public static double[] getParametersAtNewRefPoint(double[] newRefPoint,
                                                  org.lcsim.event.TrackState trkState)

Change reference point of helix (following L3 Internal Note 1666.)

Parameters:

newRefPoint - - The new reference point in XY

getParametersAtNewRefPoint

public static double[] getParametersAtNewRefPoint(double[] newRefPoint,
                                                  HpsHelicalTrackFit helicalTrackFit)

Change reference point of helix (following L3 Internal Note 1666.)

Parameters:

newRefPoint - - The new reference point in XY

getParametersAtNewRefPoint

public static double[] getParametersAtNewRefPoint(double[] newRefPoint,
                                                  double[] __refPoint,
                                                  double[] parameters)

Change reference point of helix (following L3 Internal Note 1666.)

Parameters:

newRefPoint - - The new reference point in XY

getCovarianceAtNewRefPoint

public static SymmetricMatrix getCovarianceAtNewRefPoint(double[] newRefPoint,
                                                         double[] __refPoint,
                                                         double[] par,
                                                         SymmetricMatrix helixcov)

extrapolateHelixToXPlane

public static Hep3Vector extrapolateHelixToXPlane(org.lcsim.fit.helicaltrack.HelicalTrackFit htf,
                                                  double x)

Extrapolate helix to a position along the x-axis. Re-use HelixUtils.

Parameters:

htf -

x -

Returns:

the position along the x-axis

getPhi

public static double getPhi(org.lcsim.event.Track track,
                            Hep3Vector position)

getX0

public static double getX0(org.lcsim.event.Track track)

getR

public static double getR(org.lcsim.event.Track track)

getY0

public static double getY0(org.lcsim.event.Track track)

getDoca

public static double getDoca(org.lcsim.event.Track track)

getPhi0

public static double getPhi0(org.lcsim.event.Track track)

getZ0

public static double getZ0(org.lcsim.event.Track track)

getTanLambda

public static double getTanLambda(org.lcsim.event.Track track)

getSinTheta

public static double getSinTheta(org.lcsim.event.Track track)

getCosTheta

public static double getCosTheta(org.lcsim.event.Track track)

getPhi

public static double getPhi(org.lcsim.event.TrackState track,
                            Hep3Vector position)

getX0

public static double getX0(org.lcsim.event.TrackState track)

getX0

public static double getX0(double doca,
                           double phi0)

getR

public static double getR(org.lcsim.event.TrackState track)

getY0

public static double getY0(org.lcsim.event.TrackState track)

getY0

public static double getY0(double doca,
                           double phi0)

getDoca

public static double getDoca(org.lcsim.event.TrackState track)

getPhi0

public static double getPhi0(org.lcsim.event.TrackState track)

getZ0

public static double getZ0(org.lcsim.event.TrackState track)

getPoint

public static Hep3Vector getPoint(double doca,
                                  double phi0,
                                  double z0)

getTanLambda

public static double getTanLambda(org.lcsim.event.TrackState track)

getSinTheta

public static double getSinTheta(org.lcsim.event.TrackState track)

getCosTheta

public static double getCosTheta(org.lcsim.event.TrackState track)

getCharge

public static int getCharge(org.lcsim.event.Track track)

getHelixPlaneIntercept

public static Hep3Vector getHelixPlaneIntercept(org.lcsim.fit.helicaltrack.HelicalTrackFit helfit,
                                                Hep3Vector unit_vec_normal_to_plane,
                                                Hep3Vector point_on_plane,
                                                double bfield)

Calculate the point of interception between the helix and a plane in space. Uses an iterative procedure. This function makes assumptions on the sign and convecntion of the B-field. Be careful.

Parameters:

helfit - - helix

unit_vec_normal_to_plane - - unit vector normal to the plane

point_on_plane - - point on the plane

bfield - - magnetic field value

Returns:

point at intercept

getHelixPlaneIntercept

public static Hep3Vector getHelixPlaneIntercept(org.lcsim.fit.helicaltrack.HelicalTrackFit helfit,
                                                Hep3Vector unit_vec_normal_to_plane,
                                                Hep3Vector point_on_plane,
                                                double bfield,
                                                double initial_s)

getHelixPlaneIntercept

public static Hep3Vector getHelixPlaneIntercept(org.lcsim.fit.helicaltrack.HelicalTrackFit helfit,
                                                HelicalTrackStripGbl strip,
                                                double bfield)

Calculate the point of interception between the helix and a plane in space. Uses an iterative procedure.

Parameters:

helfit - - helix

strip - - strip cluster that will define the plane

bfield - - magnetic field value

Returns:

point at intercept

getTrackPositionAtHarp

@Deprecated
public static Hep3Vector getTrackPositionAtHarp(org.lcsim.event.Track track)

Deprecated.

Get position of a track extrapolated to the HARP in the HPS test run 2012

Parameters:

track -

Returns:

position at HARP

getTrackExtrapAtEcal

public static org.lcsim.event.base.BaseTrackState getTrackExtrapAtEcal(org.lcsim.event.Track track,
                                                                       org.lcsim.geometry.FieldMap fieldMap)

Get position of a track extrapolated to the ECAL face in the HPS test run 2012

Parameters:

track -

Returns:

position at ECAL

getTrackExtrapAtEcal

public static org.lcsim.event.base.BaseTrackState getTrackExtrapAtEcal(org.lcsim.event.TrackState track,
                                                                       org.lcsim.geometry.FieldMap fieldMap)

getTrackExtrapAtEcalRK

public static org.lcsim.event.base.BaseTrackState getTrackExtrapAtEcalRK(org.lcsim.event.TrackState ts,
                                                                         org.lcsim.geometry.FieldMap fM)

getTrackExtrapAtEcalRK

public static org.lcsim.event.base.BaseTrackState getTrackExtrapAtEcalRK(org.lcsim.event.TrackState ts,
                                                                         org.lcsim.geometry.FieldMap fM,
                                                                         double stepSize)

getTrackExtrapAtEcalRK

public static org.lcsim.event.base.BaseTrackState getTrackExtrapAtEcalRK(org.lcsim.event.Track trk,
                                                                         org.lcsim.geometry.FieldMap fM,
                                                                         double stepSize)

getTrackExtrapAtEcalRK

public static org.lcsim.event.base.BaseTrackState getTrackExtrapAtEcalRK(org.lcsim.event.Track trk,
                                                                         org.lcsim.geometry.FieldMap fM)

extrapolateTrack

@Deprecated
public static Hep3Vector extrapolateTrack(org.lcsim.event.Track track,
                                                      double z)

Deprecated.

Extrapolate track to given position. For backwards compatibility.

Parameters:

track - - to be extrapolated

z -

Returns:

extrapolated position

getTrackPositionAtEcal

@Deprecated
public static Hep3Vector getTrackPositionAtEcal(org.lcsim.event.Track track)

Deprecated.

getTrackPositionAtEcal

@Deprecated
public static Hep3Vector getTrackPositionAtEcal(org.lcsim.event.TrackState track)

Deprecated.

extrapolateTrack

@Deprecated
public static Hep3Vector extrapolateTrack(org.lcsim.event.TrackState track,
                                                      double z)

Deprecated.

Extrapolate track to given position.

Parameters:

track - - to be extrapolated

z -

Returns:

extrapolated position

extrapolateTrack

public static Hep3Vector extrapolateTrack(org.lcsim.event.Track track,
                                          double z,
                                          org.lcsim.geometry.Detector detector)

Extrapolate track to given position, using dipole position from geometry.

Parameters:

track - - position along the x-axis of the helix in lcsim coordinates

Returns:

extrapolated position

extrapolateTrack

@Deprecated
public static Hep3Vector extrapolateTrack(org.lcsim.fit.helicaltrack.HelicalTrackFit helix,
                                                      double z)

Deprecated.

Extrapolate helix to given position

Parameters:

helix - - to be extrapolated

z - - position along the x-axis of the helix in lcsim coordiantes

Returns:

the extrapolated position

getHelixPlanePositionIter

@Deprecated
public static Hep3Vector getHelixPlanePositionIter(org.lcsim.fit.helicaltrack.HelicalTrackFit helix,
                                                               Hep3Vector origin,
                                                               Hep3Vector normal,
                                                               double eps)

Deprecated.

Parameters:

helix - input helix object

origin - of the plane to intercept

normal - of the plane to intercept

eps - criteria on the distance to the plane before stopping iteration

Returns:

position in space at the intercept of the plane

getHelixPosAtX

@Deprecated
private static Hep3Vector getHelixPosAtX(org.lcsim.fit.helicaltrack.HelicalTrackFit helix,
                                                     double x)

Deprecated.

findTriangleArea

public static double findTriangleArea(double x0,
                                      double y0,
                                      double x1,
                                      double y1,
                                      double x2,
                                      double y2)

sensorContainsTrack

public static boolean sensorContainsTrack(Hep3Vector trackPosition,
                                          org.lcsim.detector.tracker.silicon.SiSensor sensor)

getHitsInTopBottom

public static int[] getHitsInTopBottom(org.lcsim.event.Track track)

isTopTrack

public static boolean isTopTrack(org.lcsim.event.Track track,
                                 int minhits)

isBottomTrack

public static boolean isBottomTrack(org.lcsim.event.Track track,
                                    int minhits)

isTopOrBottomTrack

public static int isTopOrBottomTrack(org.lcsim.event.Track track,
                                     int minhits)

hasTopBotHit

public static boolean hasTopBotHit(org.lcsim.event.Track track)

isSharedHit

public static boolean isSharedHit(org.lcsim.event.TrackerHit hit,
                                  List<org.lcsim.event.Track> othertracks)

isSharedHit

public static boolean isSharedHit(org.lcsim.event.TrackerHit hit,
                                  org.lcsim.event.Track track)

numberOfSharedHits

public static int numberOfSharedHits(org.lcsim.event.Track track,
                                     List<org.lcsim.event.Track> tracklist)

Number of shared 3D hits between tracks

Parameters:

track -

tracklist -

Returns:

number of 3D hits shared on a track

numberOfSharedHits

public static int numberOfSharedHits(org.lcsim.event.Track track1,
                                     org.lcsim.event.Track track2)

Number of shared 3D hits between two tracks

Parameters:

track -

tracklist -

Returns:

number of 3D hits shared between two tracks

mostSharedHitTrack

public static org.lcsim.event.Track mostSharedHitTrack(org.lcsim.event.Track track,
                                                       List<org.lcsim.event.Track> tracklist)

Number of shared 3D hits between tracks

Parameters:

track -

tracklist -

Returns:

the track associated with the most shared hits

hasSharedHits

public static boolean hasSharedHits(org.lcsim.event.Track track,
                                    List<org.lcsim.event.Track> tracklist)

cut

public static void cut(int[] cuts,
                       EventQuality.Cut bit)

isGoodTrack

public static boolean isGoodTrack(org.lcsim.event.Track track,
                                  List<org.lcsim.event.Track> tracklist,
                                  EventQuality.Quality trk_quality)

passTrackSelections

public static int passTrackSelections(org.lcsim.event.Track track,
                                      List<org.lcsim.event.Track> tracklist,
                                      EventQuality.Quality trk_quality)

isTopTrack

public static boolean isTopTrack(org.lcsim.fit.helicaltrack.HelicalTrackFit htf)

isBottomTrack

public static boolean isBottomTrack(org.lcsim.fit.helicaltrack.HelicalTrackFit htf)

getHTF

public static org.lcsim.fit.helicaltrack.HelicalTrackFit getHTF(org.lcsim.event.MCParticle mcp,
                                                                double Bz)

Transform MCParticle into a Helix object. Note that it produces the helix parameters at nominal x=0 and assumes that there is no field at x<0

Parameters:

mcp - MC particle to be transformed

Returns:

HelicalTrackFit object based on the MC particle

getHTF

public static org.lcsim.fit.helicaltrack.HelicalTrackFit getHTF(org.lcsim.event.MCParticle mcp,
                                                                Hep3Vector origin,
                                                                double Bz)

Transform MCParticle into a HelicalTrackFit object. Note that it produces the HelicalTrackFit parameters at nominal reference point at origin and assumes that there is no field at x<0

Parameters:

mcp - - MC particle to be transformed

origin - - origin to be used for the track

Returns:

HelicalTrackFit object based on the MC particle

getHTF

public static org.lcsim.fit.helicaltrack.HelicalTrackFit getHTF(org.lcsim.event.Track track)

getHTF

public static org.lcsim.fit.helicaltrack.HelicalTrackFit getHTF(double[] par)

getHTF

public static org.lcsim.fit.helicaltrack.HelicalTrackFit getHTF(org.lcsim.event.TrackState state)

getMatchedTruthParticle

public static org.lcsim.event.MCParticle getMatchedTruthParticle(org.lcsim.event.Track track)

makeSeedTrackFromBaseTrack

public static org.lcsim.recon.tracking.seedtracker.SeedTrack makeSeedTrackFromBaseTrack(org.lcsim.event.Track track)

makeHelicalTrackHitFromTrackerHit

public static org.lcsim.fit.helicaltrack.HelicalTrackHit makeHelicalTrackHitFromTrackerHit(org.lcsim.event.TrackerHit hit)

getHitToStripsTable

public static org.lcsim.event.RelationalTable getHitToStripsTable(org.lcsim.event.EventHeader event)

getHitToRotatedTable

public static org.lcsim.event.RelationalTable getHitToRotatedTable(org.lcsim.event.EventHeader event)

getTrackTime

public static double getTrackTime(org.lcsim.event.Track track,
                                  org.lcsim.event.RelationalTable hitToStrips,
                                  org.lcsim.event.RelationalTable hitToRotated)

getTrackTimeSD

public static double getTrackTimeSD(org.lcsim.event.Track track,
                                    org.lcsim.event.RelationalTable hitToStrips,
                                    org.lcsim.event.RelationalTable hitToRotated)

getStripHits

public static List<org.lcsim.event.TrackerHit> getStripHits(org.lcsim.event.Track track,
                                                            org.lcsim.event.RelationalTable hitToStrips,
                                                            org.lcsim.event.RelationalTable hitToRotated)

sortHits

public static List<org.lcsim.event.TrackerHit> sortHits(Collection<org.lcsim.event.TrackerHit> hits)

numberOfSharedStrips

public static int numberOfSharedStrips(org.lcsim.event.Track track1,
                                       org.lcsim.event.Track track2,
                                       org.lcsim.event.RelationalTable hitToStrips,
                                       org.lcsim.event.RelationalTable hitToRotated)

Number of 2D hits shared between tracks

Parameters:

track1 -

track2 -

hitToStrips -

hitToRotated -

Returns:

hasSharedStrips

public static boolean hasSharedStrips(org.lcsim.event.Track track1,
                                      org.lcsim.event.Track track2,
                                      org.lcsim.event.RelationalTable hitToStrips,
                                      org.lcsim.event.RelationalTable hitToRotated)

Tells if there are shared 2D hits between tracks

Parameters:

track1 -

track2 -

hitToStrips -

hitToRotated -

Returns:

getLayer

public static int getLayer(org.lcsim.event.TrackerHit strip)

getIsolation

public static double getIsolation(org.lcsim.event.TrackerHit strip,
                                  org.lcsim.event.TrackerHit otherStrip,
                                  org.lcsim.event.RelationalTable hitToStrips,
                                  org.lcsim.event.RelationalTable hitToRotated)

Compute strip isolation, defined as the minimum distance to another strip in the same sensor. Strips are only checked if they formed valid crosses with the other strip in the cross (passing time and tolerance cuts).

Parameters:

strip - The strip whose isolation is being calculated.

otherStrip - The other strip in the stereo hit.

hitToStrips -

hitToRotated -

Returns:

Double_MAX_VALUE if no other strips found.

getIsolations

public static Double[] getIsolations(org.lcsim.event.Track trk,
                                     org.lcsim.event.RelationalTable hitToStrips,
                                     org.lcsim.event.RelationalTable hitToRotated,
                                     int layers)

Make an array with isolations for all 12 strip layers. If the track doesn't use a layer, the isolation is null; if there is no other strip hit in a layer, the isolation is Double.MAX_VALUE.

Parameters:

trk -

hitToStrips -

hitToRotated -

Returns:

isolations for all 12 strip layers

getIsolations

public static Double[] getIsolations(org.lcsim.event.Track trk,
                                     org.lcsim.event.RelationalTable hitToStrips,
                                     org.lcsim.event.RelationalTable hitToRotated)

extrapolateTrackUsingFieldMap

public static org.lcsim.event.base.BaseTrackState extrapolateTrackUsingFieldMap(org.lcsim.event.Track track,
                                                                                double startPositionX,
                                                                                double endPositionX,
                                                                                double stepSize,
                                                                                org.lcsim.geometry.FieldMap fieldMap)

Backward compatibility function for extrapolateTrackUsingFieldMap .

extrapolateTrackUsingFieldMap

public static org.lcsim.event.base.BaseTrackState extrapolateTrackUsingFieldMap(org.lcsim.event.TrackState track,
                                                                                double startPositionX,
                                                                                double endPosition,
                                                                                double stepSize,
                                                                                org.lcsim.geometry.FieldMap fieldMap)

Iteratively extrapolates a track to a specified value of x (z in detector frame) using the full 3D field map.

Parameters:

track - The Track object to extrapolate.

startPositionX - The position from which to start the extrapolation from. The track will be extrapolated to this point using a constant field.

endPositionX - The position to extrapolate the track to.

stepSize - The step size determining how far a track will be extrapolated after every iteration.

fieldMap - The 3D field map

Returns:

A TrackState at the final extrapolation point. Note that the "Tracking" frame is used for the reference point coordinate system.

extrapolateTrackUsingFieldMapRK

public static Pair<Hep3Vector,Hep3Vector> extrapolateTrackUsingFieldMapRK(org.lcsim.event.TrackState ts,
                                                                          Hep3Vector startPosition,
                                                                          double distanceZ,
                                                                          double stepSize,
                                                                          org.lcsim.geometry.FieldMap fM)

extrapolateTrackUsingFieldMap

public static org.lcsim.event.base.BaseTrackState extrapolateTrackUsingFieldMap(org.lcsim.event.TrackState track,
                                                                                double startPositionX,
                                                                                double endPosition,
                                                                                double stepSize,
                                                                                double epsilon,
                                                                                org.lcsim.geometry.FieldMap fieldMap)

calculatePhi

public static double calculatePhi(double x,
                                  double y,
                                  double xc,
                                  double yc,
                                  double sign)

calculatePhi

public static double calculatePhi(double px,
                                  double py)

calculateTanLambda

public static double calculateTanLambda(double pz,
                                        double p)

calculateCurvature

public static double calculateCurvature(double p,
                                        double q,
                                        double B)

getTrajectory

public static org.lcsim.util.swim.Trajectory getTrajectory(Hep3Vector p0,
                                                           org.lcsim.spacegeom.SpacePoint r0,
                                                           double q,
                                                           double B)

Parameters:

p0 -

r0 -

q -

B -

Returns:

the created trajectory

getTrackStateAtLocation

public static org.lcsim.event.TrackState getTrackStateAtLocation(org.lcsim.event.Track trk,
                                                                 int location)

Port of Track.getTrackState(int location) from the C++ LCIO API.

Parameters:

trk - A track.

location - A TrackState location constant

Returns:

The first matching TrackState; null if none is found.

getTrackStateAtECal

public static org.lcsim.event.TrackState getTrackStateAtECal(org.lcsim.event.Track trk)

getBField

public static Hep3Vector getBField(org.lcsim.geometry.Detector detector)

Sec

public static double Sec(double arg)

Power

public static double Power(double arg,
                           double pow)

Cos

public static double Cos(double arg)

Sin

public static double Sin(double arg)

Tan

public static double Tan(double arg)

ArcTan

public static double ArcTan(double x,
                            double y)

getMomentum

public static Hep3Vector getMomentum(double omega,
                                     double phi0,
                                     double tanL,
                                     double magneticField)

writeTrajectory

public static void writeTrajectory(Hep3Vector momentum,
                                   Hep3Vector position,
                                   int q,
                                   double bField,
                                   String outFile)
                            throws FileNotFoundException

Throws:

FileNotFoundException

globalToSensor

public static Hep3Vector globalToSensor(Hep3Vector trkpos,
                                        HpsSiSensor sensor)