Inheritance diagram for OutsideInMuonSeeder:

Public Member Functions
	OutsideInMuonSeeder (const edm::ParameterSet &iConfig)

void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override

	~OutsideInMuonSeeder () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
void	doDebug (const reco::Track &tk) 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

Private Attributes
const bool	debug_
	Dump deug information. More...

const double	errorRescaling_
	How much to rescale errors from STA. More...

edm::ESHandle< Chi2MeasurementEstimatorBase >	estimator_

const std::string	estimatorName_

const bool	fromVertex_
	Do inside-out. More...

edm::ESHandle< GlobalTrackingGeometry >	geometry_

const int	hitsToTry_
	How many hits to try on same layer. More...

const int	layersToTry_
	How many layers to try. More...

edm::ESHandle< MagneticField >	magfield_

const float	maxEtaForTOB_

const std::string	measurementTrackerName_

edm::EDGetTokenT< MeasurementTrackerEvent >	measurementTrackerTag_

const float	minEtaForTEC_

edm::ESHandle< Propagator >	muonPropagator_

const std::string	muonPropagatorName_

StringCutObjectSelector< reco::Muon >	selector_
	Muon selection. More...

edm::EDGetTokenT< edm::View< reco::Muon > >	src_
	Labels for input collections. More...

edm::ESHandle< Propagator >	trackerPropagator_

const std::string	trackerPropagatorName_

edm::ESHandle< TrajectoryStateUpdator >	updator_

const std::string	updatorName_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 39 of file OutsideInMuonSeeder.cc.

Constructor & Destructor Documentation

◆ OutsideInMuonSeeder()

OutsideInMuonSeeder::OutsideInMuonSeeder ( const edm::ParameterSet & iConfig )

explicit

Definition at line 94 of file OutsideInMuonSeeder.cc.

     : src_(consumes<edm::View<reco::Muon>>(iConfig.getParameter<edm::InputTag>("src"))),
       selector_(iConfig.existsAs<std::string>("cut") ? iConfig.getParameter<std::string>("cut") : "", true),
       layersToTry_(iConfig.getParameter<int32_t>("layersToTry")),
       hitsToTry_(iConfig.getParameter<int32_t>("hitsToTry")),
       fromVertex_(iConfig.getParameter<bool>("fromVertex")),
       errorRescaling_(iConfig.getParameter<double>("errorRescaleFactor")),
       trackerPropagatorName_(iConfig.getParameter<std::string>("trackerPropagator")),
       muonPropagatorName_(iConfig.getParameter<std::string>("muonPropagator")),
       measurementTrackerTag_(consumes<MeasurementTrackerEvent>(edm::InputTag("MeasurementTrackerEvent"))),
       estimatorName_(iConfig.getParameter<std::string>("hitCollector")),
       updatorName_("KFUpdator"),
       minEtaForTEC_(iConfig.getParameter<double>("minEtaForTEC")),
       maxEtaForTOB_(iConfig.getParameter<double>("maxEtaForTOB")),
       debug_(iConfig.getUntrackedParameter<bool>("debug", false)) {
   produces<std::vector<TrajectorySeed>>();
 }

◆ ~OutsideInMuonSeeder()

OutsideInMuonSeeder::~OutsideInMuonSeeder ( )

inlineoverride

Definition at line 42 of file OutsideInMuonSeeder.cc.

42 {}

Member Function Documentation

◆ doDebug()

void OutsideInMuonSeeder::doDebug ( const reco::Track & tk ) const

private

Definition at line 297 of file OutsideInMuonSeeder.cc.

                                                            {
   TrajectoryStateOnSurface tsos = trajectoryStateTransform::innerStateOnSurface(tk, *geometry_, &*magfield_);
   std::unique_ptr<Propagator> pmuon_cloned = SetPropagationDirection(*muonPropagator_, alongMomentum);
   for (unsigned int i = 0; i < tk.recHitsSize(); ++i) {
     const TrackingRecHit *hit = &*tk.recHit(i);
     const GeomDet *det = geometry_->idToDet(hit->geographicalId());
     if (det == nullptr)
       continue;
     if (i != 0)
       tsos = pmuon_cloned->propagate(tsos, det->surface());
     if (!tsos.isValid())
       continue;
     LogDebug("OutsideInMuonSeeder") << "  state " << i << " at x = " << tsos.globalPosition()
                                     << ", p = " << tsos.globalMomentum() << std::endl;
     if (hit->isValid()) {
       LogDebug("OutsideInMuonSeeder") << "         valid   rechit on detid " << hit->geographicalId().rawId()
                                       << std::endl;
     } else {
       LogDebug("OutsideInMuonSeeder") << "         invalid rechit on detid " << hit->geographicalId().rawId()
                                       << std::endl;
     }
   }
 }

References alongMomentum, geometry_, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), mps_fire::i, GlobalTrackingGeometry::idToDet(), trajectoryStateTransform::innerStateOnSurface(), TrajectoryStateOnSurface::isValid(), LogDebug, magfield_, muonPropagator_, reco::Track::recHit(), reco::Track::recHitsSize(), SetPropagationDirection(), and GeomDet::surface().

Referenced by produce().

◆ doLayer()

int OutsideInMuonSeeder::doLayer	(	const GeometricSearchDet &	layer,
		const TrajectoryStateOnSurface &	state,
		std::vector< TrajectorySeed > &	out,
		const Propagator &	muon_propagator,
		const Propagator &	tracker_propagator,
		const MeasurementTrackerEvent &	mte
	)		const

private

Definition at line 230 of file OutsideInMuonSeeder.cc.

                                                                                           {
   TrajectoryStateOnSurface onLayer(state);
   onLayer.rescaleError(errorRescaling_);
   std::vector<GeometricSearchDet::DetWithState> dets;
   layer.compatibleDetsV(onLayer, muon_propagator, *estimator_, dets);
  
   if (debug_) {
     LogDebug("OutsideInMuonSeeder") << "Query on layer around x = " << onLayer.globalPosition()
                                     << " with local pos error " << sqrt(onLayer.localError().positionError().xx())
                                     << " ,  " << sqrt(onLayer.localError().positionError().yy()) << " ,  "
                                     << " returned " << dets.size() << " compatible detectors" << std::endl;
   }
  
   std::vector<TrajectoryMeasurement> meas;
   for (std::vector<GeometricSearchDet::DetWithState>::const_iterator it = dets.begin(), ed = dets.end(); it != ed;
        ++it) {
     MeasurementDetWithData det = measurementTracker.idToDet(it->first->geographicalId());
     if (det.isNull()) {
       std::cerr << "BOGUS detid " << it->first->geographicalId().rawId() << std::endl;
       continue;
     }
     if (!it->second.isValid())
       continue;
     std::vector<TrajectoryMeasurement> mymeas =
         det.fastMeasurements(it->second, state, tracker_propagator, *estimator_);
     if (debug_)
       LogDebug("OutsideInMuonSeeder") << "Query on detector " << it->first->geographicalId().rawId() << " returned "
                                       << mymeas.size() << " measurements." << std::endl;
     for (std::vector<TrajectoryMeasurement>::const_iterator it2 = mymeas.begin(), ed2 = mymeas.end(); it2 != ed2;
          ++it2) {
       if (it2->recHit()->isValid())
         meas.push_back(*it2);
     }
   }
   int found = 0;
   std::sort(meas.begin(), meas.end(), TrajMeasLessEstim());
   for (std::vector<TrajectoryMeasurement>::const_iterator it2 = meas.begin(), ed2 = meas.end(); it2 != ed2; ++it2) {
     if (debug_) {
       LogDebug("OutsideInMuonSeeder") << "  inspecting Hit with chi2 = " << it2->estimate() << std::endl;
       LogDebug("OutsideInMuonSeeder") << "        track state     " << it2->forwardPredictedState().globalPosition()
                                       << std::endl;
       LogDebug("OutsideInMuonSeeder") << "        rechit position " << it2->recHit()->globalPosition() << std::endl;
     }
     TrajectoryStateOnSurface updated = updator_->update(it2->forwardPredictedState(), *it2->recHit());
     if (updated.isValid()) {
       if (debug_)
         LogDebug("OutsideInMuonSeeder") << "          --> updated state: x = " << updated.globalPosition()
                                         << ", p = " << updated.globalMomentum() << std::endl;
       edm::OwnVector<TrackingRecHit> seedHits;
       seedHits.push_back(*it2->recHit()->hit());
       PTrajectoryStateOnDet const &pstate =
           trajectoryStateTransform::persistentState(updated, it2->recHit()->geographicalId().rawId());
       TrajectorySeed seed(pstate, std::move(seedHits), oppositeToMomentum);
       out.push_back(seed);
       found++;
       if (found == hitsToTry_)
         break;
     }
   }
   return found;
 }

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, debug_, errorRescaling_, estimator_, MeasurementDetWithData::fastMeasurements(), newFWLiteAna::found, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), hitsToTry_, MeasurementDetWithData::isNull(), TrajectoryStateOnSurface::isValid(), phase1PixelTopology::layer, TrajectoryStateOnSurface::localError(), LogDebug, HLTSiStripMonitoring_cff::measurementTracker, eostools::move(), oppositeToMomentum, MillePedeFileConverter_cfg::out, trajectoryStateTransform::persistentState(), LocalTrajectoryError::positionError(), edm::OwnVector< T, P >::push_back(), TrajectoryStateOnSurface::rescaleError(), fileCollector::seed, jetUpdater_cfi::sort, mathSSE::sqrt(), TrajectoryStateUpdator::update(), updator_, LocalError::xx(), and LocalError::yy().

Referenced by produce().

◆ produce()

void OutsideInMuonSeeder::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 112 of file OutsideInMuonSeeder.cc.

                                                                                {
   using namespace edm;
   using namespace std;
  
   iSetup.get<IdealMagneticFieldRecord>().get(magfield_);
   iSetup.get<TrackingComponentsRecord>().get(trackerPropagatorName_, trackerPropagator_);
   iSetup.get<TrackingComponentsRecord>().get(muonPropagatorName_, muonPropagator_);
   iSetup.get<GlobalTrackingGeometryRecord>().get(geometry_);
   iSetup.get<TrackingComponentsRecord>().get(estimatorName_, estimator_);
   iSetup.get<TrackingComponentsRecord>().get(updatorName_, updator_);
  
   Handle<MeasurementTrackerEvent> measurementTracker;
   iEvent.getByToken(measurementTrackerTag_, measurementTracker);
  
   ESHandle<TrackerGeometry> tmpTkGeometry;
   iSetup.get<TrackerDigiGeometryRecord>().get(tmpTkGeometry);
  
   Handle<View<reco::Muon>> src;
   iEvent.getByToken(src_, src);
  
   auto out = std::make_unique<std::vector<TrajectorySeed>>();
  
   for (auto const &mu : *src) {
     if (mu.outerTrack().isNull() || !selector_(mu))
       continue;
     if (debug_ && mu.innerTrack().isNonnull())
       doDebug(*mu.innerTrack());
  
     // Better clone here and not directly into doLayer to avoid
     // useless clone/destroy operations to set, in the end, the
     // very same direction every single time.
     std::unique_ptr<Propagator> pmuon_cloned =
         SetPropagationDirection(*muonPropagator_, fromVertex_ ? alongMomentum : oppositeToMomentum);
     std::unique_ptr<Propagator> ptracker_cloned = SetPropagationDirection(*trackerPropagator_, alongMomentum);
  
     int sizeBefore = out->size();
     if (debug_)
       LogDebug("OutsideInMuonSeeder") << "\n\n\nSeeding for muon of pt " << mu.pt() << ", eta " << mu.eta() << ", phi "
                                       << mu.phi() << std::endl;
     const reco::Track &tk = *mu.outerTrack();
  
     TrajectoryStateOnSurface state =
         fromVertex_ ? TrajectoryStateOnSurface(trajectoryStateTransform::initialFreeState(tk, magfield_.product()))
                     : trajectoryStateTransform::innerStateOnSurface(tk, *geometry_, magfield_.product());
  
     if (std::abs(tk.eta()) < maxEtaForTOB_) {
       std::vector<BarrelDetLayer const *> const &tob = measurementTracker->geometricSearchTracker()->tobLayers();
       int found = 0;
       int iLayer = tob.size();
       if (iLayer == 0)
         LogError("OutsideInMuonSeeder") << "TOB has no layers.";
  
       for (auto it = tob.rbegin(), ed = tob.rend(); it != ed; ++it, --iLayer) {
         if (debug_)
           LogDebug("OutsideInMuonSeeder") << "\n ==== Trying TOB " << iLayer << " ====" << std::endl;
         if (doLayer(**it, state, *out, *(pmuon_cloned.get()), *(ptracker_cloned.get()), *measurementTracker)) {
           if (++found == layersToTry_)
             break;
         }
       }
     }
     if (tk.eta() > minEtaForTEC_) {
       const auto &forwLayers = tmpTkGeometry->isThere(GeomDetEnumerators::P2OTEC)
                                    ? measurementTracker->geometricSearchTracker()->posTidLayers()
                                    : measurementTracker->geometricSearchTracker()->posTecLayers();
       if (tmpTkGeometry->isThere(GeomDetEnumerators::P2OTEC)) {
         LogDebug("OutsideInMuonSeeder") << "\n We are using the Phase2 Outer Tracker (defined as a TID+). ";
       }
       LogTrace("OutsideInMuonSeeder") << "\n ==== TEC+ tot layers " << forwLayers.size() << " ====" << std::endl;
       int found = 0;
       int iLayer = forwLayers.size();
       if (iLayer == 0)
         LogError("OutsideInMuonSeeder") << "TEC+ has no layers.";
  
       if (debug_)
         LogDebug("OutsideInMuonSeeder") << "\n ==== Tot layers " << forwLayers.size() << " ====" << std::endl;
       for (auto it = forwLayers.rbegin(), ed = forwLayers.rend(); it != ed; ++it, --iLayer) {
         if (debug_)
           LogDebug("OutsideInMuonSeeder") << "\n ==== Trying Forward Layer +" << +iLayer << " ====" << std::endl;
         if (doLayer(**it, state, *out, *(pmuon_cloned.get()), *(ptracker_cloned.get()), *measurementTracker)) {
           if (++found == layersToTry_)
             break;
         }
       }
     }
     if (tk.eta() < -minEtaForTEC_) {
       const auto &forwLayers = tmpTkGeometry->isThere(GeomDetEnumerators::P2OTEC)
                                    ? measurementTracker->geometricSearchTracker()->negTidLayers()
                                    : measurementTracker->geometricSearchTracker()->negTecLayers();
       if (tmpTkGeometry->isThere(GeomDetEnumerators::P2OTEC)) {
         LogDebug("OutsideInMuonSeeder") << "\n We are using the Phase2 Outer Tracker (defined as a TID-). ";
       }
       LogTrace("OutsideInMuonSeeder") << "\n ==== TEC- tot layers " << forwLayers.size() << " ====" << std::endl;
       int found = 0;
       int iLayer = forwLayers.size();
       if (iLayer == 0)
         LogError("OutsideInMuonSeeder") << "TEC- has no layers.";
  
       if (debug_)
         LogDebug("OutsideInMuonSeeder") << "\n ==== Tot layers " << forwLayers.size() << " ====" << std::endl;
       for (auto it = forwLayers.rbegin(), ed = forwLayers.rend(); it != ed; ++it, --iLayer) {
         if (debug_)
           LogDebug("OutsideInMuonSeeder") << "\n ==== Trying Forward Layer -" << -iLayer << " ====" << std::endl;
         if (doLayer(**it, state, *out, *(pmuon_cloned.get()), *(ptracker_cloned.get()), *measurementTracker)) {
           if (++found == layersToTry_)
             break;
         }
       }
     }
     if (debug_)
       LogDebug("OutsideInMuonSeeder") << "Outcome of seeding for muon of pt " << mu.pt() << ", eta " << mu.eta()
                                       << ", phi " << mu.phi() << ": found " << (out->size() - sizeBefore) << " seeds."
                                       << std::endl;
   }
  
   iEvent.put(std::move(out));
 }

References funct::abs(), alongMomentum, debug_, doDebug(), doLayer(), estimator_, estimatorName_, reco::TrackBase::eta(), newFWLiteAna::found, fromVertex_, geometry_, edm::EventSetup::get(), get, iEvent, trajectoryStateTransform::initialFreeState(), trajectoryStateTransform::innerStateOnSurface(), TrackerGeometry::isThere(), layersToTry_, LogDebug, LogTrace, magfield_, maxEtaForTOB_, HLTSiStripMonitoring_cff::measurementTracker, measurementTrackerTag_, minEtaForTEC_, eostools::move(), amptDefaultParameters_cff::mu, muonPropagator_, muonPropagatorName_, oppositeToMomentum, MillePedeFileConverter_cfg::out, GeomDetEnumerators::P2OTEC, edm::ESHandle< T >::product(), selector_, SetPropagationDirection(), TrackRefitter_38T_cff::src, src_, trackerPropagator_, trackerPropagatorName_, updator_, and updatorName_.

Member Data Documentation

◆ debug_

const bool OutsideInMuonSeeder::debug_

private

Dump deug information.

Definition at line 83 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

◆ errorRescaling_

const double OutsideInMuonSeeder::errorRescaling_

private

How much to rescale errors from STA.

Definition at line 63 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

◆ estimator_

edm::ESHandle<Chi2MeasurementEstimatorBase> OutsideInMuonSeeder::estimator_

private

Definition at line 79 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

◆ estimatorName_

const std::string OutsideInMuonSeeder::estimatorName_

private

Definition at line 69 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ fromVertex_

const bool OutsideInMuonSeeder::fromVertex_

private

Do inside-out.

Definition at line 60 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ geometry_

edm::ESHandle<GlobalTrackingGeometry> OutsideInMuonSeeder::geometry_

private

Definition at line 78 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

◆ hitsToTry_

const int OutsideInMuonSeeder::hitsToTry_

private

How many hits to try on same layer.

Definition at line 57 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

◆ layersToTry_

const int OutsideInMuonSeeder::layersToTry_

private

How many layers to try.

Definition at line 54 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ magfield_

edm::ESHandle<MagneticField> OutsideInMuonSeeder::magfield_

private

Definition at line 75 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

◆ maxEtaForTOB_

const float OutsideInMuonSeeder::maxEtaForTOB_

private

Definition at line 73 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ measurementTrackerName_

const std::string OutsideInMuonSeeder::measurementTrackerName_

private

Definition at line 68 of file OutsideInMuonSeeder.cc.

◆ measurementTrackerTag_

edm::EDGetTokenT<MeasurementTrackerEvent> OutsideInMuonSeeder::measurementTrackerTag_

private

Definition at line 67 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ minEtaForTEC_

const float OutsideInMuonSeeder::minEtaForTEC_

private

Definition at line 72 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ muonPropagator_

edm::ESHandle<Propagator> OutsideInMuonSeeder::muonPropagator_

private

Definition at line 76 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

◆ muonPropagatorName_

const std::string OutsideInMuonSeeder::muonPropagatorName_

private

Definition at line 66 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ selector_

StringCutObjectSelector<reco::Muon> OutsideInMuonSeeder::selector_

private

Muon selection.

Definition at line 51 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ src_

edm::EDGetTokenT<edm::View<reco::Muon> > OutsideInMuonSeeder::src_

private

Labels for input collections.

Definition at line 48 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ trackerPropagator_

edm::ESHandle<Propagator> OutsideInMuonSeeder::trackerPropagator_

private

Definition at line 77 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ trackerPropagatorName_

const std::string OutsideInMuonSeeder::trackerPropagatorName_

private

Definition at line 65 of file OutsideInMuonSeeder.cc.

Referenced by produce().

◆ updator_

edm::ESHandle<TrajectoryStateUpdator> OutsideInMuonSeeder::updator_

private

Definition at line 80 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

◆ updatorName_

const std::string OutsideInMuonSeeder::updatorName_

private

Definition at line 70 of file OutsideInMuonSeeder.cc.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ OutsideInMuonSeeder()

◆ ~OutsideInMuonSeeder()

Member Function Documentation

◆ doDebug()

◆ doLayer()

◆ produce()

Member Data Documentation

◆ debug_

◆ errorRescaling_

◆ estimator_

◆ estimatorName_

◆ fromVertex_

◆ geometry_

◆ hitsToTry_

◆ layersToTry_

◆ magfield_

◆ maxEtaForTOB_

◆ measurementTrackerName_

◆ measurementTrackerTag_

◆ minEtaForTEC_

◆ muonPropagator_

◆ muonPropagatorName_

◆ selector_

◆ src_

◆ trackerPropagator_

◆ trackerPropagatorName_

◆ updator_

◆ updatorName_