40 std::vector<const TrackingRecHit *> hits;
44 hits.push_back(&**it);
46 DetId id(hits.back()->geographicalId());
47 std::cout <<
" subdet " <<
id.subdetId() <<
" layer " <<
theTrkTopo->layer(
id) <<
" valid " << hits.back()->
isValid() <<
" detid: " << id() << std::endl;
53 DetId lastId(hits.back()->geographicalId());
55 unsigned int lastLayer =
theTrkTopo->layer(lastId);
60 if (thisSubdet > lastSubdet || (thisSubdet == lastSubdet &&
theTrkTopo->layer(
id) > lastLayer)) {
62 if (
debug_)
std::cout <<
" adding subdet " <<
id.subdetId() <<
" layer " <<
theTrkTopo->layer(
id) <<
" valid " << hit->
isValid() <<
" detid: " << id() << std::endl;
64 if (
debug_)
std::cout <<
" skipping subdet " << thisSubdet <<
" layer " <<
theTrkTopo->layer(
id) <<
" valid " << hit->
isValid() <<
" detid: " << id() << std::endl;
68 size_t nHitsFirstTrack = hits.size();
75 if (
debug_)
std::cout <<
" subdet " <<
id.subdetId() <<
" layer " <<
theTrkTopo->layer(
id) <<
" valid " << valid <<
" detid: " << id() << std::endl;
78 ++iHit;
if (iHit > nHitsFirstTrack)
break;
80 if (
id.subdetId() != id2.
subdetId())
continue;
83 if (
debug_)
std::cout <<
" discared as duplicate of other hit" << id() << std::endl;
87 if (
debug_)
std::cout <<
" replacing old invalid hit on detid " << id2() << std::endl;
88 hit2 = hit; shared =
true;
break;
90 if (
debug_)
std::cout <<
" discared as additional hit on layer that already contains hit with detid " << id() << std::endl;
103 unsigned int nhits = hits.
size();
120 std::vector<TransientTrackingRecHit::RecHitPointer> ttrh(nhits);
121 for (
unsigned int i = 0;
i < nhits; ++
i) ttrh[
i] =
theBuilder->build(hits[
i]);
122 std::sort(ttrh.begin(), ttrh.end(), GlobalMomentumSort(
v));
edm::ESHandle< TrackerTopology > theTrkTopo
size_t recHitsSize() const
Get number of RecHits. (Warning, this includes invalid hits, which are not physical hits)...
virtual bool sharesInput(const TrackingRecHit *other, SharedInputType what) const
const Vector & momentum() const
track momentum vector
const math::XYZPoint & outerPosition() const
position of the outermost hit
const math::XYZPoint & innerPosition() const
position of the innermost hit
unsigned int outerDetId() const
DetId of the detector on which surface the outermost state is located.
trackingRecHit_iterator recHitsBegin() const
Iterator to first hit on the track.
int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum) ...
std::shared_ptr< TrackingRecHit const > RecHitPointer
const math::XYZVector & outerMomentum() const
momentum vector at the outermost hit position
XYZVectorD XYZVector
spatial vector with cartesian internal representation
edm::ESHandle< TransientTrackingRecHitBuilder > theBuilder
const math::XYZVector & innerMomentum() const
momentum vector at the innermost hit position
DetId geographicalId() const
unsigned int innerDetId() const
DetId of the detector on which surface the innermost state is located.
edm::ESHandle< MagneticField > theMagField
edm::ESHandle< TrackerGeometry > theGeometry
TrackingRecHitCollection::base::const_iterator trackingRecHit_iterator
iterator over a vector of reference to TrackingRecHit in the same collection
trackingRecHit_iterator recHitsEnd() const
Iterator to last hit on the track.