PoPs.cc File Reference

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include "PoPs.h"
#include "PoPs_private.h"

Include dependency graph for PoPs.cc:

Go to the source code of this file.

Classes
struct	unitConversions_s

Macros
#define	incrementalSize   1000
#define	MeV2eV   1e6
#define	MeV2keV   1e3
#define	AMU2MeV   931.494028
#define	AMU2eV   ( MeV2eV * 931.494028 )
#define	K2MeV   8.6173856922566752e-11
#define	K2eV   ( MeV2eV * K2MeV )

Typedefs
typedef struct unitConversions_s	unitConversions

Functions
static int	PoPs_particleProperIndex (int index)
static int	PoPs_sortedParticleIndex (char const *name)
static int	unitsDB_release (void)
const char *	PoPs_version (void)
int	PoPs_versionMajor (void)
int	PoPs_versionMinor (void)
int	PoPs_versionPatchLevel (void)
int	PoPs_register (void)
int	PoPs_readDatabase (statusMessageReporting *smr, char const *fileName)
int	PoPs_release (statusMessageReporting *smr)
int	PoPs_releasePrivate (statusMessageReporting *)
PoP *	PoPs_addParticleIfNeeded (statusMessageReporting *smr, PoP *pop)
PoP *	PoPs_copyAddParticleIfNeeded (statusMessageReporting *smr, PoP *pop)
PoP *	PoPs_addAliasIfNeeded (statusMessageReporting *smr, char const *name, char const *alias)
int	PoPs_numberOfParticle (void)
int	PoPs_particleIndex (char const *name)
int	PoPs_particleIndex_smr (statusMessageReporting *smr, char const *name, char const *file, int line, char const *func)
double	PoPs_getMassInUnitOf (statusMessageReporting *smr, char const *name, char const *unit)
char const *	PoPs_getName_atIndex (statusMessageReporting *smr, int index)
double	PoPs_getMassInUnitOf_atIndex (statusMessageReporting *smr, int index, char const *unit)
enum PoPs_genre	PoPs_getGenre (statusMessageReporting *smr, char const *name)
enum PoPs_genre	PoPs_getGenre_atIndex (statusMessageReporting *smr, int index)
int	PoPs_getZ_A_l (statusMessageReporting *smr, char const *name, int *Z, int *A, int *l)
int	PoPs_getZ_A_l_atIndex (statusMessageReporting *smr, int index, int *Z, int *A, int *l)
int	PoPs_hasNucleus (statusMessageReporting *smr, char const *name, int protonIsNucleus)
int	PoPs_hasNucleus_atIndex (statusMessageReporting *smr, int index, int protonIsNucleus)
char const *	PoPs_getAtomsName (statusMessageReporting *smr, char const *name)
char const *	PoPs_getAtomsName_atIndex (statusMessageReporting *smr, int index)
int	PoPs_getAtomsIndex (statusMessageReporting *smr, char const *name)
int	PoPs_getAtomsIndex_atIndex (statusMessageReporting *smr, int index)
PoP *	PoPs_getParticle_atIndex (int index)
char const *	PoPs_genreTokenToString (enum PoPs_genre genre)
void	PoPs_print (int sorted)
void	PoPs_write (FILE *f, int sorted)
PoP *	PoP_new (statusMessageReporting *smr)
int	PoP_initialize (statusMessageReporting *, PoP *pop)
int	PoP_release (PoP *pop)
PoP *	PoP_free (PoP *pop)
int	PoP_copyParticle (statusMessageReporting *smr, PoP *desc, PoP *src)
PoP *	PoP_makeParticle (statusMessageReporting *smr, enum PoPs_genre genre, char const *name, double mass, char const *massUnit)
int	PoP_setZ_A_l (statusMessageReporting *, PoP *pop, int Z, int A, int l)
int	PoP_getIndex (PoP *pop)
char const *	PoP_getName (PoP *pop)
double	PoP_getMassInUnitOf (statusMessageReporting *smr, PoP *pop, char const *unit)
PoP *	PoP_makeAlias (statusMessageReporting *smr, char const *name, char const *alias)
char const *	unitsDB_addUnitIfNeeded (statusMessageReporting *smr, char const *unit)
int	unitsDB_index (statusMessageReporting *, char const *unit)
char const *	unitsDB_stringFromIndex (statusMessageReporting *smr, int index)
int	PoPs_unitConversionRatio (char const *_from, char const *_to, double *ratio)

Variables
int	PoPs_smr_ID = smr_unknownID
static int	referenceCount = 0
static char	versionStr [64] = ""
static unitsDB	unitsRoot = { 0, 0, NULL }
static PoPs	popsRoot = { 0, 0, NULL, NULL }
static unitConversions	conversions []
static char const *	PoPs_genreStrings [] = { "invalid", "unknown", "alias", "photon", "lepton", "quark", "meson", "baryon", "nucleus", "atom" }

Macro Definition Documentation

#define AMU2eV   ( MeV2eV * 931.494028 )

Definition at line 23 of file PoPs.cc.

#define AMU2MeV   931.494028

Definition at line 22 of file PoPs.cc.

#define incrementalSize   1000

Definition at line 18 of file PoPs.cc.

#define K2eV   ( MeV2eV * K2MeV )

Definition at line 25 of file PoPs.cc.

#define K2MeV   8.6173856922566752e-11

Definition at line 24 of file PoPs.cc.

#define MeV2eV   1e6

Definition at line 20 of file PoPs.cc.

#define MeV2keV   1e3

Definition at line 21 of file PoPs.cc.

Typedef Documentation

typedef struct unitConversions_s unitConversions

Definition at line 27 of file PoPs.cc.

Function Documentation

int PoP_copyParticle	(	statusMessageReporting *	smr,
		PoP *	desc,
		PoP *	src
	)

Definition at line 567 of file PoPs.cc.

                                                                         {

    desc->index = -1;
    desc->properIndex = src->properIndex;
    desc->aliasIndex = src->aliasIndex;
    desc->genre = src->genre;
    if( ( desc->name = smr_allocateCopyString2( smr, src->name, "desc->name" ) ) == NULL ) return( 1 );
    desc->Z = src->Z;
    desc->A = src->A;
    desc->l = src->l;
    desc->mass = src->mass;
    desc->massUnit = src->massUnit;

    return( 0 );
}

Here is the caller graph for this function:

PoP* PoP_free

(

PoP *

pop

)

Definition at line 558 of file PoPs.cc.

                          {

    PoP_release( pop );
    smr_freeMemory( (void **) &pop );
    return( NULL );
}

Here is the call graph for this function:

Here is the caller graph for this function:

int PoP_getIndex

(

PoP *

pop

)

Definition at line 612 of file PoPs.cc.

                             {

    return( pop->index );
}

double PoP_getMassInUnitOf	(	statusMessageReporting *	smr,
		PoP *	pop,
		char const *	unit
	)

Definition at line 626 of file PoPs.cc.

                                                                                      {

    double mass = -1., ratio;
    /* PoP *pop_ = pop;*/

    /*if( pop->genre == PoPs_genre_alias ) pop_ = popsRoot.pops[PoPs_particleProperIndex( pop->index )];*/
    if( PoPs_unitConversionRatio( pop->massUnit, unit, &ratio ) != 0 ) {
        smr_setReportError2( smr, PoPs_smr_ID, PoPs_errorToken_badUnitConversion, "could not convert unit '%s' to '%s'", pop->massUnit, unit ); }
    else {
        mass = pop->mass * ratio;
    }

    return( mass );
}

Here is the call graph for this function:

Here is the caller graph for this function:

char const* PoP_getName

(

PoP *

pop

)

Definition at line 619 of file PoPs.cc.

                                    {

    return( pop->name );
}

int PoP_initialize	(	statusMessageReporting *	,
		PoP *	pop
	)

Definition at line 533 of file PoPs.cc.

                                                           {

    pop->index = -1;
    pop->properIndex = -1;
    pop->aliasIndex = -1;
    pop->genre = PoPs_genre_unknown;
    pop->name = NULL;
    pop->Z = 0;
    pop->A = 0;
    pop->mass = 0.0;
    pop->massUnit = NULL;
    return( 0 );
}

Here is the caller graph for this function:

PoP* PoP_makeAlias	(	statusMessageReporting *	smr,
		char const *	name,
		char const *	alias
	)

Definition at line 647 of file PoPs.cc.

                                                                                       {

    int properIndex = PoPs_particleIndex( name ), aliasIndex = PoPs_particleIndex( alias );
    PoP *pop;

    if( properIndex < 0 ) {
        smr_setReportError2( smr, PoPs_smr_ID, PoPs_errorToken_badName, "proper particle '%s' not in PoPs for alias '%s'", name, alias );
        return( NULL );
    }
    if( aliasIndex >= 0 ) {     /* alias has already been defined. */
        PoP *truePop = popsRoot.pops[aliasIndex];

        for( pop = truePop; strcmp( alias, pop->name ); pop = popsRoot.pops[aliasIndex] ) aliasIndex = pop->aliasIndex;
        if( pop->genre != PoPs_genre_alias ) {
            smr_setReportError2( smr, PoPs_smr_ID, PoPs_errorToken_badName, "particle '%s' already in PoPs and not an alias", alias );
            return( NULL );
        }
        if( pop->properIndex != properIndex ) {
            smr_setReportError2( smr, PoPs_smr_ID, PoPs_errorToken_badName, "particle '%s' already an alias for '%s', cannot re-alias to '%s'", 
                alias, truePop->name, name );
            return( NULL );
        } }
    else {
        if( ( pop = PoP_new( smr ) ) == NULL ) return( NULL );
        if( ( pop->name = smr_allocateCopyString2( smr, alias, "name" ) ) == NULL ) {
            PoP_free( pop );
            return( NULL );
        }
        pop->properIndex = properIndex;
        pop->genre = PoPs_genre_alias;
    }
    return( pop );
}

Here is the call graph for this function:

Here is the caller graph for this function:

PoP* PoP_makeParticle	(	statusMessageReporting *	smr,
		enum PoPs_genre	genre,
		char const *	name,
		double	mass,
		char const *	massUnit
	)

Definition at line 585 of file PoPs.cc.

                                                                                                                                 {

    PoP *pop;
    
    if( ( pop = PoP_new( smr ) ) == NULL ) return( NULL );
    if( ( pop->name = smr_allocateCopyString2( smr, name, "name" ) ) == NULL ) {
        PoP_free( pop );
        return( NULL );
    }
    pop->genre = genre;
    pop->mass = mass;
    if( ( pop->massUnit = unitsDB_addUnitIfNeeded( smr, massUnit ) ) == NULL ) pop = PoP_free( pop );
    return( pop );
}

Here is the call graph for this function:

PoP* PoP_new

(

statusMessageReporting *

smr

)

Definition at line 522 of file PoPs.cc.

                                            {

    PoP *pop;

    if( ( pop = (PoP *) smr_malloc2( smr, sizeof( PoP ), 0, "pop" ) ) == NULL ) return( NULL );
    if( PoP_initialize( smr, pop ) != 0 ) pop = PoP_free( pop );
    return( pop );
}

Here is the call graph for this function:

Here is the caller graph for this function:

int PoP_release

(

PoP *

pop

)

Definition at line 549 of file PoPs.cc.

                            {

    if( pop->name != NULL ) smr_freeMemory( (void **) &(pop->name ) );
    PoP_initialize( NULL, pop );                                /* Make it clean in case someone trys to use if. */
    return( 0 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

int PoP_setZ_A_l	(	statusMessageReporting *	,
		PoP *	pop,
		int	Z,
		int	A,
		int	l
	)

Definition at line 602 of file PoPs.cc.

                                                                              {

    pop->Z = Z;
    pop->A = A;
    pop->l = l;
    return( 0 );
}

Here is the call graph for this function:

PoP* PoPs_addAliasIfNeeded	(	statusMessageReporting *	smr,
		char const *	name,
		char const *	alias
	)

Definition at line 179 of file PoPs.cc.

                                                                                               {

    PoP *pop = PoP_makeAlias( smr, name, alias );

    if( pop != NULL ) {
        if( pop->index < 0 ) {
            if( PoPs_addParticleIfNeeded( smr, pop ) == NULL ) {
                PoP_free( pop );
                return( NULL );
            }
        }
    }
    
    return( pop );
}

Here is the call graph for this function:

Here is the caller graph for this function:

PoP* PoPs_addParticleIfNeeded	(	statusMessageReporting *	smr,
		PoP *	pop
	)

Definition at line 113 of file PoPs.cc.

                                                                       {
/*
    If particle with name pop->name is already in popsRoot, returns the pointer to the existing particle.
    A NULL is returned if adding particle to popsRoot fails.
*/
    int i, index = PoPs_sortedParticleIndex( pop->name );

    if( index >= 0 ) return( popsRoot.pops[PoPs_particleProperIndex( popsRoot.sorted[index]->index )] );
    if( popsRoot.numberOfParticles == popsRoot.allocated ) {
        int size = popsRoot.allocated + incrementalSize;
        PoP **sorted, **pops = (PoP **) smr_malloc2( smr, 2 * size * sizeof( PoPs * ), 0, "pops" );

        if( pops == NULL ) return( NULL );
        sorted = &(pops[size]);
        for( i = 0; i < popsRoot.numberOfParticles; i++ ) {
            pops[i] = popsRoot.pops[i];
            sorted[i] = popsRoot.sorted[i];
        }
        smr_freeMemory( (void **) &(popsRoot.pops) );
        popsRoot.pops = pops;
        popsRoot.sorted = sorted;
        popsRoot.allocated = size;
    }
    popsRoot.pops[popsRoot.numberOfParticles] = pop;
    index = -index - 1;
    for( i = popsRoot.numberOfParticles; i > index; i-- ) popsRoot.sorted[i] = popsRoot.sorted[i-1];
    popsRoot.sorted[index] = pop;
    pop->index = popsRoot.numberOfParticles;
    popsRoot.numberOfParticles++;
    if( pop->genre == PoPs_genre_alias ) {  /* Add pop->index to end of list of particles aliased by pop->properIndex. */
        PoP *pop2;

        for( pop2 = popsRoot.pops[pop->properIndex]; pop2->aliasIndex >= 0; pop2 = popsRoot.pops[pop2->aliasIndex] ) ;
        pop2->aliasIndex = pop->index;
    }
    return( pop );
}

Here is the call graph for this function:

Here is the caller graph for this function:

PoP* PoPs_copyAddParticleIfNeeded	(	statusMessageReporting *	smr,
		PoP *	pop
	)

Definition at line 153 of file PoPs.cc.

                                                                           {
/*
    If particle with name pop->name is already in popsRoot, return the address of the existing particle.
    If particle is not in popsRoot then copy particle to a new 'PoP *', add the copied PoP to popsRoot and return its address.
    A NULL is return if particle coping fails or adding particle to popsRoot fails.
*/

    int index = PoPs_particleIndex( pop->name );
    PoP *newPoP;

    if( index >= 0 ) return( popsRoot.pops[index] );

    if( ( newPoP = (PoP *) smr_malloc2( smr, sizeof( PoP ), 0, "newPoP" ) ) == NULL ) return( NULL );
    if( PoP_copyParticle( smr, newPoP, pop ) ) {
        smr_freeMemory( (void **) &newPoP );
        return( NULL );
    }
    if( PoPs_addParticleIfNeeded( smr, newPoP ) == NULL ) {
        PoP_free( newPoP );
        return( NULL );
    }
    return( newPoP );
}

Here is the call graph for this function:

char const* PoPs_genreTokenToString

(

enum PoPs_genre

genre

)

Definition at line 446 of file PoPs.cc.

                                                             {

    if( genre < PoPs_genre_invalid ) return( NULL );
    if( genre > PoPs_genre_atom ) return( NULL );
    return( PoPs_genreStrings[genre] );
}

Here is the caller graph for this function:

int PoPs_getAtomsIndex	(	statusMessageReporting *	smr,
		char const *	name
	)

Definition at line 397 of file PoPs.cc.

                                                                        {

    int index = PoPs_particleIndex_smr( smr, name, __FILE__, __LINE__, __func__ );

    if( index < 0 ) return( index );
    return( PoPs_getAtomsIndex_atIndex( smr, index ) );
}

Here is the call graph for this function:

int PoPs_getAtomsIndex_atIndex	(	statusMessageReporting *	smr,
		int	index
	)

Definition at line 407 of file PoPs.cc.

                                                                         {

    char const *p = NULL;

    if( ( index < 0 ) || ( index >= popsRoot.numberOfParticles ) ) {
        smr_setReportError2( smr, PoPs_smr_ID, PoPs_errorToken_badIndex, "index %d not in PoPs", index );
        return( -1 );
    }

    if( popsRoot.pops[index]->genre == PoPs_genre_atom ) return( index );

    if(      strcmp( "p", popsRoot.pops[index]->name ) == 0 ) {
        p = "H1"; }
    else {
        if( popsRoot.pops[index]->genre != PoPs_genre_nucleus ) return( -1 );
        else if( strcmp( "h2", popsRoot.pops[index]->name ) == 0 ) {
            p = "H2"; }
        else if( strcmp( "h3", popsRoot.pops[index]->name ) == 0 ) {
            p = "H3"; }
        else if( strcmp( "he3", popsRoot.pops[index]->name ) == 0 ) {
            p = "He3"; }
        else if( strcmp( "he4", popsRoot.pops[index]->name ) == 0 ) {
            p = "He4";
        }
    }
    if( p != NULL ) return( PoPs_particleIndex_smr( smr, p, __FILE__, __LINE__, __func__ ) );
    return( -1 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

char const* PoPs_getAtomsName	(	statusMessageReporting *	smr,
		char const *	name
	)

Definition at line 377 of file PoPs.cc.

                                                                               {

    int index = PoPs_particleIndex_smr( smr, name, __FILE__, __LINE__, __func__ );

    if( index < 0 ) return( NULL );
    return( PoPs_getAtomsName_atIndex( smr, index ) );
}

Here is the call graph for this function:

char const* PoPs_getAtomsName_atIndex	(	statusMessageReporting *	smr,
		int	index
	)

Definition at line 387 of file PoPs.cc.

                                                                                {

    int atomIndex = PoPs_getAtomsIndex_atIndex( smr, index );

    if( atomIndex < 0 ) return( NULL );
    return( popsRoot.pops[atomIndex]->name );
}

Here is the call graph for this function:

Here is the caller graph for this function:

enum PoPs_genre PoPs_getGenre	(	statusMessageReporting *	smr,
		char const *	name
	)

Definition at line 301 of file PoPs.cc.

                                                                               {

    int index = PoPs_particleIndex_smr( smr, name, __FILE__, __LINE__, __func__ );

    if( index < 0 ) return( PoPs_genre_invalid );
    return( popsRoot.pops[index]->genre );
}

Here is the call graph for this function:

enum PoPs_genre PoPs_getGenre_atIndex	(	statusMessageReporting *	smr,
		int	index
	)

Definition at line 311 of file PoPs.cc.

                                                                                {

    enum PoPs_genre genre = PoPs_genre_invalid;

    if( ( index < 0 ) || ( index >= popsRoot.numberOfParticles ) ) {
        smr_setReportError2( smr, PoPs_smr_ID, PoPs_errorToken_badIndex, "index %d not in PoPs", index ); }
    else {
        genre = popsRoot.pops[index]->genre;
    }
    return( genre );
}

double PoPs_getMassInUnitOf	(	statusMessageReporting *	smr,
		char const *	name,
		char const *	unit
	)

Definition at line 265 of file PoPs.cc.

                                                                                               {

    int index = PoPs_particleIndex_smr( smr, name, __FILE__, __LINE__, __func__ );

    if( index < 0 ) return( -1. );
    return( PoPs_getMassInUnitOf_atIndex( smr, index, unit ) );
}

Here is the call graph for this function:

double PoPs_getMassInUnitOf_atIndex	(	statusMessageReporting *	smr,
		int	index,
		char const *	unit
	)

Definition at line 286 of file PoPs.cc.

                                                                                                {

    double mass = -1.;

    if( ( index < 0 ) || ( index >= popsRoot.numberOfParticles ) ) {
        smr_setReportError2( smr, PoPs_smr_ID, PoPs_errorToken_badIndex, "index %d not in PoPs", index ); }
    else {
        mass = PoP_getMassInUnitOf( smr, popsRoot.pops[index], unit );
    }

    return( mass );
}

Here is the call graph for this function:

Here is the caller graph for this function:

char const* PoPs_getName_atIndex	(	statusMessageReporting *	smr,
		int	index
	)

Definition at line 275 of file PoPs.cc.

                                                                           {

    if( ( index < 0 ) || ( index >= popsRoot.numberOfParticles ) ) {
        smr_setReportError2( smr, PoPs_smr_ID, PoPs_errorToken_badIndex, "index %d not in PoPs", index );
        return( NULL );
    }
    return( popsRoot.pops[index]->name );
}

Here is the caller graph for this function:

PoP* PoPs_getParticle_atIndex

(

int

index

)

Definition at line 438 of file PoPs.cc.

                                           {

    if( ( index < 0 ) || ( index >= popsRoot.numberOfParticles ) ) return( NULL );
    return( popsRoot.pops[index] );
}

int PoPs_getZ_A_l	(	statusMessageReporting *	smr,
		char const *	name,
		int *	Z,
		int *	A,
		int *	l
	)

Definition at line 325 of file PoPs.cc.

                                                                                           {

    int index = PoPs_particleIndex_smr( smr, name, __FILE__, __LINE__, __func__ );

    if( index < 0 ) return( -1 );
    return( PoPs_getZ_A_l_atIndex( smr, index, Z, A, l ) );
}

Here is the call graph for this function:

int PoPs_getZ_A_l_atIndex	(	statusMessageReporting *	smr,
		int	index,
		int *	Z,
		int *	A,
		int *	l
	)

Definition at line 335 of file PoPs.cc.

                                                                                            {

    if( ( index < 0 ) || ( index >= popsRoot.numberOfParticles ) ) {
        smr_setReportError2( smr, PoPs_smr_ID, PoPs_errorToken_badIndex, "index %d not in PoPs", index );
        return( -1 );
    }
    *Z = popsRoot.pops[index]->Z;
    *A = popsRoot.pops[index]->A;
    *l = 0;
    return( 0 );
}

Here is the caller graph for this function:

int PoPs_hasNucleus	(	statusMessageReporting *	smr,
		char const *	name,
		int	protonIsNucleus
	)

Definition at line 349 of file PoPs.cc.

                                                                                          {

    int index = PoPs_particleIndex_smr( smr, name, __FILE__, __LINE__, __func__ );

    if( index < 0 ) return( -1 );
    return( PoPs_hasNucleus_atIndex( smr, index, protonIsNucleus ) );
}

Here is the call graph for this function:

Here is the caller graph for this function:

int PoPs_hasNucleus_atIndex	(	statusMessageReporting *	smr,
		int	index,
		int	protonIsNucleus
	)

Definition at line 359 of file PoPs.cc.

                                                                                           {
/*
*   If an error is encountered, a negative value is returned. A value greater than 0 means the particle
*   contains a nucleus (is an atom, ion or nucleus). Otherwise, a 0 is returned.
*/
    if( ( index < 0 ) || ( index >= popsRoot.numberOfParticles ) ) {
        smr_setReportError2( smr, PoPs_smr_ID, PoPs_errorToken_badIndex, "index %d not in PoPs", index );
        return( -1 );
    }
    if( ( popsRoot.pops[index]->genre == PoPs_genre_nucleus ) || ( popsRoot.pops[index]->genre == PoPs_genre_atom ) ) return( 1 );
    if( protonIsNucleus ) {
        if( strcmp( "p", popsRoot.pops[index]->name ) == 0 ) return( 1 );
    }
    return( 0 );
}

Here is the caller graph for this function:

int PoPs_numberOfParticle

(

void

)

Definition at line 197 of file PoPs.cc.

                                  {

    return( popsRoot.numberOfParticles );
}

int PoPs_particleIndex

(

char const *

name

)

Definition at line 204 of file PoPs.cc.

                                           {
/*
    A negative number is return if particle is not in popsRoot. Else, the Id of the real (not aliased) particle is returned.
*/
    int index = PoPs_sortedParticleIndex( name );

    if( index >= 0 ) index = PoPs_particleProperIndex( popsRoot.sorted[index]->index );
    return( index );
}

Here is the call graph for this function:

Here is the caller graph for this function:

int PoPs_particleIndex_smr	(	statusMessageReporting *	smr,
		char const *	name,
		char const *	file,
		int	line,
		char const *	func
	)

Definition at line 216 of file PoPs.cc.

                                                                                                                          {

    int index = PoPs_particleIndex( name );

    if( index < 0 )
        smr_setReportError( smr, NULL, file, line, func, PoPs_smr_ID, PoPs_errorToken_badName, "particle '%s' not in PoPs", name );
    return( index );
}

Here is the call graph for this function:

Here is the caller graph for this function:

static int PoPs_particleProperIndex ( int  index )

static

Definition at line 227 of file PoPs.cc.

                                                 {

    while( popsRoot.pops[index]->properIndex >= 0 ) index = popsRoot.pops[index]->properIndex;  /* For alias particles. */  // Loop checking, 11.05.2015, T. Koi
    return( index );
}

Here is the caller graph for this function:

void PoPs_print

(

int

sorted

)

Definition at line 455 of file PoPs.cc.

                              {

    PoPs_write( stdout, sorted );
}

Here is the call graph for this function:

int PoPs_readDatabase	(	statusMessageReporting *	smr,
		char const *	fileName
	)

Definition at line 81 of file PoPs.cc.

                                                                           {

    return( PoPs_particleReadDatabase( smr, fileName ) );
}

Here is the call graph for this function:

int PoPs_register

(

void

)

Definition at line 73 of file PoPs.cc.

                          {

    if( referenceCount < 0 ) return( -1 );
    return( ++referenceCount );
}

int PoPs_release

(

statusMessageReporting *

smr

)

Definition at line 88 of file PoPs.cc.

                                                {

    referenceCount--;
    if( referenceCount != 0 ) return( referenceCount );
    PoPs_releasePrivate( smr );
    return( 0 );
}

Here is the call graph for this function:

int PoPs_releasePrivate

(

statusMessageReporting *

)

Definition at line 98 of file PoPs.cc.

                                                      {

    int i;

    for( i = 0; i < popsRoot.numberOfParticles; i++ ) PoP_free( popsRoot.pops[i] );
    smr_freeMemory( (void **) &(popsRoot.pops) );
    popsRoot.sorted = NULL;
    popsRoot.numberOfParticles = 0;
    popsRoot.allocated = 0;
    unitsDB_release( );
    return( 0 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

static int PoPs_sortedParticleIndex ( char const *  name )

static

Definition at line 235 of file PoPs.cc.

                                                        {
/*
    If name is a particle in popsRoot, its index in the sorted list is returned; otherwise,
    a negative number is returned. For a particle not found, its index would be -returnValue + 1 if added;
*/
    int low = 0, mid, high = popsRoot.numberOfParticles, iCmp;

    if( high == 0 ) return( -1 );
    while( ( high - low ) > 1 ) {
        mid = ( low + high ) >> 1;
        iCmp = strcmp( name, popsRoot.sorted[mid]->name );
        if( iCmp == 0 ) return( mid );
        if( iCmp > 0 ) {
            low = mid; }
        else {
            high = mid;
        }
    }  // Loop checking, 11.05.2015, T. Koi
    if( high == 1 ) {           /* First point is not checked as loop exits when ( high = 1 ) - ( low = 0 ) <= 1 ). */
        if( !strcmp( name, popsRoot.sorted[0]->name ) ) return( 0 );            /* First name is a match. */
        if( strcmp( name, popsRoot.sorted[0]->name ) < 0 ) return( -1 );        /* name is less than first name. */
    }
    if( high < popsRoot.numberOfParticles ) {
        if( strcmp( name, popsRoot.sorted[high]->name ) == 0 ) return( high );
    }
    return( -high - 1 );
}

Here is the caller graph for this function:

int PoPs_unitConversionRatio	(	char const *	_from,
		char const *	_to,
		double *	ratio
	)

Definition at line 748 of file PoPs.cc.

                                                                                  {

    int i, n = sizeof( conversions ) / sizeof( conversions[0] );

    *ratio = 1.;
    if( strcmp( _from, _to ) == 0 ) return( 0 );
    for( i = 0; i < n; i++ ) {
        if( strcmp( conversions[i]._from, _from ) == 0 ) {
            if( strcmp( conversions[i]._to, _to ) == 0 ) {
                *ratio = conversions[i].ratio;
                return( 0 );
            } }
        else if( strcmp( conversions[i]._to, _from ) == 0 ) {
            if( strcmp( conversions[i]._from, _to ) == 0 ) {
                *ratio = 1. / conversions[i].ratio;
                return( 0 );
            }
        }
    }
    return( 1 );
}

Here is the caller graph for this function:

const char* PoPs_version

(

void

)

Definition at line 59 of file PoPs.cc.

                                 {

    if( versionStr[0] == 0 ) sprintf( versionStr, "PoPs version %d.%d.%d", POPS_VERSION_MAJOR, POPS_VERSION_MINOR, POPS_VERSION_PATCHLEVEL );
    return( versionStr );
}

int PoPs_versionMajor

(

void

)

Definition at line 67 of file PoPs.cc.

{ return( POPS_VERSION_MAJOR ); }

int PoPs_versionMinor

(

void

)

Definition at line 68 of file PoPs.cc.

{ return( POPS_VERSION_MINOR ); }

int PoPs_versionPatchLevel

(

void

)

Definition at line 69 of file PoPs.cc.

{ return( POPS_VERSION_PATCHLEVEL ); }

void PoPs_write	(	FILE *	f,
		int	sorted
	)

Definition at line 462 of file PoPs.cc.

                                       {

    int i1, properIndex;
    PoP *pop;

    fprintf( f, "Mass units: number of units = %d\n", unitsRoot.numberOfUnits );
    for( i1 = 0; i1 < unitsRoot.numberOfUnits; i1++ ) {
        fprintf( f, " %s", unitsRoot.unsorted[i1] );
    }
    fprintf( f, "\n\n" );
    fprintf( f, "Particles: number of particles = %d\n", popsRoot.numberOfParticles );
    fprintf( f, " name                      index   genre            mass             hasNucleus    alias info\n" );
    fprintf( f, "                                                                           Z   A l\n" );
    fprintf( f, " --------------------------------------------------------------------------------------------\n" );
    for( i1 = 0; i1 < popsRoot.numberOfParticles; i1++ ) {
        if( sorted ) {
            pop = popsRoot.sorted[i1]; }
        else {
            pop = popsRoot.pops[i1];
        }
        properIndex = PoPs_particleProperIndex( pop->index );
        fprintf( f, " %-24s %6d   %-10s %15.8e %-6s", pop->name, pop->index, PoPs_genreTokenToString( pop->genre ), 
            popsRoot.pops[properIndex]->mass, popsRoot.pops[properIndex]->massUnit );
        if( PoPs_hasNucleus( NULL, pop->name, 0 ) ) {
            fprintf( f, " T" ); }
        else {
            fprintf( f, "  " );
        }
        if( PoPs_hasNucleus( NULL, pop->name, 1 ) ) {
            fprintf( f, " T" ); }
        else {
            fprintf( f, "  " );
        }
        if( pop->Z + pop->A > 0 ) {
            fprintf( f, " %3d %3d", pop->Z, pop->A );
            if( pop->l > 0 ) {
                fprintf( f, " %d", pop->l ); }
            else {
                fprintf( f, "  " );
            } }
        else {
            fprintf( f, "          " );
        }
        if( pop->genre == PoPs_genre_alias ) {
            fprintf( f, " %s (%d)", popsRoot.pops[properIndex]->name, popsRoot.pops[properIndex]->index ); }
        else {
            int aliasIndex;

            for( aliasIndex = pop->aliasIndex; aliasIndex >= 0; aliasIndex = popsRoot.pops[aliasIndex]->aliasIndex ) fprintf( f, " %d", aliasIndex );
        }
        fprintf( f, "\n" );
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

char const* unitsDB_addUnitIfNeeded	(	statusMessageReporting *	smr,
		char const *	unit
	)

Definition at line 700 of file PoPs.cc.

                                                                                     {

    int i;

    for( i = 0; i < unitsRoot.numberOfUnits; i++ ) {
        if( strcmp(  unit, unitsRoot.unsorted[i] ) == 0 ) return( unitsRoot.unsorted[i] );
    }
    if( unitsRoot.numberOfUnits == unitsRoot.allocated ) {
        int size = unitsRoot.allocated + 20;
        char const **unsorted = (char const **) smr_malloc2( smr, size * sizeof( char * ), 0, "unsorted" );

        if( unsorted == NULL ) return( NULL );
        for( i = 0; i < unitsRoot.numberOfUnits; i++ ) unsorted[i] = unitsRoot.unsorted[i];
        smr_freeMemory( (void **) &(unitsRoot.unsorted) );
        unitsRoot.unsorted = unsorted;
        unitsRoot.allocated = size;
    }
    if( ( unitsRoot.unsorted[unitsRoot.numberOfUnits] = smr_allocateCopyString2( smr, unit, "unitsRoot.unsorted[unitsRoot.numberOfUnits]" ) ) == NULL ) 
        return( NULL );
    unitsRoot.numberOfUnits++;
    return( unitsRoot.unsorted[unitsRoot.numberOfUnits - 1] );
}

Here is the call graph for this function:

Here is the caller graph for this function:

int unitsDB_index	(	statusMessageReporting *	,
		char const *	unit
	)

Definition at line 725 of file PoPs.cc.

                                                                  {

    int i;

    for( i = 0; i < unitsRoot.numberOfUnits; i++ ) {
        if( !strcmp( unit, unitsRoot.unsorted[i] ) ) return( i );
    }
    return( -1 );
}

static int unitsDB_release ( void  )

static

Definition at line 687 of file PoPs.cc.

                                   {

    int i;

    for( i = 0; i < unitsRoot.numberOfUnits; i++ ) smr_freeMemory( (void **) &(unitsRoot.unsorted[i]) );
    smr_freeMemory( (void **) &(unitsRoot.unsorted) );
    unitsRoot.numberOfUnits = 0;
    unitsRoot.allocated = 0;
    return( 0 );
}

Here is the call graph for this function:

Here is the caller graph for this function:

char const* unitsDB_stringFromIndex	(	statusMessageReporting *	smr,
		int	index
	)

Definition at line 737 of file PoPs.cc.

                                                                              {

    if( ( index < 0 ) || ( index >= unitsRoot.numberOfUnits ) ) {
        smr_setReportError2( smr, PoPs_smr_ID, 1, "index = %d out of baounds [0 to %d)", index, unitsRoot.numberOfUnits );
        return( NULL );
    }
    return( unitsRoot.unsorted[index] );
}

Variable Documentation

unitConversions conversions[]

static

Initial value:

= { { "amu", "eV/c**2", AMU2eV }, { "amu", "MeV/c**2", AMU2MeV }, { "MeV/c**2", "eV/c**2", MeV2eV }, 
    { "MeV", "eV", MeV2eV }, { "MeV", "keV", MeV2keV }, { "K", "MeV", K2MeV }, { "K", "eV", K2eV } }

Definition at line 48 of file PoPs.cc.

char const* PoPs_genreStrings[] = { "invalid", "unknown", "alias", "photon", "lepton", "quark", "meson", "baryon", "nucleus", "atom" }

static

Definition at line 51 of file PoPs.cc.

int PoPs_smr_ID = smr_unknownID

Definition at line 35 of file PoPs.cc.

PoPs popsRoot = { 0, 0, NULL, NULL }

static

Definition at line 43 of file PoPs.cc.

int referenceCount = 0

static

Definition at line 36 of file PoPs.cc.

unitsDB unitsRoot = { 0, 0, NULL }

static

Definition at line 42 of file PoPs.cc.

char versionStr[64] = ""

static

Definition at line 37 of file PoPs.cc.