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

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

Constructor & Destructor Documentation

◆ MkFitOutputConverter()

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

explicit

Definition at line 117 of file MkFitOutputConverter.cc.

     : 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")} {}

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 381 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 Chi2MeasurementEstimator_cfi::Chi2MeasurementEstimator, 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 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 223 of file MkFitOutputConverter.cc.

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

References a, funct::abs(), b, backwardFit(), backwardFitInCMSSW_, HLT_FULL_cff::candidates, MkFitOutputWrapper::candidateTracks(), TrackingRecHit::clone(), MkFitHitIndexMap::clusterIndex(), convertInnermostState(), submitPVResolutionJobs::err, Exception, MkFitOutputWrapper::fitTracks(), MkFitHitIndexMap::hitPtr(), mps_fire::i, MkFitHitIndexMap::MkFitHit::index(), PixelPluginsPhase0_cfi::isBarrel, dqmiolumiharvest::j, LogTrace, nhits, convertSQLitetoXML_cfg::output, trajectoryStateTransform::persistentState(), HLT_FULL_cff::propagatorAlong, HLT_FULL_cff::propagatorOpposite, FastTrackerRecHitMaskProducer_cfi::recHits, InitialStep_cff::seeds, 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 493 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().

◆ createDetLayers()

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.

                                                                                                        {
   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 det rawId " << detId.rawId()
                                          << " "
                                          << " 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;
 }

References AlCaHLTBitMon_QueryRunRegistry::comp, Exception, phase1PixelTopology::layer, TrackerTopology::layer(), or, PosEndcap, TrackerTopology::side(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, PbPb_ZMuSkimMuonDPG_cff::tracker, and UNLIKELY.

Referenced by produce().

◆ fillDescriptions()

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

static

Definition at line 135 of file MkFitOutputConverter.cc.

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

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

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

References convertCandidates(), createDetLayers(), Exception, geomToken_, edm::EventSetup::getData(), InitialStep_cff::hitsSeeds, hitsSeedsToken_, iEvent, mfToken_, mteToken_, propagatorAlongToken_, propagatorOppositeToken_, putSeedStopInfoToken_, putTrackCandidateToken_, InitialStep_cff::seeds, seedToken_, tracksToken_, ttopoToken_, and ttrhBuilderToken_.

Member Data Documentation

◆ backwardFitInCMSSW_

bool MkFitOutputConverter::backwardFitInCMSSW_

private

Definition at line 114 of file MkFitOutputConverter.cc.

Referenced by convertCandidates().

◆ geomToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> MkFitOutputConverter::geomToken_

private

Definition at line 103 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ hitsSeedsToken_

edm::EDGetTokenT<MkFitInputWrapper> MkFitOutputConverter::hitsSeedsToken_

private

Definition at line 99 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ mfToken_

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

private

Definition at line 107 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ mteToken_

edm::EDGetTokenT<MeasurementTrackerEvent> MkFitOutputConverter::mteToken_

private

Definition at line 102 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ propagatorAlongName_

std::string MkFitOutputConverter::propagatorAlongName_

private

Definition at line 112 of file MkFitOutputConverter.cc.

◆ propagatorAlongToken_

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

private

Definition at line 104 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ propagatorOppositeName_

std::string MkFitOutputConverter::propagatorOppositeName_

private

Definition at line 113 of file MkFitOutputConverter.cc.

◆ propagatorOppositeToken_

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

private

Definition at line 105 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ putSeedStopInfoToken_

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

private

Definition at line 110 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ putTrackCandidateToken_

edm::EDPutTokenT<TrackCandidateCollection> MkFitOutputConverter::putTrackCandidateToken_

private

Definition at line 109 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ seedToken_

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

private

Definition at line 101 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ tracksToken_

edm::EDGetTokenT<MkFitOutputWrapper> MkFitOutputConverter::tracksToken_

private

Definition at line 100 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ ttopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> MkFitOutputConverter::ttopoToken_

private

Definition at line 106 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ ttrhBuilderName_

std::string MkFitOutputConverter::ttrhBuilderName_

private

Definition at line 111 of file MkFitOutputConverter.cc.

◆ ttrhBuilderToken_

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

private

Definition at line 108 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()

◆ createDetLayers()

◆ fillDescriptions()

◆ produce()

Member Data Documentation

◆ backwardFitInCMSSW_

◆ geomToken_

◆ hitsSeedsToken_

◆ mfToken_

◆ mteToken_

◆ propagatorAlongName_

◆ propagatorAlongToken_

◆ propagatorOppositeName_

◆ propagatorOppositeToken_

◆ putSeedStopInfoToken_

◆ putTrackCandidateToken_

◆ seedToken_

◆ tracksToken_

◆ ttopoToken_

◆ ttrhBuilderName_

◆ ttrhBuilderToken_