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...
const Chi2MeasurementEstimatorBase *	estimator_
const edm::ESGetToken< Chi2MeasurementEstimatorBase, TrackingComponentsRecord >	estimatorToken_
const bool	fromVertex_
	Do inside-out. More...
const GlobalTrackingGeometry *	geometry_
const edm::ESGetToken< GlobalTrackingGeometry, GlobalTrackingGeometryRecord >	geometryToken_
const int	hitsToTry_
	How many hits to try on same layer. More...
const int	layersToTry_
	How many layers to try. More...
const MagneticField *	magfield_
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magfieldToken_
float const	maxEtaForTOB_
const edm::EDGetTokenT< MeasurementTrackerEvent >	measurementTrackerTag_
float const	minEtaForTEC_
const Propagator *	muonPropagator_
const edm::ESGetToken< Propagator, TrackingComponentsRecord >	muonPropagatorToken_
StringCutObjectSelector< reco::Muon >	selector_
	Muon selection. More...
edm::EDGetTokenT< edm::View< reco::Muon > >	src_
	Labels for input collections. More...
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeometryToken_
const edm::ESGetToken< Propagator, TrackingComponentsRecord >	trackerPropagatorToken_
const TrajectoryStateUpdator *	updator_
const edm::ESGetToken< TrajectoryStateUpdator, TrackingComponentsRecord >	updatorToken_

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 38 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")),
      trackerPropagatorToken_(esConsumes(edm::ESInputTag("", iConfig.getParameter<std::string>("trackerPropagator")))),
      muonPropagatorToken_(esConsumes(edm::ESInputTag("", iConfig.getParameter<std::string>("muonPropagator")))),
      measurementTrackerTag_(consumes<MeasurementTrackerEvent>(edm::InputTag("MeasurementTrackerEvent"))),
      estimatorToken_(esConsumes(edm::ESInputTag("", iConfig.getParameter<std::string>("hitCollector")))),
      updatorToken_(esConsumes(edm::ESInputTag("", "KFUpdator"))),
      magfieldToken_(esConsumes()),
      geometryToken_(esConsumes()),
      tkGeometryToken_(esConsumes()),
      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 41 of file OutsideInMuonSeeder.cc.

{}

Member Function Documentation

doDebug()

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

private

Definition at line 299 of file OutsideInMuonSeeder.cc.

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

                                                           {
  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;
    }
  }
}

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 232 of file OutsideInMuonSeeder.cc.

References DMR_cfg::cerr, debug_, errorRescaling_, estimator_, MeasurementDetWithData::fastMeasurements(), newFWLiteAna::found, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), hitsToTry_, MeasurementDetWithData::isNull(), TrajectoryStateOnSurface::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, nano_mu_digi_cff::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().

                                                                                          {
  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;
}

produce()

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

override

Definition at line 115 of file OutsideInMuonSeeder.cc.

References funct::abs(), alongMomentum, debug_, doDebug(), doLayer(), estimator_, estimatorToken_, reco::TrackBase::eta(), newFWLiteAna::found, fromVertex_, geometry_, geometryToken_, edm::EventSetup::getData(), iEvent, trajectoryStateTransform::initialFreeState(), trajectoryStateTransform::innerStateOnSurface(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, layersToTry_, LogDebug, LogTrace, magfield_, magfieldToken_, maxEtaForTOB_, HLTSiStripMonitoring_cff::measurementTracker, measurementTrackerTag_, minEtaForTEC_, eostools::move(), amptDefaultParameters_cff::mu, muonPropagator_, muonPropagatorToken_, oppositeToMomentum, MillePedeFileConverter_cfg::out, GeomDetEnumerators::P2OTEC, selector_, SetPropagationDirection(), TrackRefitter_38T_cff::src, src_, tkGeometryToken_, outInSeedsFromStandaloneMuons_cfi::trackerPropagator, trackerPropagatorToken_, updator_, and updatorToken_.

                                                                               {
  using namespace edm;
  using namespace std;

  magfield_ = &iSetup.getData(magfieldToken_);
  auto const &trackerPropagator = iSetup.getData(trackerPropagatorToken_);
  muonPropagator_ = &iSetup.getData(muonPropagatorToken_);
  geometry_ = &iSetup.getData(geometryToken_);
  estimator_ = &iSetup.getData(estimatorToken_);
  updator_ = &iSetup.getData(updatorToken_);

  Handle<MeasurementTrackerEvent> measurementTracker;
  iEvent.getByToken(measurementTrackerTag_, measurementTracker);

  const auto &tmpTkGeometry = iSetup.getData(tkGeometryToken_);

  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_))
                    : trajectoryStateTransform::innerStateOnSurface(tk, *geometry_, magfield_);

    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));
}

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 62 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

estimator_

const Chi2MeasurementEstimatorBase* OutsideInMuonSeeder::estimator_

private

Definition at line 79 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

estimatorToken_

const edm::ESGetToken<Chi2MeasurementEstimatorBase, TrackingComponentsRecord> OutsideInMuonSeeder::estimatorToken_

private

Definition at line 67 of file OutsideInMuonSeeder.cc.

Referenced by produce().

fromVertex_

const bool OutsideInMuonSeeder::fromVertex_

private

Do inside-out.

Definition at line 59 of file OutsideInMuonSeeder.cc.

Referenced by produce().

geometry_

const GlobalTrackingGeometry* OutsideInMuonSeeder::geometry_

private

Definition at line 78 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

geometryToken_

const edm::ESGetToken<GlobalTrackingGeometry, GlobalTrackingGeometryRecord> OutsideInMuonSeeder::geometryToken_

private

Definition at line 70 of file OutsideInMuonSeeder.cc.

Referenced by produce().

hitsToTry_

const int OutsideInMuonSeeder::hitsToTry_

private

How many hits to try on same layer.

Definition at line 56 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

layersToTry_

const int OutsideInMuonSeeder::layersToTry_

private

How many layers to try.

Definition at line 53 of file OutsideInMuonSeeder.cc.

Referenced by produce().

magfield_

const MagneticField* OutsideInMuonSeeder::magfield_

private

Definition at line 76 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

magfieldToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> OutsideInMuonSeeder::magfieldToken_

private

Definition at line 69 of file OutsideInMuonSeeder.cc.

Referenced by produce().

maxEtaForTOB_

float const OutsideInMuonSeeder::maxEtaForTOB_

private

Definition at line 74 of file OutsideInMuonSeeder.cc.

Referenced by produce().

measurementTrackerTag_

const edm::EDGetTokenT<MeasurementTrackerEvent> OutsideInMuonSeeder::measurementTrackerTag_

private

Definition at line 66 of file OutsideInMuonSeeder.cc.

Referenced by produce().

minEtaForTEC_

float const OutsideInMuonSeeder::minEtaForTEC_

private

Definition at line 73 of file OutsideInMuonSeeder.cc.

Referenced by produce().

muonPropagator_

const Propagator* OutsideInMuonSeeder::muonPropagator_

private

Definition at line 77 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

muonPropagatorToken_

const edm::ESGetToken<Propagator, TrackingComponentsRecord> OutsideInMuonSeeder::muonPropagatorToken_

private

Definition at line 65 of file OutsideInMuonSeeder.cc.

Referenced by produce().

selector_

StringCutObjectSelector<reco::Muon> OutsideInMuonSeeder::selector_

private

Muon selection.

Definition at line 50 of file OutsideInMuonSeeder.cc.

Referenced by produce().

src_

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

private

Labels for input collections.

Definition at line 47 of file OutsideInMuonSeeder.cc.

Referenced by produce().

tkGeometryToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> OutsideInMuonSeeder::tkGeometryToken_

private

Definition at line 71 of file OutsideInMuonSeeder.cc.

Referenced by produce().

trackerPropagatorToken_

const edm::ESGetToken<Propagator, TrackingComponentsRecord> OutsideInMuonSeeder::trackerPropagatorToken_

private

Definition at line 64 of file OutsideInMuonSeeder.cc.

Referenced by produce().

updator_

const TrajectoryStateUpdator* OutsideInMuonSeeder::updator_

private

Definition at line 80 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

updatorToken_

const edm::ESGetToken<TrajectoryStateUpdator, TrackingComponentsRecord> OutsideInMuonSeeder::updatorToken_

private

Definition at line 68 of file OutsideInMuonSeeder.cc.

Referenced by produce().