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

inlineoverride

Definition at line 41 of file OutsideInMuonSeeder.cc.

41 {}

Member Function Documentation

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

private

Definition at line 299 of file OutsideInMuonSeeder.cc.

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

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::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(), fileCollector::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 115 of file OutsideInMuonSeeder.cc.

References funct::abs(), alongMomentum, debug_, doDebug(), doLayer(), estimator_, estimatorToken_, reco::TrackBase::eta(), newFWLiteAna::found, fromVertex_, geometry_, geometryToken_, edm::Event::getByToken(), edm::EventSetup::getData(), trajectoryStateTransform::initialFreeState(), trajectoryStateTransform::innerStateOnSurface(), layersToTry_, LogDebug, LogTrace, magfield_, magfieldToken_, maxEtaForTOB_, HLT_FULL_cff::measurementTracker, measurementTrackerTag_, minEtaForTEC_, eostools::move(), RPCpg::mu, muonPropagator_, muonPropagatorToken_, oppositeToMomentum, submitPVResolutionJobs::out, GeomDetEnumerators::P2OTEC, edm::Event::put(), selector_, SetPropagationDirection(), alcazmumu_cfi::src, src_, tkGeometryToken_, 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

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

Referenced by doLayer().

const Chi2MeasurementEstimatorBase* OutsideInMuonSeeder::estimator_

private

Definition at line 79 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

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

private

Definition at line 67 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const bool OutsideInMuonSeeder::fromVertex_

private

Do inside-out.

Definition at line 59 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const GlobalTrackingGeometry* OutsideInMuonSeeder::geometry_

private

Definition at line 78 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

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

private

Definition at line 70 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const int OutsideInMuonSeeder::hitsToTry_

private

How many hits to try on same layer.

Definition at line 56 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

const int OutsideInMuonSeeder::layersToTry_

private

How many layers to try.

Definition at line 53 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const MagneticField* OutsideInMuonSeeder::magfield_

private

Definition at line 76 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

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

private

Definition at line 69 of file OutsideInMuonSeeder.cc.

Referenced by produce().

float const OutsideInMuonSeeder::maxEtaForTOB_

private

Definition at line 74 of file OutsideInMuonSeeder.cc.

Referenced by produce().

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

private

Definition at line 66 of file OutsideInMuonSeeder.cc.

Referenced by produce().

float const OutsideInMuonSeeder::minEtaForTEC_

private

Definition at line 73 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const Propagator* OutsideInMuonSeeder::muonPropagator_

private

Definition at line 77 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

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

private

Definition at line 65 of file OutsideInMuonSeeder.cc.

Referenced by produce().

StringCutObjectSelector<reco::Muon> OutsideInMuonSeeder::selector_

private

Muon selection.

Definition at line 50 of file OutsideInMuonSeeder.cc.

Referenced by produce().

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

private

Labels for input collections.

Definition at line 47 of file OutsideInMuonSeeder.cc.

Referenced by produce().

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

private

Definition at line 71 of file OutsideInMuonSeeder.cc.

Referenced by produce().

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

private

Definition at line 64 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const TrajectoryStateUpdator* OutsideInMuonSeeder::updator_

private

Definition at line 80 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

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

private

Definition at line 68 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