#include <quantumnumbers.h>

A container for storing conserved values.

QuantumNumbers can be used to compare all values that should be conserved during the evolution, in particular for the whole system at fixed intervals, or for a subset of particles before- and after a DecayAction / ScatterAction / Action.

QuantumNumbers can also be used to store, and retrieve, the total value of the quantum numbers in e.g. a ScatterAction, before distributing it to the new particles.

Currently, momentum conservation (including energy conservation), charge-, isospin3-, (net-) strangeness-, (net-) bottomness-, (net-) charmness- and (net-) baryon number conservation are checked. (It should be noted, or repeated, that only the net quantities are conserved, hence this is what is stored and compared with this class.)

Usage

QuantumNumbers before(particlelist);
do_something_with(particlelist);
QuantumNumbers after(particlelist);
printf("%s", before.report_deviations(after).c_str());
if (before != after) {
  throw std::runtime_error(before.report_deviations(after));
}

Definition at line 54 of file quantumnumbers.h.

Public Member Functions
	QuantumNumbers ()
	Construct QuantumNumbers collection with all fields 0. More...

	QuantumNumbers (const FourVector &m, const int q, const int i3, const int s, const int c, const int b, const int B)

	QuantumNumbers (const Particles &particles)
	Construct QuantumNumbers collection from the conserved quantities found in a set of particles. More...

	QuantumNumbers (const std::vector< Particles > &ensembles)
	Construct QuantumNumbers from several sets of particles. More...

	QuantumNumbers (const ParticleList &part)
	Construct QuantumNumbers collection from a particle list. More...

void	add_values (const ParticleData &p)
	Add the quantum numbers of a single particle to the collection. More...

FourVector	momentum () const

int	charge () const

int	isospin3 () const

int	strangeness () const

int	charmness () const

int	bottomness () const

int	baryon_number () const

bool	operator== (const QuantumNumbers &rhs) const

bool	operator!= (const QuantumNumbers &rhs) const
	Logical complement of QuantumNumbers::operator==. More...

QuantumNumbers	operator- (const QuantumNumbers &rhs) const

std::string	report_deviations (const std::vector< Particles > &ensembles) const
	Checks if the current particle list has still the same values and reports about differences. More...

std::string	report_deviations (const QuantumNumbers &rhs) const
	Reports on deviations between two QuantumNumbers collections. More...

Private Attributes
FourVector	momentum_
	Total momentum four-vector [GeV]. More...

int	charge_
	Total charge. More...

int	isospin3_
	Total isospin-3. More...

int	strangeness_
	Total strangeness. More...

int	charmness_
	Total charm. More...

int	bottomness_
	Total bottom. More...

int	baryon_number_
	Total baryon number. More...

Constructor & Destructor Documentation

◆ QuantumNumbers() [1/5]

smash::QuantumNumbers::QuantumNumbers ( )

inline

Construct QuantumNumbers collection with all fields 0.

Definition at line 57 of file quantumnumbers.h.

       : momentum_(0., 0., 0., 0.),
         charge_(0),
         isospin3_(0),
         strangeness_(0),
         charmness_(0),
         bottomness_(0),
         baryon_number_(0) {}

◆ QuantumNumbers() [2/5]

smash::QuantumNumbers::QuantumNumbers	(	const FourVector &	m,
		const int	q,
		const int	i3,
		const int	s,
		const int	c,
		const int	b,
		const int	B
	)

inline

Returns: Constructed QuantumNumbers.

Parameters

[in]	m	Momentum FourVector [GeV]
[in]	q	Charge
[in]	i3	Isospin
[in]	s	Strangeness
[in]	c	Charmness
[in]	b	Bottomness
[in]	B	Baryon number

Definition at line 76 of file quantumnumbers.h.

       : momentum_(m),
         charge_(q),
         isospin3_(i3),
         strangeness_(s),
         charmness_(c),
         bottomness_(b),
         baryon_number_(B) {}

◆ QuantumNumbers() [3/5]

smash::QuantumNumbers::QuantumNumbers ( const Particles & particles )

inlineexplicit

Construct QuantumNumbers collection from the conserved quantities found in a set of particles.

Parameters

[in] particles set of particles for which quantum numbers are calculated and constructed

Returns: Constructed object.

Definition at line 93 of file quantumnumbers.h.

                                                       : QuantumNumbers() {
     for (const ParticleData& data : particles) {
       add_values(data);
     }
   }

◆ QuantumNumbers() [4/5]

smash::QuantumNumbers::QuantumNumbers ( const std::vector< Particles > & ensembles )

inlineexplicit

Construct QuantumNumbers from several sets of particles.

Parameters

[in] ensembles Current list of particles in all ensembles

Returns: Constructed object.

Definition at line 104 of file quantumnumbers.h.

       : QuantumNumbers() {
     for (const Particles& particles : ensembles) {
       for (const ParticleData& data : particles) {
         add_values(data);
       }
     }
   }

◆ QuantumNumbers() [5/5]

smash::QuantumNumbers::QuantumNumbers ( const ParticleList & part )

inlineexplicit

Construct QuantumNumbers collection from a particle list.

Parameters

[in] part list of particles for which quantum numbers are calculated and constructed

Returns: Constructed object.

Definition at line 119 of file quantumnumbers.h.

                                                     : QuantumNumbers() {
     for (const auto& p : part) {
       add_values(p);
     }
   }

Member Function Documentation

◆ add_values()

void smash::QuantumNumbers::add_values ( const ParticleData & p )

inline

Add the quantum numbers of a single particle to the collection.

Parameters

[in] p particle whose quantum number is added to the collection

Definition at line 129 of file quantumnumbers.h.

                                          {
     momentum_ += p.momentum();
     charge_ += p.pdgcode().charge();
     isospin3_ += p.pdgcode().isospin3();
     strangeness_ += p.pdgcode().strangeness();
     charmness_ += p.pdgcode().charmness();
     bottomness_ += p.pdgcode().bottomness();
     baryon_number_ += p.pdgcode().baryon_number();
   }

◆ momentum()

FourVector smash::QuantumNumbers::momentum ( ) const

inline

Returns: The total momentum four-vector. \(P^\mu = \sum_{i \in \mbox{particles}} (E_i, \mathbf{p}_i)\) [GeV]

See also: QuantumNumbers::momentum_; ParticleData::momentum()

Definition at line 146 of file quantumnumbers.h.

146 { return momentum_; }

◆ charge()

int smash::QuantumNumbers::charge ( ) const

inline

Returns: The total electric charge. \(Q = \sum_{i \in \mbox{particles}} q_i\)

See also: QuantumNumbers::charge_; PdgCode::charge()

Definition at line 155 of file quantumnumbers.h.

155 { return charge_; }

◆ isospin3()

int smash::QuantumNumbers::isospin3 ( ) const

inline

Returns: Twice the total isospin-3 component. \(I = \sum_{i \in \mbox{particles}} 2{I_3}_i\)

See also: QuantumNumbers::isospin3_; PdgCode::isospin3()

Definition at line 163 of file quantumnumbers.h.

163 { return isospin3_; }

◆ strangeness()

int smash::QuantumNumbers::strangeness ( ) const

inline

Returns: The total strangeness. \(S = \sum_{i \in \mbox{particles}} S_i\)

See also: QuantumNumbers::strangeness_; PdgCode::strangeness()

Definition at line 171 of file quantumnumbers.h.

171 { return strangeness_; }

◆ charmness()

int smash::QuantumNumbers::charmness ( ) const

inline

Returns: The total charm. \(C = \sum_{i \in \mbox{particles}} C_i\)

See also: QuantumNumbers::charmness_; PdgCode::charmness()

Definition at line 179 of file quantumnumbers.h.

179 { return charmness_; }

◆ bottomness()

int smash::QuantumNumbers::bottomness ( ) const

inline

Returns: The total bottom. \(b = \sum_{i \in \mbox{particles}} b_i\)

See also: QuantumNumbers::bottomness_; PdgCode::bottomness()

Definition at line 187 of file quantumnumbers.h.

187 { return bottomness_; }

◆ baryon_number()

int smash::QuantumNumbers::baryon_number ( ) const

inline

Returns: The total baryon number. \(B = \sum_{i \in \mbox{particles}} B_i\)

See also: QuantumNumbers::baryon_number_; PdgCode::baryon_number()

Definition at line 195 of file quantumnumbers.h.

195 { return baryon_number_; }

◆ operator==()

bool smash::QuantumNumbers::operator== ( const QuantumNumbers & rhs ) const

inline

Returns: true if all members compare true.

Parameters

rhs	Right-hand side.

In the comparison of FourVectors, a little leeway is built-in, so this does not rely on an exact comparison of floating point values.

See also: FourVector::operator==

Definition at line 206 of file quantumnumbers.h.

                                                    {
     // clang-format off
     // invariant mass of the current momentum
     const double m0_invariant = std::sqrt(std::fabs(momentum_.sqr()));
     // invariant mass of the momentum to be compared
     const double m1_invariant = std::sqrt(std::fabs(rhs.momentum_.sqr()));
     // average invariant mass
     const double avg_m_inv = 0.5 * (m0_invariant + m1_invariant);
     // difference in the four-momentum
     FourVector momentum_diff = momentum_ - rhs.momentum_;
     /* Check whether momentum is conserved.
      * The momentum is considered to be conserved when the deviation
      * is small enough compared to the given invariant mass. */
     bool momentum_conserved =
         (std::fabs(momentum_diff.x0()) <= small_number * avg_m_inv &&
          std::fabs(momentum_diff.x1()) <= small_number * avg_m_inv &&
          std::fabs(momentum_diff.x2()) <= small_number * avg_m_inv &&
          std::fabs(momentum_diff.x3()) <= small_number * avg_m_inv);
     return (momentum_conserved &&
             charge_ == rhs.charge_ &&
             isospin3_ == rhs.isospin3_ &&
             strangeness_ == rhs.strangeness_ &&
             charmness_ == rhs.charmness_ &&
             bottomness_ == rhs.bottomness_ &&
             baryon_number_ == rhs.baryon_number_);
     // clang-format on
   }

◆ operator!=()

bool smash::QuantumNumbers::operator!= ( const QuantumNumbers & rhs ) const

inline

Logical complement of QuantumNumbers::operator==.

Definition at line 234 of file quantumnumbers.h.

234 { return !(*this == rhs); }

◆ operator-()

QuantumNumbers smash::QuantumNumbers::operator- ( const QuantumNumbers & rhs ) const

inline

Returns: Entry-wise difference of two sets of QuantumNumbers.

Parameters

rhs	Right-hand side.

If everything is conserved, all entries of the result should be zero.

Definition at line 243 of file quantumnumbers.h.

                                                             {
     return {momentum_ - rhs.momentum_,          charge_ - rhs.charge_,
             isospin3_ - rhs.isospin3_,          strangeness_ - rhs.strangeness_,
             charmness_ - rhs.charmness_,        bottomness_ - rhs.bottomness_,
             baryon_number_ - rhs.baryon_number_};
   }

◆ report_deviations() [1/2]

std::string smash::QuantumNumbers::report_deviations ( const std::vector< Particles > & ensembles ) const

inline

Checks if the current particle list has still the same values and reports about differences.

Parameters

[in] ensembles Current list of particles in all ensembles, whose quantum number is compared

Returns: String reporting the deviations.

See also: QuantumNumbers::report_deviations(const QuantumNumbers&) const

Definition at line 259 of file quantumnumbers.h.

                                                                            {
     QuantumNumbers current_values(ensembles);
     return report_deviations(current_values);
   }

◆ report_deviations() [2/2]

std::string smash::QuantumNumbers::report_deviations ( const QuantumNumbers & rhs ) const

Reports on deviations between two QuantumNumbers collections.

Parameters

[in] rhs Other QuantumNumbers collection.

Returns: A string with information about the differences.

If there are no differences, the returned string is empty; else, a descriptive warning message is returned, e.g.

Conservation law violations detected (old vs. new)
Deviation in Charge:
 164 vs. 163
Deviation in Isospin 3:
 -88 vs. -90

Definition at line 18 of file quantumnumbers.cc.

                                                                            {
   if (rhs == *this) {
     return "";
   }
   std::stringstream error_msg;
   if (momentum_ != rhs.momentum_) {
     error_msg << "Deviation in Four-Momentum:\n" << std::scientific;
   }
   /* programmer's note: here, I'd like to simultaneously loop over an
    * integer (for the output; so that we know which component is
    * faulty) and both the current and rhs's momentum four-vector. If
    * there is a better way to do this, feel free to implement.
    *
    * I chose mu < 4 as the breaking condition out of the vague feeling
    * that comparing integers may be faster than accessing the
    * iterators. */
   int mu = 0;
   for (auto here_iter = momentum_.cbegin(), rhs_iter = rhs.momentum_.cbegin();
        mu < 4; ++here_iter, ++rhs_iter, ++mu) {
     if (!almost_equal_physics(*here_iter, *rhs_iter)) {
       error_msg << " P_" << mu << ": " << *here_iter << " vs. " << *rhs_iter
                 << "; Δ = " << (*here_iter - *rhs_iter) << "\n";
     }
   }
   if (charge_ != rhs.charge_) {
     error_msg << "Deviation in Charge:\n " << charge_ << " vs. " << rhs.charge_
               << "\n";
   }
   if (isospin3_ != rhs.isospin3_) {
     error_msg << "Deviation in Isospin 3:\n " << isospin3_ << " vs. "
               << rhs.isospin3_ << "\n";
   }
   if (strangeness_ != rhs.strangeness_) {
     error_msg << "Deviation in Strangeness:\n " << strangeness_ << " vs. "
               << rhs.strangeness_ << "\n";
   }
   if (charmness_ != rhs.charmness_) {
     error_msg << "Deviation in Charmness:\n " << charmness_ << " vs. "
               << rhs.charmness_ << "\n";
   }
   if (bottomness_ != rhs.bottomness_) {
     error_msg << "Deviation in Bottomness:\n " << bottomness_ << " vs. "
               << rhs.bottomness_ << "\n";
   }
   if (baryon_number_ != rhs.baryon_number_) {
     error_msg << "Deviation in Baryon Number:\n " << baryon_number_ << " vs. "
               << rhs.baryon_number_ << "\n";
   }
   return error_msg.str();
 }

Member Data Documentation

◆ momentum_

FourVector smash::QuantumNumbers::momentum_

private

Total momentum four-vector [GeV].

See also: QuantumNumbers::momentum()

Definition at line 290 of file quantumnumbers.h.

◆ charge_

int smash::QuantumNumbers::charge_

private

Total charge.

See also: QuantumNumbers::charge()

Definition at line 297 of file quantumnumbers.h.

◆ isospin3_

int smash::QuantumNumbers::isospin3_

private

Total isospin-3.

See also: QuantumNumbers::isospin3()

Definition at line 303 of file quantumnumbers.h.

◆ strangeness_

int smash::QuantumNumbers::strangeness_

private

Total strangeness.

See also: QuantumNumbers::strangeness()

Definition at line 309 of file quantumnumbers.h.

◆ charmness_

int smash::QuantumNumbers::charmness_

private

Total charm.

See also: QuantumNumbers::charmness()

Definition at line 315 of file quantumnumbers.h.

◆ bottomness_

int smash::QuantumNumbers::bottomness_

private

Total bottom.

See also: QuantumNumbers::bottomness()

Definition at line 321 of file quantumnumbers.h.

◆ baryon_number_

int smash::QuantumNumbers::baryon_number_

private

Total baryon number.

See also: QuantumNumbers::baryon_number()

Definition at line 327 of file quantumnumbers.h.

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

src/include/smash/quantumnumbers.h
src/quantumnumbers.cc

Usage

Public Member Functions

Private Attributes

Constructor & Destructor Documentation

◆ QuantumNumbers() [1/5]

◆ QuantumNumbers() [2/5]

◆ QuantumNumbers() [3/5]

◆ QuantumNumbers() [4/5]

◆ QuantumNumbers() [5/5]

Member Function Documentation

◆ add_values()

◆ momentum()

◆ charge()

◆ isospin3()

◆ strangeness()

◆ charmness()

◆ bottomness()

◆ baryon_number()

◆ operator==()

◆ operator!=()

◆ operator-()

◆ report_deviations() [1/2]

◆ report_deviations() [2/2]

Member Data Documentation

◆ momentum_

◆ charge_

◆ isospin3_

◆ strangeness_

◆ charmness_

◆ bottomness_

◆ baryon_number_