87 std::vector<TrajectorySeed> &
out,
97 layersToTry_(iConfig.getParameter<int32_t>(
"layersToTry")),
98 hitsToTry_(iConfig.getParameter<int32_t>(
"hitsToTry")),
109 produces<std::vector<TrajectorySeed>>();
132 auto out = std::make_unique<std::vector<TrajectorySeed>>();
134 for (
auto const &
mu : *src) {
137 if (
debug_ &&
mu.innerTrack().isNonnull())
143 std::unique_ptr<Propagator> pmuon_cloned =
147 int sizeBefore =
out->size();
149 LogDebug(
"OutsideInMuonSeeder") <<
"\n\n\nSeeding for muon of pt " <<
mu.pt() <<
", eta " <<
mu.eta() <<
", phi " 150 <<
mu.phi() << std::endl;
158 std::vector<BarrelDetLayer const *>
const &tob = measurementTracker->geometricSearchTracker()->tobLayers();
160 int iLayer = tob.size();
162 LogError(
"OutsideInMuonSeeder") <<
"TOB has no layers.";
164 for (
auto it = tob.rbegin(), ed = tob.rend(); it != ed; ++it, --iLayer) {
166 LogDebug(
"OutsideInMuonSeeder") <<
"\n ==== Trying TOB " << iLayer <<
" ====" << std::endl;
167 if (
doLayer(**it, state, *
out, *(pmuon_cloned.get()), *(ptracker_cloned.get()), *measurementTracker)) {
175 ? measurementTracker->geometricSearchTracker()->posTidLayers()
176 : measurementTracker->geometricSearchTracker()->posTecLayers();
178 LogDebug(
"OutsideInMuonSeeder") <<
"\n We are using the Phase2 Outer Tracker (defined as a TID+). ";
180 LogTrace(
"OutsideInMuonSeeder") <<
"\n ==== TEC+ tot layers " << forwLayers.size() <<
" ====" << std::endl;
182 int iLayer = forwLayers.size();
184 LogError(
"OutsideInMuonSeeder") <<
"TEC+ has no layers.";
187 LogDebug(
"OutsideInMuonSeeder") <<
"\n ==== Tot layers " << forwLayers.size() <<
" ====" << std::endl;
188 for (
auto it = forwLayers.rbegin(), ed = forwLayers.rend(); it != ed; ++it, --iLayer) {
190 LogDebug(
"OutsideInMuonSeeder") <<
"\n ==== Trying Forward Layer +" << +iLayer <<
" ====" << std::endl;
191 if (
doLayer(**it, state, *
out, *(pmuon_cloned.get()), *(ptracker_cloned.get()), *measurementTracker)) {
199 ? measurementTracker->geometricSearchTracker()->negTidLayers()
200 : measurementTracker->geometricSearchTracker()->negTecLayers();
202 LogDebug(
"OutsideInMuonSeeder") <<
"\n We are using the Phase2 Outer Tracker (defined as a TID-). ";
204 LogTrace(
"OutsideInMuonSeeder") <<
"\n ==== TEC- tot layers " << forwLayers.size() <<
" ====" << std::endl;
206 int iLayer = forwLayers.size();
208 LogError(
"OutsideInMuonSeeder") <<
"TEC- has no layers.";
211 LogDebug(
"OutsideInMuonSeeder") <<
"\n ==== Tot layers " << forwLayers.size() <<
" ====" << std::endl;
212 for (
auto it = forwLayers.rbegin(), ed = forwLayers.rend(); it != ed; ++it, --iLayer) {
214 LogDebug(
"OutsideInMuonSeeder") <<
"\n ==== Trying Forward Layer -" << -iLayer <<
" ====" << std::endl;
215 if (
doLayer(**it, state, *
out, *(pmuon_cloned.get()), *(ptracker_cloned.get()), *measurementTracker)) {
222 LogDebug(
"OutsideInMuonSeeder") <<
"Outcome of seeding for muon of pt " <<
mu.pt() <<
", eta " <<
mu.eta()
223 <<
", phi " <<
mu.phi() <<
": found " << (
out->size() - sizeBefore) <<
" seeds." 232 std::vector<TrajectorySeed> &
out,
238 std::vector<GeometricSearchDet::DetWithState> dets;
245 <<
" returned " << dets.size() <<
" compatible detectors" << std::endl;
248 std::vector<TrajectoryMeasurement> meas;
249 for (std::vector<GeometricSearchDet::DetWithState>::const_iterator it = dets.begin(), ed = dets.end(); it != ed;
253 std::cerr <<
"BOGUS detid " << it->first->geographicalId().rawId() << std::endl;
256 if (!it->second.isValid())
258 std::vector<TrajectoryMeasurement> mymeas =
261 LogDebug(
"OutsideInMuonSeeder") <<
"Query on detector " << it->first->geographicalId().rawId() <<
" returned " 262 << mymeas.size() <<
" measurements." << std::endl;
263 for (std::vector<TrajectoryMeasurement>::const_iterator it2 = mymeas.begin(), ed2 = mymeas.end(); it2 != ed2;
265 if (it2->recHit()->isValid())
266 meas.push_back(*it2);
271 for (std::vector<TrajectoryMeasurement>::const_iterator it2 = meas.begin(), ed2 = meas.end(); it2 != ed2; ++it2) {
273 LogDebug(
"OutsideInMuonSeeder") <<
" inspecting Hit with chi2 = " << it2->estimate() << std::endl;
274 LogDebug(
"OutsideInMuonSeeder") <<
" track state " << it2->forwardPredictedState().globalPosition()
276 LogDebug(
"OutsideInMuonSeeder") <<
" rechit position " << it2->recHit()->globalPosition() << std::endl;
284 seedHits.
push_back(*it2->recHit()->hit());
306 tsos = pmuon_cloned->propagate(tsos, det->
surface());
void rescaleError(double factor)
MeasurementDetWithData idToDet(const DetId &id) const
Previous MeasurementDetSystem interface.
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
const std::string estimatorName_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
size_t recHitsSize() const
Get number of RecHits. (Warning, this includes invalid hits, which are not physical hits)...
edm::ESHandle< Chi2MeasurementEstimatorBase > estimator_
edm::ESHandle< MagneticField > magfield_
constexpr uint32_t rawId() const
get the raw id
edm::ESHandle< GlobalTrackingGeometry > geometry_
float const maxEtaForTOB_
GlobalPoint globalPosition() const
edm::ESHandle< Propagator > trackerPropagator_
const bool fromVertex_
Do inside-out.
edm::ESHandle< TrajectoryStateUpdator > updator_
const Plane & surface() const
The nominal surface of the GeomDet.
LocalError positionError() const
int doLayer(const GeometricSearchDet &layer, const TrajectoryStateOnSurface &state, std::vector< TrajectorySeed > &out, const Propagator &muon_propagator, const Propagator &tracker_propagator, const MeasurementTrackerEvent &mte) const
virtual TrajectoryStateOnSurface update(const TrajectoryStateOnSurface &, const TrackingRecHit &) const =0
bool isThere(GeomDetEnumerators::SubDetector subdet) const
std::unique_ptr< Propagator > SetPropagationDirection(Propagator const &iprop, PropagationDirection dir)
const bool debug_
Dump deug information.
#define DEFINE_FWK_MODULE(type)
double eta() const
pseudorapidity of momentum vector
const double errorRescaling_
How much to rescale errors from STA.
std::vector< TrajectoryMeasurement > fastMeasurements(const TrajectoryStateOnSurface &stateOnThisDet, const TrajectoryStateOnSurface &tsos2, const Propagator &prop, const MeasurementEstimator &est) const
Abs< T >::type abs(const T &t)
const std::string trackerPropagatorName_
const std::string updatorName_
void produce(edm::Event &iEvent, const edm::EventSetup &iSetup) override
const LocalTrajectoryError & localError() const
edm::ESHandle< Propagator > muonPropagator_
StringCutObjectSelector< reco::Muon > selector_
Muon selection.
edm::EDGetTokenT< edm::View< reco::Muon > > src_
Labels for input collections.
void doDebug(const reco::Track &tk) const
const int layersToTry_
How many layers to try.
GlobalVector globalMomentum() const
TrackingRecHitRef recHit(size_t i) const
Get i-th hit on the track.
const GeomDet * idToDet(DetId) const override
const std::string measurementTrackerName_
const int hitsToTry_
How many hits to try on same layer.
edm::EDGetTokenT< MeasurementTrackerEvent > measurementTrackerTag_
DetId geographicalId() const
float const minEtaForTEC_
T const * product() const
const std::string muonPropagatorName_
virtual void compatibleDetsV(const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetWithState > &result) const
OutsideInMuonSeeder(const edm::ParameterSet &iConfig)
~OutsideInMuonSeeder() override