OutsideInMuonSeeder Class Reference

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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::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.

{}

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().