Traj2TrackHits Class Reference

#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 18 of file Traj2TrackHits.h.

Member Typedef Documentation

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

Definition at line 25 of file Traj2TrackHits.h.

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

Definition at line 24 of file Traj2TrackHits.h.

Constructor & Destructor Documentation

Traj2TrackHits::Traj2TrackHits ( )

inline

Definition at line 30 of file Traj2TrackHits.h.

{}

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

inline

Definition at line 32 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 113 of file Traj2TrackHits.h.

References makeMuonMisalignmentScenario::endcap, GeomDetType::isEndcap(), StripClusterParameterEstimator::localParameters(), hpstanc_transforms::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 101 of file Traj2TrackHits.h.

References TrackingRecHit::clone(), TrackingRecHit::det(), TrackingRecHit::detUnit(), makeMuonMisalignmentScenario::endcap, BaseTrackerRecHit::firstClusterRef(), BaseTrackerRecHit::localPosition(), BaseTrackerRecHit::localPositionError(), hpstanc_transforms::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 63 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 72 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 36 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 51 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 78 of file Traj2TrackHits.h.

References build(), clone(), MillePedeFileConverter_cfg::e, BaseTrackerRecHit::firstClusterRef(), 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) | ( 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 122 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 27 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 21 of file Traj2TrackHits.h.

bool Traj2TrackHits::removeNoDet =true

private

Definition at line 22 of file Traj2TrackHits.h.

const StripClusterParameterEstimator* Traj2TrackHits::theCPE = 0

private

Definition at line 20 of file Traj2TrackHits.h.