#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

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)

	~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

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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)

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 45 of file MuonSeedTrack.h.

Constructor & Destructor Documentation

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

explicit

constructor with config

Definition at line 39 of file MuonSeedTrack.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), MuonServiceProxy_cff::MuonServiceProxy, MuonUpdatorAtVertex_cff::MuonUpdatorAtVertex, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   // service parameters
   ParameterSet serviceParameters = pset.getParameter<ParameterSet>("ServiceParameters");
   theService = new MuonServiceProxy(serviceParameters);
   
   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 ( )

override

destructor

Definition at line 62 of file MuonSeedTrack.cc.

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

Member Function Documentation

void MuonSeedTrack::beginJob ( void )

overrideprivatevirtual

pre-job booking

empty method

Reimplemented from edm::EDProducer.

Definition at line 110 of file MuonSeedTrack.cc.

111 {

112 }

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 170 of file MuonSeedTrack.cc.

References FreeTrajectoryState::charge(), FreeTrajectoryState::curvilinearError(), debug, MuonPatternRecoDumper::dumpFTS(), TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::globalPosition(), TrackerBounds::isInside(), LogTrace, metname, FreeTrajectoryState::momentum(), AlCaHLTBitMon_ParallelJobs::p, FreeTrajectoryState::position(), HiIsolationCommonParameters_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 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);
 }

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 226 of file MuonSeedTrack.cc.

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

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

void MuonSeedTrack::endJob ( void )

overrideprivatevirtual

post-job

empty method

Reimplemented from edm::EDProducer.

Definition at line 119 of file MuonSeedTrack.cc.

119 {

120 }

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

private

get the TrajectorySeed's TrajectoryStateOnSurface

Get the TrajectoryStateOnSurface from the TrajectorySeed

Definition at line 126 of file MuonSeedTrack.cc.

References PTrajectoryStateOnDet::detId(), TrajectorySeed::startingState(), 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;
 
 }

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

overrideprivate

construct proto-tracks

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

Definition at line 76 of file MuonSeedTrack.cc.

References eostools::move(), HiIsolationCommonParameters_cff::track, and findElectronsInSiStrips_cfi::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));
 
 }

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

inlineprivate

set the Branch alias

Definition at line 67 of file MuonSeedTrack.h.

References egammaCTFFinalFitWithMaterial_cff::alias.

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

Member Data Documentation

std::string MuonSeedTrack::theAlias

private

the Branch alias

Definition at line 87 of file MuonSeedTrack.h.

bool MuonSeedTrack::theAllowNoVtxFlag

private

Definition at line 97 of file MuonSeedTrack.h.

edm::InputTag MuonSeedTrack::theSeedsLabel

private

the propagator used for extracting TSOS from seed

the TrajectorySeed label

Definition at line 93 of file MuonSeedTrack.h.

edm::EDGetTokenT<TrajectorySeedCollection> MuonSeedTrack::theSeedsToken

private

Definition at line 94 of file MuonSeedTrack.h.

MuonServiceProxy* MuonSeedTrack::theService

private

the event setup proxy, it takes care the services update

Definition at line 81 of file MuonSeedTrack.h.

MuonUpdatorAtVertex* MuonSeedTrack::theUpdatorAtVtx

private

the class used for updating a Trajectory State at veretex

Definition at line 84 of file MuonSeedTrack.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation