MuonSeedTrack Class Reference

#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::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () 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
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	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)

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...
double	computeNDOF (const TrajectorySeed &) const
	compute the TrajectorySeed's degree of freedom 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::EDProducer
typedef EDProducer	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
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
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>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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

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.

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

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

~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::EDProducer.

Definition at line 104 of file MuonSeedTrack.cc.

{}

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.

                                                                                     {
  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 ndof = computeNDOF(seed);
  double dummyChi2 = 1.0;
 
  reco::Track track(
      dummyChi2, dummyNDOF, persistentPCA, persistentMomentum, ftsAtVtx.charge(), ftsAtVtx.curvilinearError());
 
  return pair<bool, reco::Track>(true, track);
}

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

computeNDOF()

double MuonSeedTrack::computeNDOF ( const TrajectorySeed &  trajectory ) const

private

compute the TrajectorySeed's degree of freedom

Calculates number of degrees of freedom for the TrajectorySeed

Definition at line 204 of file MuonSeedTrack.cc.

                                                                        {
  const string metname = "MuonSeedTrack";
 
  BasicTrajectorySeed::const_iterator recHits1 = (trajectory.recHits().first);
  BasicTrajectorySeed::const_iterator recHits2 = (trajectory.recHits().second);
 
  double ndof = 0.;
 
  if ((*recHits1).isValid())
    ndof += (*recHits1).dimension();
  if ((*recHits2).isValid())
    ndof += (*recHits2).dimension();
 
  //const Trajectory::RecHitContainer transRecHits = trajectory.recHits();
  //for(Trajectory::RecHitContainer::const_iterator rechit = transRecHits.begin();
  //  rechit != transRecHits.end(); ++rechit)
  //if ((*rechit)->isValid()) ndof += (*rechit)->dimension();
 
  // FIXME! in case of Boff is dof - 4
  return max(ndof - 5., 0.);
}

References SiStripPI::max, metname, ndof, and TrajectorySeed::recHits().

endJob()

void MuonSeedTrack::endJob ( void  )

overrideprivatevirtual

post-job

empty method

Reimplemented from edm::EDProducer.

Definition at line 109 of file MuonSeedTrack.cc.

{}

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.

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

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

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::EDProducer.

Definition at line 74 of file MuonSeedTrack.cc.

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

References eostools::move(), InitialStep_cff::seeds, HLT_2018_cff::track, and duplicaterechits_cfi::trackCollection.

setAlias()

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

inlineprivate

set the Branch alias

Definition at line 68 of file MuonSeedTrack.h.

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

References SiStripOfflineCRack_cfg::alias, and theAlias.

Member Data Documentation

theAlias

std::string MuonSeedTrack::theAlias

private

the Branch alias

Definition at line 88 of file MuonSeedTrack.h.

Referenced by setAlias().

theAllowNoVtxFlag

bool MuonSeedTrack::theAllowNoVtxFlag

private

Definition at line 98 of file MuonSeedTrack.h.

theSeedsLabel

edm::InputTag MuonSeedTrack::theSeedsLabel

private

the propagator used for extracting TSOS from seed

the TrajectorySeed label

Definition at line 94 of file MuonSeedTrack.h.

theSeedsToken

edm::EDGetTokenT<TrajectorySeedCollection> MuonSeedTrack::theSeedsToken

private

Definition at line 95 of file MuonSeedTrack.h.

theService

MuonServiceProxy* MuonSeedTrack::theService

private

the event setup proxy, it takes care the services update

Definition at line 82 of file MuonSeedTrack.h.

theUpdatorAtVtx

MuonUpdatorAtVertex* MuonSeedTrack::theUpdatorAtVtx

private

the class used for updating a Trajectory State at veretex

Definition at line 85 of file MuonSeedTrack.h.