#include <TrackAlgoCompareUtil.h>

Inheritance diagram for TrackAlgoCompareUtil:

Public Member Functions
	TrackAlgoCompareUtil (const edm::ParameterSet &)

	~TrackAlgoCompareUtil () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

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

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override

void	SetTrackingParticleD0Dz (TrackingParticleRef tp, const reco::BeamSpot &bs, const MagneticField *bf, TPtoRecoTrack &TPRT) const

void	SetTrackingParticleD0Dz (TrackingParticleRef tp, const reco::BeamSpot &bs, const MagneticField *bf, RecoTracktoTP &RTTP) const

Private Attributes
edm::EDGetTokenT< reco::RecoToSimCollection >	associatormap_algoA_recoToSim

edm::EDGetTokenT< reco::SimToRecoCollection >	associatormap_algoA_simToReco

edm::EDGetTokenT< reco::RecoToSimCollection >	associatormap_algoB_recoToSim

edm::EDGetTokenT< reco::SimToRecoCollection >	associatormap_algoB_simToReco

edm::EDGetTokenT< reco::TrackToTrackingParticleAssociator >	assocLabel_algoA

edm::EDGetTokenT< reco::TrackToTrackingParticleAssociator >	assocLabel_algoB

edm::EDGetTokenT< reco::BeamSpot >	beamSpotLabel

edm::EDGetTokenT< TrackingParticleCollection >	trackingParticleLabel_effic

edm::EDGetTokenT< TrackingParticleCollection >	trackingParticleLabel_fakes

edm::EDGetTokenT< edm::View< reco::Track > >	trackLabel_algoA

edm::EDGetTokenT< edm::View< reco::Track > >	trackLabel_algoB

const bool	UseAssociators

const bool	UseVertex

edm::EDGetTokenT< reco::VertexCollection >	vertexLabel_algoA

edm::EDGetTokenT< reco::VertexCollection >	vertexLabel_algoB

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

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Definition at line 53 of file TrackAlgoCompareUtil.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 8 of file TrackAlgoCompareUtil.cc.

References associatormap_algoA_recoToSim, associatormap_algoA_simToReco, associatormap_algoB_recoToSim, associatormap_algoB_simToReco, assocLabel_algoA, assocLabel_algoB, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), AlCaHLTBitMon_QueryRunRegistry::string, UseAssociators, UseVertex, vertexLabel_algoA, and vertexLabel_algoB.

                                                                         :
   trackLabel_algoA(consumes<View<reco::Track>>(iConfig.getParameter<edm::InputTag>("trackLabel_algoA"))),
   trackLabel_algoB(consumes<View<reco::Track>>(iConfig.getParameter<edm::InputTag>("trackLabel_algoB"))),
   trackingParticleLabel_fakes(consumes<TrackingParticleCollection>(iConfig.getParameter<edm::InputTag>("trackingParticleLabel_fakes"))),
   trackingParticleLabel_effic(consumes<TrackingParticleCollection>(iConfig.getParameter<edm::InputTag>("trackingParticleLabel_effic"))),
   beamSpotLabel(consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpotLabel"))),
   UseAssociators(iConfig.getParameter< bool >("UseAssociators")),
   UseVertex(iConfig.getParameter< bool >("UseVertex"))
 {
     //now do what ever other initialization is needed
   if(UseVertex) {
     vertexLabel_algoA = consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("vertexLabel_algoA"));
     vertexLabel_algoB = consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("vertexLabel_algoB"));
   }
 
   if(UseAssociators) {
     assocLabel_algoA = consumes<reco::TrackToTrackingParticleAssociator>(iConfig.getUntrackedParameter<std::string>("assocLabel_algoA", "trackAssociatorByHits"));
     assocLabel_algoB = consumes<reco::TrackToTrackingParticleAssociator>(iConfig.getUntrackedParameter<std::string>("assocLabel_algoB", "trackAssociatorByHits"));
   }
   else {
     edm::InputTag algoA = iConfig.getParameter< edm::InputTag >("associatormap_algoA");
     edm::InputTag algoB = iConfig.getParameter< edm::InputTag >("associatormap_algoB");
 
     associatormap_algoA_recoToSim = consumes<reco::RecoToSimCollection>(algoA);
     associatormap_algoB_recoToSim = consumes<reco::RecoToSimCollection>(algoB);
     associatormap_algoA_simToReco = consumes<reco::SimToRecoCollection>(algoA);
     associatormap_algoB_simToReco = consumes<reco::SimToRecoCollection>(algoB);
   }
   
     produces<RecoTracktoTPCollection>("AlgoA");
     produces<RecoTracktoTPCollection>("AlgoB");
     produces<TPtoRecoTrackCollection>("TP");
 }

TrackAlgoCompareUtil::~TrackAlgoCompareUtil ( )

override

Definition at line 43 of file TrackAlgoCompareUtil.cc.

44 {

45 }

Member Function Documentation

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

overrideprivate

Definition at line 50 of file TrackAlgoCompareUtil.cc.

References reco::TrackToTrackingParticleAssociator::associateRecoToSim(), reco::TrackToTrackingParticleAssociator::associateSimToReco(), associatormap_algoA_recoToSim, associatormap_algoA_simToReco, associatormap_algoB_recoToSim, associatormap_algoB_simToReco, assocLabel_algoA, assocLabel_algoB, ecalDrivenElectronSeedsParameters_cff::beamSpot, beamSpotLabel, edm::AssociationMap< Tag >::end(), edm::AssociationMap< Tag >::find(), edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, seedCreatorFromRegionConsecutiveHitsEDProducer_cff::magneticField, eostools::move(), reco::BeamSpot::position(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), edm::Event::put(), RecoTracktoTP::SetBeamSpot(), TPtoRecoTrack::SetBeamSpot(), RecoTracktoTP::SetRecoTrack(), TPtoRecoTrack::SetRecoTrack_AlgoA(), TPtoRecoTrack::SetRecoTrack_AlgoB(), RecoTracktoTP::SetRecoVertex(), TPtoRecoTrack::SetRecoVertex_AlgoA(), TPtoRecoTrack::SetRecoVertex_AlgoB(), RecoTracktoTP::SetShared(), TPtoRecoTrack::SetShared_AlgoA(), TPtoRecoTrack::SetShared_AlgoB(), TPtoRecoTrack::SetTrackingParticle(), RecoTracktoTP::SetTrackingParticle(), SetTrackingParticleD0Dz(), trackingParticleLabel_effic, trackingParticleLabel_fakes, trackLabel_algoA, trackLabel_algoB, UseAssociators, UseVertex, vertexLabel_algoA, and vertexLabel_algoB.

 {
      // create output collection instance
     auto outputAlgoA = std::make_unique<RecoTracktoTPCollection>();
     auto outputAlgoB = std::make_unique<RecoTracktoTPCollection>();
     auto outputTP = std::make_unique<TPtoRecoTrackCollection>();
   
     // Get Inputs
     edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
     iEvent.getByToken(beamSpotLabel, recoBeamSpotHandle);
     reco::BeamSpot beamSpot = *recoBeamSpotHandle; 
   
     edm::Handle<View<reco::Track> > trackCollAlgoA;
     iEvent.getByToken(trackLabel_algoA, trackCollAlgoA);
   
     edm::Handle< View<reco::Track> > trackCollAlgoB;
     iEvent.getByToken(trackLabel_algoB, trackCollAlgoB);
   
     edm::Handle<TrackingParticleCollection> trackingParticleCollFakes;
     iEvent.getByToken(trackingParticleLabel_fakes, trackingParticleCollFakes);
   
     edm::Handle<TrackingParticleCollection> trackingParticleCollEffic;
     iEvent.getByToken(trackingParticleLabel_effic, trackingParticleCollEffic);
   
     edm::Handle<reco::VertexCollection> vertexCollAlgoA;
     edm::Handle<reco::VertexCollection> vertexCollAlgoB;
     if(UseVertex) 
     {
         iEvent.getByToken(vertexLabel_algoA, vertexCollAlgoA);
         iEvent.getByToken(vertexLabel_algoB, vertexCollAlgoB);
     }
   
     // call the associator functions:
     reco::RecoToSimCollection recSimColl_AlgoA; 
     reco::RecoToSimCollection recSimColl_AlgoB;
   
     reco::SimToRecoCollection simRecColl_AlgoA;
     reco::SimToRecoCollection simRecColl_AlgoB;
 
     if(UseAssociators)
     {
         edm::Handle<reco::TrackToTrackingParticleAssociator> theAssociator_algoA;
         iEvent.getByToken(assocLabel_algoA, theAssociator_algoA);
   
         edm::Handle<reco::TrackToTrackingParticleAssociator> theAssociator_algoB;
         iEvent.getByToken(assocLabel_algoB, theAssociator_algoB);
   
         recSimColl_AlgoA = theAssociator_algoA->associateRecoToSim(trackCollAlgoA, trackingParticleCollFakes);
         recSimColl_AlgoB = theAssociator_algoB->associateRecoToSim(trackCollAlgoB, trackingParticleCollFakes);
 
         simRecColl_AlgoA = theAssociator_algoA->associateSimToReco(trackCollAlgoA, trackingParticleCollEffic);
         simRecColl_AlgoB = theAssociator_algoB->associateSimToReco(trackCollAlgoB, trackingParticleCollEffic);
     }
     else
     {
         Handle<reco::RecoToSimCollection > recotosimCollectionH_AlgoA;
         iEvent.getByToken(associatormap_algoA_recoToSim,recotosimCollectionH_AlgoA);
         recSimColl_AlgoA  = *(recotosimCollectionH_AlgoA.product());
         
         Handle<reco::RecoToSimCollection > recotosimCollectionH_AlgoB;
         iEvent.getByToken(associatormap_algoB_recoToSim,recotosimCollectionH_AlgoB);
         recSimColl_AlgoB  = *(recotosimCollectionH_AlgoB.product());
         
         Handle<reco::SimToRecoCollection > simtorecoCollectionH_AlgoA;
         iEvent.getByToken(associatormap_algoA_simToReco, simtorecoCollectionH_AlgoA);
         simRecColl_AlgoA = *(simtorecoCollectionH_AlgoA.product());
 
         Handle<reco::SimToRecoCollection > simtorecoCollectionH_AlgoB;
         iEvent.getByToken(associatormap_algoB_simToReco, simtorecoCollectionH_AlgoB);
         simRecColl_AlgoB = *(simtorecoCollectionH_AlgoB.product());
     }
     
     // define the vector of references to trackingParticleColl associated with a given reco::Track
     std::vector<std::pair<TrackingParticleRef, double> > associatedTrackingParticles;
   
     // define the vector of references to trackColl associated with a given TrackingParticle
     std::vector<std::pair<reco::TrackBaseRef, double> > associatedRecoTracks;
   
     // Get the magnetic field data from the event (used to calculate the point of closest TrackingParticle)
     edm::ESHandle<MagneticField> theMagneticField;
     iSetup.get<IdealMagneticFieldRecord>().get(theMagneticField);
     const MagneticField *magneticField = theMagneticField.product();
     
     // fill collection algoA
     for(View<reco::Track>::size_type i = 0; i < trackCollAlgoA->size(); ++i)
     {
         // get recoTrack algo A
         reco::TrackBaseRef recoTrack(trackCollAlgoA, i);
         RecoTracktoTP  recoTracktoTP;
         recoTracktoTP.SetRecoTrack(recoTrack);
         recoTracktoTP.SetBeamSpot(beamSpot.position());
         
         // get the associated trackingParticle
         if(recSimColl_AlgoA.find(recoTrack) != recSimColl_AlgoA.end())
         {
             associatedTrackingParticles = recSimColl_AlgoA[recoTrack];
             recoTracktoTP.SetTrackingParticle( associatedTrackingParticles.begin()->first );
             recoTracktoTP.SetShared( associatedTrackingParticles.begin()->second );
             SetTrackingParticleD0Dz(associatedTrackingParticles.begin()->first, beamSpot, magneticField, recoTracktoTP);
         }
         else
         {           
             recoTracktoTP.SetTrackingParticle(TrackingParticleRef());
             recoTracktoTP.SetShared(-1.0);
         }       
         
         // get the reco primary vertex info
         if(UseVertex && !vertexCollAlgoA->empty())  
         {
             recoTracktoTP.SetRecoVertex( reco::VertexRef(vertexCollAlgoA, 0) );
         }
         else
         {
             recoTracktoTP.SetRecoVertex( reco::VertexRef() );
         }   
        
         outputAlgoA->push_back(recoTracktoTP);
     }
   
   
     // fill collection algoB
     for(reco::TrackCollection::size_type i = 0; i < trackCollAlgoB->size(); ++i)
     {
         // get recoTrack algo B
         reco::TrackBaseRef recoTrack(trackCollAlgoB, i);
         RecoTracktoTP  recoTracktoTP;
         recoTracktoTP.SetRecoTrack(recoTrack);
         recoTracktoTP.SetBeamSpot(beamSpot.position());
       
         // get the associated trackingParticle
         if(recSimColl_AlgoB.find(recoTrack) != recSimColl_AlgoB.end())
         {
             associatedTrackingParticles = recSimColl_AlgoB[recoTrack];
             recoTracktoTP.SetTrackingParticle( associatedTrackingParticles.begin()->first );
             recoTracktoTP.SetShared( associatedTrackingParticles.begin()->second );
             SetTrackingParticleD0Dz(associatedTrackingParticles.begin()->first, beamSpot, magneticField, recoTracktoTP);
         }
         else
         {           
             recoTracktoTP.SetTrackingParticle(TrackingParticleRef());
             recoTracktoTP.SetShared(-1.0);
         }       
         
         // get the reco primary vertex info
         if(UseVertex && !vertexCollAlgoB->empty())  
         {
             recoTracktoTP.SetRecoVertex( reco::VertexRef(vertexCollAlgoB, 0) );
         }
         else
         {
             recoTracktoTP.SetRecoVertex( reco::VertexRef() );
         }
         
         outputAlgoB->push_back(recoTracktoTP);
     }
   
   
     for(TrackingParticleCollection::size_type i = 0; i < trackingParticleCollEffic->size(); ++i)
     {
         // initialize the trackingParticle (sim) info
         TrackingParticleRef tparticle(trackingParticleCollEffic, i);
         TPtoRecoTrack  tptoRecoTrack;
         tptoRecoTrack.SetBeamSpot(beamSpot.position());
         tptoRecoTrack.SetTrackingParticle(tparticle);
         SetTrackingParticleD0Dz(tparticle, beamSpot,  magneticField, tptoRecoTrack);
       
         // get the assocated recoTrack algoA
         if(simRecColl_AlgoA.find(tparticle) != simRecColl_AlgoA.end())
         {
             associatedRecoTracks = simRecColl_AlgoA[tparticle];
             tptoRecoTrack.SetRecoTrack_AlgoA(associatedRecoTracks.begin()->first );
             tptoRecoTrack.SetShared_AlgoA(associatedRecoTracks.begin()->second );
         }
         else
         {
             tptoRecoTrack.SetRecoTrack_AlgoA(reco::TrackBaseRef());
             tptoRecoTrack.SetShared_AlgoA(-1.0);
         }
         
         // get the recoVertex algo A
         if(UseVertex && !vertexCollAlgoA->empty())
         {
             tptoRecoTrack.SetRecoVertex_AlgoA( reco::VertexRef(vertexCollAlgoA, 0) );
         }
         else
         {
             tptoRecoTrack.SetRecoVertex_AlgoA( reco::VertexRef() );
         }
       
         // get the assocated recoTrack algoB
         if(simRecColl_AlgoB.find(tparticle) != simRecColl_AlgoB.end())
         {
             associatedRecoTracks = simRecColl_AlgoB[tparticle];
             tptoRecoTrack.SetRecoTrack_AlgoB(associatedRecoTracks.begin()->first );
             tptoRecoTrack.SetShared_AlgoB(associatedRecoTracks.begin()->second );
         }
         else
         {
             tptoRecoTrack.SetRecoTrack_AlgoB(reco::TrackBaseRef());
             tptoRecoTrack.SetShared_AlgoB(-1.0);
         }
         // get the recoVertex algo B
         if(UseVertex && !vertexCollAlgoB->empty())
         {
             tptoRecoTrack.SetRecoVertex_AlgoB( reco::VertexRef(vertexCollAlgoB, 0) );
         }
         else
         {
             tptoRecoTrack.SetRecoVertex_AlgoB( reco::VertexRef() );
         }
       
         outputTP->push_back(tptoRecoTrack);
     }
 
 
     // put the collection in the event record
     iEvent.put(std::move(outputAlgoA), "AlgoA");
     iEvent.put(std::move(outputAlgoB), "AlgoB");
     iEvent.put(std::move(outputTP), "TP");
 }

void TrackAlgoCompareUtil::SetTrackingParticleD0Dz	(	TrackingParticleRef	tp,
		const reco::BeamSpot &	bs,
		const MagneticField *	bf,
		TPtoRecoTrack &	TPRT
	)		const

private

Definition at line 272 of file TrackAlgoCompareUtil.cc.

References TrajectoryStateClosestToBeamLine::isValid(), FreeTrajectoryState::momentum(), AlCaHLTBitMon_ParallelJobs::p, FreeTrajectoryState::position(), TPtoRecoTrack::SetTrackingParticleMomentumPCA(), TPtoRecoTrack::SetTrackingParticlePCA(), and TrajectoryStateClosestToBeamLine::trackStateAtPCA().

Referenced by produce().

 {
     GlobalPoint trackingParticleVertex( tp->vertex().x(), tp->vertex().y(), tp->vertex().z() );
     GlobalVector trackingParticleP3(tp->g4Track_begin()->momentum().x(),
                                     tp->g4Track_begin()->momentum().y(),
                                     tp->g4Track_begin()->momentum().z() );
     TrackCharge trackingParticleCharge(tp->charge());
   
     FreeTrajectoryState ftsAtProduction( trackingParticleVertex, trackingParticleP3, trackingParticleCharge, bf );
     TSCBLBuilderNoMaterial tscblBuilder;
     TrajectoryStateClosestToBeamLine tsAtClosestApproach = tscblBuilder(ftsAtProduction, bs);  //as in TrackProducerAlgorithm
   
     if(tsAtClosestApproach.isValid())
     {
         GlobalPoint v1 = tsAtClosestApproach.trackStateAtPCA().position();
         GlobalVector p = tsAtClosestApproach.trackStateAtPCA().momentum();
 
         TPRT.SetTrackingParticleMomentumPCA(p);
         TPRT.SetTrackingParticlePCA(v1);
     }
     else
     {
         TPRT.SetTrackingParticleMomentumPCA(GlobalVector(-9999.0, -9999.0, -9999.0));
         TPRT.SetTrackingParticlePCA(GlobalPoint(-9999.0, -9999.0, -9999.0));
     }
 }

void TrackAlgoCompareUtil::SetTrackingParticleD0Dz	(	TrackingParticleRef	tp,
		const reco::BeamSpot &	bs,
		const MagneticField *	bf,
		RecoTracktoTP &	RTTP
	)		const

private

Definition at line 300 of file TrackAlgoCompareUtil.cc.

References TrajectoryStateClosestToBeamLine::isValid(), FreeTrajectoryState::momentum(), AlCaHLTBitMon_ParallelJobs::p, FreeTrajectoryState::position(), RecoTracktoTP::SetTrackingParticleMomentumPCA(), RecoTracktoTP::SetTrackingParticlePCA(), and TrajectoryStateClosestToBeamLine::trackStateAtPCA().

 {
     GlobalPoint trackingParticleVertex( tp->vertex().x(), tp->vertex().y(), tp->vertex().z() );
     GlobalVector trackingParticleP3(tp->g4Track_begin()->momentum().x(),
                                     tp->g4Track_begin()->momentum().y(),
                                     tp->g4Track_begin()->momentum().z() );
     TrackCharge trackingParticleCharge(tp->charge());
 
     FreeTrajectoryState ftsAtProduction( trackingParticleVertex, trackingParticleP3, trackingParticleCharge, bf );
     TSCBLBuilderNoMaterial tscblBuilder;
     TrajectoryStateClosestToBeamLine tsAtClosestApproach = tscblBuilder(ftsAtProduction, bs);  //as in TrackProducerAlgorithm
 
     if(tsAtClosestApproach.isValid())
     {
         GlobalPoint v1 = tsAtClosestApproach.trackStateAtPCA().position();
         GlobalVector p = tsAtClosestApproach.trackStateAtPCA().momentum();
 
         RTTP.SetTrackingParticleMomentumPCA(p);
         RTTP.SetTrackingParticlePCA(v1);
     }
     else
     {
         RTTP.SetTrackingParticleMomentumPCA(GlobalVector(-9999.0, -9999.0, -9999.0));
         RTTP.SetTrackingParticlePCA(GlobalPoint(-9999.0, -9999.0, -9999.0));
     }
 }

Member Data Documentation

edm::EDGetTokenT<reco::RecoToSimCollection> TrackAlgoCompareUtil::associatormap_algoA_recoToSim

private

Definition at line 75 of file TrackAlgoCompareUtil.h.

Referenced by produce(), and TrackAlgoCompareUtil().

edm::EDGetTokenT<reco::SimToRecoCollection> TrackAlgoCompareUtil::associatormap_algoA_simToReco

private

Definition at line 77 of file TrackAlgoCompareUtil.h.

Referenced by produce(), and TrackAlgoCompareUtil().

edm::EDGetTokenT<reco::RecoToSimCollection> TrackAlgoCompareUtil::associatormap_algoB_recoToSim

private

Definition at line 76 of file TrackAlgoCompareUtil.h.

Referenced by produce(), and TrackAlgoCompareUtil().

edm::EDGetTokenT<reco::SimToRecoCollection> TrackAlgoCompareUtil::associatormap_algoB_simToReco

private

Definition at line 78 of file TrackAlgoCompareUtil.h.

Referenced by produce(), and TrackAlgoCompareUtil().

edm::EDGetTokenT<reco::TrackToTrackingParticleAssociator> TrackAlgoCompareUtil::assocLabel_algoA

private

Definition at line 79 of file TrackAlgoCompareUtil.h.

Referenced by produce(), and TrackAlgoCompareUtil().

edm::EDGetTokenT<reco::TrackToTrackingParticleAssociator> TrackAlgoCompareUtil::assocLabel_algoB

private

Definition at line 80 of file TrackAlgoCompareUtil.h.

Referenced by produce(), and TrackAlgoCompareUtil().

edm::EDGetTokenT<reco::BeamSpot> TrackAlgoCompareUtil::beamSpotLabel

private

Definition at line 74 of file TrackAlgoCompareUtil.h.

Referenced by produce().

edm::EDGetTokenT<TrackingParticleCollection> TrackAlgoCompareUtil::trackingParticleLabel_effic

private

Definition at line 71 of file TrackAlgoCompareUtil.h.

Referenced by produce().

edm::EDGetTokenT<TrackingParticleCollection> TrackAlgoCompareUtil::trackingParticleLabel_fakes

private

Definition at line 70 of file TrackAlgoCompareUtil.h.

Referenced by produce().

edm::EDGetTokenT<edm::View<reco::Track> > TrackAlgoCompareUtil::trackLabel_algoA

private

Definition at line 68 of file TrackAlgoCompareUtil.h.

Referenced by produce().

edm::EDGetTokenT<edm::View<reco::Track> > TrackAlgoCompareUtil::trackLabel_algoB

private

Definition at line 69 of file TrackAlgoCompareUtil.h.

Referenced by produce().

const bool TrackAlgoCompareUtil::UseAssociators

private

Definition at line 81 of file TrackAlgoCompareUtil.h.

Referenced by produce(), and TrackAlgoCompareUtil().

const bool TrackAlgoCompareUtil::UseVertex

private

Definition at line 82 of file TrackAlgoCompareUtil.h.

Referenced by produce(), and TrackAlgoCompareUtil().

edm::EDGetTokenT<reco::VertexCollection> TrackAlgoCompareUtil::vertexLabel_algoA

private

Definition at line 72 of file TrackAlgoCompareUtil.h.

Referenced by produce(), and TrackAlgoCompareUtil().

edm::EDGetTokenT<reco::VertexCollection> TrackAlgoCompareUtil::vertexLabel_algoB

private

Definition at line 73 of file TrackAlgoCompareUtil.h.

Referenced by produce(), and TrackAlgoCompareUtil().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation