CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
OutsideInMuonSeeder Class Referencefinal
Inheritance diagram for OutsideInMuonSeeder:
edm::stream::EDProducer<> edm::stream::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

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
 
- Public Member Functions inherited from edm::stream::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDProducerBase () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
virtual ~ProducerBase () noexcept(false)
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

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< Chi2MeasurementEstimatorBaseestimator_
 
const std::string estimatorName_
 
const bool fromVertex_
 Do inside-out. More...
 
edm::ESHandle< GlobalTrackingGeometrygeometry_
 
const int hitsToTry_
 How many hits to try on same layer. More...
 
const int layersToTry_
 How many layers to try. More...
 
edm::ESHandle< MagneticFieldmagfield_
 
float const maxEtaForTOB_
 
const std::string measurementTrackerName_
 
edm::EDGetTokenT< MeasurementTrackerEventmeasurementTrackerTag_
 
float const minEtaForTEC_
 
edm::ESHandle< PropagatormuonPropagator_
 
const std::string muonPropagatorName_
 
StringCutObjectSelector< reco::Muonselector_
 Muon selection. More...
 
edm::EDGetTokenT< edm::View< reco::Muon > > src_
 Labels for input collections. More...
 
edm::ESHandle< PropagatortrackerPropagator_
 
const std::string trackerPropagatorName_
 
edm::ESHandle< TrajectoryStateUpdatorupdator_
 
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, GlobalCacheLuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCacheRunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 
- Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- 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 40 of file OutsideInMuonSeeder.cc.

Constructor & Destructor Documentation

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

Definition at line 96 of file OutsideInMuonSeeder.cc.

96  :
98  selector_(iConfig.existsAs<std::string>("cut") ? iConfig.getParameter<std::string>("cut") : "", true),
99  layersToTry_(iConfig.getParameter<int32_t>("layersToTry")),
100  hitsToTry_(iConfig.getParameter<int32_t>("hitsToTry")),
101  fromVertex_(iConfig.getParameter<bool>("fromVertex")),
102  errorRescaling_(iConfig.getParameter<double>("errorRescaleFactor")),
103  trackerPropagatorName_(iConfig.getParameter<std::string>("trackerPropagator")),
104  muonPropagatorName_(iConfig.getParameter<std::string>("muonPropagator")),
105  measurementTrackerTag_(consumes<MeasurementTrackerEvent>(edm::InputTag("MeasurementTrackerEvent"))),
106  estimatorName_(iConfig.getParameter<std::string>("hitCollector")),
107  updatorName_("KFUpdator"),
108  minEtaForTEC_(iConfig.getParameter<double>("minEtaForTEC")),
109  maxEtaForTOB_(iConfig.getParameter<double>("maxEtaForTOB")),
110  debug_(iConfig.getUntrackedParameter<bool>("debug",false))
111 {
112  produces<std::vector<TrajectorySeed> >();
113 }
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
const std::string estimatorName_
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:186
const bool fromVertex_
Do inside-out.
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
const bool debug_
Dump deug information.
const double errorRescaling_
How much to rescale errors from STA.
const std::string trackerPropagatorName_
const std::string updatorName_
StringCutObjectSelector< reco::Muon > selector_
Muon selection.
edm::EDGetTokenT< edm::View< reco::Muon > > src_
Labels for input collections.
const int layersToTry_
How many layers to try.
const int hitsToTry_
How many hits to try on same layer.
edm::EDGetTokenT< MeasurementTrackerEvent > measurementTrackerTag_
const std::string muonPropagatorName_
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().

277  {
279  std::unique_ptr<Propagator> pmuon_cloned = SetPropagationDirection(*muonPropagator_, alongMomentum);
280  for (unsigned int i = 0; i < tk.recHitsSize(); ++i) {
281  const TrackingRecHit *hit = &*tk.recHit(i);
282  const GeomDet *det = geometry_->idToDet(hit->geographicalId());
283  if (det == nullptr) continue;
284  if (i != 0) tsos = pmuon_cloned->propagate(tsos, det->surface());
285  if (!tsos.isValid()) continue;
286  LogDebug("OutsideInMuonSeeder") << " state " << i << " at x = " << tsos.globalPosition() << ", p = " << tsos.globalMomentum() << std::endl;
287  if (hit->isValid()) {
288  LogDebug("OutsideInMuonSeeder") << " valid rechit on detid " << hit->geographicalId().rawId() << std::endl;
289  } else {
290  LogDebug("OutsideInMuonSeeder") << " invalid rechit on detid " << hit->geographicalId().rawId() << std::endl;
291  }
292  }
293 }
#define LogDebug(id)
size_t recHitsSize() const
Get number of RecHits. (Warning, this includes invalid hits, which are not physical hits)...
Definition: Track.h:119
edm::ESHandle< MagneticField > magfield_
edm::ESHandle< GlobalTrackingGeometry > geometry_
GlobalPoint globalPosition() const
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:42
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
std::unique_ptr< Propagator > SetPropagationDirection(Propagator const &iprop, PropagationDirection dir)
edm::ESHandle< Propagator > muonPropagator_
bool isValid() const
GlobalVector globalMomentum() const
TrackingRecHitRef recHit(size_t i) const
Get i-th hit on the track.
Definition: Track.h:114
const GeomDet * idToDet(DetId) const override
DetId geographicalId() const
TrajectoryStateOnSurface innerStateOnSurface(const reco::Track &tk, const TrackingGeometry &geom, const MagneticField *field, bool withErr=true)
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().

231  {
232  TrajectoryStateOnSurface onLayer(state);
233  onLayer.rescaleError(errorRescaling_);
234  std::vector< GeometricSearchDet::DetWithState > dets;
235  layer.compatibleDetsV(onLayer, muon_propagator, *estimator_, dets);
236 
237  if (debug_) {
238  LogDebug("OutsideInMuonSeeder") << "Query on layer around x = " << onLayer.globalPosition() <<
239  " with local pos error " << sqrt(onLayer.localError().positionError().xx()) << " , " << sqrt(onLayer.localError().positionError().yy()) << " , " <<
240  " returned " << dets.size() << " compatible detectors" << std::endl;
241  }
242 
243  std::vector<TrajectoryMeasurement> meas;
244  for (std::vector<GeometricSearchDet::DetWithState>::const_iterator it = dets.begin(), ed = dets.end(); it != ed; ++it) {
245  MeasurementDetWithData det = measurementTracker.idToDet(it->first->geographicalId());
246  if (det.isNull()) { std::cerr << "BOGUS detid " << it->first->geographicalId().rawId() << std::endl; continue; }
247  if (!it->second.isValid()) continue;
248  std::vector < TrajectoryMeasurement > mymeas = det.fastMeasurements(it->second, state, tracker_propagator, *estimator_);
249  if (debug_) LogDebug("OutsideInMuonSeeder") << "Query on detector " << it->first->geographicalId().rawId() << " returned " << mymeas.size() << " measurements." << std::endl;
250  for (std::vector<TrajectoryMeasurement>::const_iterator it2 = mymeas.begin(), ed2 = mymeas.end(); it2 != ed2; ++it2) {
251  if (it2->recHit()->isValid()) meas.push_back(*it2);
252  }
253  }
254  int found = 0;
255  std::sort(meas.begin(), meas.end(), TrajMeasLessEstim());
256  for (std::vector<TrajectoryMeasurement>::const_iterator it2 = meas.begin(), ed2 = meas.end(); it2 != ed2; ++it2) {
257  if (debug_) {
258  LogDebug("OutsideInMuonSeeder") << " inspecting Hit with chi2 = " << it2->estimate() << std::endl;
259  LogDebug("OutsideInMuonSeeder") << " track state " << it2->forwardPredictedState().globalPosition() << std::endl;
260  LogDebug("OutsideInMuonSeeder") << " rechit position " << it2->recHit()->globalPosition() << std::endl;
261  }
262  TrajectoryStateOnSurface updated = updator_->update(it2->forwardPredictedState(), *it2->recHit());
263  if (updated.isValid()) {
264  if (debug_) LogDebug("OutsideInMuonSeeder") << " --> updated state: x = " << updated.globalPosition() << ", p = " << updated.globalMomentum() << std::endl;
266  seedHits.push_back(*it2->recHit()->hit());
267  PTrajectoryStateOnDet const & pstate = trajectoryStateTransform::persistentState(updated, it2->recHit()->geographicalId().rawId());
268  TrajectorySeed seed(pstate, std::move(seedHits), oppositeToMomentum);
269  out.push_back(seed);
270  found++; if (found == hitsToTry_) break;
271  }
272  }
273  return found;
274 }
#define LogDebug(id)
edm::ESHandle< Chi2MeasurementEstimatorBase > estimator_
PTrajectoryStateOnDet persistentState(const TrajectoryStateOnSurface &ts, unsigned int detid)
GlobalPoint globalPosition() const
edm::ESHandle< TrajectoryStateUpdator > updator_
virtual TrajectoryStateOnSurface update(const TrajectoryStateOnSurface &, const TrackingRecHit &) const =0
void push_back(D *&d)
Definition: OwnVector.h:290
const bool debug_
Dump deug information.
const double errorRescaling_
How much to rescale errors from STA.
T sqrt(T t)
Definition: SSEVec.h:18
std::vector< TrajectoryMeasurement > fastMeasurements(const TrajectoryStateOnSurface &stateOnThisDet, const TrajectoryStateOnSurface &tsos2, const Propagator &prop, const MeasurementEstimator &est) const
GlobalVector globalMomentum() const
const int hitsToTry_
How many hits to try on same layer.
virtual void compatibleDetsV(const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetWithState > &result) const
def move(src, dest)
Definition: eostools.py:510
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().

116  {
117  using namespace edm;
118  using namespace std;
119 
120  iSetup.get<IdealMagneticFieldRecord>().get(magfield_);
126 
127  Handle<MeasurementTrackerEvent> measurementTracker;
128  iEvent.getByToken(measurementTrackerTag_, measurementTracker);
129 
130  ESHandle<TrackerGeometry> tmpTkGeometry;
131  iSetup.get<TrackerDigiGeometryRecord>().get(tmpTkGeometry);
132 
134  iEvent.getByToken(src_, src);
135 
136 
137  auto out = std::make_unique<std::vector<TrajectorySeed>>();
138 
139  for (auto const & mu : *src) {
140  if (mu.outerTrack().isNull() || !selector_(mu)) continue;
141  if (debug_ && mu.innerTrack().isNonnull()) doDebug(*mu.innerTrack());
142 
143  // Better clone here and not directly into doLayer to avoid
144  // useless clone/destroy operations to set, in the end, the
145  // very same direction every single time.
146  std::unique_ptr<Propagator> pmuon_cloned = SetPropagationDirection(*muonPropagator_,
148  std::unique_ptr<Propagator> ptracker_cloned = SetPropagationDirection(*trackerPropagator_, alongMomentum);
149 
150  int sizeBefore = out->size();
151  if (debug_) LogDebug("OutsideInMuonSeeder") << "\n\n\nSeeding for muon of pt " << mu.pt() << ", eta " << mu.eta() << ", phi " << mu.phi() << std::endl;
152  const reco::Track &tk = *mu.outerTrack();
153 
157 
158  if (std::abs(tk.eta()) < maxEtaForTOB_) {
159  std::vector< BarrelDetLayer const* > const & tob = measurementTracker->geometricSearchTracker()->tobLayers();
160  int found = 0;
161  int iLayer = tob.size();
162  if(iLayer==0) LogError("OutsideInMuonSeeder") << "TOB has no layers." ;
163 
164  for (auto it = tob.rbegin(), ed = tob.rend(); it != ed; ++it, --iLayer) {
165  if (debug_) LogDebug("OutsideInMuonSeeder") << "\n ==== Trying TOB " << iLayer << " ====" << std::endl;
166  if (doLayer(**it, state, *out,
167  *(pmuon_cloned.get()),
168  *(ptracker_cloned.get()),
169  *measurementTracker)) {
170  if (++found == layersToTry_) break;
171  }
172  }
173  }
174  if (tk.eta() > minEtaForTEC_) {
175  const auto& forwLayers = tmpTkGeometry->isThere(GeomDetEnumerators::P2OTEC) ?
176  measurementTracker->geometricSearchTracker()->posTidLayers() : measurementTracker->geometricSearchTracker()->posTecLayers();
177  if (tmpTkGeometry->isThere(GeomDetEnumerators::P2OTEC)) {
178  LogDebug("OutsideInMuonSeeder") << "\n We are using the Phase2 Outer Tracker (defined as a TID+). ";
179  }
180  LogTrace("OutsideInMuonSeeder") << "\n ==== TEC+ tot layers " << forwLayers.size() << " ====" << std::endl;
181  int found = 0;
182  int iLayer = forwLayers.size();
183  if(iLayer==0) LogError("OutsideInMuonSeeder") << "TEC+ has no layers." ;
184 
185  if (debug_) LogDebug("OutsideInMuonSeeder") << "\n ==== Tot layers " << forwLayers.size() << " ====" << std::endl;
186  for (auto it = forwLayers.rbegin(), ed = forwLayers.rend(); it != ed; ++it, --iLayer) {
187  if (debug_) LogDebug("OutsideInMuonSeeder") << "\n ==== Trying Forward Layer +" << +iLayer << " ====" << std::endl;
188  if (doLayer(**it, state, *out,
189  *(pmuon_cloned.get()),
190  *(ptracker_cloned.get()),
191  *measurementTracker)) {
192  if (++found == layersToTry_) break;
193  }
194  }
195  }
196  if (tk.eta() < -minEtaForTEC_) {
197  const auto& forwLayers = tmpTkGeometry->isThere(GeomDetEnumerators::P2OTEC) ?
198  measurementTracker->geometricSearchTracker()->negTidLayers() : measurementTracker->geometricSearchTracker()->negTecLayers();
199  if (tmpTkGeometry->isThere(GeomDetEnumerators::P2OTEC)) {
200  LogDebug("OutsideInMuonSeeder") << "\n We are using the Phase2 Outer Tracker (defined as a TID-). ";
201  }
202  LogTrace("OutsideInMuonSeeder") << "\n ==== TEC- tot layers " << forwLayers.size() << " ====" << std::endl;
203  int found = 0;
204  int iLayer = forwLayers.size();
205  if(iLayer==0) LogError("OutsideInMuonSeeder") << "TEC- has no layers." ;
206 
207  if (debug_) LogDebug("OutsideInMuonSeeder") << "\n ==== Tot layers " << forwLayers.size() << " ====" << std::endl;
208  for (auto it = forwLayers.rbegin(), ed = forwLayers.rend(); it != ed; ++it, --iLayer) {
209  if (debug_) LogDebug("OutsideInMuonSeeder") << "\n ==== Trying Forward Layer -" << -iLayer << " ====" << std::endl;
210  if (doLayer(**it, state, *out,
211  *(pmuon_cloned.get()),
212  *(ptracker_cloned.get()),
213  *measurementTracker)) {
214  if (++found == layersToTry_) break;
215  }
216  }
217  }
218  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;
219 
220  }
221 
222  iEvent.put(std::move(out));
223 }
#define LogDebug(id)
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:127
const std::string estimatorName_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:508
edm::ESHandle< Chi2MeasurementEstimatorBase > estimator_
edm::ESHandle< MagneticField > magfield_
edm::ESHandle< GlobalTrackingGeometry > geometry_
edm::ESHandle< Propagator > trackerPropagator_
const bool fromVertex_
Do inside-out.
edm::ESHandle< TrajectoryStateUpdator > updator_
int doLayer(const GeometricSearchDet &layer, const TrajectoryStateOnSurface &state, std::vector< TrajectorySeed > &out, const Propagator &muon_propagator, const Propagator &tracker_propagator, const MeasurementTrackerEvent &mte) const
bool isThere(GeomDetEnumerators::SubDetector subdet) const
std::unique_ptr< Propagator > SetPropagationDirection(Propagator const &iprop, PropagationDirection dir)
const bool debug_
Dump deug information.
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:651
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
const std::string trackerPropagatorName_
const std::string updatorName_
const int mu
Definition: Constants.h:22
#define LogTrace(id)
edm::ESHandle< Propagator > muonPropagator_
const T & get() const
Definition: EventSetup.h:55
StringCutObjectSelector< reco::Muon > selector_
Muon selection.
edm::EDGetTokenT< edm::View< reco::Muon > > src_
Labels for input collections.
void doDebug(const reco::Track &tk) const
const int layersToTry_
How many layers to try.
HLT enums.
edm::EDGetTokenT< MeasurementTrackerEvent > measurementTrackerTag_
FreeTrajectoryState initialFreeState(const reco::Track &tk, const MagneticField *field, bool withErr=true)
T const * product() const
Definition: ESHandle.h:86
const std::string muonPropagatorName_
def move(src, dest)
Definition: eostools.py:510
TrajectoryStateOnSurface innerStateOnSurface(const reco::Track &tk, const TrackingGeometry &geom, const MagneticField *field, bool withErr=true)

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