recob::DumpPFParticles Class Reference

Prints the content of all the ParticleFlow particles on screen. More...

Inheritance diagram for recob::DumpPFParticles:

Classes
struct	Config

Public Types
using	Parameters = art::EDAnalyzer::Table< Config >

Public Member Functions
	DumpPFParticles (Parameters const &config)
	Default constructor. More...
virtual void	analyze (const art::Event &evt) override
	Does the printing. More...

Private Member Functions
void	MakePFParticleGraph (art::Event const &event, art::ValidHandle< std::vector< recob::PFParticle >> const &handle) const

Static Private Member Functions
static std::string	DotFileName (art::EventID const &evtID, art::Provenance const &prodInfo)

Private Attributes
art::InputTag	fInputTag
	input tag of the PFParticle product More...
std::string	fOutputCategory
	category for LogInfo output More...
bool	fPrintHexFloats
	whether to print floats in base 16 More...
unsigned int	fMaxDepth
	maximum generation to print (0: only primaries) More...
bool	fMakeEventGraphs
	whether to create one DOT file per event More...

Detailed Description

Prints the content of all the ParticleFlow particles on screen.

This analyser prints the content of all the ParticleFlow particles into the LogInfo/LogVerbatim stream.

Configuration parameters

• PFModuleLabel (art::InputTag, required): label of the producer used to create the recob::PFParticle collection to be dumped
• OutputCategory (string, default: "DumpPFParticles"): the category used for the output (useful for filtering)
• PrintHexFloats (boolean, default: false): print all the floating point numbers in base 16
• MaxDepth (unsigned int, optional): if specified, at most this number of particle generations will be printed; 1 means printing only primaries and their daughters, 0 only primaries. If not specified, no limit will be applied. This is useful for buggy PFParticles with circular references.
• MakeParticleGraphs (boolean, default: false): creates a DOT file for each event, with a graph of PFParticle relations; each file is named as: ProcessName_ModuleLabel_InstanceName_Run::_Subrun::_Event::_particles.dot, where the the input label elements refer to the data product being plotted.

Particle connection graphs

When MakeParticleGraphs configuration option is activated, a file is created for each event, that contains the particle flow tree in GraphViz format. The GraphViz dot command can be used to render it into a PDF, SVG, EPS or one of the many supported bitmap formats. The typical command to use is:

dot -Tpdf -oPMTrk.pdf PMTrk.dot

A bash command to convert all files into a OutputFormat format:

OutputFormat='pdf'
for DotFile in *.dot ; do
  OutputFile="${DotFile%.dot}.${OutputFormat}"
  [[ "$OutputFile" -ot "$DotFile" ]] || continue # up to date already
  echo "${DotFile} => ${OutputFile} ..."
  dot -T"$OutputFormat" -o"$OutputFile" "$DotFile" || break
done

which will also skip files already converted.

The output shows one cell ("node") per particle. The format of the node follows these prescriptions:

• the label has the particle ID number prepended by a hash character (#)
• if the particle has a PDG ID, that also appears in the label (either the name of the corresponding particle, or, if unknown, just the PDG ID number)
• if the particle is primary, it is rendered in bold font
• if the particle is referred by other particles, but it is not present ("ghost particle"), its border is red and dashed

The relations between particles in the flow are represented by connecting lines ("edges"). Connection information is redundant: the parent particle should have the daughter in the daughter list, and the daughter should have the parent particle referenced as such. Since the connection is usually from two sources, there are usually two arrow heads, each one close to the particle that provides information on that connection; all arrow heads point from parent to daughter.

• when the information of daughter and parent is consistent, a black line with two arrow heads both pointing to the daughter is shown
• when the parent is ghost, the arrow head close to the daughter is hollow; ghost particles have no arrow heads close to them
• when the daughter is ghost, the arrow head close to the parent is hollow; ghost particles have no arrow heads close to them

If you are trying to interpret an existing diagram, the following list is more direct to the point. Nodes: represent particles (see above for the label content)

• bold label: primary particle
• red, dashed border: "ghost particle" (missing but referenced by others)
• other: just a particle

Connecting lines ("edges"):

• all arrow heads point from parent to daughter
• black with two full arrow heads: regular parent to daughter
• black with a single inward empty arrow head: the particle close to the arrow claims the particle pointed by the arrow as a daughter, but there is no information on that daughter (ghost daughter)
• black with a single outward empty arrow head: the particle at the tip of the arrow claims to be daughter of the other particle, but there is no information on that parent (ghost parent)
• red, outward arrow: the daughter (at the tip of the only arrow) claims the other particle as parent, but that parent does not recognise it as daughter
• orange, inward arrow: the parent (close to the only arrow head) claims the other particle as daughter, but that daighter does not recognise it as parent

Definition at line 131 of file DumpPFParticles_module.cc.

Member Typedef Documentation

using recob::DumpPFParticles::Parameters = art::EDAnalyzer::Table<Config>

Definition at line 168 of file DumpPFParticles_module.cc.

Constructor & Destructor Documentation

recob::DumpPFParticles::DumpPFParticles ( Parameters const &  config )

explicit

Default constructor.

Definition at line 1118 of file DumpPFParticles_module.cc.

    : EDAnalyzer(config)
    , fInputTag(config().PFModuleLabel())
    , fOutputCategory(config().OutputCategory())
    , fPrintHexFloats(config().PrintHexFloats())
    , fMaxDepth(std::numeric_limits<unsigned int>::max())
    , fMakeEventGraphs(config().MakeParticleGraphs())
    {
      // here we are handling the optional configuration key as it had just a
      // default value
      if (!config().MaxDepth(fMaxDepth))
        fMaxDepth = std::numeric_limits<unsigned int>::max();
    }

Member Function Documentation

void recob::DumpPFParticles::analyze ( const art::Event &  evt )

overridevirtual

Does the printing.

Definition at line 1174 of file DumpPFParticles_module.cc.

                                                   {

    //
    // collect all the available information
    //
    // fetch the data to be dumped on screen
    art::ValidHandle<std::vector<recob::PFParticle>> PFParticles
      = evt.getValidHandle<std::vector<recob::PFParticle>>(fInputTag);

    if (fMakeEventGraphs)
      MakePFParticleGraph(evt, PFParticles);

    art::FindOne<recob::Vertex> const ParticleVertices
      (PFParticles, evt, fInputTag);
    art::FindMany<recob::Track> const ParticleTracks
      (PFParticles, evt, fInputTag);
    art::FindMany<recob::Cluster> const ParticleClusters
      (PFParticles, evt, fInputTag);
    art::FindMany<recob::Seed> const ParticleSeeds
      (PFParticles, evt, fInputTag);
    art::FindMany<recob::SpacePoint> const ParticleSpacePoints
      (PFParticles, evt, fInputTag);
    art::FindMany<recob::PCAxis> const ParticlePCAxes
      (PFParticles, evt, fInputTag);

    size_t const nParticles = PFParticles->size();
    mf::LogVerbatim(fOutputCategory) << "Event " << evt.id()
      << " contains " << nParticles << " particles from '"
      << fInputTag.encode() << "'";

    // prepare the dumper
    ParticleDumper::PrintOptions_t options;
    options.hexFloats = fPrintHexFloats;
    options.maxDepth = fMaxDepth;
    options.streamName = fOutputCategory;
    ParticleDumper dumper(*PFParticles, options);
    if (ParticleVertices.isValid()) dumper.SetVertices(&ParticleVertices);
    else mf::LogPrint("DumpPFParticles") << "WARNING: vertex information not available";
    if (ParticleTracks.isValid()) dumper.SetTracks(&ParticleTracks);
    else mf::LogPrint("DumpPFParticles") << "WARNING: track information not available";
    if (ParticleClusters.isValid()) dumper.SetClusters(&ParticleClusters);
    else mf::LogPrint("DumpPFParticles") << "WARNING: cluster information not available";
    if (ParticleSeeds.isValid()) dumper.SetSeeds(&ParticleSeeds);
    else mf::LogPrint("DumpPFParticles") << "WARNING: seed information not avaialble";
    if (ParticleSpacePoints.isValid())
      dumper.SetSpacePoints(&ParticleSpacePoints);
    else {
      mf::LogPrint("DumpPFParticles")
        << "WARNING: space point information not available";
    }
    if (ParticlePCAxes.isValid())
      dumper.SetPCAxes(&ParticlePCAxes);
    else {
      mf::LogPrint("DumpPFParticles")
        << "WARNING: principal component axis not available";
    }
    dumper.DumpAllParticles("  ");

    mf::LogVerbatim(fOutputCategory) << "\n"; // two empty lines

  } // DumpPFParticles::analyze()

std::string recob::DumpPFParticles::DotFileName ( art::EventID const &  evtID, art::Provenance const &  prodInfo  )

staticprivate

Definition at line 1135 of file DumpPFParticles_module.cc.

  {
    return prodInfo.processName()
      + '_' + prodInfo.moduleLabel()
      + '_' + prodInfo.productInstanceName()
      + "_Run" + std::to_string(evtID.run())
      + "_Subrun" + std::to_string(evtID.subRun())
      + "_Event" + std::to_string(evtID.event())
      + "_particles.dot";
  } // DumpPFParticles::DotFileName()

void recob::DumpPFParticles::MakePFParticleGraph ( art::Event const &  event, art::ValidHandle< std::vector< recob::PFParticle >> const &  handle  ) const

private

Definition at line 1148 of file DumpPFParticles_module.cc.

          {
    art::EventID const eventID = event.id();
    std::string fileName = DotFileName(eventID, *(handle.provenance()));
    std::ofstream outFile(fileName); // overwrite by default

    outFile
      <<   "// " << fileName
      << "\n// "
      << "\n// Created for run " << eventID.run()
               << " subrun " << eventID.subRun()
               << " event " << eventID.event()
      << "\n// "
      << "\n// dump of " << handle->size() << " particles"
      << "\n// "
      << std::endl;

    PFParticleGraphMaker graphMaker;
    graphMaker.MakeGraph(outFile, *handle);

  } // DumpPFParticles::MakePFParticleGraph()

Member Data Documentation

art::InputTag recob::DumpPFParticles::fInputTag

private

input tag of the PFParticle product

Definition at line 178 of file DumpPFParticles_module.cc.

bool recob::DumpPFParticles::fMakeEventGraphs

private

whether to create one DOT file per event

Definition at line 182 of file DumpPFParticles_module.cc.

unsigned int recob::DumpPFParticles::fMaxDepth

private

maximum generation to print (0: only primaries)

Definition at line 181 of file DumpPFParticles_module.cc.

std::string recob::DumpPFParticles::fOutputCategory

private

category for LogInfo output

Definition at line 179 of file DumpPFParticles_module.cc.

bool recob::DumpPFParticles::fPrintHexFloats

private

whether to print floats in base 16

Definition at line 180 of file DumpPFParticles_module.cc.

---

The documentation for this class was generated from the following file:

• DumpPFParticles_module.cc