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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () 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_

float const	maxEtaForTOB_

const std::string	measurementTrackerName_

edm::EDGetTokenT< MeasurementTrackerEvent >	measurementTrackerTag_

float const	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 ( 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 ( )

inlineoverride

Definition at line 42 of file OutsideInMuonSeeder.cc.

References iEvent, and produce().

42 {}

Member Function Documentation

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

private

Definition at line 297 of file OutsideInMuonSeeder.cc.

References alongMomentum, DEFINE_FWK_MODULE, TrackingRecHit::geographicalId(), geometry_, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), mps_fire::i, GlobalTrackingGeometry::idToDet(), trajectoryStateTransform::innerStateOnSurface(), TrajectoryStateOnSurface::isValid(), TrackingRecHit::isValid(), LogDebug, magfield_, muonPropagator_, DetId::rawId(), 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;
     }
   }
 }

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.

References beam_dqm_sourceclient-live_cfg::cerr, GeometricSearchDet::compatibleDetsV(), debug_, errorRescaling_, estimator_, MeasurementDetWithData::fastMeasurements(), newFWLiteAna::found, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), hitsToTry_, MeasurementTrackerEvent::idToDet(), MeasurementDetWithData::isNull(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), LogDebug, eostools::move(), oppositeToMomentum, trajectoryStateTransform::persistentState(), LocalTrajectoryError::positionError(), edm::OwnVector< T, P >::push_back(), TrajectoryStateOnSurface::rescaleError(), SurveyInfoScenario_cff::seed, 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;
 }

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

override

Definition at line 112 of file OutsideInMuonSeeder.cc.

References funct::abs(), alongMomentum, debug_, doDebug(), doLayer(), estimator_, estimatorName_, reco::TrackBase::eta(), newFWLiteAna::found, fromVertex_, geometry_, edm::EventSetup::get(), edm::Event::getByToken(), 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(), edm::Event::put(), selector_, SetPropagationDirection(), TrackRefitter_38T_cff::src, src_, trackerPropagator_, trackerPropagatorName_, updator_, and updatorName_.

Referenced by ~OutsideInMuonSeeder().

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

Member Data Documentation

const bool OutsideInMuonSeeder::debug_

private

Dump deug information.

Definition at line 83 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

const double OutsideInMuonSeeder::errorRescaling_

private

How much to rescale errors from STA.

Definition at line 63 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

edm::ESHandle<Chi2MeasurementEstimatorBase> OutsideInMuonSeeder::estimator_

private

Definition at line 79 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

const std::string OutsideInMuonSeeder::estimatorName_

private

Definition at line 69 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const bool OutsideInMuonSeeder::fromVertex_

private

Do inside-out.

Definition at line 60 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<GlobalTrackingGeometry> OutsideInMuonSeeder::geometry_

private

Definition at line 78 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

const int OutsideInMuonSeeder::hitsToTry_

private

How many hits to try on same layer.

Definition at line 57 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

const int OutsideInMuonSeeder::layersToTry_

private

How many layers to try.

Definition at line 54 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<MagneticField> OutsideInMuonSeeder::magfield_

private

Definition at line 75 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

float const OutsideInMuonSeeder::maxEtaForTOB_

private

Definition at line 73 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const std::string OutsideInMuonSeeder::measurementTrackerName_

private

Definition at line 68 of file OutsideInMuonSeeder.cc.

edm::EDGetTokenT<MeasurementTrackerEvent> OutsideInMuonSeeder::measurementTrackerTag_

private

Definition at line 67 of file OutsideInMuonSeeder.cc.

Referenced by produce().

float const OutsideInMuonSeeder::minEtaForTEC_

private

Definition at line 72 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<Propagator> OutsideInMuonSeeder::muonPropagator_

private

Definition at line 76 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

const std::string OutsideInMuonSeeder::muonPropagatorName_

private

Definition at line 66 of file OutsideInMuonSeeder.cc.

Referenced by produce().

StringCutObjectSelector<reco::Muon> OutsideInMuonSeeder::selector_

private

Muon selection.

Definition at line 51 of file OutsideInMuonSeeder.cc.

Referenced by produce().

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

private

Labels for input collections.

Definition at line 48 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<Propagator> OutsideInMuonSeeder::trackerPropagator_

private

Definition at line 77 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const std::string OutsideInMuonSeeder::trackerPropagatorName_

private

Definition at line 65 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<TrajectoryStateUpdator> OutsideInMuonSeeder::updator_

private

Definition at line 80 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

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

Member Function Documentation

Member Data Documentation