#include <MuonSeedTrack.h>

Inheritance diagram for MuonSeedTrack:

Public Member Functions
	MuonSeedTrack (const edm::ParameterSet &)
	constructor with config More...

	~MuonSeedTrack () override
	destructor More...

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

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

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

Private Member Functions
void	beginJob () override
	pre-job booking More...

std::pair< bool, reco::Track >	buildTrackAtPCA (const TrajectorySeed &) const
	Build a track at the PCA WITHOUT any vertex constriant. More...

void	endJob () override
	post-job More...

TrajectoryStateOnSurface	getSeedTSOS (const TrajectorySeed &seed) const
	get the TrajectorySeed's TrajectoryStateOnSurface More...

void	produce (edm::Event &, const edm::EventSetup &) override
	construct proto-tracks More...

void	setAlias (std::string alias)
	set the Branch alias More...

Private Attributes
std::string	theAlias
	the Branch alias More...

bool	theAllowNoVtxFlag

edm::InputTag	theSeedsLabel
	the propagator used for extracting TSOS from seed More...

edm::EDGetTokenT< TrajectorySeedCollection >	theSeedsToken

MuonServiceProxy *	theService
	the event setup proxy, it takes care the services update More...

MuonUpdatorAtVertex *	theUpdatorAtVtx
	the class used for updating a Trajectory State at veretex More...

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Make a (fake) reco::Track from a TrajectorySeed. The (fake) track can be associated to a TrackingParticle (SimTrack) using SimTracker/TrackAssociation. The association can then be used in Validation packages, such as Validation/RecoMuon/MuonTrackValidator.cc

Author: Adam Everett Purdue University

Definition at line 46 of file MuonSeedTrack.h.

Constructor & Destructor Documentation

◆ MuonSeedTrack()

MuonSeedTrack::MuonSeedTrack ( const edm::ParameterSet & pset )

explicit

constructor with config

Definition at line 39 of file MuonSeedTrack.cc.

References MuonServiceProxy_cff::MuonServiceProxy, MuonUpdatorAtVertex_cff::MuonUpdatorAtVertex, muonDTDigis_cfi::pset, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                         {
   // service parameters
   ParameterSet serviceParameters = pset.getParameter<ParameterSet>("ServiceParameters");
   theService = new MuonServiceProxy(serviceParameters, consumesCollector());
 
   ParameterSet updatorPar = pset.getParameter<ParameterSet>("MuonUpdatorAtVertexParameters");
   //theSeedPropagatorName = updatorPar.getParameter<string>("Propagator");
   theSeedsLabel = pset.getParameter<InputTag>("MuonSeed");
   theSeedsToken = consumes<TrajectorySeedCollection>(theSeedsLabel);
   theUpdatorAtVtx = new MuonUpdatorAtVertex(updatorPar, theService);
 
   theAllowNoVtxFlag = pset.getUntrackedParameter<bool>("AllowNoVertex", false);
 
   //register products
   setAlias(pset.getParameter<std::string>("@module_label"));
   produces<reco::TrackCollection>().setBranchAlias(theAlias + "Tracks");
 }

◆ ~MuonSeedTrack()

MuonSeedTrack::~MuonSeedTrack ( )

override

destructor

Definition at line 60 of file MuonSeedTrack.cc.

                               {
   if (theService)
     delete theService;
   if (theUpdatorAtVtx)
     delete theUpdatorAtVtx;
 }

Member Function Documentation

◆ beginJob()

void MuonSeedTrack::beginJob ( void )

overrideprivatevirtual

pre-job booking

empty method

Reimplemented from edm::one::EDProducerBase.

Definition at line 104 of file MuonSeedTrack.cc.

104 {}

◆ buildTrackAtPCA()

pair< bool, reco::Track > MuonSeedTrack::buildTrackAtPCA ( const TrajectorySeed & seed ) const

private

Build a track at the PCA WITHOUT any vertex constriant.

First calls getSeedTSOS, then propagates to the vertex. Creates a reco::Track from the propagated initial FreeTrajectoryState.

Definition at line 155 of file MuonSeedTrack.cc.

References FreeTrajectoryState::charge(), FreeTrajectoryState::curvilinearError(), debug, TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::globalPosition(), TrackerBounds::isInside(), LogTrace, metname, FreeTrajectoryState::momentum(), AlCaHLTBitMon_ParallelJobs::p, FreeTrajectoryState::position(), fileCollector::seed, HLT_2023v12_cff::track, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                      {
   const string metname = "MuonSeedTrack";
 
   MuonPatternRecoDumper debug;
 
   TSOS seedTSOS = getSeedTSOS(seed);
   // This is needed to extrapolate the tsos at vertex
   LogTrace(metname) << "Propagate to PCA...";
   pair<bool, FreeTrajectoryState> extrapolationResult = theUpdatorAtVtx->propagateToNominalLine(seedTSOS);
   FreeTrajectoryState ftsAtVtx;
 
   if (extrapolationResult.first) {
     ftsAtVtx = extrapolationResult.second;
   } else {
     if (TrackerBounds::isInside(seedTSOS.globalPosition())) {
       LogWarning(metname) << "Track in the Tracker: taking the innermost state instead of the state at PCA";
       ftsAtVtx = *seedTSOS.freeState();
     } else {
       if (theAllowNoVtxFlag) {
         LogWarning(metname) << "Propagation to PCA failed, taking the innermost state instead of the state at PCA";
         ftsAtVtx = *seedTSOS.freeState();
       } else {
         LogWarning(metname) << "Stand Alone track: this track will be rejected";
         return pair<bool, reco::Track>(false, reco::Track());
       }
     }
   }
 
   LogTrace(metname) << "TSOS after the extrapolation at vtx";
   LogTrace(metname) << debug.dumpFTS(ftsAtVtx);
 
   GlobalPoint pca = ftsAtVtx.position();
   math::XYZPoint persistentPCA(pca.x(), pca.y(), pca.z());
   GlobalVector p = ftsAtVtx.momentum();
   math::XYZVector persistentMomentum(p.x(), p.y(), p.z());
 
   double dummyNDOF = 1.0;
   double dummyChi2 = 1.0;
 
   reco::Track track(
       dummyChi2, dummyNDOF, persistentPCA, persistentMomentum, ftsAtVtx.charge(), ftsAtVtx.curvilinearError());
 
   return pair<bool, reco::Track>(true, track);
 }

◆ endJob()

void MuonSeedTrack::endJob ( void )

overrideprivatevirtual

post-job

empty method

Reimplemented from edm::one::EDProducerBase.

Definition at line 109 of file MuonSeedTrack.cc.

109 {}

◆ getSeedTSOS()

TrajectoryStateOnSurface MuonSeedTrack::getSeedTSOS ( const TrajectorySeed & seed ) const

private

get the TrajectorySeed's TrajectoryStateOnSurface

Get the TrajectoryStateOnSurface from the TrajectorySeed

Definition at line 114 of file MuonSeedTrack.cc.

References PTrajectoryStateOnDet::detId(), fileCollector::seed, and trajectoryStateTransform::transientState().

                                                                                     {
   // Get the Trajectory State on Det (persistent version of a TSOS) from the seed
   PTrajectoryStateOnDet pTSOD = seed.startingState();
 
   // Transform it in a TrajectoryStateOnSurface
 
   DetId seedDetId(pTSOD.detId());
 
   const GeomDet* gdet = theService->trackingGeometry()->idToDet(seedDetId);
 
   TrajectoryStateOnSurface initialState =
       trajectoryStateTransform::transientState(pTSOD, &(gdet->surface()), &*theService->magneticField());
 
   /*
   // Get the layer on which the seed relies
   const DetLayer *initialLayer = theService->detLayerGeometry()->idToLayer( seedDetId );
 
   PropagationDirection detLayerOrder = oppositeToMomentum;
 
   // ask for compatible layers
   vector<const DetLayer*> detLayers;
   detLayers = initialLayer->compatibleLayers( *initialState.freeState(),detLayerOrder);
   
   TrajectoryStateOnSurface result = initialState;
   if(detLayers.size()){
     const DetLayer* finalLayer = detLayers.back();
     const TrajectoryStateOnSurface propagatedState = theService->propagator(theSeedPropagatorName)->propagate(initialState, finalLayer->surface());
     if(propagatedState.isValid())
       result = propagatedState;
   }
   
   return result;
   */
 
   return initialState;
 }

◆ produce()

void MuonSeedTrack::produce	(	edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overrideprivatevirtual

construct proto-tracks

For each seed, make a (fake) reco::Track

Implements edm::one::EDProducerBase.

Definition at line 74 of file MuonSeedTrack.cc.

References options_cfi::eventSetup, eostools::move(), HLT_FULL_cff::seeds, HLT_2023v12_cff::track, and JetHT_cfg::trackCollection.

                                                                             {
   using namespace edm;
 
   // Update the services
   theService->update(eventSetup);
 
   // the track collectios; they will be loaded in the event
   unique_ptr<reco::TrackCollection> trackCollection(new reco::TrackCollection());
   // ... and its reference into the event
   reco::TrackRefProd trackCollectionRefProd = event.getRefBeforePut<reco::TrackCollection>();
 
   Handle<TrajectorySeedCollection> seeds;
   event.getByToken(theSeedsToken, seeds);
 
   for (TrajectorySeedCollection::const_iterator iSeed = seeds->begin(); iSeed != seeds->end(); iSeed++) {
     pair<bool, reco::Track> resultOfTrackExtrapAtPCA = buildTrackAtPCA(*iSeed);
     if (!resultOfTrackExtrapAtPCA.first)
       continue;
     // take the "bare" track at PCA
     reco::Track& track = resultOfTrackExtrapAtPCA.second;
     // fill the TrackCollection
     trackCollection->push_back(track);
   }
 
   event.put(std::move(trackCollection));
 }

◆ setAlias()

void MuonSeedTrack::setAlias ( std::string alias )

inlineprivate

set the Branch alias

Definition at line 68 of file MuonSeedTrack.h.

References SiStripOfflineCRack_cfg::alias, and theAlias.

                                  {
     alias.erase(alias.size() - 1, alias.size());
     theAlias = alias;
   }

Member Data Documentation

◆ theAlias

std::string MuonSeedTrack::theAlias

private

the Branch alias

Definition at line 85 of file MuonSeedTrack.h.

Referenced by setAlias().

◆ theAllowNoVtxFlag

bool MuonSeedTrack::theAllowNoVtxFlag

private

Definition at line 95 of file MuonSeedTrack.h.

◆ theSeedsLabel

edm::InputTag MuonSeedTrack::theSeedsLabel

private

the propagator used for extracting TSOS from seed

the TrajectorySeed label

Definition at line 91 of file MuonSeedTrack.h.

◆ theSeedsToken

edm::EDGetTokenT<TrajectorySeedCollection> MuonSeedTrack::theSeedsToken

private

Definition at line 92 of file MuonSeedTrack.h.

◆ theService

MuonServiceProxy* MuonSeedTrack::theService

private

the event setup proxy, it takes care the services update

Definition at line 79 of file MuonSeedTrack.h.

◆ theUpdatorAtVtx

MuonUpdatorAtVertex* MuonSeedTrack::theUpdatorAtVtx

private

the class used for updating a Trajectory State at veretex

Definition at line 82 of file MuonSeedTrack.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MuonSeedTrack()

◆ ~MuonSeedTrack()

Member Function Documentation

◆ beginJob()

◆ buildTrackAtPCA()

◆ endJob()

◆ getSeedTSOS()

◆ produce()

◆ setAlias()

Member Data Documentation

◆ theAlias

◆ theAllowNoVtxFlag

◆ theSeedsLabel

◆ theSeedsToken

◆ theService

◆ theUpdatorAtVtx