TrackHistoryAnalyzer Class Reference

Inheritance diagram for TrackHistoryAnalyzer:

Public Member Functions
	TrackHistoryAnalyzer (const edm::ParameterSet &)
	~TrackHistoryAnalyzer () override
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
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::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	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
TrackClassifier	classifier_
edm::ESHandle< ParticleDataTable >	pdt_
std::size_t	totalTracks_
edm::InputTag	trackProducer_

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
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 = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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...
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

Definition at line 35 of file TrackHistoryAnalyzer.cc.

Constructor & Destructor Documentation

TrackHistoryAnalyzer()

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

explicit

Definition at line 65 of file TrackHistoryAnalyzer.cc.

                                                                        : classifier_(config, consumesCollector()) {
  trackProducer_ = config.getUntrackedParameter<edm::InputTag>("trackProducer");
  consumes<edm::View<reco::Track>>(trackProducer_);
}

References trackProducer_.

~TrackHistoryAnalyzer()

TrackHistoryAnalyzer::~TrackHistoryAnalyzer ( )

override

Definition at line 70 of file TrackHistoryAnalyzer.cc.

{}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 72 of file TrackHistoryAnalyzer.cc.

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

References classifier_, gather_cfg::cout, TrackClassifier::evaluate(), TrackCategories::Fake, genParticles2HepMC_cfi::genParticles, HistoryBase::genParticleTrail(), HistoryBase::genVertexTrail(), TrackClassifier::history(), TrackCategories::is(), TrackClassifier::newEvent(), particleString(), HiggsValidation_cfi::pdg_id, EgammaValidation_cff::pdgId, singleTopDQM_cfi::setup, HistoryBase::simParticleTrail(), HistoryBase::simVertexTrail(), duplicaterechits_cfi::trackCollection, trackProducer_, and vertexString().

beginJob()

void TrackHistoryAnalyzer::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 148 of file TrackHistoryAnalyzer.cc.

{ totalTracks_ = 0; }

References totalTracks_.

beginRun()

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

privatevirtual

Definition at line 143 of file TrackHistoryAnalyzer.cc.

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

References pdt_, and singleTopDQM_cfi::setup.

particleString()

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

private

Definition at line 150 of file TrackHistoryAnalyzer.cc.

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

References LHEGenericFilter_cfi::ParticleID, PbPb_ZMuSkimMuonDPG_cff::particleType, EgammaValidation_cff::pdgId, and pdt_.

Referenced by analyze().

vertexString() [1/2]

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

private

Definition at line 169 of file TrackHistoryAnalyzer.cc.

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

References recoMuon::in, dqmiolumiharvest::j, MillePedeFileConverter_cfg::out, LHEGenericFilter_cfi::ParticleID, PbPb_ZMuSkimMuonDPG_cff::particleType, EgammaValidation_cff::pdgId, and pdt_.

Referenced by analyze().

vertexString() [2/2]

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

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

References recoMuon::in, dqmiolumiharvest::j, MillePedeFileConverter_cfg::out, LHEGenericFilter_cfi::ParticleID, PbPb_ZMuSkimMuonDPG_cff::particleType, and pdt_.

Member Data Documentation

classifier_

TrackClassifier TrackHistoryAnalyzer::classifier_

private

Definition at line 55 of file TrackHistoryAnalyzer.cc.

Referenced by analyze().

pdt_

edm::ESHandle<ParticleDataTable> TrackHistoryAnalyzer::pdt_

private

Definition at line 51 of file TrackHistoryAnalyzer.cc.

Referenced by beginRun(), particleString(), and vertexString().

totalTracks_

std::size_t TrackHistoryAnalyzer::totalTracks_

private

Definition at line 49 of file TrackHistoryAnalyzer.cc.

Referenced by beginJob().

trackProducer_

edm::InputTag TrackHistoryAnalyzer::trackProducer_

private

Definition at line 47 of file TrackHistoryAnalyzer.cc.

Referenced by analyze(), and TrackHistoryAnalyzer().