Reconstructed track at TOP. More...

#include <TOPTrack.h>

Classes
struct	AssumedEmission
	assumed photon emission point in local frame More...

struct	SelectedHit
	selected photon hit from TOPDigits More...

struct	TrackAngles
	Sine and cosine of track polar and azimuthal angles at assumed photon emission. More...

Public Member Functions
	TOPTrack ()
	Default constructor.

	TOPTrack (const Track &track, const std::string &digitsName="", const Const::ChargedStable &chargedStable=Const::pion)
	Constructor from mdst track - isValid() must be checked before using the object.

	TOPTrack (const ExtHit *extHit, const std::string &digitsName="")
	Constructor from extrapolated track hit - isValid() must be checked before using the object.

bool	overrideTransformation (const ROOT::Math::Transform3D &transform)
	Overrides transformation from local to nominal frame, which is by default obtained from DB.

bool	isValid () const
	Checks if track is successfully constructed.

int	getModuleID () const
	Returns slot ID.

double	getMomentumMag () const
	Returns momentum magnitude (extrapolated to TOP)

double	getTransverseMomentum () const
	Returns transverse momentum (at POCA)

double	getCharge () const
	Returns charge.

double	getTrackLength () const
	Returns track length from IP to the average position of photon emission within quartz.

double	getLengthInQuartz () const
	Returns track length within quartz.

double	getBeta (const Const::ChargedStable &particle, double overrideMass=0) const
	Returns particle beta.

const TOPTrack::AssumedEmission &	getEmissionPoint (double dL=0) const
	Returns assumed photon emission position and track direction.

double	getTOF (const Const::ChargedStable &particle, double dL=0, double overrideMass=0) const
	Returns time-of-flight from IP to photon emission position.

const TOP::HelixSwimmer &	getHelix () const
	Returns helix Helix is given in nominal slot frame and with reference position at average photon emission.

const Track *	getTrack () const
	Returns mdst track.

const ExtHit *	getExtHit () const
	Returns extrapolated hit (track entrance to the bar)

const MCParticle *	getMCParticle () const
	Returns MC particle assigned to this track (if any)

int	getPDGCode () const
	Returns PDG code of associated MCParticle (returns 0 if none)

const TOPBarHit *	getBarHit () const
	Returns bar hit of MC particle assigned to this track (if any)

const std::vector< SelectedHit > &	getSelectedHits () const
	Returns selected photon hits from TOPDigits belonging to the slot ID.

double	getBkgRate () const
	Returns estimated background hit rate.

bool	isScanRequired (unsigned col, double time, double wid) const
	Checks if scan method of YScanner is needed to construct PDF for a given pixel column.

Private Member Functions
void	set (const Track &track, const std::string &digitsName, const Const::ChargedStable &chargedStable)
	Sets the object (called by constructors)

bool	setHelix (const ROOT::Math::Transform3D &transform)
	Sets helix (helix is given in nominal frame)

bool	xsecPrism (std::vector< double > &lengths, std::vector< ROOT::Math::XYZPoint > &positions, const RaytracerBase::Prism &prism, const ROOT::Math::Transform3D &transform)
	Calculates intersection of trajectory with prism.

Private Attributes
int	m_moduleID = 0
	slot ID

double	m_momentum = 0
	track momentum magnitude at TOP

double	m_pT = 0
	transverse momentum at POCA

double	m_charge = 0
	track charge in units of elementary charge

double	m_TOFLength = 0
	trajectory length corresponding to TOF of extrapolated hit

double	m_trackLength = 0
	trajectory length from IP to average photon emission point

double	m_length = 0
	trajectory length within quartz

TOP::HelixSwimmer	m_helix
	trajectory helix in nominal slot frame

DBObjPtr< TOPCalModuleAlignment >	m_alignment
	module alignment constants

DBObjPtr< TOPFrontEndSetting >	m_feSetting
	front-end settings

const Track *	m_track = 0
	mdst track

const ExtHit *	m_extHit = 0
	extrapolated hit

const MCParticle *	m_mcParticle = 0
	MC particle.

const TOPBarHit *	m_barHit = 0
	bar hit

bool	m_valid = false
	true for properly defined track

std::vector< SelectedHit >	m_selectedHits
	selected photon hits from TOPDigits belonging to this slot ID

double	m_bkgRate = 0
	estimated background hit rate

std::unordered_multimap< unsigned, const SelectedHit * >	m_columnHits
	selected hits mapped to pixel columns

std::map< double, TOPTrack::AssumedEmission >	m_emissionPoints
	assumed emission points in module local frame

Detailed Description

Reconstructed track at TOP.

Definition at line 39 of file TOPTrack.h.

Constructor & Destructor Documentation

◆ TOPTrack() [1/3]

TOPTrack ( )

inline

Default constructor.

Definition at line 102 of file TOPTrack.h.

103 {}

◆ TOPTrack() [2/3]

TOPTrack	(	const Track &	track,
		const std::string &	digitsName = `""`,
		const Const::ChargedStable &	chargedStable = `Const::pion`
	)

explicit

Constructor from mdst track - isValid() must be checked before using the object.

Parameters

track	mdst track
digitsName	name of TOPDigits collection
chargedStable	hypothesis used in mdst track extrapolation

Definition at line 40 of file TOPTrack.cc.

    {
      // find extrapolated track hit at TOP
 
      Const::EDetector myDetID = Const::EDetector::TOP;
      const auto* geo = TOPGeometryPar::Instance()->getGeometry();
      int numModules = geo->getNumModules();
      int pdgCode = std::abs(chargedStable.getPDGCode());
 
      RelationVector<ExtHit> extHits = track.getRelationsWith<ExtHit>();
      double tmin = 1e10; // some large time
      for (const auto& extHit : extHits) {
        if (std::abs(extHit.getPdgCode()) != pdgCode) continue;
        if (extHit.getDetectorID() != myDetID) continue;
        if (extHit.getCopyID() < 1 or extHit.getCopyID() > numModules) continue;
        if (extHit.getTOF() < tmin) {
          tmin = extHit.getTOF();
          m_extHit = &extHit;
        }
      }
      if (not m_extHit) return;
 
      // set the object
 
      set(track, digitsName, chargedStable);
    }

◆ TOPTrack() [3/3]

TOPTrack	(	const ExtHit *	extHit,
		const std::string &	digitsName = `""`
	)

explicit

Constructor from extrapolated track hit - isValid() must be checked before using the object.

Parameters

extHit	extrapolated track hit
digitsName	name of TOPDigits collection

Definition at line 68 of file TOPTrack.cc.

    {
      if (not extHit) return;
      m_extHit = extHit;
 
      const auto* track = extHit->getRelated<Track>();
      if (not track) {
        B2ERROR("TOPTrack: no related Track found for valid ExtHit");
        return;
      }
 
      auto chargedStable = Const::chargedStableSet.find(std::abs(extHit->getPdgCode()));
      if (chargedStable == Const::invalidParticle) {
        B2ERROR("TOPTrack: extrapolation hypothesis of ExtHit is not ChargedStable");
        return;
      }
 
      // set the object
 
      set(*track, digitsName, chargedStable);
    }

Member Function Documentation

◆ getBarHit()

const TOPBarHit * getBarHit ( ) const

inline

Returns bar hit of MC particle assigned to this track (if any)

Returns: bar hit or NULL pointer

Definition at line 238 of file TOPTrack.h.

238{return m_barHit;}

Belle2::TOP::TOPTrack::m_barHit

const TOPBarHit * m_barHit

bar hit

Definition: TOPTrack.h:303

◆ getBeta()

double getBeta	(	const Const::ChargedStable &	particle,
		double	overrideMass = `0`
	)		const

inline

Returns particle beta.

Parameters

particle	particle mass hypothesis
overrideMass	alternative mass value to be used to calculate beta. Ignored if <= 0.

Returns: particle beta

Definition at line 316 of file TOPTrack.h.

    {
      double mass = chargedStable.getMass();
      if (overrideMass > 0)
        mass = overrideMass;
      return m_momentum / sqrt(m_momentum * m_momentum + mass * mass);
    }

◆ getBkgRate()

double getBkgRate ( ) const

inline

Returns estimated background hit rate.

Returns: background hit rate per module

Definition at line 250 of file TOPTrack.h.

250{return m_bkgRate;}

Belle2::TOP::TOPTrack::m_bkgRate

double m_bkgRate

estimated background hit rate

Definition: TOPTrack.h:307

◆ getCharge()

double getCharge ( ) const

inline

Returns charge.

Returns: charge in units of elementary charge

Definition at line 161 of file TOPTrack.h.

161{return m_charge;}

Belle2::TOP::TOPTrack::m_charge

double m_charge

track charge in units of elementary charge

Definition: TOPTrack.h:292

◆ getEmissionPoint()

const TOPTrack::AssumedEmission & getEmissionPoint ( double dL = 0 ) const

Returns assumed photon emission position and track direction.

Parameters

dL	length difference to the average emission position

Returns: assumed photon emission position and track direction in local frame

Definition at line 357 of file TOPTrack.cc.

    {
      if (m_emissionPoints.size() > 1000) m_emissionPoints.clear(); // prevent blow-up
 
      auto& emissionPoint = m_emissionPoints[dL];
      if (not emissionPoint.isSet) {
        emissionPoint.position = m_helix.getPosition(dL);
        emissionPoint.trackAngles = TrackAngles(m_helix.getDirection(dL));
        emissionPoint.isSet = true;
      }
      return emissionPoint;
    }

◆ getExtHit()

const ExtHit * getExtHit ( ) const

inline

Returns extrapolated hit (track entrance to the bar)

Returns: extrapolated hit

Definition at line 216 of file TOPTrack.h.

216{return m_extHit;}

◆ getHelix()

const TOP::HelixSwimmer & getHelix ( ) const

inline

Returns helix Helix is given in nominal slot frame and with reference position at average photon emission.

Returns: helix

Definition at line 204 of file TOPTrack.h.

204{return m_helix;}

◆ getLengthInQuartz()

double getLengthInQuartz ( ) const

inline

Returns track length within quartz.

Returns: track length within quartz

Definition at line 173 of file TOPTrack.h.

173{return m_length;}

Belle2::TOP::TOPTrack::m_length

double m_length

trajectory length within quartz

Definition: TOPTrack.h:295

◆ getMCParticle()

const MCParticle * getMCParticle ( ) const

inline

Returns MC particle assigned to this track (if any)

Returns: MC particle or NULL pointer

Definition at line 222 of file TOPTrack.h.

222{return m_mcParticle;}

Belle2::TOP::TOPTrack::m_mcParticle

const MCParticle * m_mcParticle

MC particle.

Definition: TOPTrack.h:302

◆ getModuleID()

int getModuleID ( ) const

inline

Returns slot ID.

Returns: slot ID

Definition at line 143 of file TOPTrack.h.

143{return m_moduleID;}

Belle2::TOP::TOPTrack::m_moduleID

int m_moduleID

slot ID

Definition: TOPTrack.h:289

◆ getMomentumMag()

double getMomentumMag ( ) const

inline

Returns momentum magnitude (extrapolated to TOP)

Returns: momentum magnitude

Definition at line 149 of file TOPTrack.h.

149{return m_momentum;}

◆ getPDGCode()

int getPDGCode ( ) const

inline

Returns PDG code of associated MCParticle (returns 0 if none)

Returns: PDG code or 0

Definition at line 228 of file TOPTrack.h.

      {
        if (m_mcParticle) return m_mcParticle->getPDG();
        return 0;
      }

◆ getSelectedHits()

const std::vector< SelectedHit > & getSelectedHits ( ) const

inline

Returns selected photon hits from TOPDigits belonging to the slot ID.

Returns: selected photon hits

Definition at line 244 of file TOPTrack.h.

244{return m_selectedHits;}

Belle2::TOP::TOPTrack::m_selectedHits

std::vector< SelectedHit > m_selectedHits

selected photon hits from TOPDigits belonging to this slot ID

Definition: TOPTrack.h:306

◆ getTOF()

double getTOF	(	const Const::ChargedStable &	particle,
		double	dL = `0`,
		double	overrideMass = `0`
	)		const

inline

Returns time-of-flight from IP to photon emission position.

Parameters

particle	particle mass hypothesis
dL	length difference to the average emission position
overrideMass	alternative mass value to be used to calculate beta. Ignored if <= 0.

Returns: time-of-flight

Definition at line 324 of file TOPTrack.h.

    {
      return (m_trackLength + dL) / getBeta(chargedStable, overrideMass) / Const::speedOfLight;
    }

◆ getTrack()

const Track * getTrack ( ) const

inline

Returns mdst track.

Returns: mdst track

Definition at line 210 of file TOPTrack.h.

210{return m_track;}

Belle2::TOP::TOPTrack::m_track

const Track * m_track

mdst track

Definition: TOPTrack.h:300

◆ getTrackLength()

double getTrackLength ( ) const

inline

Returns track length from IP to the average position of photon emission within quartz.

Returns: track length from IP

Definition at line 167 of file TOPTrack.h.

167{return m_trackLength;}

◆ getTransverseMomentum()

double getTransverseMomentum ( ) const

inline

Returns transverse momentum (at POCA)

Returns: transverse momentum

Definition at line 155 of file TOPTrack.h.

155{return m_pT;}

Belle2::TOP::TOPTrack::m_pT

double m_pT

transverse momentum at POCA

Definition: TOPTrack.h:291

◆ isScanRequired()

bool isScanRequired	(	unsigned	col,
		double	time,
		double	wid
	)		const

Checks if scan method of YScanner is needed to construct PDF for a given pixel column.

Parameters

col	pixel column (0-based)
time	PDF peak time
wid	PDF peak width squared

Returns: true if at least one of the detected photons is within the PDF peak convoluted with TTS.

Definition at line 370 of file TOPTrack.cc.

    {
      const auto& tts = TOPGeometryPar::Instance()->getGeometry()->getTTS(0); // PMT independent TTS should be fine here
      const auto& range = m_columnHits.equal_range(col);
      for (auto it = range.first; it != range.second; ++it) {
        const auto hit = it->second;
        for (const auto& gaus : tts.getTTS()) {
          double sigsq = wid + pow(gaus.sigma, 2) + pow(hit->timeErr, 2);
          double x = pow(hit->time - time - gaus.position, 2) / sigsq;
          if (x < 10) return true;
        }
      }
 
      return false;
    }

◆ isValid()

bool isValid ( ) const

inline

Checks if track is successfully constructed.

Returns: true if successful

Definition at line 137 of file TOPTrack.h.

137{return m_valid;}

Belle2::TOP::TOPTrack::m_valid

bool m_valid

true for properly defined track

Definition: TOPTrack.h:304

◆ overrideTransformation()

bool overrideTransformation ( const ROOT::Math::Transform3D & transform )

inline

Overrides transformation from local to nominal frame, which is by default obtained from DB.

Needed for module alignment.

Parameters

transform transformation from local to nominal frame

Returns: true if quartz bar intersection still exists

Definition at line 127 of file TOPTrack.h.

      {
        m_valid = setHelix(transform);
        return m_valid;
      }

◆ set()

void set	(	const Track &	track,
		const std::string &	digitsName,
		const Const::ChargedStable &	chargedStable
	)

private

Sets the object (called by constructors)

Parameters

track	mdst track
digitsName	name of TOPDigits collection
chargedStable	hypothesis used in mdst track extrapolation

Definition at line 91 of file TOPTrack.cc.

    {
      // require fitResult
 
      const auto* fitResult = track.getTrackFitResultWithClosestMass(chargedStable);
      if (not fitResult) {
        B2ERROR("TOPTrack: no TrackFitResult available for PDGCode = " << chargedStable.getPDGCode());
        return;
      }
      m_pT = fitResult->getTransverseMomentum();
 
      // require hits in CDC
 
      if (fitResult->getHitPatternCDC().getNHits() == 0) return;
 
      // set pointers
 
      m_track = &track;
      m_mcParticle = track.getRelated<MCParticle>();
      if (m_mcParticle) {
        const auto barHits = m_mcParticle->getRelationsWith<TOPBarHit>();
        for (const auto& barHit : barHits) {
          if (barHit.getModuleID() == m_extHit->getCopyID()) m_barHit = &barHit;
        }
      }
 
      // set track parameters and helix
 
      m_moduleID = m_extHit->getCopyID();
      m_momentum = m_extHit->getMomentum().R();
      m_charge = fitResult->getChargeSign();
      m_TOFLength = m_extHit->getTOF() * Const::speedOfLight * getBeta(chargedStable);
      m_valid = setHelix(m_alignment->getTransformation(m_moduleID));
      if (not m_valid) return;
 
      // selection of photon hits belonging to this track
 
      unsigned numHitsOtherSlots = 0;
      StoreArray<TOPDigit> digits(digitsName);
      for (const auto& digit : digits) {
        if (digit.getHitQuality() != TOPDigit::c_Good) continue;
        if (digit.getModuleID() == m_moduleID) {
          m_selectedHits.push_back(SelectedHit(digit.getPixelID(), digit.getTime(), digit.getTimeError()));
        } else {
          numHitsOtherSlots++;
        }
      }
 
      // background rate estimation (TODO to be improved ...)
 
      const auto* geo = TOPGeometryPar::Instance()->getGeometry();
      const auto& tdc = geo->getNominalTDC();
      double timeWindow = m_feSetting->getReadoutWindows() * tdc.getSyncTimeBase() / static_cast<double>(TOPNominalTDC::c_syncWindows);
 
      const auto& backgroundPDFs = TOPRecoManager::getBackgroundPDFs();
      unsigned k = m_moduleID - 1;
      double effi = (k < backgroundPDFs.size()) ? backgroundPDFs[k].getEfficiency() : 0;
      double effiSum = 0;
      for (const auto& bkg : backgroundPDFs) effiSum += bkg.getEfficiency();
      m_bkgRate = (effiSum > effi) ? numHitsOtherSlots * effi / (effiSum - effi) / timeWindow : 0;
 
      // selected photon hits mapped to pixel columns
 
      const auto* yScanner = TOPRecoManager::getYScanner(m_moduleID);
      if (not yScanner) {
        B2ERROR("TOPTrack: YScanner for given module not found (must be a bug!)" << LogVar("slot", m_moduleID));
        m_valid = false;
        return;
      }
      unsigned numCols = yScanner->getPixelPositions().getNumPixelColumns();
      for (const auto& hit : m_selectedHits) {
        unsigned col = (hit.pixelID - 1) % numCols;
        m_columnHits.emplace(col, &hit);
      }
 
      m_valid = effi > 0; // no sense to provide PID for the track if the module is fully inefficient
    }

◆ setHelix()

bool setHelix ( const ROOT::Math::Transform3D & transform )

private

Sets helix (helix is given in nominal frame)

Parameters

transform transformation from local to nominal frame

Returns: true if quartz bar intersection exists

Definition at line 170 of file TOPTrack.cc.

    {
      m_emissionPoints.clear();
 
      // helix in module nominal frame (z-axis still parallel to magnetic field)
 
      const auto* geo = TOPGeometryPar::Instance()->getGeometry();
      const auto& module = geo->getModule(m_moduleID);
      auto globalPosition = m_extHit->getPosition();
      auto position = module.pointGlobalToNominal(static_cast<XYZPoint>(globalPosition));
      auto momentum = module.momentumGlobalToNominal(m_extHit->getMomentum());
      double Bz = BFieldManager::getField(globalPosition).Z();
      m_helix.set(position, momentum, m_charge, Bz);
      m_helix.setTransformation(transform);
 
      // geometry data
 
      const RaytracerBase::BarSegment bar(module);
      const RaytracerBase::Mirror mirror(module);
 
      // bar surfaces in module nominal frame
 
      std::vector<XYZPoint> points;
      std::vector<XYZVector> normals;
      points.push_back(transform * XYZPoint(0, -bar.B / 2, 0)); // lower
      normals.push_back(transform * XYZVector(0, -1, 0));
 
      points.push_back(transform * XYZPoint(0, bar.B / 2, 0)); // upper
      normals.push_back(transform * XYZVector(0, 1, 0));
 
      points.push_back(transform * XYZPoint(-bar.A / 2, 0, 0)); // left side
      normals.push_back(transform * XYZVector(-1, 0, 0));
 
      points.push_back(transform * XYZPoint(bar.A / 2, 0, 0)); // right side
      normals.push_back(transform * XYZVector(1, 0, 0));
 
      // intersection with quartz bar
 
      std::vector<double> lengths; // w.r.t extHit position
      std::vector<XYZPoint> positions; // in module local frame
      for (size_t i = 0; i < 2; i++) {
        double t = m_helix.getDistanceToPlane(points[i], normals[i]);
        if (isnan(t)) return false;
        auto r = m_helix.getPosition(t);
        if (std::abs(r.X()) > bar.A / 2) {
          auto k = (r.X() > 0) ? 3 : 2;
          t = m_helix.getDistanceToPlane(points[k], normals[k]);
          if (isnan(t)) return false;
          r = m_helix.getPosition(t);
          if (r.Z() >= bar.zL and std::abs(r.Y()) > bar.B / 2) return false;
        }
        lengths.push_back(t);
        positions.push_back(r);
      }
 
      // crossing prism?
 
      if (positions[0].Z() < bar.zL or positions[1].Z() < bar.zL) {
        const RaytracerBase::Prism prism(module);
        bool ok = xsecPrism(lengths, positions, prism, transform);
        if (not ok) return false;
      }
 
      // crossing mirror surface?
 
      std::vector<bool> outOfBar;
      XYZPoint rc(mirror.xc, mirror.yc, mirror.zc);
      double Rsq = mirror.R * mirror.R;
      for (const auto& r : positions) {
        outOfBar.push_back((r - rc).Mag2() > Rsq);
      }
      if (outOfBar[0] and outOfBar[1]) return false;
 
      if (outOfBar[0] or outOfBar[1]) { // track crosses mirror surface, where?
        if (outOfBar[0]) std::reverse(lengths.begin(), lengths.end());
        double t1 = lengths[0];
        double t2 = lengths[1];
        for (int i = 0; i < 20; i++) {
          double t = (t1 + t2) / 2;
          auto r = m_helix.getPosition(t);
          if ((r - rc).Mag2() > Rsq) {
            t2 = t;
          } else {
            t1 = t;
          }
        }
        lengths[1] = (t1 + t2) / 2;
      }
 
      // set track length in quartz and full length from IP to the average emission point
 
      m_length = std::abs(lengths[1] - lengths[0]);
      double length = (lengths[0] + lengths[1]) / 2;
      m_trackLength = m_TOFLength + length;
 
      // finally move reference position to average emission point
 
      m_helix.moveReferencePosition(length);
 
      return m_length > bar.B / 2; // require minimal track lenght inside quartz (more than half of bar thickness)
    }

◆ xsecPrism()

bool xsecPrism	(	std::vector< double > &	lengths,
		std::vector< ROOT::Math::XYZPoint > &	positions,
		const RaytracerBase::Prism &	prism,
		const ROOT::Math::Transform3D &	transform
	)

private

Calculates intersection of trajectory with prism.

Parameters

lengths	trajectory lengths relative to extHit position of intersection points [in/out]
positions	positions of intersection points in module local frame [in/out]
prism	prism geometry data
transform	transformation from local to nominal frame

Returns: true if intersection exists

Definition at line 273 of file TOPTrack.cc.

    {
      std::vector<XYZPoint> points;
      std::vector<XYZVector> normals;
 
      points.push_back(transform * XYZPoint(0, prism.yDown, 0)); // lower-most surface
      normals.push_back(transform * XYZVector(0, -1, 0));
 
      points.push_back(transform * XYZPoint(0, prism.yUp, 0)); // upper surface
      normals.push_back(transform * XYZVector(0, 1, 0));
 
      points.push_back(transform * XYZPoint(-prism.A / 2, 0, 0)); // left-side surface
      normals.push_back(transform * XYZVector(-1, 0, 0));
 
      points.push_back(transform * XYZPoint(prism.A / 2, 0, 0)); // right-side surface
      normals.push_back(transform * XYZVector(1, 0, 0));
 
      for (size_t i = 0; i < 2; i++) {
        if (positions[i].Z() < prism.zR) {
          double t = m_helix.getDistanceToPlane(points[i], normals[i]);
          if (isnan(t)) return false;
          auto r = m_helix.getPosition(t);
          if (i == 0 and r.Z() > prism.zFlat) { // intersection with slanted surface -> find it using bisection
            auto point = transform * XYZPoint(0, -prism.B / 2, 0);
            double t1 = m_helix.getDistanceToPlane(point, normals[0]);
            if (isnan(t1)) return false;
            double t2 = t;
            for (int iter = 0; iter < 20; iter++) {
              t = (t1 + t2) / 2;
              r = m_helix.getPosition(t);
              double ySlanted = prism.yDown + prism.slope * (r.Z() - prism.zFlat);
              if (r.Y() < ySlanted) t2 = t;
              else t1 = t;
            }
            t = (t1 + t2) / 2;
          }
          if (std::abs(r.X()) > prism.A / 2) { // intersection on the side surface
            auto k = (r.X() > 0) ? 3 : 2;
            t = m_helix.getDistanceToPlane(points[k], normals[k]);
            if (isnan(t)) return false;
            r = m_helix.getPosition(t);
            if (r.Z() < prism.zR) { // yes, it's on the prism side
              if (r.Y() > prism.yUp or r.Y() < prism.yDown) return false;
              double ySlanted = prism.yDown + prism.slope * (r.Z() - prism.zFlat);
              if (r.Y() < ySlanted) return false;
            } else { // no, it's on the prism entrance but outside the bar exit window -> find it using bisection
              double t1 = lengths[i];
              double t2 = t;
              for (int iter = 0; iter < 20; iter++) {
                t = (t1 + t2) / 2;
                r = m_helix.getPosition(t);
                if (r.Z() < prism.zR) t1 = t;
                else t2 = t;
              }
              t = (t1 + t2) / 2;
            }
          }
          lengths[i] = t;
          positions[i] = m_helix.getPosition(t);
        }
      }
 
      if (positions[0].Z() < prism.zL and positions[1].Z() < prism.zL) return false;
 
      if (positions[0].Z() < prism.zL or positions[1].Z() < prism.zL) { // intersection is on exit window
        int i0 = (positions[0].Z() < prism.zL) ? 0 : 1;
        double t1 = lengths[i0];
        double t2 = lengths[(i0 + 1) % 2];
        for (int iter = 0; iter < 20; iter++) {
          double t = (t1 + t2) / 2;
          auto r = m_helix.getPosition(t);
          if (r.Z() < prism.zL) t1 = t;
          else t2 = t;
        }
        double t = (t1 + t2) / 2;
        lengths[i0] = t;
        positions[i0] = m_helix.getPosition(t);
      }
 
      return true;
    }

Member Data Documentation

◆ m_alignment

DBObjPtr<TOPCalModuleAlignment> m_alignment

private

module alignment constants

Definition at line 297 of file TOPTrack.h.

◆ m_barHit

const TOPBarHit* m_barHit = 0

private

bar hit

Definition at line 303 of file TOPTrack.h.

◆ m_bkgRate

double m_bkgRate = 0

private

estimated background hit rate

Definition at line 307 of file TOPTrack.h.

◆ m_charge

double m_charge = 0

private

track charge in units of elementary charge

Definition at line 292 of file TOPTrack.h.

◆ m_columnHits

std::unordered_multimap<unsigned, const SelectedHit*> m_columnHits

private

selected hits mapped to pixel columns

Definition at line 308 of file TOPTrack.h.

◆ m_emissionPoints

std::map<double, TOPTrack::AssumedEmission> m_emissionPoints

mutableprivate

assumed emission points in module local frame

Definition at line 311 of file TOPTrack.h.

◆ m_extHit

const ExtHit* m_extHit = 0

private

extrapolated hit

Definition at line 301 of file TOPTrack.h.

◆ m_feSetting

DBObjPtr<TOPFrontEndSetting> m_feSetting

private

front-end settings

Definition at line 298 of file TOPTrack.h.

◆ m_helix

TOP::HelixSwimmer m_helix

private

trajectory helix in nominal slot frame

Definition at line 296 of file TOPTrack.h.

◆ m_length

double m_length = 0

private

trajectory length within quartz

Definition at line 295 of file TOPTrack.h.

◆ m_mcParticle

const MCParticle* m_mcParticle = 0

private

MC particle.

Definition at line 302 of file TOPTrack.h.

◆ m_moduleID

int m_moduleID = 0

private

slot ID

Definition at line 289 of file TOPTrack.h.

◆ m_momentum

double m_momentum = 0

private

track momentum magnitude at TOP

Definition at line 290 of file TOPTrack.h.

◆ m_pT

double m_pT = 0

private

transverse momentum at POCA

Definition at line 291 of file TOPTrack.h.

◆ m_selectedHits

std::vector<SelectedHit> m_selectedHits

private

selected photon hits from TOPDigits belonging to this slot ID

Definition at line 306 of file TOPTrack.h.

◆ m_TOFLength

double m_TOFLength = 0

private

trajectory length corresponding to TOF of extrapolated hit

Definition at line 293 of file TOPTrack.h.

◆ m_track

const Track* m_track = 0

private

mdst track

Definition at line 300 of file TOPTrack.h.

◆ m_trackLength

double m_trackLength = 0

private

trajectory length from IP to average photon emission point

Definition at line 294 of file TOPTrack.h.

◆ m_valid

bool m_valid = false

private

true for properly defined track

Definition at line 304 of file TOPTrack.h.

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

top/reconstruction_cpp/include/TOPTrack.h
top/reconstruction_cpp/src/TOPTrack.cc

Classes

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TOPTrack() [1/3]

◆ TOPTrack() [2/3]

◆ TOPTrack() [3/3]

Member Function Documentation

◆ getBarHit()

◆ getBeta()

◆ getBkgRate()

◆ getCharge()

◆ getEmissionPoint()

◆ getExtHit()

◆ getHelix()

◆ getLengthInQuartz()

◆ getMCParticle()

◆ getModuleID()

◆ getMomentumMag()

◆ getPDGCode()

◆ getSelectedHits()

◆ getTOF()

◆ getTrack()

◆ getTrackLength()

◆ getTransverseMomentum()

◆ isScanRequired()

◆ isValid()

◆ overrideTransformation()

◆ set()

◆ setHelix()

◆ xsecPrism()

Member Data Documentation

◆ m_alignment

◆ m_barHit

◆ m_bkgRate

◆ m_charge

◆ m_columnHits

◆ m_emissionPoints

◆ m_extHit

◆ m_feSetting

◆ m_helix

◆ m_length

◆ m_mcParticle

◆ m_moduleID

◆ m_momentum

◆ m_pT

◆ m_selectedHits

◆ m_TOFLength

◆ m_track

◆ m_trackLength

◆ m_valid