Inheritance diagram for MkFitOutputConverter:

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

	~MkFitOutputConverter () override=default

Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::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)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
std::pair< TrajectoryStateOnSurface, const GeomDet * >	backwardFit (const FreeTrajectoryState &fts, const edm::OwnVector< TrackingRecHit > &hits, const Propagator &propagatorAlong, const Propagator &propagatorOpposite, const TkClonerImpl &hitCloner, bool lastHitWasInvalid, bool lastHitWasChanged) const

TrackCandidateCollection	convertCandidates (const MkFitOutputWrapper &mkFitOutput, const MkFitHitIndexMap &hitIndexMap, const edm::View< TrajectorySeed > &seeds, const TrackerGeometry &geom, const MagneticField &mf, const Propagator &propagatorAlong, const Propagator &propagatorOpposite, const TkClonerImpl &hitCloner, const std::vector< const DetLayer * > &detLayers, const mkfit::TrackVec &mkFitSeeds) const

std::pair< TrajectoryStateOnSurface, const GeomDet * >	convertInnermostState (const FreeTrajectoryState &fts, const edm::OwnVector< TrackingRecHit > &hits, const Propagator &propagatorAlong, const Propagator &propagatorOpposite) const

std::vector< const DetLayer * >	createDetLayers (const mkfit::LayerNumberConverter &lnc, const GeometricSearchTracker &tracker, const TrackerTopology &ttopo) const

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

Private Attributes
bool	backwardFitInCMSSW_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_

edm::EDGetTokenT< MkFitInputWrapper >	hitsSeedsToken_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	mfToken_

edm::EDGetTokenT< MeasurementTrackerEvent >	mteToken_

std::string	propagatorAlongName_

edm::ESGetToken< Propagator, TrackingComponentsRecord >	propagatorAlongToken_

std::string	propagatorOppositeName_

edm::ESGetToken< Propagator, TrackingComponentsRecord >	propagatorOppositeToken_

edm::EDPutTokenT< std::vector< SeedStopInfo > >	putSeedStopInfoToken_

edm::EDPutTokenT< TrackCandidateCollection >	putTrackCandidateToken_

edm::EDGetTokenT< edm::View< TrajectorySeed > >	seedToken_

edm::EDGetTokenT< MkFitOutputWrapper >	tracksToken_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	ttopoToken_

std::string	ttrhBuilderName_

edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	ttrhBuilderToken_

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	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

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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 61 of file MkFitOutputConverter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 117 of file MkFitOutputConverter.cc.

References backwardFitInCMSSW_, geomToken_, edm::ParameterSet::getParameter(), mfToken_, mteToken_, propagatorAlongToken_, propagatorOppositeToken_, putSeedStopInfoToken_, putTrackCandidateToken_, seedToken_, tracksToken_, ttopoToken_, and ttrhBuilderToken_.

     : hitsSeedsToken_{consumes<MkFitInputWrapper>(iConfig.getParameter<edm::InputTag>("hitsSeeds"))},
       tracksToken_{consumes<MkFitOutputWrapper>(iConfig.getParameter<edm::InputTag>("tracks"))},
       seedToken_{consumes<edm::View<TrajectorySeed>>(iConfig.getParameter<edm::InputTag>("seeds"))},
       mteToken_{consumes<MeasurementTrackerEvent>(iConfig.getParameter<edm::InputTag>("measurementTrackerEvent"))},
       geomToken_{esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>()},
       propagatorAlongToken_{
           esConsumes<Propagator, TrackingComponentsRecord>(iConfig.getParameter<edm::ESInputTag>("propagatorAlong"))},
       propagatorOppositeToken_{esConsumes<Propagator, TrackingComponentsRecord>(
           iConfig.getParameter<edm::ESInputTag>("propagatorOpposite"))},
       ttopoToken_{esConsumes<TrackerTopology, TrackerTopologyRcd>()},
       mfToken_{esConsumes<MagneticField, IdealMagneticFieldRecord>()},
       ttrhBuilderToken_{esConsumes<TransientTrackingRecHitBuilder, TransientRecHitRecord>(
           iConfig.getParameter<edm::ESInputTag>("ttrhBuilder"))},
       putTrackCandidateToken_{produces<TrackCandidateCollection>()},
       putSeedStopInfoToken_{produces<std::vector<SeedStopInfo>>()},
       backwardFitInCMSSW_{iConfig.getParameter<bool>("backwardFitInCMSSW")} {}

MkFitOutputConverter::~MkFitOutputConverter ( )

overridedefault

Member Function Documentation

std::pair< TrajectoryStateOnSurface, const GeomDet * > MkFitOutputConverter::backwardFit	(	const FreeTrajectoryState &	fts,
		const edm::OwnVector< TrackingRecHit > &	hits,
		const Propagator &	propagatorAlong,
		const Propagator &	propagatorOpposite,
		const TkClonerImpl &	hitCloner,
		bool	lastHitWasInvalid,
		bool	lastHitWasChanged
	)		const

private

Definition at line 380 of file MkFitOutputConverter.cc.

References Chi2MeasurementEstimator_cfi::Chi2MeasurementEstimator, mps_fire::i, PixelPluginsPhase0_cfi::isBarrel, Trajectory::isValid(), Trajectory::lastMeasurement(), LogDebug, LogTrace, Propagator::magneticField(), oppositeToMomentum, FreeTrajectoryState::position(), Propagator::propagateWithPath(), HLT_2018_cff::propagatorAlong, HLT_2018_cff::propagatorOpposite, TrajectoryStateOnSurface::rescaleError(), edm::OwnVector< T, P >::size(), TrajectoryFitter::standard, and std::swap().

Referenced by convertCandidates().

                                   {
   // First filter valid hits as in TransientInitialStateEstimator
   TransientTrackingRecHit::ConstRecHitContainer firstHits;
 
   for (int i = hits.size() - 1; i >= 0; --i) {
     if (hits[i].det()) {
       // TransientTrackingRecHit::ConstRecHitContainer has shared_ptr,
       // and it is passed to backFitter below so it is really needed
       // to keep the interface. Since we keep the ownership in hits,
       // let's disable the deleter.
       firstHits.emplace_back(&(hits[i]), edm::do_nothing_deleter{});
     }
   }
 
   // Then propagate along to the surface of the last hit to get a TSOS
   const auto& lastHitSurface = firstHits.front()->det()->surface();
 
   const Propagator* tryFirst = &propagatorAlong;
   const Propagator* trySecond = &propagatorOpposite;
   if (lastHitWasInvalid || lastHitWasChanged) {
     LogTrace("MkFitOutputConverter") << "Propagating first opposite, then along, because lastHitWasInvalid? "
                                      << lastHitWasInvalid << " or lastHitWasChanged? " << lastHitWasChanged;
     std::swap(tryFirst, trySecond);
   } else {
     const auto lastHitSubdet = firstHits.front()->geographicalId().subdetId();
     const auto& surfacePos = lastHitSurface.position();
     const auto& lastHitPos = firstHits.front()->globalPosition();
     bool doSwitch = false;
     if (isBarrel(lastHitSubdet)) {
       doSwitch = (surfacePos.perp2() < lastHitPos.perp2());
     } else {
       doSwitch = (surfacePos.z() < lastHitPos.z());
     }
     if (doSwitch) {
       LogTrace("MkFitOutputConverter")
           << "Propagating first opposite, then along, because surface is inner than the hit; surface perp2 "
           << surfacePos.perp() << " hit " << lastHitPos.perp2() << " surface z " << surfacePos.z() << " hit "
           << lastHitPos.z();
 
       std::swap(tryFirst, trySecond);
     }
   }
 
   auto tsosDouble = tryFirst->propagateWithPath(fts, lastHitSurface);
   if (!tsosDouble.first.isValid()) {
     LogDebug("MkFitOutputConverter") << "Propagating to startingState failed, trying in another direction next";
     tsosDouble = trySecond->propagateWithPath(fts, lastHitSurface);
   }
   auto& startingState = tsosDouble.first;
 
   if (!startingState.isValid()) {
     edm::LogWarning("MkFitOutputConverter")
         << "startingState is not valid, FTS was\n"
         << fts << " last hit surface surface:"
         << "\n position " << lastHitSurface.position() << "\n phiSpan " << lastHitSurface.phiSpan().first << ","
         << lastHitSurface.phiSpan().first << "\n rSpan " << lastHitSurface.rSpan().first << ","
         << lastHitSurface.rSpan().first << "\n zSpan " << lastHitSurface.zSpan().first << ","
         << lastHitSurface.zSpan().first;
     return std::pair<TrajectoryStateOnSurface, const GeomDet*>();
   }
 
   // Then return back to the logic from TransientInitialStateEstimator
   startingState.rescaleError(100.);
 
   // avoid cloning
   KFUpdator const aKFUpdator;
   Chi2MeasurementEstimator const aChi2MeasurementEstimator(100., 3);
   KFTrajectoryFitter backFitter(
       &propagatorAlong, &aKFUpdator, &aChi2MeasurementEstimator, firstHits.size(), nullptr, &hitCloner);
 
   // assume for now that the propagation in mkfit always alongMomentum
   PropagationDirection backFitDirection = oppositeToMomentum;
 
   // only direction matters in this context
   TrajectorySeed fakeSeed(PTrajectoryStateOnDet(), edm::OwnVector<TrackingRecHit>(), backFitDirection);
 
   // ignore loopers for now
   Trajectory fitres = backFitter.fitOne(fakeSeed, firstHits, startingState, TrajectoryFitter::standard);
 
   LogDebug("MkFitOutputConverter") << "using a backward fit of :" << firstHits.size() << " hits, starting from:\n"
                                    << startingState << " to get the estimate of the initial state of the track.";
 
   if (!fitres.isValid()) {
     edm::LogWarning("MkFitOutputConverter") << "FitTester: first hits fit failed";
     return std::pair<TrajectoryStateOnSurface, const GeomDet*>();
   }
 
   TrajectoryMeasurement const& firstMeas = fitres.lastMeasurement();
 
   // magnetic field can be different!
   TrajectoryStateOnSurface firstState(firstMeas.updatedState().localParameters(),
                                       firstMeas.updatedState().localError(),
                                       firstMeas.updatedState().surface(),
                                       propagatorAlong.magneticField());
 
   firstState.rescaleError(100.);
 
   LogDebug("MkFitOutputConverter") << "the initial state is found to be:\n:" << firstState
                                    << "\n it's field pointer is: " << firstState.magneticField()
                                    << "\n the pointer from the state of the back fit was: "
                                    << firstMeas.updatedState().magneticField();
 
   return std::make_pair(firstState, firstMeas.recHit()->det());
 }

TrackCandidateCollection MkFitOutputConverter::convertCandidates	(	const MkFitOutputWrapper &	mkFitOutput,
		const MkFitHitIndexMap &	hitIndexMap,
		const edm::View< TrajectorySeed > &	seeds,
		const TrackerGeometry &	geom,
		const MagneticField &	mf,
		const Propagator &	propagatorAlong,
		const Propagator &	propagatorOpposite,
		const TkClonerImpl &	hitCloner,
		const std::vector< const DetLayer * > &	detLayers,
		const mkfit::TrackVec &	mkFitSeeds
	)		const

private

Definition at line 222 of file MkFitOutputConverter.cc.

References a, funct::abs(), edm::View< T >::at(), b, edm::OwnVector< T, P >::back(), backwardFit(), backwardFitInCMSSW_, HLT_2018_cff::candidates, MkFitOutputWrapper::candidateTracks(), TrackingRecHit::clone(), MkFitHitIndexMap::clusterIndex(), convertInnermostState(), runTheMatrix::err, Exception, MkFitOutputWrapper::fitTracks(), TrackingRecHit::geographicalId(), TrackingRecHit::globalPosition(), MkFitHitIndexMap::hitPtr(), mps_fire::i, MkFitHitIndexMap::MkFitHit::index(), PixelPluginsPhase0_cfi::isBarrel, dqmiolumiharvest::j, LogTrace, PV3DBase< T, PVType, FrameType >::mag2(), nhits, convertSQLitetoXML_cfg::output, trajectoryStateTransform::persistentState(), edm::OwnVector< T, P >::push_back(), DetId::rawId(), FastTrackerRecHitMaskProducer_cfi::recHits, edm::View< T >::refAt(), edm::OwnVector< T, P >::sort(), UNINITIALIZED, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by produce().

                                                                                                         {
   TrackCandidateCollection output;
   const auto& candidates = backwardFitInCMSSW_ ? mkFitOutput.candidateTracks() : mkFitOutput.fitTracks();
   output.reserve(candidates.size());
 
   LogTrace("MkFitOutputConverter") << "Number of candidates " << mkFitOutput.candidateTracks().size();
 
   int candIndex = -1;
   for (const auto& cand : candidates) {
     ++candIndex;
     LogTrace("MkFitOutputConverter") << "Candidate " << candIndex << " pT " << cand.pT() << " eta " << cand.momEta()
                                      << " phi " << cand.momPhi() << " chi2 " << cand.chi2();
 
     // hits
     edm::OwnVector<TrackingRecHit> recHits;
     // nTotalHits() gives sum of valid hits (nFoundHits()) and
     // invalid/missing hits (up to a maximum of 32 inside mkFit,
     // restriction to be lifted in the future)
     const int nhits = cand.nTotalHits();
     bool lastHitInvalid = false;
     for (int i = 0; i < nhits; ++i) {
       const auto& hitOnTrack = cand.getHitOnTrack(i);
       LogTrace("MkFitOutputConverter") << " hit on layer " << hitOnTrack.layer << " index " << hitOnTrack.index;
       if (hitOnTrack.index < 0) {
         // See index-desc.txt file in mkFit for description of negative values
         //
         // In order to use the regular InvalidTrackingRecHit I'd need
         // a GeomDet (and "unfortunately" that is needed in
         // TrackProducer).
         //
         // I guess we could take the track state and propagate it to
         // each layer to find the actual module the track crosses, and
         // check whether it is active or not to be able to mark
         // inactive hits
         const auto* detLayer = detLayers.at(hitOnTrack.layer);
         if (detLayer == nullptr) {
           throw cms::Exception("LogicError") << "DetLayer for layer index " << hitOnTrack.layer << " is null!";
         }
         // In principle an InvalidTrackingRecHitNoDet could be
         // inserted here, but it seems that it is best to deal with
         // them in the TrackProducer.
         lastHitInvalid = true;
       } else {
         recHits.push_back(hitIndexMap.hitPtr(MkFitHitIndexMap::MkFitHit{hitOnTrack.index, hitOnTrack.layer})->clone());
         LogTrace("MkFitOutputConverter") << "  pos " << recHits.back().globalPosition().x() << " "
                                          << recHits.back().globalPosition().y() << " "
                                          << recHits.back().globalPosition().z() << " mag2 "
                                          << recHits.back().globalPosition().mag2() << " detid "
                                          << recHits.back().geographicalId().rawId() << " cluster "
                                          << hitIndexMap.clusterIndex(
                                                 MkFitHitIndexMap::MkFitHit{hitOnTrack.index, hitOnTrack.layer});
         lastHitInvalid = false;
       }
     }
 
     const auto lastHitId = recHits.back().geographicalId();
 
     // MkFit hits are *not* in the order of propagation, sort by 3D radius for now (as we don't have loopers)
     // TODO: Improve the sorting (extract keys? maybe even bubble sort would work well as the hits are almost in the correct order)
     recHits.sort([](const auto& a, const auto& b) {
       const auto asub = a.geographicalId().subdetId();
       const auto bsub = b.geographicalId().subdetId();
       if (asub != bsub) {
         // Subdetector order (BPix, FPix, TIB, TID, TOB, TEC) corresponds also the navigation
         return asub < bsub;
       }
 
       const auto& apos = a.globalPosition();
       const auto& bpos = b.globalPosition();
 
       if (isBarrel(asub)) {
         return apos.perp2() < bpos.perp2();
       }
       return std::abs(apos.z()) < std::abs(bpos.z());
     });
 
     const bool lastHitChanged = (recHits.back().geographicalId() != lastHitId);  // TODO: make use of the bools
 
     // seed
     const auto seedIndex = cand.label();
     LogTrace("MkFitOutputConverter") << " from seed " << seedIndex << " seed hits";
     const auto& mkseed = mkFitSeeds.at(cand.label());
     for (int i = 0; i < mkseed.nTotalHits(); ++i) {
       const auto& hitOnTrack = mkseed.getHitOnTrack(i);
       LogTrace("MkFitOutputConverter") << "  hit on layer " << hitOnTrack.layer << " index " << hitOnTrack.index;
       // sanity check for now
       const auto& candHitOnTrack = cand.getHitOnTrack(i);
       if (hitOnTrack.layer != candHitOnTrack.layer) {
         throw cms::Exception("LogicError")
             << "Candidate " << candIndex << " from seed " << seedIndex << " hit " << i
             << " has different layer in candidate (" << candHitOnTrack.layer << ") and seed (" << hitOnTrack.layer
             << ")."
             << " Hit indices are " << candHitOnTrack.index << " and " << hitOnTrack.index << ", respectively";
       }
       if (hitOnTrack.index != candHitOnTrack.index) {
         throw cms::Exception("LogicError") << "Candidate " << candIndex << " from seed " << seedIndex << " hit " << i
                                            << " has different hit index in candidate (" << candHitOnTrack.index
                                            << ") and seed (" << hitOnTrack.index << ") on layer " << hitOnTrack.layer;
       }
     }
 
     // state
     auto state = cand.state();  // copy because have to modify
     state.convertFromCCSToCartesian();
     const auto& param = state.parameters;
     const auto& err = state.errors;
     AlgebraicSymMatrix66 cov;
     for (int i = 0; i < 6; ++i) {
       for (int j = i; j < 6; ++j) {
         cov[i][j] = err.At(i, j);
       }
     }
 
     auto fts = FreeTrajectoryState(
         GlobalTrajectoryParameters(
             GlobalPoint(param[0], param[1], param[2]), GlobalVector(param[3], param[4], param[5]), state.charge, &mf),
         CartesianTrajectoryError(cov));
     if (!fts.curvilinearError().posDef()) {
       edm::LogWarning("MkFitOutputConverter") << "Curvilinear error not pos-def\n"
                                               << fts.curvilinearError().matrix() << "\noriginal 6x6 covariance matrix\n"
                                               << cov << "\ncandidate ignored";
       continue;
     }
 
     auto tsosDet =
         backwardFitInCMSSW_
             ? backwardFit(fts, recHits, propagatorAlong, propagatorOpposite, hitCloner, lastHitInvalid, lastHitChanged)
             : convertInnermostState(fts, recHits, propagatorAlong, propagatorOpposite);
     if (!tsosDet.first.isValid()) {
       edm::LogWarning("MkFitOutputConverter")
           << "Backward fit of candidate " << candIndex << " failed, ignoring the candidate";
       continue;
     }
 
     // convert to persistent, from CkfTrackCandidateMakerBase
     auto pstate = trajectoryStateTransform::persistentState(tsosDet.first, tsosDet.second->geographicalId().rawId());
 
     output.emplace_back(
         recHits,
         seeds.at(seedIndex),
         pstate,
         seeds.refAt(seedIndex),
         0,                                               // mkFit does not produce loopers, so set nLoops=0
         static_cast<uint8_t>(StopReason::UNINITIALIZED)  // TODO: ignore details of stopping reason as well for now
     );
   }
   return output;
 }

std::pair< TrajectoryStateOnSurface, const GeomDet * > MkFitOutputConverter::convertInnermostState	(	const FreeTrajectoryState &	fts,
		const edm::OwnVector< TrackingRecHit > &	hits,
		const Propagator &	propagatorAlong,
		const Propagator &	propagatorOpposite
	)		const

private

Definition at line 492 of file MkFitOutputConverter.cc.

References DEFINE_FWK_MODULE, Exception, LogDebug, and Propagator::propagateWithPath().

Referenced by convertCandidates().

                                                 {
   auto det = hits[0].det();
   if (det == nullptr) {
     throw cms::Exception("LogicError") << "Got nullptr from the first hit det()";
   }
 
   const auto& firstHitSurface = det->surface();
 
   auto tsosDouble = propagatorAlong.propagateWithPath(fts, firstHitSurface);
   if (!tsosDouble.first.isValid()) {
     LogDebug("MkFitOutputConverter") << "Propagating to startingState along momentum failed, trying opposite next";
     tsosDouble = propagatorOpposite.propagateWithPath(fts, firstHitSurface);
   }
 
   return std::make_pair(tsosDouble.first, det);
 }

std::vector< const DetLayer * > MkFitOutputConverter::createDetLayers	(	const mkfit::LayerNumberConverter &	lnc,
		const GeometricSearchTracker &	tracker,
		const TrackerTopology &	ttopo
	)		const

private

Definition at line 181 of file MkFitOutputConverter.cc.

References GeometricSearchTracker::allLayers(), AlCaHLTBitMon_QueryRunRegistry::comp, constexpr, Exception, TrackerTopology::layer(), or, PosEndcap, TrackerTopology::side(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, and UNLIKELY.

Referenced by produce().

                                                                                                        {
   std::vector<const DetLayer*> dets(lnc.nLayers(), nullptr);
 
   auto isPlusSide = [&ttopo](const DetId& detid) {
     return ttopo.side(detid) == static_cast<unsigned>(TrackerDetSide::PosEndcap);
   };
   auto setDet = [&lnc, &dets, &isPlusSide](
                     const int subdet, const int layer, const int isStereo, const DetId& detId, const DetLayer* lay) {
     const int index = lnc.convertLayerNumber(subdet, layer, false, isStereo, isPlusSide(detId));
     if (index < 0 or static_cast<unsigned>(index) >= dets.size()) {
       throw cms::Exception("LogicError") << "Invalid mkFit layer index " << index << " for DetId " << detId
                                          << " subdet " << subdet << " layer " << layer << " isStereo " << isStereo;
     }
     dets[index] = lay;
   };
   constexpr int monoLayer = 0;
   constexpr int stereoLayer = 1;
   for (const DetLayer* lay : tracker.allLayers()) {
     const auto& comp = lay->basicComponents();
     if (UNLIKELY(comp.empty())) {
       throw cms::Exception("LogicError") << "Got a tracker layer (subdet " << lay->subDetector()
                                          << ") with empty basicComponents.";
     }
     // First component is enough for layer and side information
     const auto& detId = comp.front()->geographicalId();
     const auto subdet = detId.subdetId();
     const auto layer = ttopo.layer(detId);
 
     // TODO: mono/stereo structure is still hardcoded for phase0/1 strip tracker
     setDet(subdet, layer, monoLayer, detId, lay);
     if (((subdet == StripSubdetector::TIB or subdet == StripSubdetector::TOB) and (layer == 1 or layer == 2)) or
         subdet == StripSubdetector::TID or subdet == StripSubdetector::TEC) {
       setDet(subdet, layer, stereoLayer, detId, lay);
     }
   }
 
   return dets;
 }

void MkFitOutputConverter::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 135 of file MkFitOutputConverter.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), and edm::ParameterDescriptionNode::setComment().

                                                                                       {
   edm::ParameterSetDescription desc;
 
   desc.add("hitsSeeds", edm::InputTag{"mkFitInputConverter"});
   desc.add("tracks", edm::InputTag{"mkFitProducer"});
   desc.add("seeds", edm::InputTag{"initialStepSeeds"});
   desc.add("measurementTrackerEvent", edm::InputTag{"MeasurementTrackerEvent"});
   desc.add("ttrhBuilder", edm::ESInputTag{"", "WithTrackAngle"});
   desc.add("propagatorAlong", edm::ESInputTag{"", "PropagatorWithMaterial"});
   desc.add("propagatorOpposite", edm::ESInputTag{"", "PropagatorWithMaterialOpposite"});
   desc.add("backwardFitInCMSSW", false)
       ->setComment("Do backward fit (to innermost hit) in CMSSW (true) or mkFit (false)");
 
   descriptions.addWithDefaultLabel(desc);
 }

void MkFitOutputConverter::produce	(	edm::StreamID	iID,
		edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		const

overrideprivate

Definition at line 151 of file MkFitOutputConverter.cc.

References watchdog::const, convertCandidates(), createDetLayers(), edm::Event::emplace(), Exception, geomToken_, edm::Event::get(), edm::EventSetup::getData(), hitsSeedsToken_, mfToken_, mteToken_, propagatorAlongToken_, propagatorOppositeToken_, putSeedStopInfoToken_, putTrackCandidateToken_, seedToken_, tracksToken_, ttopoToken_, and ttrhBuilderToken_.

                                                                                                        {
   const auto& seeds = iEvent.get(seedToken_);
   const auto& hitsSeeds = iEvent.get(hitsSeedsToken_);
   const auto& mte = iEvent.get(mteToken_);
 
   const auto& ttrhBuilder = iSetup.getData(ttrhBuilderToken_);
   const auto* tkBuilder = dynamic_cast<TkTransientTrackingRecHitBuilder const*>(&ttrhBuilder);
   if (!tkBuilder) {
     throw cms::Exception("LogicError") << "TTRHBuilder must be of type TkTransientTrackingRecHitBuilder";
   }
 
   // Convert mkfit presentation back to CMSSW
   const auto detlayers =
       createDetLayers(hitsSeeds.layerNumberConverter(), *(mte.geometricSearchTracker()), iSetup.getData(ttopoToken_));
   iEvent.emplace(putTrackCandidateToken_,
                  convertCandidates(iEvent.get(tracksToken_),
                                    hitsSeeds.hitIndexMap(),
                                    seeds,
                                    iSetup.getData(geomToken_),
                                    iSetup.getData(mfToken_),
                                    iSetup.getData(propagatorAlongToken_),
                                    iSetup.getData(propagatorOppositeToken_),
                                    tkBuilder->cloner(),
                                    detlayers,
                                    hitsSeeds.seeds()));
 
   // TODO: SeedStopInfo is currently unfilled
   iEvent.emplace(putSeedStopInfoToken_, seeds.size());
 }