MCGIDI_pop.cc File Reference

#include <string.h>
#include <PoPs.h>
#include "MCGIDI.h"

Go to the source code of this file.

Functions
MCGIDI_POPs *	MCGIDI_POPs_new (statusMessageReporting *smr, int size)
int	MCGIDI_POPs_initial (statusMessageReporting *, MCGIDI_POPs *pops, int size)
void *	MCGIDI_POPs_free (MCGIDI_POPs *pops)
int	MCGIDI_POPs_release (MCGIDI_POPs *pops)
MCGIDI_POP *	MCGIDI_POPs_addParticleIfNeeded (statusMessageReporting *smr, MCGIDI_POPs *pops, char const *name, double mass_MeV, double level_MeV, MCGIDI_POP *parent, int globalParticle)
int	MCGIDI_POPs_findParticleIndex (MCGIDI_POPs *pops, char const *name)
MCGIDI_POP *	MCGIDI_POPs_findParticle (MCGIDI_POPs *pops, char const *name)
void	MCGIDI_POPs_writeSortedList (MCGIDI_POPs *pops, FILE *f)
void	MCGIDI_POPs_printSortedList (MCGIDI_POPs *pops)
MCGIDI_POP *	MCGIDI_POP_new (statusMessageReporting *smr, char const *name, double mass_MeV, double level_MeV, MCGIDI_POP *parent)
MCGIDI_POP *	MCGIDI_POP_free (MCGIDI_POP *pop)
MCGIDI_POP *	MCGIDI_POP_release (MCGIDI_POP *pop)
double	MCGIDI_POP_getMass_MeV (MCGIDI_POP *pop)

Function Documentation

MCGIDI_POP_free()

MCGIDI_POP * MCGIDI_POP_free

(

MCGIDI_POP *

pop

)

Definition at line 191 of file MCGIDI_pop.cc.

                                               {
 
    if( pop == NULL ) return( NULL );
    MCGIDI_POP_release( pop );
    smr_freeMemory( (void **) &pop );
    return( NULL );
}

Referenced by MCGIDI_POPs_release().

MCGIDI_POP_getMass_MeV()

double MCGIDI_POP_getMass_MeV

(

MCGIDI_POP *

pop

)

Definition at line 212 of file MCGIDI_pop.cc.

                                                 {
 
    return( pop->mass_MeV );
}

Referenced by MCGIDI_product_getMass_MeV(), MCGIDI_target_heated_getProjectileMass_MeV(), and MCGIDI_target_heated_getTargetMass_MeV().

MCGIDI_POP_new()

MCGIDI_POP * MCGIDI_POP_new	(	statusMessageReporting *	smr,
		char const *	name,
		double	mass_MeV,
		double	level_MeV,
		MCGIDI_POP *	parent )

Definition at line 165 of file MCGIDI_pop.cc.

                                                                                                                                   {
 
    int Z, A, m, level;
    MCGIDI_POP *pop = (MCGIDI_POP *) smr_malloc2( smr, sizeof( MCGIDI_POP ), 0, "pop" );
 
    if( pop == NULL ) return( NULL );
    pop->next = NULL;
    pop->parent = parent;
    if( ( pop->name = smr_allocateCopyString2( smr, name, "pop->name" ) ) == NULL ) {
        smr_freeMemory( (void **) &pop );
        return( NULL );
    }
    MCGIDI_miscNameToZAm( smr, name, &Z, &A, &m, &level );
    pop->Z = Z;
    pop->A = A;
    pop->level = level;
    pop->m = m;
    pop->mass_MeV = mass_MeV;
    pop->level_MeV = level_MeV;
    pop->numberOfGammaBranchs = 0;
    pop->gammas = NULL;
    return( pop );
}

Referenced by MCGIDI_POPs_addParticleIfNeeded().

MCGIDI_POP_release()

MCGIDI_POP * MCGIDI_POP_release

(

MCGIDI_POP *

pop

)

Definition at line 201 of file MCGIDI_pop.cc.

                                                  {
 
    if( pop == NULL ) return( NULL );
    smr_freeMemory( (void **) &(pop->name) );
    pop->numberOfGammaBranchs = 0;
    if( pop->gammas != NULL ) smr_freeMemory( (void **) &(pop->gammas) );
    return( NULL );
}

Referenced by MCGIDI_POP_free().

MCGIDI_POPs_addParticleIfNeeded()

MCGIDI_POP * MCGIDI_POPs_addParticleIfNeeded	(	statusMessageReporting *	smr,
		MCGIDI_POPs *	pops,
		char const *	name,
		double	mass_MeV,
		double	level_MeV,
		MCGIDI_POP *	parent,
		int	globalParticle )

Definition at line 67 of file MCGIDI_pop.cc.

                                                                   {
 
    int i, index;
    MCGIDI_POP *pop;
 
    if( ( index = MCGIDI_POPs_findParticleIndex( pops, name ) ) >= 0 ) return( pops->sorted[index] );
    if( pops->size == pops->numberOfPOPs ) {
        int size = pops->size + pops->increment;
        MCGIDI_POP **sorted = (MCGIDI_POP **) smr_malloc2( smr, size * sizeof( MCGIDI_POP * ), 0, "sorted" );
 
        if( sorted == NULL ) return( NULL );
        for( i = 0; i < pops->numberOfPOPs; i++ ) sorted[i] = pops->sorted[i];
        smr_freeMemory( (void **) &(pops->sorted) );
        pops->sorted = sorted;
        pops->size = size;
    }
    index = -index - 1;
    if( ( pop = MCGIDI_POP_new( smr, name, mass_MeV, level_MeV, parent ) ) == NULL ) return( NULL );
    for( i = pops->numberOfPOPs; i > index; i-- ) pops->sorted[i] = pops->sorted[i-1];
    pops->sorted[index] = pop;
    if( pops->first == NULL ) {
        pops->first = pop; }
    else {
        pops->last->next = pop;
    }
    pops->last = pop;
    pops->numberOfPOPs++;
    pop->globalPoPsIndex = -1;
    if( globalParticle ) {
        if( ( pop->globalPoPsIndex = lPoPs_addParticleIfNeeded( smr, name, "LLNL" ) ) < 0 ) return( NULL );
    }
    return( pop );
}

MCGIDI_POPs_findParticle()

MCGIDI_POP * MCGIDI_POPs_findParticle	(	MCGIDI_POPs *	pops,
		char const *	name )

Definition at line 133 of file MCGIDI_pop.cc.

                                                                            {
 
    int index = MCGIDI_POPs_findParticleIndex( pops, name );
 
    if( index < 0 ) return( NULL );
    return( pops->sorted[index] );
}

Referenced by MCGIDI_particle_getInternalID(), MCGIDI_product_parseFromTOM(), and MCGIDI_target_heated_read().

MCGIDI_POPs_findParticleIndex()

int MCGIDI_POPs_findParticleIndex	(	MCGIDI_POPs *	pops,
		char const *	name )

Definition at line 104 of file MCGIDI_pop.cc.

                                                                         {
 
    int iCmp = 0, min = 0, mid = 0, max = pops->numberOfPOPs;
 
    if( max == 0 ) return( -1 );
    while( ( max - min ) > 1 ) {
        mid = ( min + max ) / 2;
        iCmp = strcmp( name, pops->sorted[mid]->name );
        if( iCmp == 0 ) return( mid );
        if( iCmp < 0 ) {
            max = mid; }
        else {
            min = mid;
        }
    }  // Loop checking, 11.05.2015, T. Koi
    if( max == 1 ) {            /* First point is not checked as loop exits when ( max = 1 ) - ( min = 0 ) !> 1 ). */
        if( strcmp( name, pops->sorted[0]->name ) == 0 ) return( 0 );
    }
    if( max < pops->numberOfPOPs ) {
        if( strcmp( name, pops->sorted[max]->name ) == 0 ) return( max );
    }
    if( max == 1 ) {
        if( strcmp( name, pops->sorted[0]->name ) < 0 ) return( -1 );
    }
    return( -max - 1 );
}

Referenced by MCGIDI_POPs_addParticleIfNeeded(), and MCGIDI_POPs_findParticle().

MCGIDI_POPs_free()

void * MCGIDI_POPs_free

(

MCGIDI_POPs *

pops

)

Definition at line 41 of file MCGIDI_pop.cc.

                                            {
 
    if( pops == NULL ) return( NULL );
    MCGIDI_POPs_release( pops );
    smr_freeMemory( (void **) &pops );
    return( NULL );
}

MCGIDI_POPs_initial()

int MCGIDI_POPs_initial	(	statusMessageReporting *	smr,
		MCGIDI_POPs *	pops,
		int	size )

Definition at line 30 of file MCGIDI_pop.cc.

                                                                                   {
 
    memset( pops, 0, sizeof( MCGIDI_POPs ) );
    if( size < 10 ) size = 10;
    pops->increment = size;
    
    return( 0 );
}

Referenced by MCGIDI_POPs_new(), MCGIDI_POPs_release(), and MCGIDI_target_heated_initialize().

MCGIDI_POPs_new()

MCGIDI_POPs * MCGIDI_POPs_new	(	statusMessageReporting *	smr,
		int	size )

Definition at line 19 of file MCGIDI_pop.cc.

                                                                      {
 
    MCGIDI_POPs *pops;
 
    if( ( pops = (MCGIDI_POPs *) smr_malloc2( smr, sizeof( MCGIDI_POPs ), 0, "pops->sorted" ) ) == NULL ) return( NULL );
    if( MCGIDI_POPs_initial( smr, pops, size ) ) smr_freeMemory( (void **) &pops );
    return( pops );
}

MCGIDI_POPs_printSortedList()

void MCGIDI_POPs_printSortedList

(

MCGIDI_POPs *

pops

)

Definition at line 153 of file MCGIDI_pop.cc.

                                                      {
 
    MCGIDI_POPs_writeSortedList( pops, stdout );
}

MCGIDI_POPs_release()

int MCGIDI_POPs_release

(

MCGIDI_POPs *

pops

)

Definition at line 51 of file MCGIDI_pop.cc.

                                             {
 
    MCGIDI_POP *pop, *next;
 
    if( pops == NULL ) return( 0 );
    for( pop = pops->first; pop != NULL; pop = next ) {
        next = pop->next;
        MCGIDI_POP_free( pop );
    }
    smr_freeMemory( (void **) &(pops->sorted) );
    MCGIDI_POPs_initial( NULL, pops, 0 );
    return( 0 );
}

Referenced by MCGIDI_POPs_free(), and MCGIDI_target_heated_release().

MCGIDI_POPs_writeSortedList()

void MCGIDI_POPs_writeSortedList	(	MCGIDI_POPs *	pops,
		FILE *	f )

Definition at line 143 of file MCGIDI_pop.cc.

                                                               {
 
    int i;
 
    fprintf( f, "POPs Information: n = %d\n", pops->numberOfPOPs );
    for( i = 0; i < pops->numberOfPOPs; i++ ) fprintf( f, "    %-20s  %e\n", pops->sorted[i]->name, pops->sorted[i]->mass_MeV );
}

Referenced by MCGIDI_POPs_printSortedList().