Inheritance diagram for TrackHistoryAnalyzer:

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

	~TrackHistoryAnalyzer ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

virtual void	beginJob ()

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

edm::ESHandle< ParticleDataTable >	pdt_

std::size_t	totalTracks_

edm::InputTag	trackProducer_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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

CurrentProcessingContext const *	currentContext () const

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 35 of file TrackHistoryAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 74 of file TrackHistoryAnalyzer.cc.

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

                                                                         : classifier_(config)
 {
     trackProducer_ = config.getUntrackedParameter<edm::InputTag> ( "trackProducer" );
 }

TrackHistoryAnalyzer::~TrackHistoryAnalyzer ( )

Definition at line 80 of file TrackHistoryAnalyzer.cc.

80 { }

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 83 of file TrackHistoryAnalyzer.cc.

References classifier_, gather_cfg::cout, TrackClassifier::evaluate(), TrackCategories::Fake, genParticleCandidates2GenParticles_cfi::genParticles, HistoryBase::genParticleTrail(), HistoryBase::genVertexTrail(), TrackClassifier::history(), getHLTprescales::index, TrackCategories::is(), TrackClassifier::newEvent(), particleString(), benchmark_cfg::pdgId, HistoryBase::simParticleTrail(), HistoryBase::simVertexTrail(), trackProducer_, and vertexString().

 {
     // Track collection
     edm::Handle<edm::View<reco::Track> > trackCollection;
     event.getByLabel(trackProducer_, trackCollection);
 
     // Set the classifier for a new event
     classifier_.newEvent(event, setup);
 
     // Get a constant reference to the track history associated to the classifier
     TrackHistory const & tracer = classifier_.history();
 
     // Loop over the track collection.
     for (std::size_t index = 0; index < trackCollection->size(); index++)
     {
         std::cout << std::endl << "History for track #" << index << " : " << std::endl;
 
         // Classify the track and detect for fakes
         if ( ! classifier_.evaluate( reco::TrackBaseRef(trackCollection, index) ).is(TrackClassifier::Fake) )
         {
             // Get the list of TrackingParticles associated to
             TrackHistory::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
             TrackHistory::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
             TrackHistory::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
             TrackHistory::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 track" << std::endl;
 
         std::cout << "  track categories : " << classifier_;
         std::cout << std::endl;
     }
 }

void TrackHistoryAnalyzer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 182 of file TrackHistoryAnalyzer.cc.

References totalTracks_.

 {
     totalTracks_ = 0;
 }

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

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 175 of file TrackHistoryAnalyzer.cc.

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

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

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

private

Definition at line 188 of file TrackHistoryAnalyzer.cc.

References RecoTau_DiTaus_pt_20-420_cfg::ParticleID, benchmark_cfg::pdgId, pdt_, and evf::utils::pid.

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 TrackHistoryAnalyzer::vertexString	(	const TrackingParticleRefVector &	in,
		const TrackingParticleRefVector &	out
	)		const

private

Definition at line 211 of file TrackHistoryAnalyzer.cc.

References j, RecoTau_DiTaus_pt_20-420_cfg::ParticleID, benchmark_cfg::pdgId, pdt_, evf::utils::pid, 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 TrackHistoryAnalyzer::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 268 of file TrackHistoryAnalyzer.cc.

References recoMuon::in, j, dbtoconf::out, RecoTau_DiTaus_pt_20-420_cfg::ParticleID, pdt_, and evf::utils::pid.

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