Inheritance diagram for VertexHistoryAnalyzer:

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

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

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

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	beginJob () override

virtual void	beginRun (const edm::Run &, const edm::EventSetup &)

std::string	particleString (int) const

std::string	vertexString (const TrackingParticleRefVector &, const TrackingParticleRefVector &) const

std::string	vertexString (HepMC::GenVertex::particles_in_const_iterator, HepMC::GenVertex::particles_in_const_iterator, HepMC::GenVertex::particles_out_const_iterator, HepMC::GenVertex::particles_out_const_iterator) const

Private Attributes
VertexClassifier	classifier_

edm::ESHandle< ParticleDataTable >	pdt_

edm::InputTag	vertexProducer_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
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)

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 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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

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 34 of file VertexHistoryAnalyzer.cc.

Constructor & Destructor Documentation

VertexHistoryAnalyzer::VertexHistoryAnalyzer ( const edm::ParameterSet & config )

explicit

Definition at line 61 of file VertexHistoryAnalyzer.cc.

References edm::ParameterSet::getUntrackedParameter(), and vertexProducer_.

     : classifier_(config, consumesCollector()) {
   vertexProducer_ = config.getUntrackedParameter<edm::InputTag>("vertexProducer");
   consumes<edm::View<reco::Vertex>>(vertexProducer_);
 }

Member Function Documentation

void VertexHistoryAnalyzer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 67 of file VertexHistoryAnalyzer.cc.

References classifier_, gather_cfg::cout, VertexClassifier::evaluate(), VertexCategories::Fake, genParticleCandidates2GenParticles_cfi::genParticles, HistoryBase::genParticleTrail(), HistoryBase::genVertexTrail(), VertexClassifier::history(), VertexCategories::is(), VertexClassifier::newEvent(), particleString(), HistoryBase::simParticleTrail(), HistoryBase::simVertexTrail(), GoodVertex_cfg::vertexCollection, vertexProducer_, and vertexString().

                                                                                      {
   // Set the classifier for a new event
   classifier_.newEvent(event, setup);
 
   // Vertex collection
   edm::Handle<edm::View<reco::Vertex>> vertexCollection;
   event.getByLabel(vertexProducer_, vertexCollection);
 
   // Get a constant reference to the track history associated to the classifier
   VertexHistory const &tracer = classifier_.history();
 
   // Loop over the track collection.
   for (std::size_t index = 0; index < vertexCollection->size(); index++) {
     std::cout << std::endl << "History for vertex #" << index << " : " << std::endl;
 
     // Classify the track and detect for fakes
     if (!classifier_.evaluate(reco::VertexBaseRef(vertexCollection, index)).is(VertexClassifier::Fake)) {
       // Get the list of TrackingParticles associated to
       VertexHistory::SimParticleTrail simParticles(tracer.simParticleTrail());
 
       // Loop over all simParticles
       for (std::size_t hindex = 0; hindex < simParticles.size(); hindex++) {
         std::cout << "  simParticles [" << hindex << "] : " << particleString(simParticles[hindex]->pdgId())
                   << std::endl;
       }
 
       // Get the list of TrackingVertexes associated to
       VertexHistory::SimVertexTrail simVertexes(tracer.simVertexTrail());
 
       // Loop over all simVertexes
       if (!simVertexes.empty()) {
         for (std::size_t hindex = 0; hindex < simVertexes.size(); hindex++) {
           std::cout << "  simVertex    [" << hindex << "] : "
                     << vertexString(simVertexes[hindex]->sourceTracks(), simVertexes[hindex]->daughterTracks())
                     << std::endl;
         }
       } else
         std::cout << "  simVertex no found" << std::endl;
 
       // Get the list of GenParticles associated to
       VertexHistory::GenParticleTrail genParticles(tracer.genParticleTrail());
 
       // Loop over all genParticles
       for (std::size_t hindex = 0; hindex < genParticles.size(); hindex++) {
         std::cout << "  genParticles [" << hindex << "] : " << particleString(genParticles[hindex]->pdg_id())
                   << std::endl;
       }
 
       // Get the list of TrackingVertexes associated to
       VertexHistory::GenVertexTrail genVertexes(tracer.genVertexTrail());
 
       // Loop over all simVertexes
       if (!genVertexes.empty()) {
         for (std::size_t hindex = 0; hindex < genVertexes.size(); hindex++) {
           std::cout << "  genVertex    [" << hindex << "] : "
                     << vertexString(genVertexes[hindex]->particles_in_const_begin(),
                                     genVertexes[hindex]->particles_in_const_end(),
                                     genVertexes[hindex]->particles_out_const_begin(),
                                     genVertexes[hindex]->particles_out_const_end())
                     << std::endl;
         }
       } else
         std::cout << "  genVertex no found" << std::endl;
     } else
       std::cout << "  fake vertex" << std::endl;
 
     std::cout << "  vertex categories : " << classifier_;
     std::cout << std::endl;
   }
 }

void VertexHistoryAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 143 of file VertexHistoryAnalyzer.cc.

143 {}

void VertexHistoryAnalyzer::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

privatevirtual

Definition at line 138 of file VertexHistoryAnalyzer.cc.

References edm::EventSetup::getData(), and pdt_.

                                                                                   {
   // Get the particles table.
   setup.getData(pdt_);
 }

std::string VertexHistoryAnalyzer::particleString ( int pdgId ) const

private

Definition at line 145 of file VertexHistoryAnalyzer.cc.

References HLT_FULL_cff::particleType, and pdt_.

Referenced by analyze().

                                                                {
   ParticleData const *pid;
 
   std::ostringstream vDescription;
 
   HepPDT::ParticleID particleType(pdgId);
 
   if (particleType.isValid()) {
     pid = pdt_->particle(particleType);
     if (pid)
       vDescription << pid->name();
     else
       vDescription << pdgId;
   } else
     vDescription << pdgId;
 
   return vDescription.str();
 }

std::string VertexHistoryAnalyzer::vertexString	(	const TrackingParticleRefVector &	in,
		const TrackingParticleRefVector &	out
	)		const

private

Definition at line 164 of file VertexHistoryAnalyzer.cc.

References dqmiolumiharvest::j, HLT_FULL_cff::particleType, pdt_, and edm::RefVector< C, T, F >::size().

Referenced by analyze().

                                                                                             {
   ParticleData const *pid;
 
   std::ostringstream vDescription;
 
   for (std::size_t j = 0; j < in.size(); j++) {
     if (!j)
       vDescription << "(";
 
     HepPDT::ParticleID particleType(in[j]->pdgId());
 
     if (particleType.isValid()) {
       pid = pdt_->particle(particleType);
       if (pid)
         vDescription << pid->name();
       else
         vDescription << in[j]->pdgId();
     } else
       vDescription << in[j]->pdgId();
 
     if (j == in.size() - 1)
       vDescription << ")";
     else
       vDescription << ",";
   }
 
   vDescription << "->";
 
   for (std::size_t j = 0; j < out.size(); j++) {
     if (!j)
       vDescription << "(";
 
     HepPDT::ParticleID particleType(out[j]->pdgId());
 
     if (particleType.isValid()) {
       pid = pdt_->particle(particleType);
       if (pid)
         vDescription << pid->name();
       else
         vDescription << out[j]->pdgId();
     } else
       vDescription << out[j]->pdgId();
 
     if (j == out.size() - 1)
       vDescription << ")";
     else
       vDescription << ",";
   }
 
   return vDescription.str();
 }

std::string VertexHistoryAnalyzer::vertexString	(	HepMC::GenVertex::particles_in_const_iterator	in_begin,
		HepMC::GenVertex::particles_in_const_iterator	in_end,
		HepMC::GenVertex::particles_out_const_iterator	out_begin,
		HepMC::GenVertex::particles_out_const_iterator	out_end
	)		const

private

Definition at line 217 of file VertexHistoryAnalyzer.cc.

References recoMuon::in, dqmiolumiharvest::j, submitPVResolutionJobs::out, HLT_FULL_cff::particleType, and pdt_.

                                                                                                           {
   ParticleData const *pid;
 
   std::ostringstream vDescription;
 
   std::size_t j = 0;
 
   HepMC::GenVertex::particles_in_const_iterator in, itmp;
 
   for (in = in_begin; in != in_end; in++, j++) {
     if (!j)
       vDescription << "(";
 
     HepPDT::ParticleID particleType((*in)->pdg_id());
 
     if (particleType.isValid()) {
       pid = pdt_->particle(particleType);
       if (pid)
         vDescription << pid->name();
       else
         vDescription << (*in)->pdg_id();
     } else
       vDescription << (*in)->pdg_id();
 
     itmp = in;
 
     if (++itmp == in_end)
       vDescription << ")";
     else
       vDescription << ",";
   }
 
   vDescription << "->";
   j = 0;
 
   HepMC::GenVertex::particles_out_const_iterator out, otmp;
 
   for (out = out_begin; out != out_end; out++, j++) {
     if (!j)
       vDescription << "(";
 
     HepPDT::ParticleID particleType((*out)->pdg_id());
 
     if (particleType.isValid()) {
       pid = pdt_->particle(particleType);
       if (pid)
         vDescription << pid->name();
       else
         vDescription << (*out)->pdg_id();
     } else
       vDescription << (*out)->pdg_id();
 
     otmp = out;
 
     if (++otmp == out_end)
       vDescription << ")";
     else
       vDescription << ",";
   }
 
   return vDescription.str();
 }