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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

bool	wantsStreamLuminosityBlocks () const noexcept final

bool	wantsStreamRuns () const noexcept 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESResolverIndex > 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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, 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::ESRecordsToProductResolverIndices 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, const bool isPhase1, bool lastHitWasInvalid, bool lastHitWasChanged) const

std::vector< float >	computeDNNs (TrackCandidateCollection const &tkCC, const std::vector< TrajectoryStateOnSurface > &states, const reco::BeamSpot bs, const reco::VertexCollection vertices, const tensorflow::Session *session, const std::vector< float > &chi2, const bool rescaledError) 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 MkFitGeometry &mkFitGeom, const TrackerTopology &tTopo, const TkClonerImpl &hitCloner, const std::vector< const DetLayer > &detLayers, const mkfit::TrackVec &mkFitSeeds, const reco::BeamSpot bs, const reco::VertexCollection vertices, const tensorflow::Session session) 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 int	algo_

const bool	algoCandSelection_

const float	algoCandWorkingPoint_

const int	bsize_

const edm::EDGetTokenT< reco::BeamSpot >	bsToken_

const bool	doErrorRescale_

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 edm::ESGetToken< Propagator, TrackingComponentsRecord >	propagatorAlongToken_

const edm::ESGetToken< Propagator, TrackingComponentsRecord >	propagatorOppositeToken_

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

const edm::EDPutTokenT< TrackCandidateCollection >	putTrackCandidateToken_

const float	qualityMaxInvPt_

const float	qualityMaxPosErrSq_

const float	qualityMaxRsq_

const float	qualityMaxZ_

const float	qualityMinTheta_

const bool	qualitySignPt_

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

const edm::EDGetTokenT< MkFitClusterIndexToHit >	stripClusterIndexToHitToken_

const std::string	tfDnnLabel_

const edm::ESGetToken< TfGraphDefWrapper, TfGraphRecord >	tfDnnToken_

const edm::EDGetTokenT< MkFitOutputWrapper >	tracksToken_

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

const edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >	ttrhBuilderToken_

const edm::EDGetTokenT< reco::VertexCollection >	verticesToken_

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
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)

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

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

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

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

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

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

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

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

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

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

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

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

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 70 of file MkFitOutputConverter.cc.

Constructor & Destructor Documentation

◆ MkFitOutputConverter()

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

explicit

Definition at line 153 of file MkFitOutputConverter.cc.

References edm::ParameterSet::getParameter().

     : 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>()},
       tTopoToken_{esConsumes<TrackerTopology, TrackerTopologyRcd>()},
       putTrackCandidateToken_{produces<TrackCandidateCollection>()},
       putSeedStopInfoToken_{produces<std::vector<SeedStopInfo>>()},
       qualityMaxInvPt_{float(iConfig.getParameter<double>("qualityMaxInvPt"))},
       qualityMinTheta_{float(iConfig.getParameter<double>("qualityMinTheta"))},
       qualityMaxRsq_{float(pow(iConfig.getParameter<double>("qualityMaxR"), 2))},
       qualityMaxZ_{float(iConfig.getParameter<double>("qualityMaxZ"))},
       qualityMaxPosErrSq_{float(pow(iConfig.getParameter<double>("qualityMaxPosErr"), 2))},
       qualitySignPt_{iConfig.getParameter<bool>("qualitySignPt")},
       doErrorRescale_{iConfig.getParameter<bool>("doErrorRescale")},
       algo_{reco::TrackBase::algoByName(
           TString(iConfig.getParameter<edm::InputTag>("seeds").label()).ReplaceAll("Seeds", "").Data())},
       algoCandSelection_{bool(iConfig.getParameter<bool>("candMVASel"))},
       algoCandWorkingPoint_{float(iConfig.getParameter<double>("candWP"))},
       bsize_{int(iConfig.getParameter<int>("batchSize"))},
       bsToken_(algoCandSelection_ ? consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"))
                                   : edm::EDGetTokenT<reco::BeamSpot>()),
       verticesToken_(algoCandSelection_ ? consumes<reco::VertexCollection>(edm::InputTag("firstStepPrimaryVertices"))
                                         : edm::EDGetTokenT<reco::VertexCollection>()),
       tfDnnLabel_(iConfig.getParameter<std::string>("tfDnnLabel")),
       tfDnnToken_(esConsumes(edm::ESInputTag("", tfDnnLabel_))) {}

◆ ~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,
		const bool	isPhase1,
		bool	lastHitWasInvalid,
		bool	lastHitWasChanged
	)		const

private

Definition at line 534 of file MkFitOutputConverter.cc.

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

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 {
     bool doSwitch = false;
     const auto& surfacePos = lastHitSurface.position();
     const auto& lastHitPos = firstHits.front()->globalPosition();
     if (isPhase1) {
       const auto lastHitSubdet = firstHits.front()->geographicalId().subdetId();
       if (isPhase1Barrel(lastHitSubdet)) {
         doSwitch = (surfacePos.perp2() < lastHitPos.perp2());
       } else {
         doSwitch = (surfacePos.z() < lastHitPos.z());
       }
     } else {
       const auto lastHitSubdet = firstHits.front()->det()->subDetector();
       if (GeomDetEnumerators::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());
 }

◆ computeDNNs()

std::vector< float > MkFitOutputConverter::computeDNNs	(	TrackCandidateCollection const &	tkCC,
		const std::vector< TrajectoryStateOnSurface > &	states,
		const reco::BeamSpot *	bs,
		const reco::VertexCollection *	vertices,
		const tensorflow::Session *	session,
		const std::vector< float > &	chi2,
		const bool	rescaledError
	)		const

private

Definition at line 677 of file MkFitOutputConverter.cc.

References funct::abs(), algo_, cms::cuda::bs, bsize_, isoTrack_cff::chi2, f, input1, input2, PixelMapPlotter::inputs, createfilelist::int, TrajectoryStateClosestToBeamLine::isValid(), WZElectronSkims53X_cff::max, PixelVertexMonitor_cff::ndof, nt, convertSQLitetoXML_cfg::output, PatBasicFWLiteJetAnalyzer_Selector_cfg::outputs, AlCaHLTBitMon_ParallelJobs::p, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, rpcPointValidation_cfi::recHit, tensorflow::run(), run_AlCaRecoTriggerBitsUpdateWorkflow::session, nano_mu_digi_cff::strip, TrajectoryStateClosestToBeamLine::trackStateAtPCA(), trk, V0Monitor_cff::v0, bphysicsOniaDQM_cfi::vertex, and AlignmentTracksFromVertexSelector_cfi::vertices.

Referenced by convertCandidates().

                                                                                      {
   int size_in = (int)tkCC.size();
   int nbatches = size_in / bsize_;
 
   std::vector<float> output(size_in, 0);
 
   TSCBLBuilderNoMaterial tscblBuilder;
 
   tensorflow::Tensor input1(tensorflow::DT_FLOAT, {bsize_, 29});
   tensorflow::Tensor input2(tensorflow::DT_FLOAT, {bsize_, 1});
 
   for (auto nb = 0; nb < nbatches + 1; nb++) {
     std::vector<bool> invalidProp(bsize_, false);
 
     for (auto nt = 0; nt < bsize_; nt++) {
       int itrack = nt + bsize_ * nb;
       if (itrack >= size_in)
         continue;
 
       auto const& tkC = tkCC.at(itrack);
 
       TrajectoryStateOnSurface state = states.at(itrack);
 
       if (rescaledError)
         state.rescaleError(1 / 100.f);
 
       TrajectoryStateClosestToBeamLine tsAtClosestApproachTrackCand =
           tscblBuilder(*state.freeState(), *bs);  //as in TrackProducerAlgorithm
 
       if (!(tsAtClosestApproachTrackCand.isValid())) {
         edm::LogVerbatim("TrackBuilding") << "TrajectoryStateClosestToBeamLine not valid";
         invalidProp[nt] = true;
         continue;
       }
 
       auto const& stateAtPCA = tsAtClosestApproachTrackCand.trackStateAtPCA();
       auto v0 = stateAtPCA.position();
       auto p = stateAtPCA.momentum();
       math::XYZPoint pos(v0.x(), v0.y(), v0.z());
       math::XYZVector mom(p.x(), p.y(), p.z());
 
       //pseudo track for access to easy methods
       reco::Track trk(0, 0, pos, mom, stateAtPCA.charge(), stateAtPCA.curvilinearError());
 
       // get best vertex
       float dzmin = std::numeric_limits<float>::max();
       float dxy_zmin = 0;
 
       for (auto const& vertex : *vertices) {
         if (std::abs(trk.dz(vertex.position())) < dzmin) {
           dzmin = trk.dz(vertex.position());
           dxy_zmin = trk.dxy(vertex.position());
         }
       }
 
       // loop over the RecHits
       int ndof = 0;
       int pix = 0;
       int strip = 0;
       for (auto const& recHit : tkC.recHits()) {
         ndof += recHit.dimension();
         auto const subdet = recHit.geographicalId().subdetId();
         if (subdet == PixelSubdetector::PixelBarrel || subdet == PixelSubdetector::PixelEndcap)
           pix++;
         else
           strip++;
       }
       ndof = ndof - 5;
 
       input1.matrix<float>()(nt, 0) = trk.pt();  //using inner track only
       input1.matrix<float>()(nt, 1) = p.x();
       input1.matrix<float>()(nt, 2) = p.y();
       input1.matrix<float>()(nt, 3) = p.z();
       input1.matrix<float>()(nt, 4) = p.perp();
       input1.matrix<float>()(nt, 5) = p.x();
       input1.matrix<float>()(nt, 6) = p.y();
       input1.matrix<float>()(nt, 7) = p.z();
       input1.matrix<float>()(nt, 8) = p.perp();
       input1.matrix<float>()(nt, 9) = trk.ptError();
       input1.matrix<float>()(nt, 10) = dxy_zmin;
       input1.matrix<float>()(nt, 11) = dzmin;
       input1.matrix<float>()(nt, 12) = trk.dxy(bs->position());
       input1.matrix<float>()(nt, 13) = trk.dz(bs->position());
       input1.matrix<float>()(nt, 14) = trk.dxyError();
       input1.matrix<float>()(nt, 15) = trk.dzError();
       input1.matrix<float>()(nt, 16) = ndof > 0 ? chi2[itrack] / ndof : chi2[itrack] * 1e6;
       input1.matrix<float>()(nt, 17) = trk.eta();
       input1.matrix<float>()(nt, 18) = trk.phi();
       input1.matrix<float>()(nt, 19) = trk.etaError();
       input1.matrix<float>()(nt, 20) = trk.phiError();
       input1.matrix<float>()(nt, 21) = pix;
       input1.matrix<float>()(nt, 22) = strip;
       input1.matrix<float>()(nt, 23) = ndof;
       input1.matrix<float>()(nt, 24) = 0;
       input1.matrix<float>()(nt, 25) = 0;
       input1.matrix<float>()(nt, 26) = 0;
       input1.matrix<float>()(nt, 27) = 0;
       input1.matrix<float>()(nt, 28) = 0;
 
       input2.matrix<float>()(nt, 0) = algo_;
     }
 
     //inputs finalized
     tensorflow::NamedTensorList inputs;
     inputs.resize(2);
     inputs[0] = tensorflow::NamedTensor("x", input1);
     inputs[1] = tensorflow::NamedTensor("y", input2);
 
     //eval and rescale
     std::vector<tensorflow::Tensor> outputs;
     tensorflow::run(session, inputs, {"Identity"}, &outputs);
 
     for (auto nt = 0; nt < bsize_; nt++) {
       int itrack = nt + bsize_ * nb;
       if (itrack >= size_in)
         continue;
 
       float out0 = 2.0 * outputs[0].matrix<float>()(nt, 0) - 1.0;
       if (invalidProp[nt])
         out0 = -1;
 
       output[itrack] = out0;
     }
   }
 
   return output;
 }

◆ 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 MkFitGeometry &	mkFitGeom,
		const TrackerTopology &	tTopo,
		const TkClonerImpl &	hitCloner,
		const std::vector< const DetLayer *> &	detLayers,
		const mkfit::TrackVec &	mkFitSeeds,
		const reco::BeamSpot *	bs,
		const reco::VertexCollection *	vertices,
		const tensorflow::Session *	session
	)		const

private

Definition at line 267 of file MkFitOutputConverter.cc.

References a, funct::abs(), algoCandSelection_, algoCandWorkingPoint_, b, backwardFit(), cms::cuda::bs, HLT_2024v14_cff::candidates, isoTrack_cff::chi2, computeDNNs(), Options::const, convertInnermostState(), doErrorRescale_, submitPVResolutionJobs::err, InitialStepPreSplitting_cff::eventOfHits, Exception, hfClusterShapes_cfi::hits, MkFitClusterIndexToHit::hits(), mps_fire::i, GeomDetEnumerators::isBarrel(), pixelClusterTagInfos_cfi::isPhase1, MkFitGeometry::isPhase1(), fastTrackerRecHitType::isPixel(), dqmiolumiharvest::j, LogTrace, M_PI, WZElectronSkims53X_cff::max, TrackingDataMCValidation_Standalone_cff::nhits, convertSQLitetoXML_cfg::output, trajectoryStateTransform::persistentState(), MkFitOutputWrapper::propagatedToFirstLayer(), L1TMuonDQMOffline_cfi::propagatorAlong, L1TMuonDQMOffline_cfi::propagatorOpposite, qualityMaxInvPt_, qualityMaxPosErrSq_, qualityMaxRsq_, qualityMaxZ_, qualityMinTheta_, qualitySignPt_, FastTrackerRecHitMaskProducer_cfi::recHits, offlineSlimmedPrimaryVertices_cfi::score, HLT_2024v14_cff::seeds, run_AlCaRecoTriggerBitsUpdateWorkflow::session, StripSubdetector::TOB, TrackerTopology::tobSide(), MkFitOutputWrapper::tracks(), UNINITIALIZED, and AlignmentTracksFromVertexSelector_cfi::vertices.

Referenced by produce().

                                                                                                          {
   TrackCandidateCollection output;
   const auto& candidates = mkFitOutput.tracks();
   output.reserve(candidates.size());
 
   LogTrace("MkFitOutputConverter") << "Number of candidates " << candidates.size();
 
   std::vector<float> chi2;
   std::vector<TrajectoryStateOnSurface> states;
   chi2.reserve(candidates.size());
   states.reserve(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();
 
     // state: check for basic quality first
     if (cand.state().invpT() > qualityMaxInvPt_ || (qualitySignPt_ && cand.state().invpT() < 0) ||
         cand.state().theta() < qualityMinTheta_ || (M_PI - cand.state().theta()) < qualityMinTheta_ ||
         cand.state().posRsq() > qualityMaxRsq_ || std::abs(cand.state().z()) > qualityMaxZ_ ||
         (cand.state().errors.At(0, 0) + cand.state().errors.At(1, 1) + cand.state().errors.At(2, 2)) >
             qualityMaxPosErrSq_) {
       edm::LogInfo("MkFitOutputConverter")
           << "Candidate " << candIndex << " failed state quality checks" << cand.state().parameters;
       continue;
     }
 
     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::LogInfo("MkFitOutputConverter")
           << "Curvilinear error not pos-def\n"
           << fts.curvilinearError().matrix() << "\ncandidate " << candIndex << "ignored";
       continue;
     }
 
     //Sylvester's criterion, start from the smaller submatrix size
     double det = 0;
     if ((!fts.curvilinearError().matrix().Sub<AlgebraicSymMatrix22>(0, 0).Det(det)) || det < 0) {
       edm::LogInfo("MkFitOutputConverter")
           << "Fail pos-def check sub2.det for candidate " << candIndex << " with fts " << fts;
       continue;
     } else if ((!fts.curvilinearError().matrix().Sub<AlgebraicSymMatrix33>(0, 0).Det(det)) || det < 0) {
       edm::LogInfo("MkFitOutputConverter")
           << "Fail pos-def check sub3.det for candidate " << candIndex << " with fts " << fts;
       continue;
     } else if ((!fts.curvilinearError().matrix().Sub<AlgebraicSymMatrix44>(0, 0).Det(det)) || det < 0) {
       edm::LogInfo("MkFitOutputConverter")
           << "Fail pos-def check sub4.det for candidate " << candIndex << " with fts " << fts;
       continue;
     } else if ((!fts.curvilinearError().matrix().Det2(det)) || det < 0) {
       edm::LogInfo("MkFitOutputConverter")
           << "Fail pos-def check det for candidate " << candIndex << " with fts " << fts;
       continue;
     }
 
     // hits
     edm::OwnVector<TrackingRecHit> recHits;
     // nTotalHits() gives sum of valid hits (nFoundHits()) and invalid/missing hits.
     const int nhits = cand.nTotalHits();
     //std::cout << candIndex << ": " << nhits << " " << cand.nFoundHits() << std::endl;
     bool lastHitInvalid = false;
     const auto isPhase1 = mkFitGeom.isPhase1();
     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_pixel();
         auto const& hits = isPixel ? pixelClusterIndexToHit.hits() : stripClusterIndexToHit.hits();
 
         auto const& thit = static_cast<BaseTrackerRecHit const&>(*hits[hitOnTrack.index]);
         if (isPhase1) {
           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()));
           }
         } else {
           if (thit.firstClusterRef().isPixel()) {
             recHits.push_back(hits[hitOnTrack.index]->clone());
           } else if (thit.firstClusterRef().isPhase2()) {
             recHits.push_back(std::make_unique<Phase2TrackerRecHit1D>(
                 thit.localPosition(),
                 LocalError(thit.localPositionError().xx(), 0.f, std::numeric_limits<float>::max()),
                 *thit.det(),
                 thit.firstClusterRef().cluster_phase2OT()));
           }
         }
         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([&tTopo, &isPhase1](const auto& a, const auto& b) {
       //const GeomDetEnumerators::SubDetector asub = a.det()->subDetector();
       //const GeomDetEnumerators::SubDetector bsub = b.det()->subDetector();
       //const auto& apos = a.globalPosition();
       //const auto& bpos = b.globalPosition();
       // For Phase-1, can rely on subdetector index
       if (isPhase1) {
         const auto asub_ph1 = a.geographicalId().subdetId();
         const auto bsub_ph1 = b.geographicalId().subdetId();
         const auto& apos_ph1 = a.globalPosition();
         const auto& bpos_ph1 = b.globalPosition();
         if (asub_ph1 != bsub_ph1) {
           // Subdetector order (BPix, FPix, TIB, TID, TOB, TEC) corresponds also the navigation
           return asub_ph1 < bsub_ph1;
         } else {
           //if (GeomDetEnumerators::isBarrel(asub)) {
           if (isPhase1Barrel(asub_ph1)) {
             return apos_ph1.perp2() < bpos_ph1.perp2();
           } else {
             return std::abs(apos_ph1.z()) < std::abs(bpos_ph1.z());
           }
         }
       }
 
       // For Phase-2, can not rely uniquely on subdetector index
       const GeomDetEnumerators::SubDetector asub = a.det()->subDetector();
       const GeomDetEnumerators::SubDetector bsub = b.det()->subDetector();
       const auto& apos = a.globalPosition();
       const auto& bpos = b.globalPosition();
       const auto aid = a.geographicalId().rawId();
       const auto bid = b.geographicalId().rawId();
       const auto asubid = a.geographicalId().subdetId();
       const auto bsubid = b.geographicalId().subdetId();
       if (GeomDetEnumerators::isBarrel(asub) || GeomDetEnumerators::isBarrel(bsub)) {
         // For barrel tilted modules, or in case (only) one of the two modules is barrel, use 3D position
         if ((asubid == StripSubdetector::TOB && tTopo.tobSide(aid) < 3) ||
             (bsubid == StripSubdetector::TOB && tTopo.tobSide(bid) < 3) ||
             !(GeomDetEnumerators::isBarrel(asub) && GeomDetEnumerators::isBarrel(bsub))) {
           return apos.mag2() < bpos.mag2();
         }
         // For fully barrel comparisons and no tilt, use 2D position
         else {
           return apos.perp2() < bpos.perp2();
         }
       }
       // For fully endcap comparisons, use z position
       else {
         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";
 
     // Rescale candidate error if candidate is already propagated to first layer,
     // to be consistent with TransientInitialStateEstimator::innerState used in CkfTrackCandidateMakerBase
     // Error is only rescaled for candidates propagated to first layer;
     // otherwise, candidates undergo backwardFit where error is already rescaled
     if (mkFitOutput.propagatedToFirstLayer() && doErrorRescale_)
       fts.rescaleError(100.);
     auto tsosDet = mkFitOutput.propagatedToFirstLayer()
                        ? convertInnermostState(fts, recHits, propagatorAlong, propagatorOpposite)
                        : backwardFit(fts,
                                      recHits,
                                      propagatorAlong,
                                      propagatorOpposite,
                                      hitCloner,
                                      isPhase1,
                                      lastHitInvalid,
                                      lastHitChanged);
     if (!tsosDet.first.isValid()) {
       edm::LogInfo("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
     );
 
     chi2.push_back(cand.chi2());
     states.push_back(tsosDet.first);
   }
 
   if (algoCandSelection_) {
     const auto& dnnScores = computeDNNs(
         output, states, bs, vertices, session, chi2, mkFitOutput.propagatedToFirstLayer() && doErrorRescale_);
 
     TrackCandidateCollection reducedOutput;
     reducedOutput.reserve(output.size());
     int scoreIndex = 0;
     for (const auto& score : dnnScores) {
       if (score > algoCandWorkingPoint_)
         reducedOutput.push_back(output[scoreIndex]);
       scoreIndex++;
     }
 
     output.swap(reducedOutput);
   }
 
   return output;
 }

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

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

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

◆ fillDescriptions()

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

static

Definition at line 190 of file MkFitOutputConverter.cc.

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

                                                                                       {
   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"});
 
   desc.add<double>("qualityMaxInvPt", 100)->setComment("max(1/pt) for converted tracks");
   desc.add<double>("qualityMinTheta", 0.01)->setComment("lower bound on theta (or pi-theta) for converted tracks");
   desc.add<double>("qualityMaxR", 120)->setComment("max(R) for the state position for converted tracks");
   desc.add<double>("qualityMaxZ", 280)->setComment("max(|Z|) for the state position for converted tracks");
   desc.add<double>("qualityMaxPosErr", 100)->setComment("max position error for converted tracks");
   desc.add<bool>("qualitySignPt", true)->setComment("check sign of 1/pt for converted tracks");
 
   desc.add<bool>("doErrorRescale", true)->setComment("rescale candidate error before final fit");
 
   desc.add<std::string>("tfDnnLabel", "trackSelectionTf");
 
   desc.add<bool>("candMVASel", false)->setComment("flag used to trigger MVA selection at cand level");
   desc.add<double>("candWP", 0)->setComment("MVA selection at cand level working point");
   desc.add<int>("batchSize", 16)->setComment("batch size for cand DNN evaluation");
 
   descriptions.addWithDefaultLabel(desc);
 }

◆ produce()

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

overrideprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 221 of file MkFitOutputConverter.cc.

References algoCandSelection_, ConfigBuilder::beamspot, bsToken_, Options::const, convertCandidates(), eventOfHitsToken_, Exception, edm::EventSetup::getData(), iEvent, mfToken_, mkFitGeomToken_, mkfitSeedToken_, pixelClusterIndexToHitToken_, propagatorAlongToken_, propagatorOppositeToken_, putSeedStopInfoToken_, putTrackCandidateToken_, HLT_2024v14_cff::seeds, seedToken_, run_AlCaRecoTriggerBitsUpdateWorkflow::session, stripClusterIndexToHitToken_, tfDnnToken_, tracksToken_, tTopoToken_, ttrhBuilderToken_, AlignmentTracksFromVertexSelector_cfi::vertices, and verticesToken_.

                                                                                                        {
   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_);
   // primary vertices under the algo_ because in initialStepPreSplitting there are no firstStepPrimaryVertices
   // beamspot as well since the producer can be used in hlt
   const reco::VertexCollection* vertices = nullptr;
   const reco::BeamSpot* beamspot = nullptr;
   if (algoCandSelection_) {
     vertices = &iEvent.get(verticesToken_);
     beamspot = &iEvent.get(bsToken_);
   }
 
   const tensorflow::Session* session = nullptr;
   if (algoCandSelection_)
     session = iSetup.getData(tfDnnToken_).getSession();
 
   // 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_),
                                    iSetup.getData(mkFitGeomToken_),
                                    iSetup.getData(tTopoToken_),
                                    tkBuilder->cloner(),
                                    mkFitGeom.detLayers(),
                                    mkfitSeeds.seeds(),
                                    beamspot,
                                    vertices,
                                    session));
 
   // TODO: SeedStopInfo is currently unfilled
   iEvent.emplace(putSeedStopInfoToken_, seeds.size());
 }

Member Data Documentation

◆ algo_

const int MkFitOutputConverter::algo_

private

Definition at line 143 of file MkFitOutputConverter.cc.

Referenced by computeDNNs().

◆ algoCandSelection_

const bool MkFitOutputConverter::algoCandSelection_

private

Definition at line 144 of file MkFitOutputConverter.cc.

Referenced by convertCandidates(), and produce().

◆ algoCandWorkingPoint_

const float MkFitOutputConverter::algoCandWorkingPoint_

private

Definition at line 145 of file MkFitOutputConverter.cc.

Referenced by convertCandidates().

◆ bsize_

const int MkFitOutputConverter::bsize_

private

Definition at line 146 of file MkFitOutputConverter.cc.

Referenced by computeDNNs().

◆ bsToken_

const edm::EDGetTokenT<reco::BeamSpot> MkFitOutputConverter::bsToken_

private

Definition at line 147 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ doErrorRescale_

const bool MkFitOutputConverter::doErrorRescale_

private

Definition at line 141 of file MkFitOutputConverter.cc.

Referenced by convertCandidates().

◆ eventOfHitsToken_

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

private

Definition at line 119 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ mfToken_

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

private

Definition at line 127 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ mkFitGeomToken_

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

private

Definition at line 129 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ mkfitSeedToken_

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

private

Definition at line 122 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ pixelClusterIndexToHitToken_

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

private

Definition at line 120 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ propagatorAlongToken_

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

private

Definition at line 125 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ propagatorOppositeToken_

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

private

Definition at line 126 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ putSeedStopInfoToken_

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

private

Definition at line 132 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ putTrackCandidateToken_

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

private

Definition at line 131 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ qualityMaxInvPt_

const float MkFitOutputConverter::qualityMaxInvPt_

private

Definition at line 134 of file MkFitOutputConverter.cc.

Referenced by convertCandidates().

◆ qualityMaxPosErrSq_

const float MkFitOutputConverter::qualityMaxPosErrSq_

private

Definition at line 138 of file MkFitOutputConverter.cc.

Referenced by convertCandidates().

◆ qualityMaxRsq_

const float MkFitOutputConverter::qualityMaxRsq_

private

Definition at line 136 of file MkFitOutputConverter.cc.

Referenced by convertCandidates().

◆ qualityMaxZ_

const float MkFitOutputConverter::qualityMaxZ_

private

Definition at line 137 of file MkFitOutputConverter.cc.

Referenced by convertCandidates().

◆ qualityMinTheta_

const float MkFitOutputConverter::qualityMinTheta_

private

Definition at line 135 of file MkFitOutputConverter.cc.

Referenced by convertCandidates().

◆ qualitySignPt_

const bool MkFitOutputConverter::qualitySignPt_

private

Definition at line 139 of file MkFitOutputConverter.cc.

Referenced by convertCandidates().

◆ seedToken_

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

private

Definition at line 124 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ stripClusterIndexToHitToken_

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

private

Definition at line 121 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ tfDnnLabel_

const std::string MkFitOutputConverter::tfDnnLabel_

private

Definition at line 149 of file MkFitOutputConverter.cc.

◆ tfDnnToken_

const edm::ESGetToken<TfGraphDefWrapper, TfGraphRecord> MkFitOutputConverter::tfDnnToken_

private

Definition at line 150 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ tracksToken_

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

private

Definition at line 123 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ tTopoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> MkFitOutputConverter::tTopoToken_

private

Definition at line 130 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ ttrhBuilderToken_

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

private

Definition at line 128 of file MkFitOutputConverter.cc.

Referenced by produce().

◆ verticesToken_

const edm::EDGetTokenT<reco::VertexCollection> MkFitOutputConverter::verticesToken_

private

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

◆ computeDNNs()

◆ convertCandidates()

◆ convertInnermostState()

◆ fillDescriptions()

◆ produce()

Member Data Documentation

◆ algo_

◆ algoCandSelection_

◆ algoCandWorkingPoint_

◆ bsize_

◆ bsToken_

◆ doErrorRescale_

◆ eventOfHitsToken_

◆ mfToken_

◆ mkFitGeomToken_

◆ mkfitSeedToken_

◆ pixelClusterIndexToHitToken_

◆ propagatorAlongToken_

◆ propagatorOppositeToken_

◆ putSeedStopInfoToken_

◆ putTrackCandidateToken_

◆ qualityMaxInvPt_

◆ qualityMaxPosErrSq_

◆ qualityMaxRsq_

◆ qualityMaxZ_

◆ qualityMinTheta_

◆ qualitySignPt_

◆ seedToken_

◆ stripClusterIndexToHitToken_

◆ tfDnnLabel_

◆ tfDnnToken_

◆ tracksToken_

◆ tTopoToken_

◆ ttrhBuilderToken_

◆ verticesToken_