Functions
size_t	index (const Gaudi::ParticleProperty property, const Gaudi::Interfaces::IParticlePropertySvc service)
	helper utility for mapping of Gaudi::ParticleProperty object into non-negative integral sequential identifier

size_t	index (const Gaudi::ParticleID &pid, const Gaudi::Interfaces::IParticlePropertySvc *service)
	helper utility for mapping of Gaudi::ParticleID object into non-negative integral sequential identifier

const Gaudi::ParticleProperty *	particle (const size_t index, const Gaudi::Interfaces::IParticlePropertySvc *service)
	the inverse mapping of the integer sequential number onto Gaudi::ParticleID object

const Gaudi::ParticleID	particleID (const size_t index, const Gaudi::Interfaces::IParticlePropertySvc *service)
	the inverse mapping of the integer sequential number onto Gaudi::ParticleID object

const Gaudi::ParticleProperty *	byPythiaID (const int pythia, const Gaudi::Interfaces::IParticlePropertySvc *svc)
	mapping by pythiaID

const Gaudi::ParticleProperty *	byEvtGenName (const std::string &evtGen, const Gaudi::Interfaces::IParticlePropertySvc *svc)
	mapping by EvtGen-name

template<class INPUT, class PREDICATE, class OUTPUT>
OUTPUT	get (INPUT first, INPUT last, const PREDICATE &cut, OUTPUT output)
	get all particle properties which satisfy the certain criteria

template<class PREDICATE, class OUTPUT>
OUTPUT	get (const Gaudi::Interfaces::IParticlePropertySvc *service, const PREDICATE &cut, OUTPUT output)
	get all particle properties which satisfy the certain criteria

Gaudi::Interfaces::IParticlePropertySvc::ParticleProperties	allProperties (const Gaudi::Interfaces::IParticlePropertySvc *service)
	get all the properties at once

GAUDI_API std::ostream &	printAsTable_ (const std::vector< const Gaudi::ParticleProperty * > &particles, std::ostream &stream, const Gaudi::Interfaces::IParticlePropertySvc *service=0)
	print a list of properties in a form of the table

GAUDI_API std::string	printAsTable (const std::vector< const Gaudi::ParticleProperty * > &particles, const Gaudi::Interfaces::IParticlePropertySvc *service=0)
	print a list of properties in a form of the table

GAUDI_API MsgStream &	printAsTable (const std::vector< const Gaudi::ParticleProperty * > &particles, MsgStream &stream, const Gaudi::Interfaces::IParticlePropertySvc *service=0)
	print a list of properties in a form of the table

GAUDI_API std::ostream &	printAsTable_ (const std::vector< Gaudi::ParticleID > &particles, std::ostream &stream, const Gaudi::Interfaces::IParticlePropertySvc *service=0)
	print properties in a form of the table

template<class C_, class A_>
std::ostream &	printAsTable_ (const std::set< Gaudi::ParticleID, C_, A_ > &particles, std::ostream &stream, const Gaudi::Interfaces::IParticlePropertySvc *service=0)

GAUDI_API std::string	printAsTable (const std::vector< Gaudi::ParticleID > &particles, const Gaudi::Interfaces::IParticlePropertySvc *service=0)
	print properties in a form of the table

template<class C_, class A_>
std::string	printAsTable (const std::set< Gaudi::ParticleID, C_, A_ > &particles, const Gaudi::Interfaces::IParticlePropertySvc *service=0)

GAUDI_API MsgStream &	printAsTable (const std::vector< Gaudi::ParticleID > &particles, MsgStream &stream, const Gaudi::Interfaces::IParticlePropertySvc *service=0)
	print properties in a form of the table

template<class C_, class A_>
MsgStream &	printAsTable (const std::set< Gaudi::ParticleID, C_, A_ > &particles, MsgStream &stream, const Gaudi::Interfaces::IParticlePropertySvc *service=0)

Function Documentation

◆ allProperties()

Gaudi::Interfaces::IParticlePropertySvc::ParticleProperties Gaudi::ParticleProperties::allProperties ( const Gaudi::Interfaces::IParticlePropertySvc * service )

get all the properties at once

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 201 of file IParticlePropertySvc.cpp.

                                                                                           {
  return service ? Gaudi::Interfaces::IParticlePropertySvc::ParticleProperties( service->begin(), service->end() )
                 : Gaudi::Interfaces::IParticlePropertySvc::ParticleProperties();
}

◆ byEvtGenName()

const Gaudi::ParticleProperty * Gaudi::ParticleProperties::byEvtGenName	(	const std::string &	evtGen,
		const Gaudi::Interfaces::IParticlePropertySvc *	svc )

mapping by EvtGen-name

const std::string& evtGen
 
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
const Gaudi::ParticleProperty* pp = byEvtGenName ( pythiaID , svc ) ;

Attention: the method is not very efficient and should not be abused

See also: Gaudi::ParticleProperties::particle

Parameters

evtGen	the particle naem in EvtGen-generator
svc	pointer to particle property service

Returns: the particle property for the given EvtGen-name

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 180 of file IParticlePropertySvc.cpp.

                                                                                          {
  if ( !svc ) { return nullptr; }
  // to be more efficient:
  // 1) try to use the regular name (fast, logarithmic search)
  const Gaudi::ParticleProperty* pp = svc->find( evtGen );
  // 2) check the proper evtgen name
  if ( pp && evtGen == pp->evtGen() ) { return pp; }
  // 3) use the regular (linear search)
  auto begin = svc->begin();
  auto end   = svc->end();
  auto found = std::find_if( begin, end, [&]( const Gaudi::ParticleProperty* pp ) { return pp->evtGen() == evtGen; } );
  //
  return found != end ? *found : nullptr;
}

◆ byPythiaID()

const Gaudi::ParticleProperty * Gaudi::ParticleProperties::byPythiaID	(	const int	pythia,
		const Gaudi::Interfaces::IParticlePropertySvc *	svc )

mapping by pythiaID

const int pythiaID = ... ;
 
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
const Gaudi::ParticleProeprty* pp = byPythiaID( pythiaID , svc ) ;

Attention: the method is not very efficient and should not be abused

See also: Gaudi::ParticleProperties::particle

Parameters

pythia	pythia identifier
svc	pointer to particle property service

Returns: the particle property for the given pythia ID

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 144 of file IParticlePropertySvc.cpp.

                                                                                                      {
  if ( !svc ) { return nullptr; }
  // to be efficient
  // 1) try to use PDG-ID (fast, logarithmic search)
  const Gaudi::ParticleProperty* pp = svc->find( Gaudi::ParticleID( pythia ) );
  // 2) check the proper pythia ID
  if ( pp && pythia == pp->pythiaID() ) { return pp; }
  // 3) use the regular (linear search)
  auto begin = svc->begin();
  auto end   = svc->end();
  auto found =
      std::find_if( begin, end, [&]( const Gaudi::ParticleProperty* pp ) { return pp->pythiaID() == pythia; } );
  return found != end ? *found : nullptr;
}

◆ get() [1/2]

template<class PREDICATE, class OUTPUT>

OUTPUT Gaudi::ParticleProperties::get	(	const Gaudi::Interfaces::IParticlePropertySvc *	service,
		const PREDICATE &	cut,
		OUTPUT	output )

get all particle properties which satisfy the certain criteria

e.g. select all particles with 'invalid' Pythia ID = 0 :

#include "boost/lambda/lambda.hpp"
#include "boost/lambda/bind.hpp"
 
typedef Gaudi::Interfaces::IParticlePropertySvc::ParticleProperties Vector ;
using namespace boost::lambda ;
 
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
// create the output vector:
Vector output ;
// use the function
Gaudi::ParticleProperties::get
  ( svc ,
    bind ( &Gaudi::ParticleProperty::pythiaID , _1 ) == 0 , // predicate
    std::back_inserter( output ) ) ; // output

or select all particles with 'invalid' EvtGen = "unknown" :

#include "boost/lambda/lambda.hpp"
#include "boost/lambda/bind.hpp"
 
typedef Gaudi::Interfaces::IParticlePropertySvc::ParticleProperties Vector ;
using namespace boost::lambda ;
 
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
// create the output vector:
Vector output ;
// use the function
Gaudi::ParticleProperties::get
  ( svc ,
    bind ( &Gaudi::ParticleProperty::evtGen , _1 ) == "unknown" , // predicate
    std::back_inserter( output ) ) ; // output

Select all leptons:

#include "boost/lambda/lambda.hpp"
#include "boost/lambda/bind.hpp"
 
typedef Gaudi::Interfaces::IParticlePropertySvc::ParticleProperties Vector ;
using namespace boost::lambda ;
 
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
// create the output vector:
Vector leptons ;
// use the function
Gaudi::ParticleProperties::get
  ( svc ,
    bind ( &Gaudi::ParticleID::isLepton ,
    bind ( &Gaudi::ParticleProperty::particleID , _1 ) ) ,
    std::back_inserter ( leptons ) ) ; // output

Parameters

first	begin-iterator of input sequence of particle properties
end	end-iterator of input sequence of particle properties
cut	the predicate
output	the output iterator

Returns: the updated position of output iterator

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 476 of file IParticlePropertySvc.h.

                                                                                                          {
      if ( 0 == service ) { return output; }
      return service->get( cut, output );
    }

◆ get() [2/2]

template<class INPUT, class PREDICATE, class OUTPUT>

OUTPUT Gaudi::ParticleProperties::get	(	INPUT	first,
		INPUT	last,
		const PREDICATE &	cut,
		OUTPUT	output )

get all particle properties which satisfy the certain criteria

e.g. select all particles with 'invalid' Pythia ID = 0 :

#include "boost/lambda/lambda.hpp"
#include "boost/lambda/bind.hpp"
 
typedef Gaudi::Interfaces::IParticlePropertySvc::ParticleProperties Vector ;
using namespace boost::lambda ;
 
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
// create the output vector:
Vector output ;
// use the function
Gaudi::ParticleProperties::get
       ( svc->begin() ,  // begin-sequence
         svc->end  () ,  // end-sequence
         bind ( &Gaudi::ParticleProperty::pythiaID , _1 ) == 0 , // predicate
         std::back_inserter( output ) ) ; // output

or select all particles with 'invalid' EvtGen = "unknown" :

#include "boost/lambda/lambda.hpp"
#include "boost/lambda/bind.hpp"
 
typedef Gaudi::Interfaces::IParticlePropertySvc::ParticleProperties Vector ;
using namespace boost::lambda ;
 
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
// create the output vector:
Vector output ;
// use the function
Gaudi::ParticleProperties::get
       ( svc->begin() ,  // begin-sequence
         svc->end  () ,  // end-sequence
         bind ( &Gaudi::ParticleProperty::evtGen , _1 ) == "unknown" , // predicate
         std::back_inserter( output ) ) ; // output

select leptons:

#include "boost/lambda/lambda.hpp"
#include "boost/lambda/bind.hpp"
...
...
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
typedef Gaudi::Interfaces::IParticlePropertySvc::ParticleProperties Vector ;
using namespace boost::lambda ;
 
// create the output vector:
Vector leptons ;
// use the function
Gaudi::ParticleProperties::get
       ( svc->begin() ,  // begin-sequence
         svc->end  () ,  // end-sequence
         bind ( &Gaudi::ParticleID::isLepton , bind ( &Gaudi::ParticleProperty::particleID ,_1 ) ) ,
         std::back_inserter ( lepton ) ) ; // output

Essentially it is just a missing std::copy_if STL-algorithm

Parameters

first	begin-iterator of input sequence of particle properties
last	end-iterator of input sequence of particle properties
cut	the predicate
output	the output iterator

Returns: the updated position of output iterator

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 388 of file IParticlePropertySvc.h.

                                                                               {
      for ( ; first != last; ++first ) {
        if ( cut( *first ) ) {
          *output = *first;
          ++output;
        }
      }
      return output;
    }

◆ index() [1/2]

size_t Gaudi::ParticleProperties::index	(	const Gaudi::ParticleID &	pid,
		const Gaudi::Interfaces::IParticlePropertySvc *	service )

helper utility for mapping of Gaudi::ParticleID object into non-negative integral sequential identifier

const Gaudi::ParticlePropertySvc* svc = ... ;
 
const Gaudi::ParticleID& pid = ... ;
 
size_t index = index ( pid , svc ) ;
 
if ( !index )
 {
   // error here: conversion is not possible
 }

This appears to be useful operation, but since it is "pure technical" it does not appear as interface method.

For invalid/missing PID and/or service 0 is returned. The valid result is always satisfy the condition: index <= service->size()

Parameters

pid	the object to be mapped
service	the service

Returns: the sequential non-negative index

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 67 of file IParticlePropertySvc.cpp.

                                                                                              {
  if ( !service ) { return 0; }
  const Gaudi::ParticleProperty* pp = service->find( pid );
  return pp ? Gaudi::ParticleProperties::index( pp, service ) : 0;
}

◆ index() [2/2]

size_t Gaudi::ParticleProperties::index	(	const Gaudi::ParticleProperty *	property,
		const Gaudi::Interfaces::IParticlePropertySvc *	service )

helper utility for mapping of Gaudi::ParticleProperty object into non-negative integral sequential identifier

const Gaudi::ParticlePropertySvc* svc = ... ;
 
const Gaudi::ParticleProperty* pp = ... ;
 
 size_t index = index ( pp , svc ) ;
 
 if ( !index )
 {
   // error here: conversion is not possible
 }

This appears to be useful operation, but since it is "pure technical" it does not appear as interface method.

For invalid/missing property and/or service 0 is returned. The valid result is always satisfy the condition: index <= service->size()

Parameters

property	the property to be mapped
service	the service

Returns: the sequential non-negative index

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 39 of file IParticlePropertySvc.cpp.

                                                                                              {
  if ( !property || !service ) { return 0; }
 
  auto first = service->begin();
  auto last  = service->end();
  // start the binary_search
  static constexpr auto cmp   = Gaudi::ParticleProperty::Compare();
  auto                  ifind = std::lower_bound( first, last, property, cmp );
  return last != ifind && !cmp( *ifind, property ) ? ( ifind - first + 1 ) : 0;
}

◆ particle()

const Gaudi::ParticleProperty * Gaudi::ParticleProperties::particle	(	const size_t	index,
		const Gaudi::Interfaces::IParticlePropertySvc *	service )

the inverse mapping of the integer sequential number onto Gaudi::ParticleID object

const Gaudi::ParticlePropertySvc* svc = ... ;
 
const size_t index = ...
 
const Gaudi::ParticleProperty* pp = particle ( index , svc ) ;
 
if ( !pp )
 {
   // error here: conversion is not possible
 }

This appears to be useful operation, but since it is "pure technical" it does not appear as interface method.

For invalid/missing PID and/or service NULL is returned.

Parameters

pid	the object to be mapped
service	the service

Returns: the sequential non-negative index

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 90 of file IParticlePropertySvc.cpp.

                                                                                                          {
  if ( 0 == index || !service ) { return nullptr; }
  // get the iterators from the service
  auto first = service->begin();
  auto last  = service->end();
  if ( index > (size_t)std::distance( first, last ) ) { return 0; }
  std::advance( first, index - 1 );
  return *first;
}

◆ particleID()

const Gaudi::ParticleID Gaudi::ParticleProperties::particleID	(	const size_t	index,
		const Gaudi::Interfaces::IParticlePropertySvc *	service )

the inverse mapping of the integer sequential number onto Gaudi::ParticleID object

const Gaudi::ParticlePropertySvc* svc = ... ;
 
const size_t index = ...
 
const Gaudi::ParticleID pid = particleID ( index , svc ) ;
 
if ( Gaudi::ParticleID() == pid )
 {
   // error here: conversion is not possible
 }

This appears to be useful operation, but since it is "pure technical" it does not appear as interface method.

For invalid/missing index and/or service Gaudi::ParticleID() is returned.

Parameters

pid	the object to be mapped
service	the service

Returns: the sequential non-negative index

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 117 of file IParticlePropertySvc.cpp.

                                                                                                            {
  if ( 0 == index || !service ) { return Gaudi::ParticleID(); }
  const Gaudi::ParticleProperty* pp = Gaudi::ParticleProperties::particle( index, service );
  return pp ? pp->particleID() : Gaudi::ParticleID();
}

◆ printAsTable() [1/6]

template<class C_, class A_>

std::string Gaudi::ParticleProperties::printAsTable	(	const std::set< Gaudi::ParticleID, C_, A_ > &	particles,
		const Gaudi::Interfaces::IParticlePropertySvc *	service = 0 )

inline

Definition at line 376 of file ParticleProperty.h.

                                                                                              {
      return printAsTable( std::vector<Gaudi::ParticleID>( particles.begin(), particles.end() ), service );
    }

◆ printAsTable() [2/6]

template<class C_, class A_>

MsgStream & Gaudi::ParticleProperties::printAsTable	(	const std::set< Gaudi::ParticleID, C_, A_ > &	particles,
		MsgStream &	stream,
		const Gaudi::Interfaces::IParticlePropertySvc *	service = 0 )

inline

Definition at line 394 of file ParticleProperty.h.

                                                                                             {
      return printAsTable( std::vector<Gaudi::ParticleID>( particles.begin(), particles.end() ), stream, service );
    }

◆ printAsTable() [3/6]

std::string Gaudi::ParticleProperties::printAsTable	(	const std::vector< const Gaudi::ParticleProperty * > &	particles,
		const Gaudi::Interfaces::IParticlePropertySvc *	service = 0 )

print a list of properties in a form of the table

Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
const std::vector<const Gaudi::ParticleProperty*>& props = ... ;
 
std::cout << "Properties" << std::endl <<
Gaudi::ParticleProperties::printAsTable ( props , svc ) << std::endl

The utility is easy to use in conjunction with Gaudi::ParticleProperties:;get utilities: e.g. get all leptons from the service and print them as table:

#include "boost/lambda/lambda.hpp"
#include "boost/lambda/bind.hpp"
 
typedef std::vector<const Gaudi::IParticleProperty*> Vector ;
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
// create the output vector:
Vector leptons ;
// use the function
Gaudi::ParticleProperties::get
       ( svc ,
         // create the predicate:
         boost::lambda::bind ( &Gaudi::ParticleID::isLepton ,
             boost::lambda::bind ( &Gaudi::ParticleProperty::particleID , boost::lambda::_1 ) ) ,
         std::back_inserter ( leptons ) ) ; // output
 
// print the leptons:
std::cout << "LEPTONS"
          << std::endl
Gaudi::ParticleProperties::printAsTable ( leptons , svc )
          << std::endl ;

See also: Gaudi::ParticleProperty; Gaudi::Interfaces::IParticlePropertySvc; Gaudi::ParticleProperties::get

Parameters

particles	the list of particle properties
service	the service to extract global information

Returns: the string output

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 319 of file ParticleProperty.cpp.

                                                                                                              {
  std::ostringstream s;
  printAsTable_( particles, s, service );
  return s.str();
}

◆ printAsTable() [4/6]

MsgStream & Gaudi::ParticleProperties::printAsTable	(	const std::vector< const Gaudi::ParticleProperty * > &	particles,
		MsgStream &	stream,
		const Gaudi::Interfaces::IParticlePropertySvc *	service = 0 )

print a list of properties in a form of the table

Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
const std::vector<const Gaudi::ParticleProperty*>& props = ... ;
 
MsgStream& log = ... ;
Gaudi::ParticleProperties::printAsTable ( props , log , svc ) ;
log << endmsg ;

The utility is easy to use in conjunction with Gaudi::ParticleProperties:;get utilities: e.g. get all leptons from the service and print them as table:

#include "boost/lambda/lambda.hpp"
#include "boost/lambda/bind.hpp"
 
typedef std::vector<const Gaudi::IParticleProperty*> Vector ;
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
// create the output vector:
Vector leptons ;
// use the function
Gaudi::ParticleProperties::get
       ( svc ,
         // create the predicate:
         boost::lambda::bind ( &Gaudi::ParticleID::isLepton ,
             boost::lambda::bind ( &Gaudi::ParticleProperty::particleID , boost::lambda::_1 ) ) ,
         std::back_inserter ( leptons ) ) ; // output
 
// print the leptons:
MsgStream& log = ... ;
log << "LEPTONS" << std::endl ;
Gaudi::ParticleProperties::printAsTable ( leptons , log , svc ) ;
log << endmsg ;

See also: Gaudi::ParticleProperty; Gaudi::Interfaces::IParticlePropertySvc; Gaudi::ParticleProperties::get

Parameters

particles	the list of particle properties
stream	the reference to the output stream
service	the service to extract global information

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 333 of file ParticleProperty.cpp.

                                                                                                             {
  if ( stream.isActive() ) { printAsTable_( particles, stream.stream(), service ); }
  return stream;
}

◆ printAsTable() [5/6]

std::string Gaudi::ParticleProperties::printAsTable	(	const std::vector< Gaudi::ParticleID > &	particles,
		const Gaudi::Interfaces::IParticlePropertySvc *	service = 0 )

print properties in a form of the table

Parameters

particles	(INPUT) list of particles
service	(INPUT) pointer to particle property service

Returns: string-representation

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@i.nosp@m.tep..nosp@m.ru

Date: 2010-01-04

Definition at line 393 of file ParticleProperty.cpp.

                                                                                                          {
  std::ostringstream s;
  printAsTable_( particles, s, service );
  return s.str();
}

◆ printAsTable() [6/6]

MsgStream & Gaudi::ParticleProperties::printAsTable	(	const std::vector< Gaudi::ParticleID > &	particles,
		MsgStream &	stream,
		const Gaudi::Interfaces::IParticlePropertySvc *	service = 0 )

print properties in a form of the table

Parameters

particles	(INPUT) list of particles
stream	(UPDATE) the stream
service	(INPUT) pointer to particle property service

Returns: the stream

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@i.nosp@m.tep..nosp@m.ru

Date: 2010-01-04

Definition at line 381 of file ParticleProperty.cpp.

                                                                                                         {
  if ( stream.isActive() ) { printAsTable_( particles, stream.stream(), service ); }
  return stream;
}

◆ printAsTable_() [1/3]

template<class C_, class A_>

std::ostream & Gaudi::ParticleProperties::printAsTable_	(	const std::set< Gaudi::ParticleID, C_, A_ > &	particles,
		std::ostream &	stream,
		const Gaudi::Interfaces::IParticlePropertySvc *	service = 0 )

inline

Definition at line 359 of file ParticleProperty.h.

                                                                                                 {
      return printAsTable_( std::vector<Gaudi::ParticleID>( particles.begin(), particles.end() ), stream, service );
    }

◆ printAsTable_() [2/3]

std::ostream & Gaudi::ParticleProperties::printAsTable_	(	const std::vector< const Gaudi::ParticleProperty * > &	particles,
		std::ostream &	stream,
		const Gaudi::Interfaces::IParticlePropertySvc *	service = 0 )

print a list of properties in a form of the table

Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
const std::vector<const Gaudi::ParticleProperty*>& props = ... ;
 
std::cout << "Properties" << std::endl ;
Gaudi::ParticleProperties::printAsTable ( props , std::cout , svc ) ;

The utility is easy to use in conjunction with Gaudi::ParticleProperties:;get utilities: e.g. get all leptons from the service and print them as table:

#include "boost/lambda/lambda.hpp"
#include "boost/lambda/bind.hpp"
 
typedef std::vector<const Gaudi::IParticleProperty*> Vector ;
const Gaudi::Interfaces::IParticlePropertySvc* svc = ... ;
 
// create the output vector:
Vector leptons ;
// use the function
Gaudi::ParticleProperties::get
       ( svc ,
         // create the predicate:
         boost::lambda::bind ( &Gaudi::ParticleID::isLepton ,
             boost::lambda::bind ( &Gaudi::ParticleProperty::particleID , boost::lambda::_1 ) ) ,
         std::back_inserter ( leptons ) ) ; // output
 
// print the leptons:
std::cout << "LEPTONS" << std::endl ;
Gaudi::ParticleProperties::printAsTable_ ( leptons , std::cout , svc ) ;

See also: Gaudi::ParticleProperty; Gaudi::Interfaces::IParticlePropertySvc; Gaudi::ParticleProperties::get

Parameters

particles	the list of particle properties
stream	the reference to the output stream
service	the service to extract global information

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@n.nosp@m.ikhe.nosp@m.f.nl

Date: 2008-08-03

Definition at line 173 of file ParticleProperty.cpp.

                                                                                                                 {
  std::string fmt, hdr;
  //
  if ( 0 != service ) {
    fmt += " | %|-4.4s| "; // index
    hdr += " | %|-4.4s| "; // index
  }
  //
  fmt += "| %|-18.18s|";  // the name
  fmt += "| %|12d| ";     // PDGID
  fmt += "| %|=4s| ";     // charge
  fmt += "| %|13.8g|";    // mass
  fmt += " %|-3s| ";      // unit
  fmt += "| %|13.8g|";    // lifetime/gamma/ctau
  fmt += " %|-3s| ";      // unit
  fmt += "| %|=10.4g| ";  // maxwidth
  fmt += "| %|=20.20s|";  // evtgen name
  fmt += " | %|=10d| |";  // pythia ID
  fmt += " %|=20.20d| |"; // antiparticle
  //
  hdr += "| %|=18.18s|";  // the name
  hdr += "| %|=12.12s| "; // PDGID
  hdr += "| %|=4.4s| ";   // charge
  hdr += "| %|=18.18s|";  // mass
  // hdr +=   "%|=5.5s|"     ;  // unit
  hdr += "| %|=18.18s|"; // lifetime/gamma/ctau
  // hdr +=   "%|=5.5s|"     ;  // unit
  hdr += "| %|=10.10g| ";   // maxwidth
  hdr += "| %|=20.20s|";    // evtgen name
  hdr += " | %|=10.10s| |"; // pythia ID
  hdr += " %|=20.20d| |";   // antiparticle
 
  boost::format header( hdr );
  header.exceptions( boost::io::all_error_bits ^ ( boost::io::too_many_args_bit | boost::io::too_few_args_bit ) );
 
  if ( service ) { header % "#"; }
  header % "Name" % "PdgID" % "Q" % "Mass" % "(c*)Tau/Gamma" % "MaxWidth" % "EvtGen" % "PythiaID" % "Antiparticle";
 
  const std::string hdrLine = header.str();
  const std::string dashes  = " " + std::string( hdrLine.size() - 1, '-' );
  stream << dashes << std::endl << hdrLine << std::endl << dashes << std::endl;
 
  for ( auto ip = particles.begin(); particles.end() != ip; ++ip ) {
    const Gaudi::ParticleProperty* pp = *ip;
    if ( 0 == pp ) { continue; }
    boost::format line( fmt );
    line.exceptions( boost::io::all_error_bits ^ ( boost::io::too_many_args_bit | boost::io::too_few_args_bit ) );
    //
    if ( 0 != service ) {
      size_t index = Gaudi::ParticleProperties::index( pp, service );
      if ( 0 == index ) {
        line % "?";
      } else {
        line % index;
      }
    }
    //
    line % pp->name()               // name
        % pp->particleID().pid();   // PDG-ID
                                    //
#ifdef __INTEL_COMPILER             // Disable ICC remark
#  pragma warning( disable : 1572 ) // Floating-point equality and inequality comparisons are unreliable
#  pragma warning( push )
#endif
    if ( 0 == pp->charge() ) {
      line % "0";
    } else if ( 0.3 == pp->charge() ) {
      line % "+1/3";
    } else if ( -0.3 == pp->charge() ) {
      line % "-1/3";
    } else if ( 0.7 == pp->charge() ) {
      line % "+2/3";
    } else if ( -0.7 == pp->charge() ) {
      line % "-2/3";
    } else {
      line % pp->charge();
    }
#ifdef __INTEL_COMPILER // Re-enable ICC remark
#  pragma warning( pop )
#endif
    //
    // mass
    if ( pp->mass() < 1 * Gaudi::Units::keV ) {
      line % ( pp->mass() / Gaudi::Units::eV ) % "eV";
    } else if ( pp->mass() < 1 * Gaudi::Units::MeV ) {
      line % ( pp->mass() / Gaudi::Units::keV ) % "keV";
    } else if ( pp->mass() < 1 * Gaudi::Units::GeV ) {
      line % ( pp->mass() / Gaudi::Units::MeV ) % "MeV";
    } else if ( pp->mass() < 1 * Gaudi::Units::TeV ) {
      line % ( pp->mass() / Gaudi::Units::GeV ) % "GeV";
    } else {
      line % ( pp->mass() / Gaudi::Units::TeV ) % "TeV";
    }
    // lifetime/width/ctau
    const double ctau  = pp->ctau();
    const double width = pp->width();
    const double ltime = pp->lifetime();
    if ( 3600 * Gaudi::Units::second < ltime ) {
      line % "stable" % "";
    } else if ( 1 * Gaudi::Units::second < ltime ) {
      line % ( ltime / Gaudi::Units::second ) % "s";
    } else if ( 1 * Gaudi::Units::meter < ctau ) {
      line % ( ctau / Gaudi::Units::meter ) % "m";
    } else if ( 1 * Gaudi::Units::cm < ctau ) {
      line % ( ctau / Gaudi::Units::cm ) % "cm";
    } else if ( 0.5 * Gaudi::Units::mm < ctau ) {
      line % ( ctau / Gaudi::Units::mm ) % "mm";
    } else if ( 0.1 * Gaudi::Units::micrometer < ctau ) {
      line % ( ctau / Gaudi::Units::micrometer ) % "um";
    } else if ( 1 * Gaudi::Units::GeV < width ) {
      line % ( width / Gaudi::Units::GeV ) % "GeV";
    } else if ( 1 * Gaudi::Units::MeV < width ) {
      line % ( width / Gaudi::Units::MeV ) % "MeV";
    } else if ( 1 * Gaudi::Units::keV < width ) {
      line % ( width / Gaudi::Units::keV ) % "keV";
    } else if ( 1 * Gaudi::Units::eV < width ) {
      line % ( width / Gaudi::Units::eV ) % "eV";
    } else {
      line % ( ltime / Gaudi::Units::ns ) % "ns";
    }
    //
    line % pp->maxWidth() % pp->evtGen() % pp->pythia();
    //
    const Gaudi::ParticleProperty* const anti = pp->antiParticle();
    //
    if ( 0 == anti ) {
      line % "-";
    } else if ( anti == pp ) {
      line % "self-cc";
    } else {
      line % anti->name();
    }
    //
    stream << line << std::endl;
  }
  return stream << dashes << std::endl;
}

◆ printAsTable_() [3/3]

std::ostream & Gaudi::ParticleProperties::printAsTable_	(	const std::vector< Gaudi::ParticleID > &	particles,
		std::ostream &	stream,
		const Gaudi::Interfaces::IParticlePropertySvc *	service = 0 )

print properties in a form of the table

Parameters

particles	(INPUT) list of particles
stream	(UPDATE) the stream
service	(INPUT) pointer to particle property service

Returns: the stream

Author: Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@i.nosp@m.tep..nosp@m.ru

Date: 2010-01-04

Definition at line 358 of file ParticleProperty.cpp.

                                                                                                             {
  //
  if ( 0 == service ) { return Gaudi::Utils::toStream( particles, stream ); }
 
  // convert list of ParticleIDs to list fo Particle properties:
  std::vector<const Gaudi::ParticleProperty*> props;
  for ( std::vector<Gaudi::ParticleID>::const_iterator ipid = particles.begin(); particles.end() != ipid; ++ipid ) {
    const Gaudi::ParticleProperty* pp = service->find( *ipid );
    if ( 0 == pp ) { continue; }
    props.push_back( pp );
  }
  return printAsTable_( props, stream, service );
}

Functions

Function Documentation

◆ allProperties()

◆ byEvtGenName()

◆ byPythiaID()

◆ get() [1/2]

◆ get() [2/2]

◆ index() [1/2]

◆ index() [2/2]

◆ particle()

◆ particleID()

◆ printAsTable() [1/6]

◆ printAsTable() [2/6]

◆ printAsTable() [3/6]

◆ printAsTable() [4/6]

◆ printAsTable() [5/6]

◆ printAsTable() [6/6]

◆ printAsTable_() [1/3]

◆ printAsTable_() [2/3]

◆ printAsTable_() [3/3]