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

bool	hasAbilityToProduceInRuns () 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 40 of file OutsideInMuonSeeder.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 96 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 43 of file OutsideInMuonSeeder.cc.

References iEvent, and produce().

43 { }

Member Function Documentation

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

private

Definition at line 277 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 226 of file OutsideInMuonSeeder.cc.

References MessageLogger_cfi::cerr, GeometricSearchDet::compatibleDetsV(), debug_, errorRescaling_, estimator_, MeasurementDetWithData::fastMeasurements(), runEdmFileComparison::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 116 of file OutsideInMuonSeeder.cc.

References funct::abs(), alongMomentum, debug_, doDebug(), doLayer(), estimator_, estimatorName_, reco::TrackBase::eta(), runEdmFileComparison::found, fromVertex_, geometry_, edm::EventSetup::get(), edm::Event::getByToken(), trajectoryStateTransform::initialFreeState(), trajectoryStateTransform::innerStateOnSurface(), TrackerGeometry::isThere(), layersToTry_, LogDebug, LogTrace, magfield_, maxEtaForTOB_, measurementTrackerTag_, minEtaForTEC_, eostools::move(), RPCpg::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 84 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

const double OutsideInMuonSeeder::errorRescaling_

private

How much to rescale errors from STA.

Definition at line 64 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

edm::ESHandle<Chi2MeasurementEstimatorBase> OutsideInMuonSeeder::estimator_

private

Definition at line 80 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

const std::string OutsideInMuonSeeder::estimatorName_

private

Definition at line 70 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const bool OutsideInMuonSeeder::fromVertex_

private

Do inside-out.

Definition at line 61 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<GlobalTrackingGeometry> OutsideInMuonSeeder::geometry_

private

Definition at line 79 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 58 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

const int OutsideInMuonSeeder::layersToTry_

private

How many layers to try.

Definition at line 55 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<MagneticField> OutsideInMuonSeeder::magfield_

private

Definition at line 76 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

float const OutsideInMuonSeeder::maxEtaForTOB_

private

Definition at line 74 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const std::string OutsideInMuonSeeder::measurementTrackerName_

private

Definition at line 69 of file OutsideInMuonSeeder.cc.

edm::EDGetTokenT<MeasurementTrackerEvent> OutsideInMuonSeeder::measurementTrackerTag_

private

Definition at line 68 of file OutsideInMuonSeeder.cc.

Referenced by produce().

float const OutsideInMuonSeeder::minEtaForTEC_

private

Definition at line 73 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<Propagator> OutsideInMuonSeeder::muonPropagator_

private

Definition at line 77 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

const std::string OutsideInMuonSeeder::muonPropagatorName_

private

Definition at line 67 of file OutsideInMuonSeeder.cc.

Referenced by produce().

StringCutObjectSelector<reco::Muon> OutsideInMuonSeeder::selector_

private

Muon selection.

Definition at line 52 of file OutsideInMuonSeeder.cc.

Referenced by produce().

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

private

Labels for input collections.

Definition at line 49 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<Propagator> OutsideInMuonSeeder::trackerPropagator_

private

Definition at line 78 of file OutsideInMuonSeeder.cc.

Referenced by produce().

const std::string OutsideInMuonSeeder::trackerPropagatorName_

private

Definition at line 66 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<TrajectoryStateUpdator> OutsideInMuonSeeder::updator_

private

Definition at line 81 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

const std::string OutsideInMuonSeeder::updatorName_

private

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