Inheritance diagram for OutsideInMuonSeeder:

Public Member Functions
	OutsideInMuonSeeder (const edm::ParameterSet &iConfig)

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

virtual	~OutsideInMuonSeeder ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	doDebug (const reco::Track &tk) const

int	doLayer (const GeometricSearchDet &layer, const TrajectoryStateOnSurface &state, std::vector< TrajectorySeed > &out, const MeasurementTrackerEvent &mte) const

Private Attributes
bool	debug_
	Dump deug information. More...

Plane::PlanePointer	dummyPlane_
	Surface used to make a TSOS at the PCA to the beamline. More...

double	errorRescaling_
	How much to rescale errors from STA. More...

edm::ESHandle < Chi2MeasurementEstimatorBase >	estimator_

std::string	estimatorName_

bool	fromVertex_
	Do inside-out. More...

edm::ESHandle < GlobalTrackingGeometry >	geometry_

int	hitsToTry_
	How many hits to try on same layer. More...

int	layersToTry_
	How many layers to try. More...

edm::ESHandle< MagneticField >	magfield_

double	maxEtaForTOB_

std::string	measurementTrackerName_

edm::EDGetTokenT < MeasurementTrackerEvent >	measurementTrackerTag_

double	minEtaForTEC_

edm::ESHandle< Propagator >	muonPropagator_

std::string	muonPropagatorName_

StringCutObjectSelector < reco::Muon >	selector_
	Muon selection. More...

edm::InputTag	src_
	Labels for input collections. More...

edm::ESHandle< Propagator >	trackerPropagator_

std::string	trackerPropagatorName_

edm::ESHandle < TrajectoryStateUpdator >	updator_

std::string	updatorName_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 44 of file OutsideInMuonSeeder.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 97 of file OutsideInMuonSeeder.cc.

References updatorName_.

                                                                         :
     src_(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")),
     minEtaForTEC_(iConfig.getParameter<double>("minEtaForTEC")),
     maxEtaForTOB_(iConfig.getParameter<double>("maxEtaForTOB")),
     debug_(iConfig.getUntrackedParameter<bool>("debug",false)),
     dummyPlane_(Plane::build(Plane::PositionType(), Plane::RotationType()))
 {
     produces<std::vector<TrajectorySeed> >(); 
     updatorName_ = "KFUpdator";
 }

virtual OutsideInMuonSeeder::~OutsideInMuonSeeder ( )

inlinevirtual

Definition at line 47 of file OutsideInMuonSeeder.cc.

47 { }

Member Function Documentation

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

private

Definition at line 242 of file OutsideInMuonSeeder.cc.

References alongMomentum, gather_cfg::cout, TrackingRecHit::geographicalId(), geometry_, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), i, trajectoryStateTransform::innerStateOnSurface(), TrackingRecHit::isValid(), TrajectoryStateOnSurface::isValid(), magfield_, muonPropagator_, DetId::rawId(), reco::Track::recHit(), reco::Track::recHitsSize(), and GeomDet::surface().

Referenced by produce().

                                                       {
     TrajectoryStateOnSurface tsos = trajectoryStateTransform::innerStateOnSurface(tk, *geometry_, &*magfield_);
     muonPropagator_->setPropagationDirection(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 == 0) continue;
         if (i != 0) tsos = muonPropagator_->propagate(tsos, det->surface());
         if (!tsos.isValid()) continue;
         std::cout << "  state " << i << " at x = " << tsos.globalPosition() << ", p = " << tsos.globalMomentum() << std::endl;
         if (hit->isValid()) {
             std::cout << "         valid   rechit on detid " << hit->geographicalId().rawId() << std::endl;
         } else {
             std::cout << "         invalid rechit on detid " << hit->geographicalId().rawId() << std::endl;
         }
     }
 }

int OutsideInMuonSeeder::doLayer	(	const GeometricSearchDet &	layer,
		const TrajectoryStateOnSurface &	state,
		std::vector< TrajectorySeed > &	out,
		const MeasurementTrackerEvent &	mte
	)		const

private

Definition at line 196 of file OutsideInMuonSeeder.cc.

References dtNoiseDBValidation_cfg::cerr, GeometricSearchDet::compatibleDetsV(), gather_cfg::cout, debug_, errorRescaling_, estimator_, MeasurementDetWithData::fastMeasurements(), newFWLiteAna::found, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), hitsToTry_, MeasurementTrackerEvent::idToDet(), MeasurementDetWithData::isNull(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), muonPropagator_, oppositeToMomentum, trajectoryStateTransform::persistentState(), LocalTrajectoryError::positionError(), edm::OwnVector< T, P >::push_back(), TrajectoryStateOnSurface::rescaleError(), python.multivaluedict::sort(), mathSSE::sqrt(), trackerPropagator_, updator_, LocalError::xx(), and LocalError::yy().

Referenced by produce().

                                                                                                                                                                                             {
     TrajectoryStateOnSurface onLayer(state);
     onLayer.rescaleError(errorRescaling_);
     std::vector< GeometricSearchDet::DetWithState > dets;
     layer.compatibleDetsV(onLayer, *muonPropagator_, *estimator_, dets); 
 
     if (debug_) {
         std::cout << "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, *trackerPropagator_, *estimator_);
         if (debug_) std::cout << "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_) {
             std::cout << "  inspecting Hit with chi2 = " << it2->estimate() << std::endl;
             std::cout << "        track state     " << it2->forwardPredictedState().globalPosition() << std::endl; 
             std::cout << "        rechit position " << it2->recHit()->globalPosition() << std::endl; 
         }
         TrajectoryStateOnSurface updated = updator_->update(it2->forwardPredictedState(), *it2->recHit());
         if (updated.isValid()) {
             if (debug_) std::cout << "          --> 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
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 118 of file OutsideInMuonSeeder.cc.

References funct::abs(), alongMomentum, gather_cfg::cout, debug_, doDebug(), doLayer(), dummyPlane_, estimator_, estimatorName_, reco::TrackBase::eta(), reco::LeafCandidate::eta(), newFWLiteAna::found, fromVertex_, geometry_, edm::EventSetup::get(), edm::Event::getByLabel(), edm::Event::getByToken(), trajectoryStateTransform::initialFreeState(), trajectoryStateTransform::innerStateOnSurface(), reco::Muon::innerTrack(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::isNull(), layersToTry_, magfield_, maxEtaForTOB_, measurementTrackerTag_, minEtaForTEC_, RPCpg::mu, muonPropagator_, muonPropagatorName_, oppositeToMomentum, dbtoconf::out, reco::Muon::outerTrack(), reco::LeafCandidate::phi(), FreeTrajectoryState::position(), edm::ESHandle< class >::product(), reco::LeafCandidate::pt(), edm::Event::put(), selector_, alcazmumu_cfi::src, src_, trackerPropagator_, trackerPropagatorName_, updator_, and updatorName_.

                                                                             {
     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);
 
     Handle<View<reco::Muon> > src;
     iEvent.getByLabel(src_, src);
 
     
     auto_ptr<vector<TrajectorySeed> > out(new vector<TrajectorySeed>());
 
     for (View<reco::Muon>::const_iterator it = src->begin(), ed = src->end(); it != ed; ++it) {
         const reco::Muon &mu = *it;
         if (mu.outerTrack().isNull() || !selector_(mu)) continue;
         if (debug_ && mu.innerTrack().isNonnull()) doDebug(*mu.innerTrack());
 
         muonPropagator_->setPropagationDirection(fromVertex_ ? alongMomentum : oppositeToMomentum);
         trackerPropagator_->setPropagationDirection(alongMomentum);
 
         int sizeBefore = out->size();
         if (debug_) std::cout << "\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;
         if (fromVertex_) {
             FreeTrajectoryState fstate = trajectoryStateTransform::initialFreeState(tk, magfield_.product());
             dummyPlane_->move(fstate.position() - dummyPlane_->position());
             state = TrajectoryStateOnSurface(fstate, *dummyPlane_);
         } else {
             state = trajectoryStateTransform::innerStateOnSurface(tk, *geometry_, magfield_.product());
         }
         if (std::abs(tk.eta()) < maxEtaForTOB_) {
             std::vector< BarrelDetLayer * > const & tob = measurementTracker->geometricSearchTracker()->tobLayers();
             int iLayer = 6, found = 0;
             for (std::vector<BarrelDetLayer *>::const_reverse_iterator it = tob.rbegin(), ed = tob.rend(); it != ed; ++it, --iLayer) {
                 if (debug_) std::cout << "\n ==== Trying TOB " << iLayer << " ====" << std::endl;
                 if (doLayer(**it, state, *out, *measurementTracker)) {
                     if (++found == layersToTry_) break;
                 }
             }
         }
         if (tk.eta() > minEtaForTEC_) {
             int iLayer = 9, found = 0;
             std::vector< ForwardDetLayer * > const & tec = measurementTracker->geometricSearchTracker()->posTecLayers();
             for (std::vector<ForwardDetLayer *>::const_reverse_iterator it = tec.rbegin(), ed = tec.rend(); it != ed; ++it, --iLayer) {
                 if (debug_) std::cout << "\n ==== Trying TEC " << +iLayer << " ====" << std::endl;
                 if (doLayer(**it, state, *out, *measurementTracker)) {
                     if (++found == layersToTry_) break;
                 }
             }
         }
         if (tk.eta() < -minEtaForTEC_) {
             int iLayer = 9, found = 0;
             std::vector< ForwardDetLayer * > const & tec = measurementTracker->geometricSearchTracker()->negTecLayers();
             for (std::vector<ForwardDetLayer *>::const_reverse_iterator it = tec.rbegin(), ed = tec.rend(); it != ed; ++it, --iLayer) {
                 if (debug_) std::cout << "\n ==== Trying TEC " << -iLayer << " ====" << std::endl;
                 if (doLayer(**it, state, *out, *measurementTracker)) {
                     if (++found == layersToTry_) break;
                 }
             }
         }
         if (debug_) std::cout << "Outcome of seeding for muon of pt " << mu.pt() << ", eta " << mu.eta() << ", phi " << mu.phi() << ": found " << (out->size() - sizeBefore) << " seeds."<< std::endl;
         
     }
 
     iEvent.put(out);
 }

Member Data Documentation

bool OutsideInMuonSeeder::debug_

private

Dump deug information.

Definition at line 87 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

Plane::PlanePointer OutsideInMuonSeeder::dummyPlane_

private

Surface used to make a TSOS at the PCA to the beamline.

Definition at line 90 of file OutsideInMuonSeeder.cc.

Referenced by produce().

double OutsideInMuonSeeder::errorRescaling_

private

How much to rescale errors from STA.

Definition at line 68 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

edm::ESHandle<Chi2MeasurementEstimatorBase> OutsideInMuonSeeder::estimator_

private

Definition at line 83 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

std::string OutsideInMuonSeeder::estimatorName_

private

Definition at line 74 of file OutsideInMuonSeeder.cc.

Referenced by produce().

bool OutsideInMuonSeeder::fromVertex_

private

Do inside-out.

Definition at line 65 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<GlobalTrackingGeometry> OutsideInMuonSeeder::geometry_

private

Definition at line 82 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

int OutsideInMuonSeeder::hitsToTry_

private

How many hits to try on same layer.

Definition at line 62 of file OutsideInMuonSeeder.cc.

Referenced by doLayer().

int OutsideInMuonSeeder::layersToTry_

private

How many layers to try.

Definition at line 59 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<MagneticField> OutsideInMuonSeeder::magfield_

private

Definition at line 79 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), and produce().

double OutsideInMuonSeeder::maxEtaForTOB_

private

Definition at line 77 of file OutsideInMuonSeeder.cc.

Referenced by produce().

std::string OutsideInMuonSeeder::measurementTrackerName_

private

Definition at line 73 of file OutsideInMuonSeeder.cc.

edm::EDGetTokenT<MeasurementTrackerEvent> OutsideInMuonSeeder::measurementTrackerTag_

private

Definition at line 72 of file OutsideInMuonSeeder.cc.

Referenced by produce().

double OutsideInMuonSeeder::minEtaForTEC_

private

Definition at line 77 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<Propagator> OutsideInMuonSeeder::muonPropagator_

private

Definition at line 80 of file OutsideInMuonSeeder.cc.

Referenced by doDebug(), doLayer(), and produce().

std::string OutsideInMuonSeeder::muonPropagatorName_

private

Definition at line 71 of file OutsideInMuonSeeder.cc.

Referenced by produce().

StringCutObjectSelector<reco::Muon> OutsideInMuonSeeder::selector_

private

Muon selection.

Definition at line 56 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::InputTag OutsideInMuonSeeder::src_

private

Labels for input collections.

Definition at line 53 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<Propagator> OutsideInMuonSeeder::trackerPropagator_

private

Definition at line 81 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

std::string OutsideInMuonSeeder::trackerPropagatorName_

private

Definition at line 70 of file OutsideInMuonSeeder.cc.

Referenced by produce().

edm::ESHandle<TrajectoryStateUpdator> OutsideInMuonSeeder::updator_

private

Definition at line 84 of file OutsideInMuonSeeder.cc.

Referenced by doLayer(), and produce().

std::string OutsideInMuonSeeder::updatorName_

private

Definition at line 75 of file OutsideInMuonSeeder.cc.

Referenced by OutsideInMuonSeeder(), and produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation