mkfit::MkFitter Class Reference

#include <MkFitter.h>

Inheritance diagram for mkfit::MkFitter:

Public Member Functions
void	checkAlignment ()
int	countInvalidHits (int itrack, int end_hit) const
int	countInvalidHits (int itrack) const
int	countValidHits (int itrack, int end_hit) const
int	countValidHits (int itrack) const
void	fitTracksWithInterSlurp (const std::vector< HitVec > &layersohits, const PropagationFlags &pflags, int N_proc)
void	inputSeedsTracksAndHits (const std::vector< Track > &seeds, const std::vector< Track > &tracks, const std::vector< HitVec > &layerHits, int beg, int end)
void	inputTracksAndHitIdx (const std::vector< Track > &tracks, int beg, int end, bool inputProp)
void	inputTracksAndHitIdx (const std::vector< std::vector< Track > > &tracks, const std::vector< std::pair< int, int > > &idxs, int beg, int end, bool inputProp)
void	inputTracksAndHits (const std::vector< Track > &tracks, const std::vector< HitVec > &layerHits, int beg, int end)
void	inputTracksAndHits (const std::vector< Track > &tracks, const std::vector< LayerOfHits > &layerHits, int beg, int end)
void	inputTracksForFit (const std::vector< Track > &tracks, int beg, int end)
	MkFitter ()
void	outputFittedTracks (std::vector< Track > &tracks, int beg, int end) const
void	outputFittedTracksAndHitIdx (std::vector< Track > &tracks, int beg, int end, bool outputProp) const
void	outputPropagatedTracks (std::vector< Track > &tracks, int beg, int end) const
void	outputTracks (std::vector< Track > &tracks, int beg, int end, int iCP) const
void	printPt (int idx)
MPlexLS &	refErr0 ()
MPlexLV &	refPar0 ()
void	setNhits (int newnhits)
void	slurpInTracksAndHits (const std::vector< Track > &tracks, const std::vector< HitVec > &layerHits, int beg, int end)
Public Member Functions inherited from mkfit::MkBase
void	clearFailFlag ()
float	getPar (int itrack, int i, int par) const
	MkBase ()
void	propagateTracksToHitR (const MPlexHV &par, const int N_proc, const PropagationFlags &pf, const MPlexQI *noMatEffPtr=nullptr)
void	propagateTracksToHitZ (const MPlexHV &par, const int N_proc, const PropagationFlags &pf, const MPlexQI *noMatEffPtr=nullptr)
void	propagateTracksToPCAZ (const int N_proc, const PropagationFlags &pf)
void	propagateTracksToR (float r, const int N_proc, const PropagationFlags &pf)
void	propagateTracksToZ (float z, const int N_proc, const PropagationFlags &pf)
float	radius (int itrack, int i) const
float	radiusSqr (int itrack, int i) const

Private Attributes
MPlexQI	m_CandIdx
MPlexQF	m_Chi2
MPlexQHoT	m_HoTArr [Config::nMaxTrkHits]
MPlexQI	m_Label
MPlexHS	m_msErr [Config::nMaxTrkHits]
MPlexHV	m_msPar [Config::nMaxTrkHits]
int	m_Nhits
MPlexQI	m_SeedIdx
MPlexHitIdx	m_XHitArr
MPlexQI	m_XHitSize

Additional Inherited Members
Static Public Attributes inherited from mkfit::MkBase
static constexpr int	iC = 0
static constexpr int	iP = 1
Protected Attributes inherited from mkfit::MkBase
MPlexQI	m_Chg
MPlexLS	m_Err [2]
MPlexQI	m_FailFlag
MPlexLV	m_Par [2]

Detailed Description

Definition at line 16 of file MkFitter.h.

Constructor & Destructor Documentation

MkFitter()

mkfit::MkFitter::MkFitter ( )

inline

Definition at line 18 of file MkFitter.h.

: m_Nhits(0) {}

Member Function Documentation

checkAlignment()

void mkfit::MkFitter::checkAlignment

(

)

Definition at line 13 of file MkFitter.cc.

References Matriplex::align_check(), mkfit::MkBase::m_Err, m_msErr, m_msPar, and mkfit::MkBase::m_Par.

                                {
    printf("MkFitter alignment check:\n");
    Matriplex::align_check("  m_Err[0]   =", &m_Err[0].fArray[0]);
    Matriplex::align_check("  m_Err[1]   =", &m_Err[1].fArray[0]);
    Matriplex::align_check("  m_Par[0]   =", &m_Par[0].fArray[0]);
    Matriplex::align_check("  m_Par[1]   =", &m_Par[1].fArray[0]);
    Matriplex::align_check("  m_msErr[0] =", &m_msErr[0].fArray[0]);
    Matriplex::align_check("  m_msPar[0] =", &m_msPar[0].fArray[0]);
  }

countInvalidHits() [1/2]

int mkfit::MkFitter::countInvalidHits	(	int	itrack,
		int	end_hit
	)		const

Definition at line 38 of file MkFitter.cc.

References m_HoTArr, and mps_fire::result.

                                                              {
    int result = 0;
    for (int hi = 0; hi < end_hit; ++hi) {
      // XXXX MT: Should also count -2 hits as invalid?
      if (m_HoTArr[hi](itrack, 0, 0).index == -1)
        result++;
    }
    return result;
  }

countInvalidHits() [2/2]

int mkfit::MkFitter::countInvalidHits ( int  itrack ) const

inline

Definition at line 35 of file MkFitter.h.

References countInvalidHits(), and m_Nhits.

Referenced by countInvalidHits().

{ return countInvalidHits(itrack, m_Nhits); }

countValidHits() [1/2]

int mkfit::MkFitter::countValidHits	(	int	itrack,
		int	end_hit
	)		const

Definition at line 29 of file MkFitter.cc.

References m_HoTArr, and mps_fire::result.

                                                            {
    int result = 0;
    for (int hi = 0; hi < end_hit; ++hi) {
      if (m_HoTArr[hi](itrack, 0, 0).index >= 0)
        result++;
    }
    return result;
  }

countValidHits() [2/2]

int mkfit::MkFitter::countValidHits ( int  itrack ) const

inline

Definition at line 34 of file MkFitter.h.

References countValidHits(), and m_Nhits.

Referenced by countValidHits().

{ return countValidHits(itrack, m_Nhits); }

fitTracksWithInterSlurp()

void mkfit::MkFitter::fitTracksWithInterSlurp	(	const std::vector< HitVec > &	layersohits,
		const PropagationFlags &	pflags,
		int	N_proc
	)

Definition at line 302 of file MkFitter.cc.

References mkfit::MatriplexErrParPackerSlurpIn< T, D >::addInput(), mkfit::MatriplexPackerSlurpIn< D >::addNullInput(), data, mps_fire::i, mkfit::MkBase::iC, cuy::ii, mkfit::MkBase::iP, mkfit::kalmanUpdate(), mkfit::MkBase::m_Err, m_HoTArr, m_msErr, m_msPar, m_Nhits, mkfit::MkBase::m_Par, mkfit::MatriplexErrParPackerSlurpIn< T, D >::pack(), and mkfit::MkBase::propagateTracksToHitR().

                                                           {
    // XXXX This has potential issues hits coming from different layers!
    // Expected to only work reliably with barrel (consecutive layers from 0 -> m_Nhits)
    // and with hits present on every layer for every track.

    // Loops over layers and:
    // a) slurps in hit parameters;
    // b) propagates and updates tracks

    for (int ii = 0; ii < m_Nhits; ++ii) {
      // XXXX Assuming hit index corresponds to layer!
      MatriplexHitPacker mhp(layersohits[ii].data());

      for (int i = 0; i < N_proc; ++i) {
        const int hidx = m_HoTArr[ii](i, 0, 0).index;
        const int hlyr = m_HoTArr[ii](i, 0, 0).layer;

        // XXXXMT4K What to do with hidx < 0 ????
        // This could solve the unbalanced fit.
        // Or, if the hidx is the "universal" missing hit, it could just work.
        // Say, hidx = 0 ... grr ... but then we don't know it is missing.

        if (hidx < 0 || hlyr < 0) {
          mhp.addNullInput();
        } else {
          mhp.addInput(layersohits[hlyr][hidx]);
        }
      }

      mhp.pack(m_msErr[0], m_msPar[0]);

      propagateTracksToHitR(m_msPar[0], N_proc, pflags);

      kalmanUpdate(m_Err[iP], m_Par[iP], m_msErr[0], m_msPar[0], m_Err[iC], m_Par[iC], N_proc);
    }
  }

inputSeedsTracksAndHits()

void mkfit::MkFitter::inputSeedsTracksAndHits	(	const std::vector< Track > &	seeds,
		const std::vector< Track > &	tracks,
		const std::vector< HitVec > &	layerHits,
		int	beg,
		int	end
	)

Definition at line 222 of file MkFitter.cc.

References mkfit::TrackBase::charge(), mkfit::TrackBase::chi2(), mkfit::TrackBase::errors(), mkfit::Track::getHitIdx(), mkfit::Track::getHitOnTrack(), mps_fire::i, mkfit::MkBase::iC, mkfit::TrackBase::label(), mkfit::MkBase::m_Chg, m_Chi2, mkfit::MkBase::m_Err, m_HoTArr, m_Label, m_msErr, m_msPar, m_Nhits, mkfit::MkBase::m_Par, mkfit::TrackBase::parameters(), HLT_FULL_cff::seeds, and pwdgSkimBPark_cfi::tracks.

                                                  {
    // Assign track parameters to initial state and copy hit values in.

    // This might not be true for the last chunk!
    // assert(end - beg == NN);

    int itrack;
    for (int i = beg; i < end; ++i) {
      itrack = i - beg;

      const Track& see = seeds[i];

      m_Label(itrack, 0, 0) = see.label();
      if (see.label() < 0)
        continue;

      m_Err[iC].copyIn(itrack, see.errors().Array());
      m_Par[iC].copyIn(itrack, see.parameters().Array());

      m_Chg(itrack, 0, 0) = see.charge();
      m_Chi2(itrack, 0, 0) = see.chi2();

      const Track& trk = tracks[see.label()];

      // CopyIn seems fast enough, but indirections are quite slow.
      for (int hi = 0; hi < m_Nhits; ++hi) {
        m_HoTArr[hi](itrack, 0, 0) = trk.getHitOnTrack(hi);

        const int hidx = trk.getHitIdx(hi);
        if (hidx < 0)
          continue;  //fixme, check if this is harmless

        const Hit& hit = layerHits[hi][hidx];
        m_msErr[hi].copyIn(itrack, hit.errArray());
        m_msPar[hi].copyIn(itrack, hit.posArray());
      }
    }
  }

inputTracksAndHitIdx() [1/2]

void mkfit::MkFitter::inputTracksAndHitIdx	(	const std::vector< Track > &	tracks,
		int	beg,
		int	end,
		bool	inputProp
	)

Definition at line 165 of file MkFitter.cc.

References mkfit::TrackBase::charge(), mkfit::TrackBase::chi2(), mkfit::TrackBase::errors(), mkfit::Track::getHitOnTrack(), mps_fire::i, mkfit::MkBase::iC, mkfit::MkBase::iP, mkfit::TrackBase::label(), mkfit::MkBase::m_Chg, m_Chi2, mkfit::MkBase::m_Err, m_HoTArr, m_Label, m_Nhits, mkfit::MkBase::m_Par, mkfit::TrackBase::parameters(), and pwdgSkimBPark_cfi::tracks.

                                                                                                      {
    // Assign track parameters to initial state and copy hit values in.

    // This might not be true for the last chunk!
    // assert(end - beg == NN);

    const int iI = inputProp ? iP : iC;

    int itrack = 0;
    for (int i = beg; i < end; ++i, ++itrack) {
      const Track& trk = tracks[i];

      m_Err[iI].copyIn(itrack, trk.errors().Array());
      m_Par[iI].copyIn(itrack, trk.parameters().Array());

      m_Chg(itrack, 0, 0) = trk.charge();
      m_Chi2(itrack, 0, 0) = trk.chi2();
      m_Label(itrack, 0, 0) = trk.label();

      for (int hi = 0; hi < m_Nhits; ++hi) {
        m_HoTArr[hi](itrack, 0, 0) = trk.getHitOnTrack(hi);
      }
    }
  }

inputTracksAndHitIdx() [2/2]

void mkfit::MkFitter::inputTracksAndHitIdx	(	const std::vector< std::vector< Track > > &	tracks,
		const std::vector< std::pair< int, int > > &	idxs,
		int	beg,
		int	end,
		bool	inputProp
	)

Definition at line 190 of file MkFitter.cc.

References mkfit::TrackBase::charge(), mkfit::TrackBase::chi2(), mkfit::TrackBase::errors(), mkfit::Track::getHitOnTrack(), mps_fire::i, mkfit::MkBase::iC, mkfit::MkBase::iP, mkfit::TrackBase::label(), m_CandIdx, mkfit::MkBase::m_Chg, m_Chi2, mkfit::MkBase::m_Err, m_HoTArr, m_Label, m_Nhits, mkfit::MkBase::m_Par, m_SeedIdx, mkfit::TrackBase::parameters(), and pwdgSkimBPark_cfi::tracks.

                                                      {
    // Assign track parameters to initial state and copy hit values in.

    // This might not be true for the last chunk!
    // assert(end - beg == NN);

    const int iI = inputProp ? iP : iC;

    int itrack = 0;
    for (int i = beg; i < end; ++i, ++itrack) {
      const Track& trk = tracks[idxs[i].first][idxs[i].second];

      m_Label(itrack, 0, 0) = trk.label();
      m_SeedIdx(itrack, 0, 0) = idxs[i].first;
      m_CandIdx(itrack, 0, 0) = idxs[i].second;

      m_Err[iI].copyIn(itrack, trk.errors().Array());
      m_Par[iI].copyIn(itrack, trk.parameters().Array());

      m_Chg(itrack, 0, 0) = trk.charge();
      m_Chi2(itrack, 0, 0) = trk.chi2();

      for (int hi = 0; hi < m_Nhits; ++hi) {
        m_HoTArr[hi](itrack, 0, 0) = trk.getHitOnTrack(hi);
      }
    }
  }

inputTracksAndHits() [1/2]

void mkfit::MkFitter::inputTracksAndHits	(	const std::vector< Track > &	tracks,
		const std::vector< HitVec > &	layerHits,
		int	beg,
		int	end
	)

Definition at line 50 of file MkFitter.cc.

References mkfit::TrackBase::charge(), mkfit::TrackBase::chi2(), mkfit::TrackBase::errors(), mkfit::Track::getHitIdx(), mkfit::Track::getHitOnTrack(), mps_fire::i, mkfit::MkBase::iC, mkfit::TrackBase::label(), mkfit::MkBase::m_Chg, m_Chi2, mkfit::MkBase::m_Err, m_HoTArr, m_Label, m_msErr, m_msPar, m_Nhits, mkfit::MkBase::m_Par, mkfit::TrackBase::parameters(), and pwdgSkimBPark_cfi::tracks.

                                             {
    // Assign track parameters to initial state and copy hit values in.

    // This might not be true for the last chunk!
    // assert(end - beg == NN);

    int itrack = 0;

    for (int i = beg; i < end; ++i, ++itrack) {
      const Track& trk = tracks[i];

      m_Err[iC].copyIn(itrack, trk.errors().Array());
      m_Par[iC].copyIn(itrack, trk.parameters().Array());

      m_Chg(itrack, 0, 0) = trk.charge();
      m_Chi2(itrack, 0, 0) = trk.chi2();
      m_Label(itrack, 0, 0) = trk.label();

      // CopyIn seems fast enough, but indirections are quite slow.
      for (int hi = 0; hi < m_Nhits; ++hi) {
        m_HoTArr[hi](itrack, 0, 0) = trk.getHitOnTrack(hi);

        const int hidx = trk.getHitIdx(hi);
        if (hidx < 0)
          continue;

        const Hit& hit = layerHits[hi][hidx];
        m_msErr[hi].copyIn(itrack, hit.errArray());
        m_msPar[hi].copyIn(itrack, hit.posArray());
      }
    }
  }

inputTracksAndHits() [2/2]

void mkfit::MkFitter::inputTracksAndHits	(	const std::vector< Track > &	tracks,
		const std::vector< LayerOfHits > &	layerHits,
		int	beg,
		int	end
	)

Definition at line 86 of file MkFitter.cc.

References mkfit::TrackBase::charge(), mkfit::TrackBase::chi2(), mkfit::TrackBase::errors(), mkfit::Track::getHitIdx(), mkfit::Track::getHitLyr(), mkfit::Track::getHitOnTrack(), mps_fire::i, mkfit::MkBase::iC, mkfit::TrackBase::label(), mkfit::MkBase::m_Chg, m_Chi2, mkfit::MkBase::m_Err, m_HoTArr, m_Label, m_msErr, m_msPar, m_Nhits, mkfit::MkBase::m_Par, mkfit::TrackBase::parameters(), and pwdgSkimBPark_cfi::tracks.

                                             {
    // Assign track parameters to initial state and copy hit values in.

    // This might not be true for the last chunk!
    // assert(end - beg == NN);

    int itrack;

    for (int i = beg; i < end; ++i) {
      itrack = i - beg;
      const Track& trk = tracks[i];

      m_Label(itrack, 0, 0) = trk.label();

      m_Err[iC].copyIn(itrack, trk.errors().Array());
      m_Par[iC].copyIn(itrack, trk.parameters().Array());

      m_Chg(itrack, 0, 0) = trk.charge();
      m_Chi2(itrack, 0, 0) = trk.chi2();

      // CopyIn seems fast enough, but indirections are quite slow.
      for (int hi = 0; hi < m_Nhits; ++hi) {
        const int hidx = trk.getHitIdx(hi);
        const int hlyr = trk.getHitLyr(hi);
        const Hit& hit = layerHits[hlyr].refHit(hidx);

        m_msErr[hi].copyIn(itrack, hit.errArray());
        m_msPar[hi].copyIn(itrack, hit.posArray());

        m_HoTArr[hi](itrack, 0, 0) = trk.getHitOnTrack(hi);
      }
    }
  }

inputTracksForFit()

void mkfit::MkFitter::inputTracksForFit	(	const std::vector< Track > &	tracks,
		int	beg,
		int	end
	)

Definition at line 269 of file MkFitter.cc.

References mkfit::MatriplexPackerSlurpIn< D >::addInput(), mkfit::MatriplexErrParPackerSlurpIn< T, D >::addInput(), mkfit::TrackBase::charge(), mkfit::TrackBase::chi2(), mkfit::Track::getHitsOnTrackArray(), mps_fire::i, mkfit::MkBase::iC, mkfit::TrackBase::label(), mkfit::MkBase::m_Chg, m_Chi2, mkfit::MkBase::m_Err, m_HoTArr, m_Label, mkfit::MkBase::m_Par, mkfit::Config::nLayers, mkfit::MatriplexPackerSlurpIn< D >::pack(), mkfit::MatriplexErrParPackerSlurpIn< T, D >::pack(), and pwdgSkimBPark_cfi::tracks.

                                                                                   {
    // Loads track parameters and hit indices.

    // XXXXMT4K has Config::nLayers: How many hits do we read in?
    // Check for max? Expect an argument?
    // What to do with invalid hits? Skip?

    // XXXX MT Here the same idx array WAS used for slurping in of tracks and
    // hots. With this, two index arrays are built, one within each packer.

    MatriplexTrackPacker mtp(&tracks[beg]);
    MatriplexHoTPacker mhotp(tracks[beg].getHitsOnTrackArray());

    int itrack = 0;

    for (int i = beg; i < end; ++i, ++itrack) {
      const Track& trk = tracks[i];

      m_Chg(itrack, 0, 0) = trk.charge();
      m_Chi2(itrack, 0, 0) = trk.chi2();
      m_Label(itrack, 0, 0) = trk.label();

      mtp.addInput(trk);

      mhotp.addInput(*trk.getHitsOnTrackArray());
    }

    mtp.pack(m_Err[iC], m_Par[iC]);
    for (int ll = 0; ll < Config::nLayers; ++ll) {
      mhotp.pack(m_HoTArr[ll], ll);
    }
  }

outputFittedTracks()

void mkfit::MkFitter::outputFittedTracks ( std::vector< Track > &  tracks, int  beg, int  end  ) const

inline

Definition at line 60 of file MkFitter.h.

References mkfit::MkBase::iC, outputTracks(), and pwdgSkimBPark_cfi::tracks.

                                                                              {
      return outputTracks(tracks, beg, end, iC);
    }

outputFittedTracksAndHitIdx()

void mkfit::MkFitter::outputFittedTracksAndHitIdx	(	std::vector< Track > &	tracks,
		int	beg,
		int	end,
		bool	outputProp
	)		const

Definition at line 362 of file MkFitter.cc.

References mps_fire::i, mkfit::MkBase::iC, mkfit::MkBase::iP, mkfit::MkBase::m_Chg, m_Chi2, mkfit::MkBase::m_Err, m_HoTArr, m_Label, m_Nhits, mkfit::MkBase::m_Par, and pwdgSkimBPark_cfi::tracks.

                                                                                                              {
    // Copies last track parameters (updated) into Track objects and up to m_Nhits.
    // The tracks vector should be resized to allow direct copying.

    const int iO = outputProp ? iP : iC;

    int itrack = 0;
    for (int i = beg; i < end; ++i, ++itrack) {
      m_Err[iO].copyOut(itrack, tracks[i].errors_nc().Array());
      m_Par[iO].copyOut(itrack, tracks[i].parameters_nc().Array());

      tracks[i].setCharge(m_Chg(itrack, 0, 0));
      tracks[i].setChi2(m_Chi2(itrack, 0, 0));
      tracks[i].setLabel(m_Label(itrack, 0, 0));

      // QQQQ Could do resize and std::copy, as in MkFinder::copy_out(), but
      // we do not know the correct N_found_hits.
      tracks[i].resetHits();
      tracks[i].reserveHits(m_Nhits);
      for (int hi = 0; hi < m_Nhits; ++hi) {
        tracks[i].addHitIdx(m_HoTArr[hi](itrack, 0, 0), 0.);
      }
    }
  }

outputPropagatedTracks()

void mkfit::MkFitter::outputPropagatedTracks ( std::vector< Track > &  tracks, int  beg, int  end  ) const

inline

Definition at line 64 of file MkFitter.h.

References mkfit::MkBase::iP, outputTracks(), and pwdgSkimBPark_cfi::tracks.

                                                                                  {
      return outputTracks(tracks, beg, end, iP);
    }

outputTracks()

void mkfit::MkFitter::outputTracks	(	std::vector< Track > &	tracks,
		int	beg,
		int	end,
		int	iCP
	)		const

Definition at line 345 of file MkFitter.cc.

References mps_fire::i, mkfit::MkBase::m_Chg, m_Chi2, mkfit::MkBase::m_Err, m_Label, mkfit::MkBase::m_Par, and pwdgSkimBPark_cfi::tracks.

Referenced by outputFittedTracks(), and outputPropagatedTracks().

                                                                                       {
    // Copies last track parameters (updated) into Track objects.
    // The tracks vector should be resized to allow direct copying.

    int itrack = 0;
    for (int i = beg; i < end; ++i, ++itrack) {
      m_Err[iCP].copyOut(itrack, tracks[i].errors_nc().Array());
      m_Par[iCP].copyOut(itrack, tracks[i].parameters_nc().Array());

      tracks[i].setCharge(m_Chg(itrack, 0, 0));

      // XXXXX chi2 is not set (also not in SMatrix fit, it seems)
      tracks[i].setChi2(m_Chi2(itrack, 0, 0));
      tracks[i].setLabel(m_Label(itrack, 0, 0));
    }
  }

printPt()

void mkfit::MkFitter::printPt

(

int

idx

)

Definition at line 23 of file MkFitter.cc.

References Matriplex::hypot(), mps_fire::i, heavyIonCSV_trainingSettings::idx, mkfit::MkBase::m_Par, and mkfit::NN.

                                {
    for (int i = 0; i < NN; ++i) {
      printf("%5.2f  ", std::hypot(m_Par[idx].At(i, 3, 0), m_Par[idx].At(i, 4, 0)));
    }
  }

refErr0()

MPlexLS& mkfit::MkFitter::refErr0 ( )

inline

Definition at line 21 of file MkFitter.h.

References mkfit::MkBase::m_Err.

{ return m_Err[0]; }

refPar0()

MPlexLV& mkfit::MkFitter::refPar0 ( )

inline

Definition at line 22 of file MkFitter.h.

References mkfit::MkBase::m_Par.

{ return m_Par[0]; }

setNhits()

void mkfit::MkFitter::setNhits ( int  newnhits )

inline

Definition at line 30 of file MkFitter.h.

References m_Nhits, SiStripPI::min, and mkfit::Config::nMaxTrkHits.

{ m_Nhits = std::min(newnhits, Config::nMaxTrkHits - 1); }

slurpInTracksAndHits()

void mkfit::MkFitter::slurpInTracksAndHits	(	const std::vector< Track > &	tracks,
		const std::vector< HitVec > &	layerHits,
		int	beg,
		int	end
	)

Definition at line 123 of file MkFitter.cc.

References mkfit::MatriplexErrParPackerSlurpIn< T, D >::addInput(), mkfit::TrackBase::charge(), mkfit::TrackBase::chi2(), data, mps_fire::i, mkfit::MkBase::iC, mkfit::TrackBase::label(), mkfit::MkBase::m_Chg, m_Chi2, mkfit::MkBase::m_Err, m_HoTArr, m_Label, m_msErr, m_msPar, m_Nhits, mkfit::MkBase::m_Par, mkfit::MatriplexErrParPackerSlurpIn< T, D >::pack(), and pwdgSkimBPark_cfi::tracks.

                                               {
    // Assign track parameters to initial state and copy hit values in.

    // This might not be true for the last chunk!
    // assert(end - beg == NN);

    MatriplexTrackPacker mtp(&tracks[beg]);

    for (int i = beg; i < end; ++i) {
      int itrack = i - beg;
      const Track& trk = tracks[i];

      m_Label(itrack, 0, 0) = trk.label();

      mtp.addInput(trk);

      m_Chg(itrack, 0, 0) = trk.charge();
      m_Chi2(itrack, 0, 0) = trk.chi2();
    }

    mtp.pack(m_Err[iC], m_Par[iC]);

    // CopyIn seems fast enough, but indirections are quite slow.
    for (int hi = 0; hi < m_Nhits; ++hi) {
      MatriplexHitPacker mhp(layerHits[hi].data());

      for (int i = beg; i < end; ++i) {
        const int hidx = tracks[i].getHitIdx(hi);
        const Hit& hit = layerHits[hi][hidx];

        m_HoTArr[hi](i - beg, 0, 0) = tracks[i].getHitOnTrack(hi);

        mhp.addInput(hit);
      }

      mhp.pack(m_msErr[hi], m_msPar[hi]);
    }
  }

Member Data Documentation

m_CandIdx

MPlexQI mkfit::MkFitter::m_CandIdx

private

Definition at line 80 of file MkFitter.h.

Referenced by inputTracksAndHitIdx().

m_Chi2

MPlexQF mkfit::MkFitter::m_Chi2

private

Definition at line 73 of file MkFitter.h.

Referenced by inputSeedsTracksAndHits(), inputTracksAndHitIdx(), inputTracksAndHits(), inputTracksForFit(), outputFittedTracksAndHitIdx(), outputTracks(), and slurpInTracksAndHits().

m_HoTArr

MPlexQHoT mkfit::MkFitter::m_HoTArr[Config::nMaxTrkHits]

private

Definition at line 82 of file MkFitter.h.

Referenced by countInvalidHits(), countValidHits(), fitTracksWithInterSlurp(), inputSeedsTracksAndHits(), inputTracksAndHitIdx(), inputTracksAndHits(), inputTracksForFit(), outputFittedTracksAndHitIdx(), and slurpInTracksAndHits().

m_Label

MPlexQI mkfit::MkFitter::m_Label

private

Definition at line 78 of file MkFitter.h.

Referenced by inputSeedsTracksAndHits(), inputTracksAndHitIdx(), inputTracksAndHits(), inputTracksForFit(), outputFittedTracksAndHitIdx(), outputTracks(), and slurpInTracksAndHits().

m_msErr

MPlexHS mkfit::MkFitter::m_msErr[Config::nMaxTrkHits]

private

Definition at line 75 of file MkFitter.h.

Referenced by checkAlignment(), fitTracksWithInterSlurp(), inputSeedsTracksAndHits(), inputTracksAndHits(), and slurpInTracksAndHits().

m_msPar

MPlexHV mkfit::MkFitter::m_msPar[Config::nMaxTrkHits]

private

Definition at line 76 of file MkFitter.h.

Referenced by checkAlignment(), fitTracksWithInterSlurp(), inputSeedsTracksAndHits(), inputTracksAndHits(), and slurpInTracksAndHits().

m_Nhits

int mkfit::MkFitter::m_Nhits

private

Definition at line 88 of file MkFitter.h.

Referenced by countInvalidHits(), countValidHits(), fitTracksWithInterSlurp(), inputSeedsTracksAndHits(), inputTracksAndHitIdx(), inputTracksAndHits(), outputFittedTracksAndHitIdx(), setNhits(), and slurpInTracksAndHits().

m_SeedIdx

MPlexQI mkfit::MkFitter::m_SeedIdx

private

Definition at line 79 of file MkFitter.h.

Referenced by inputTracksAndHitIdx().

m_XHitArr

MPlexHitIdx mkfit::MkFitter::m_XHitArr

private

Definition at line 86 of file MkFitter.h.

m_XHitSize

MPlexQI mkfit::MkFitter::m_XHitSize

private

Definition at line 85 of file MkFitter.h.