ParticleDecayDrawer Class Reference

Inheritance diagram for ParticleDecayDrawer:

Public Member Functions
	ParticleDecayDrawer (const edm::ParameterSet &)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () 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::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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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
bool	accept (const reco::Candidate &, const std::list< const reco::Candidate * > &) const
	accept candidate More...
void	analyze (const edm::Event &, const edm::EventSetup &) override
std::string	decay (const reco::Candidate &, std::list< const reco::Candidate * > &) const
bool	hasValidDaughters (const reco::Candidate &) const
	has valid daughters in the chain More...
std::string	printP4 (const reco::Candidate &) const
	print 4 momenta More...
bool	select (const reco::Candidate &) const
	select candidate More...

Private Attributes
edm::ESHandle< ParticleDataTable >	pdt_
edm::ESGetToken< ParticleDataTable, edm::DefaultRecord >	pdtToken_
bool	printP4_
	print parameters More...
bool	printPtEtaPhi_
bool	printVertex_
edm::EDGetTokenT< edm::View< reco::Candidate > >	srcToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
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 13 of file ParticleDecayDrawer.cc.

Constructor & Destructor Documentation

ParticleDecayDrawer()

ParticleDecayDrawer::ParticleDecayDrawer

(

const edm::ParameterSet &

cfg

)

Definition at line 47 of file ParticleDecayDrawer.cc.

    : srcToken_(consumes<edm::View<reco::Candidate> >(cfg.getParameter<InputTag>("src"))),
      pdtToken_(esConsumes<ParticleDataTable, edm::DefaultRecord>()),
      printP4_(cfg.getUntrackedParameter<bool>("printP4", false)),
      printPtEtaPhi_(cfg.getUntrackedParameter<bool>("printPtEtaPhi", false)),
      printVertex_(cfg.getUntrackedParameter<bool>("printVertex", false)) {}

Member Function Documentation

accept()

bool ParticleDecayDrawer::accept ( const reco::Candidate &  , const std::list< const reco::Candidate * > &    ) const

private

accept candidate

Definition at line 54 of file ParticleDecayDrawer.cc.

                                                                                                  {
  if (find(skip.begin(), skip.end(), &c) != skip.end())
    return false;
  return select(c);
}

References c, spr::find(), select(), and optionsL1T::skip.

Referenced by decay(), and esMonitoring.FDJsonServer::handle_accept().

analyze()

void ParticleDecayDrawer::analyze ( const edm::Event &  event, const edm::EventSetup &  es  )

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 70 of file ParticleDecayDrawer.cc.

                                                                          {
  pdt_ = es.getHandle(pdtToken_);
  Handle<View<Candidate> > particles;
  event.getByToken(srcToken_, particles);
  list<const Candidate *> skip;
  vector<const Candidate *> nodes, moms;
  for (View<Candidate>::const_iterator p = particles->begin(); p != particles->end(); ++p) {
    if (p->numberOfMothers() > 1) {
      if (select(*p)) {
        skip.push_back(&*p);
        nodes.push_back(&*p);
        for (size_t j = 0; j < p->numberOfMothers(); ++j) {
          const Candidate *mom = p->mother(j);
          const Candidate *grandMom;
          while ((grandMom = mom->mother()) != nullptr)
            mom = grandMom;
          if (select(*mom)) {
            moms.push_back(mom);
          }
        }
      }
    }
  }
  cout << "-- decay: --" << endl;
  if (!moms.empty()) {
    if (moms.size() > 1)
      for (size_t m = 0; m < moms.size(); ++m) {
        string dec = decay(*moms[m], skip);
        if (!dec.empty())
          cout << "{ " << dec << " } ";
      }
    else
      cout << decay(*moms[0], skip);
  }
  if (!nodes.empty()) {
    cout << "-> ";
    if (nodes.size() > 1) {
      for (size_t n = 0; n < nodes.size(); ++n) {
        skip.remove(nodes[n]);
        string dec = decay(*nodes[n], skip);
        if (!dec.empty()) {
          if (dec.find("->", 0) != string::npos)
            cout << " ( " << dec << " )";
          else
            cout << " " << dec;
        }
      }
    } else {
      skip.remove(nodes[0]);
      cout << decay(*nodes[0], skip);
    }
  }
  cout << endl;
}

References gather_cfg::cout, TauDecayModes::dec, decay(), edm::EventSetup::getHandle(), dqmiolumiharvest::j, visualization-live-secondInstance_cfg::m, reco::Candidate::mother(), dqmiodumpmetadata::n, class-composition::nodes, AlCaHLTBitMon_ParallelJobs::p, ecalTrigSettings_cff::particles, pdt_, pdtToken_, select(), optionsL1T::skip, and srcToken_.

decay()

string ParticleDecayDrawer::decay ( const reco::Candidate &  , std::list< const reco::Candidate * > &    ) const

private

Definition at line 136 of file ParticleDecayDrawer.cc.

                                                                                         {
  string out;
  if (find(skip.begin(), skip.end(), &c) != skip.end())
    return out;
  skip.push_back(&c);
 
  int id = c.pdgId();
  const ParticleData *pd = pdt_->particle(id);
  assert(pd != nullptr);
  out += (pd->name() + printP4(c));
 
  size_t validDau = 0, ndau = c.numberOfDaughters();
  for (size_t i = 0; i < ndau; ++i)
    if (accept(*c.daughter(i), skip))
      ++validDau;
  if (validDau == 0)
    return out;
 
  out += " ->";
 
  for (size_t i = 0; i < ndau; ++i) {
    const Candidate *d = c.daughter(i);
    if (accept(*d, skip)) {
      string dec = decay(*d, skip);
      if (dec.find("->", 0) != string::npos)
        out += (" ( " + dec + " )");
      else
        out += (" " + dec);
    }
  }
  return out;
}

References accept(), cms::cuda::assert(), c, ztail::d, TauDecayModes::dec, spr::find(), mps_fire::i, MillePedeFileConverter_cfg::out, pdt_, printP4(), and optionsL1T::skip.

Referenced by analyze().

hasValidDaughters()

bool ParticleDecayDrawer::hasValidDaughters ( const reco::Candidate &  c ) const

private

has valid daughters in the chain

Definition at line 62 of file ParticleDecayDrawer.cc.

                                                                        {
  size_t ndau = c.numberOfDaughters();
  for (size_t i = 0; i < ndau; ++i)
    if (select(*c.daughter(i)))
      return true;
  return false;
}

References c, mps_fire::i, and select().

printP4()

string ParticleDecayDrawer::printP4 ( const reco::Candidate &  c ) const

private

print 4 momenta

Definition at line 125 of file ParticleDecayDrawer.cc.

                                                            {
  ostringstream cout;
  if (printP4_)
    cout << " (" << c.px() << ", " << c.py() << ", " << c.pz() << "; " << c.energy() << ")";
  if (printPtEtaPhi_)
    cout << " [" << c.pt() << ": " << c.eta() << ", " << c.phi() << "]";
  if (printVertex_)
    cout << " {" << c.vx() << ", " << c.vy() << ", " << c.vz() << "}";
  return cout.str();
}

References c, gather_cfg::cout, printP4_, printPtEtaPhi_, and printVertex_.

Referenced by decay().

select()

bool ParticleDecayDrawer::select ( const reco::Candidate &  c ) const

private

select candidate

Definition at line 60 of file ParticleDecayDrawer.cc.

{ return c.status() == 3; }

References c.

Referenced by accept(), analyze(), and hasValidDaughters().

Member Data Documentation

pdt_

edm::ESHandle<ParticleDataTable> ParticleDecayDrawer::pdt_

private

Definition at line 21 of file ParticleDecayDrawer.cc.

Referenced by analyze(), and decay().

pdtToken_

edm::ESGetToken<ParticleDataTable, edm::DefaultRecord> ParticleDecayDrawer::pdtToken_

private

Definition at line 22 of file ParticleDecayDrawer.cc.

Referenced by analyze().

printP4_

bool ParticleDecayDrawer::printP4_

private

print parameters

Definition at line 24 of file ParticleDecayDrawer.cc.

Referenced by printP4().

printPtEtaPhi_

bool ParticleDecayDrawer::printPtEtaPhi_

private

Definition at line 24 of file ParticleDecayDrawer.cc.

Referenced by printP4().

printVertex_

bool ParticleDecayDrawer::printVertex_

private

Definition at line 24 of file ParticleDecayDrawer.cc.

Referenced by printP4().

srcToken_

edm::EDGetTokenT<edm::View<reco::Candidate> > ParticleDecayDrawer::srcToken_

private

Definition at line 19 of file ParticleDecayDrawer.cc.

Referenced by analyze().