#include <Traj2TrackHits.h>

Public Types
using	Chi2sFive = std::vector< unsigned char >

using	TrajParams = std::vector< LocalTrajectoryParameters >

Public Member Functions
void	operator() (Trajectory const &traj, TrackingRecHitCollection &hits, bool splitting) const

void	operator() (Trajectory const &traj, TrackingRecHitCollection &hits, TrajParams &trajParams, Chi2sFive &chi2s) const

	Traj2TrackHits ()

	Traj2TrackHits (const TransientTrackingRecHitBuilder *builder, bool ikeepOrder, bool noNoDet=true)

Static Public Member Functions
static unsigned char	toChi2x5 (float chi2)

Private Member Functions
BaseTrackerRecHit *	build (GeomDetUnit const &idet, OmniClusterRef const &clus) const

TrackingRecHit *	clone (BaseTrackerRecHit const &hit2D) const

template<typename HI >
void	copy (HI itm, HI e, TrackingRecHitCollection &hits, TrajParams &trajParams, Chi2sFive &chi2s) const

template<typename HI >
void	copy (HI itm, HI e, TrackingRecHitCollection &hits) const

template<typename HI >
void	split (HI itm, HI e, TrackingRecHitCollection &hits, bool along) const

void	split (TrajectoryMeasurement const &itm, SiStripMatchedRecHit2D const &mhit, TrackingRecHitCollection &hits, float zdir) const

Private Attributes
bool	keepOrder = false

bool	removeNoDet = true

const StripClusterParameterEstimator *	theCPE = 0

Detailed Description

Definition at line 16 of file Traj2TrackHits.h.

Member Typedef Documentation

using Traj2TrackHits::Chi2sFive = std::vector<unsigned char>

Definition at line 23 of file Traj2TrackHits.h.

using Traj2TrackHits::TrajParams = std::vector<LocalTrajectoryParameters>

Definition at line 22 of file Traj2TrackHits.h.

Constructor & Destructor Documentation

Traj2TrackHits::Traj2TrackHits ( )

inline

Definition at line 31 of file Traj2TrackHits.h.

31 {}

Traj2TrackHits::Traj2TrackHits	(	const TransientTrackingRecHitBuilder *	builder,
		bool	ikeepOrder,
		bool	noNoDet = `true`
	)

inline

Definition at line 33 of file Traj2TrackHits.h.

       : theCPE(static_cast<TkTransientTrackingRecHitBuilder const *>(builder)->stripClusterParameterEstimator()),
         keepOrder(ikeepOrder),
         removeNoDet(noNoDet) {}

Member Function Documentation

BaseTrackerRecHit* Traj2TrackHits::build	(	GeomDetUnit const &	idet,
		OmniClusterRef const &	clus
	)		const

inlineprivate

Definition at line 131 of file Traj2TrackHits.h.

References makeMuonMisalignmentScenario::endcap, GeomDetType::isEndcap(), StripClusterParameterEstimator::localParameters(), SiStripPI::max, OmniClusterRef::stripCluster(), and GeomDet::type().

Referenced by split().

                                                                                       {
     //Use 2D SiStripRecHit in endcap
     bool endcap = idet.type().isEndcap();
     auto &&lv = theCPE->localParameters(clus.stripCluster(), idet);
     if (endcap)
       return new SiStripRecHit2D(lv.first, lv.second, idet, clus);
     return new SiStripRecHit1D(
         lv.first, LocalError(lv.second.xx(), 0.f, std::numeric_limits<float>::max()), idet, clus);
   }

TrackingRecHit* Traj2TrackHits::clone ( BaseTrackerRecHit const & hit2D ) const

inlineprivate

Definition at line 119 of file Traj2TrackHits.h.

References TrackingRecHit::clone(), TrackingRecHit::det(), TrackingRecHit::detUnit(), makeMuonMisalignmentScenario::endcap, BaseTrackerRecHit::firstClusterRef(), BaseTrackerRecHit::localPosition(), BaseTrackerRecHit::localPositionError(), SiStripPI::max, and LocalError::xx().

Referenced by split().

                                                               {
     auto const &detU = *hit2D.detUnit();
     //Use 2D SiStripRecHit in endcap
     bool endcap = detU.type().isEndcap();
     if (endcap)
       return hit2D.clone();
     return new SiStripRecHit1D(hit2D.localPosition(),
                                LocalError(hit2D.localPositionError().xx(), 0.f, std::numeric_limits<float>::max()),
                                *hit2D.det(),
                                hit2D.firstClusterRef());
   }

template<typename HI >

void Traj2TrackHits::copy	(	HI	itm,
		HI	e,
		TrackingRecHitCollection &	hits,
		TrajParams &	trajParams,
		Chi2sFive &	chi2s
	)		const

inlineprivate

Definition at line 74 of file Traj2TrackHits.h.

References MillePedeFileConverter_cfg::e, edm::OwnVector< T, P >::push_back(), and toChi2x5().

Referenced by operator()().

                                                                                                           {
     for (; itm != e; ++itm)
       if ((!removeNoDet) | ((*itm).recHitR().det() != nullptr)) {
         hits.push_back((*itm).recHitR().cloneHit());
         trajParams.push_back((*itm).updatedState().localParameters());
         chi2s.push_back(toChi2x5((*itm).estimate()));
       }
   }

template<typename HI >

void Traj2TrackHits::copy	(	HI	itm,
		HI	e,
		TrackingRecHitCollection &	hits
	)		const

inlineprivate

Definition at line 84 of file Traj2TrackHits.h.

References MillePedeFileConverter_cfg::e, and edm::OwnVector< T, P >::push_back().

                                                                 {
     for (; itm != e; ++itm)
       if ((!removeNoDet) | ((*itm).recHitR().det() != nullptr))
         hits.push_back((*itm).recHitR().cloneHit());
   }

void Traj2TrackHits::operator()	(	Trajectory const &	traj,
		TrackingRecHitCollection &	hits,
		bool	splitting
	)		const

inline

Definition at line 38 of file Traj2TrackHits.h.

References alongMomentum, copy(), Trajectory::direction(), hfClusterShapes_cfi::hits, Trajectory::measurements(), edm::OwnVector< T, P >::reserve(), edm::OwnVector< T, P >::shrink_to_fit(), and split().

                                                                                                 {
     // ---  NOTA BENE: the convention is to sort hits and measurements "along the momentum".
     bool along = traj.direction() == alongMomentum;
     auto const &meas = traj.measurements();
     hits.reserve(splitting ? 2 * meas.size() : meas.size());
     if (!splitting) {
       if (keepOrder | along)
         copy(meas.begin(), meas.end(), hits);
       else
         copy(meas.rbegin(), meas.rend(), hits);
       return;
     }
     if (keepOrder | along)
       split(meas.begin(), meas.end(), hits, along);
     else
       split(meas.rbegin(), meas.rend(), hits, along);
     hits.shrink_to_fit();
   }

void Traj2TrackHits::operator()	(	Trajectory const &	traj,
		TrackingRecHitCollection &	hits,
		TrajParams &	trajParams,
		Chi2sFive &	chi2s
	)		const

inline

Definition at line 57 of file Traj2TrackHits.h.

References alongMomentum, copy(), Trajectory::direction(), hfClusterShapes_cfi::hits, and Trajectory::measurements().

                                           {
     // ---  NOTA BENE: the convention is to sort hits and measurements "along the momentum".
     bool along = traj.direction() == alongMomentum;
     auto const &meas = traj.measurements();
     trajParams.reserve(meas.size());
     chi2s.reserve(meas.size());
     if (keepOrder | along)
       copy(meas.begin(), meas.end(), hits, trajParams, chi2s);
     else
       copy(meas.rbegin(), meas.rend(), hits, trajParams, chi2s);
   }

template<typename HI >

void Traj2TrackHits::split	(	HI	itm,
		HI	e,
		TrackingRecHitCollection &	hits,
		bool	along
	)		const

inlineprivate

Definition at line 91 of file Traj2TrackHits.h.

References build(), clone(), MillePedeFileConverter_cfg::e, BaseTrackerRecHit::firstClusterRef(), trackerHitRTTI::isNotFromCluster(), trackerHitRTTI::isUndef(), and edm::OwnVector< T, P >::push_back().

Referenced by operator()().

                                                                              {
     for (; itm != e; ++itm) {
       auto const &hit = *(*itm).recHit()->hit();
       if ((removeNoDet) & ((*itm).recHitR().det() == nullptr))
         continue;
       if (trackerHitRTTI::isUndef(hit) | trackerHitRTTI::isNotFromCluster(hit) | (hit.dimension() != 2)) {
         hits.push_back(hit.clone());
         continue;
       }
       auto const &thit = static_cast<BaseTrackerRecHit const &>(hit);
       auto const &clus = thit.firstClusterRef();
       if (clus.isPixel())
         hits.push_back(hit.clone());
       else if (clus.isPhase2())
         hits.push_back(hit.clone());
       else if (thit.isMatched()) {
         auto zdir = itm->updatedState().localDirection().z();
         if (keepOrder & (!along))
           zdir = -zdir;
         split(*itm, static_cast<SiStripMatchedRecHit2D const &>(thit), hits, zdir);
       } else if (thit.isProjected()) {
         auto detU = static_cast<ProjectedSiStripRecHit2D const &>(thit).originalDet();
         hits.push_back(build(*detU, clus));
       } else
         hits.push_back(clone(thit));
     }
   }

void Traj2TrackHits::split	(	TrajectoryMeasurement const &	itm,
		SiStripMatchedRecHit2D const &	mhit,
		TrackingRecHitCollection &	hits,
		float	zdir
	)		const

inlineprivate

Definition at line 141 of file Traj2TrackHits.h.

References build(), TrackingRecHit::det(), SiStripMatchedRecHit2D::monoClusterRef(), GluedGeomDet::monoDet(), GeomDet::position(), edm::OwnVector< T, P >::push_back(), hitfit::scalar(), SiStripMatchedRecHit2D::stereoClusterRef(), GluedGeomDet::stereoDet(), TrajectoryStateOnSurface::surface(), GloballyPositioned< T >::toLocal(), TrajectoryMeasurement::updatedState(), and PV3DBase< T, PVType, FrameType >::z().

                                {
     const GluedGeomDet *gdet = static_cast<const GluedGeomDet *>(mhit.det());
 
     auto hitM = build(*gdet->monoDet(), mhit.monoClusterRef());
     auto hitS = build(*gdet->stereoDet(), mhit.stereoClusterRef());
 
     // we should find a faster way
     LocalPoint firstLocalPos = itm.updatedState().surface().toLocal(gdet->monoDet()->position());
     LocalPoint secondLocalPos = itm.updatedState().surface().toLocal(gdet->stereoDet()->position());
     LocalVector Delta = secondLocalPos - firstLocalPos;
     float scalar = Delta.z() * zdir;
     // hit along the direction
     if (scalar < 0) {
       hits.push_back(hitS);
       hits.push_back(hitM);
     } else {
       hits.push_back(hitM);
       hits.push_back(hitS);
     }
   }

static unsigned char Traj2TrackHits::toChi2x5 ( float chi2 )

inlinestatic

Definition at line 25 of file Traj2TrackHits.h.

References f, and min().

Referenced by copy(), and HITrackClusterRemover::process().

                                             {
     float tc = std::round(5.f * chi2);
     return std::min(tc, 255.f);
   }

Member Data Documentation

bool Traj2TrackHits::keepOrder = false

private

Definition at line 19 of file Traj2TrackHits.h.

bool Traj2TrackHits::removeNoDet = true

private

Definition at line 20 of file Traj2TrackHits.h.

const StripClusterParameterEstimator* Traj2TrackHits::theCPE = 0

private

Definition at line 18 of file Traj2TrackHits.h.

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation