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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

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 mkfit::EventOfHits &eventOfHits, const MkFitClusterIndexToHit &pixelClusterIndexToHit, const MkFitClusterIndexToHit &stripClusterIndexToHit, const edm::View< TrajectorySeed > &seeds, 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

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

Private Attributes
const edm::EDGetTokenT< MkFitEventOfHits >	eventOfHitsToken_

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	mfToken_

const edm::ESGetToken< MkFitGeometry, TrackerRecoGeometryRecord >	mkFitGeomToken_

const edm::EDGetTokenT< MkFitSeedWrapper >	mkfitSeedToken_

const edm::EDGetTokenT< MkFitClusterIndexToHit >	pixelClusterIndexToHitToken_

const std::string	propagatorAlongName_

const edm::ESGetToken< Propagator, TrackingComponentsRecord >	propagatorAlongToken_

const std::string	propagatorOppositeName_

const edm::ESGetToken< Propagator, TrackingComponentsRecord >	propagatorOppositeToken_

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

const edm::EDPutTokenT< TrackCandidateCollection >	putTrackCandidateToken_

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

const edm::EDGetTokenT< MkFitClusterIndexToHit >	stripClusterIndexToHitToken_

const edm::EDGetTokenT< MkFitOutputWrapper >	tracksToken_

const std::string	ttrhBuilderName_

const 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
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

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

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

Detailed Description

Definition at line 58 of file MkFitOutputConverter.cc.

Constructor & Destructor Documentation

◆ MkFitOutputConverter()

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

explicit

Definition at line 111 of file MkFitOutputConverter.cc.

     : eventOfHitsToken_{consumes<MkFitEventOfHits>(iConfig.getParameter<edm::InputTag>("mkFitEventOfHits"))},
       pixelClusterIndexToHitToken_{consumes(iConfig.getParameter<edm::InputTag>("mkFitPixelHits"))},
       stripClusterIndexToHitToken_{consumes(iConfig.getParameter<edm::InputTag>("mkFitStripHits"))},
       mkfitSeedToken_{consumes<MkFitSeedWrapper>(iConfig.getParameter<edm::InputTag>("mkFitSeeds"))},
       tracksToken_{consumes<MkFitOutputWrapper>(iConfig.getParameter<edm::InputTag>("tracks"))},
       seedToken_{consumes<edm::View<TrajectorySeed>>(iConfig.getParameter<edm::InputTag>("seeds"))},
       propagatorAlongToken_{
           esConsumes<Propagator, TrackingComponentsRecord>(iConfig.getParameter<edm::ESInputTag>("propagatorAlong"))},
       propagatorOppositeToken_{esConsumes<Propagator, TrackingComponentsRecord>(
           iConfig.getParameter<edm::ESInputTag>("propagatorOpposite"))},
       mfToken_{esConsumes<MagneticField, IdealMagneticFieldRecord>()},
       ttrhBuilderToken_{esConsumes<TransientTrackingRecHitBuilder, TransientRecHitRecord>(
           iConfig.getParameter<edm::ESInputTag>("ttrhBuilder"))},
       mkFitGeomToken_{esConsumes<MkFitGeometry, TrackerRecoGeometryRecord>()},
       putTrackCandidateToken_{produces<TrackCandidateCollection>()},
       putSeedStopInfoToken_{produces<std::vector<SeedStopInfo>>()} {}

References edm::ParameterSet::getParameter().

◆ ~MkFitOutputConverter()

MkFitOutputConverter::~MkFitOutputConverter ( )

overridedefault

Member Function Documentation

◆ backwardFit()

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

                                   {
   // 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());
 }

References hfClusterShapes_cfi::hits, mps_fire::i, PixelPluginsPhase0_cfi::isBarrel, Trajectory::isValid(), Trajectory::lastMeasurement(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), LogDebug, LogTrace, TrajectoryStateOnSurface::magneticField(), oppositeToMomentum, FreeTrajectoryState::position(), Propagator::propagateWithPath(), HLT_FULL_cff::propagatorAlong, HLT_FULL_cff::propagatorOpposite, TrajectoryMeasurement::recHit(), TrajectoryStateOnSurface::rescaleError(), TrajectoryFitter::standard, TrajectoryStateOnSurface::surface(), std::swap(), and TrajectoryMeasurement::updatedState().

Referenced by convertCandidates().

◆ convertCandidates()

TrackCandidateCollection MkFitOutputConverter::convertCandidates	(	const MkFitOutputWrapper &	mkFitOutput,
		const mkfit::EventOfHits &	eventOfHits,
		const MkFitClusterIndexToHit &	pixelClusterIndexToHit,
		const MkFitClusterIndexToHit &	stripClusterIndexToHit,
		const edm::View< TrajectorySeed > &	seeds,
		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 174 of file MkFitOutputConverter.cc.

                                                                                                         {
   TrackCandidateCollection output;
   const auto& candidates = mkFitOutput.tracks();
   output.reserve(candidates.size());
  
   LogTrace("MkFitOutputConverter") << "Number of candidates " << candidates.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 {
         auto const isPixel = eventOfHits[hitOnTrack.layer].is_pix_lyr();
         auto const& hits = isPixel ? pixelClusterIndexToHit.hits() : stripClusterIndexToHit.hits();
  
         auto const& thit = static_cast<BaseTrackerRecHit const&>(*hits[hitOnTrack.index]);
         if (thit.firstClusterRef().isPixel() || thit.detUnit()->type().isEndcap()) {
           recHits.push_back(hits[hitOnTrack.index]->clone());
         } else {
           recHits.push_back(std::make_unique<SiStripRecHit1D>(
               thit.localPosition(),
               LocalError(thit.localPositionError().xx(), 0.f, std::numeric_limits<float>::max()),
               *thit.det(),
               thit.firstClusterRef()));
         }
         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 " << hitOnTrack.index;
         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";
  
     // state
     auto state = cand.state();  // copy because have to modify
     state.convertFromCCSToGlbCurvilinear();
     const auto& param = state.parameters;
     const auto& err = state.errors;
     AlgebraicSymMatrix55 cov;
     for (int i = 0; i < 5; ++i) {
       for (int j = i; j < 5; ++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),
         CurvilinearTrajectoryError(cov));
     if (!fts.curvilinearError().posDef()) {
       edm::LogWarning("MkFitOutputConverter") << "Curvilinear error not pos-def\n"
                                               << fts.curvilinearError().matrix() << "\ncandidate ignored";
       continue;
     }
  
     auto tsosDet =
         mkFitOutput.propagatedToFirstLayer()
             ? convertInnermostState(fts, recHits, propagatorAlong, propagatorOpposite)
             : backwardFit(fts, recHits, propagatorAlong, propagatorOpposite, hitCloner, lastHitInvalid, lastHitChanged);
     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;
 }

References a, funct::abs(), b, backwardFit(), HLT_FULL_cff::candidates, convertInnermostState(), submitPVResolutionJobs::err, InitialStepPreSplitting_cff::eventOfHits, Exception, hfClusterShapes_cfi::hits, MkFitClusterIndexToHit::hits(), mps_fire::i, PixelPluginsPhase0_cfi::isBarrel, fastTrackerRecHitType::isPixel(), dqmiolumiharvest::j, LogTrace, SiStripPI::max, nhits, convertSQLitetoXML_cfg::output, trajectoryStateTransform::persistentState(), MkFitOutputWrapper::propagatedToFirstLayer(), HLT_FULL_cff::propagatorAlong, HLT_FULL_cff::propagatorOpposite, FastTrackerRecHitMaskProducer_cfi::recHits, DetachedQuadStep_cff::seeds, MkFitOutputWrapper::tracks(), and UNINITIALIZED.

Referenced by produce().

◆ convertInnermostState()

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

                                                 {
   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);
 }

References Exception, hfClusterShapes_cfi::hits, LogDebug, HLT_FULL_cff::propagatorAlong, and HLT_FULL_cff::propagatorOpposite.

Referenced by convertCandidates().

◆ fillDescriptions()

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

static

Definition at line 129 of file MkFitOutputConverter.cc.

                                                                                       {
   edm::ParameterSetDescription desc;
  
   desc.add("mkFitEventOfHits", edm::InputTag{"mkFitEventOfHits"});
   desc.add("mkFitPixelHits", edm::InputTag{"mkFitSiPixelHits"});
   desc.add("mkFitStripHits", edm::InputTag{"mkFitSiStripHits"});
   desc.add("mkFitSeeds", edm::InputTag{"mkFitSeedConverter"});
   desc.add("tracks", edm::InputTag{"mkFitProducer"});
   desc.add("seeds", edm::InputTag{"initialStepSeeds"});
   desc.add("ttrhBuilder", edm::ESInputTag{"", "WithTrackAngle"});
   desc.add("propagatorAlong", edm::ESInputTag{"", "PropagatorWithMaterial"});
   desc.add("propagatorOpposite", edm::ESInputTag{"", "PropagatorWithMaterialOpposite"});
  
   descriptions.addWithDefaultLabel(desc);
 }

References edm::ConfigurationDescriptions::addWithDefaultLabel(), and submitPVResolutionJobs::desc.

◆ produce()

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

overrideprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 145 of file MkFitOutputConverter.cc.

                                                                                                        {
   const auto& seeds = iEvent.get(seedToken_);
   const auto& mkfitSeeds = iEvent.get(mkfitSeedToken_);
  
   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";
   }
   const auto& mkFitGeom = iSetup.getData(mkFitGeomToken_);
  
   // Convert mkfit presentation back to CMSSW
   iEvent.emplace(putTrackCandidateToken_,
                  convertCandidates(iEvent.get(tracksToken_),
                                    iEvent.get(eventOfHitsToken_).get(),
                                    iEvent.get(pixelClusterIndexToHitToken_),
                                    iEvent.get(stripClusterIndexToHitToken_),
                                    seeds,
                                    iSetup.getData(mfToken_),
                                    iSetup.getData(propagatorAlongToken_),
                                    iSetup.getData(propagatorOppositeToken_),
                                    tkBuilder->cloner(),
                                    mkFitGeom.detLayers(),
                                    mkfitSeeds.seeds()));
  
   // TODO: SeedStopInfo is currently unfilled
   iEvent.emplace(putSeedStopInfoToken_, seeds.size());
 }

References convertCandidates(), eventOfHitsToken_, Exception, edm::EventSetup::getData(), iEvent, mfToken_, mkFitGeomToken_, mkfitSeedToken_, pixelClusterIndexToHitToken_, propagatorAlongToken_, propagatorOppositeToken_, putSeedStopInfoToken_, putTrackCandidateToken_, DetachedQuadStep_cff::seeds, seedToken_, stripClusterIndexToHitToken_, tracksToken_, and ttrhBuilderToken_.

Member Data Documentation

◆ eventOfHitsToken_

const edm::EDGetTokenT<MkFitEventOfHits> MkFitOutputConverter::eventOfHitsToken_

private

Definition at line 93 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ mfToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> MkFitOutputConverter::mfToken_

private

Definition at line 101 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ mkFitGeomToken_

const edm::ESGetToken<MkFitGeometry, TrackerRecoGeometryRecord> MkFitOutputConverter::mkFitGeomToken_

private

Definition at line 103 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ mkfitSeedToken_

const edm::EDGetTokenT<MkFitSeedWrapper> MkFitOutputConverter::mkfitSeedToken_

private

Definition at line 96 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ pixelClusterIndexToHitToken_

const edm::EDGetTokenT<MkFitClusterIndexToHit> MkFitOutputConverter::pixelClusterIndexToHitToken_

private

Definition at line 94 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ propagatorAlongName_

const std::string MkFitOutputConverter::propagatorAlongName_

private

Definition at line 107 of file MkFitOutputConverter.cc.

◆ propagatorAlongToken_

const edm::ESGetToken<Propagator, TrackingComponentsRecord> MkFitOutputConverter::propagatorAlongToken_

private

Definition at line 99 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ propagatorOppositeName_

const std::string MkFitOutputConverter::propagatorOppositeName_

private

Definition at line 108 of file MkFitOutputConverter.cc.

◆ propagatorOppositeToken_

const edm::ESGetToken<Propagator, TrackingComponentsRecord> MkFitOutputConverter::propagatorOppositeToken_

private

Definition at line 100 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ putSeedStopInfoToken_

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

private

Definition at line 105 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ putTrackCandidateToken_

const edm::EDPutTokenT<TrackCandidateCollection> MkFitOutputConverter::putTrackCandidateToken_

private

Definition at line 104 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ seedToken_

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

private

Definition at line 98 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ stripClusterIndexToHitToken_

const edm::EDGetTokenT<MkFitClusterIndexToHit> MkFitOutputConverter::stripClusterIndexToHitToken_

private

Definition at line 95 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ tracksToken_

const edm::EDGetTokenT<MkFitOutputWrapper> MkFitOutputConverter::tracksToken_

private

Definition at line 97 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ ttrhBuilderName_

const std::string MkFitOutputConverter::ttrhBuilderName_

private

Definition at line 106 of file MkFitOutputConverter.cc.

◆ ttrhBuilderToken_

const edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> MkFitOutputConverter::ttrhBuilderToken_

private

Definition at line 102 of file MkFitOutputConverter.cc.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MkFitOutputConverter()

◆ ~MkFitOutputConverter()

Member Function Documentation

◆ backwardFit()

◆ convertCandidates()

◆ convertInnermostState()

◆ fillDescriptions()

◆ produce()

Member Data Documentation

◆ eventOfHitsToken_

◆ mfToken_

◆ mkFitGeomToken_

◆ mkfitSeedToken_

◆ pixelClusterIndexToHitToken_

◆ propagatorAlongName_

◆ propagatorAlongToken_

◆ propagatorOppositeName_

◆ propagatorOppositeToken_

◆ putSeedStopInfoToken_

◆ putTrackCandidateToken_

◆ seedToken_

◆ stripClusterIndexToHitToken_

◆ tracksToken_

◆ ttrhBuilderName_

◆ ttrhBuilderToken_