Functions to dump Monte Carlo object information into a stream. More...

Functions

template<typename Stream >
void	DumpMCParticle (Stream &&out, simb::MCParticle const &particle, std::string indent, std::string firstIndent)
	Dumps the content of the specified particle in the output stream. More...

template<typename Stream >
void	DumpMCParticle (Stream &&out, simb::MCParticle const &particle, std::string indent="")
	Dumps the content of the specified particle in the output stream. More...


template<typename Stream >
void	DumpMCParticleTrajectory (Stream &&out, simb::MCTrajectory const &trajectory, unsigned int pointsPerLine, std::string indent)
	Dumps the specified particle trajectory into the output stream. More...

template<typename Stream >
void	DumpMCParticleTrajectory (Stream &&out, simb::MCTrajectory const &trajectory)
	Dumps the specified particle trajectory into the output stream. More...


template<typename Stream >
void	DumpMCNeutrino (Stream &&out, simb::MCNeutrino const &neutrino, std::string indent, std::string firstIndent)
	Dumps the content of the specified neutrino in the output stream. More...

template<typename Stream >
void	DumpMCNeutrino (Stream &&out, simb::MCNeutrino const &neutrino, std::string indent="")
	Dumps the content of the specified neutrino in the output stream. More...


template<typename Stream >
void	DumpMCTruth (Stream &&out, simb::MCTruth const &truth, unsigned int pointsPerLine, std::string indent, std::string firstIndent)
	Dumps the content of the specified MC truth in the output stream. More...

template<typename Stream >
void	DumpMCTruth (Stream &&out, simb::MCTruth const &truth, unsigned int pointsPerLine, std::string indent="")
	Dumps the content of the specified MC truth in the output stream. More...

template<typename Stream >
void	DumpMCTruth (Stream &&out, simb::MCTruth const &truth, std::string indent, std::string firstIndent)
	Dumps the content of the specified MC truth in the output stream. More...

template<typename Stream >
void	DumpMCTruth (Stream &&out, simb::MCTruth const &truth, std::string indent="")
	Dumps the content of the specified MC truth in the output stream. More...


template<typename Stream >
void	DumpGTruth (Stream &&out, simb::GTruth const &truth, std::string indent, std::string firstIndent)
	Dumps the content of the GENIE truth in the output stream. More...

template<typename Stream >
void	DumpGTruth (Stream &&out, simb::GTruth const &truth, std::string indent="")
	Dumps the content of the GENIE truth in the output stream. More...

Detailed Description

Functions to dump Monte Carlo object information into a stream.

Function Documentation

template<typename Stream >

void sim::dump::DumpGTruth	(	Stream &&	out,
		simb::GTruth const &	truth,
		std::string	indent,
		std::string	firstIndent
	)

Dumps the content of the GENIE truth in the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
truth	the truth information to be dumped
indent	base indentation string (default: none)
firstIndent	string to be used for indentation of the first line (default: as `indent`)

The indent string is prepended to every line of output, except for the first one, where firstIndent is used.

The output starts on the current line, and the last line is NOT broken.

Definition at line 377 of file MCDumpers.h.

References simb::GTruth::fDiffXsec, simb::GTruth::fFShadSystP4, simb::GTruth::fGint, simb::GTruth::fgQ2, simb::GTruth::fgq2, simb::GTruth::fGscatter, simb::GTruth::fgT, simb::GTruth::fgW, simb::GTruth::fgX, simb::GTruth::fgY, simb::GTruth::fHitNucP4, simb::GTruth::fIsCharm, simb::GTruth::fIsSeaQuark, simb::GTruth::fNumNeutron, simb::GTruth::fNumPi0, simb::GTruth::fNumPiMinus, simb::GTruth::fNumPiPlus, simb::GTruth::fNumProton, simb::GTruth::fprobability, simb::GTruth::fProbeP4, simb::GTruth::fProbePDG, simb::GTruth::fResNum, simb::GTruth::ftgtA, simb::GTruth::ftgtPDG, simb::GTruth::ftgtZ, simb::GTruth::fVertex, simb::GTruth::fweight, simb::GTruth::fXsec, and sim::ParticleName().

Referenced by sim::DumpGTruth::analyze(), DumpGTruth(), DumpMCTruth(), and evgen::GENIEGen::produce().

 {
 
   unsigned int const nCharged = truth.fNumPiPlus + truth.fNumPiMinus + truth.fNumProton;
   unsigned int const nNeutral = truth.fNumPi0 + truth.fNumNeutron;
   unsigned int const nPions = truth.fNumPiPlus + truth.fNumPiMinus + truth.fNumPi0;
   unsigned int const nNucleons = truth.fNumProton + truth.fNumNeutron;
   unsigned int const nTotalParticles = nCharged + nNeutral;
 
   out << firstIndent << "interaction code: " << truth.fGint
       << ", neutrino scattering code: " << truth.fGscatter << " at " << truth.fVertex << "\n"
       << indent << "probe: " << ParticleName(truth.fProbePDG) << " with cp=" << truth.fProbeP4
       << " hit nucleon with cp=" << truth.fHitNucP4 << " GeV"
       << " (" << (truth.fIsSeaQuark ? "" : "not a ") << "sea quark)"
       << " in target: " << ParticleName(truth.ftgtPDG) << " (Z: " << truth.ftgtZ
       << ", A: " << truth.ftgtA << ")"
       << "\n"
       << indent << "event interaction weight (genie internal): " << truth.fweight
       << ", interaction probability: " << truth.fprobability << ", cross section: " << truth.fXsec
       << ", differential cross section: " << truth.fDiffXsec << "\n"
       << indent << "particles after reaction, before FSI: " << truth.fNumPiPlus << " pi+"
       << ", " << truth.fNumPiMinus << " pi-"
       << ", " << truth.fNumPi0 << " pi0"
       << ", " << truth.fNumProton << " p/pbar"
       << ", " << truth.fNumNeutron << " n/nbar"
       << "\n"
       << indent << "  total " << nTotalParticles
       << " particles after reaction before FSI"
          ": "
       << nCharged << "/" << nNeutral
       << " charged/neutral"
          ", "
       << nPions << " pions, " << nNucleons << " nucleons"
       << "\n"
       << indent << "process " << (truth.fIsCharm ? "with" : "without") << " charmed hadron";
   if (truth.fResNum == -1)
     out << ", no resonance";
   else
     out << ", resonance: #" << truth.fResNum;
   out << "\n"
       << indent << "internal (on shell) genie kinematics: Q^2: " << truth.fgQ2 << " GeV^2"
       << " q^2: " << truth.fgq2 << " GeV^2"
       << ", w: " << truth.fgW << " GeV^2"
       << ", t: " << truth.fgT << " GeV^2"
       << ", x: " << truth.fgX << ", y: " << truth.fgY << "\n"
       << indent << "FShadSyst: " << truth.fFShadSystP4;
 
 } // sim::DumpGTruth::DumpTruth()

template<typename Stream >

void sim::dump::DumpGTruth	(	Stream &&	out,
		simb::GTruth const &	truth,
		std::string	indent = `""`
	)

Dumps the content of the GENIE truth in the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
truth	the truth information to be dumped
indent	base indentation string (default: none)
firstIndent	string to be used for indentation of the first line (default: as `indent`)

The indent string is prepended to every line of output, except for the first one, where firstIndent is used.

The output starts on the current line, and the last line is NOT broken.

Definition at line 209 of file MCDumpers.h.

References DumpGTruth().

     {
       DumpGTruth(std::forward<Stream>(out), truth, indent, indent);
     }

template<typename Stream >

void sim::dump::DumpMCNeutrino	(	Stream &&	out,
		simb::MCNeutrino const &	neutrino,
		std::string	indent,
		std::string	firstIndent
	)

Dumps the content of the specified neutrino in the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
neutrino	the neutrino to be dumped
indent	base indentation string (default: none)
firstIndent	string to be used for indentation of the first line (default: as `indent`)

The indent string is prepended to every line of output, except for the first one, where firstIndent is used.

The output starts on the current line, and the last line is NOT broken.

Definition at line 310 of file MCDumpers.h.

Referenced by DumpMCNeutrino(), DumpMCParticleTrajectory(), and DumpMCTruth().

 {
 
   out << firstIndent << "from " << TruthCCNCname(nu.CCNC()) << ", "
       << TruthInteractionTypeName(nu.InteractionType()) << ", mode: " << nu.Mode() << " ("
       << TruthReactionMode(nu.Mode()) << ")" << '\n'
       << indent << "target: " << nu.Target() << " (" << ParticleName(nu.Target()) << ")";
   if (nu.HitNuc() != 0) {
     out << ", hit nucleon: " << nu.HitNuc() << " (" << ParticleName(nu.HitNuc()) << ")";
   }
   if (nu.HitQuark() != 0) {
     out << ", hit quark: " << nu.HitQuark() << " (" << ParticleName(nu.HitQuark()) << ")";
   }
   out << '\n'
       << indent << "x=" << nu.X() << " y=" << nu.Y() << " w=" << nu.W() << " Q^2=" << nu.QSqr()
       << " GeV^2; theta=" << nu.Theta() << " rad pT=" << nu.Pt() << " GeV/c";
   out << '\n' << indent << "neutrino: ";
   DumpMCParticle(std::forward<Stream>(out), nu.Nu(), indent + "  ", "");
   out << '\n' << indent << "outgoing lepton: ";
   DumpMCParticle(std::forward<Stream>(out), nu.Lepton(), indent + "  ", "");
 
 } // sim::dump::DumpMCNeutrino()

template<typename Stream >

void sim::dump::DumpMCNeutrino	(	Stream &&	out,
		simb::MCNeutrino const &	neutrino,
		std::string	indent = `""`
	)

Dumps the content of the specified neutrino in the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
neutrino	the neutrino to be dumped
indent	base indentation string (default: none)
firstIndent	string to be used for indentation of the first line (default: as `indent`)

The indent string is prepended to every line of output, except for the first one, where firstIndent is used.

The output starts on the current line, and the last line is NOT broken.

Definition at line 122 of file MCDumpers.h.

References DumpMCNeutrino(), and DumpMCTruth().

     {
       DumpMCNeutrino(std::forward<Stream>(out), neutrino, indent, indent);
     }

template<typename Stream >

void sim::dump::DumpMCParticle	(	Stream &&	out,
		simb::MCParticle const &	particle,
		std::string	indent,
		std::string	firstIndent
	)

Dumps the content of the specified particle in the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
particle	the particle to be dumped
indent	base indentation string (default: none)
firstIndent	indentation of first line (default: as `indent`)

The indent string is prepended to every line of output, except for the first one, where firstIndent is used.

The output starts on the current line, and the last line is NOT broken.

Definition at line 226 of file MCDumpers.h.

Referenced by DumpMCNeutrino(), DumpMCParticle(), sim::DumpMCParticles::DumpMCParticle(), DumpMCTruth(), and larg4::LArG4::produce().

 {
   out << firstIndent << "ID=" << particle.TrackId() << ": " << ParticleName(particle.PdgCode())
       << " mass=" << particle.Mass() << " GeV/c2 "
       << " status=" << particle.StatusCode() << " (" << ParticleStatusName(particle.StatusCode())
       << ")";
   if (particle.Weight() != 1.0) out << " weight=" << particle.Weight();
   if (particle.Rescatter() != simb::MCParticle::s_uninitialized) {
     out << " rescattered (" << particle.Rescatter() << ")";
   }
   out << "; generator vertex " << particle.GetGvtx();
   out << "\n"
       << indent << "created via " << (particle.Process().empty() ? "magics" : particle.Process());
   if (particle.Mother() == 0)
     out << " by the gods";
   else
     out << " from ID=" << particle.Mother();
 
   const unsigned int nDaughters = particle.NumberDaughters();
   const unsigned int nPoints = particle.NumberTrajectoryPoints();
   if (nPoints > 0) {
     TLorentzVector const& start = particle.Position();
     TLorentzVector const& start_mom = particle.Momentum();
     out << " at " << start << " cm with momentum " << start_mom << " GeV/c";
   }
   if (particle.Polarization().Mag2() != 0.) {
     out << " with polarization " << lar::dump::vector3D(particle.Polarization());
   }
   // does this particle seem to stop? (by decay, or because of extended traj.)
   if ((nPoints > 1) || (nDaughters > 0)) {
     out << "\n"
         << indent << ((nDaughters > 0) ? "ends" : "stops") << " by "
         << (particle.EndProcess().empty() ? "magics" : particle.EndProcess());
     if (nPoints > 1) {
       TLorentzVector const& stop = particle.EndPosition();
       TLorentzVector const& stop_mom = particle.EndMomentum();
       out << " at " << stop << " cm with momentum " << stop_mom << " GeV/c";
     }
     if (nDaughters > 0) {
       out << " into ";
       if (nDaughters == 1)
         out << "particle ID=" << particle.FirstDaughter();
       else {
         out << nDaughters << " particles from ID=" << particle.FirstDaughter()
             << " to ID=" << particle.LastDaughter();
       }
     } // if daughters
   }   // if it stops
   if (nPoints > 1) {
     simb::MCTrajectory const& traj = particle.Trajectory();
     out << "\n"
         << indent << "comes with a trajectory " << traj.TotalLength() << " cm long in " << nPoints
         << " points";
   } // if has points
 
 } // sim::dump::DumpMCParticle()

template<typename Stream >

void sim::dump::DumpMCParticle	(	Stream &&	out,
		simb::MCParticle const &	particle,
		std::string	indent = `""`
	)

Dumps the content of the specified particle in the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
particle	the particle to be dumped
indent	base indentation string (default: none)
firstIndent	indentation of first line (default: as `indent`)

The indent string is prepended to every line of output, except for the first one, where firstIndent is used.

The output starts on the current line, and the last line is NOT broken.

Definition at line 59 of file MCDumpers.h.

References DumpMCParticle(), and DumpMCParticleTrajectory().

     {
       DumpMCParticle(std::forward<Stream>(out), particle, indent, indent);
     }

template<typename Stream >

void sim::dump::DumpMCParticleTrajectory	(	Stream &&	out,
		simb::MCTrajectory const &	trajectory,
		unsigned int	pointsPerLine,
		std::string	indent
	)

Dumps the specified particle trajectory into the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
trajectory	the particle trajectory to be dumped
pointsPerLine	number of points dumped per line (default: all)
indent	base indentation string (default: none)

All points of the specified Monte Carlo trajectory are printed on screen, pointsPerLine on each line. The points are printed starting on a new line, and each line is applied the same indentation string (indent). As an exception, if pointsPerLine is not specified, all points are printed on the current line and indent is ignored.

The last line of the output is NOT broken.

Definition at line 288 of file MCDumpers.h.

Referenced by DumpMCParticle(), sim::DumpMCParticles::DumpMCParticle(), DumpMCParticleTrajectory(), and DumpMCTruth().

 {
   unsigned int page = 0;
   for (auto const& pair : trajectory) {
     if ((pointsPerLine > 0) && (page-- == 0)) {
       out << "\n" << indent << "  ";
       page = pointsPerLine - 1;
     }
     else
       out << " -- ";
 
     TLorentzVector const& pos = pair.first;
     out << pos;
   } // for trajectory points
 
 } // sim::dump::DumpMCParticleTrajectory()

template<typename Stream >

void sim::dump::DumpMCParticleTrajectory	(	Stream &&	out,
		simb::MCTrajectory const &	trajectory
	)

Dumps the specified particle trajectory into the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
trajectory	the particle trajectory to be dumped
pointsPerLine	number of points dumped per line (default: all)
indent	base indentation string (default: none)

All points of the specified Monte Carlo trajectory are printed on screen, pointsPerLine on each line. The points are printed starting on a new line, and each line is applied the same indentation string (indent). As an exception, if pointsPerLine is not specified, all points are printed on the current line and indent is ignored.

The last line of the output is NOT broken.

Definition at line 92 of file MCDumpers.h.

References DumpMCNeutrino(), and DumpMCParticleTrajectory().

     {
       DumpMCParticleTrajectory(std::forward<Stream>(out), trajectory, 0U, "");
     }

template<typename Stream >

void sim::dump::DumpMCTruth	(	Stream &&	out,
		simb::MCTruth const &	truth,
		unsigned int	pointsPerLine,
		std::string	indent,
		std::string	firstIndent
	)

Dumps the content of the specified MC truth in the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
truth	the truth information to be dumped
pointsPerLine	(optional) number of trajectory points dumped per line (default: `0`, no points shown)
indent	base indentation string (default: none)
firstIndent	string to be used for indentation of the first line (default: same as `indent`)

The indent string is prepended to every line of output, except for the first one, where firstIndent is used.

The argument pointsPerLine regulates the dump of trajectory points from all the particles in the record (except the ones stored in the neutrino object). Setting it to 0, or leaving it out, will suppress the dump of particle trajectories completely. There is no option to reduce the number of printed trajectory points: it's just all or none.

The output starts on the current line, and the last line is NOT broken.

Definition at line 338 of file MCDumpers.h.

References DumpMCNeutrino(), DumpMCParticle(), DumpMCParticleTrajectory(), for(), simb::MCGeneratorInfo::generator, simb::MCTruth::GeneratorInfo(), simb::MCTruth::GetNeutrino(), simb::MCTruth::GetParticle(), simb::MCTruth::NeutrinoSet(), simb::MCTruth::NParticles(), simb::MCParticle::NumberTrajectoryPoints(), simb::MCTruth::Origin(), simb::MCParticle::Trajectory(), and sim::TruthOriginName().

Referenced by sim::DumpMCTruth::analyze(), DumpMCNeutrino(), DumpMCTruth(), evgen::GENIEGen::produce(), and larg4::LArG4::produce().

 {
   // add a dumper of the generator info
   auto genInfo = truth.GeneratorInfo();
   out << "Generator ID: " << (int)genInfo.generator << "\n";
   out << "Generator Version: " << genInfo.generatorVersion << "\n";
   out << "Generator configuration string map.  Pairs separated by colons. "
       << "\n";
   for (const auto& gci : genInfo.generatorConfig) {
     out << "   " << gci.first << " : " << gci.second << "\n";
   }
 
   unsigned int const nParticles = truth.NParticles();
   out << firstIndent << nParticles << " particles from " << TruthOriginName(truth.Origin());
   if (truth.NeutrinoSet()) {
     out << '\n' << indent << "neutrino information: ";
     DumpMCNeutrino(std::forward<Stream>(out), truth.GetNeutrino(), indent + "  ", "");
   }
   for (unsigned int i = 0; i < nParticles; ++i) {
     out << '\n' << indent << "[#" << i << "] ";
     simb::MCParticle const& particle = truth.GetParticle(i);
     DumpMCParticle(std::forward<Stream>(out), particle, indent + "  ", "");
 
     const unsigned int nPoints = particle.NumberTrajectoryPoints();
     if ((nPoints > 0) && (pointsPerLine > 0)) {
       out << ":";
       DumpMCParticleTrajectory(
         std::forward<Stream>(out), particle.Trajectory(), pointsPerLine, indent + "    ");
     } // if has points
   }   // for all particles
 
 } // sim::dump::DumpMCTruth()

template<typename Stream >

void sim::dump::DumpMCTruth	(	Stream &&	out,
		simb::MCTruth const &	truth,
		unsigned int	pointsPerLine,
		std::string	indent = `""`
	)

Dumps the content of the specified MC truth in the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
truth	the truth information to be dumped
pointsPerLine	(optional) number of trajectory points dumped per line (default: `0`, no points shown)
indent	base indentation string (default: none)
firstIndent	string to be used for indentation of the first line (default: same as `indent`)

The indent string is prepended to every line of output, except for the first one, where firstIndent is used.

The argument pointsPerLine regulates the dump of trajectory points from all the particles in the record (except the ones stored in the neutrino object). Setting it to 0, or leaving it out, will suppress the dump of particle trajectories completely. There is no option to reduce the number of printed trajectory points: it's just all or none.

The output starts on the current line, and the last line is NOT broken.

Definition at line 161 of file MCDumpers.h.

References DumpMCTruth().

     {
       DumpMCTruth(std::forward<Stream>(out), truth, pointsPerLine, indent, indent);
     }

template<typename Stream >

void sim::dump::DumpMCTruth	(	Stream &&	out,
		simb::MCTruth const &	truth,
		std::string	indent,
		std::string	firstIndent
	)

Dumps the content of the specified MC truth in the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
truth	the truth information to be dumped
pointsPerLine	(optional) number of trajectory points dumped per line (default: `0`, no points shown)
indent	base indentation string (default: none)
firstIndent	string to be used for indentation of the first line (default: same as `indent`)

The indent string is prepended to every line of output, except for the first one, where firstIndent is used.

The argument pointsPerLine regulates the dump of trajectory points from all the particles in the record (except the ones stored in the neutrino object). Setting it to 0, or leaving it out, will suppress the dump of particle trajectories completely. There is no option to reduce the number of printed trajectory points: it's just all or none.

The output starts on the current line, and the last line is NOT broken.

Definition at line 170 of file MCDumpers.h.

References DumpMCTruth().

     {
       DumpMCTruth(std::forward<Stream>(out), truth, 0, indent, firstIndent);
     }

template<typename Stream >

void sim::dump::DumpMCTruth	(	Stream &&	out,
		simb::MCTruth const &	truth,
		std::string	indent = `""`
	)

Dumps the content of the specified MC truth in the output stream.

Template Parameters

Stream the type of output stream

Parameters

out	the output stream
truth	the truth information to be dumped
pointsPerLine	(optional) number of trajectory points dumped per line (default: `0`, no points shown)
indent	base indentation string (default: none)
firstIndent	string to be used for indentation of the first line (default: same as `indent`)

The indent string is prepended to every line of output, except for the first one, where firstIndent is used.

The argument pointsPerLine regulates the dump of trajectory points from all the particles in the record (except the ones stored in the neutrino object). Setting it to 0, or leaving it out, will suppress the dump of particle trajectories completely. There is no option to reduce the number of printed trajectory points: it's just all or none.

The output starts on the current line, and the last line is NOT broken.

Definition at line 179 of file MCDumpers.h.

References DumpGTruth(), and DumpMCTruth().

     {
       DumpMCTruth(std::forward<Stream>(out), truth, indent, indent);
     }