127 auto trajectoryCollection = std::make_unique<std::vector<Trajectory>>();
134 for (reco::TrackCollection::const_iterator
track = tracks->begin();
track != tracks->end(); ++
track) {
140 std::vector<TrajectoryMeasurement::ConstRecHitPointer> transHits;
141 std::vector<TrajectoryStateOnSurface> TSOSes;
144 DetId id = (*hit)->geographicalId();
160 transHits.push_back(hitPtr);
161 TSOSes.push_back(
TrajectoryStateOnSurface(localTrajectoryParameters, localTrajectoryError, layerSurface, &*magneticField));
165 assert(clonedHits.
size() == transHits.size());
166 assert(transHits.size() == TSOSes.size());
169 if (!clonedHits.
empty()) {
175 std::vector<TrajectoryMeasurement::ConstRecHitPointer>::const_iterator transHitPtr = transHits.begin();
176 std::vector<TrajectoryStateOnSurface>::const_iterator
TSOS = TSOSes.begin();
177 for (; clonedHit != clonedHits.
end(); ++clonedHit, ++transHitPtr, ++
TSOS) {
181 trajectoryCollection->push_back(trajectory);
185 reference_map[trajCounter] = trackCounter;
190 unsigned int numTrajectories = trajectoryCollection->size();
196 auto trajTrackMap = std::make_unique<TrajTrackAssociationCollection>();
198 for (trajCounter = 0; trajCounter < numTrajectories; trajCounter++) {
RecHitPointer build(const TrackingRecHit *p, edm::ESHandle< GlobalTrackingGeometry > trackingGeometry) const
Call the MuonTransientTrackingRecHit::specificBuild.
std::remove_cv< typename std::remove_reference< argument_type >::type >::type key_type
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
const GeomDet * idToDet(DetId) const override
const Plane & surface() const
The nominal surface of the GeomDet.
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
TrajectoryStateOnSurface TSOS
ROOT::Math::SVector< double, 5 > AlgebraicVector5
TrackingRecHit::ConstRecHitPointer ConstRecHitPointer
TrackingRecHitCollection::base::const_iterator trackingRecHit_iterator
iterator over a vector of reference to TrackingRecHit in the same collection