ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
MCGIDI_particle.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file MCGIDI_particle.cc
1 /*
2 # <<BEGIN-copyright>>
3 # <<END-copyright>>
4 */
5 #include <string.h>
6 #include "MCGIDI.h"
7 
8 #if defined __cplusplus
9 namespace GIDI {
10 using namespace GIDI;
11 #endif
12 
13 static int initialSizeOfList = 1000, incrementalSizeOfList = 1000;
17 /*
18 ************************************************************
19 */
21 
22  MCGIDI_particle *particle = (MCGIDI_particle *) smr_malloc2( smr, sizeof( MCGIDI_particle ), 0, "particle" );
23 
24  if( particle == NULL ) return( NULL );
25  MCGIDI_particle_initialize( smr, particle );
26  return( particle );
27 }
28 /*
29 ************************************************************
30 */
32 
33  memset( particle, 0, sizeof( MCGIDI_particle ) );
34  return( 0 );
35 }
36 /*
37 ************************************************************
38 */
40 
41  int i, j;
43 
44  for( i = 0, p = particleSortedList; i < numberOfParticles; i++, p++ ) {
45  if( *p == particle ) {
46  numberOfParticles--;
47  for( j = i; j < numberOfParticles; j++, p++ ) *p = p[1];
48  break;
49  }
50  }
51  if( particle == particleListEnd ) particleListEnd = particle->prior;
52  if( particle == particleList ) particleList = particle->next;
53  if( particle->prior != NULL ) particle->prior->next = particle->next;
54  if( particle->next != NULL ) particle->next->prior = particle->prior;
55  MCGIDI_particle_release( smr, particle );
56  smr_freeMemory( (void **) &particle );
57  return( NULL );
58 }
59 /*
60 ************************************************************
61 */
63 
64  smr_freeMemory( (void **) &(particle->name) );
65  return( 0 );
66 }
67 /*
68 ************************************************************
69 */
71 
72  while( particleList != NULL ) MCGIDI_particle_free( smr, particleList ); // Loop checking, 11.06.2015, T. Koi
74  return( 0 );
75 }
76 /*
77 ************************************************************
78 */
80 
81  int i, iCmp, min, mid, max, Z, A, m, level;
83  MCGIDI_POP *pop;
84 
85  iCmp = 0;
86  min = mid = 0;
87  max = numberOfParticles;
88  while( min != max ) { // Loop checking, 11.06.2015, T. Koi
89  mid = ( min + max ) / 2;
90  iCmp = strcmp( name, particleSortedList[mid]->name );
91  if( iCmp == 0 ) return( particleSortedList[mid] );
92  if( iCmp < 0 ) {
93  max = mid - 1;
94  if( mid == 0 ) max = 0; }
95  else {
96  min = mid + 1;
97  if( min > max ) min = max;
98  }
99  }
100  mid = min;
101  if( numberOfParticles > 0 ) {
102  iCmp = strcmp( name, particleSortedList[mid]->name );
103  if( iCmp == 0 ) return( particleSortedList[mid] );
104  if( ( iCmp < 0 ) && ( mid != 0 ) ) {
105  mid--;
106  iCmp = strcmp( name, particleSortedList[mid]->name );
107  }
108  }
109 
110  if( ( particle = MCGIDI_particle_new( smr ) ) == NULL ) return( NULL );
111  if( ( particle->name = smr_allocateCopyString( smr, name, "particle->name", __FILE__, __LINE__, __func__ ) ) == NULL ) goto err;
112  if( MCGIDI_miscNameToZAm( smr, name, &Z, &A, &m, &level ) != 0 ) goto err;
113  particle->prior = NULL;
114  particle->next = NULL;
115  particle->Z = Z;
116  particle->A = A;
117  particle->m = m;
118  if( ( pop = MCGIDI_POPs_findParticle( pops, name ) ) == NULL ) { /* This should not happend. */
119  particle->mass_MeV = MCGIDI_AMU2MeV * MCGIDI_particleMass_AMU( smr, name ); }
120  else {
121  particle->mass_MeV = pop->mass_MeV;
122  }
123  if( !smr_isOk( smr ) ) goto err;
124 
126  if( sizeOfParticleSortedList == 0 ) {
128  else {
130  }
132  "particleSortedList" ) ) == NULL ) goto err;
133  }
134 
135  if( particleList == NULL ) {
136  particle->ordinal = 0;
138  else {
139  particle->ordinal = particleListEnd->ordinal + 1;
140  particle->prior = particleListEnd;
143  }
144 
145  if( ( mid != 0 ) || ( iCmp > 0 ) ) mid++;
146  for( i = numberOfParticles; i > mid; i-- ) particleSortedList[i] = particleSortedList[i-1];
149 
150  return( particle );
151 
152 err:
153  MCGIDI_particle_free( smr, particle );
154  return( NULL );
155 }
156 /*
157 ************************************************************
158 */
160 
161  int i;
162 
163  for( i = 0; i < numberOfParticles; i++ ) printf( "%s\n", particleSortedList[i]->name );
164  return( 0 );
165 }
166 
167 #if defined __cplusplus
168 }
169 #endif
170